The bosebush
program
produces OA( 2q^2, k, q, 2 )
, k <= 2q+1
, for powers of 2, q=2^r
.
createBoseBush(q, ncol, bRandom = TRUE)
the number of symbols in the array
number of parameters or columns
should the array be randomized
an orthogonal array
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 )
.
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.
Other methods to create orthogonal arrays [createBush()], [createBose()], [createAddelKemp()], [createAddelKemp3()], [createAddelKempN()], [createBusht()], [createBoseBushl()]
A <- createBoseBush(4, 3, FALSE)
B <- createBoseBush(8, 3, TRUE)