如何修正ggplot代码，为林分重量数据添加置信区间误差棒、散点及抖动图？

问题需求 #

针对有无防护林的林分分组展示物种总重量的均值，误差棒代表均值的置信区间，同时添加散点及抖动图，修正现有代码错误以达到目标效果。

原代码存在的问题 #

• 分组逻辑错误：group_by("Yes")/group_by("No")使用字符串而非原始变量名，导致分组无效
• 置信区间计算硬编码样本量（12、8），易出错且不灵活
• 冗余拆分数据集：没必要分别生成mean_sd_yes/mean_sd_no，可统一处理
• 散点数据处理完全错误：points_df_yes/points_df_no的写法不符合语法规范
• ggplot初始映射错误：使用不存在的mean_sd数据集，y轴映射为标准差而非均值

修正后的完整代码 #

library(tidyverse)

# 统一计算分组的均值、置信区间
summary_data <- data %>%
  group_by(`Type of plot with shelterbelt`) %>%
  summarise(
    mean_weight = mean(Weight, na.rm = TRUE),
    sd_weight = sd(Weight, na.rm = TRUE),
    n = n(),
    ci_lower = mean_weight - qt(0.975, n - 1) * sd_weight / sqrt(n),
    ci_upper = mean_weight + qt(0.975, n - 1) * sd_weight / sqrt(n)
  )

# 绘制图形
ggplot() +
  # 添加抖动散点（原始数据）
  geom_jitter(data = data, 
              aes(x = `Type of plot with shelterbelt`, y = Weight),
              width = 0.15, alpha = 0.6, color = "#636363") +
  # 添加置信区间误差棒
  geom_errorbar(data = summary_data,
                aes(x = `Type of plot with shelterbelt`, 
                    ymin = ci_lower, ymax = ci_upper),
                width = 0.2, color = "#333333") +
  # 添加均值点
  geom_point(data = summary_data,
             aes(x = `Type of plot with shelterbelt`, y = mean_weight),
             size = 3, color = "#E6550D", fill = "#E6550D", shape = 21) +
  # 图形美化
  labs(x = "是否有防护林", y = "物种总重量") +
  theme_bw()

代码说明 #

1. 数据汇总：使用group_by按防护林类型分组，自动计算每组样本量n，基于t分布计算置信区间，避免硬编码样本量
2. 图层顺序：先绘制抖动散点（底层），再叠加误差棒和均值点（上层），保证视觉层次清晰
3. 抖动散点：width控制横向抖动幅度，alpha调整透明度避免点重叠，color设置点的颜色
4. 均值点：使用带填充的圆形（shape=21）突出显示均值，与散点形成区分
5. 图形美化：使用theme_bw()生成简洁的白色背景主题，添加清晰的坐标轴标签

内容的提问来源于stack exchange，提问作者Madelein Victor