Abstract:Objectives: The development of clinical prediction models is often impeded by the occurrence of missing values in the predictors. Various methods for imputing missing values before modeling have been proposed. Some of them are based on variants of multiple imputations by chained equations, while others are based on single imputation. These methods may include elements of flexible modeling or machine learning algorithms, and for some of them user-friendly software packages are available. The aim of this study was to investigate by simulation if some of these methods consistently outperform others in performance measures of clinical prediction models. Study design and setting: We simulated development and validation cohorts by mimicking observed distributions of predictors and outcome variable of a real data set. In the development cohorts, missing predictor values were created in 36 scenarios defined by the missingness mechanism and proportion of noncomplete cases. We applied three imputation algorithms that were available in R software (R Foundation for Statistical Computing, Vienna, Austria): mice, aregImpute, and missForest. These algorithms differed in their use of linear or flexible models, or random forests, the way of sampling from the predictive posterior distribution, and the generation of a single or multiple imputed data set. For multiple imputation, we also investigated the impact of the number of imputations. Logistic regression models were fitted with the simulated development cohorts before (full data analysis) and after missing value generation (complete case analysis), and with the imputed data. Prognostic model performance was measured by the scaled Brier score, c-statistic, calibration intercept and slope, and by the mean absolute prediction error evaluated in validation cohorts without missing values. Performance of full data analysis was considered as ideal. Results: None of the imputation methods achieved the model's predictive accuracy that would be obtained in case of no missingness. In general, complete case analysis yielded the worst performance, and deviation from ideal performance increased with increasing percentage of missingness and decreasing sample size. Across all scenarios and performance measures, aregImpute and mice, both with 100 imputations, resulted in highest predictive accuracy. Surprisingly, aregImpute outperformed full data analysis in achieving calibration slopes very close to one across all scenarios and outcome models. The increase of mice's performance with 100 compared to five imputations was only marginal. The differences between the imputation methods decreased with increasing sample sizes and decreasing proportion of noncomplete cases. Conclusion: In our simulation study, model calibration was more affected by the choice of the imputation method than model discrimination. While differences in model performance after using imputation methods were generally small, multiple imputation methods as mice and aregImpute that can handle linear or nonlinear associations between predictors and outcome are an attractive and reliable choice in most situations.

Auxiliary variables in multiple imputation in regression with missing X: a warning against including too many in small sample research

Analyses Using Multiple Imputation Need to Consider Missing Data in Auxiliary Variables

Multiple imputation using auxiliary imputation variables that only predict missingness can increase bias due to data missing not at random

A comparison of strategies for selecting auxiliary variables for multiple imputation

Multiple Imputation When Variables Exceed Observations: An Overview of Challenges and Solutions

Multiple imputation for missing data: fully conditional specification versus multivariate normal imputation

Random forest analysis and lasso regression outperform traditional methods in identifying missing data auxiliary variables when the MAR mechanism is nonlinear (p.s. Stop using Little's MCAR test)

A Causal View on Bias in Missing Data Imputation: The Impact of Evil Auxiliary Variables on Norming of Test Scores

Comparison of methods for imputing limited-range variables: a simulation study

Regression with missing Ys: An improved strategy for analyzing multiply imputed data

Imputation for prediction: beware of diminishing returns

Multiple imputation in data that grow over time: A comparison of three strategies

The `Why' behind including `Y' in your imputation model

Multiple Imputation for Multilevel Data with Continuous and Binary Variables

Random features models: a way to study the success of naive imputation

The performance of prognostic models depended on the choice of missing value imputation algorithm: a simulation study

An information criterion for auxiliary variable selection in incomplete data analysis

On the Relation between Prediction and Imputation Accuracy under Missing Covariates

Using auxiliary data for binomial parameter estimation with nonignorable nonresponse

How handling missing data may impact conclusions: A comparison of six different imputation methods for categorical questionnaire data

How many imputations do you need? A two-stage calculation using a quadratic rule