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.

Between columns :<br />

P …y:j y† 2 ˆ Pc<br />

jˆ1<br />

C 2 j =r T 2 =N:<br />

The Total SSq is similarly calculated as<br />

P …yij y† 2 ˆ P y 2 ij T 2 =N,<br />

and the Residual SSq may be obta<strong>in</strong>ed by subtraction:<br />

Residual SSq ˆ Total SSq Between-Rows SSq Between-Columns SSq:<br />

…9:5†<br />

The analysis so far is purely a consequence of algebraic identities. The relationships<br />

given above are true irrespective of the validity of the model. We now<br />

complete the analysis of variance by some steps which depend for their validity<br />

on that of the model. First, the degrees of freedom (DF) are allotted as shown <strong>in</strong><br />

Table 9.2. Those for rows and columns follow from the one-way analysis; if the<br />

only classification had been <strong>in</strong>to rows, for example, the first l<strong>in</strong>e of Table 9.2 would<br />

have been shown as Between groups and the SSq shown <strong>in</strong> Table 9.2 as Between<br />

columns and Residual would have added to form the With<strong>in</strong>-Groups SSq. With<br />

r 1 and c 1 as degrees of freedom for rows and columns, respectively, and<br />

N 1 for the Total SSq, the DF for Residual SSq follow by subtraction:<br />

…rc 1† …r 1† …c 1† ˆrc r c ‡ 1 ˆ…r 1†…c 1†:<br />

The mean squares (MSq) for rows, columns and residual are obta<strong>in</strong>ed <strong>in</strong> each<br />

case by the formula MSq ˆ SSq/DF, and those for rows and columns may each<br />

be tested aga<strong>in</strong>st the Residual MSq, s 2 , as shown <strong>in</strong> Table 9.2. The test for rows,<br />

for <strong>in</strong>stance, has the follow<strong>in</strong>g justification. On the null hypothesis (which we<br />

shall call HR) that all the row constants ai <strong>in</strong> (9.1) are equal (and therefore equal<br />

to zero, s<strong>in</strong>ce P ai ˆ 0), both s 2 R and s2 are unbiased estimates of s 2 .IfHR is not<br />

true, so that the ai differ, s 2 R has expectation greater than s2 whereas s 2 is still an<br />

unbiased estimate of s 2 . Hence FR tends to be greater than 1, and sufficiently<br />

high values <strong>in</strong>dicate a significant departure from HR. This test is valid whatever<br />

values the b j take, s<strong>in</strong>ce add<strong>in</strong>g a constant on to all the read<strong>in</strong>gs <strong>in</strong> a particular<br />

column has no effect on either s 2 R or s2 .<br />

Table 9.2 Two-way analysis of variance table.<br />

SSq DF MSq VR<br />

Between rows<br />

P<br />

i R2i =c T 2 =N r 1 s2 R FR ˆ s2 R =s2<br />

Between columns P<br />

j C2 j =r T 2 =N c 1 s2 C FC ˆ s2 C =s2<br />

Residual By subtraction …r 1†…c 1† s2 Total<br />

P<br />

y2 ij T 2 =N rc 1…ˆ N 1†<br />

i, j<br />

9.2 Two-way analysis of variance: randomized blocks 241

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

Saved successfully!

Ooh no, something went wrong!