R: Fangfehler in `nls`

Question

Ich passe einige exponentielle Daten mit nls.

Der Code Ich verwende ist:

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

expFit als

definiert ist

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

Das funktioniert gut für die meisten meiner Daten, für die die Startparameter (100, -3 und 0) gut funktionieren. Manchmal habe ich jedoch Daten, die mit diesen Parametern nicht gut gehen, und ich bekomme Fehler von nls (z. B. "singular gradient" oder solche Dinge). Wie "erfasse" ich diese Fehler?

Ich habe versucht, so etwas wie

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

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

Aber das wird nicht funktionieren zu tun, weil nls scheint die Ausführung und den Code nach nls zu stoppen nicht ausführen ...

Irgendwelche Ideen?

Dank nico

Answer 1

ich diesen Trick in der Regel verwenden:

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 

}