Создать контрастную матрицу дизайна, используя автоматическую функцию
Я хотел бы создать функцию, которая возвращает матрицу контрастности дизайна факторов с количеством факторов, уровней каждого фактора в качестве входных данных. Я мой пример, у меня есть:
1-й фактор
treat <- gl(3, 20, labels = paste("t", 1:3, sep=""))
2-й фактор
imp <- rep(gl(2, 10, labels = c("yes", "no")), 3)
создать дизайн матрицы контрастов для "имп"
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]
Теперь я спрашиваю, как я могу начать создание матрицы Imp и Treat и всех комбинаций Imp и Treat, например, учитывая Коэффициент =2, Уровни = c(3,2) и Nsize=60?