To examine the effect of missing data, multiple imputation by chained equations was used to impute 20 complete datasets on which the analysis was repeated and the results pooled. Predictive mean matching with five nearest neighbours was used for continuous variables and logistic regression for binary variables. Plots were examined to assess convergence and plausibility of estimates.