r - 如何将推算数据(w/aregImpute)添加到数据帧？

我的问题是如何将估算的数据添加到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，因此您无需提取完整的数据帧即可执行分析)