Re: [R] metaplot

From: cheba meier <cheba.meier_at_googlemail.com>
Date: Fri, 08 Apr 2011 10:31:31 +0100

Thank you for your reply.

> forest(log(rr), ci.lb=lci, ci.ub=uci, xlab="Log Relative Risk")
Error in sei^2 : 'sei' is missing

Regards,
Cheba

2011/4/7 Viechtbauer Wolfgang (STAT) <
wolfgang.viechtbauer_at_maastrichtuniversity.nl>

> Actually, Scott suggested ggplot2, which is available on CRAN:
>
> http://cran.r-project.org/web/packages/ggplot2/index.html
>
> Moreover, for the metafor package, you do not need "(a,b,c,d) or standard
> errors for the effect estimates".
>
> rr <- c(0.61, 0.35, 1.16, 0.94, 1.16, 0.46, 0.99, 1.05, 1.15, 0.27)
> lci <- c(-1.69, -1.27, -0.24, -0.54, -1.01, -1.8,-0.89, -1.38, -1.05,-2.04)
> uci <- c(0.70, -0.82,0.54,0.41,1.30,0.24,0.87,1.48,1.33,-0.57)
>
> The lci and uci values are obviously the log-transformed values of the
> upper and lower CI bounds for the RR. So, you could do:
>
> library(metafor)
> forest(log(rr), ci.lb=lci, ci.ub=uci, xlab="Log Relative Risk")
>
> or:
>
> forest(log(rr), ci.lb=lci, ci.ub=uci, xlab="Relative Risk", atransf=exp)
>
> help(forest.default) will tell you more about using the forest() function.
>
> 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
> ________________________________________
> From: r-help-bounces_at_r-project.org [r-help-bounces_at_r-project.org] On
> Behalf Of cheba meier [cheba.meier_at_googlemail.com]
> Sent: Thursday, April 07, 2011 4:11 PM
> To: Scott Chamberlain
> Cc: R-help_at_r-project.org; tlumley_at_uw.edu
> Subject: Re: [R] metaplot
>
> Many thanks for your reply! the Package ‘ggplot’ was removed from the CRAN
> repository. Formerly available versions can be obtained from the
> archive<http://cran.r-project.org/src/contrib/Archive/ggplot
> >(Unfortunately,
> I am a windows user). The metafor and meta.summaries()
> still need the (a,b,c,d) or standard errors for the effect estimates! I
> have
> the following data
>
> rr <- c(0.61, 0.35, 1.16, 0.94, 1.16, 0.46, 0.99, 1.05, 1.15, 0.27)
> lci <- c(-1.69, -1.27, -0.24, -0.54, -1.01, -1.8,-0.89, -1.38, -1.05,-2.04)
> uci <- c(0.70, -0.82,0.54,0.41,1.30,0.24,0.87,1.48,1.33,-0.57)
> rrdata <- data.frame(rr,lci,uci)
>
> I can compute standard errors from my data? but the results show different
> confidence interval?
>
> Many thanks in advance,
> Cheba
>
>
> 2011/4/6 Scott Chamberlain <scttchamberlain4_at_gmail.com>
>
> > What about the metafor package?
> >
> >
> > Or just create your own plot.
> >
> > For example, using ggplot2 package:
> >
> > limits <- aes(ymax = OR + (OR - 95%LCI), ymin = OR - (OR - 95%LCI))
> > ggplot(dataframe, aes(x = Study.Name, y = OR)) + geom_point() +
> > geom_errobar(limits)
> >
> > Best,
> > Scott
> >
> > On Wednesday, April 6, 2011 at 11:53 AM, cheba meier wrote:
> >
> > Dear all,
> >
> > I have a four variable: Stuy.Name, OR, 95%LCI and 95%UCI and I would like
> > to
> > create a meta analysis plot. I can't use meta.MH function in metaplot
> > because I do not have
> > n.trt, n.ctrl, col.trt, col.ctrl are not available! Is there an
> alternative
> > way to do it?
> >
> > Many thanks in advance,
> > Cheba
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > 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.
>

        [[alternative HTML version deleted]]



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 Fri 08 Apr 2011 - 09:34:43 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 Fri 08 Apr 2011 - 10:00:27 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