What version of the metafor package are you using? The option of directly providing the CI bounds was added in version 1.5-0, so make sure you are using the newest version.

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<http://www.wvbauer.com/>

________________________________
From: cheba meier [cheba.meier@googlemail.com]
Sent: Friday, April 08, 2011 11:31 AM
To: Viechtbauer Wolfgang (STAT)
Cc: R-help@r-project.org
Subject: Re: [R] metaplot

Thank you for your reply.

> forest(log(rr), ci.lb<http://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@maastrichtuniversity.nl<mailto:wolfgang.viechtbauer@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<http://ci.lb>=lci, ci.ub=uci, xlab="Log Relative Risk")

or:

forest(log(rr), ci.lb<http://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@r-project.org<mailto:r-help-bounces@r-project.org> [r-help-bounces@r-project.org<mailto:r-help-bounces@r-project.org>] On Behalf Of cheba meier [cheba.meier@googlemail.com<mailto:cheba.meier@googlemail.com>]
Sent: Thursday, April 07, 2011 4:11 PM
To: Scott Chamberlain
Cc: R-help@r-project.org<mailto:R-help@r-project.org>; tlumley@uw.edu<mailto:tlumley@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@gmail.com<mailto:scttchamberlain4@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<mailto: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]]

