[R] questions about meta-analysis

[David Winsemius]

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