相关疑难解决方法(0)

lm()中的poly():原始与正交之间的差异

我有

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=Traw=F.

我看到产生的系数与poly(...,raw=T)那些系数相匹配~age+I(age^2)+I(age^3)+I(age^4).

如果我想使用系数来"手动"(不使用predict()函数)获得预测,那么我应该注意什么?我该如何解释正交多项式的系数poly().

r

25
推荐指数
1
解决办法
2万
查看次数

标签 统计

r ×1