我有
library(ISLR)
attach(Wage)
# Polynomial Regression and Step Functions
fit=lm(wage~poly(age,4),data=Wage)
coef(summary(fit))
fit2=lm(wage~poly(age,4,raw=T),data=Wage)
coef(summary(fit2))
plot(age, wage)
lines(20:350, predict(fit, newdata = data.frame(age=20:350)), lwd=3, col="darkred")
lines(20:350, predict(fit2, newdata = data.frame(age=20:350)), lwd=3, col="darkred")
Run Code Online (Sandbox Code Playgroud)
预测线似乎是相同的,但为什么系数如此不同?你如何解读为此他们raw=T和raw=F.
我看到产生的系数与poly(...,raw=T)那些系数相匹配~age+I(age^2)+I(age^3)+I(age^4).
如果我想使用系数来"手动"(不使用predict()函数)获得预测,那么我应该注意什么?我该如何解释正交多项式的系数poly().
r ×1