[R] forest + igraph ?

Viechtbauer Wolfgang (STAT) wolfgang.viechtbauer at maastrichtuniversity.nl
Tue Apr 12 12:24:02 CEST 2011


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



More information about the R-help mailing list