The object "a" has log.estimates and selog.estimates:
> str(a)
List of 10
$ logOR : num [1:7] -0.548 -1.804 -1.404 -0.357 -1.055 ...
$ selogOR : num [1:7] 0.22 1.11 0.611 0.372 0.715 ...
$ logMH : num -0.632
$ selogMH : num 0.16
snipped remainder of output.
You can also look at the the code of meta.MH:
meta.MH
On Jun 27, 2009, at 1:02 AM, sdzhangping wrote:
Dear R users:
In the example of meta-analysis (cochrane, package rmeta), I can not
found the p-value of Test for overall effect, and some other indices
(Z, I, weight and et al). How can I get the these indices listed?
library(rmeta)
data(cochrane)
cochrane
name ev.trt n.trt ev.ctrl n.ctrl
1 Auckland 36 532 60 538
2 Block 1 69 5 61
3 Doran 4 81 11 63
4 Gamsu 14 131 20 137
5 Morrison 3 67 7 59
6 Papageorgiou 1 71 7 75
7 Tauesch 8 56 10 71
a=meta.MH(n.trt,n.ctrl,ev.trt,ev.ctrl,names=name,data=cochrane)
summary(a)
Fixed effects ( Mantel-Haenszel ) meta-analysis
Call: meta.MH(ntrt = n.trt, nctrl = n.ctrl, ptrt = ev.trt, pctrl =
ev.ctrl,
names = name, data = cochrane)
------------------------------------
OR (lower 95% upper)
Auckland 0.58 0.38 0.89
Block 0.16 0.02 1.45
Doran 0.25 0.07 0.81
Gamsu 0.70 0.34 1.45
Morrison 0.35 0.09 1.41
Papageorgiou 0.14 0.02 1.16
Tauesch 1.02 0.37 2.77
------------------------------------
Mantel-Haenszel OR =0.53 95% CI ( 0.39,0.73 ) (where is Z and p-
value ?)
Test for heterogeneity: X^2( 6 ) = 6.9 ( p-value 0.3303 )
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
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.
