作者ppp1987 (ppp)
看板R_Language
標題[問題] glmmPQL無法收斂問題
時間Mon Nov 21 10:52:00 2016
[問題類型]:
程式諮詢(我想用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
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.161.23
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1479696722.A.679.html
※ 編輯: ppp1987 (140.109.161.23), 11/21/2016 10:55:41
→ andrew43: 會不會某些ID之下沒有重覆? 11/25 07:00
→ andrew43: 我不肯定答案,只是幫瞎猜 11/25 07:01