tests/NewKobi1.R

In planor: Generation of Regular Factorial Designs

library("planor")

#---------------------------------------------------------------------------
# EXAMPLES FROM THE PLANOR MANUAL
#---------------------------------------------------------------------------
# Exemple 1 page 7 : SEARCH FOR ALL SOLUTIONS
#---------------------------------------------------------------------------
cat("\n")
cat("***************** EXEMPLE 1 PAGE 7 *****************\n")
cat("\n")
cat("Four 3-level treatment factors and one 3-level block factor\n")
cat("Model: bloc+(A+B+C+D)^2\n")
cat("Estimate: A+B+C+D\n")
cat("N=3^3\n")
cat("\n")
# 
cat("*** RUN ***\n")
cat("\n")
F1 <- planor.factors( factors=c(LETTERS[1:4], "Bl"),
                     nlevels=rep(3,5),
                     block=~Bl )
M1 <- planor.model( model=~Bl+(A+B+C+D)^2, estimate=~A+B+C+D )
K1 <- planor.designkey(factors=F1, model=M1, nunits=3^3,
                       base=~A+B+C, max.sol=Inf, verbose=TRUE)

summary(K1[1])
alias(K1[1])
a <- K1[[1]];
print(class(a@.Data[[1]]))
print(alias(a@.Data[[1]]))
print(summary(a@.Data[[1]]))

summary(K1[34])
alias(K1[34])