前提
以下のスクリプトは正常に動作します
b0 = 0.5143477
、b1 = 0.1021859
と、それっぽい解が得られています
R
1# sample.R 2x <- -50:50 3a0 <- 0.5 4a1 <- 0.1 5y <- a0 + a1*x + rnorm(101) 6 7f <- y ~ b0 + b1*x 8res <- nls(f, start=c(b0=0, b1=0), trace=T) 9yy <- predict(res)
$ Rscript sample.R 1012.126 : 0 0 88.9638 : 0.5143477 0.1021859
一方で、こちらはエラーが出ます
b0 = 0.5
、b1 = 0.1
と、ちゃんと解が得られているのに!
R
1# sample.R 2x <- -50:50 3a0 <- 0.5 4a1 <- 0.1 5y <- a0 + a1*x # 誤差の項を取り払った 6 7f <- y ~ b0 + b1*x 8res <- nls(f, start=c(b0=0, b1=0), trace=T) 9yy <- predict(res)
$ Rscript sample.R 883.75 : 0 0 2.656736e-28 : 0.5 0.1 1.087765e-30 : 0.5 0.1 〜〜省略〜〜 1.087765e-30 : 0.5 0.1 1.087765e-30 : 0.5 0.1 Error in nls(f, start = c(b0 = 0, b1 = 0), trace = T) : number of iterations exceeded maximum of 50 Execution halted
質問したいこと
この挙動の理由が知りたいです(そんなこと聞いて何になると言われればそれまでですが...)
補足情報(FW/ツールのバージョンなど)
公式サイトや質問サイトを閲覧しましたが、現状は誤差が0
にならないように工夫するしかなさそうです
R
のバージョン 3.4.4
あなたの回答
tips
プレビュー