如何为含多分组的ggplot线性回归散点图添加各组置信区间

问题背景 #

已完成硕士论文数据的线性回归分析，绘制了x轴含3个分类（对应3个窝数）、3个不同处理组的散点图，需要为每个处理组的每个x轴分类添加95%置信区间，现有绘图代码如下：

ClutchvsAge = ggplot(All_Pippi, aes(Clutch, Age, colour = Treatment)) + 
 geom_point() + geom_jitter(width=0.2) +
 theme(axis.text=element_text(size=11),
       axis.title=element_text(size=12)) + 
   labs(x= "Clutch Number", y = "Age") + 
   scale_x_continuous(breaks=c(1,2,3)) + 
   scale_color_manual(values=c("#228B22", "#E69F00", "#56B4E9"))

原始绘图效果：

需求是使用geom_errorbar实现置信区间绘制，不清楚如何给ymin、ymax参数赋值每个处理-窝数组合对应的置信区间。

实现方案 #

步骤1：计算每个分组的置信区间统计量 #

需要先按处理组和窝数两个维度对原始数据分组，计算每组的Age均值、标准误，再推导置信区间上下限，示例代码（依赖dplyr包）：

library(dplyr)
ci_summary <- All_Pippi %>%
  # 按两个维度分组，匹配你需要的交叉组合统计逻辑
  group_by(Treatment, Clutch) %>%
  summarise(
    age_mean = mean(Age, na.rm = TRUE),
    # 计算标准误
    se = sd(Age, na.rm = TRUE)/sqrt(n()),
    # 95%置信区间上下限，也可直接用t.test(Age)$conf.int[1]和[2]获取更精准的结果
    ci_min = age_mean - 1.96 * se,
    ci_max = age_mean + 1.96 * se
  )

步骤2：修改绘图代码新增误差条图层 #

在原有绘图代码中新增geom_errorbar层，指定使用刚才生成的汇总数据集，将对应的置信区间上下限赋值给ymin和ymax即可，同时可调整位置参数避免不同处理组的误差条重叠：

library(ggplot2)
ClutchvsAge <- ggplot(All_Pippi, aes(x = Clutch, y = Age, colour = Treatment)) + 
  # 调整抖动逻辑，避免不同处理组的散点重叠
  geom_jitter(position = position_jitterdodge(jitter.width = 0.1, dodge.width = 0.2)) +
  # 添加置信区间误差条
  geom_errorbar(data = ci_summary,
                aes(y = age_mean, ymin = ci_min, ymax = ci_max),
                width = 0.15, position = position_dodge(width = 0.2)) +
  # 可选添加每组均值点，可读性更强
  geom_point(data = ci_summary, aes(y = age_mean),
             size = 3, position = position_dodge(width = 0.2)) +
  theme(axis.text = element_text(size = 11),
        axis.title = element_text(size = 12)) + 
  labs(x = "Clutch Number", y = "Age") + 
  scale_x_continuous(breaks = c(1,2,3)) + 
  scale_color_manual(values = c("#228B22", "#E69F00", "#56B4E9"))

# 输出绘图结果
print(ClutchvsAge)

注意事项 #

• 分组汇总时必须同时包含处理组和窝数两个维度，才能得到每个交叉组合的置信区间
• 如果需要调整置信区间水平，比如改为99%，只需要把1.96替换为对应的分位数（99%对应2.58）即可
• position_dodge的宽度可以根据实际显示效果调整，避免不同处理组的元素重叠

内容的提问来源于stack exchange，提问作者Anna Olsson