2017-01-17 40 views
1

我的問題是我怎麼去有關添加的估算數據到quakes.missing數據幀?如何將估算數據(w/aregImpute)添加到數據框?

我在下面創建了一個可重現的示例。

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 
+1

檢查出'impute_quakes $ imputed'。 –

+0

謝謝@ChirayuChamoli。我將嘗試使用apply函數添加數據。 –

+0

你怎樣「將估算數據添加到quake.missing」數據框?你想垂直(rbind)還是水平(rbind)添加數據?或者,您是否僅僅想提取多重估算數據集? – Ahmedov

回答

1

既然你有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,所以你不需要拉出完整的數據幀,以執行分析)