r - 如何将 tau^2 添加到 metafor R 包中的 Forest Plot？

Question

我正在使用 metafor 包在 R 中进行元分析。我按照http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups上的教程进行操作

如何将 tau^2 添加到“所有研究的 RE 模型”？

library(metafor)

### load BCG vaccine dataset
data(dat.bcg)

res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="OR",
       slab=paste(author, year, sep=", "), method="REML")

###Plot Forest Plot
forest(res)

### Add Text
text(-5, -1, pos=4, cex=0.75, bquote(paste("RE Model for All Studies (Q 
= ",.(formatC(res$QE, digits=2, format="f")), ", df = ", .(res$k - 
res$p),", p = ", .(formatC(res$QEp, digits=2, format="f")), "; ", I^2, " = ",.(formatC(res$I2, digits=1, format="f")), "%)")))

在此先感谢，C.

Answer 1

这将打印一个带有上标 2 的希腊 tau，后跟 res 对象中 tau2 项的两位小数：

text(-5, -0.5, pos=4, cex=0.75, bquote(
paste("RE Model for All Studies (Q = ",
.(formatC(res$QE, digits=2, format="f")), 
", df = ", .(res$k - res$p),", p = ", 
.(formatC(res$QEp, digits=2, format="f")),
 "; ", I^2, " = ",
.(formatC(res$I2, digits=1, format="f")), 
 "%)", "; ", tau^2 == 
.(formatC(res$tau2, digits=2, format="f")))))

我将文本向上移动了半行，这样它就不会覆盖摘要 OR-diamond。您应该了解您实际上是在使用plotmath粘贴功能，这意味着您可以更改tau为Tau如果您想要一个大写的希腊语 Tau。阅读?plotmath。