15

我刚刚将 lme4 更新到 1.0-4 版本,当我运行 lmer() 之前正在收敛的混合效果模型时,现在会打印以下警告:

Warning message:
In (function (fn, par, lower = rep.int(-Inf, n), upper = rep.int(Inf,  :
  failure to converge in 10000 evaluations

所以,我想尝试增加迭代次数,看看我是否可以解决这个问题。(我必须说我不知道​​是什么导致了警告,因为消息的第一部分听起来有点不透明)。无论如何,我在文档中读到现在我应该使用lmerControl(),但我无法实现它。有人可以给我一个具体的例子来说明你如何做到具体吗?(帮助文件没有帮助)。这是我的模型:

m <- lmer(RT ~ Factor1*Factor2 + (0+Factor1+Factor2|Subject) + (1|Subject)  + (1|Item) + (0+Factor1+Factor2|Item), data= data)

非常感谢!

4

1 回答 1

23

lmerControl函数允许您选择优化器并将控制参数传递给它。控制迭代次数或评估次数的参数因函数而异(如 帮助页面中所述lmerControl)。默认优化器是“Nelder_Mead”,对于该优化器选择,可以通过在“optCtrl”参数列表中指定“maxfun”来更改最大评估次数:

m <- lmer(RT ~ Factor1*Factor2 + (0+Factor1+Factor2|Subject) + 
               (1|Subject)  + (1|Item) + (0+Factor1+Factor2|Item),
          data= data, control=lmerControl(optCtrl=list(maxfun=20000) ) )

这并不能保证会达到收敛。(我的经验是默认最大值通常就足够了。)很有可能您的数据不足以支持模型的复杂性,或者模型的构造不正确以用于研究设计。

并感谢@NBrouwer 将这个建议扩展到glmerwith的说明glmControl

于 2013-10-20T16:26:58.487 回答