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.

94 GENERALIZED LINEAR MODELS<br />

Table 3.7. Table for Problem 3.10 on Cancer Remission<br />

Standard Likelihood Ratio 95% Chi-<br />

Parameter Estimate Error Confidence Limits Square Pr > ChiSq<br />

Intercept −2.3178 0.7795 −4.0114 −0.9084 8.84 0.0029<br />

LI 0.0878 0.0328 0.0275 0.1575 7.19 0.0073<br />

as independent Poisson variates having means μA and μB. Consider the model<br />

log μ = α + βx, where x = 1 for treatment B and x = 0 for treatment A.<br />

a. Show that β = log μB − log μA = log(μB/μA) and e β = μb/μA.<br />

b. Fit the model. Report the prediction equation and interpret ˆβ.<br />

c. Test H0: μA = μB by conducting the Wald or likelihood-ratio test of<br />

H0: β = 0. Interpret.<br />

d. Construct a 95% confidence interval for μB/μA. [Hint: Construct one for<br />

β = log(μB/μA) and then exponentiate.]<br />

3.12 Refer <strong>to</strong> Problem 3.11. The wafers are also classified by thickness of silicon<br />

coating (z = 0, low; z = 1, high). The first five imperfection counts reported<br />

for each treatment refer <strong>to</strong> z = 0 and the last five refer <strong>to</strong> z = 1. Analyze these<br />

data, making inferences about the effects of treatment type and of thickness of<br />

coating.<br />

3.13 Access the horseshoe crab data of Table 3.2 at www.stat.ufl.edu/∼aa/introcda/appendix.html.<br />

a. Using x = weight and Y = number of satellites, fit a Poisson loglinear<br />

model. Report the prediction equation.<br />

b. Estimate the mean of Y for female crabs of average weight 2.44 kg.<br />

c. Use ˆβ <strong>to</strong> describe the weight effect. Construct a 95% confidence interval<br />

for β and for the multiplicative effect of a 1 kg increase.<br />

d. Conduct a Wald test of the hypothesis that the mean of Y is independent of<br />

weight. Interpret.<br />

e. Conduct a likelihood-ratio test about the weight effect. Interpret.<br />

3.14 Refer <strong>to</strong> the previous exercise. Allow overdispersion by fitting the negative<br />

binomial loglinear model.<br />

a. Report the prediction equation and the estimate of the dispersion parameter<br />

and its SE. Is there evidence that this model gives a better fit than the<br />

Poisson model?<br />

b. Construct a 95% confidence interval for β. Compare it with the one in (c)<br />

in the previous exercise. Interpret, and explain why the interval is wider<br />

with the negative binomial model.

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

Saved successfully!

Ooh no, something went wrong!