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

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

您的函數參數。