회귀모델에 대한 여러가지 가정을 테스트해본다.
library(UsingR)
잔차에 대한 표준편차를 구하는 여러 방법들
fit <- lm(child ~ parent, galton)
# 회귀 모델에 의한 표준편차
summary(fit)$sigma
## [1] 2.238547
# sd(y) * sqrt(1 - cor(x, y)^2)
sd(galton$child) * sqrt(1 - cor(galton$child, galton$parent)^2)
## [1] 2.237339
# sqrt(잔차 제곱합 / (n - 1))
sqrt(sum(fit$residuals^2) /927)
## [1] 2.237339
# sd(잔차)
sd(fit$residuals)
## [1] 2.237339
# sqrt(잔차 분산)
sqrt(var(fit$residuals))
## [1] 2.237339