DM 'LOG; CLEAR; OUT; CLEAR;';

ODS LISTING;
ODS PRINTER PDF file='C:\COURSES\ST578\SAS\TWO3.PDF';
OPTIONS NODATE NONUMBER;

OPTIONS PS=54 LS=76 NODATE NONUMBER;

DATA IN;
   DO C = -1 TO 1 BY 2;
   DO B = -1 TO 1 BY 2;
   DO A = -1 TO 1 BY 2;
   DO REP = 1 TO 3;
      INPUT Y @@; 
      X1=A;  X2=B;  X3=C; 
      X1X2 = X1*X2;  X1X3 = X1*X3;  X2X3 = X2*X3;
      OUTPUT;
   END; END; END; END;
LINES;
22 31 25   32 43 29   35 34 50   55 47 46   
44 45 38   40 37 36   60 50 54   39 41 47
;

PROC GLM DATA=IN PLOTS=NONE;
  CLASS A B C;
  MODEL Y = A B C / SS3;
  MEANS A B C;
  ESTIMATE 'A' A -1 1;
  ESTIMATE 'B' B -1 1;
  ESTIMATE 'C' C -1 1;
TITLE 'MODEL 1: ADDITIVE FIXED EFFECTS MODEL';
/*
PROC REG DATA=IN PLOTS=(DIAGNOSTICS);
     MODEL Y = X1 X2 X3 / LACKFIT VIF;
TITLE 'MODEL 2: FIRST ORDER REGRESSION MODEL';
*/
PROC GLM DATA=IN PLOTS=NONE;
  CLASS A B C;
  MODEL Y = A|B|C@2 / SS3 ;
  MEANS A|B|C@2;
  ESTIMATE 'A' A -1 1;
  ESTIMATE 'B' B -1 1;
  ESTIMATE 'C' C -1 1;
  ESTIMATE 'A*B' A*B 1 -1 -1 1 / DIVISOR=2;
  ESTIMATE 'A*C' A*C 1 -1 -1 1 / DIVISOR=2;
  ESTIMATE 'B*C' B*C 1 -1 -1 1 / DIVISOR=2;
* ESTIMATE 'A*B*C' A*B*C -1 1 1 -1 1 -1 -1 1 ;
TITLE 'MODEL 3: INTERACTION FIXED EFFECTS MODEL';
/*
PROC REG DATA=IN PLOTS=(DIAGNOSTICS);
     MODEL Y = X1 X2 X3 X1X2 X1X3 X2X3 / LACKFIT VIF;
TITLE 'MODEL 4: INTERACTION REGRESSION MODEL';
*/
RUN;