gpt4 book ai didi

r - 如何在r中的多重比较中将字符串变量传递给multcomp::glht函数

转载 作者:行者123 更新时间:2023-12-02 08:23:51 24 4
gpt4 key购买 nike

我正在进行单向方差分析和事后多重比较。以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)) ...但没有运气..非常感谢您的帮助。

最佳答案

我们可以使用do.call方法

aov1 <- do.call("aov", list(formula = as.formula(paste('mpg~', var)), data = quote(mtcars)))
out2 <- summary(multcomp::glht(aov1, linfct = do.call(mcp, setNames(list("Dunnet"), var))))

检查OP帖子中的输出

out1 <- summary(multcomp::glht(aov, linfct=mcp(cyl='Dunnet')))
all.equal(aov, aov1)
#[1] TRUE

all.equal(out1, out2)
#[1] TRUE
<小时/>

上面的内容可以封装在一个函数中

f1 <- function(dat, Var){
form1 <- formula(paste('mpg~', Var))
model <- aov(form1, data = dat)
model$call$formula <- eval(form1)
model$call$data <- substitute(dat)
summary(multcomp::glht(model, linfct = do.call(mcp, setNames(list("Dunnet"), Var))))

}

f1(mtcars, var)

关于r - 如何在r中的多重比较中将字符串变量传递给multcomp::glht函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/48124429/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com