使用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))
?
答案 0 :(得分:0)
您可以使用该功能提取每个输入数据集的残差 res = sapply(fit [[4]],residuals)
和一个选项(虽然我不认为它已在任何地方发布)是平均多个插补数据集上的残差, 申请(RES,1,平均值)