14.03.2014 Views

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

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.

Chapter 3 Fitting St<strong>and</strong>ard Least Squares Models 131<br />

Examples with Statistical Details<br />

Analysis of Covariance with Separate Slopes<br />

This example is a continuation of the Drug.jmp example presented in the previous section. The example<br />

uses data from Snedecor <strong>and</strong> Cochran (1967, p. 422). A one-way analysis of variance for a variable called<br />

Drug, shows a difference in the mean response among the levels a, d, <strong>and</strong> f, with a significance probability of<br />

0.03.<br />

The lack of fit test for the model with main effect Drug <strong>and</strong> covariate x is not significant. However, for the<br />

sake of illustration, this example includes the main effects <strong>and</strong> the Drug*x effect. This model tests whether<br />

the regression on the covariate has separate slopes for different Drug levels.<br />

This specification adds two columns to the linear model (call them x 4i <strong>and</strong> x 5i ) that allow the slopes for the<br />

covariate to be different for each Drug level. The new variables are formed by multiplying the dummy<br />

variables for Drug by the covariate values, giving the following formula:<br />

y i<br />

= β 0<br />

+ β 1<br />

x 1i<br />

+ β 2<br />

x 2i<br />

+ β 3<br />

x 3i<br />

+ β 4<br />

x 4i<br />

+ β 5<br />

x 5i<br />

+ ε i<br />

Table 3.17 on page 131, shows the coding of this Analysis of Covariance with Separate Slopes. (The mean<br />

of X is 10.7333.)<br />

Table 3.17 Coding of Analysis of Covariance with Separate Slopes<br />

Regressor Effect Values<br />

x 1 Drug[a] +1 if a, 0 if d, –1 if f<br />

x 2 Drug[d] 0 if a, +1 if d, –1 if f<br />

x 3 x the values of x<br />

x 4 Drug[a]*(x - 10.733) +x – 10.7333 if a, 0 if d, –(x – 10.7333) if f<br />

x 5 Drug[d]*(x - 10.733) 0 if a, +x – 10.7333 if d, –(x – 10.7333) if f<br />

A portion of the report is shown in Figure 3.56. The Regression Plot shows fitted lines with different slopes.<br />

The Effect Tests report gives a p-value for the interaction of 0.56. This is not significant statistically,<br />

indicating the model does not need to have different slopes.

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

Saved successfully!

Ooh no, something went wrong!