マルチプルシリーズのNLSを見つける

Question

モデルにフィットしたい/ 70シリーズを使って機能的なフォームを探したい。私はそれらのそれぞれのための30の値があります。プロットは傾向が線形ではないので、NLSを使用しようとしていることを示しています。

これが正しいアプローチであるかどうかはわかりません。私は以前これを実装していません。私はモデル化するために何か他のものを使用すべきでしょうか？次のエラーが表示されています。

Error in nls(y ~ a * exp(b * x), start = list(a = a_start, b = b_start)) :
singular gradient

xの値は1:30です。これは70回繰り返されます。 Yの値は0と1の間で変化します。私は71のベクトル（長さ30）を使うべきですか、それらを2つの列に集めるべきですか？現時点では、私は2つの列にすべてを集めました。

データは次のようになります。

---

a_start<-0.2 
b_start<-2*log(2)/a_start 

私の質問が明確でない場合は、さらに明確化のためにお問い合わせください。どんな助力も高く評価されます。

Answer 1

問題は、あなたの数式は、あなたのデータを正確に反映した形ではないということです。 y〜a * exp（b * x）は決してあなたのデータのようには見えないので、意味を成さない鋭いグラデーションが得られます。

これを試してみてください：

#some made up data, since none provided 
x <- runif(1000) * 30 
y <- -0.9 * exp(-0.5 * x) + 1.1 + rnorm(length(x), sd=0.1) 

fit <- nls(y ~ a * exp(b * x) + c, start = list(a = -1, b=-.1, c=1)) 

plot(x,y) 
lines(x[order(x)], predict(fit)[order(x)], col="red", lwd=4)