R: bắt lỗi trong `nls`

Question

Tôi đang lắp một số dữ liệu theo hàm mũ bằng cách sử dụng nls.

Mã Tôi đang sử dụng là:

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

expFit được định nghĩa là

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

này hoạt động tốt cho hầu hết các dữ liệu của tôi, mà các thông số bắt đầu được cung cấp (100, -3 và 0) hoạt động tốt. Đôi khi, mặc dù, tôi có dữ liệu không hoạt động tốt với các thông số đó và tôi gặp lỗi từ nls (ví dụ: "số ít" hoặc những thứ tương tự). Làm thế nào để "bắt" các lỗi này?

tôi đã cố gắng để làm một cái gì đó giống như

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

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

Nhưng điều này sẽ không làm việc vì nls dường như dừng lại việc thực hiện và các mã sau nls sẽ không thực hiện ...

Bất kỳ ý tưởng?

Cảm ơn nico

Answer 1

Tôi thường sử dụng thủ thuật này:

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 

}