当将遗传算法与 lme4 一起使用时，glmulti 无限期运行

我在 R 中使用 glmulti 进行模型平均。我的模型中有大约 10 个变量，使得详尽的筛选不切实际 - 因此我需要使用遗传算法 (GA)（调用：method =“g”）。

我需要包含随机效应，因此我使用 glmulti 作为 lme4 的包装器。此处提供了执行此操作的方法http://www.inside-r.org/packages/cran/glmulti/docs/glmulti http://www.inside-r.org/packages/cran/glmulti/docs/glmultiglmulti 包中还包含一个 pdf 文件，其中包含更多详细信息。问题是，当告诉 glmulti 在这个设置中使用 GA 时，它会无限期地运行，即使在找到最佳模型之后也是如此。

这是从 glmulti 包中包含的 pdf 中获取的示例：

library(lme4)
library(glmulti)

# create a function for glmulti to act as a wrapper for lmer:
lmer.glmulti <- function (formula, data, random = "", ...) {
lmer(paste(deparse(formula), random), data = data, REML=F, ...)
}

# set some random variables:
y = runif(30,0,10) # mock dependent variable
a = runif(30) # dummy covariate
b = runif(30) # another dummy covariate
c = runif(30) # an another one
x = as.factor(round(runif(30),1))# dummy grouping factor

# run exhaustive screening with lmer:
bab <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, random = "+(1|x)")

这很好用。问题是当我告诉它使用遗传算法时：

babs <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, random = "+(1|x)", method = "g")

它只是无限期地继续运行，并且 AIC 不会改变：

...

After 19550 generations:
Best model: y~1
Crit= 161.038899734164
Mean crit= 164.13629335762
Change in best IC: 0 / Change in mean IC: 0

After 19560 generations:
Best model: y~1
Crit= 161.038899734164
Mean crit= 164.13629335762
Change in best IC: 0 / Change in mean IC: 0

After 19570 generations:
Best model: y~1
Crit= 161.038899734164
Mean crit= 164.13629335762

... etc.

我尝试过使用告诉 glmulti 何时停止的调用（deltaB = 0，deltaM = 0.01，conseq = 6），但似乎没有任何效果。我认为问题一定出在设置函数上（？）。这可能是非常明显的事情，但我是 R 新手，我无法解决它。

任何对此的帮助将不胜感激。

我从包维护者那里收到了解决方案。问题是探索的模型数量是由参数决定的confsetsize。默认值为 100。

根据?glmulti，这个论证是：

要查找的模型的数量，即返回的置信集的大小。

解决办法是设置confsetsize使其小于或等于模型总数。

从OP中没有停止的例子开始：

babs <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, 
                random = "+(1|x)", method = "g")

glmulti将使用以下方法确定候选模型的总数method = "d"

babs <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, 
                random = "+(1|x)", method = "d")



Initialization...
TASK: Diagnostic of candidate set.
Sample size: 30
0 factor(s).
3 covariate(s).
...
Your candidate set contains 64 models.

因此，设置confsetsize小于或等于 64 将产生所需的行为。

babs <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, 
                random = "+(1|x)", method = "g", confsetsize = 64)

然而，对于小型模型，使用穷举搜索可能就足够了（方法=“h”）：

babs <- glmulti(y~a*b*c, level = 2, fitfunc = lmer.glmulti, 
                random = "+(1|x)", method = "h")