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.viechtba...@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-boun...@r-project.org [r-help-boun...@r-project.org] On > Behalf Of cheba meier [cheba.me...@googlemail.com] > Sent: Thursday, April 07, 2011 4:11 PM > To: Scott Chamberlain > Cc: R-help@r-project.org; tlum...@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 <scttchamberla...@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@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@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.
