ALPHA: Program for Constructing α-designs

  1. Introduction
  2. Using ALPHA
  3. Output
  4. Examples
  5. References

Introduction

Patterson & Williams (1976a) introduced a new class of resolvable incomplete block design called α-designs. Since their introduction, α-designs have become popular among designers of experiments due to two main reasons:

The construction of an α-design of size (r,k,s) begins with an r×k array α with elements in set of residues mod s. Each row of α is used to generate s-1 further rows by cyclic substitution. Denote the intermediate rs×k array by α*. Now add s to the second column of α*, 2s to the third column of α* and so on. The rows of the resulting array are blocks of the α-design. Each set of rows generated from the same row of α forms a complete replication. An α-design with pairs of treatments appearing in blocks either 0 or 1 times is referred to as α(0,1)-design and either 0, 1 or 2 times is referred to as α(0,1,2)-design. Chapter 4 of John & Williams (1995) gives an excellent summary of resolvable incomplete block designs including α-designs. In this note this book is abbreviated as JW.

ALPHA is a Gendex module for constructing optimal or near-optimal α-design. ALPHA handles up to 10,000 treatments. Designs constructed by ALPHA can be used as column components of resolvable row-column designs (see Nguyen & Williams 1993). ALPHA uses a 2-stage optimization process. Each stage of the optimization process uses an algorithm similar to the cyclic-coordinate exchange algorithm described in Nguyen (2002) . The algorithm uses extensively some theoretical results described in Patterson & Williams (1976b).

Using ALPHA

Let's assume all Gendex class files are in the directory c:\gendex and suppose you want to construct an α-design of size (r,k,s)=(2,4,5) (Example 1). At the working directory, type the following command at the Command Prompt (case is important):

java -cp c:\gendex ALPHA

The ALPHA window will pop up. Enter the the values 2, 4 and 5 in the (r,k,s) combination fields, the ALPHA window will become:

Now click START, ALPHA will start running and after try 1, the plan of the constructed design for this try pops up in the ALPHA output window (as the ratio E/U reaches 1) and ALPHA stops:

The START button has changed to RESET button. If you click the RESET button, the output will disappear and you can use ALPHA for a new design problem. Note that the default random seed is the one obtained from the system clock and the default number of tries is 1000. You can change these default values if you wish to.

Output

The result of the best try is displayed in the ALPHA output window and is also saved in the file ALPHA.htm in the working directory. This file can be read by a browser such as IE or Google Chrome. Information for this try includes:

  1. Try number;
  2. The number of iterations;
  3. E, the efficiency factor of this design;
  4. The ratio E/U where U is the upper bound of an IBD. U=min(UJ,UWP) for resolvable IBDs. UWP is the bound of Williams & Patterson (1977) good for any resolvable IBDs with v≥b. UJ is the bound of Jarrett (1989) good for any regular graph design (RGD). An RGD is an IBD with concurrences differing by at most 1. The program automatically stops if this ratio reaches 1.
  5. The distribution of the concurrences of this design;
  6. The design plan and the associated random seed;
  7. The time in seconds ALPHA used to construct this design;

An additional output of a ALPHA session is the file form.htm. This file contains the randomized plan of the constructed design.

Examples

  1. An α-design of size (r,k,s)=(3,4,3) (http://designcomputing.net/gendex/alpha/a1.html).
  2. An α-design of size (r,k,s)=(4,4,6) (http://designcomputing.net/gendex/alpha/a2.html).
  3. An α-design of size (r,k,s)=(2,3,4) (http://designcomputing.net/gendex/alpha/a3.html).
  4. An α-design of size (r,k,s)=(2,9,6) (http://designcomputing.net/gendex/alpha/a4.html).
  5. An α-design of size (r,k,s)=(3,5,8) (http://designcomputing.net/gendex/alpha/a5.html).
  6. An α-design of size (r,k,s)=(3,6,4) (http://designcomputing.net/gendex/alpha/a6.html).
  7. An α-design of size (r,k,s)=(2,30,33) (http://designcomputing.net/gendex/alpha/a7.html).

Notes:

References

Jarrett, R.G. (1989). A review of bounds for the efficiency factor of block designs. Austral. J. Statist. 31, 118-129.
John, J.A. & Williams E.R. (1995). Cyclic designs and computer-generated designs. New York: Chapman & Hall.
Nguyen, N-K (2002) A modified cyclic-coordinate exchange algorithm as illustrated by the construction of minimum-point second-order designs. Advances in Statistics, Combinatorics and Related Areas. Edited by C. Gulati, Y-X Lin, S. Mishra, J. Rayner. World Scientific Publishing Co. Pty. Ltd., 205-210.
Nguyen, N-K. and Williams, E.R. (1993). An algorithm for constructing optimal resolvable row-column designs. Austral. J. Statist. 35, 363-370.
Patterson, H.D. & Williams, E.R. (1976a). A new class of resolvable incomplete block designs. Biometrika 63, 83-92.
Patterson, H.D. & Williams, E.R. (1976b) Some theoretical results on general block designs. In Proceedings of the 5th British Combinatorial Conference. Congressus Numeratium XV, 489-496, Utilitas Mathematica, Winnipeg.
Williams & Patterson (1977) Upper bound for efficiency factors in block designs. Austral. J. Statist. 19, 194-201.

©2000-2014 Design Computing