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

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（优化器中，par =开始，fn =函数（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))

您的函数参数。