同一个回归的不同结果

Question

为什么我从

得到不同的结果

summary(lm(mpg~horsepower + I(horsepower^2),data = Auto))$coef

和

summary(lm(mpg∼poly(horsepower,2) ,data=Auto))$coef

PS：我正在练习 ISLR

的实验

Answer 1

poly 默认使用正交多项式。如果您使用 poly(..., 2, raw = TRUE) 它将使用原始多项式，在这种情况下结果是相同的。

如果您使用默认的正交多项式，那么尽管它以不同的方式对模型进行参数化，但模型仍然会给出相同的预测。也就是说，fitted(lm(...)) 对于您的两个模型都是相同的。

library(ISLR)

fo1 <- mpg ~ horsepower + I(horsepower ^ 2)
fo2 <- mpg ~ poly(horsepower, 2)
fo3 <- mpg ~ poly(horsepower, 2, raw = TRUE)

fm1 <- lm(fo1, Auto)
fm2 <- lm(fo2, Auto)
fm3 <- lm(fo3, Auto)

all.equal(coef(summary(fm1)), coef(summary(fm3)), check.attributes = FALSE)
## [1] TRUE

all.equal(fitted(fm1), fitted(fm2))
## [1] TRUE

all.equal(fitted(fm1), fitted(fm3))
## [1] TRUE

这里的讨论可能会有帮助：What does the R function `poly` really do?

同一个回归的不同结果

Different result from the same regression

r

glm

lm