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

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（优化器，PAR =启动，FN =函数（x）的DD（mkpar（npar1，：从bobyqa 收敛代码1：bobyqa - 功能评价的最大数目超过

当我尝试：

confint(Model1, parm="beta_") 

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

这些拦截是否正常，或者以前的警告是否使我谨慎接受那些CI？

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

Answer 1

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

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

给你的函数的参数。