createBoseBushl: Create an orthogonal array using the Bose-Bush algorithm with...

In bertcarnell/lhs: Latin Hypercube Samples

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

Description

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

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 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)

bertcarnell/lhs documentation built on July 4, 2024, 3:55 p.m.