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.

8.7 Comparison of several variances<br />

The one-way analysis of variance (§8.1) is a generalization of the two-sample t<br />

test (§4.3). Occasionally one requires a generalization of the F test (used, as <strong>in</strong><br />

§5.1, for the comparison of two variances) to the situation where more than two<br />

estimates of variance are to be compared. In a one-way analysis of variance, for<br />

example, the primary purpose is to compare means, but one might wish to test<br />

the significance of differences between variances, both for the <strong>in</strong>tr<strong>in</strong>sic <strong>in</strong>terest of<br />

this comparison and also because the analysis of variance <strong>in</strong>volves an assumption<br />

that the group variances are equal.<br />

Suppose there are k estimates of variance, s2 i , hav<strong>in</strong>g possibly different<br />

degrees of freedom, ni. (If the ith group conta<strong>in</strong>s ni observations, ni ˆ ni 1.)<br />

On the assumption that the observations are randomly selected from normal<br />

distributions, an approximate significance test due to Bartlett (1937) consists <strong>in</strong><br />

calculat<strong>in</strong>g<br />

s 2 ˆ P nis2 P<br />

i = ni<br />

and<br />

1<br />

C ˆ 1 ‡<br />

3…k 1†<br />

8.7 Comparison of several variances 233<br />

M ˆ P … ni†ln<br />

s 2 P ni ln s2 i<br />

X 1<br />

ni<br />

1<br />

P ni<br />

and referr<strong>in</strong>g M=C to the x2 …k 1† distribution. Here `ln' refers to the natural<br />

logarithm (see p. 126). The quantity C is likely to be near 1 and need be<br />

calculated only <strong>in</strong> marg<strong>in</strong>al cases. Worked examples are given by Snedecor and<br />

Cochran (1989, §13.10).<br />

Bartlett's test is perhaps less useful than might be thought, for two reasons.<br />

First, like the F test, it is rather sensitive to non-normality. Secondly, with<br />

have to differ very considerably<br />

samples of moderate size, the true variances s2 i<br />

before there is a reasonable chance of obta<strong>in</strong><strong>in</strong>g a significant test result. To put<br />

this po<strong>in</strong>t another way, even if M=C is non-significant, the estimated s2 i may<br />

differ substantially, and so may the true s2 i . If possible <strong>in</strong>equality <strong>in</strong> the s2 i is<br />

important, it may therefore be wise to assume it even if the test result is nonsignificant.<br />

In some situations moderate <strong>in</strong>equality <strong>in</strong> the s2 i will not matter very<br />

much, so aga<strong>in</strong> the significance test is not relevant.<br />

An alternative test which is less <strong>in</strong>fluenced by non-normality is due to Levene<br />

(1960). In this test the deviations of each value from its group mean, or median,<br />

are calculated and negative deviations changed to positive, that is, the absolute<br />

deviations are used. Then a test of the equality of the mean values of the<br />

absolute deviations over the groups us<strong>in</strong>g a one-way analysis of variance (§8.1)<br />

is carried out. S<strong>in</strong>ce the mean value of the absolute deviation is proportional to

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

Saved successfully!

Ooh no, something went wrong!