使用 confint() 估計 GLMM 的 CI 時出錯

Question

我有一組裝有二元響應變量和一組連續變量的 GLMM，我想獲得每個模型的置信區間。 我一直在使用confint()函數，使用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

Answer 1

按照以下說明解決了問題：

錯誤消息表明，在分析期間，優化器發現了一個比最初返回的假設最小偏差解決方案明顯更好的擬合值（由“devtol”參數表征）。 如果您想忽略這一點，您可以提高 'devtol' 參數（當前設置為保守的 1e-9 ...）——但是，非單調配置文件也警告您配置文件可能有問題.

來自https://stat.ethz.ch/pipermail/r-sig-mixed-models/2014q3/022394.html

我使用了confint.merMod包中的lme4 ，並將“devtol”參數提升到 1e-8，這對我的模型不起作用，然后提升到 1e-7。 有了這個值，它起作用了