2017-12-22 19 views
0

여러 가지 S 자형 커브를 피팅하는 데 문제가 있습니다. 더 나은 피팅을위한 약간의 도움은 정말 감사하겠습니다.nls 열악한 곡선 피팅

나는 세 가지 모델을 시도하고있다 : 무게

4 개 매개 변수

nls(y ~ a + (k-a) /(1 + (x/c)^e)^m, start=list(a=a_start, k=k_start, c=c_start, m=m_start, e=e_start), control = list(maxiter = 50000, minFactor=1/2000, warnOnly=T), algorithm = "port") 

5 매개 변수

nls(y ~ a + (k-a) /(1 + exp(-1/c*(x-e)))^m, start=list(a=a_start, k=k_start, c=c_start, m=m_start, e=e_start), trace=T, lower=c(a=-10000, k=0, c=0, m=0, e=0), upper=c(a=1000, k=100000, c=50, m=5, e=50), control = list(maxiter = 50000, minFactor=1/2000, warnOnly=T)) 

5 매개 변수를

nls(y ~ a + (k-a) /(1 + exp(-1/c*(x-e)))^m, weights = (1/residuals(fit)^2), start=list(a=a_start, k=k_start, c=c_start, m=m_start, e=e_start), trace=T, lower=c(a=-10000, k=0, c=0, m=0, e=0), upper=c(a=1000, k=100000, c=50, m=5, e=50), control = list(maxiter = 50000, minFactor=1/2000, warnOnly=T)) 

4PL 좋은 피팅 enter image description here

5pl 좋은 피팅 enter image description here

가난한 피팅 enter image description here

시그 모이 기능 enter image description here

+0

나는 이러한 데이터 포인트를 S 자형 곡선으로 설명 할 수있는 방법을 추측 할 수 없습니다. 시그 모이 드 위쪽의 데이터가 없습니다. – Roland

+0

죄송합니다. 불량 피팅을 표시하려면 간격을 확대하십시오. 이제 전체 데이터 포인트 예제의 하단에 추가했습니다. – Ribas

답변

0

나는 부분적으로 내 문제를 해결했다. 저는 성장 초기에 좋은 관심을 보이기 때문에 그 지역에서 더 많은 중요성을 부여하기 위해 맞춤식 가중치를 사용했습니다.