* Constructing Central Composite Designs in SAS;

/*********************************************
The following code will list available central
composite designs in k factors, where 2<=k<=8:

%adxgen;
%adxff;
%adxcc;
%adxpcc(k);

Try this code for 5 factors:
*********************************************/

%adxgen;
%adxff;
%adxcc;
%adxpcc(5);

/*********************************************
Here is the output (found in the output window).
Compare this to what you find in Table 7.3, p. 307 of the text:


>     Number of
>    Runs in the
>     Factorial        Number of       Axial         Total Number
>      Portion       Center Points    Extreme          of runs
>    -----------   ----------------   -------   --------------------
>1.       16       10                 2.0000     36
>2.       16        7 = ( 1*6) +  1   2.0000     33 = ( 1* 22) +  11
>3.       32       16                 2.3784     58
>4.       32       12 = ( 4*2) +  4   2.3664     54 = ( 4* 10) +  14
>
>
>                      %adxccd() parameters to construct:
>                  -----------------------------------------
>              1.  %adxccd(*data set name*,5,16,10,2.0000)
>              2.  %adxccd(*data set name*,5,16,6/1,2.0000,2)
>              3.  %adxccd(*data set name*,5,32,16,2.3784)
>              4.  %adxccd(*data set name*,5,32,2/4,2.3664,5)
>
>
>  For blocked designs, equations give
>
>               Number of   Number in each      Number in
>     Total = ( factorial *   factorial   )  +   axial
>                blocks         block            block
>


The part of the output labeled "%adxccd() parameters to construct:"
is SAS code to construct the indicated CCD and output the design to
the data set *data set name*. Let's try design 2, a 2^(5-1) run in two
blocks, and see what we get.
*********************************************/
%adxccd(design2,5,16,6/1,2.0000,2);
proc print;run;