0

使用该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))

4

1 回答 1

0

您可以使用函数 res= sapply(fit[[4]],residuals) 提取每个输入数据集的残差

一个选项(尽管我认为它没有在任何地方发布)是对多个估算数据集的残差进行平均,apply(res,1,mean)

于 2017-08-07T14:19:37.507 回答