[问题类型]: 程式谘询(我想用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://webptt.com/cn.aspx?n=bbs/R_Language/M.1479696722.A.679.html ※ 编辑: ppp1987 (140.109.161.23), 11/21/2016 10:55:41