<<目次へ

多項式モデルの回帰分析

> x <- c(12, 15, 18, 21)
> y <- c(56, 46, 59, 70)
> x2 <- x^2
> result <- lm(y ~ x2 + x) #2次式の場合
> result

Call:
lm(formula = y ~ x2 + x)

Coefficients:
(Intercept)           x2            x  
   179.7500       0.5833     -17.4167  

> x3 <- x^3
> result2 <- lm(y ~ x3 + x2 + x) #3次式の場合
> reuslt2

Call:
lm(formula = y ~ x3 + x2 + x)

Coefficients:
(Intercept)           x3           x2            x  
   826.0000      -0.1543       8.2222    -140.6111  

#データのプロット（曲線の重ね描き）
> plot(x,y,col="blue",pch=19)
> xx <- seq(12,21,0.01)
> points(xx, 179.7500 + 0.5833*xx^2 - 17.4167*xx, pch=".") #2次曲線のプロット
> points(xx, 826.0000 - 0.1543*xx^3  + 8.2222*xx^2 - 140.6111*xx, pch=".",col="pink")　#3次曲線のプロット