基于同一数据框为ggplot2箱线图添加双Y轴的技术问询

嘿，我完全懂你的痛点！当某组数据的量级远超过其他组时，单Y轴箱线图直接把小数值的特征给挤没了，根本没法看清细节。下面我就手把手教你用ggplot2实现双Y轴箱线图，同时也给你一个更稳妥的替代方案～

我先构造和你描述一致的模拟数据，方便后续演示：

library(ggplot2)

set.seed(123) # 固定随机种子，保证结果可复现
df <- data.frame(
  Species = rep(c("物种1", "物种2"), each = 30),
  Appendage = rep(c("附肢A", "附肢B", "附肢C"), each = 10, times = 2),
  SA_V = c(
    rnorm(20, mean = 6, sd = 1.2),  # 附肢A和B的SA:V（低量级）
    rnorm(20, mean = 45, sd = 4)    # 附肢C的SA:V（远高于前两者）
  )
)

先跑一下单轴代码，你就能看到你遇到的问题：附肢C的箱线图把另外两个的箱线挤得几乎看不见

ggplot(df, aes(x = Species, y = SA_V, fill = Appendage)) +
  geom_boxplot(position = position_dodge(0.8)) +
  labs(title = "单Y轴SA:V箱线图", x = "物种", y = "表面积体积比") +
  theme_minimal()

核心思路是：把高量级的附肢数据按比例缩小，适配低量级的轴范围，再通过第二个Y轴把数值还原回去。具体代码如下：

# 1. 定义转换系数：把附肢C的SA:V缩小到和其他附肢同量级（这里选8，你可以根据自己的数据调整）
scale_factor <- 8

# 2. 绘制双Y轴箱线图
ggplot() +
  # 第一个图层：绘制低量级的附肢A和B，用原始Y轴
  geom_boxplot(
    data = subset(df, Appendage != "附肢C"),
    aes(x = Species, y = SA_V, fill = Appendage),
    position = position_dodge(0.8), # 保证同一物种的箱线对齐
    width = 0.7
  ) +
  # 第二个图层：绘制高量级的附肢C，用转换后的Y值
  geom_boxplot(
    data = subset(df, Appendage == "附肢C"),
    aes(x = Species, y = SA_V / scale_factor, fill = Appendage),
    position = position_dodge(0.8),
    width = 0.7
  ) +
  # 3. 设置双Y轴：左侧是原始低量级刻度，右侧是还原后的高量级刻度
  scale_y_continuous(
    name = "SA:V（附肢A/B）",
    sec.axis = sec_axis(
      ~ . * scale_factor, # 把缩小的数值还原回去
      name = "SA:V（附肢C）"
    )
  ) +
  # 4. 自定义颜色，让不同附肢更易区分
  scale_fill_manual(
    values = c("附肢A" = "#1f77b4", "附肢B" = "#ff7f0e", "附肢C" = "#2ca02c")
  ) +
  labs(title = "双Y轴SA:V箱线图", x = "物种") +
  theme_minimal() +
  # 5. 优化视觉：把右侧Y轴的文字改成和附肢C一样的颜色，方便对应
  theme(
    legend.position = "top",
    axis.title.y.right = element_text(color = "#2ca02c"),
    axis.text.y.right = element_text(color = "#2ca02c")
  )

其实ggplot2官方并不太推荐双Y轴（容易误导读者对数据关系的判断），如果可以的话，分面展示是更清晰的选择：

ggplot(df, aes(x = Species, y = SA_V, fill = Appendage)) +
  geom_boxplot() +
  facet_wrap(~ Appendage, scales = "free_y") # 每个子图用独立的Y轴
  labs(title = "分面展示SA:V箱线图", x = "物种", y = "表面积体积比") +
  theme_minimal()

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