我想在散点图中添加一个二次拟合的拟合线,但是点的排序在某种程度上搞砸了。
attach(mtcars)
plot(hp, mpg)
fit <- lm(mpg ~ hp + I(hp^2))
summary(fit)
res <- data.frame(cbind(mpg, fitted(fit), hp))
with(res, plot(hp, mpg))
with(res, lines(hp, V2))
这样可以在整个地方绘制线条,而不是通过散点图绘制。我确信这很简单,但我有点难过。
答案 0 :(得分:8)
绘制线条时,所有点都按接收顺序连接。您希望在连接点之前对hp
值进行排序
res <- data.frame(cbind(mpg, fitted(fit), hp))
res <- res[order(hp), ]
with(res, plot(hp, mpg))
with(res, lines(hp, V2))
获取
此外,为了获得更平滑的线条,您可以考虑预测除了您观察到的hp
值之外的其他点。适合您的模型后,您可以
php <- seq(min(hp), max(hp), length.out=100)
p <- predict(fit, newdata=data.frame(hp=php))
plot(hp, mpg)
lines(php, p)