如何在ggplot2分面图中添加多层X轴标签？

问题说明 #

要在ggplot2分面柱状图中实现多层X轴：

• xaxis1对应每个X轴刻度显示
• xaxis2仅在对应分组上方显示一次
  使用annotate添加文本时触发参数长度不匹配错误。

示例数据 #

df <- data.frame(START = rep(seq(1:10),2),
                 SAMPLE = rep(c("A", "B"), each=10),
                 yaxis = runif(20, min = 1, max = 10),
                 xaxis1 = rep(c("E1", "E1", "E2", "E2", "E1", "E2", "E3", "E4", "E1", "E2"),2),
                 xaxis2 = rep(c("G1", "G1", "G1", "G1", "G2", "G2", "G2", "G2", "G3", "G3"),2))

错误原因 #

annotate要求x、y、label三个参数长度完全一致。报错代码中第二个annotate的x是长度20的序列，但label用了unique(df$xaxis2)（长度3），两者长度不匹配；同时直接用START作为x轴，分面后每个面板x轴范围重复1-10，也会导致多层轴位置混乱。

解决方案：基于分组中心的多层轴实现 #

可以通过预处理数据计算分组中心，结合geom_text适配分面场景，代码如下：

library(ggplot2)
library(dplyr)

# 预处理：计算每个分面内xaxis2分组的x轴中心位置
df_summary <- df %>%
  group_by(SAMPLE, xaxis2) %>%
  summarise(x_center = mean(START), .groups = "drop")

ggplot(data = df, aes(x = START, y = yaxis, fill = yaxis)) +
  geom_col(width = 1, position = "identity") +
  # 添加xaxis1标签（对应每个刻度）
  geom_text(aes(label = xaxis1), y = -0.5, size = 4, vjust = 0) +
  # 添加xaxis2标签（每个分组仅显示一次）
  geom_text(data = df_summary, aes(x = x_center, label = xaxis2), 
            y = -1.5, size = 6, vjust = 0) +
  # 分面设置
  facet_wrap(SAMPLE ~ ., ncol = 1, strip.position = "left") +
  # 调整x轴刻度
  scale_x_continuous(breaks = 1:10) +
  # 允许标签超出绘图区域，避免截断
  coord_cartesian(clip = "off") +
  # 调整主题，隐藏原x轴元素
  theme(
    plot.margin = margin(b = 40, unit = "pt"),
    axis.title.x = element_blank(),
    axis.text.x = element_blank(),
    axis.ticks.x = element_blank()
  )

代码说明 #

1. 预处理数据时按分面变量SAMPLE和分组变量xaxis2聚合，计算每个分组的x轴中心位置，确保标签显示在分组正上方
2. 用geom_text替代annotate，自动适配分面逻辑，避免长度不匹配问题
3. 通过coord_cartesian(clip = "off")和plot.margin调整，防止底部的多层轴标签被截断
4. 隐藏原x轴的刻度和文本，改用自定义的多层标签实现需求

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