14.03.2014 Views

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Appendix A Statistical Details 657<br />

The Factor Models<br />

Nominal Factors<br />

Nominal factors are transformed into indicator variables for the design matrix. <strong>SAS</strong> GLM constructs an<br />

indicator column for each nominal level. JMP constructs the same indicator columns for each nominal level<br />

except the last level. When the last nominal level occurs, a one is subtracted from all the other columns of<br />

the factor. For example, consider a nominal factor A with three levels coded for GLM <strong>and</strong> for JMP as shown<br />

below.<br />

Table A.1 Nominal Factor A<br />

GLM<br />

JMP<br />

A A1 A2 A3 A13 A23<br />

A1 1 0 0 1 0<br />

A2 0 1 0 0 1<br />

A3 0 0 1 –1 –1<br />

In GLM, the linear model design matrix has linear dependencies among the columns, <strong>and</strong> the least squares<br />

solution employs a generalized inverse. The solution chosen happens to be such that the A3 parameter is set<br />

to zero.<br />

In JMP, the linear model design matrix is coded so that it achieves full rank unless there are missing cells or<br />

other incidental collinearity. The parameter for the A effect for the last level is the negative sum of the other<br />

levels, which makes the parameters sum to zero over all the effect levels.<br />

Interpretation of Parameters<br />

Note: The parameter for a nominal level is interpreted as the differences in the predicted response for that<br />

level from the average predicted response over all levels.<br />

The design column for a factor level is constructed as the zero-one indicator of that factor level minus the<br />

indicator of the last level. This is the coding that leads to the parameter interpretation above.<br />

Table A.2 Interpreting Parameters<br />

JMP Parameter Report How to Interpret Design Column Coding<br />

Intercept mean over all levels 1´<br />

A[1]<br />

A[2]<br />

α 1<br />

– 1 ⁄ 3( α 1<br />

+ α 2<br />

+ α 3<br />

)<br />

α 2<br />

– 1 ⁄ 3( α 1<br />

+ α 2<br />

+ α 3<br />

)<br />

(A==1) – (A==3)<br />

(A==2) – (A==3)

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

Saved successfully!

Ooh no, something went wrong!