[R] Error messages of nlme

Grathwohl,Dominik,LAUSANNE,NRC/NT dominik.grathwohl at rdls.nestle.com
Tue Sep 10 08:25:14 CEST 2002


Hello, everyone,

I would like to know the meaning of the error messages of nlme
in order to get converging fits. The application is a simple 
non-linear model to a dataset of increasing Apo B. 
I use:

platform i386-pc-mingw32
arch     i386           
os       mingw32        
system   i386, mingw32  
status                  
major    1              
minor    5.1            
year     2002           
month    06             
day      17             
language R

The non-linear function:

mu2 <- function(k, A, B0, t) {
	c0 <- exp(B0)
	i  <- exp(k)
	c1 <- exp(A)
	c0 * exp(-i*t) + c1/i * (1-exp(-i*t))
}

For example this error message:

> m2e.1 <- nlme(y ~ mu2(k, A, B0, t), data=thelp, 
+ fixed = list(k + A + B0 ~ 1),
+ random = list(id = B0 ~ 1, f = B0 ~ 1),
+ start=c(k=-1.7, A=0.86, B0=1.44))
Error: Singularity in backsolve at level 0, block 1

or:

> m2e.15 <- update(m2e.13, corr = corExp())
Error in as.vector(dist(1:length(x))) : couldn't find function "dist"

What have I to do to get proper fits? Nevertheless, I get some nice fits 
with other models. From my experience with this models, they are relatively 
robust against the choose of starting values. 
Can somebody explain me the meaning of the error messages?

Kind regards,

Dominik
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list