I recently used aregImpute to handle missing values in the dataset used for developing a risk prediction model and I computed pooled estimates by fit.mult.impute. However, I would like to extract imputed data to have a “fixed” dataset for reproducible analyses without the need for running a new imputation process. I extracted the matrix containing the final imputations from the aregImpute object and made a conversion to data.frame. The resulting object had no missing values but when I ran the model (same fitter) I got inconsistent results. Do you have any suggestions?