使用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）