如何在y轴对数化的ggplot图中添加矩形且不改变y轴刻度

核心思路 #

要避免矩形影响y轴刻度，关键是让矩形的y范围不参与ggplot的自动刻度计算，或者使用和原始数据匹配的对数空间极值，同时锁定y轴范围。

可行方案 #

方案1：用annotate()配合对数转换后的y极值 #

先从原始数据提取y的实际极值，转换为对数空间后传入annotate，同时用coord_cartesian()锁定y轴范围，避免自动扩展。

示例代码：

# 计算原始数据的y极值（对应你代码里的RelativeAbundance+0.001）
y_min <- min(mydata$RelativeAbundance + 0.001)
y_max <- max(mydata$RelativeAbundance + 0.001)

# 转换为log10空间的值（匹配scale_y_log10的转换规则）
log_y_min <- log10(y_min)
log_y_max <- log10(y_max)

p <- ggplot(mydata, aes(x=SampleDate, y=RelativeAbundance+0.001, colour = Taxa)) + 
  scale_y_log10() +
  geom_smooth(method="gam", se=F) +
  # 添加矩形，使用对数空间的y值
  annotate("rect", 
           xmin = as.Date("2020-01-01"),  # 替换为你的起始日期
           xmax = as.Date("2020-12-31"),  # 替换为你的结束日期
           ymin = log_y_min, 
           ymax = log_y_max,
           alpha = 0.2, fill = "grey") +
  # 锁定y轴范围与原始数据一致
  coord_cartesian(ylim = c(y_min, y_max))

方案2：用annotation_custom()配合rectGrob() #

annotation_custom()的元素不会参与坐标轴范围计算，只要避开0值，使用数据对应的对数极值即可：

library(grid)

# 将日期转换为ggplot内部使用的数值格式
x_start <- as.numeric(as.Date("2020-01-01"))
x_end <- as.numeric(as.Date("2020-12-31"))

# 获取原始y的极值
y_min <- min(mydata$RelativeAbundance + 0.001)
y_max <- max(mydata$RelativeAbundance + 0.001)

p <- ggplot(mydata, aes(x=SampleDate, y=RelativeAbundance+0.001, colour = Taxa)) + 
  scale_y_log10() +
  geom_smooth(method="gam", se=F) +
  annotation_custom(
    grob = rectGrob(gp = gpar(fill = "grey", alpha = 0.2)),
    xmin = x_start,
    xmax = x_end,
    ymin = log10(y_min),
    ymax = log10(y_max)
  )

方案3：用geom_rect()并设置inherit.aes=FALSE+锁定y轴 #

确保geom_rect不继承主图的映射关系，同时使用对数空间的y值，配合coord_cartesian锁定范围：

y_min <- min(mydata$RelativeAbundance + 0.001)
y_max <- max(mydata$RelativeAbundance + 0.001)

p <- ggplot(mydata, aes(x=SampleDate, y=RelativeAbundance+0.001, colour = Taxa)) + 
  scale_y_log10() +
  geom_smooth(method="gam", se=F) +
  geom_rect(
    inherit.aes = FALSE,
    aes(xmin = as.Date("2020-01-01"), xmax = as.Date("2020-12-31"),
        ymin = log10(y_min), ymax = log10(y_max)),
    fill = "grey", alpha = 0.2
  ) +
  coord_cartesian(ylim = c(y_min, y_max))

关键注意点 #

• 所有y值要匹配scale_y_log10()的转换规则，即使用log10处理后的数值，因为ggplot在对数刻度下是基于转换后的空间绘制元素；
• coord_cartesian()可以强制锁定y轴范围，避免矩形导致刻度自动扩展；
• annotation_custom()是最稳妥的方法之一，因为它的元素不会影响坐标轴范围的计算。

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