Engine Mapping Problem

 /****************************************************************/
 /*          S A S   S A M P L E   L I B R A R Y                 */
 /*                                                              */
 /*    NAME: OPTEX4                                              */
 /*   TITLE: Engine Mapping Problem                              */
 /* PRODUCT: QC                                                  */
 /*  SYSTEM: ALL                                                 */
 /*    KEYS: Optimal Designs,                                    */
 /*   PROCS: FREQ OPTEX PLAN                                     */
 /*    DATA:                                                     */
 /*                                                              */
 /*     REF: SAS/QC Software:  Usage and Reference, Version 6,   */
 /*          First Edition, Volume 1 and Volume 2                */
 /*                                                              */
 /*          Vance, L.C. (1986), "Computer Construction of       */
 /*          Experimental Designs," General Motors Research      */
 /*          Report GMR-5411, General Motors Laboratories,       */
 /*          Warren, Michigan.                                   */
 /*    MISC:                                                     */
 /*                                                              */
 /****************************************************************/

options ps=60;
proc plan seed=46062;
   factors af=4 ordered egr=6 ordered sa=8 ordered / noprint;
   output out=a af  nvals=(15,16,17,18)
                egr nvals=(.020,.177,.377,.566,.921,1.117)
                sa  nvals=(10,16,22,28,34,40,46,52);

title1 "Engine Mapping Problem:";
title2 "First Attempt With The DETMAX Method";
proc optex seed=61552;
   model af|egr|sa@2 af*af egr*egr sa*sa;
   generate n=50 method=detmax;
run;

title1 "Engine Mapping Problem:";
title2 "First Attempt With The Sequential Algorithm";
proc optex seed=33805;
   model af|egr|sa@2 af*af egr*egr sa*sa;
   generate n=50 method=sequential;
run;

title1 "Engine Mapping Problem:";
title2 "Attempt To Improve On The Sequential Design?";
proc optex seed=33805;
   model af|egr|sa@2 af*af egr*egr sa*sa;
   generate n=50 method=detmax initdesign=sequential;
run;

title1 "Engine Mapping Problem:";
title2 "Initializing With A Replicated Saturated Design";
proc optex seed=33805;
    model af|egr|sa@2 af*af egr*egr sa*sa;
    generate n=10 method=detmax;
    output out=b;

data c; set b; do i=1 to 5; output; end; drop i;

proc optex data=a seed=33805;
    model af|egr|sa@2 af*af egr*egr sa*sa;
    generate n=50 method=detmax initdesign=c;
run;

title1 "Engine Mapping Problem:";
title2 "Searching for an Optimal 25-point Design";
proc optex data=a seed=36926;
    model af|egr|sa@2 af*af egr*egr sa*sa;
    generate n=25 method=detmax;
    output out=b;
run;

title2 "Augmenting the Best 25-point Design";
proc optex data=a seed=37034;
    model af|egr|sa@2 af*af egr*egr sa*sa;
    generate n=50 method=detmax augment=b;
run;

title1 "Engine Mapping Problem:";
title2 "The Sequential Design";
proc optex data=a seed=33805 noprint;
   model af|egr|sa@2 af*af egr*egr sa*sa;
   generate n=50 method=sequential;
   output out=b;
proc freq;
   table af egr sa / nocum;
run;

proc plan seed=46145;
   factors af=4 ordered egr=4 ordered sa=4 ordered / noprint;
   output out=a af  nvals=(  15,  16,  17,   18)
                egr nvals=(.020,.377,.566,1.117)
                sa  nvals=(  10,  28,  34,   52);

title1 "Engine Mapping Problem:";
title2 "A Smaller Candidate Set";
proc optex seed=61552;
   model af|egr|sa@2 af*af egr*egr sa*sa;
   generate n=50 method=detmax;
run;

title;