[R] How to do a meta-analysis plot

Gabor Grothendieck ggrothendieck at gmail.com
Mon Aug 25 20:50:01 CEST 2008


See ?metaplot

On Mon, Aug 25, 2008 at 2:44 PM, Jorge Ivan Velez
<jorgeivanvelez at gmail.com> wrote:
> Dear R-list,
>
> I'd like to do a meta-analysis plot similar to
>
> install.packages('rmeta')
> require(rmeta)
> data(catheter)
> a <- meta.MH(n.trt, n.ctrl, col.trt, col.ctrl, data=catheter,
>             names=Name, subset=c(13,6,5,3,7,12,4,11,1,8,10,2))
> summary(a)
> plot(a)
>
>
> (see attached file) by using my own OR (Odds Ratio) and 95% Confidence
> Interval data set, which looks like
>
> mydata=data.frame(OR=c(2.04545454545, 1.10434782609, 1.22588104401,
> 1.14102564103,
> 1.20579245527, 1.375, 1.16535433071),
> L95=c(1.22839621997, 0.858106819302, 1.0964802088, 0.841934120955,
> 0.969786886818, 1.01498537023, 0.919391492382),
> U95=c(3.40546755139, 1.42122051928, 1.37055308613, 1.54632513827,
> 1.49917372998, 1.86258857302, 1.47707220868)
> )
> rownames(mydata)=c(paste("Study",1:6,sep=""),'Summary')
> mydata
>
>
> My problem is that I don't have the raw data as rmeta _requires_ and, even
> when I have my data set in the _same_ (?) format that summary(a), when I
> tried plot(mydata) it doesn't work. Another approach I used was to change
> the class of my object but it didn't work either. I'm running XP SP2 on a
> 2.4 GHz Intel-Core 2 Duo processor and my R-session info is the following:
>
> R version 2.7.2 RC (2008-08-18 r46388)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
> States.1252;LC_MONETARY=English_United
> States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
>
> attached base packages:
> [1] grid      stats     graphics  grDevices utils     datasets  methods
> base
>
> other attached packages:
> [1] rmeta_2.14  RODBC_1.2-3
>
> loaded via a namespace (and not attached):
> [1] tools_2.7.2
>
>
> I would greatly appreciate any ideas about how should I proceed.
>
> Thanks in advance,
>
>
> Jorge
>
> ______________________________________________
> 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