Problem Note 8764: PROC GLM may not compute the correct error sums of squares when solving
the Expected Mean Squares for RANDOM effects
If the TEST option is requested on the RANDOM statement and
if the reported ERROR term for a random effect test is
"almost all" MS(ERROR), then PROC GLM and PROC ANOVA will incorrectly
report the ERROR DF, ERROR SS, ERROR MS, F-Value, PR > F under the Tests
of Hypothesis for Mixed Model Analysis of Variance table.
The exact criterion for "almost all" is if
sum_i (l(i) * MS(i))**2 / DF(i) < 1e-8
where: l(i) is the coefficient of the mean square of the non-error
components
MS(i) is the mean square for the ith component, and
DF(i) is the degrees of freedom for the ith component
In the example below there is only one non-error component:
Error: 0.0001*MS(A*B) + 0.9999*MS(Error)
with
Source DF Type III SS Mean Square F Value Pr > F
A*B 11 0.01666472 0.00151497 1.21 0.3138
so
sum_i (l[i] * MS(i))**2 / DF(i)
~ (0.0001 * 0.00151497)**2 / 11 < 1e-8.
There is no circumvention to this problem in PROC GLM or PROC ANOVA.
However, you can use PROC MIXED with the METHOD=TYPE3 option on the
PROC statement to obtain the correct results. In general, all RANDOM
effects will go on the RANDOM statement (and not on the MODEL statement)
in PROC MIXED.
A Technical Support hot fix for Release 8.2 TSLEVEL TS2M0 for this
issue is available at:
http://www.sas.com/techsup/download/hotfix/82_sbcs_prod_list.html#008764
For customers running SAS with Asian Language Support (DBCS), this
hot fix should be downloaded from:
http://www.sas.com/techsup/download/hotfix/82_dbcs_prod_list.html#008764
Operating System and Release Information
| SAS System | SAS/STAT | HP-UX | 8 TS M0 | |
| Microsoft Windows XP Professional | 8.2 TS2M0 | 9.1 TS1M0 |
| Microsoft Windows 95/98 | 8 TS M0 | |
| Microsoft Windows NT Workstation | 8 TS M0 | 9.1 TS1M0 |
| Windows Millennium Edition (Me) | 8.2 TS2M0 | |
| Microsoft Windows 2000 Server | 8 TS M0 | 9.1 TS1M0 |
| Microsoft® Windows® for 64-Bit Itanium-based Systems | 9 TS M0 | 9.1 TS1M0 |
| Microsoft Windows 2000 Datacenter Server | 8 TS M0 | 9.1 TS1M0 |
| Microsoft Windows 2000 Professional | 8 TS M0 | 9.1 TS1M0 |
| OpenVMS VAX | 8 TS M0 | |
| Microsoft Windows 2000 Advanced Server | 8 TS M0 | 9.1 TS1M0 |
| 64-bit Enabled Solaris | 8 TS M0 | |
| Solaris | 8 TS M0 | |
| IRIX | 8 TS M0 | |
| z/OS | 8 TS M0 | 9.1 TS1M0 |
| OS/2 | 8 TS M0 | |
| ABI+ for Intel Architecture | 8 TS M0 | |
| Linux | 8.2 TS2M0 | 9.1 TS1M0 |
| 64-bit Enabled HP-UX | 8 TS M0 | |
| HP-UX IPF | 9 TS M0 | 9.1 TS1M0 |
| OpenVMS Alpha | 8 TS M0 | 9.1 TS1M0 |
| CMS | 8 TS M0 | |
| AIX | 8 TS M0 | |
| Tru64 UNIX | 8 TS M0 | 9.1 TS1M0 |
| 64-bit Enabled AIX | 8 TS M0 | |
*
For software releases that are not yet generally available, the Fixed
Release is the software release in which the problem is planned to be
fixed.
| Type: | Problem Note |
| Priority: | alert |
| Date Modified: | 2002-11-20 12:57:04 |
| Date Created: | 2002-10-30 09:34:16 |