SOAs | R Documentation |
takes an OA(n,m,s,t) and creates an SOA(n,m',s^t',t') with t'<=t.
SOAs(
oa,
t = 3,
m = NULL,
noptim.rounds = 1,
noptim.repeats = 1,
optimize = TRUE,
dmethod = "manhattan",
p = 50
)
oa |
matrix or data.frame that contains an ingoing symmetric OA. Levels must be denoted as 0 to s-1 or as 1 to s. |
t |
the strength the SOA should have, can be 2, 3, 4, or 5. Must not
be larger than the strength of |
m |
the requested number of columns (see details for permitted numbers of columns) |
noptim.rounds |
the number of optimization rounds for each independent restart |
noptim.repeats |
the number of independent restarts of optimizations with |
optimize |
logical, default |
dmethod |
method for the calculation of |
p |
p for |
The resulting SOA will have at most m' columns in s^t levels and will be of
strength t. m'(m, t) is a function of the number of columns of oa
(denoted as m) and the strength t: m'(m,2)=m, m'(m,3)=m-1, m'(m,4)=floor(m/2),
m'(m,5)=floor((m-1)/2).
Suitable OAs for argument oa
can e.g. be constructed with OA creation functions
from package lhs or can be obtained from arrays listed in R package DoE.base.
matrix of class SOA
with the attributes that are listed below. All attributes can be accessed using function attributes
, or individual attributes can be accessed using function attr
. These are the attributes:
the type of array (SOA
or OSOA
)
character string that gives the strength
the phi_p value (smaller=better)
logical indicating whether optimization was applied
matrix that lists the id numbers of the permutations used
optional element, when optimization was conducted: the overall permutation list to which the numbers in permlist refer
the call that created the object
Ulrike Groemping
For full detail, see SOAs-package
.
He and Tang (2013)
Weng (2014)
aus <- SOAs(DoE.base::L27.3.4, optimize=FALSE) ## t=3 is the default
dim(aus)
soacheck2D(aus, s=3, el=3) ## check for 2*
soacheck3D(aus, s=3, el=3) ## check for 3
aus2 <- SOAs(DoE.base::L27.3.4, t=2, optimize=FALSE)
## t can be smaller than the array strength
## --> more columns with fewer levels each
dim(aus2)
soacheck2D(aus2, s=3, el=2, t=2) # check for 2
soacheck3D(aus2, s=3, el=2) # t=3 is the default (check for 3-)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.