Create an orthogonal array using the Bose-Bush algorithm with alternate strength >= 3.

The bosebushl program produces OA( lambda*q^2, k, q, 2 ), k <= lambda*q+1, for prime powers q and lambda > 1. Both q and lambda must be powers of the same prime.

Usage

createBoseBushl(q, ncol, lambda, bRandom = TRUE)

Arguments

q

the number of symbols in the array

ncol

number of parameters or columns

lambda

the lambda of the BoseBush algorithm

bRandom

should the array be randomized

Value

an orthogonal array

Details

From Owen: An orthogonal array A is a matrix of n rows, k columns with every element being one of q symbols 0,...,q-1. The array has strength t if, in every n by t submatrix, the q^t possible distinct rows, all appear the same number of times. This number is the index of the array, commonly denoted lambda. Clearly, lambda*q^t=n. The notation for such an array is OA( n, k, q, t ).

References

Owen, Art. Orthogonal Arrays for: Computer Experiments, Visualizations, and Integration in high dimensions. https://lib.stat.cmu.edu/designs/oa.c. 1994 R.C. Bose and K.A. Bush (1952) Annals of Mathematical Statistics, Vol 23 pp 508-524.

See also

Other methods to create orthogonal arrays [createBoseBush()], [createBose()], [createBush()], [createAddelKemp()], [createAddelKemp3()], [createAddelKempN()], [createBusht()]

Examples

A <- createBoseBushl(3, 3, 3, TRUE)
B <- createBoseBushl(4, 4, 16, TRUE)