tests/Kobi3.R In planor: Generation of Regular Factorial Designs

```library("planor")

#---------------------------------------------------------------------------
# EXAMPLES FROM THE PLANOR MANUAL
#---------------------------------------------------------------------------
# Exemple 3 page 16
#---------------------------------------------------------------------------
cat("\n")
cat("***************** EXEMPLE 3 PAGE 16 *****************\n")
cat("\n")
cat("Five treatment factors at 2 levels \n")
cat("robotFac \n")
cat("Block structure: 4 plates / (2 rows x 4 columns) \n")
cat("Model: nsoil + qsoil + Rug + conc + Tact \n")
cat("Estimate: idem \n")
cat("Hierarchy 1: conc constant in plate \n")
cat("Hierarchy 2: Tact constant in plate \n")
cat("Hierarchy 3: nsoil constant in plate x column \n")
cat("Hierarchy 4: qsoil constant in plate x column \n")
cat("N=32 units\n")
cat("\n")
#

cat("*** RUN ***\n")
# ptmtotal <- proc.time()
cat("\n")
robotFac <- planor.factors( factors=list(
conc=c(1,3),
Tact=c(15,30),
nsoil=c("curd","Saint-Paulin"),
qsoil=c("0.01g","0.10g"),
Rug=c(0.25,0.73),
plate=1:4,
row=1:2,
col=1:4),
block=~plate+row+col,
hierarchy=list(~conc/plate,
~Tact/plate,
~nsoil/(plate*col),
~qsoil/(plate*col)))
# Obligation for the moment of a complicated model formulation:
# - all factors in the design must appear in the model-estimate formulae
#   (it is not sufficient to declare them as basic factors)
robotMod <- planor.model( model=~nsoil * qsoil * Rug * conc * Tact,
listofmodels=list(c(~plate+row+col+Rug,~Rug)) )
# ptm <- proc.time()
robotKey <- planor.designkey(factors=robotFac, model=robotMod,
nunits=32,
base=~plate+row+col, max.sol=1)
# cat("TEMPS planor.designkey", proc.time()-ptm,"\n")
cat("\n")
# ptm <- proc.time()
robotDes <- planor.design(robotKey)
# cat("TEMPS planor.design", proc.time()-ptm,"\n")
print(robotDes@design)
cat("\n")
# Increase print amount:
options(planor.max.print=100)
# ptm <- proc.time()
summary(robotKey, show="dtbw")
# cat("TEMPS summary(pick)", proc.time()-ptm,"\n")
# cat("TEMPS total3A", proc.time()-ptmtotal,"\n")
cat("\n")
# --------------------------------------------------

cat("\n ROBOT1B \n")
cat("Block structure: 2 plates / (2 rows x 4 columns) \n")
cat("Model: nsoil + qsoil + Rug + conc \n")
cat("Estimate: idem \n")
cat("Hierarchy 1: conc constant in plate \n")
cat("Hierarchy 2: Tact constant in plate \n")
cat("Hierarchy 3: nsoil constant in plate x column \n")
cat("Hierarchy 4: qsoil constant in plate x column \n")
cat("N=16 units\n")
cat("\n")
#
cat("*** RUN ***\n")
# ptmtotal <- proc.time()
cat("\n")
F3B <- planor.factors( factors=list(
conc=c(1,3),
Tact=c(15,30),
nsoil=c("curd","Saint-Paulin"),
qsoil=c("0.01g","0.10g"),
Rug=c(0.25,0.73),
plate=1:4,
row=1:2,
col=1:4),
block=~plate+row+col,
hierarchy=list(~conc/plate,
~Tact/plate,
~nsoil/(plate*col),
~qsoil/(plate*col)))
# Obligation for the moment of a complicated model formulation:
# - all factors in the design must appear in the model-estimate formulae
#   (it is not sufficient to declare them as basic factors)
M3B <- planor.model( model=~nsoil * qsoil * Rug * conc * Tact,
listofmodels=list(c(~plate+row+col,~Rug)) )
# ptm <- proc.time()
K3B <- planor.designkey(factors=F3B, model=M3B, nunits=32,
base=~plate+row+col, max.sol=1)
# cat("TEMPS planor.designkey", proc.time()-ptm,"\n")
cat("\n")
# ptm <- proc.time()
P3B <- planor.design(key=K3B, select=1)
# cat("TEMPS planor.design", proc.time()-ptm,"\n")
# ptm <- proc.time()
resum3B <-summary(K3B[1], show="dtbw")
# cat("TEMPS summary(pick)", proc.time()-ptm,"\n")
# cat("TEMPS total", proc.time()-ptmtotal,"\n")
```

Try the planor package in your browser

Any scripts or data that you put into this service are public.

planor documentation built on March 19, 2020, 1:06 a.m.