01.06.2013 Views

Statistical Methods in Medical Research 4ed

Statistical Methods in Medical Research 4ed

Statistical Methods in Medical Research 4ed

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

216 Comparison of several groups<br />

should be regarded as heterogeneous, <strong>in</strong> the sense that sampl<strong>in</strong>g variation is<br />

unlikely to expla<strong>in</strong> the differences among them.<br />

We now consider a method of wide generality for tackl<strong>in</strong>g other problems of<br />

this sort. It depends on certa<strong>in</strong> assumptions, which may or may not be precisely<br />

true <strong>in</strong> any particular application, but it is nevertheless very useful <strong>in</strong> provid<strong>in</strong>g<br />

approximate solutions to many problems. The general situation is that data are<br />

available <strong>in</strong> k groups, each provid<strong>in</strong>g an estimate of some parameter. We wish<br />

first to test whether there is evidence of heterogeneity between the estimates, and<br />

then <strong>in</strong> the absence of heterogeneity to obta<strong>in</strong> a s<strong>in</strong>gle estimate of the parameter<br />

from the whole data set.<br />

Suppose we observe k quantities, Y1, Y2, ..., Yk, and we know that Yi is<br />

N…mi, Vi†, where the means mi are unknown but the variances Vi are known. To<br />

test the hypothesis that all the mi are equal to some specified value m, we could<br />

calculate a weighted sum of squares<br />

G0 ˆ P<br />

‰…Yi m† 2 =ViŠ ˆ P<br />

wi…Yi m† 2 , …8:12†<br />

i<br />

where wi ˆ 1=Vi. The quantity G0 is the sum of squares of k standardized normal<br />

deviates and, from (5.2), it follows the x 2 …k† distribution.<br />

However, to test the hypothesis of homogeneity we do not usually wish to<br />

specify the value m. Let us replace m by the weighted mean<br />

and calculate<br />

P<br />

Y ˆ<br />

i wiYi<br />

P<br />

i wi<br />

i<br />

, …8:13†<br />

G ˆ P<br />

wi…Yi Y† 2 : …8:14†<br />

i<br />

The quantity wi is called a weight: note that it is the reciprocal of the variance Vi,<br />

so if, for example, Y1 is more precise than Y2, then V1 < V2 and w1 > w2 and Y1<br />

will be given the higher weight.<br />

A little algebra gives the alternative formula:<br />

G ˆ P<br />

i<br />

wiY 2<br />

i<br />

… P<br />

i wiYi† 2 = P<br />

i wi: …8:15†<br />

Note that if all the wi ˆ 1, (8.14) is the usual sum of squares about the mean, Y<br />

becomes the usual unweighted mean and (8.15) becomes the usual short-cut<br />

formula (2.3).<br />

It can be shown that, on the null hypothesis of homogeneity, G is distributed<br />

as x2 …k 1† . Replac<strong>in</strong>g m by Y has resulted <strong>in</strong> the loss of one degree of freedom.<br />

High values of G <strong>in</strong>dicate evidence aga<strong>in</strong>st homogeneity.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!