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 Design
Run the code above in your browser using DataLab