R 使用自动功能创建设计对比矩阵

R 使用自动功能创建设计对比矩阵,r,functional-programming,compare-contrast,R,Functional Programming,Compare Contrast,我想创建一个函数,返回设计对比矩阵的因素与数量的因素,在每个因素的水平作为输入。以我为例,我有: 第一因素 treat查看?展开.grid请多给infominer一些评论 treat <- gl(3, 20, labels = paste("t", 1:3, sep="")) imp <- rep(gl(2, 10, labels = c("yes", "no")), 3) contrasts(imp) <- c(-1, 1) Imp <- model.matrix(

我想创建一个函数,返回设计对比矩阵的因素与数量的因素,在每个因素的水平作为输入。以我为例,我有:

第一因素
treat查看
?展开.grid
请多给infominer一些评论
treat <- gl(3, 20, labels = paste("t", 1:3, sep=""))
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]