vp.full = facDesign(k = 3)                  #returns a 2^3 full factorial design
response(vp.full) = rnorm(2^3)              #generate some random response
summary(vp.full)                            #summary of the full factorial design (especially no defining relation)
#------------
vp.rep = facDesign(k = 2, replicates = 3, centerCube = 4)  #returns a full factorial design with 3 replications per factor combination and 4 center points
names(vp.rep) = c("Name 1", "Name 2") #set names
units(vp.rep) = c("min", "F")        #set units
lows(vp.rep) = c(20, 40, 60)          #set low and high factor values
highs(vp.rep) = c(40, 60, 80)         
summary(vp.rep)                       #summary of the replicated full factorial DesignRun the code above in your browser using DataLab