# Full extended orthonormalised candidate set for two 4-level factors,
# full quadratic polynomial model as primary model and all three-order terms as potential.
K<-2; Levels <- rep(list(1:4),K)
cand.trt <- candidate_trt_set(Levels, K)
cand.full <- candidate_set_full(cand.trt, K)
prime.terms <- colnames(cand.full)[2:7]
poten.terms <- colnames(cand.full)[8:11]
Parameters <- c(1, rep(1,K), rep(1,K), K*(K-1)/2)
candidate_set_orth(cand.full, prime.terms, poten.terms)
Run the code above in your browser using DataLab