A Brief Gosset Tutorial


The first step in generating a design with Gosset is to write a program describing the experimental situation. Let's assume we want to generate an I-Optimal design for 8 factors with a full quadratic model. We also want to have three trials more than the minimum required. The Gosset program for this is,

10 range X1 X2 X3 X4 X5 X6 X7 X8 -1 1
20 model (1+X1+X2+X3+X4+X5+X6+X7+X8)^2

This simple program specifies the situation. Notice that each line requires a line number. The commands will be executed in numerical order. The command "range" specifies the low and high levels for continuous factors. In this case, all of the factors from X1 to X8 will range from -1 to 1.

The "model" command requires a little more explanation. Notice that if you multiply (1+X1+X2+X3+X4+X5+X6+X7+X8) *(1+X1+X2+X3+X4+X5+X6+X7+X8) you get 1+X1+X2+...+X8_X1X2+X1X3+...+X7X8+X1^2+X2^2+...+X8^2 - the full quadratic model with all 2-factor interactions, but with the b coefficients missing. Gosset knows you want the b coefficients. This is a shorthand for model specification. If you wish to subtract terms from the model, simply subtract them at the end of line 20 above.

Now that we have the program, we need to compile it. This is accomplished by simply typing "c" and pressing return. Gosset will provide a list of terms in your model. If the list looks like you expected, you are doing things right.

The next step is to compute the "moments matrix." This is a mathematical detail handled by Gosset. All we need to do is type

"m n=100000"

and press return. When Gosset is done, it will let us know.

Last, we have to generate the actual design. To do this, we type

"d n=100 extra=3 tiny=1e-5 type=I"

and press return. "n=100" tells Gosset to generate 100 designs and keep the best. "extra=3" tells Gosset to include 3 more trials in the design than the minimum required. "tiny=1e-5" tells Gosset the smallest difference to consider to be real. "type=I" tells Gosset to generate an I-Optimal design.

We will see a list of the progress as Gosset works. When it is finished, we type "interp" to see the design.