Re: [R] Help on glmmPQL

From: <Lorenz.Gygax_at_fat.admin.ch>
Date: Thu 01 Jun 2006 - 15:50:26 EST


See help file of glmmPQL (the first thing you should be doing):

glmmPQL(fixed, random, family, data, correlation, weights,

        control, niter = 10, verbose = TRUE, ...)
                 ^^^^^^^^^^

Arguments
...
niter maximum number of iterations.
...

Thus, you can increase the value of niter to say 25. In my experience it is likely that the model fit is problematic if 10 iterations are not enough. Perhaps if you are lucky you were close and only some few more iterations are needed. If that is not the case, you may have too few occurrences of your response or too many explanatory variables.

Regards, Lorenz

> -----Original Message-----
> From: r-help-bounces@stat.math.ethz.ch
> [mailto:r-help-bounces@stat.math.ethz.ch]On Behalf Of
> Francisco Redelico
> Sent: Wednesday, May 31, 2006 4:54 PM
> To: r-help@stat.math.ethz.ch
> Subject: [R] Help on glmmPQL
>
>
> Dear R-Users,
>
> I am trying to fit a generalized linear mixed models
> using glmmPQL,
>
> mymodel.glmm = glmmPQL(efect~ time+ cov1+ cov2,
> random = ~ fact1+ fact2+ fact3+
> fact1*time + fact2*time+ fact3*- cov1 - cov2- time|
> subject, family=poisson)
>
> I get this error message:
>
> Error in lme.formula(fixed = zz ~ time + cov1 + cov2 ,
> random = ~fact1 + :
> iteration limit reached without convergence (9)
>
> I used
>
> meControl( maxIter = 10000, msMaxIter=10000,
> tolerance=1e-4, nlmStepMax=10000, msVerbose = TRUE)
>
> But it does not work.
>
> Could you tell me what is the error about? What does
> the number in brackets mean?
>
> Thanks in advance
>
> Francisco



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jun 01 17:06:31 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 01 Jun 2006 - 18:10:28 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.