The optimization is done with function NeighbourcalcUniversal.R
OSOAarbitrary(oa, el = 3, m = NULL, permlist = NULL, random = TRUE)an orthogonal array
matrix (preferred) or data.frame
must be 2 or 3. el=3: Li et al; el=2: Zhou and Tang
desired number of columns. Defaults to the maximum possible.
a list of length m of lists of length 1 of permutations of the levels
a logical (a random draw of permutations is used if permlist is NULL)
Li et al. Zhou and Tang