错误: ord < 26L不正确

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时(我需要使用/需要使用自适应高斯hermite近似，而不是laplace近似)。

这是我现在运行我的模型时所遇到的错误，我用了30个点

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个点)，没有发现任何错误，并收到了与我最初在R中运行模型时没有错误的输出(就在上周左右！)。

我不知道为什么lme4现在给我错误的代码，我已经成功运行了很长时间。lme4最近更新了吗?是否有与更新相关的bug？任何人所能给予的任何帮助都将是美妙的。

我试着卸载和重新安装R，RStudio和lme4，没有改变错误信息.

Answer 1

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

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

(对不起，我本应该评论一下的，但我的许可不允许)