如何用ggplot绘制叠加90%与95%置信区间的系数图？

第一步：计算90%置信区间 #

置信区间的核心公式是 系数 ± 临界z值×标准误，90%双侧置信区间对应的z临界值是qnorm(0.95)（双侧各留5%的尾部，取标准正态分布的95%分位数）。直接在现有数据框中新增90% CI的上下限列：

方法1：用dplyr（推荐，代码更简洁）

library(dplyr)

z_90 <- qnorm(0.95)
test <- test %>%
  mutate(ci_lower_90 = coef - z_90 * stderr,
         ci_upper_90 = coef + z_90 * stderr)

方法2：基础R（无需额外包）

z_90 <- qnorm(0.95)
test$ci_lower_90 <- test$coef - z_90 * test$stderr
test$ci_upper_90 <- test$coef + z_90 * test$stderr

第二步：叠加双CI的可视化 #

用ggplot分层绘制：先画半透明的95% CI，再画更醒目的90% CI，最后保留系数点，确保90% CI在视觉上更突出：

library(ggplot2)

ggplot(test, aes(x = coef, y = var)) +
  # 95% CI：半透明灰色，作为背景层
  geom_errorbarh(aes(xmin = ci_lower, xmax = ci_upper),
                 height = 0.2,  # 调整竖线高度，避免不同变量的误差线重叠
                 size = 1.5,
                 alpha = 0.3,
                 color = "gray50") +
  # 90% CI：黑色实线，叠加在上方
  geom_errorbarh(aes(xmin = ci_lower_90, xmax = ci_upper_90),
                 height = 0.2,
                 size = 1,
                 color = "black") +
  # 系数点：黑色实心点，明确标注系数位置
  geom_point(size = 3, color = "black") +
  theme_minimal()

关键细节说明 #

• 图层顺序决定视觉优先级：先画95% CI（底层），再画90% CI（上层），避免90% CI被遮挡
• 用alpha参数降低95% CI的透明度，既保留完整信息，又不让它抢90% CI的焦点
• 调整height参数可以让不同变量的误差线竖线更清晰，提升可读性

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