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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Table 19.6 Subject-years at risk ( y) of asbestos workers and death rates (d ) per 100 000 for lung<br />

cancer <strong>in</strong> men <strong>in</strong> England and Wales.<br />

Period<br />

1941±50 1951±60 1961±70 1971±80<br />

Age group y d y d y d y d<br />

25±29 82 2 62 1 13 1 3 1<br />

30±34 148 3 273 4 156 3 43 2<br />

35±39 74 9 446 10 435 8 141 6<br />

40±44 41 21 395 25 677 22 290 16<br />

45±49 33 46 229 58 749 54 485 46<br />

50±54 23 78 172 124 590 119 642 106<br />

55±59 14 112 158 216 399 226 621 201<br />

60±64 11 137 109 294 288 370 479 346<br />

65±69 4 137 78 343 185 508 273 530<br />

70±74 0 107 47 325 124 562 151 651<br />

75‡ 0 86 16 270 58 518 72 756<br />

The expected number of deaths due to lung cancer is<br />

…82 2 ‡ 62 1 ‡ 13 1 ‡ ...‡ 72 756† 10 5 ˆ 13 8:<br />

The observed number was 67, so the SMR is 67=13 8 ˆ 4 9. Us<strong>in</strong>g the methods given <strong>in</strong><br />

§5.2, a test of the null hypothesis that the asbestos workers experienced national death<br />

rates is<br />

z ˆ…67 13 8†= 1<br />

p 3 8 ˆ 14 3 …P < 0 001†,<br />

and approximate 95% confidence limits for the SMR are 3 8 and 6 2.<br />

The method has usually been applied to compare observed and expected<br />

mortality with<strong>in</strong> s<strong>in</strong>gle groups but may be extended to compare the SMR<br />

between different subgroups or, more generally, to take account of covariates<br />

recorded for each <strong>in</strong>dividual, by express<strong>in</strong>g the SMR as a proportional-hazards<br />

regression model (Berry, 1983)Ðthat is, by the use of a generalized l<strong>in</strong>ear model<br />

(§14.4). For subgroup i, if mi is the cumulative hazard from the reference<br />

population and g i the proportional-hazards multiplier, then m i, the expected<br />

number of deaths, is given by<br />

ln m i ˆ ln mi ‡ ln g i:<br />

If g i is modelled <strong>in</strong> terms of a set of covariates by<br />

then<br />

ln g i ˆ b 0 ‡ b 1x i1 ‡ b 2x i2 ‡ ...‡ bpxip<br />

19.7 Subject-years method 687

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

Saved successfully!

Ooh no, something went wrong!