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.merMod包中的lme4，并将'devtol'参数提升为1e-8（对我的模型不起作用），然后提升为1e-7。有了这个值，它就起作用了

使用confint（）

1 个答案: