13.11.2012 Views

Introduction to Categorical Data Analysis

Introduction to Categorical Data Analysis

Introduction to Categorical Data Analysis

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

300 RANDOM EFFECTS: GENERALIZED LINEAR MIXED MODELS<br />

(yi1 = 1,yi2 = 1), whereas with a large negative ui it is common <strong>to</strong> see outcomes<br />

(yi1 = 0,yi2 = 0). When a high proportion of cases have these outcomes, the association<br />

between the repeated responses is positive. Greater association results from<br />

greater heterogeneity (i.e., larger σ ).<br />

10.1.3 Example: Sacrifices for the Environment Revisited<br />

Table 10.1 shows a 2 × 2 table from the General Social Survey, analyzed originally in<br />

Chapter 8. Subjects were asked whether, <strong>to</strong> help the environment, they were willing<br />

<strong>to</strong> (1) raise taxes, (2) accept a cut in living standards. The ML fit of model (10.3),<br />

treating {ui} as normal, yields ˆβ = 0.210 (SE = 0.130), with ˆσ = 2.85. For a given<br />

subject, the estimated odds of a “yes” response on accepting higher taxes equal<br />

exp(0.210) = 1.23 times the odds of a “yes” response on accepting a lower standard<br />

of living.<br />

Table 10.1. Opinions Relating <strong>to</strong> Environment<br />

Cut Living Standards<br />

Pay Higher Taxes Yes No Total<br />

Yes 227 132 359<br />

No 107 678 785<br />

Total 334 810 1144<br />

The relatively large ˆσ value of 2.85 reflects a strong association between the two<br />

responses. In fact, Table 10.1 has a sample odds ratio of 10.9. Whenever the sample<br />

log odds ratio in such a table is nonnegative, as it usually is, the ML estimate of β with<br />

this random effects approach is identical <strong>to</strong> the conditional ML estimate from treating<br />

{αi} in model (10.2) as fixed effects. Section 8.2.3 presented this conditional ML<br />

approach. For these data the conditional ML estimate is ˆβ = log(132/107) = 0.210,<br />

with SE =[(1/107) + (1/132)] 1/2 = 0.130.<br />

10.1.4 Differing Effects in Conditional Models and Marginal Models<br />

As Sections 9.1.3 and 8.2.2 discussed, parameters in GLMMs and marginal models<br />

have different interpretations. The parameters in GLMMs have conditional (clusterspecific)<br />

intepretations, given the random effect. By contrast, effects in marginal<br />

models are averaged over all clusters (i.e., population-averaged), and so those effects<br />

do not refer <strong>to</strong> a comparison at a fixed value of a random effect.<br />

Section 8.2.2 noted that the cluster-specific model (10.2) applies naturally <strong>to</strong> the<br />

data as displayed in a separate partial table for each cluster, displaying the two matched<br />

responses. For the survey data on the environmental issues, each subject is a cluster<br />

and has their own table. The first row shows the response on taxes (a 1 in the first<br />

column for “yes” or in the second column for “no”), and the second row shows the

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

Saved successfully!

Ooh no, something went wrong!