ggplot2按分组仅展示单向误差棒（依组别调整方向）的实现方法

核心是在不手动计算分组统计量的前提下，自动识别每个时间点下的均值最高分组，动态隐藏误差棒的重叠侧：最高分组仅保留上误差棒，其余分组仅保留下误差棒，同时保留原有的正态置信区间统计逻辑，适配多值变量批量处理需求。

第一步：加载依赖包与自动标记分组 #

library(ggplot2)
library(dplyr)

# 示例数据集
Subject <- c(1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3)
Condition <- c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C")
Time <- c(1,1,1,2,2,2,1,1,1,2,2,2,1,1,1,2,2,2)
Value1 <- c(600,550,450,300,325,250,610,545,453,323,299,280,575,560,475,100,140,85)
DF1 <- data.frame (Subject, Condition, Time, Value1)

# 自动标记每个时间点的均值最高分组，无需手动计算
DF1 <- DF1 %>%
  group_by(Time, Condition) %>%
  mutate(group_mean = mean(Value1)) %>%
  group_by(Time) %>%
  mutate(is_top_group = group_mean == max(group_mean)) %>%
  ungroup()

第二步：修改后的ggplot绘图代码 #

仅需调整误差棒部分的stat_summary调用，分两次绘制单侧误差棒即可：

PublicationPlot <- ggplot(DF1, aes(Time, Value1, shape = Condition))
PublicationPlot + 
  stat_summary(fun = mean,
               geom = "point",
               size = 2,
               aes(group = Condition)) +
  stat_summary(fun = mean, 
               geom = "line",
               aes(group = Condition,
                   linetype = Condition)) +
  # 绘制最高分组的上误差棒：ymin设为均值，仅显示上半段
  stat_summary(fun.data = mean_cl_normal,
               geom = "errorbar",
               width = 0.075,
               aes(group = Condition,
                   ymin = after_stat(ifelse(DF1$is_top_group[group], y, ymin))),
               data = ~ filter(.x, is_top_group)) +
  # 绘制其余分组的下误差棒：ymax设为均值，仅显示下半段
  stat_summary(fun.data = mean_cl_normal,
               geom = "errorbar",
               width = 0.075,
               aes(group = Condition,
                   ymax = after_stat(ifelse(!DF1$is_top_group[group], y, ymax))),
               data = ~ filter(.x, !is_top_group)) +
  xlab("Measurement Times") +
  ylab("Value 1 (Units)") +
  theme(panel.grid.major = element_blank(),
        panel.grid.minor = element_blank(),
        panel.background = element_blank(),
        axis.line = element_line(color = "black"),
        legend.key = element_rect(fill = "white"),
        axis.title.x = element_text(size = 15),
        axis.text.x = element_text(size = 13),
        axis.title.y = element_text(size = 15),
        axis.text.y = element_text(size = 13),
        legend.title = element_text(size = 12), 
        legend.text = element_text(size = 12))

如果你有9个Value变量，只需用tidyr::pivot_longer将宽格式数据转为长格式，增加区分不同指标的分组列，后续在绘图时可以通过分面或颜色映射区分不同指标，标记最高分组的逻辑无需修改即可批量复用。

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