我想将二次拟合的拟合线添加到 scatteprlot,但点的排序不知何故搞砸了。

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))

这会在整个地方绘制线条,而不是通过散点图进行平滑拟合。我相信这很简单,但我有点难住了。

最佳答案

当你绘制一条线时,所有的点都按照它们被接收的顺序连接起来。看起来您想在连接点之前对 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)

关于r - R线图中点的排序,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/25730676/

10-12 18:09