使用R语言meta包forest函数绘制OR森林图的置信区间异常问题

核心原因 #

你直接将原始比值比(OR)传入metagen的TE参数，但当sm="OR"时，TE要求传入的是log转换后的OR值（logOR）。函数默认会把你输入的TE当作logOR处理，后续的SE和置信区间计算全基于这个错误的前提，自然和你输入的lower/upper对不上。

解决方案1：修正输入格式（meta包） #

先把OR、置信上下限都转换为对数尺度，再传入metagen，函数会直接使用你提供的置信区间，不会重新计算SE：

library(meta)

# 原始数据
OR_values <- c(0.58, 0.04, 0.26, 0.92, 0.3)
lower_OR <- c(0.21, 0.02, 0.16, 0.64, 0.0)
upper_OR <- c(1.59, 0.08, 0.42, 1.32, 1.2)

# 转换为对数尺度（注意：log(0)无意义，把0替换为极小值0.001）
logTE <- log(OR_values)
loglower <- log(ifelse(lower_OR == 0, 0.001, lower_OR))
logupper <- log(upper_OR)

# 构建meta对象
m <- metagen(
  TE = logTE,
  lower = loglower,
  upper = logupper,
  overall = FALSE, random = FALSE, fixed = FALSE,
  sm = "OR"
)

# 绘制森林图
forest(m)

解决方案2：手动传入计算好的SE（meta包） #

按公式SE = (log(upper) - log(lower))/3.92计算对数尺度下的标准误，直接传给seTE参数，函数会跳过自动计算步骤：

# 计算对数尺度的SE
logSE <- (logupper - loglower)/3.92

# 传入TE和seTE
m <- metagen(
  TE = logTE,
  seTE = logSE,
  overall = FALSE, random = FALSE, fixed = FALSE,
  sm = "OR"
)

forest(m)

更适合非Meta场景的替代方案：forestplot包 #

如果只是做OR和置信区间的可视化，不需要Meta分析的统计计算，forestplot包更简洁直接，无需处理log转换的细节：

library(forestplot)

# 整理数据（替换0为极小值）
plot_data <- data.frame(
  研究 = paste0("研究", 1:5),
  OR = OR_values,
  下限 = ifelse(lower_OR == 0, 0.001, lower_OR),
  上限 = upper_OR
)

# 绘制森林图
forestplot(
  labeltext = plot_data$研究,
  mean = plot_data$OR,
  lower = plot_data$下限,
  upper = plot_data$上限,
  xlog = TRUE, # OR用对数轴更符合统计展示逻辑
  title = "比值比(OR)与95%置信区间",
  xticks = c(0.1, 0.5, 1, 2, 5)
)

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