****************************************************************;
*** GENERATE DESIGN FOR WELCH ET AL (1990) JQT ARTICLE MODEL ***;
****************************************************************;

OPTION PS=60 LS=64 NODATE NONUMBER;

TITLE 'JQT ARTICLE DESIGN';

DATA IN;
     DO W1=-1 TO 1 BY 1;
     DO W2=-1 TO 1 BY 1;
     DO W3=-1 TO 1 BY 1;
     DO W4=-1 TO 1 BY 1; 
     DO W5=-1 TO 1 BY 1; 
     DO W6=-1 TO 1 BY 1; 
     DO D = 1 TO 5 BY 1;
        OUTPUT; 
     END; END; END; END; END; END; END;

PROC OPTEX DATA=IN SEED=3029751;
  CLASS D;
  EXAMINE DESIGN;
  MODEL W1|W2|W3|W4@2 W5|W6 W1*W1 W2*W2 W3*W3 W4*W4 W5*W5 W6*W6
        D D*W1 D*W2 D*W3 D*W4 D*W5 D*W6;
  GENERATE METHOD=DETMAX N=60 ITER=50;
TITLE 'GENERATE DESIGN FOR JQT ARTICLE MODEL';

PROC OPTEX DATA=IN SEED=3029751;
  CLASS D;
  EXAMINE DESIGN;
  MODEL W1|W2|W3|W4@2 W5|W6 W1*W1 W2*W2 W3*W3 W4*W4 W5*W5 W6*W6
        D D*W1 D*W2 D*W3 D*W4 D*W5 D*W6 ;
  GENERATE METHOD=DETMAX N=60 ITER=50;
  OUTPUT OUT=GDSGN NUMBER = GBEST;
TITLE2 'BEST G-EFFICIENT DESIGN AMONG 50 BEST D-EFFICIENT DESIGNS';
PROC PRINT DATA=GDSGN;
RUN;