커스텀 pert 퍼츠 커브에 맞춰 줄 수 있습니다.

Question

사용자 정의 갬 파트 (Gompertz) 방정식에 데이터를 맞추려고합니다. NLS (식는, 데이터 = GPdata는, = 목록 (K = 1), 알고리즘 = "기본"을 시작으로

오류 :

은 아마 이것을 알고 단일 기울기를 다음과 같이 그러나, 나는 오류가 발생 내 시작 값이 수렴에 충분하지 않다는 것을 의미하지만 추측 이외의 다른 방법이 있습니까? 또한 nls2 함수 (무차별 대입)를 사용하여 가능한 K 값의 거대한 벡터를 확장했지만 동일한 오류가 발생했습니다. 여기

샘플 데이터를 포함, 내 코드입니다.

t = c(0, 2.4, 4, 5.9, 8.3, 10.2, 12, 14.5, 16.1, 17.8, 19, 21.1, 23, 24.7, 26.4, 28.3, 
    30.4, 32.3, 34.2, 36) 
Cells = c(14619994708, 18945074477, 19999450160, 23461507263, 22803949639, 19400861775, 
     19806214226, 17390826226, 17390826226, 18186961935, 14066488593, 
     13285017004, 11058253613, 8085853098, 9168523935, 8435042020, 
     8560592341, 7734294268, 7720887693, 7503367578) 

GPdata = data.frame(t, Cells) 

Ninf = (1 * 10^6) * (1/0.001) 
Eq <- Cells ~ exp(log(Ninf) * (1 - exp(-K * t))) 
Gompstart = nls(Eq, data = GPdata, start = list(K = 1), 
     algorithm = "default", trace = TRUE) 

도움이나 의견을 크게 주시면 감사하겠습니다.

감사합니다.

Answer 1

1) t 및 K 모든 것을 참고 :

LHS >= min(LSH) = min(Cells) = 7503367578 > 1e9 = Ninf = max(RHS) >= RHS 

따라서 LHS> RHS 모든 t 및 K 그래서 모델이 작동하지 않습니다에 대한.

2) 지수는 상승 및 하락 현상을 모델링하지 않습니다.