如何为ggplot2分面图的每个子图设置独立Y轴及自定义最大刻度

问题背景 #

需要为每个物种生成按Region分面的柱状图，核心要求：

• 每个分面使用独立Y轴（因不同区域记录数差异较大）
• Y轴最大刻度需对应该区域的rounded_max_count值（避免柱状图顶部贴合分面边界）

示例数据结构 #

# 基础数据框
max_counts <- data.frame(
  Region = c("Eco1", "Eco1", "Eco2", "Eco2"),
  Date = c(1, 2, 1, 2),
  max_count = c(10, 15, 8, 12),
  rounded_max_count = c(10, 15, 10, 15)
)

# 复制生成多物种数据
species_list <- c("Sp1", "Sp2", "Sp3")
max_counts_list <- list()
for (species in species_list) {
  max_counts_species <- max_counts
  max_counts_species$Species <- species
  max_counts_list[[length(max_counts_list) + 1]] <- max_counts_species
}

# 合并为总数据框
max_counts_df <- do.call(rbind, max_counts_list)

原尝试代码及问题 #

原代码无法实现每个分面Y轴自动匹配对应区域的rounded_max_count：

# 原尝试代码
plot <- ggplot(max_counts, aes(x = Date, y = max_count, fill = Region)) +
  geom_bar(stat = "identity") +
  scale_fill_manual(values = mycolors) +  
  facet_wrap(~ Region, ncol = 1, scales = "free_y") +  
  scale_x_continuous(breaks = seq(min(filtered_ecoregions$Date), max(filtered_ecoregions$Date), by = 1)) +
  scale_y_continuous(breaks = seq(0, max(max_counts$rounded_max_count), by = 5), expand = c(0, 0.01)) +
  theme_bw()

问题：要么所有分面共用统一Y轴，要么启用独立Y轴后，无法显示对应区域的最大刻度值。

解决方案 #

方法1：使用ggh4x包（推荐） #

ggh4x的facetted_pos_scales功能支持为每个分面单独设置坐标轴刻度，完美适配需求：

1. 安装并加载依赖包：

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

2. 预计算每个区域的Y轴最大刻度值：

region_y_limits <- max_counts_df %>%
  group_by(Region) %>%
  summarise(max_y = max(rounded_max_count)) %>%
  tibble::deframe() # 转换为以Region为名称、max_y为值的向量

3. 循环生成每个物种的分面图：

mycolors <- c("Eco1" = "#1f77b4", "Eco2" = "#ff7f0e") # 自定义填充色
plot_list <- list()

for (sp in unique(max_counts_df$Species)) {
  # 筛选当前物种的数据
  sp_data <- filter(max_counts_df, Species == sp)
  
  # 为每个区域生成专属Y轴刻度配置
  y_scales <- lapply(names(region_y_limits), function(reg) {
    scale_y_continuous(
      breaks = seq(0, region_y_limits[reg], by = 5),
      expand = c(0, 0.01),
      limits = c(0, region_y_limits[reg])
    )
  })
  names(y_scales) <- names(region_y_limits)
  
  # 绘制图形
  p <- ggplot(sp_data, aes(x = Date, y = max_count, fill = Region)) +
    geom_bar(stat = "identity") +
    scale_fill_manual(values = mycolors) +
    facet_wrap(~ Region, ncol = 1, scales = "free_y") +
    scale_x_continuous(breaks = unique(sp_data$Date)) +
    facetted_pos_scales(y = y_scales) + # 绑定分面与对应Y轴配置
    theme_bw() +
    labs(title = paste("物种:", sp))
  
  plot_list[[sp]] <- p
}

# 查看指定物种的图，例如Sp1
plot_list[["Sp1"]]

方法2：不依赖额外包（手动拆分拼接） #

如果不想安装第三方包，可以将每个区域的图单独绘制后拼接：

library(gridExtra)
library(dplyr)

plot_list <- list()
mycolors <- c("Eco1" = "#1f77b4", "Eco2" = "#ff7f0e")

for (sp in unique(max_counts_df$Species)) {
  sp_data <- filter(max_counts_df, Species == sp)
  region_plots <- list()
  
  for (reg in unique(sp_data$Region)) {
    reg_data <- filter(sp_data, Region == reg)
    max_y <- max(reg_data$rounded_max_count)
    
    # 绘制单个区域的图
    p <- ggplot(reg_data, aes(x = Date, y = max_count, fill = Region)) +
      geom_bar(stat = "identity") +
      scale_fill_manual(values = mycolors) +
      scale_x_continuous(breaks = unique(reg_data$Date)) +
      scale_y_continuous(breaks = seq(0, max_y, by = 5), expand = c(0, 0.01), limits = c(0, max_y)) +
      labs(title = paste("物种:", sp, "| 区域:", reg)) +
      theme_bw() +
      theme(legend.position = "none")
    
    region_plots[[reg]] <- p
  }
  
  # 拼接当前物种的所有区域图
  combined_p <- grid.arrange(grobs = region_plots, ncol = 1)
  plot_list[[sp]] <- combined_p
}

# 查看指定物种的图，例如Sp2
plot_list[["Sp2"]]

关键说明 #

• 方法1基于原生ggplot分面逻辑，通过ggh4x扩展实现自定义刻度，代码更简洁易维护
• 两种方法均确保每个分面的Y轴最大刻度严格匹配对应区域的rounded_max_count，同时保持独立Y轴
• 循环中针对每个物种单独筛选数据，避免跨物种数据干扰

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