****************************************************************;
*** D-EFFICIENT PROC OPTEX DESIGN -- K=3, 4 CENTER PTS, N=22 ***;
****************************************************************;

OPTION PS=60 LS=64 NODATE NONUMBER;

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

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3|X4@2 X1*X1 X2*X2 X3*X3 X4*X4;
  GENERATE METHOD=FEDOROV CRITERION=D ITER=32 AUGMENT=MUST N=22;
  OUTPUT OUT=DSGN1;
TITLE 'D-EFFICIENT PROC OPTEX DESIGN -- K=3, 4 CENTER PTS, N=22';

PROC PRINT DATA=DSGN1;

DATA MUST; INPUT X1 X2 X3 X4 @@; CARDS;
0 0 0 0   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;
     DO X4=-1 TO 1 BY 1; OUTPUT; END; END; END; END;
TITLE;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3|X4@2 X1*X1 X2*X2 X3*X3 X4*X4;
  GENERATE METHOD=FEDOROV CRITERION=D AUGMENT=MUST N=22;
TITLE 'D-EFFICIENT PROC OPTEX DESIGN -- K=3, 4 CENTER PTS, N=22';
TITLE2 'NOTE THE CHANGE IN AVE. PREDICTION VARIANCE USING 3^K POINTS';
RUN;