r - 如何将因子的参考值传递给glmulti？

Question

当使用 glmulti 包来适应多个嵌套模型时，我想传递一个因子的参考值（x2：“A”，“B”）。安装普通 glm 时，我会使用relevel(x1,"B"). 这在 glmulti 中不起作用，其中候选变量作为文本传递。

我的尝试：

multimodel<- glmulti(y="response", xr= c("x1", "relevel(x2,'A')", "x3", "x4"), 
   data=mydata, level=2, method="g", crit="aicc")

错误信息：

Initialization...
Error in .subset2(x, i, exact = exact) : 
attempt to select less than one element

关于如何解决这个问题的任何建议？

score 0 · Accepted Answer

没关系。我刚刚意识到因子的重新排序可以在分析之前直接在数据框中mydata$x2<- relevel(mydata$x2,"B")完成，这完全消除了我使用glmulti时遇到的具体问题。我懒得在发布 Q 之前不发现这一点。

原来的调用也可以作为普通的glm进行，这样也消除了这个问题：

multimodel<- glmulti(response ~ x1 + relevel(x2,"A") + x3 + x4,
   data=mydata, level=2, crit="aicc")

1 回答 1