如何使用ggplot2调整堆叠柱状图的误差棒位置

问题根因 #

误差棒错位的核心原因是堆叠柱状图的柱子默认按position = "stack"累加高度，但误差棒没有配置堆叠位置参数，默认按原始分组均值的绝对高度绘制，没有叠加下层柱子的高度。

方案1：最简代码修改（无需预处理数据） #

直接给误差棒的stat_summary添加position = position_stack(vjust = 0.5)参数，让误差棒跟随柱子自动堆叠对齐：

ggplot(test_1, aes(ID,  value , fill = type)) +
  theme_bw() +
  stat_summary(
    geom = "col",
    fun = "mean",
    position = "stack",
    width = 0.6,
    col = "black",
    size = 1
  ) +
  stat_summary(
    geom = "errorbar",
    fun.min = function(x) mean(x) - sd(x),
    fun.max = function(x) mean(x) + sd(x),
    width = 0.1,
    size = 1,
    # 新增堆叠位置参数，vjust=0.5让误差棒对齐对应柱子的中心高度
    position = position_stack(vjust = 0.5)
  )

方案2：手动预计算统计量（灵活度更高，适合复杂自定义需求） #

如果需要自定义误差逻辑或者更精准控制位置，可以先汇总每个分组的均值、标准差，再计算堆叠偏移量后绘图：

1. 预处理统计量：

library(dplyr)
library(ggplot2)

# 计算每个ID+type分组的统计量和堆叠偏移高度
test_1_stat <- test_1 %>%
  group_by(ID, type) %>%
  summarise(
    mean_val = mean(value),
    sd_val = sd(value),
    .groups = "drop"
  ) %>%
  group_by(ID) %>%
  # 按type的因子顺序排序，计算堆叠的底部偏移量
  arrange(type, .by_group = TRUE) %>%
  mutate(
    # 下层柱子的总高度
    lower_offset = cumsum(lag(mean_val, default = 0)),
    # 当前柱子的中心高度（误差棒的基准高度）
    bar_center = lower_offset + mean_val,
    err_min = bar_center - sd_val,
    err_max = bar_center + sd_val
  )

2. 基于预处理数据绘图：

ggplot(test_1_stat, aes(x = ID, fill = type)) +
  theme_bw() +
  geom_col(aes(y = mean_val), position = "stack", width = 0.6, col = "black", size = 1) +
  geom_errorbar(aes(ymin = err_min, ymax = err_max, y = bar_center), width = 0.1, size = 1)

两种方案都可以实现误差棒和对应堆叠柱子的精准对齐，方案1修改量最小，方案2适合需要进一步自定义绘图逻辑的场景。

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