## 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;

```