如何将估算数据(w / aregImpute)添加到数据框?

How to add imputed data (w/ aregImpute) to a data frame?

我的问题是如何将推算数据添加到 quakes.missing 数据框?

我在下面创建了一个可重现的示例。

library(Hmisc)
library(missForest) #load packages

data("quakes") 
quakes

quakes.missing <- prodNA(quakes, noNA = 0.1) #create missing values

summary(is.na(quakes.missing)) #confirm that data is missing

impute_quakes <- aregImpute(~ lat + long + depth + mag + stations, data = quakes.missing, n.impute = 5)

impute_quakes

因为你有 5 个插补,你会有 5 个完整的数据框,你可以用这样的函数把它们拉出来:

fill_data <- function(impute = impute_quakes, data = quakes.missing, im = 1) {
  cbind.data.frame(impute.transcan(x = impute, 
                                   imputation = im, 
                                   data = data, 
                                   list.out = TRUE, 
                                   pr = FALSE))
 }
full_dat1 <- fill_data(im = 1)
full_dat2 <- fill_data(im = 2)
...

(另外,我相信你知道,但是 Hmisc 也有一个很棒的功能 fit.mult.impute 所以你不需要为了执行分析而拉出完整的数据框)