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.

128 Analys<strong>in</strong>g means and proportions<br />

where A…c† is the value of a which, with the observed marg<strong>in</strong>al totals, would<br />

give an odds ratio of c, and var (a; c) is the variance of a for that value of c.<br />

That is,<br />

A…d a ‡ A†<br />

ˆ c …4:28†<br />

…a ‡ b A†…a ‡ c A†<br />

and<br />

var…a; c† ' 1<br />

A ‡<br />

1<br />

a ‡ c A ‡<br />

1<br />

a ‡ b A ‡<br />

1<br />

1<br />

: …4:29†<br />

d a ‡ A<br />

It is tedious to solve (4.27), but the calculation can readily be set up as a<br />

spreadsheet calculation and solved quickly by trial and error. Us<strong>in</strong>g a trial<br />

value for A, c is obta<strong>in</strong>ed from (4.28) and var…a; c† from (4.29). These are then<br />

substituted <strong>in</strong> the left-hand side of (4.27). The aim is to choose values of A so that<br />

(4.27) gives 1 96. This is achieved by try<strong>in</strong>g different values and iterat<strong>in</strong>g. See<br />

also Breslow and Day (1980, §4.3) or Fleiss (1981, §5.6).<br />

Example 4.11<br />

Liddell et al. (1984) reported on a case±control study <strong>in</strong>vestigat<strong>in</strong>g the association of<br />

bronchial carc<strong>in</strong>oma and asbestos exposure <strong>in</strong> the Canadian chrysolite m<strong>in</strong>es and mills.<br />

The data were as follows:<br />

The calculations are:<br />

Asbestos exposure Lung cancer Controls<br />

Exposed 148 (a) 372 (c)<br />

Not exposed 75 (b) 343 (d)<br />

OR ˆ…148<br />

ln OR ˆ 0 599<br />

343†=…75 372† ˆ1 82<br />

r<br />

SE…ln OR† ˆ<br />

1 1 1 1<br />

‡ ‡ ‡<br />

148 75 372 343<br />

p<br />

ˆ 0 02569 ˆ 0 160:<br />

95% confidence <strong>in</strong>terval for ln OR ˆ 0 599 1 96 0 160<br />

ˆ 0 285 to 0 913,<br />

95% confidence <strong>in</strong>terval for OR ˆ exp…0 285† to exp …0 913†<br />

ˆ 1 33 to 2 49:<br />

Cornfield's method, us<strong>in</strong>g (4.27) ± (4.29), gives the same limits of 1 33 and 2 49.<br />

2 2 tables and x 2 tests<br />

An alternative way of display<strong>in</strong>g the data of Example 4.10 is shown <strong>in</strong> Table 4.5.<br />

This is called a 2 2, or sometimes a fourfold, cont<strong>in</strong>gency table. The total

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

Saved successfully!

Ooh no, something went wrong!