07-02-2014 06:16 AM
How can we obtain using PROC MIXED the "Tests of Hypotheses for Mixed Model Analysis of Variance" table like the one which is obtained using PROC GLM (TEST option in RANDOM statement)? We have by default only the "Type 3 Tests of Fixed Effects".
Thanking you in advance for your answers.
07-02-2014 02:24 PM
To get variance component estimates use COVTEST option of the PROC MIXED statement: SAS/STAT(R) 13.1 User's Guide
proc mixed data=yourdata covtest;
Note that COVTEST option produces the asymptotic standard errors and Wald Z-tests for the covariance parameter estimates.
Hope this woud be helpful.
07-02-2014 02:39 PM
You cannot get tests like those as GLM depends on partitioning sums of squares and expected mean squares for its tests using random effects. MIXED uses a (RE)ML method to construct quadratic forms that depend on maximizing the likelihood of the observed values, given the fixed effects of the model and conditional on the random effect estimates. 's response will help if you are interested in the covariance parameter estimates only. The Wald tests of significance for the parameters are woefully inadequate in most cases due to lack of sample size. If you have sufficient data, using PROC GLIMMIX with Gauss-Hermitian quadrature methods will improve the accuracy of the estimates and the (still) asymptotic errors.