R中的并排森林地块

Question

所以这是一些通用数据。 我目前使用metafor包和par()函数来并排绘制这些图。 但是，我想理想地摆脱第二个图中的研究列和两个图中的日志RR估计。 我希望两个森林地块也能更紧密地联系在一起。 很抱歉以前没有提供数据 - 我是新来的。

library(metafor)
par(mfrow=c(1,2))

### load BCG vaccine data
data(dat.bcg)
### calculate log relative risks and corresponding sampling variances
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
### default forest plot of the observed log relative risks
forest(dat$yi, dat$vi)
### default forest plot of the observed log relative risks for second plot
forest(dat$yi, dat$vi)

===============================

有谁知道如何在r中制作一个并排的森林情节？ 我希望左侧有子组名，然后是右侧的两个森林图（彼此相邻）有点像这样http://www.nature.com/ng/journal/v43/n8/images_article/ng .882-F1.jpg但更简单。

我已尝试使用metafor包的par()函数，但它不起作用。

Answer 1

这让你非常接近：

library(metafor)
data(dat.bcg)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)

par(mfrow=c(1,2))
par(mar=c(5,4,1,1))
forest(dat$yi, dat$vi, annotate=FALSE, cex=.8, at=seq(-3,2,1), digits=1, xlim=c(-6,2))
text(0, 15, "Figure 1", cex=.8, font=2)
par(mar=c(5,3,1,2))
forest(dat$yi, dat$vi, annotate=FALSE, slab=rep("",length(dat$yi)), cex=.8, at=seq(-3,2,1), digits=1, xlim=c(-5,3))
text(0, 15, "Figure 2", cex=.8, font=2)