[R] How to do a meta-analysis plot

Michael Dewey info at aghmed.fsnet.co.uk
Tue Aug 26 13:16:36 CEST 2008


At 19:44 25/08/2008, Jorge Ivan Velez wrote:
>Dear R-list,
>
>I'd like to do a meta-analysis plot similar to

Since these plots are known as forest plots
?forestplot
might help you.
As it says you have to do a bit more work but you do get much more flexibility


>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
>
>Content-Type: application/pdf; name="example MH.pdf"
>X-Attachment-Id: f_fkbfhzhb0
>Content-Disposition: attachment; filename="example MH.pdf"

Michael Dewey
http://www.aghmed.fsnet.co.uk



More information about the R-help mailing list