随机截距和随机斜率模型中的误差

Question

我在 R 中使用 lme4 包并尝试适应随机斜率和随机截距模型。 如果有人能在我运行随机斜率和随机截距模型时帮助我理解这个错误以及如何处理这个错误，我会非常有帮助：

mdl17<-lmer(yld.res ~ brk + (1+brk|state),data=data1,REML="FALSE")
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
   Model failed to converge with max|grad| = 1.84098 (tol = 0.002, component 3)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
   Model failed to converge: degenerate  Hessian with 1 negative eigenvalues

非常感谢

Answer 1

正如@CarlWitthoft 所说，警告会通知您模型收敛失败。它可以定义不足或过度定义。您还应该检查您尝试拟合的模型是否有意义。此外，您绝对应该通过运行update.packages("lme4") 或（如果您想更新所有软件包）update.packages() 来安装最新版本的lme4 软件包。当前的lme4 包包含它自己的关于收敛警告的页面。通过library(lme4) 加载它后，您可以通过发出?convergence 访问该手册页。在那里，您将找到有关您的警告消息的其他信息。 （请注意，未来版本可能不一定包含此页面。）作为附加措施，您可以搜索如何使用不同的优化器。 lme4 的新版本，例如导入可以使用的nloptr。

您还应该查看 lme4 的 Github 页面：https://github.com/lme4/lme4/，以防您遇到令您担心的警告或错误消息。对于您的具体情况，有一些可能相关的信息：

“lme4 (e.g. 1.1-6) 的最新版本给出了错误的收敛警告。在 r-sig-mixed-models 上有一个摘要帖子。如果您收到有关 max|grad| 的警告但模型通过了此测试：

dd <- fit@optinfo$derivs
with(dd,max(abs(solve(Hessian,gradient)))<2e-3)

那么您看到的是误报警告，并且该问题将在以后的版本中消失(1.1-7 and up)。"^[1]

如果Github页面包含相关信息（例如警告或错误特定于CRAN上的当前版本），您应该考虑从lme4安装lme4的当前master并检查是否有问题使用时消失：

(1) 首先安装devtools:install.packages("devtools") 并将其附加到您的命名空间library(devtools)。

(2) 然后使用：install_github("lme4/lme4", dependencies = TRUE) 从 Github 安装最新的 lme4 包 master。 （如果您遇到一个错误，表明尽管dependencies = TRUE 构建小插图失败，您应该将build_vignettes = FALSE 传递给您对install_github() 的调用。）

---

这里是收敛警告摘要的链接：http://thread.gmane.org/gmane.comp.lang.r.lme4.devel/11893（检索时间：2014-07-16T10:04+02:00）