OPTIONS NODATE NONUMBER PS=60 LS=72;

*******************************************************;
***  ANALYSIS OF VARIANCE FOR A LATIN SQUARE DESIGN ***;
*******************************************************;
*** FROM HICKS PAGE 77 ***;

DATA latin; INPUT car $ position brand $ wear @@; cards;
  I 1 C 12      I 2 B 14      I 3 A 17      I 4 D 13
 II 1 D 11     II 2 C 12     II 3 B 14     II 4 A 14
III 1 A 13    III 2 D 11    III 3 C 10    III 4 B 13
 IV 1 B  8     IV 2 A 13     IV 3 D  9     IV 4 C  9
;
PROC GLM DATA=latin;
     CLASS car position brand;
     MODEL wear = car position brand / SS3;
     MEANS position car;
     OUTPUT OUT=diag P=pred R=resid;
TITLE 'ANOVA RESULTS FOR LATIN SQUARE TIRE BRAND EXAMPLE';

PROC PRINT DATA=diag;