ggplot中使用自定义缩放geom_rect后保留原始X轴的方法

核心思路 #

无需修改数据缩放逻辑，通过自定义坐标轴变换或手动映射刻度与标签实现：既保留图形居中的缩放效果，又让X轴显示原始数值。

方案一：自定义坐标轴变换函数 #

适配批量处理场景，直接在scale_x_continuous中指定变换规则，无需额外修改数据列。

示例代码 #

假设你的数据集包含species、season、depth、abundance、biomass字段，代码逻辑如下：

library(ggplot2)

# 定义自定义缩放变换：实现原始值与缩放值的双向转换
custom_scale_trans <- function(max_abs) {
  trans_new(
    name = "custom_centered",
    transform = function(x) x / max_abs,  # 原始值转缩放值
    inverse = function(x) x * max_abs,    # 缩放值转原始值（用于坐标轴标签）
    breaks = function(x) pretty(x * max_abs) / max_abs  # 基于原始值生成刻度位置
  )
}

# 批量绘图循环
for (sp in unique(df$species)) {
  sp_data <- subset(df, species == sp)
  # 计算当前物种的最大绝对值（用于缩放基准）
  max_abs <- max(abs(sp_data$abundance), abs(sp_data$biomass), na.rm = TRUE)
  
  p <- ggplot(sp_data, aes(y = depth)) +
    # 绘制丰度（负方向矩形）
    geom_rect(aes(xmin = abundance, xmax = 0, ymin = depth - 0.5, ymax = depth + 0.5), 
              fill = "#4292c6", alpha = 0.7) +
    # 绘制生物量（正方向矩形）
    geom_rect(aes(xmin = 0, xmax = biomass, ymin = depth - 0.5, ymax = depth + 0.5), 
              fill = "#238b45", alpha = 0.7) +
    # 应用自定义变换，X轴显示原始值
    scale_x_continuous(
      breaks = pretty(c(-max_abs, max_abs)),  # 原始值刻度点
      labels = function(x) x,                 # 直接显示原始数值
      trans = custom_scale_trans(max_abs),
      limits = c(-1, 1)                       # 缩放后的固定范围，保证居中
    ) +
    scale_y_reverse() +  # 深度轴从上到下递增，反转Y轴符合常规展示逻辑
    facet_wrap(~season) +
    labs(x = "丰度（负值）/ 生物量（正值）", y = "深度") +
    theme_bw()
  
  ggsave(paste0(sp, "_depth_dist.png"), p, width = 10, height = 8)
}

方案二：手动映射刻度与标签 #

若自定义变换出现兼容问题，可先缩放数据，再手动将缩放后的刻度位置对应到原始值标签。

示例代码 #

library(ggplot2)

# 批量绘图循环
for (sp in unique(df$species)) {
  sp_data <- subset(df, species == sp)
  max_abs <- max(abs(sp_data$abundance), abs(sp_data$biomass), na.rm = TRUE)
  
  # 提前计算缩放后的数据列
  sp_data$abundance_scaled <- sp_data$abundance / max_abs
  sp_data$biomass_scaled <- sp_data$biomass / max_abs
  
  # 生成原始值刻度和对应的缩放后位置
  original_breaks <- pretty(c(-max_abs, max_abs))
  scaled_breaks <- original_breaks / max_abs
  
  p <- ggplot(sp_data, aes(y = depth)) +
    geom_rect(aes(xmin = abundance_scaled, xmax = 0, ymin = depth - 0.5, ymax = depth + 0.5), 
              fill = "#4292c6", alpha = 0.7) +
    geom_rect(aes(xmin = 0, xmax = biomass_scaled, ymin = depth - 0.5, ymax = depth + 0.5), 
              fill = "#238b45", alpha = 0.7) +
    scale_x_continuous(
      breaks = scaled_breaks,    # 缩放后的刻度位置
      labels = original_breaks,  # 显示原始值标签
      limits = c(-1, 1)
    ) +
    scale_y_reverse() +
    facet_wrap(~season) +
    labs(x = "丰度（负值）/ 生物量（正值）", y = "深度") +
    theme_bw()
  
  ggsave(paste0(sp, "_depth_dist.png"), p, width = 10, height = 8)
}

关键注意事项 #

• 两种方案均保留了原有的缩放居中逻辑，仅调整坐标轴显示方式
• 方案一的自定义变换更贴合ggplot语法，批量处理时代码更简洁
• 计算max_abs时需添加na.rm = TRUE，避免缺失值导致出错
• 可调整pretty()函数参数，控制刻度的数量与间隔，适配不同物种的数据范围

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