IBD: Program for Constructing Optimal and Near-Optimal Incomplete Block Designs

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

Introduction

An incomplete block design (IBD) of size (v,k,r) is an arrangement of v treatments set out in blocks of size k (≤v) such that each treatment is replicated r times. We confine to binary IBDs, i.e. IBDs in which no treatment occurs more than once in a block. An IBD is said to be a balanced IBD, or BIBD, if every pair of treatments occurs in exactly λ blocks. IBDs in Examples 2, 4 and 9 are examples of BIBDs. An IBD is said to be t-resolvable if its blocks can be divided into subsets and each subset is an IBD of size (v,k,t). IBDs in Examples 1, 6, 7, 8 and 9 are examples of 1-resolvable IBDs or resolvable IBDs. IBDs in Examples 2 and 10 are examples of 2-resolvable IBDs. The books of John (1980) and John & Williams (1995) are recommended references on IBDs.

IBD is a Gendex module for constructing optimal or near-optimal IBDs (resolvable and non-resolvable). The optimality criterion and the algorithm used in bib is discussed in Nguyen (1993, 1994). Designs produced by bib are comparable to α-designs of Patterson & Williams (1976) and generalized cyclic designs of Hall & Jarrett (1981) in terms of the efficiency factor E of the design. These designs, in turn, can be used as column components of row-column designs (Nguyen & Williams 1993; Nguyen 1997).

Using IBD

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

java -cp c:\gendex ibd

The IBD window will pop up. Enter the the number of (v,k,r)=(6,2,4), the IBD window will become:

Now click START, since the combination (v,k,r)=(6,2,4) allows resolvable IBD, the following window will pop up:

Click Yes, the IBD window will become:

and the OUTPUT window showing the constructed IBD of size (v,k,r)=(6,2,4) will pop up. Note that the START button has been changed to the STOP one. If you close the pop-up window, the STOP button will become a RESET one. If you click this RESET button, the output will disappear and you can use IBD for a new design problem. Also note that the default random seed is the one obtained from the system clock and the default number of tries is 100. You can change these default values if you wish to.

Output

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

  1. try number;
  2. the random seed used;
  3. the number of iterations;
  4. E, the efficiency factor of this design.
  5. the ratio E/U where U is the upper bound of an IBD. U=UJ for non-resolvable IBDs and 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). RGDs are IBDs with concurrences differing by at most 1. IBDs in Examples 1-6 and 8-10 are examples of RGDs. The program automatically stops if this ratio reaches 1.
  6. the distribution of the concurrences of this design;
  7. the design plan;
  8. the time in seconds bib used to construct this design.

An additional output of a IBD session is the file form.htm. The following is the content of file form.htm of Example 1.

Rep   Block Plot  Treat
1     1     1     0
1     1     2     3
1     2     1     5
1     2     2     1
1     3     1     2
1     3     2     4
2     1     1     0
2     1     2     4
2     2     1     1
2     2     2     3
2     3     1     2
2     3     2     5
3     1     1     2
3     1     2     3
3     2     1     4
3     2     2     5
3     3     1     0
3     3     2     1
4     1     1     4
4     1     2     3
4     2     1     0
4     2     2     5
4     3     1     2
4     3     2     1

Examples

  1. An resolvable IBD of size (v,k,r)=(6,2,4) (http://designcomputing.net/gendex/ibd/b1.html).
  2. A 2-resolvable BIBD of size (v,k,r)=(6,4,10). (http://designcomputing.net/gendex/ibd/b2.html).
  3. An IBD of size (v,k,r)=(14,5,10) (http://designcomputing.net/gendex/ibd/b3.html).
  4. A BIBD of size (v,k,r)=(15,5,14) (http://designcomputing.net/gendex/ibd/b4.html).
  5. An IBD of size (v,k,r)=(60,9,3) (http://designcomputing.net/gendex/ibd/b5.html).
  6. A resolvable IBD of size (v,k,r)=(30,5,4) (http://designcomputing.net/gendex/ibd/b6.html).
  7. A resolvable IBD of size (v,k,r)=(36,6,4) (http://designcomputing.net/gendex/ibd/b7.html).
  8. A resolvable IBD of size (v,k,r)=(98,7,2) (http://designcomputing.net/gendex/ibd/b8.html).
  9. A resolvable BIBD of size (v,k,r)=(15,3,7) (http://designcomputing.net/gendex/ibd/b9.html).
  10. A 2-resolvable IBD of size (v,k,r)=(21,6,10) (http://designcomputing.net/gendex/ibd/b10.html).

Notes:

References

Box, G.E.P, Hunter, W.G. & Hunter, J.S. (1978) Statistics for experimenters. New York: John Wiley.
Clatworthy, W.H. (1973) Tables of two-associates-class partially balanced designs. Appl. Math. Ser. 63. National Bureau of Standards, Washington.
Hall, W.B. & Jarrett, R.G. (1981) Non-resolvable incomplete block designs with few replications. Biometrika 68, 617-627.
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. (1987) Cyclic designs and computer-generated designs. New York: Chapman & Hall.
John, P.W.M. (1980) Incomplete block designs. New York: Marcel Dekker, Inc.
Nguyen, N-K. (1993) An algorithm for constructing optimal resolvable block designs. Commun. Statist. B 22, 911-923.
Nguyen, N-K. (1994) Construction of optimal incomplete block designs by computer. Technometrics 36, 300-307.
Nguyen, N-K. and Williams, E.R. (1993) An algorithm for constructing optimal resolvable row-column designs. Austral. J. Statist. 35, 363-370.
Nguyen, N-K. (1997) Construction of optimal row-column designs by computer. Computing Science & Statistics 28, 471-475.
Patterson, H.D. & Williams, E.R. (1976) A class of resolvable incomplete block designs. Biometrika 63, 83-92.
Williams & Patterson (1977) Upper bound for efficiency factors in block designs. Austral. J. Statist. 19, 194-201.

©2008 Design Computing