Find centralized, trusted content and collaborate around the technologies you use most.
Teams
Q&A for work
Connect and share knowledge within a single location that is structured and easy to search.
使用该mice软件包,我们如何检查合并分析的残差平方和?
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))?
fit
pool(fit)
lm
residuals(pool(fit))
您可以使用函数 res= sapply(fit[[4]],residuals) 提取每个输入数据集的残差
一个选项(尽管我认为它没有在任何地方发布)是对多个估算数据集的残差进行平均,apply(res,1,mean)