Multivariable flexible modelling for estimating complete, smoothed life tables for sub-national populations