我正在尝试比较包含我的数据集中所有变量的各种混合效应逻辑（或某些结果伽玛）模型的AIC（或AICc）值，可以在此处下载其简化版本：https://drive.google.com/file/d/1YO17J7Dx1cFD0Wf3fNGe-a37ccTKyWNp/view?usp=sharing。但是我是glmulti的新手，仅使用lme4一个月左右，所以我认为我做错了什么：
我已经建立了包含所有变量的初始模型，如下所示：
data_reprex <- read_csv("reprex_glmulti_data.csv")
data_reprex <- within(data_reprex, {
Dog <- factor(Dog)
Box <-factor(Box)
Sex <- factor(Sex)
Group <- factor(Group)
Breedtype <- factor(Breedtype)
OwnerSeverity <- factor(OwnerSeverity, levels = c("None", "Mild", "Moderate", "Severe"))
})
outcome_model <- glmer(SuccessBinary ~ Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity +
Group + Group*Interval + Group*Box + Group*Trial + Group*Age + Group*Sex + Group*Breedtype
+ (1 | Dog), data = data_reprex, family = binomial, nAGQ=100)
该模型运行良好（除了我已经研究过的一些警告之外）。但是我想为模型运行glmulti并使用这些变量的每种组合来确定最佳模型，并根据我在网上找到的各种示例尝试了几种不同形式的语法，但是这些都不起作用（我也尝试过更改标准和“ confsetsize”设置）：
attempt1 <- glmulti(SuccessBinary ~ Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity +
Group + (1 | Dog), data=data_reprex,
level=1, fitfunction=glmer, family= binomial, crit="aicc", confsetsize=150)
attempt2 <- glmulti(SuccessBinary ~ Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity +
Group + (1 | Dog), data=data_reprex,
level=1, fitfunction=glmer, crit="aicc", confsetsize=150)
attempt3 <- glmulti(outcome_model, level=2, fitfunction=glmer, crit=AICc)
attempt4 <- glmulti(outcome_model, level=2, crit=AICc)
错误消息包括：
Improper call of glmulti.
Error in .subset2(x, i, exact = exact) :
attempt to select less than one element in get1index
有时我还会收到警告消息，例如：
In Ops.factor(Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity, :
‘+’ not meaningful for factors
In Ops.factor(Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity + : ‘|’ not meaningful for factors
我假设这是一个语法问题，我需要更改代码的结构，因为我使用的是Logistic glmer，而不是像我在网上找到的大多数示例一样使用glm，并且从文档中我认为glmulti应该兼容使用来自lme4的模型。请有人能建议我如何组织它以便运行吗？ （此外，是否有一种相对简单的方法来仅包括与变量“ Group”的交互，而不包括所有其他变量的交互？）
编辑：作为最后的尝试，我也尝试使用此处提供的建议（尽管我认为我可能有些困惑）：glmulti and liner mixed models使用以下代码，但这也不起作用：
glmer2.glmulti <- setMethod('getfit', 'merMod', function(object, ...) {
summ=summary(object)$coef
summ1=summ[,1:2]
if (length(dimnames(summ)[[1]])==1) {
summ1=matrix(summ1, nr=1, dimnames=list(c("(Intercept)"),c("Estimate","Std. Error")))
}
cbind(summ1, df=rep(10000,length(fixef(object))))
})
attempt5 <- glmulti(SuccessBinary ~ Interval + Box + Trial + Age + Sex + Breedtype + OwnerSeverity +
Group + (1 | Dog), data=data_reprex,
level=1, fitfunction=glmer2.glmulti, family=binomial, crit="aicc", confsetsize=150)
提前致谢！
最佳答案
我很困惑，也尝试将glmulti应用于装配了glmmTMB的混合效果模型，并且得到了与您非常相似的错误消息。我没有设法使它与glmulti一起使用，但是看来MuMIn包here中的函数dredge对我有用。也许您可以尝试一下？ （以防您尚未找到解决方案）。