我正在做单向方差分析和事后多重比较。以 mtcars 数据集为例:
mtcars$cyl <- as.factor(mtcars$cyl)
aov<- aov(mpg~cyl, data=mtcars)
summary(multcomp::glht(aov, linfct=mcp(cyl='Dunnet')))
但是,我不想将变量硬编码为 cyl。所以我创建了一个变量 var='cyl':
var <- 'cyl'
aov <- aov(as.formula(paste('mpg~', var)), data=mtcars)
summary(multcomp::glht(aov, linfct=mcp( var='Dunnet')))
我收到错误消息:
Error in mcp2matrix(model, linfct = linfct) : Variable(s) ‘var’ have been specified in ‘linfct’ but cannot be found in ‘model’!
我认为问题来自在 mcp 函数中传递 var 。我怎样才能解决这个问题?我试过: as.name(var) , eval(quote(var))... 但没有运气.. 非常感谢您的帮助。