在ggplot2的scale="free"分面图中设置Y轴固定下限

我完全懂你遇到的麻烦——带着异常值的分面柱状图，想让每个分面的Y轴都从固定值（比如100）开始，同时让上限根据每个分面的数据自动计算，还不想一个个单独绘制每个Taxa。你之前试的几种方法要么被忽略，要么直接把柱子搞没了，这是因为ggplot默认的scale_y_continuous(limits)会过滤超出范围的数据，而expand_limits在scale="free_y"的分面下会被覆盖，coord_cartesian的语法你用错了而且也不兼容free_y分面。

下面给你两种可行的解决方案，都不需要单独绘制每个Taxa：

方法一：用ggh4x包自定义分面轴缩放（推荐） #

ggh4x是专门增强ggplot分面功能的包，能轻松为每个分面单独设置坐标轴范围：

首先安装并加载包（如果还没装的话）：

install.packages("ggh4x")
library(ggplot2)
library(ggh4x)
library(dplyr)

然后用以下代码处理你的数据和绘图：

# 你的原始数据集
df <- data.frame(
  Taxa=c(rep("A", 4), rep("B", 4)),
  Day=as.factor(c(1,1,2,2,1,1,2,2)),
  Treatment=c(rep(c("t1","t2"),4)),
  Values=c(100,110,120,130, 120, 130,150, 1800),
  SD = c(10,11,12,13,12,13,14,200)
)

# 提前计算每个Taxa的Y轴范围：固定下限100，上限取该组(Values+SD)的最大值并留10%余量
y_limits <- df %>%
  group_by(Taxa) %>%
  summarise(
    ymin = 100,
    ymax = max(Values + SD) * 1.1
  )

# 绘制分面图并自定义每个分面的Y轴
p1 <- ggplot(df, aes(x=Day, ymax=Values+SD, ymin=Values-SD, y=Values, fill=Treatment))+
  geom_bar(stat="identity", position=position_dodge(width=0.7), col="black", width=0.7) +
  facet_wrap(~Taxa, nrow=2) + # 这里不用scale="free_y"，因为我们自定义了缩放
  scale_y_continuous(labels = scales::scientific) +
  geom_errorbar(position=position_dodge(width=0.7), col="black", lty="solid", width=0.3) +
  # 为每个分面设置Y轴范围
  facetted_pos_scales(
    y = list(
      Taxa == "A" ~ scale_y_continuous(limits = c(y_limits$ymin[y_limits$Taxa=="A"], y_limits$ymax[y_limits$Taxa=="A"]), labels = scales::scientific),
      Taxa == "B" ~ scale_y_continuous(limits = c(y_limits$ymin[y_limits$Taxa=="B"], y_limits$ymax[y_limits$Taxa=="B"]), labels = scales::scientific)
    )
  )

p1

这个方法的好处是代码清晰，而且能完美保留所有数据（包括异常值），每个分面的Y轴下限固定，上限自动适配该组数据。

方法二：手动修改ggplot的内部对象（无需额外包） #

如果不想安装新包，可以通过ggplot_build和ggplot_gtable直接修改绘图的内部轴设置：

library(ggplot2)
library(dplyr)

# 原始数据集和基础绘图
df <- data.frame(
  Taxa=c(rep("A", 4), rep("B", 4)),
  Day=as.factor(c(1,1,2,2,1,1,2,2)),
  Treatment=c(rep(c("t1","t2"),4)),
  Values=c(100,110,120,130, 120, 130,150, 1800),
  SD = c(10,11,12,13,12,13,14,200)
)

p1 <- ggplot(df, aes(x=Day, ymax=Values+SD, ymin=Values-SD, y=Values, fill=Treatment))+
  geom_bar(stat="identity", position=position_dodge(width=0.7), col="black", width=0.7) +
  facet_wrap(~Taxa, scale="free_y", nrow=2) +
  scale_y_continuous(labels = scales::scientific) +
  geom_errorbar(position=position_dodge(width=0.7), col="black", lty="solid", width=0.3)

# 计算每个Taxa的Y轴范围
y_limits <- df %>%
  group_by(Taxa) %>%
  summarise(
    ymin = 100,
    ymax = max(Values + SD) * 1.1
  ) %>%
  arrange(Taxa) # 确保顺序和分面一致

# 手动修改绘图对象的Y轴范围
gb <- ggplot_build(p1)
for(i in seq_along(y_limits$Taxa)){
  gb$layout$panel_scales_y[[i]]$limits <- c(y_limits$ymin[i], y_limits$ymax[i])
}
gp <- ggplot_gtable(gb)

# 显示修改后的图
plot(gp)

这个方法通过直接修改ggplot生成的内部对象，替换每个分面的Y轴范围，同样能达到固定下限、自动适配上限的效果，而且不需要额外安装包。

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