r-MICE刺激后的柱平均值和SD

Question

我已经在m=18的数据集上运行了MICE，我希望用汇总估计来确定列的新平均值和SD，因为我希望将其与我原始的未调整数据集进行比较，但不太清楚如何实现这一点。我想知道是否有人对此有什么建议？

到目前为止，我已经能够使用计算该列的每个估算数据集的平均值或sd(在下面的代码中替换sd的平均值(

library(mice)
aux_vart <- mice::quickpred(
nhanes,
mincor = 0.1
)
imp <- mice::mice(nhanes, pred = aux_vart, m = 16, meth = "pmm")
give_imp_n <- with(imp, expr = mean(bmi))

但是当我尝试将这些集合在一起时，使用：

give_pool_n <- mice::pool(give_imp_n)

我收到一个错误：

summary中的错误(.，m=n((，qbar=均值(.data$estimate(，ubar=均值(.deata$std.Error^2(，：

由.data$estimate中的错误引起：！在.data中找不到列estimate。

Answer 1

pool函数似乎需要线性回归的输出。如果你将平均值建模为仅截距回归，它似乎有效：

give_imp_n <- with(imp, expr = lm(bmi~1))
mice::pool(give_imp_n)
# Class: mipo    m = 2 
#          term m estimate      ubar        b        t dfcom       df       riv    lambda       fmi
# 1 (Intercept) 2   26.292 0.5801967 0.326432 1.069845    24 3.419396 0.8439345 0.4576814 0.6266439