09.08.2013 Views

Fundamentals of epidemiology - an evolving text - Are you looking ...

Fundamentals of epidemiology - an evolving text - Are you looking ...

Fundamentals of epidemiology - an evolving text - Are you looking ...

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.

Moreover, the logistic model, we will see, corresponds to a multiplicative model, which we saw earlier is<br />

the model that is implied by stratified <strong>an</strong>alysis based on the OR or the risk ratio. Furthermore, the<br />

coefficients that we estimate using logistic regression c<strong>an</strong> be converted into OR's, so that we now have<br />

a ratio measure <strong>of</strong> association.<br />

It is easy to discover what the logistic coefficients are. Since the logit is the logarithm <strong>of</strong> the odds, then<br />

the difference <strong>of</strong> two logits is the logarithm <strong>of</strong> <strong>an</strong> OR (because subtraction <strong>of</strong> logs corresponds to<br />

division <strong>of</strong> their arguments – see the appendix to the chapter on Measures <strong>of</strong> Frequency <strong>an</strong>d Extent).<br />

Suppose that X3 is a dichotomous (0-1) variable indicating absence (0) or presence (1) <strong>of</strong> <strong>an</strong> exposure.<br />

First write the model with the exposure "present" (X3=1), <strong>an</strong>d underneath write the model with the<br />

exposure "absent" (X3=0).<br />

logit(D=1|X1,X2,X3=1) = α + β1X1 + β2X2 + β3 (X3 =1, present)<br />

– logit(D=1|X1,X2,X3=0) = α + β1X1 + β2X2 + 0 (X3 = 0, absent)<br />

_______________________________________________________<br />

When we subtract the second model from the first, all the terms on the right are removed except the<br />

coefficient for X3. On the left, we have the (rather messy) difference <strong>of</strong> the two logits, one for X3<br />

present <strong>an</strong>d the other for X3 absent:<br />

logit(D=1|X1,X2,X3=1) – logit(D=1|X1,X2,X3=0) = β3<br />

Spelling out the logits:<br />

ln(odds(D=1|X1,X2,X3=1)) – ln(odds(D=1|X1,X2,X3=0)) = β3<br />

<strong>an</strong>d, since a difference <strong>of</strong> logarithms is the logarithm <strong>of</strong> a ratio:<br />

odds(D=1|X1,X2,X3=1)<br />

ln [ ———————————— ] = β3<br />

odds(D=1|X1,X2,X3=0)<br />

A ratio <strong>of</strong> odds is simply <strong>an</strong> OR, in this case, the OR for the disease with respect to the exposure<br />

represented by X3:<br />

ln [ OR ] = β3<br />

exp (ln [ OR ] ) = exp(β3)<br />

________________________________________________________________________________________________<br />

www.sph.unc.edu/courses/EPID 168, © Victor J. Schoenbach 13. Multicausality ― <strong>an</strong>alysis approaches ― 446<br />

rev. 10/28/1999, 11/16/2000, 4/2/2001

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

Saved successfully!

Ooh no, something went wrong!