ggplot2森林图：x轴限制+误差棒箭头与log2转换问题

问题解决与修正代码 #

核心问题分析

1. x轴刻度异常：xlim(0,16)会强制覆盖scale_x_continuous的对数转换，导致刻度回到线性刻度。同时log2(0)无意义，不能直接将下限设为0。
2. 箭头错误与误差棒截断：现有箭头从OR值开始绘制，而非截断后的误差棒末端；未正确处理超出16的误差棒截断逻辑。

完整修正代码

library(tidyverse, scales)

data <- tibble(
  SNP = c('Univariable CM','1','2','3','4','5','6','7','8','9'),
  OR = c(10.13, 5.16, 2.83, 5.92, 4.71, 6.22, 3.79, 6.64, 7.80, 7.88),
  `lower OR` = c(5.3, 1.8, 1.6, 2.8, 2.3, 2.7, 1.3, 3.5, 3.8, 3.9),
  `upper OR` = c(19,14,5,12,7,14,11,13,18,17)
) %>%
  # 提前处理因子顺序，避免重复代码
  mutate(
    SNP = fct_relevel(SNP, "Univariable CM") %>% fct_rev(),
    # 截断上限到16，标记需要箭头的行
    trunc_upper = pmin(`upper OR`, 16),
    need_arrow = `upper OR` > 16
  )

data %>%
  ggplot(aes(y = SNP)) +
  theme_classic() +
  # 绘制点估计
  geom_point(aes(x = OR), shape = 15, size = 3) +
  # 绘制截断后的误差棒
  geom_linerange(aes(xmin = `lower OR`, xmax = trunc_upper)) +
  # 绘制超出范围的箭头（需关闭坐标裁剪）
  geom_segment(
    data = filter(data, need_arrow),
    aes(x = 16, xend = 16 * 1.1, y = SNP, yend = SNP),
    arrow = arrow(length = unit(0.3, "cm"))
  ) +
  # OR=1的参考线
  geom_vline(xintercept = 1, linetype = "dashed") +
  # 正确设置log2转换的x轴
  scale_x_continuous(
    trans = 'log2',
    breaks = c(1,2,4,8,16),
    labels = trans_format('log2', math_format(2^.x)),
    limits = c(0.5, 16), # 下限用0.5而非0，避免log2(0)报错
    expand = c(0,0)
  ) +
  # 关闭坐标裁剪，让箭头显示在x轴范围外
  coord_cartesian(clip = "off") +
  labs(x = "OR", y = "") +
  theme(
    axis.line.y = element_blank(),
    axis.ticks.y = element_blank(),
    axis.text.y = element_blank(),
    axis.title.y = element_blank(),
    legend.position = 'none',
    # 右侧留空容纳箭头
    plot.margin = margin(5.5, 20, 5.5, 5.5)
  )

关键改动说明

• 因子顺序处理：提前在数据中重新排序SNP因子，减少绘图代码重复。
• 误差棒截断：用pmin()将上限截断到16，只保留可见部分的误差棒。
• 箭头绘制：筛选出需要箭头的行，从16的位置向外延伸绘制箭头，并通过coord_cartesian(clip="off")让箭头显示在x轴范围外。
• x轴修正：在scale_x_continuous内设置limits替代xlim()，使用0.5作为下限确保log2转换生效，同时保留期望的刻度值。

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