Expansive replacement for two orthogonal arrays
an orthogonal array, must be a matrix;
an orthogonal array, must be a matrix
This function mainly is meant for combining two orthogonal arrays via what Kuhfeld (2009) calls expansive replacement.
array2 is a full factorial in two factors, argument
all = TRUE
creates a list of expanded arrays obtained by permuting the second array in all ways
that may lead to combinatorially different end results. With s1 and s2 the numbers
of levels of the factors in
array2, this is the number of partitions of the
array2 into s1 equally-sized groups, multiplied with s2!^(s1-1) for
the possibilities of permuting the levels of the second factor within all but the first
level of the first factor. This functionality is primarily meant for the creation of
strength 3 arrays in combination with arrays listed in the catalogue
(see an example on expanding the 6-level factor in
The function returns an object of classes
which can be used in function
oa.design, or a list of such arrays,
all = TRUE.
Whether or not the object is an orthogonal array depends on the choice of suitable input
arrays by the user. The properties of the resulting array(s) can e.g. be inspected with functions
This package is still under development. Bug reports and feature requests are welcome.
Kuhfeld, W. (2009). Orthogonal arrays. Website courtesy of SAS Institute http://support.sas.com/techsup/technote/ts723.html.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
myL126.96.36.199.1 <- expansive.replace(L188.8.131.52.1, L184.108.40.206.1) L220.127.116.11.1.4.2_list <- expansive.replace(L18.104.22.168.2.6.1, cbind(U=rep(c(1,2),each=3), V=rep(1:3,2)), all=TRUE) ## the list of 60 resolution IV arrays can be used for design creation, ## e.g. as follows: ## Not run: ## resolution IV designs obtained from the 60 different arrays deslist <- lapply(L22.214.171.124.1.4.2_list, function(aa) oa.design(aa, nlevels=c(2,2,2,2,3,4,4), columns="min34")) table(A4s <- sapply(deslist, length4)) ## a single best design exists best <- deslist[[which(A4s < 2)]] GWLP(best) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.