[問題] glmmPQL無法收斂問題

看板R_Language作者 (ppp)時間9年前 (2016/11/21 10:52), 9年前編輯推噓0(002)
留言2則, 1人參與, 最新討論串1/1
[問題類型]: 程式諮詢(我想用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

11/25 07:00, , 1F
會不會某些ID之下沒有重覆?
11/25 07:00, 1F

11/25 07:01, , 2F
我不肯定答案,只是幫瞎猜
11/25 07:01, 2F
文章代碼(AID): #1OCc5IPv (R_Language)
文章代碼(AID): #1OCc5IPv (R_Language)