13.11.2012 Views

Introduction to Categorical Data Analysis

Introduction to Categorical Data Analysis

Introduction to Categorical Data Analysis

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.

9.4 TRANSITIONAL MODELING, GIVEN THE PAST 289<br />

Table 9.8. Child’s Respira<strong>to</strong>ry Illness by Age and Maternal<br />

Smoking<br />

No Maternal Maternal<br />

Smoking Smoking<br />

Child’s Respira<strong>to</strong>ry Illness Age 10 Age 10<br />

Age 7 Age 8 Age 9 No Yes No Yes<br />

No No No 237 10 118 6<br />

Yes 15 4 8 2<br />

Yes No 16 2 11 1<br />

Yes 7 3 6 4<br />

Yes No No 24 3 7 3<br />

Yes 3 2 3 1<br />

Yes No 6 2 4 2<br />

Yes 5 11 4 7<br />

Source: Thanks <strong>to</strong> Dr. James Ware for these data.<br />

Let yt denote the response on respira<strong>to</strong>ry illness at age t. For the regressive logistic<br />

model<br />

logit[P(Yt = 1)] =α + βyt−1 + β1s + β2t, t = 8, 9, 10<br />

each subject contributes three observations <strong>to</strong> the model fitting. The data set consists<br />

of 12 binomials, for the 2 × 3 × 2 combinations of (s, t, yt−1). For instance, for the<br />

combination (0, 8, 0), from Table 9.8 we see that y8 = 0 for 237 + 10 + 15 + 4 = 266<br />

subjects and y8 = 1 for 16 + 2 + 7 + 3 = 28 subjects.<br />

The ML fit of this regressive logistic model is<br />

logit[ ˆ<br />

P(Yt = 1)] =−0.293 + 2.210yt−1 + 0.296s − 0.243t<br />

The SE values are 0.158 for the yt−1 effect, 0.156 for the s effect, and 0.095 for<br />

the t effect. Not surprisingly, yt−1 has a strong effect – a multiplicative impact of<br />

e 2.21 = 9.1 on the odds. Given that and the child’s age, there is slight evidence of a<br />

positive effect of maternal smoking: The likelihood-ratio statistic for H0: β1 = 0is<br />

3.55 (df = 1, P = 0.06). The maternal smoking effect weakens further if we add<br />

yt−2 <strong>to</strong> the model (Problem 9.13).<br />

9.4.3 Comparisons that Control for Initial Response<br />

The transitional type of model can be especially useful for matched-pairs data. The<br />

marginal models that are the main focus of this chapter would evaluate how the marginal<br />

distributions of Y1 and Y2 depend on explana<strong>to</strong>ry variables. It is often more relevant<br />

<strong>to</strong> treat Y2 as a univariate response, evaluating effects of explana<strong>to</strong>ry variables while<br />

controlling for the initial response y1. That is the focus of a transitional model.

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

Saved successfully!

Ooh no, something went wrong!