如何通过创建数据框子集在ggplot2的geom_ribbon中指定HC与HCc曲线间的填充界限

要在HC和HCc两条曲线之间添加填充区域，核心问题是你的原始数据是长格式（每个分组的数值各占一行），而geom_ribbon需要每个x轴刻度（Time）对应的上下限数值在同一行里。我们可以先整理数据，再添加填充层，具体步骤如下：

1. 整理数据，提取HC和HCc的均值并转为宽格式 #

首先从你的high数据集中筛选出HC和HCc的记录，然后用tidyr::pivot_wider把它们转换成宽格式，这样每个Time对应的HC均值（作为上限）和HCc均值（作为下限）会在同一行：

library(tidyr)
# 筛选HC和HCc的数据，转换为宽格式
ribbon_data <- high %>%
  filter(Trat %in% c("HC", "HCc")) %>%
  pivot_wider(
    id_cols = Time,
    names_from = Trat,
    values_from = mean,
    names_prefix = "mean_"
  )

2. 修改ggplot代码，添加geom_ribbon #

在你原有的ggplot代码基础上，添加geom_ribbon层，使用上面整理好的ribbon_data，指定ymax = mean_HC（HC的均值作为上限）、ymin = mean_HCc（HCc的均值作为下限），同时设置填充样式：

library(ggplot2)

ggplot(high, aes(Time, mean, color=Trat, group = Trat))+
  # 先添加填充层（放在线的下方，避免遮挡曲线）
  geom_ribbon(
    data = ribbon_data,
    aes(x = Time, ymax = mean_HC, ymin = mean_HCc),
    fill = "gray", alpha = 0.5,
    inherit.aes = FALSE  # 不继承主ggplot的aes设置，避免颜色/分组干扰
  )+
  geom_line(aes(linetype=Trat, color=Trat), size=1.5, na.rm=T)+
  labs(x = ("Time (days)"), y = bquote("Clorophyll-a"~ (µg.L^-1))) +
  scale_x_discrete(breaks=c("T1", "T2","T3", "T4", "T5","T6","T7","T8","T9","T10"),
                   labels=c("1","2","3","4","5","6","7","8","9","10"))+
  theme_classic(base_size = 15)+
  theme(plot.title = element_text(hjust = 0.5))+
  theme(legend.position="right")+
  theme(legend.title = element_blank())+
  guides(col = guide_legend(nrow = 2))+
  scale_color_grey()+
  geom_errorbar(aes(ymin=mean-sd, ymax=mean+sd), width=1, position=position_dodge(.01))+
  guides(colour=guide_legend(nrow=3))

关键说明 #

• 把geom_ribbon放在geom_line之前，这样填充区域会显示在曲线下方，不会遮挡线条。
• 使用inherit.aes = FALSE，因为主ggplot的aes包含了color=Trat，而我们的填充层不需要这个映射，避免出现不必要的分组错误。
• 如果你还想把误差棒的范围也考虑进去（比如填充HC的误差上限到HCc的误差下限），只需要在整理ribbon_data时同时提取sd，计算mean_HC + sd_HC和mean_HCc - sd_HCc之类的数值，再对应到ymax和ymin即可。

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