使用R中的鼠標包進行多次插補的殘差平方和

Question

使用mice包，我們如何檢查合並分析的殘差平方和？

library(mice) 
imp <- mice(nhanes, seed = 24500)
fit <- with(imp, lm(chl ~ age + bmi))
pool(fit)
summary(pool(fit))

fit包含對每個輸入數據集的分析，並合並pool(fit)合並的結果。 是否有一個命令來檢查標准lm對象的殘差平方和，類似於residuals(pool(fit)) ？

Answer 1

您可以使用函數res = sapply（fit [[4]]，residuals）為每個輸入數據集提取殘差

一種選擇（盡管我不認為它已在任何地方發布）是對多個估算數據集上的殘差求平均值，apply（res，1，mean）