lines（）无法正确显示二次拟合

Question

我只是试图显示我使用lm（）生成的拟合，但是lines函数给我一个奇怪的结果，其中有多条线从一个点出来。

这是我的代码：

library(ISLR)
data(Wage)
lm.mod<-lm(wage~poly(age, 4), data=Wage)
Wage$lm.fit<-predict(lm.mod, Wage)

plot(Wage$age, Wage$wage)
lines(Wage$age, Wage$lm.fit, col="blue")

我曾尝试用dev.off（）重置剧情，但我没有运气。 我正在使用rStudio。 FWIW，如果仅使回归线性化，则该线显示得非常好，但是一旦将其变为二次方或更高（使用I（age ^ 2）或poly（）），就会得到一个奇怪的图。 而且，points（）函数可以与poly（）一起使用。

谢谢您的帮助。

Answer 1

因为您忘了先按年龄排序点，所以这些行将随机出现。 线性回归也是如此。 他之所以适用于直线，是因为沿着直线上的任何点集行驶……停留在直线上！

plot(Wage$age, Wage$wage)
lines(sort(Wage$age), Wage$lm.fit[order(Wage$age)], col = 'blue')

考虑增加线宽以获得更好的视图：

lines(sort(Wage$age), Wage$lm.fit[order(Wage$age)], col = 'blue', lwd = 3)

Answer 2

只是为了在绘制模型预测时添加另一个更通用的技巧：

经常使用的策略是创建一个新数据集（例如newdat ），其中包含一系列可能值中的预测变量的值序列。 然后使用此数据显示您的预测值。 在此数据集中，您可以很好地分散预测变量值，但是并非总是如此。 使用新的数据集，可以确保您的行代表变量范围内的均匀分布的值：

例

newdat <- data.frame(age=seq(min(Wage$age), max(Wage$age),length=1000))
newdat$pred <- predict(lm.mod, newdata=newdat)
plot(Wage$age, Wage$wage, col=8, ylab="Wage", xlab="Age")
lines(newdat$age, newdat$pred, col="blue", lwd=2)