R/createOA.R
createAddelKemp.RdThe addelkemp program produces OA( 2q^2, k, q, 2 ), k <= 2q+1,
for odd prime powers q.
createAddelKemp(q, ncol, bRandom = TRUE)
| q | the number of symbols in the array |
|---|---|
| ncol | number of parameters or columns |
| bRandom | 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 dimenstions. http://lib.stat.cmu.edu/designs/oa.c. 1994 S. Addelman and O. Kempthorne (1961) Annals of Mathematical Statistics, Vol 32 pp 1167-1176.
Other methods to create orthogonal arrays [createBoseBush()], [createBose()], [createAddelKemp3()], [createAddelKempN()], [createBusht()], [createBoseBushl()]
A <- createAddelKemp(3, 3, TRUE) B <- createAddelKemp(3, 5, FALSE)