Skip to content
Prev 256559 / 398502 Next

forest + igraph ?

You said that you do NOT want to use par(mfrow=c(2,1)). Why not?

Isn't this (below) what you want?

library(metafor)
data(dat.bcg)

windows(height=8, width=6, pointsize=10)
par(mfrow=c(2,1))
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
forest(dat$yi, dat$vi, atransf=exp)
title("Forest Plot of Relative Risks")
dat <- escalc(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
forest(dat$yi, dat$vi, atransf=exp)
title("Forest Plot of Odds Ratios")

Best,