* refinery.sas, refinery experiment, Table 15.54 (page 541);
options ls=74;
;
DATA REFINERY;
  INPUT A B C D Y;
  AB=MOD(A+B,3); AB2=MOD(A+2*B,3);
  AC=MOD(A+C,3); AC2=MOD(A+2*C,3);
  AD=MOD(A+D,3); AD2=MOD(A+2*D,3);
  BC=MOD(B+C,3); BC2=MOD(B+2*C,3);
  BD=MOD(B+D,3); BD2=MOD(B+2*D,3);
  CD=MOD(C+D,3); CD2=MOD(C+2*D,3);
  LINES;
  0 0 0 0      4.2
  0 0 1 2      5.9
  0 0 2 1      8.2
  0 1 0 2     13.1
  0 1 1 1     16.4
  0 1 2 0     30.7
  0 2 0 1      9.5
  0 2 1 0     22.2
  0 2 2 2     31.0
  1 0 0 2      7.7
  1 0 1 1     16.5
  1 0 2 0     14.3
  1 1 0 1     11.0
  1 1 1 0     29.0
  1 1 2 2     55.0
  1 2 0 0      8.5
  1 2 1 2     37.4
  1 2 2 1     66.3
  2 0 0 1     11.4
  2 0 1 0     21.1
  2 0 2 2     57.9
  2 1 0 0     13.5
  2 1 1 2     51.6
  2 1 2 1     76.5
  2 2 0 2     31.0
  2 2 1 1     74.5
  2 2 2 0     85.1 
;
PROC PRINT;
PROC GLM;
  CLASS A B C D AB AB2 AC AC2 AD AD2 BC BC2 BD BD2 CD CD2;
  MODEL Y = A B C D AB AB2 AC AC2 AD AD2 BC BC2 BD BD2 CD CD2;
  CONTRAST 'A LIN'  A -1  0  1;
  CONTRAST 'A QUAD' A  1 -2  1;
  CONTRAST 'B LIN'  B -1  0  1;
  CONTRAST 'B QUAD' B  1 -2  1;
  CONTRAST 'C LIN'  C -1  0  1;
  CONTRAST 'C QUAD' C  1 -2  1;
  CONTRAST 'D LIN'  D -1  0  1;
  CONTRAST 'D QUAD' D  1 -2  1;
PROC GLM;  CLASS A B C D;
  MODEL Y = A B C D A*B A*C A*D B*C B*D C*D;