因子の数、各因子のレベルを入力として因子の設計対比行列を返す関数を作成したいと思います。私の例では、私は持っています：
第 1 因子
treat <- gl(3, 20, labels = paste("t", 1:3, sep=""))
第 2 因子
imp <- rep(gl(2, 10, labels = c("yes", "no")), 3)
「imp」の対比の計画マトリックスを作成する
contrasts(imp) <- c(-1, 1)
Imp <- model.matrix(~ imp)[, -1]
「治療」の対比の計画行列を作成する
contrasts(treat) <- cbind(c(0,1,0),c(0,0,1))
Treat <- model.matrix(~ treat)[, -1]
次に、Factor=2、Levels = c(3,2)、Nsize=60 の場合、Imp と Treat マトリックスと Imp と Treat のすべての組み合わせを作成する方法を質問します。