[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
入門(寫過其他程式，只是對語法不熟悉)
[問題敘述]:
我現在使用MASS裡的glmmPQL來分析資料，模型長這樣
model.n7.8.gmm.4<- glmmPQL(Y~ factor(x1) + factor(x2) + factor(x3) +
x4 + x5 + x6 + x7 + x8 +factor(x9) + factor(x10),
random = ~ x5 + x6 +x7 + x8 | ID,
family=binomial(link="logit"),
correlation=corLin(form=~DayOrder), data = data.m7)
iteration 1
iteration 2
Error in lme.formula(fixed = zz ~ factor(sport.7.15) + factor(sick.7) + :
nlminb problem, convergence error code = 1
message = iteration limit reached without convergence (10)
但是出現無法收斂的訊息
就算我random改成 random = ~ x5 | ID也是出現同樣的訊息
如果我random改成只有常數項 random = ~ 1 | ID 就跑得出來
所以是模型太複雜了嗎?(資料約14000筆)
想請問有沒有解決的方式
謝謝
[關鍵字]:
glmmPQL