如何将 tau^2 添加到 metafor R 包中的森林图中?
How can I add tau^2 to Forest Plot in metafor R-package?
我正在使用 metafor 包在 R 中进行元分析。我遵循了 http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups
上的教程
如何将 tau^2 添加到 "RE Model for All studies"?
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.
这将打印一个带有上标 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
.
我正在使用 metafor 包在 R 中进行元分析。我遵循了 http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups
上的教程如何将 tau^2 添加到 "RE Model for All studies"?
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.
这将打印一个带有上标 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
.