使用混合模型计算置信区间的警告

Question

我对 R 还很陌生，但一直在使用 lmer 运行模型。现在我正在尝试计算截距和条件的置信区间。

例如：

Model1 <- lmer(Response ~ Treatment + 
               (1+Treatment|Hospital), data=Data) 
confint(Model1)

我收到以下输出和警告，例如：

              2.5 %     97.5 %
.sig01         5.6742125  6.7891223
.sig02        -0.9234187 -0.5258181
.sig03         2.7687916  8.7684621
.sigma         1.8734195  2.3471127
(Intercept)   26.7283448 28.9332127
TreatmentDrugx 12.1278192 20.4275981

警告信息： 1: 在 optwrap(optimizer, par = start, fn = function(x) dd(mkpar(npar1, : 来自 bobyqa 的收敛代码 1：bobyqa -- 函数的最大数量 > 超出评估

2: 在 optwrap(optimizer, par = start, fn = function(x) dd(mkpar(npar1, : 来自 bobyqa 的收敛代码 1：bobyqa -- 超出函数评估的最大数量

当我尝试时：

confint(Model1, parm="beta_")

我没有收到警告，但只获得了截距和治疗水平的置信区间。

这些拦截是否正常，或者之前的警告是否应该让我谨慎接受这些作为 CI？

对不起，如果这是一个愚蠢的问题，同样，不需要太多深入的统计知识或编码经验。

Answer 1

它可能不会有很大的不同，但你应该尝试添加

control=lmerControl(optCtrl=list(maxfun=10000))

到你的函数参数。