一昨日のSZ Pscの件、2次曲線のfittingにnlsは、大鉈を振るいすぎというか、nlsは非線形の場合に使うものだと、Katさんに指摘を受けました。
で、
r <- nls(y~a+b*x+c*x^2, star=c(a=1, b=1, c=1))
は、
x2 <- x*x
r <- lm(y ~ x + x2)
のほうが良いとのことでした。
再計算してみて、答えは一緒でした。
ちなみに今回は苦肉の策で、2次曲線で近似して極小時刻を求めましたが、本来は、光度変化は2次曲線にはならないので、他の方法で決めるのが望ましいです。
