考虑dat1
在这里创建:
dat1 <- data.frame(Region = rep(c("r1","r2"), each = 100),
State = rep(c("NY","MA","FL","GA"), each = 10),
Loc = rep(c("a","b","c","d","e","f","g","h"),each = 5),
ID = rep(c(1:10), each=2),
var1 = rnorm(200),
var2 = rnorm(200),
var3 = rnorm(200),
var4 = rnorm(200),
var5 = rnorm(200))
我想编写一个函数,允许我指定数据、分组变量(Region
、State
或Loc
)和变量的行名(var1
: var5
),我希望它使用这些变量进行 manova 并以整洁的格式返回结果。运行该函数将如下所示:manova.test(dat = dat1, groupvar = "Region", cols = 5:9)
该函数看起来像这样(不起作用):
manova.test <- function(dat, groupvar, cols){
var.mat <- as.matrix(dat[, cols])
group.man <- manova(cbind(var.mat) ~ groupvar, data = dat)
summary(group.man)
}
你怎么做到这一点?我对如何以这种格式指定公式感到特别困惑!