使用mice包完成多重插补后，如何按组获取chl的合并均值？

针对你用mice完成多重插补后，想要按hyp分组获取chl合并均值的需求，可以按照以下步骤操作：

1. 对每个插补数据集执行分组统计
   使用mice包的with()函数，遍历所有插补后的数据集，运行分组均值计算：

# 基于插补对象执行分组统计
imp_stats <- with(imp, expr = {
  library(dplyr) # 确保在每个插补环境中加载dplyr
  nhanes %>% 
    group_by(hyp) %>% 
    summarise(chl_mean = mean(chl, na.rm = TRUE))
})

2. 合并多重插补的结果
   用pool()函数按照Rubin规则合并多个插补数据集的统计结果，整合各插补结果的变异，得到合并后的估计值：

pooled_results <- pool(imp_stats)

3. 查看最终合并结果
   通过summary()函数提取合并后的均值估计、标准误及置信区间：

# 输出带置信区间的合并结果
summary(pooled_results, conf.int = TRUE)

说明 #

• with(imp, ...)会自动在每个插补后的数据集上执行指定的统计代码，返回每个插补结果的统计值
• pool()函数遵循Rubin多重插补合并规则，同时考虑了插补间变异和插补内变异，结果比单一未插补的统计更稳健

内容的提问来源于stack exchange，提问作者J. Doe