****************************************************;
*** SAS EFFICIENCIES FOR K=3, 4 CENTER POINT CCD ***;
****************************************************;

OPTION PS=60 LS=64 NODATE NONUMBER;

DATA MUST; INPUT X1 X2 X3 @@; CARDS;
 1 1 1    1 1 -1    1 -1 1    1 -1 -1
-1 1 1   -1 1 -1   -1 -1 1   -1 -1 -1
 1 0 0  -1 0 0   0 1 0   0 -1 0   0 0 1   0 0 -1
 0 0 0   0 0 0   0 0 0   0 0 0
;
DATA IN;
     DO X1=-1 TO 1 BY .1;
     DO X2=-1 TO 1 BY .1;
     DO X3=-1 TO 1 BY .1; OUTPUT; END; END; END;
TITLE;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3@2 X1*X1 X2*X2 X3*X3;
  GENERATE CRITERION=D AUGMENT=MUST N=18 ITER=1;
  OUTPUT OUT=DSGN1;
TITLE 'SAS EFFICIENCIES FOR K=3, 4 CENTER POINT CCD';
PROC PRINT DATA=DSGN1;
RUN;