如何用ggplot2绘制类似Stata coefplot的平滑置信区间误差条

geom_errorbar本身只支持硬边界、统一粗细的误差条绘制，无法实现平滑渐变/两端收窄的效果，可采用以下三类方案实现需求：

方案1：使用ggdist包的geom_pointinterval（最便捷） #

该几何对象原生支持自定义区间的粗细变化、透明度渐变，完美匹配coefplot的平滑置信区间效果：

library(ggplot2)
library(ggdist)
library(dplyr)

# 示例数据结构：回归结果表res_df包含变量名term、系数estimate、95%置信区间上下限conf.low/conf.high
res_df <- tibble(
  term = c("变量1", "变量2", "变量3", "控制变量1", "控制变量2"),
  estimate = c(0.2, 0.5, -0.3, 0.1, -0.15),
  conf.low = c(0.05, 0.3, -0.45, -0.02, -0.28),
  conf.high = c(0.35, 0.7, -0.15, 0.22, -0.02)
)

ggplot(res_df, aes(y = term, x = estimate, xmin = conf.low, xmax = conf.high)) +
  # interval_size_range控制区间从中心到两端的粗细变化，.width指定置信区间水平
  geom_pointinterval(point_size = 3, interval_size_range = c(0.2, 1.8), .width = 0.95) +
  labs(x = "回归系数", y = "变量名") +
  theme_bw()

如果需要渐变透明效果，可以额外添加interval_alpha参数调整。

方案2：原生ggplot2手动实现（无需额外扩展包） #

通过给置信区间分段、给不同位置的线段设置不同粗细/透明度，模拟平滑效果：

# 为每个变量的置信区间生成100个分段点，计算每个点的透明度/粗细值
res_df_expand <- res_df %>%
  group_by(term) %>%
  reframe(
    estimate = unique(estimate),
    x_pos = seq(conf.low, conf.high, length.out = 100),
    # 距离系数中心越远，透明度越高、线条越细
    alpha_val = 1 - abs(x_pos - estimate)/max(abs(x_pos - estimate)),
    size_val = 2 * (1 - abs(x_pos - estimate)/max(abs(x_pos - estimate)))
  )

ggplot(res_df_expand, aes(y = term)) +
  geom_segment(aes(x = x_pos, xend = x_pos, 
                   y = as.numeric(term) - 0.1, yend = as.numeric(term) + 0.1,
                   alpha = alpha_val, linewidth = size_val),
               show.legend = FALSE) +
  geom_point(aes(x = estimate), size = 3) +
  labs(x = "回归系数", y = "变量名") +
  theme_bw()

方案3：使用回归可视化封装包 #

modelsummary、ggeffects这类专门处理回归结果的包，内置了平滑置信区间的绘图选项，传入拟合好的回归模型即可直接出图，无需手动整理置信区间数据。

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