| createBusht {lhs} | R Documentation |
Create an orthogonal array using the Bush algorithm with alternate strength.
Description
The busht program produces OA( q^t, k, q, t ), k <= q+1, t>=3,
for prime powers q.
Usage
createBusht(q, ncol, strength, bRandom = TRUE)
Arguments
q |
the number of symbols in the array |
ncol |
number of parameters or columns |
strength |
the strength of the array to be created |
bRandom |
should the array be randomized |
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 ).
Value
an orthogonal array
References
Owen, Art. Orthogonal Arrays for: Computer Experiments, Visualizations, and Integration in high dimensions. https://lib.stat.cmu.edu/designs/oa.c. 1994 K.A. Bush (1952) Annals of Mathematical Statistics, Vol 23 pp 426-434
See Also
Other methods to create orthogonal arrays [createBoseBush()], [createBose()], [createAddelKemp()], [createAddelKemp3()], [createAddelKempN()], [createBoseBushl()]
Examples
set.seed(1234)
A <- createBusht(3, 4, 2, TRUE)
B <- createBusht(3, 4, 3, FALSE)
G <- createBusht(3, 4, 3, TRUE)