R中的并排森林地块

use*_*672 4 plot r graph

所以这是一些通用数据.我目前使用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)
Run Code Online (Sandbox Code Playgroud)

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

有谁知道如何在r中制作一个并排的森林情节?我想马上两个森林的地块对左侧的子组的名称,然后在右边(对方旁边)有点像这个http://www.nature.com/ng/journal/v43/n8/images_article/ng .882-F1.jpg 但更简单.

我已经尝试了par()metafor包的功能,但它不起作用.

小智 6

这让你非常接近:

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)
Run Code Online (Sandbox Code Playgroud)

并排的森林地块