Plot polynomial regression curve in R

by

I like to use ggplot2 for this because it’s usually very intuitive to add layers of data.

library(ggplot2)
fit <- lm(mpg ~ hp + I(hp^2), data = mtcars)
prd <- data.frame(hp = seq(from = range(mtcars$hp)[1], to = range(mtcars$hp)[2], length.out = 100))
err <- predict(fit, newdata = prd, se.fit = TRUE)

prd$lci <- err$fit - 1.96 * err$se.fit
prd$fit <- err$fit
prd$uci <- err$fit + 1.96 * err$se.fit

ggplot(prd, aes(x = hp, y = fit)) +
  theme_bw() +
  geom_line() +
  geom_smooth(aes(ymin = lci, ymax = uci), stat = "identity") +
  geom_point(data = mtcars, aes(x = hp, y = mpg))

enter image description here

More Related Contents:

  1. Messy plot when plotting predictions of a polynomial regression using lm() in R
  2. Linear Regression and group by in R
  3. How to change legend title in ggplot
  4. Position geom_text on dodged barplot
  5. Control ggplot2 legend look without affecting the plot
  6. Plot with conditional colors based on values in R [duplicate]
  7. Plot correlation matrix into a graph
  8. Increase distance between text and title on the y-axis
  9. Change the spacing of tick marks on the axis of a plot?
  10. Fast pairwise simple linear regression between variables in a data frame
  11. Align ggplot2 plots vertically
  12. How to overlay density plots in R?
  13. R: 4D plot, x, y, z, colours
  14. R Scatter Plot: symbol color represents number of overlapping points
  15. R: lm() result differs when using `weights` argument and when using manually reweighted data
  16. Is there a faster lm function
  17. How to make geom_text plot within the canvas’s bounds
  18. Change colors in ggpairs now that params is deprecated
  19. Plot two graphs in a same plot
  20. Plot 3D data in R
  21. Ordering of points in R lines plot
  22. Overlay two ggplot2 stat_density2d plots with alpha channels
  23. Showing equation of nls model with ggpmisc
  24. Plot inside a plot
  25. How do I draw gridlines using abline() that are behind the data?
  26. Adding ellipses to a principal component analysis (PCA) plot
  27. Finding where two linear fits intersect in R
  28. how to plot ellipse given a general equation in R?
  29. how to debug “factor has new levels” error for linear model and prediction [duplicate]
  30. Chronological timeline with points in time and format date

Leave a Comment