MICE插补后如何合并计算列的均值与标准差？

错误原因 #

mice::pool() 函数是为合并模型拟合结果（如线性回归、逻辑回归输出）设计的，要求输入对象必须包含estimate和std.error列。你直接用with(imp, mean(bmi))得到的只是各插补数据集的均值向量，没有这些必要字段，因此触发报错。

解决方案 #

1. 合并均值 #

MICE插补后的合并均值可直接取各插补数据集均值的算术平均值：

# 提取所有插补数据集的bmi均值
imp_means <- with(imp, mean(bmi))
# 计算合并后的均值
pooled_mean <- mean(imp_means)

2. 合并标准差（遵循Rubin's规则） #

合并标准差需要同时考虑插补内方差和插补间方差，不能直接取各插补集标准差的平均值，需按Rubin规则计算：

方法一：利用pool()函数

先构造包含estimate（均值）和std.error（标准差）的结果对象，再传入pool()：

# 提取每个插补数据集的bmi均值和标准差
imp_stats <- with(imp, {
  mean_val <- mean(bmi)
  sd_val <- sd(bmi)
  data.frame(estimate = mean_val, std.error = sd_val)
})

# 合并统计量
pooled_obj <- pool(imp_stats)
pooled_summary <- summary(pooled_obj)

# 提取合并后的均值和标准差
pooled_mean <- pooled_summary$estimate
pooled_sd <- sqrt(pooled_summary$ubar + pooled_summary$b + pooled_summary$b / imp$m)

方法二：手动计算（更直观）

m <- imp$m  # 插补数据集数量

# 提取各插补集的均值和方差
q_list <- with(imp, mean(bmi))
u_list <- with(imp, var(bmi))

# 计算插补内方差（各插补集方差的平均值）
ubar <- mean(u_list)
# 合并均值
qbar <- mean(q_list)
# 计算插补间方差
b <- var(q_list) * m / (m - 1)
# 总方差（含校正项）
total_var <- ubar + b + b/m
# 合并后的标准差
pooled_sd <- sqrt(total_var)

与原始数据集对比 #

计算完合并后的均值和标准差后，直接和原始未插补数据集的统计量对比即可：

# 原始数据集的bmi均值与标准差
original_mean <- mean(nhanes$bmi, na.rm = TRUE)
original_sd <- sd(nhanes$bmi, na.rm = TRUE)

# 输出对比结果
cat("原始数据集：均值=", original_mean, "，标准差=", original_sd, "\n")
cat("插补合并后：均值=", pooled_mean, "，标准差=", pooled_sd, "\n")

内容的提问来源于stack exchange，提问作者Seafra Barrett