I have just installed the package as
> install.packages("metafor")
--- Please select a CRAN mirror for use in this session ---
.....
downloaded 403 Kb
package 'metafor' successfully unpacked and MD5 sums checked
> library("metafor")
Loading required package: Formula
Loading required package: nlme
Loading 'metafor' package (version 1.4-0). For an overview
and introduction to the package please type: help(metafor).
Warning message:
package 'metafor' was built under R version 2.11.1
> 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)
> library(metafor)
> forest(log(rr), ci.lb=lci, ci.ub=uci, xlab="Log Relative Risk")
Error in sei^2 : 'sei' is missing
Is it just because I am using R version 2.11.0?
Regards,
Cheba
2011/4/8 Viechtbauer Wolfgang (STAT) <
wolfgang.viechtba...@maastrichtuniversity.nl>
> 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
>
> ------------------------------
> *From:* cheba meier [cheba.me...@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=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.
