- 元となる y 軸と x 軸のデータを作成
> y <- c(5,10,12,13,4,-2) > y [1] 5 10 12 13 4 -2 > x <- c(1,2,3,4,5,6) > x [1] 1 2 3 4 5 6
- 重回帰用の変数を計算
x の二乗を計算します。
> x2 <- x ** 2 > x2 [1] 1 4 9 16 25 36
- 回帰分析の実行
y = a･x^2 + b･x + c の各係数を求めます。
> lm(y ~ x + x2) Call: lm(formula = y ~ x + x2) Coefficients: (Intercept) x x2 -4.300 10.889 -1.768
よってこの場合、y = -1.768x^2 + 10.889x - 4.3、が近似であるという結果になります。