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;