Re: [R] forest + igraph ?

From: Viechtbauer Wolfgang (STAT) <wolfgang.viechtbauer_at_maastrichtuniversity.nl>
Date: Tue, 12 Apr 2011 12:24:02 +0200

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
______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Received on Tue 12 Apr 2011 - 10:31:42 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 12 Apr 2011 - 15:40:29 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive