[R] Create patterns within a plot?

Viechtbauer Wolfgang (STAT) wolfgang.viechtbauer at maastrichtuniversity.nl
Thu Mar 14 14:26:37 CET 2013


Same idea, but just using base graphics and using an x-axis with log units:

library(metafor)

data(dat.bcg)

res <- rma(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, slab=paste(dat.bcg$author, dat.bcg$year, sep=", "))

forest(res, cex=.9, at=log(c(.05, .25, 1, 4)), xlim=c(-10,7), atransf=exp)
rect(log(0.5), 0.1, log(1.5), 13.9, col="gray90", border=NA)
par(new=TRUE)
forest(res, cex=.9, at=log(c(.05, .25, 1, 4)), xlim=c(-10,7), atransf=exp)

Best,
Wolfgang

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
> On Behalf Of Jim Lemon
> Sent: Thursday, March 14, 2013 12:19
> To: Öhagen Patrik
> Cc: R hELP
> Subject: Re: [R] Create patterns within a plot?
> 
> On 03/14/2013 09:19 PM, Öhagen Patrik wrote:
> >
> > Please note that I want to use log-scale on the x-axis which makes the
> polygon function sad.
> >
> > Cheers, Patrik
> >
> >
> > -----Ursprungligt meddelande-----
> > Från: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
> För Öhagen Patrik
> > Skickat: den 14 mars 2013 08:47
> > Till: R hELP
> > Ämne: [R] Create patterns within a plot?
> >
> >
> >
> > Dear All,
> >
> > As an attempt to highlight the overall pattern in a Forest plot, I would
> like to "highlight" the area around HR=1. I cannot find any simple tools
> for painting a grey ribbon between 0.9 and 1.1. Any suggestions?
> >
> Hi Patrik,
> This is a bit indirect, but it may do what you want. The example plot is
> from the metafor package (see forest.default):
> 
> forest(dat$yi, dat$vi)
> library(plotrix)
> gradient.rect(0.9,-0.1,1.1,14,
>   reds=c(1,0.5,1),greens=c(1,0.5,1),
>   blues=c(1,0.5,1),gradient="x",border=NA)
> par(new=TRUE)
> forest(dat$yi, dat$vi)
> 
> Jim
> 
> ______________________________________________
> 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.



More information about the R-help mailing list