增加lmer新版本的迭代次数?

时间:2013-10-20 14:42:18

标签: r iteration lme4 lmer

我刚刚将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)

非常感谢!

1 个答案:

答案 0 :(得分:16)

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) ) )

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