Eu quero adicionar uma linha ajustada de um ajuste quadrático a um gráfico de dispersão, mas a ordem dos pontos está de alguma maneira confusa.

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

Isso desenha linhas por todo o lugar, em oposição ao smooh que se encaixa no gráfico de dispersão. Tenho certeza de que isso é bem direto, mas estou um pouco perplexo.