From yen@stat.wisc.edu Wed Mar  4 21:10:02 1998
Date: Wed, 4 Mar 1998 19:59:44 -0600 (CST)
From: Chen-Hsing Yen <yen@stat.wisc.edu>
To: yandell@stat.wisc.edu
Cc: Chen-Hsing Yen <yen@stat.wisc.edu>
Subject: About data of ex26.2 in M&J

Hi, Prof. Yandell
The following is the SAS program I use. Problem is that the result 
is not exactly the same as or close to the anova table in the book.
Maybe you will have better solution. 

Chen-Hsing Yen
================================================================
options nodate ls=79 ps=79;
data new;
infile 'comforta.dat' firstobs=2;
input obs envir rep sex$ cloth time score;

proc mixed;
class envir rep sex cloth time;
model score = envir rep(envir)
              sex cloth sex*cloth envir*sex envir*cloth envir*sex*cloth
              sex*cloth*rep(envir)
              time time*envir time*sex time*cloth time*sex*cloth
              time*envir*sex time*envir*cloth time*envir*sex*cloth;
random rep(envir) sex*cloth*rep(envir);
test h=envir e=rep(envir);
test h=sex cloth sex*cloth envir*sex envir*cloth envir*sex*cloth
     e=sex*cloth*rep(envir);