2017-12-22 34 views
0

私はいくつかの異なるS字曲線にフィットするのに問題があります。 より良いフィッティングのための少しの助けは本当に感謝されるでしょう。nlsひどいカーブフィッティング

私は、3つの異なるモデルをしようとしています:重みで

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字状の曲線で記述できることをどのように推測できるのかわかりません。あなたはS字状の上部からのデータを持っていません。 – Roland

+0

申し訳ありませんが、間違っているフィッティングを表示するために間隔を拡大します。今度は完全なデータポイントの例を底に追加しました。 – Ribas

答えて

0

私は部分的に私の問題を解決してきました。 私は成長の初期に良いフィッティングに関心があったので、私はその地域をより重要視するためにカスタムウェイトを使用しました。