[问题类型]:
程式咨询(我想用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