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
example MH.pdf
Description: Adobe PDF document
______________________________________________ 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.
