************************************************;
*** MIXTURE DESIGN -- QUADRATIC MODEL -- N=6 ***;
************************************************;

DATA IN;
     DO X1 = 0 TO       1 BY .1;
     DO X2 = 0 TO    1-X1 BY .1;
     DO X3 = 0 TO 1-X1-X2 BY .1;
        IF .999 < SUM(OF X1-X3) < 1.001 THEN OUTPUT;
     END; END; END;
DATA OTHER;
     X1=1/3; X2=1/3; X3=1/3; OUTPUT;  *** CENTROID ***;

     X1=1/4; X2=1/4; X3=1/2; OUTPUT;
     X1=1/4; X2=1/2; X3=1/4; OUTPUT;
     X1=1/2; X2=1/4; X3=1/4; OUTPUT;

     X1=2/3; X2=1/3; X3=0  ; OUTPUT;
     X1=2/3; X2=0  ; X3=1/3; OUTPUT;
     X1=1/3; X2=2/3; X3=0  ; OUTPUT;
     X1=1/3; X2=0  ; X3=2/3; OUTPUT;
     X1=0  ; X2=1/3; X3=2/3; OUTPUT;
     X1=0  ; X2=2/3; X3=1/3; OUTPUT;
DATA IN; SET IN OTHER;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3@2 / NOINT;
  GENERATE METHOD=DETMAX N=6;
TITLE 'MIXTURE DESIGN -- QUADRATIC MODEL -- N=6';
RUN;

*************************************************;
*** MIXTURE DESIGN -- QUADRATIC MODEL -- N=10 ***;
*************************************************;

DATA IN;
     DO X1 = 0 TO       1 BY .1;
     DO X2 = 0 TO    1-X1 BY .1;
     DO X3 = 0 TO 1-X1-X2 BY .1;
        IF .999<SUM(OF X1-X3)<1.001  THEN OUTPUT;
     END; END; END;
DATA OTHER;
     X1=1/3; X2=1/3; X3=1/3; OUTPUT;  *** CENTROID ***;

     X1=1/4; X2=1/4; X3=1/2; OUTPUT;
     X1=1/4; X2=1/2; X3=1/4; OUTPUT;
     X1=1/2; X2=1/4; X3=1/4; OUTPUT;

     X1=2/3; X2=1/3; X3=0  ; OUTPUT;
     X1=2/3; X2=0  ; X3=1/3; OUTPUT;
     X1=1/3; X2=2/3; X3=0  ; OUTPUT;
     X1=1/3; X2=0  ; X3=2/3; OUTPUT;
     X1=0  ; X2=1/3; X3=2/3; OUTPUT;
     X1=0  ; X2=2/3; X3=1/3; OUTPUT;
DATA IN; SET IN OTHER;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3@2 / NOINT;
  GENERATE METHOD=DETMAX N=10;
TITLE 'MIXTURE DESIGN -- QUADRATIC MODEL -- N=10';
RUN;

****************************************************;
*** MIXTURE DESIGN -- SPECIAL CUBIC MODEL -- N=7 ***;
****************************************************;

DATA IN;
     DO X1 = 0 TO       1 BY .1;
     DO X2 = 0 TO    1-X1 BY .1;
     DO X3 = 0 TO 1-X1-X2 BY .1;
        IF .999<SUM(OF X1-X3)<1.001  THEN OUTPUT;
     END; END; END;
DATA OTHER;
     X1=1/3; X2=1/3; X3=1/3; OUTPUT;  *** CENTROID ***;

     X1=1/4; X2=1/4; X3=1/2; OUTPUT;
     X1=1/4; X2=1/2; X3=1/4; OUTPUT;
     X1=1/2; X2=1/4; X3=1/4; OUTPUT;

     X1=2/3; X2=1/3; X3=0  ; OUTPUT;
     X1=2/3; X2=0  ; X3=1/3; OUTPUT;
     X1=1/3; X2=2/3; X3=0  ; OUTPUT;
     X1=1/3; X2=0  ; X3=2/3; OUTPUT;
     X1=0  ; X2=1/3; X3=2/3; OUTPUT;
     X1=0  ; X2=2/3; X3=1/3; OUTPUT;
DATA IN; SET IN OTHER;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3 / NOINT;
  GENERATE METHOD=DETMAX N=7;
TITLE 'MIXTURE DESIGN -- SPECIAL CUBIC MODEL -- N=17;
RUN;

*****************************************************;
*** MIXTURE DESIGN -- SPECIAL CUBIC MODEL -- N=10 ***;
*****************************************************;

DATA IN;
     DO X1 = 0 TO       1 BY .1;
     DO X2 = 0 TO    1-X1 BY .1;
     DO X3 = 0 TO 1-X1-X2 BY .1;
        IF .999<SUM(OF X1-X3)<1.001  THEN OUTPUT;
     END; END; END;
DATA OTHER;
     X1=1/3; X2=1/3; X3=1/3; OUTPUT;  *** CENTROID ***;

     X1=1/4; X2=1/4; X3=1/2; OUTPUT;
     X1=1/4; X2=1/2; X3=1/4; OUTPUT;
     X1=1/2; X2=1/4; X3=1/4; OUTPUT;

     X1=2/3; X2=1/3; X3=0  ; OUTPUT;
     X1=2/3; X2=0  ; X3=1/3; OUTPUT;
     X1=1/3; X2=2/3; X3=0  ; OUTPUT;
     X1=1/3; X2=0  ; X3=2/3; OUTPUT;
     X1=0  ; X2=1/3; X3=2/3; OUTPUT;
     X1=0  ; X2=2/3; X3=1/3; OUTPUT;
DATA IN; SET IN OTHER;

PROC OPTEX DATA=IN;
  EXAMINE DESIGN;
  MODEL X1|X2|X3 / NOINT;
  GENERATE METHOD=DETMAX N=10;
TITLE 'MIXTURE DESIGN -- SPECIAL CUBIC MODEL -- N=10';

RUN;