Uses the BLKL algorithm to construct designs for estimating response surfaces.
|Printed output required (
||Model to be fitted when the design is used; no default i.e. this option must be specified|
||How to treat the constant in the model (
||Limit for expansion of terms in the model; default 2|
||Number of units (or trials) in the design|
||Number of design points to consider for deletion; default takes
||Number of design points to consider for inclusion; default takes
||Number of times to run the algorithm; default 100|
||Maximum amount by which to perturb the design points in the adjustment algorithm; default
||Number of blocks; default 1 i.e. design not blocked|
||Saves the block factor (if any) for the design|
||Number of units in each block of the design|
||Supplies values of the blocking factor for any previous experiments that are to be included in the analysis of the results of the design|
||Lists any variates that are part of a mixture (their values must be greater than zero and sum to one)|
||Seed for random numbers used to construct the initial design; default 124195|
||Saves the determinant of the information matrix for the best design|
||Saves the mean value of the standardized variance of predictions obtained from the design over a grid of x-values|
||Saves the maximum value of the standardized variance of predictions obtained from the design over a grid of x-values|
||Number of grid points in each x-direction to use for
||Lists the variates to be investigated in the design; these need not be supplied if none of the other parameters are required|
||Lists identifiers to be used to represent squares of the x-variates in the model|
||Lists identifiers to be used to represent cubes of the x-variates in the model|
||Support points for each x-variate in the design; if these are not (all) specified, they are formed automatically|
||Supplies values of the x-variates for any previous experiments that are to be included in the analysis of the results of the design|
AFRESPONSESURFACE uses the BLKL algorithm of Atkinson & Donev (1992) to construct a design to estimate parameters of a response-surface model. The algorithm searches for a D-optimal design: that is, a design that will provide a maximum value for the determinant of the information matrix of the model parameters. The model is specified using the
TERMS option, with the
CONSTANT option indicating whether or not it is to contain the constant term (or intercept). The
FACTORIAL sets a limit on the number of variates in each model term; by default this is 2.
NUNITS option specifies the number of units in the design. If there is to be a blocking factor in the design, the
NBLOCKS option specifies its number of levels, and the
BLOCKFACTOR option saves its values. The
BLOCKSIZE option specifies the number of units to be contained in each block of the design, in a scalar (if they are all the same) or a variate. If the block sizes are fixed, the specified sizes must sum to the number of units. However, if you specify sizes that sum to a value greater that the required number of units, the algorithm will search for the optimum block sizes.
When the model is to contain squares or cubes of x-variables, you will need to specify identifiers to represent these using the parameters of the directive. (When using regression directives such as
FIT to fit the model, you can use the
POL function but this is not recognised by
AFRESPONSESURFACE.) The x-variates in the model must then all be listed by the
X parameter. The corresponding squares are listed by the
X2 parameter, and the cubes by the
After specifying the
X parameter, you can also use the
SUPPORTPOINTS parameter to specify the x-values of the points to be considered when constructing the design; if this is not specified, these support points are formed automatically. Note that the variates are all assumed to be scaled to have values between -1 and 1. However, the criterion for D-optimality is unaffected by linear transformations of the X-variables. So you can rescale afterwards in any way you like.
AFRESPONSESURFACE allows for a set of mixture variates, whose values must all be positive and which must sum to 1. The variates in the mixture are specified using the
PREVIOUSVALUES parameter can be used to supply values of the x-variates for any previous experiments that are to be included in the analysis of the results of the new experiment, or to specify points that must be included in the design. The
PREVIOUSBLOCKS option should then indicate the blocks to which these previous observations belonged.
The BLKL algorithm starts by forming an initial design by making a random selection of points from the set of support points. The
SEED option defines the seed for the random numbers used to make the selection (default 124195). The algorithm then uses an exchanges algorithm to improve the design. At each exchange, the K points with the lowest variance of prediction amongst the points of design are considered for replacement by the L points with the highest variance of prediction amongst the candidate points for inclusion in the design. The algorithm makes the best one of these exchanges, continuing until there are none that increase the determinant. The values for K and L are specified by the
NINCLUSION options respectively. The best values depend on the design parameters, including the number of model parameters and the number of residual degrees of freedom. If they are unset,
AFRESPONSESURFACE sets them to the number of units divided by 4, or 4 if this larger. The
NRUNS option can be set to request that the algorithm is run several times, with different starting designs; the default is 100. The design parameters are saved only for the best design found, but you can set option
PRINT=monitoring to print information about each attempt.
There is also be a final adjustment algorithm which can be used except when the design contains mixtures. This examines the design points one at a time to see whether the design can be improved by moving it a small amount along any x-axis. If an increase is possible, the point providing the greatest increase is moved. The process is then repeated until no improvment is possible. This phase is selected by setting the
ADJUSTMENTSTEP option to the maximum amount (e.g. 0.2) by which the point may be moved on any axis.
DETERMINANT option allows you to saves the determinant of the information matrix for the best design. An alternative way of evaluating the design is to examine the standardized variance of the predictions that would be obtained from the design at other points, not in the design. The
MEANGRID option can save the mean value of the standardized variance of prediction over a grid of x-values, and the
MAXGRID option can save the maximum value. Number of grid points in each x-direction is specified by the
NGRIDPOINTSMETHOD option (default 5).
The algorithm is described in Sections 15.6 and 15.7 of Atkinson & Donev (1992). The source code of the algorithm was provided by Alex Donev. This, and his assistance generally with this Genstat implementation, is gratefully acknowledged.
Restrictions on the
PREVIOUSVALUES parameters are ignored.
Atkinson, A.C. & Donev, A.N. (1992) Optimum Experimental Designs. Oxford University Press, Oxford.
" Examples 2:4.9.14a-c " AFRESPONSESURFACE [NDELETION=10; NINCLUSION=40; NRUNS=1000; NUNITS=17;\ TERMS=X1 * X2 * X3 * X4 + X1_2 + X2_2 + X3_2 + X4_2;\ CONSTANT=estimate; DETERMINANT=Det; MAXGRID=Dmax; MEANGRID=Dave]\ X=X1,X2,X3,X4; X2=X1_2,X2_2,X3_2,X4_2 PRINT X1,X2,X3,X4; FIELD=8; DECIMALS=3 & Det,Dmax,Dave VARIATE [NVALUES=15] S1,S2,S3 READ [PRINT=data,errors] S1,S2,S3 -0.5774 -0.5774 -0.5774 0.5774 -0.5774 -0.5774 -0.5774 0.5774 -0.5774 0.5774 0.5774 -0.5774 -0.5774 -0.5774 0.5774 0.5774 -0.5774 0.5774 -0.5774 0.5774 0.5774 0.5774 0.5774 0.5774 1.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000 : VARIATE [NVALUES=6] P1,P2,P3 READ [PRINT=data,errors] P1,P2,P3 1.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 -1.0000 : AFRESPONSESURFACE [NUNITS=16; NDELETION=3; NINCLUSION=3;\ TERMS=X1 * X2 * X3 + X1_2 + X2_2 + X3_2; CONSTANT=estimate;\ DETERMINANT=Det; MAXGRID=Dmax; MEANGRID=Dave] \ X=X1,X2,X3; X2=X1_2,X2_2,X3_2; SUPPORTPOINTS=S1,S2,S3 SORT [INDEX=X1,X2,X3] X1,X2,X3 PRINT X1,X2,X3; FIELD=8; DECIMALS=3 & Det,Dmax,Dave AFRESPONSESURFACE [NUNITS=12; NDELETION=3; NINCLUSION=3; NRUNS=500; \ TERMS=X1 + X2 + X3 + X4 + X1.X2 + X1.X3 + X2.X3 + X4_2; CONSTANT=omit; MIXTURE=X1,X2,X3;\ DETERMINANT=Det; MAXGRID=Dmax; MEANGRID=Dave]\ X=X1,X2,X3,X4; X2=*,*,*,X4_2 PRINT X1,X2,X3,X4; FIELD=8; DECIMALS=3 & Det,Dmax,Dave