1require(DoE.base) 2set.seed(9890) 3plan <- oa.design(nlevels=c(2,6,2), replications=2, repeat.only=TRUE) 4y <- rnorm(24) 5plan <- reptowide(add.response(plan,y)) 6aggregate(plan) 7 8plan <- param.design(oa.design(nlevels=c(2,6,2), factor.names=c("C1","C2","C3")), oa.design(ID=L4.2.3)) 9y <- rnorm(48) 10plan <- paramtowide(add.response(plan,y)) 11aggregate(plan) 12aggregate(plan,response="y") 13 14plan <- param.design(oa.design(nlevels=c(2,6,2), factor.names=c("C1","C2","C3")), oa.design(ID=L4.2.3)) 15plan <- paramtowide(plan) 16aggregate(plan) 17 18plan <- cross.design(oa.design(nlevels=c(2,6,2), factor.names=c("C1","C2","C3")), oa.design(ID=L4.2.3,replications=2,repeat.only=TRUE)) 19y <- rnorm(96) 20plan <- reptowide(add.response(plan,y)) 21aggregate(plan) 22aggregate(plan,response="y",FUN=sd) 23