vp.full = facDesign(k = 3) #returns a 2^3 full factorial design
vp.full = blocking(vp.full, 2) #design in 2 blocks
response(vp.full) = rnorm(2^3) #generate some random response
summary(vp.full) #summary of the full factorial design (especially no defining relation)
vp.frac = fracDesign(k = 4, gen = "D=ABC") #returns a 2^4-1 fractional factorial design. Factor D will be aliased with
response(vp.frac) = rnorm(2^(4-1)) #the three-way-interaction ABC (i.e. I = ABCD)
summary(vp.frac) #summary of the fractional factorial design
vp.rep = fracDesign(k = 3, replicates = 3, centerCube = 4) #returns a full factorial design with 3 replications per factor combination and 4 center points
summary(vp.rep) #summary of the replicated fractional factorial DesignRun the code above in your browser using DataLab