***********************************************;
*** 2^2 Design: Bit size and cutting speed. ***;
*** Montgomery (2001) Problem 6-5, page 277 ***;
***********************************************;

DATA IN;
   DO B = -1 TO 1 BY 2;
   DO A = -1 TO 1 BY 2;
   DO REP = 1 TO 4;
      INPUT Y @@; OUTPUT;
   END; END; END;
CARDS;
18.2 18.9 12.9 14.4    27.2 24.0 22.4 22.5
15.9 14.5 15.1 14.2    41.0 43.9 36.3 39.9

PROC GLM DATA=IN;
  CLASS A B;
  MODEL Y = A B A*B / SS3 ;
  MEANS A B A*B;
TITLE 'A 2**2 DESIGN EXAMPLE';
TITLE2 'FIXED EFFECTS ANALYSIS';

PROC GLM DATA=IN;
     MODEL Y = X1 X2 X1*X2 / SS3 SOLUTION;
     OUTPUT OUT=DIAG R=RESID P=PRED;
TITLE2 'REGRESSION MODEL ANALYSIS';

RUN;