計算混合模型置信區間的警告

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)) 

給你的函數的參數。