如何在ggplot2多分面图中跨分面添加显著性标记？

问题描述 #

使用ggplot2的facet_wrap创建多分面图表时，无法通过geom_signif在不同分面的对应类别（如示例中分面C=1和C=2的A=1）之间添加跨分面的显著性连线与标注。由于分面的x轴独立，无法设置正确的xmax值连接相邻分面，调整分组变量也无效，且不能调换fill与分面变量。期望实现：绘制连接分面1与分面2中A=1的连线，并在连线中间标注*。

解决方案 #

由于facet_wrap的分面是独立坐标系，geom_signif无法直接跨分面绘制元素。可以通过操作ggplot的底层gtable对象精确计算跨分面坐标，手动添加连线和标注，具体步骤如下：

1. 准备数据与基础绘图 #

首先生成可重复的示例数据，并创建基础分面图：

set.seed(123) # 保证结果可重复
library(ggplot2)
library(dplyr)
library(gtable)
library(grid)

value <- sample(1:50, 100, replace = T)
B <- as.factor(rep(c(1,2,2,1), 50))
C <- as.factor(rep(c(1,2), each=100))
A <- as.factor(rep(c(4,4,3,3,2,2,1,1), 25))

dt1 <- data.frame(A,B,C,value) %>%
  group_by(A,B,C)

# 创建基础分面图
p <- ggplot(dt1, aes(x=A, y=value, fill=B)) +
  theme_classic() +
  facet_wrap(vars(C)) +
  geom_boxplot(position='dodge', outlier.size=0) +
  geom_point(size=0.5, position=position_jitterdodge(0.2))

2. 转换为gtable并计算跨分面坐标 #

将ggplot对象转换为gtable，获取两个分面的位置信息，计算目标类别（A=1）在两个分面中的全局坐标：

# 构建绘图对象获取坐标参数
pb <- ggplot_build(p)
# 转换为gtable对象
gt <- ggplotGrob(p)

# 定位两个分面的布局位置
panel1 <- gt$layout[gt$layout$name == "panel-1-1", ]
panel2 <- gt$layout[gt$layout$name == "panel-1-2", ]

# 目标类别A=1在分面内的x坐标（因子水平对应x=1）
x_panel <- 1
# 目标y位置（对应数据中的y值）
y_target <- 55

# 转换为全局坐标（基于gtable的网格单位）
x1_global <- panel1$l + (x_panel - pb$layout$panel_params[[1]]$x.range[1]) * 
  (panel1$r - panel1$l) / diff(pb$layout$panel_params[[1]]$x.range)
x2_global <- panel2$l + (x_panel - pb$layout$panel_params[[2]]$x.range[1]) * 
  (panel2$r - panel2$l) / diff(pb$layout$panel_params[[2]]$x.range)
y_global <- panel1$t + (y_target - pb$layout$panel_params[[1]]$y.range[1]) * 
  (panel1$b - panel1$t) / diff(pb$layout$panel_params[[1]]$y.range)

3. 添加跨分面的连线与标注 #

在gtable中添加连线、短竖线和显著性标注：

# 添加跨分面的水平连线
gt <- gtable_add_grob(gt, linesGrob(x = c(x1_global, x2_global), y = c(y_global, y_global)), 
                      t = panel1$t, b = panel1$b, l = panel1$l, r = panel2$r)

# 添加两端的短竖线（连接连线到箱线图顶部）
gt <- gtable_add_grob(gt, linesGrob(x = c(x1_global, x1_global), y = c(y_global - 0.5, y_global)), 
                      t = panel1$t, b = panel1$b, l = panel1$l, r = panel1$r)
gt <- gtable_add_grob(gt, linesGrob(x = c(x2_global, x2_global), y = c(y_global - 0.5, y_global)), 
                      t = panel2$t, b = panel2$b, l = panel2$l, r = panel2$r)

# 添加中间的显著性标注
gt <- gtable_add_grob(gt, textGrob("*", x = (x1_global + x2_global)/2, y = y_global + 0.5, 
                                   gp = gpar(size=12)), 
                      t = panel1$t, b = panel1$b, l = panel1$l, r = panel2$r)

# 绘制最终图表
grid.draw(gt)

说明 #

• 该方法通过直接操作gtable对象，绕过了分面的独立坐标系限制，能精确控制跨分面元素的位置。
• 若需要添加多个跨分面标记，只需重复计算对应类别的坐标并添加相应的grob元素即可。
• 可根据实际需求调整y_target（连线高度）、标注大小size等参数。

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