Constrained Mixture Experiment

 /****************************************************************/
 /*          S A S   S A M P L E   L I B R A R Y                 */
 /*                                                              */
 /*    NAME: OPTEX5                                              */
 /*   TITLE: Constrained Mixture Experiment                      */
 /* PRODUCT: QC                                                  */
 /*  SYSTEM: ALL                                                 */
 /*    KEYS: Optimal Designs,                                    */
 /*   PROCS: OPTEX                                               */
 /*    DATA:                                                     */
 /*                                                              */
 /*     REF: SAS/QC Software:  Examples                          */
 /*    MISC: ADX Macros are stored in the AUTOCALL library       */
 /*                                                              */
 /****************************************************************/

options ps=60;

title "A Design for a Constrained Mixture Experiment";

%adxgen;
%adxmix;
%adxinit

%adxxvert(lube,add .07-.18 /compa -.3 /compb .37-.7 /compc -.15)
proc print data=lube; run;

proc optex data=lube seed=54321;
   model add|compa|compb@2 add*add compa*compa compb*compb;
   generate n=18;
   id dimen dist;
   output out=optex;

proc sort; by dimen add compa compb;
proc print; run;

* Construct Snee's design for comparison;
data snee; set lube;
   if    ( ( dimen = 0 )
       | ( ( dimen = 1 ) & ( dist > .2 ) )
       | ( ( dimen = 2 ) & ( dist > .1 ) )
       |   ( dimen = 3 )                 );
proc print;
run;

proc optex data=lube seed=27180;
   model add|compa|compb@2 add*add compa*compa compb*compb;
   generate method=sequential initdesign=snee;
run;

title;