createBusht: Create an orthogonal array using the Bush algorithm with...

In bertcarnell/lhs: Latin Hypercube Samples

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)

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