混合模型lme4中的警告信息

Question

以下是使用lme4包拟合'glmer'模型时出现的警告信息的含义是什么？

Warning messages:
1: glm.fit: fitted probabilities numerically 0 or 1 occurred 
2: In mer_finalize(ans) : false convergence (8)

我正在尝试拟合的模型如下所示：

glmer(dummy ~ constituency.coa + I(governat.part) + I(district2) + gdp.cap + lula.power + ifdm + bf.cap + year + (1 | munname), data=pool, family=binomial(link = "logit"), REML=T, verbose=T)

谢谢

- user1172558

---

这可能是一些很好的背景信息：http://www.ats.ucla.edu/stat/mult_pkg/faq/general/complete_separation_logit_models.htm - Aaron left Stack Overflow

Answer 1

对于警告2，您可以增加迭代次数，缺省值为300，当您添加更多迭代时，可以查看其是否收敛。尝试如下：

glmer(dummy ~ constituency.coa + I(governat.part) + I(district2) + gdp.cap + lula.power + ifdm + bf.cap + year + (1 | munname), data=pool, family=binomial(link = "logit"), REML=T, verbose=T, control = list(maxIter = 600))

这将把迭代次数从300增加到600，但如果这不起作用，您可以尝试更多的迭代次数。

- cconnell

---

尽管在这些情况下进行更多的迭代通常没有帮助（不过还是值得一试的）。 - Ben Bolker

Answer 2

警告1：对于一个或多个观测值，拟合值变为0或1，但在逻辑回归下这不应该发生。原因有很多；其中一个在?glm的帮助页面上讨论过，但那只是指向其他文件的指针。由于只是一个警告，可能不是问题，但它是一个警告，表示拟合存在一些问题。

警告2：我不知道确切的含义，但代码告诉你优化程序声明拟合过程已经收敛到估计值，但这种说法是错误的，拟合实际上没有真正收敛。

需要查看的一件事是是否存在可分离性问题，其中一个预测变量或线性组合可以完美地分割0和1事件。

建议您在R-SIG-Mixed邮件列表上跟进此事，那里有真正的专家可以提供进一步的帮助。您可能需要提供拟合过程的更多细节（打开详细模式）甚至数据，以便诊断问题。

- Gavin Simpson