[R] How to do a meta-analysis plot

[Jorge Ivan Velez]

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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: example MH.pdf
Type: application/pdf
Size: 4022 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080825/6b7256ee/attachment.pdf>