4

我确定我在这里遗漏了一些明显的东西；我正在尝试使用循环将输入更改为公式，然后生成glht涉及该公式的公式。（glht=一般线性假设检验）。谢谢大家！

df2 <- data.frame(f1=seq(1:5),f2=seq(4,8), y=sample(x=c(0,1),size=5,replace=TRUE))
df2 <- data.frame(lapply(df2,factor)) # convert to factors

library(multcomp)

for (i in 1:(length(df2)-1)){
  fmla <- as.formula( paste("y~",colnames(df2)[i],sep=""))
  fm1 <- glm (fmla, family = binomial("logit"), data=df2)
 }

这很有效；现在我正在尝试动态传递f1以使以下行：

  tuk <- glht(fm1, linfct=mcp(f1 ="Tukey"))

我努力了：

tuk <- glht(fm1,linfct=mcp(colnames(df2)[i] ="Tukey"))

但这给出了：

Error: unexpected '=' in " tuk <- glht(fm1,linfct=mcp(colnames(df2)[i] ="

我认为这与parse/eval. 所以我尝试了：

e1 <- paste("fm1,linfct=mcp(",colnames(df2)[i],"=Tukey)",sep="")
glht(e1)

给出：

Error in object$coefficients : $ operator is invalid for atomic vectors
Error in modelparm.default(model, ...) :  no 'coef' method for 'model' found!

然后尝试：

e1 <- parse(text=paste("fm1,linfct=mcp(",colnames(df2)[i],"=Tukey)",sep="") )
glht(e1)

给出：

Error in parse.default(text = paste("fm1,linfct=mcp(", colnames(df2)[i],  : 
 <text>:1:4: unexpected ','
1: fm1,
      ^

最后：

e1 <- as.expression(paste("fm1,linfct=mcp(",colnames(df2)[i],"=Tukey)",sep="") )
glht(e1)

给出：

Error in UseMethod("vcov") : no applicable method for 'vcov' applied to an object of class "expression"
Error in modelparm.default(model, ...) :  no 'vcov' method for 'model' found!
4

1 回答 1

4

去年在 SO 上提出并回答了一个基本相同的问题。下面的答案使用了@hadley 然后概述的策略，更详细地解释了它的工作原理。

这个问题的棘手之处在于您需要将每个列名作为唯一参数的名称mcp()传递给. 使用do.call()可能是实现这一目标的最佳方式。

你想要什么：

mcp(f1 = "Tukey")

当你有一个变量时如何得到它X="f1"：

X <- "f1"
arg <- list("Tukey")
names(arg) <- X
do.call(mcp, arg)

在您的情况下，您可以执行以下操作：

require(multcomp)
df2 <- data.frame(f1=seq(1:5),f2=seq(4,8), y=sample(x=c(0,1),
                  size=5,replace=TRUE))
df2 <- data.frame(lapply(df2,factor)) # convert to factors

nms <- head(colnames(df2), -1)
lapply(nms, 
       function(X) {
           fmla <- as.formula( paste("y~", X, sep=""))
           fm1 <- glm (fmla, family = binomial("logit"), data=df2)
           ## From here on down is the solution
           args <- list("Tukey")
           names(args) <- X
           cmp <- do.call(mcp, args)
           glht(fm1, linfct = cmp)
       })
于 2012-08-16T04:41:19.127 回答