Re: [R] how to set the starting value in lme

From: Daniel Kaschek <daniel.kaschek_at_physik.uni-freiburg.de>
Date: Wed, 16 Mar 2011 12:40:11 +0100

Hi Jia,

in order to test if the failing parameter estimation really has to do with wrong initial values, you could first simulate data with your mixed effects model. Just assume parameters, produce a data set with these parameters, add noise and look if you can estimate the parameters correctly with varying initial guesses.

On Tue, 2011-03-15 at 12:34 -0700, jia liu wrote:
> Hi,
> When I fit the mixed-effects model by lme, I got theconvergence error code = 1
> message = iteration limit reached without convergence (9). Even after I increase the iteration numbers, it still did not work. The only way I found to avoid the error message is to increase thenumber of iterations for the EM algorithm, which need to be at least 1100. But the weird thing is that the correlation between two random coefficients is always around 1, which is not true. So I am thinking the error should be related with the improper starting value. But I do not know how to set up my desired starting value in lme function. Are there anybody know it? Or met the same problem with me before and know how to avoid the error? Thank you for your help.
> Bests,
> Jia
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

-- 
Daniel Kaschek
Physikalisches Institut, Freiburg
Hermann-Herder-Str. 3
79104 Freiburg

Office:  Westbau, 02020
Phone:   +49-761-203-8531

______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Received on Wed 16 Mar 2011 - 11:44:34 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 16 Mar 2011 - 13:00:22 GMT.

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

list of date sections of archive