我有一组装有二元响应变量和一组连续变量的 GLMM，我想获得每个模型的置信区间。我一直在使用confint()功能，在 95% 并且与profile方法，如果将其应用于没有交互的模型，那么它的工作不会出现任何问题。
然而，当我申请时confint()对于具有交互作用的模型（连续*连续），我收到此错误：
m1CI <- confint(m1, level=0.95, method="profile")
zeta(shiftpar, start = opt[seqpar1][-w]) 中的错误：
分析检测到新的、较低的偏差
该模型运行没有任何问题（尽管我应用了优化器，因为某些模型存在收敛问题），这是其中之一的最终形式：
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
按照以下说明解决了问题：
该错误消息表明，在分析过程中，优化器发现
明显更好的拟合值（其特征在于
'devtol' 参数）比返回的假定最小偏差解要大
首先。您可以提高“devtol”参数（即
目前设置为保守的 1e-9 ...）如果你想忽略这一点 -
然而，非单调的配置文件也警告您有些事情
可能对个人资料感到奇怪。
From https://stat.ethz.ch/pipermail/r-sig-mixed-models/2014q3/022394.html
我用的是confint.merMod来自lme4包，并将“devtol”参数提升，首先提升到 1e-8，这对我的模型不起作用，然后提升到 1e-7。有了这个值，它就起作用了
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
