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

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 = start，fn = function（x）dd（mkpar（npar1，：      来自bobyqa的收敛码1：bobyqa - 最大函数数＆gt;评估超出

     

2：在optwrap中（优化器，par = start，fn = function（x）dd（mkpar（npar1，：     来自bobyqa的收敛代码1：bobyqa - 超出最大功能评估数

当我尝试：

confint(Model1, parm="beta_")

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

这些拦截是否可以，或者之前的警告是否会让我谨慎接受那些作为CI？

对不起，如果这是一个愚蠢的问题，再次，不必深入了解统计知识或编码经验。

Answer 1

它可能不会产生巨大的影响，但你应该尝试添加

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

到您的函数参数。