Create an orthogonal array using the Bose algorithm.

The bose program produces OA( q^2, k, q, 2 ), k <= q+1 for prime powers q.

createBose(q, ncol, bRandom = TRUE)

Arguments

q

the number of symbols in the array

ncol

number of parameters or columns

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 (1938) Sankhya Vol 3 pp 323-338

See also

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

Examples

A <- createBose(3, 3, FALSE)
B <- createBose(5, 4, TRUE)