使用confint（）估算GLMM的CI时出错

Question

我有一组GLMM配有二元响应变量和一组连续变量，我想获得每个模型的置信区间。我一直在使用confint()函数，95％和profile方法，如果它应用于没有交互的模型，它没有任何问题。

但是，当我将confint()应用于具有交互的模型（连续*连续）时，我一直收到此错误：

> m1CI <- confint(m1, level=0.95, method="profile")
Error in zeta(shiftpar, start = opt[seqpar1][-w]) : 
  profiling detected new, lower deviance

该模型运行没有任何问题（虽然我应用了一个优化器，因为有些模型存在收敛问题），这里是其中一个的最终形式：

m1 <- glmer(Use~RSr2*W+RSr3*W+RShw*W+RScon*W+RSmix*W+(1|Pack/Year),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = data0516RS, family=binomial(link="logit"))

有谁知道为什么会这样，我该如何解决？

我正在使用R版本3.4.3和lme4 1.1-17

Answer 1

按照以下说明解决了这个问题：https://stat.ethz.ch/pipermail/r-sig-mixed-models/2014q3/022394.html

我使用confint.merModfrom lme4包，并提升'devtol'参数，首先是1e-8，这对我的模型不起作用，然后到1e-7。有了这个值，它就有效了