如何使用ggplot为过长误差棒添加箭头截断并正确定位？

要实现仅给指定绿色误差棒末端添加截断箭头，可按以下步骤调整代码：

核心思路 #

1. 先设置y轴显示上限，截断超出上限的误差棒数值，避免长条遮挡其余结果
2. 筛选出需要添加箭头的目标观测，调用与误差棒完全一致的position_dodge参数保证位置对齐
3. 适配翻转后的坐标系调整箭头方向

完整可运行代码 #

library(ggplot2)
library(scales)

# 加载示例数据
dat <- 
structure(list(outcome = c("cancer", "cancer", "cancer", "cancer", 
"cancer", "cancer", "cancer", "cancer"), exposure = c("insulin", 
"insulin", "insulin", "insulin", "LDL-cholesterol", "LDL-cholesterol", 
"LDL-cholesterol", "LDL-cholesterol"), OR = c(3.925713203, 3.49400025, 
3.545796066, 8.283204105, 0.910187777, 0.902219745, 0.918412256, 
0.953392987), LCI = c(2.285899133, 1.609867994, 0.879681108, 
0.672007727, 0.817076527, 0.789469504, 0.82104523, 0.870226921
), UCI = c(6.741865347, 7.583253904, 14.29230392, 14.292304, 
1.013909666, 1.03107272, 1.027325952, 1.044507088), method = c("Inverse variance weighted", 
"Weighted median", "Weighted mode", "MR Egger", "Weighted mode", 
"MR Egger", "Weighted median", "Inverse variance weighted")), row.names = c(73L, 
74L, 75L, 76L, 139L, 140L, 141L, 142L), class = "data.frame")

# 自定义y轴显示上限，可根据需求调整
y_max <- 10
# 截断超出上限的误差棒上限值
dat$UCI_trim <- ifelse(dat$UCI > y_max, y_max, dat$UCI)
# 筛选需要加箭头的观测：这里取UCI最大的绿色误差棒记录，也可按method/exposure手动筛选
arrow_data <- dat[which.max(dat$UCI), ]

# 绘图
ggplot(dat, aes(x=exposure, y=OR, ymin=LCI, ymax=UCI_trim, col=method)) + 
  geom_linerange(position=position_dodge(width = 0.5)) +
  geom_hline(yintercept=1, lty=2) +
  geom_point(stroke = 0.5, position=position_dodge(width = 0.5)) +
  scale_x_discrete(name="暴露因素") +
  # oob_squish参数保证超出范围的数值不会被判定为缺失
  scale_y_continuous(name="比值比 (95% 置信区间)", limits = c(0, y_max), oob = oob_squish) +
  coord_flip() +
  scale_colour_manual(values=c("#D55E00","#009E73", "#56B4E9","#E69F00")) + 
  theme_bw() + 
  theme(panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(),
        text = element_text(size=11), 
        legend.position="none") + 
  # 添加箭头，和误差棒使用相同的dodge参数保证对齐
  geom_segment(data = arrow_data,
               aes(x = exposure, y = y_max - 0.3, xend = exposure, yend = y_max), 
               position = position_dodge(width = 0.5),
               size = 0.25,
               arrow = arrow(length = unit(0.13, "cm"), angle = 15))

调整说明 #

• 如果需要箭头旋转180度，仅需交换geom_segment中的y和yend参数值即可
• 可调整y_max修改y轴显示范围，调整y_max - 0.3的差值修改箭头长度
• 如果有多条需要截断的误差棒，仅需调整arrow_data的筛选逻辑批量生成箭头数据即可

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