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.

290 MODELING CORRELATED, CLUSTERED RESPONSES<br />

Consider the insomnia study of Problem 9.3.2 from the previous section. Let<br />

y1 be the initial time <strong>to</strong> fall asleep, let Y2 be the follow-up time, with explana<strong>to</strong>ry<br />

variable x defining the two treatment groups (1 = active drug, 0 = placebo). We will<br />

now treat Y2 as an ordinal response and y1 as an explana<strong>to</strong>ry variable, using scores<br />

{10, 25, 45, 75}. In the model<br />

logit[P(Y2 ≤ j)]=αj + β1x + β2y1<br />

(9.2)<br />

β1 compares the follow-up distributions for the treatments, controlling for the initial<br />

observation. This models the follow-up response (Y2), conditional on y1, rather than<br />

marginal distributions of (Y1,Y2). It’s the type of model for an ordinal response that<br />

Section 6.2 discussed. Here, the initial response y1 plays the role of an explana<strong>to</strong>ry<br />

variable.<br />

From software for ordinary cumulative logit models, the ML treatment effect<br />

estimate is ˆβ1 = 0.885 (SE = 0.246). This provides strong evidence that follow-up<br />

time <strong>to</strong> fall asleep is lower for the active drug group. For any given value for the<br />

initial response, the estimated odds of falling asleep by a particular time for the active<br />

treatment are exp(0.885) = 2.4 times those for the placebo group. Exercise 9.12<br />

considers alternative analyses for these data.<br />

9.4.4 Transitional Models Relate <strong>to</strong> Loglinear Models<br />

Effects in transitional models differ from effects in marginal models, both in magnitude<br />

and in their interpretation. The effect of a predic<strong>to</strong>r xj on Yt is conditional<br />

on yt−1 in a transitional model, but it ignores yt−1 in a marginal model. Effects in<br />

transitional models are often considerably weaker than effects in marginal models,<br />

because conditioning on a previous response attenuates the effect of a predic<strong>to</strong>r.<br />

Transitional models have connections with the loglinear models of Chapter 7,<br />

which described joint distributions. Associations in loglinear models are conditional<br />

on the other response variables. In addition, a joint distribution of (Y1,Y2,...,YT ) can<br />

be fac<strong>to</strong>red in<strong>to</strong> the distribution of Y1, the distribution of Y2 given Y1, the distribution<br />

of Y3 given Y1 and Y2, and so forth.<br />

PROBLEMS<br />

9.1 Refer <strong>to</strong> Table 7.3 on high school students’ use of alcohol, cigarettes, and<br />

marijuana. View the table as matched triplets.<br />

a. Construct the marginal distribution for each substance. Find the sample<br />

proportions of students who used (i) alcohol, (ii) cigarettes, (iii) marijuana.<br />

b. Specify a marginal model that could be fitted as a way of comparing the<br />

margins. Explain how <strong>to</strong> interpret the parameters in the model. State the<br />

hypothesis, in terms of the model parameters, that corresponds <strong>to</strong> marginal<br />

homogeneity.

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

Saved successfully!

Ooh no, something went wrong!