[R] Specifying starting values in lme (nlme package) using msScale

From: Catherine A. Holt <caholt_at_u.washington.edu>
Date: Tue, 18 Dec 2007 12:02:46 -0800 (PST)

I am using package nlme and would like to specify initial values for a linear mixed-effects model to help with convergence. I am trying to specify those initial values using the msScale option under ‘control’ in the lme() function:

lme(Y ~ X1, random= ~ X1|X2, control=list(msScale=lmeScale))

where, (as far as I understand), lmeScale is a function that can take initial values for parameters. However, I am unsure about how to input those starting values (e.g., what names lme will recognize for fixed and random effects, in what format, and if a partial list of initial values would be acceptable?).

Any advice or examples of code inputting starting values would be extremely helpful. I have been unable to find examples myself online. Note, although it may be easier to do this in the lme4 package, I would prefer to use nlme.

Thank you for your attention.
Carrie Holt, Ph.D., M.Sc., B.Sc.(Honours) University of Washington
School of Aquatic & Fishery Sciences
Box 355020
Seattle, WA 98195

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 Tue 18 Dec 2007 - 20:06:14 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 Tue 18 Dec 2007 - 20:30:19 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.