R语言森林图添加作者年份上标报错问题求助

核心问题 #

你尝试用expression()生成带上下标的文本时报错，是因为forestplot不直接支持expression向量作为标签输入，得换用它能识别的格式：HTML文本或者grid包的textGrob对象。

方法1：HTML格式文本（简单高效） #

如果你的forestplot是较新版本，直接拼接带<sup>标签的字符串即可：

library(forestplot)

# 拼接带HTML上标的标签列
m$label_sup <- paste0(m$fulltitle, "<sup>", m$refnum, "</sup>")

# 绘制森林图，开启HTML解析
forestplot(labeltext = m$label_sup,
           mean = m$mean,
           lower = m$lower,
           upper = m$upper,
           # 确保文本样式正常
           txt_gp = fpTxtGp(label = gpar(fontsize = 10)),
           # 标记所有行不是汇总行（如果有汇总行需单独调整）
           is.summary = rep(FALSE, nrow(m)))

方法2：textGrob通用方案（兼容旧版本） #

如果你的forestplot版本不支持HTML，用grid包生成带解析表达式的文本对象：

library(forestplot)
library(grid)

# 生成带上下标的textGrob列表
label_grobs <- lapply(seq(nrow(m)), function(i) {
  textGrob(
    label = paste0(m$fulltitle[i], "^", m$refnum[i]),
    gp = gpar(fontsize = 10),
    just = "left",
    # 让R解析^作为上标符号
    parse = TRUE
  )
})

# 传入grob列表绘制森林图
forestplot(labeltext = label_grobs,
           mean = m$mean,
           lower = m$lower,
           upper = m$upper)

模拟数据测试 #

如果你的数据结构如下，上面两种方法都能正常生成带Smith 2020¹格式的上标：

# 模拟你的数据框m
m <- data.frame(
  fulltitle = c("Smith 2020", "Jones 2021", "Lee 2022"),
  refnum = c(1, 2, 3),
  mean = c(0.5, 0.7, 0.6),
  lower = c(0.3, 0.5, 0.4),
  upper = c(0.7, 0.9, 0.8)
)

为什么expression会报错？ #

forestplot的labeltext只接受字符向量、矩阵或grob对象列表，expression属于R的表达式对象，不在支持范围内，所以直接传入会触发错误。而上面两种方法都是把文本转换成了forestplot能识别的格式。

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