R: errores de captura en `nls`

Question

Estoy ajustando algunos datos exponenciales usando nls.

El código que estoy usando es:

fit <- nls(y ~ expFit(times, A, tau, C), start = c(A=100, tau=-3, C=0)) 

expFit se define como

expFit <- function(t, A, tau, C) 
    { 
    expFit <- A*(exp(-t/tau))+C 
    } 

Esto funciona bien para la mayoría de mis datos, para lo cual los parámetros iniciales proporcionaron (100, -3 y 0) funcionan bien. A veces, sin embargo, tengo datos que no van bien con esos parámetros y recibo errores de nls (por ejemplo, "gradiente singular" o cosas por el estilo). ¿Cómo "atrapo" estos errores?

Me trataron de hacer algo como

fit <- NULL 
fit <- nls(...) 

if (is.null(fit)) 
    { 
    // Try nls with other starting parameters 
    } 

Pero esto no va a funcionar porque nls parece detener la ejecución y el código después de nls no ejecutará ...

¿Alguna idea?

Gracias nico

Answer 1

normalmente utilizo este truco:

params<-... # setup default params. 

while(TRUE){ 

fit<-NULL 
try(fit<-nls(...)); # does not stop in the case of error 

if(!is.null(fit))break; # if nls works, then quit from the loop 

params<-... # change the params for nls 

}