我可以在R metafor包的forest()中对齐/证明置信区间

时间:2016-06-17 20:19:29

标签: r

我有13项研究,我正在荟萃分析其中的平均差异。有没有办法证明或对齐置信区间,以便括号全部排列在森林图中?这是我当前的代码和生成的森林图的图片。

代码:

forest(result.mbw,  xlab=NA, 
       alim=c(-300,300), xlim=c(-700,500), digits=0)
text(c(-600),14, c("Study"))
text(c(200),14, c("Mean BW Difference (g) [95% CI]")) 

Forest plot

1 个答案:

答案 0 :(得分:1)

你必须使用固定宽度的字体才能实现。因此,您可以在创建绘图之前使用par(family="mono")。由于文本宽度可能更宽,因此您可能必须使用xlim函数的forest()参数来使其适合。这是一个例子:

library(metafor)
dat <- escalc(measure="MD", m1i=m1i, sd1i=sd1i, n1i=n1i, m2i=m2i, sd2i=sd2i, n2i=n2i, data=dat.normand1999)
res <- rma(yi, vi, data=dat)
par(family="mono")
forest(res, xlim=c(-180,230), cex=.8)
text(-180, 11, "Study", pos=4, font=2, cex=.8)
text( 230, 11, "Mean Difference [95% CI]", pos=2, font=2, cex=.8)

我不得不减少一点cex以使其更合适。此外,我在添加文字时使用pos=4pos=2,以便它也能正确对齐。

如果您使用的是Windows并希望使用其他字体,则可以执行以下操作:

windowsFonts(Consolas=windowsFont("Consolas"))
par(family="Consolas")

或者使用您喜欢的任何其他固定宽度字体。不确定这在MacOS或Unix / Linux下是如何工作的。