如何为geom_col添加误差棒？ggplot2绘图问题求助

首先还原你的数据集：

Village <- as.factor(c("Beng-Gaty", "Beng-Gaty", "Ohchra", "Ohchra", "Ohtrone", "Ohtrone", "Sraevly", "Sraevly"))
Prevalence_Measure <- as.factor(c("Crude", "Adjusted", "Crude", "Adjusted", "Crude", "Adjusted", "Crude", "Adjusted"))
Prevalence <- as.numeric(c(37.7, 42.3, 39.5, 43.9, 19.4, 30, 18.2, 21.2))
Lower_CI <- as.numeric(c(0.285, 0.329, 0.316, 0.360, 0.123, 0.215, 0.124, 0.152))
Upper_CI <- as.numeric(c(0.476, 0.524, 0.477, 0.524, 0.283, 0.399, 0.252, 0.288))

Village_Prevalences_Comp <- data.frame(Village, Prevalence_Measure, Prevalence, Lower_CI, Upper_CI)

你遇到的三个问题本质是两个核心错误：

• 误差棒未设置分组偏移：geom_pointrange/geom_errorbar没有和柱状图同步设置偏移规则，且未在映射中指定分组依据，导致无法识别需要分开显示的两组误差棒。
• CI数值与患病率尺度不匹配：Lower_CI和Upper_CI是0-1的比例值，但Prevalence是百分比数值（如37.7），尺度差异导致误差棒出现在y轴0附近。

修正后的完整代码 #

library(ggplot2)
library(dplyr)

# 先把置信区间转换为百分比，和患病率统一尺度
Village_Prevalences_Comp <- Village_Prevalences_Comp %>%
  mutate(Lower_CI_pct = Lower_CI * 100,
         Upper_CI_pct = Upper_CI * 100)

# 绘制带正确误差棒的柱状图
Village_Prevalences_Comp %>%
  ggplot(aes(x = Village, y = Prevalence, fill = Prevalence_Measure)) +
  geom_col(position = position_dodge(width = 0.9)) +
  # 给误差棒设置和柱状图一致的偏移宽度，确保对齐
  geom_errorbar(aes(ymin = Lower_CI_pct, ymax = Upper_CI_pct),
                position = position_dodge(width = 0.9),
                width = 0.2) +
  labs(title = "村庄患病率对比（含置信区间）",
       x = "村庄",
       y = "患病率（%）",
       fill = "测量类型") +
  theme_minimal()

关键修正点说明 #

1. 统一数值尺度：通过mutate将置信区间的比例值乘以100转为百分比，和Prevalence的单位保持一致，解决误差棒位置过低的问题。
2. 同步偏移设置：给geom_errorbar设置position = position_dodge(width = 0.9)，宽度和geom_col的默认偏移宽度一致，确保误差棒和对应柱状图精准对齐。
3. 分组映射：在全局aes中指定fill=Prevalence_Measure，让ggplot自动识别分组逻辑，无需额外手动指定group参数。

如果想用geom_pointrange替代geom_errorbar，只需替换对应图层即可，逻辑完全一致：

Village_Prevalences_Comp %>%
  ggplot(aes(x = Village, y = Prevalence, fill = Prevalence_Measure)) +
  geom_col(position = position_dodge(width = 0.9)) +
  geom_pointrange(aes(ymin = Lower_CI_pct, ymax = Upper_CI_pct),
                  position = position_dodge(width = 0.9)) +
  labs(title = "村庄患病率对比（含置信区间）",
       x = "村庄",
       y = "患病率（%）",
       fill = "测量类型") +
  theme_minimal()

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