如何在ggplot的facet_wrap中让所有分面以0为中心对齐？

问题 #

使用facet_wrap()垂直堆叠3个图表，要求所有图表中yintercept=0的虚线完全上下对齐，同时保留x轴（经coord_flip()后对应原y轴）的自由缩放能力（已通过scales="free"实现）。

解决方法 #

方法1：用ggh4x的增强分面快速实现 #

ggh4x包的facet_wrap2()提供了align="all"参数，能在保留轴自由缩放的同时，对齐所有分面的基准线（这里是0点），是最简便的方案。

首先安装并加载包：

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

修改后的绘图代码：

d_example %>% 
  ggplot(aes(x = var, y = coef, 
             ymin = ci_lower, ymax = ci_upper,
             color = race)) +
  geom_point(position = position_dodge(width = 0.7), size = 3) +
  geom_errorbar(width = 0, size = 1.1, alpha = 0.6, 
                position = position_dodge(width = 0.7)) +
  facet_wrap2(~ outcome, nrow = 3, scales = "free", align = "all") +  # 替换为facet_wrap2并设置对齐规则
  coord_flip() +
  geom_hline(yintercept = 0, linetype = "dashed", color = "black", size = .3) +
  theme_minimal()

方法2：手动计算对称轴范围（可避免依赖第三方包） #

如果不想额外安装包，可以手动计算每个分面的轴范围，让每个分面的y轴（原y轴，翻转后对应x轴）以0为中心，即范围设置为[-max_abs, max_abs]，其中max_abs是该分面所有y值（包括置信区间）的最大绝对值。

先计算每个分面的对称范围：

library(dplyr)

axis_ranges <- d_example %>%
  group_by(outcome) %>%
  summarise(
    max_abs = max(abs(c(coef, ci_lower, ci_upper))),
    .groups = "drop"
  ) %>%
  mutate(xmin = -max_abs, xmax = max_abs)

再结合ggh4x的facetted_pos_scales动态设置每个分面的轴范围：

library(ggh4x)
library(purrr)

# 生成每个分面的y轴缩放规则
scale_list <- map2(axis_ranges$xmin, axis_ranges$xmax, ~scale_y_continuous(limits = c(.x, .y)))
names(scale_list) <- axis_ranges$outcome

# 绘图
d_example %>% 
  ggplot(aes(x = var, y = coef, 
             ymin = ci_lower, ymax = ci_upper,
             color = race)) +
  geom_point(position = position_dodge(width = 0.7), size = 3) +
  geom_errorbar(width = 0, size = 1.1, alpha = 0.6, 
                position = position_dodge(width = 0.7)) +
  facet_wrap(~ outcome, nrow = 3, scales = "free") +
  coord_flip() +
  geom_hline(yintercept = 0, linetype = "dashed", color = "black", size = .3) +
  facetted_pos_scales(y = scale_list) +
  theme_minimal()

示例数据 #

d_example <- structure(list(var = c("score_A", "score_B", "score_C", "score_A", 
"score_B", "score_C", "score_A", "score_B", "score_C", "friends_A", 
"friends_B", "friends_A", "friends_B", "friends_A", "friends_B", 
"poverty_A", "poverty_B", "poverty_A", "poverty_B", "poverty_A", 
"poverty_B"), coef = c(-0.8, -0.7, 0.12, -0.7, -0.9, -0.05, -0.73, 
-1, 0.02, 0.55, -0.09, 0.51, 0.2, 0.65, 0.05, 1.8, 0.34, 1.103, 
0.077, 1.02, 0.06), race = c("white", "white", "white", "black", 
"black", "black", "hispanic", "hispanic", "hispanic", "white", 
"white", "black", "black", "hispanic", "hispanic", "white", "white", 
"black", "black", "hispanic", "hispanic"), outcome = c("score", 
"score", "score", "score", "score", "score", "score", "score", 
"score", "friends", "friends", "friends", "friends", "friends", 
"friends", "ses", "ses", "ses", "ses", "ses", "ses"), ci_lower = c(-1, 
-0.9, -0.08, -0.9, -1.1, -0.25, -0.93, -1.2, -0.18, 0.35, -0.29, 
0.31, 0, 0.45, -0.15, 1.6, 0.14, 0.903, -0.123, 0.82, -0.14), 
    ci_upper = c(-0.6, -0.5, 0.32, -0.5, -0.7, 0.15, -0.53, -0.8, 
    0.22, 0.75, 0.11, 0.71, 0.4, 0.85, 0.25, 2, 0.54, 1.303, 
    0.277, 1.22, 0.26)), class = c("spec_tbl_df", "tbl_df", "tbl", 
"data.frame"), row.names = c(NA, -21L), spec = structure(list(
    cols = list(var = structure(list(), class = c("collector_character", 
    "collector")), coef = structure(list(), class = c("collector_double", 
    "collector")), race = structure(list(), class = c("collector_character", 
    "collector")), outcome = structure(list(), class = c("collector_character", 
    "collector")), ci_lower = structure(list(), class = c("collector_double", 
    "collector")), ci_upper = structure(list(), class = c("collector_double", 
    "collector"))), default = structure(list(), class = c("collector_guess", 
    "collector")), skip = 1L), class = "col_spec"))

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