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.

682 Statistical Details Appendix A<br />

<strong>Multivariate</strong> Details<br />

Table A.17 Approximate F-statistics (Continued)<br />

Hotelling-Lawley<br />

Trace<br />

F<br />

=<br />

---------------------------------<br />

2( sn + 1)U<br />

s 2 ( 2m + s+<br />

1)<br />

s( 2m+ s + 1) 2( sn + 1)<br />

Roy’s Max Root<br />

F<br />

=<br />

Θ<br />

--------------------------------------------------<br />

( v – max( pq , ) + q)<br />

max( pq , )<br />

max( pq , ) v – max( pq , ) + q<br />

Canonical Details<br />

The canonical correlations are computed as<br />

ρ i<br />

=<br />

λ<br />

------------- i<br />

1 + λ i<br />

The canonical Y’s are calculated as<br />

Ỹ<br />

=<br />

YMV<br />

where Y is the matrix of response variables, M is the response design matrix, <strong>and</strong> V is the matrix of<br />

eigenvectors of E – 1 H . Canonical Y’s are saved for eigenvectors corresponding to eigenvalues larger than<br />

zero.<br />

The total sample centroid is computed as<br />

Gr<strong>and</strong> = yMV<br />

where V is the matrix of eigenvectors of E – 1 H .<br />

The centroid values for effects are calculated as<br />

E<br />

m = ( c' 1 x j , c' 2 x j , …,<br />

c' g x j ) wherec i = v' i<br />

-----------<br />

N–<br />

r<br />

v – 1 ⁄ 2<br />

<br />

i <br />

v i<br />

<strong>and</strong> the vs are columns of V, the eigenvector matrix of E – 1 H , x j refers to the multivariate least squares<br />

mean for the jth effect, g is the number of eigenvalues of E – 1 H greater than 0, <strong>and</strong> r is the rank of the X<br />

matrix.<br />

The centroid radii for effects are calculated as<br />

d =<br />

2<br />

χ g ( 0.95)<br />

-----------------------------<br />

LX'X ( ) – 1 L'<br />

where g is the number of eigenvalues of E – 1 H greater than 0 <strong>and</strong> the denominator L’s are from the<br />

multivariate least squares means calculations.

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

Saved successfully!

Ooh no, something went wrong!