2
我對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?
對不起,如果這是一個愚蠢的問題,再次,不需要太多深入的統計知識或編碼經驗。