[R] forest + igraph ?

[Viechtbauer Wolfgang (STAT)]

Sure. Here is an example:

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)
res <- rma(yi, vi, data=dat)
forest(res, atransf=exp)
title("Forest Plot of Relative Risks")
dat <- escalc(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
res <- rma(yi, vi, data=dat)
forest(res, atransf=exp)
title("Forest Plot of Odds Ratios")

Best,

-- 
Wolfgang Viechtbauer 
Department of Psychiatry and Neuropsychology 
School for Mental Health and Neuroscience 
Maastricht University, P.O. Box 616 
6200 MD Maastricht, The Netherlands 
Tel: +31 (43) 368-5248 
Fax: +31 (43) 368-8689 
Web: http://www.wvbauer.com 

> -----Original Message-----
> From: Samor Gandhi [mailto:samorgandhi at yahoo.com]
> Sent: Tuesday, April 12, 2011 15:59
> To: r-help at stat.math.ethz.ch; Viechtbauer Wolfgang (STAT)
> Subject: RE: [R] forest + igraph ?
> 
> Thanks, That is a nice one. Is there any option that I can plot the pooled
> estimate?
> 
> Regards,
> Samor
> 
> --- On Tue, 12/4/11, Viechtbauer Wolfgang (STAT)
> <wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
> 
> From: Viechtbauer Wolfgang (STAT)
> <wolfgang.viechtbauer at maastrichtuniversity.nl>
> Subject: RE: [R] forest + igraph ?
> To: "Samor Gandhi" <samorgandhi at yahoo.com>, "r-help at stat.math.ethz.ch" <r-
> help at stat.math.ethz.ch>
> Date: Tuesday, 12 April, 2011, 15:54
> 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,
> 
> --
> Wolfgang Viechtbauer
> Department of Psychiatry and Neuropsychology
> School for Mental Health and Neuroscience
> Maastricht University, P.O. Box 616
> 6200 MD Maastricht, The Netherlands
> Tel: +31 (43) 368-5248
> Fax: +31 (43) 368-8689
> Web: http://www.wvbauer.com
> 
> 
> > -----Original Message-----
> > From: Samor Gandhi [mailto:samorgandhi at yahoo.com]
> > Sent: Tuesday, April 12, 2011 11:52
> > To: r-help at stat.math.ethz.ch; Viechtbauer Wolfgang (STAT)
> > Subject: RE: [R] forest + igraph ?
> >
> > Thank you for your reply. I would like to have two forest plots one on
> top
> > and the other on the bottom. I am using R version 2.12.2 (32-bit)
> Windows.
> > The code you sent me still plotting two windows one after the other?
> >
> > Best wishes and many thanks,
> > Samor
> >
> > --- On Tue, 12/4/11, Viechtbauer Wolfgang (STAT)
> > <wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
> >
> > From: Viechtbauer Wolfgang (STAT)
> > <wolfgang.viechtbauer at maastrichtuniversity.nl>
> > Subject: RE: [R] forest + igraph ?
> > To: "r-help at stat.math.ethz.ch" <r-help at stat.math.ethz.ch>
> > Cc: "Samor Gandhi" <samorgandhi at yahoo.com>
> > Date: Tuesday, 12 April, 2011, 13:42
> > It is not clear (at least to me) what exactly you want. You want two
> > forest plots "in one graph" but apparently not side-by-side or one on
> top
> > and the other on the bottom. So, you want to superimpose them? How do
> you
> > want to do that without creating an illegible mess? Or do you want one
> > graph, where you can scroll through various plots? Then try this:
> >
> > library(metafor)
> > data(dat.bcg)
> >
> > if (interactive()) {
> >    windows(record=TRUE)
> >    dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg,
> > data=dat.bcg)
> >    forest(dat$yi, dat$vi)
> >    dat <- escalc(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg,
> > data=dat.bcg)
> >    forest(dat$yi, dat$vi)
> > }
> >
> > and then use PageUp and PageDown to switch between the figures.
> >
> > (I don't know what OS you are using, so windows() may not work).
> >
> > Best,
> >
> > --
> > Wolfgang Viechtbauer
> > Department of Psychiatry and Neuropsychology
> > School for Mental Health and Neuroscience
> > Maastricht University, P.O. Box 616
> > 6200 MD Maastricht, The Netherlands
> > Tel: +31 (43) 368-5248
> > Fax: +31 (43) 368-8689
> > Web: http://www.wvbauer.com
> >
> > > -----Original Message-----
> > > From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org]
> > > On Behalf Of Samor Gandhi
> > > Sent: Monday, April 11, 2011 18:25
> > > To: r-help at stat.math.ethz.ch
> > > Subject: [R] forest + igraph ?
> > >
> > > Hello,
> > >
> > > Is it possible to have two meta-plots in one graph (not
> > par(mfrow=c(2,1))?
> > > But somthing like
> > >
> > >  library(metafor)
> > >  library("igraph")
> > >
> > >  if (interactive()) {
> > >     forest(dat.Treat$RR, ci.lb=dat.Treat$lower, ci.ub=dat.Treat$upper,
> > > xlab="Relative Risk",slab=dat.Treat$ID,refline=1)
> > >     forest(dat.Control$RR, ci.lb=dat.Control$lower,
> > > ci.ub=dat.Control$upper, xlab="Relative
> > > Risk",slab=dat.Control$ID,refline=1)
> > >   }
> > >
> > > i.e. both metaplots on the same graph!
> > >
> > > Regards,
> > > Samor