如何使用ggplot2制作带整体误差棒的堆叠条形图

报错原因 #

1. 你将ymin、ymax参数写到了geom_bar的美学映射里，geom_bar本身不识别这两个参数，它们是geom_errorbar的专属参数；同时你直接在代码里引用的count被R识别为dplyr的count()函数（闭包类型），不是你预期的计数值，因此触发类型转换报错。
2. 你之前按color和clarity分组后计算标准差返回全为NA，是因为每个分组下只有唯一的1条计数记录，单值无法计算标准差，你需要先确认误差的统计来源（比如重复实验的多次计数结果）。

解决思路 #

先完成数据统计汇总，计算好每个分组的计数值、误差值、堆叠对应的y轴位置，再使用预计算的数据绘图，不要依赖ggplot的内置统计逻辑自动计数。

示例代码 #

# 加载依赖包
library(ggplot2)
library(dplyr)

# 步骤1：统计每个color+clarity分组的计数值
diamonds_summary <- diamonds %>%
  group_by(color, clarity) %>%
  summarize(count = n(), .groups = "drop") 

# 步骤2：补充误差值（此处为模拟值，实际使用时替换为你真实计算的标准差）
set.seed(123)
diamonds_summary <- diamonds_summary %>%
  mutate(sd_count = runif(n(), min = 20, max = 100))

# 步骤3：计算堆叠条形图每个分段的y轴位置，用于匹配误差棒位置
diamonds_summary <- diamonds_summary %>%
  group_by(clarity) %>%
  arrange(color, .by_group = TRUE) %>%
  mutate(
    stack_y_top = cumsum(count), # 每个堆叠段的顶部y坐标
    error_ymin = stack_y_top - sd_count,
    error_ymax = stack_y_top + sd_count
  )

# 步骤4：绘图
ggplot(diamonds_summary, aes(x = clarity, y = count, fill = color)) +
  geom_col(position = "stack") + # 用geom_col绘制预计算数值的条形图
  geom_errorbar(
    aes(ymin = error_ymin, ymax = error_ymax),
    width = 0.2, colour = "black"
  )

注意事项 #

• 如果你的数据有重复样本，需要先按样本ID、color、clarity分组统计单次重复的计数，再计算多轮重复的均值和标准差，再代入上述代码即可得到真实的误差棒。
• 如果你不需要堆叠条形图，而是分组并列的条形图，只需要将position = "stack"改为position = position_dodge(width = 0.9)，同时调整误差棒的position参数为对应值即可，不需要额外计算堆叠位置。

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