「統計学入門」（東京図書） - R入門

目次

最小二乗法による単純回帰モデルのパラメーターα、βの推定（pp.356-357）

> dtf <- read.csv("table10_1.csv", header = TRUE)
> n <- nrow(dtf)
> yy <- dtf$w
> xx <- dtf$v1
> sxx <- sum(xx)
> syy <- sum(yy)
> sxx2 <- sum(xx ^ 2)
> sxxyy <- sum(xx * yy)
> syy2 <- sum(yy ^ 2)
> xxm <- mean(xx)
> yym <- mean(yy)
> sx2 <- sum((xx - xxm) ^ 2)
> sxy <- sum((xx - xxm) * (yy - yym))
> sy2 <- sum((yy - yym) ^ 2)
> b <- sxy / sx2
> a <- yym - b * xxm
> cat(sprintf("βの最小二乗推定量 = %.2f\n", b))
βの最小二乗推定量 = 71.98
> cat(sprintf("αの最小二乗推定量 = %.2f\n", a))
αの最小二乗推定量 = -25.56