Rで回帰分析のため、回帰直線ではなく、近似曲線を求めてみます。 今回使うのは以下のようなデータ。 No,データ 1,121696 2,122012 3,122469 4,123033 5,123756 6,124454 7,125019 (略) 46,125983 47,126249 48,126643 49,127263 50,127816このデータをモデル化してplotしておきます。 x <- data[,1] y <- data[,2] model = (y~x) 明らかに5次式で表されるグラフのようですが、簡易的に三次式で近似します。 近似する場合には、まず未定係数の三次式を作成します。 y = a*x^3 + b*x^2 + c*x + dこの数式の係数を求めるために、nls関数を用います。 nlsを用いるときに、a,b,c,dのおおよその初期値を与えておくようです。(この初期値