Functional form was assessed by plotting deviance residuals against each predictor variable and assessing the locally estimated scatterplot smoothing curve. Models were run with and without imputation of missing data. Both analyses gave similar results, and only complete case results are shown.Multiple imputation by chained equations was used to impute 10 complete datasets, and results were pooled.