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