cand.list = expand.grid(Factor1 = c("A", "B", "C", "D"), Factor2 = c("I", "II", "III"), Factor3 = c("Low", "High"), Factor4 = c("Yes", "No")) set.seed(69) optFederov( ~ ., data = cand.list, nTrials = 13)