[R] lme, best model without convergence

From: Thomas Wutzler <thomas.wutzler_at_web.de>
Date: Fri 26 May 2006 - 22:26:49 EST


Dear R-help list readers,

I am fitting mixed models with the lme function of the nlme package. If I get convergence depends on how the method (ML/REM) and which (and how much) parameters will depend randomly on the cluster-variable.

How get the bist fit without convergence?

I set the parameters msVerbose and returnObject to TRUE:

lmeControl(maxIter=50000, msMaxIter=200, tolerance=1e-4, niter=50, msTol=1e-5, nlmStepMax=500,

,msVerbose=TRUE
,returnObject=TRUE

)

However, the lme-functions does not produce verbose output, nor does it return the best fit if lme is not converging. It returns only an error:
Error in lme.formula(y ~ lndbh + I(lndbh^2) + lnh + I(lnh^2), random = ~lndbh + :

         iteration limit reached without convergence (9)

Best regards
Thomas



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 Fri May 26 22:26:47 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 Tue 30 May 2006 - 04:10:19 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.