错误：ord < 26L 不是 TRUE

Question

我正在运行一个混合效果模型：

a1fit1.1=glmer(Outcome_AprWatch_Exclude ~ Herd_size_2 +  w4a + w4b + as.factor(AgeSex_4_Coded_Clean) + (1 | EncounterID), data=attack, family=binomial(link="logit"), verbose=TRUE, nAGQ=30) 

结果是二元 (0/1)，预测变量是连续的，除了 AgeSex 是一个因素。我有 303 个观察值嵌套在 173 个遭遇（这是我的随机变量）中。我已经成功使用此代码一年半了，运行了各种类似的模型，并且成功地操纵了集成点的数量（nAGQ=1-30 左右）。

由于某种原因，我现在在运行代码时开始收到错误消息，但只有当我有 nAGQ>1 时（我需要它，因为我想/需要使用自适应高斯厄米近似而不是拉普拉斯近似)。

这是我现在使用 30 个 intpoint 运行模型时遇到的错误

Error: ord < 26L is not TRUE

这是我现在收到的警告，其中 intpoint 值较低（即 2-7）

Warning message:
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
Model failed to converge with max|grad| = 0.00114606 (tol = 0.001)

同样，当我使用 nAGQ=1 运行它时，我没有收到任何错误

正如我所说，这是我多次运行的同一模型，之前没有错误。此外，我在 STATA 中检查了这个模型（具有 30 个 intpoints 的同一模型），没有发现任何错误，并且收到了与我最初在 R 中运行模型时没有错误的相同输出（仅在上周左右！ )。

我不确定为什么 lme4 现在给我已经成功运行了很长时间的代码错误。 lme4 最近是否更新过，是否存在与更新相关的错误？任何人都可以提供的任何帮助都会很棒。

仅供参考 - 我尝试卸载并重新安装 R、RStudio 和 lme4，但错误消息没有任何变化。

Answer 1

lme4 的 1.0 版本于去年问世，并对代码进行了重大更改。由于您在此之前已经开始使用lme4，因此可以想象您的 R 可能直到今天才更新包。至少这就是我上个月发生的事情。

实际上开发了一个 lme4.0 包来比较新旧lme4 的结果。这是一个链接 (http://hlplab.wordpress.com/2014/03/17/old-and-new-lme4/)。请注意，作者最近在页面顶部发布了更新。

（对不起，我应该刚刚评论，但我的权限不允许）