Add regression line (and goodness-of-fit stats) to scatterplot

Question

After reviewing other stackoverflow posts, I am attempting to add a regression line to my scatter plot with:

plot(subdata2$PeakToGone, subdata2$NO3_AVG, xlim = c(0, 70))
abline(lm(PeakToGone~NO3_AVG, data = subdata2))

However, it is not showing the line. I would also like to add the R^2, RMSE, and p-value from lm as text on the plot.

How can I add the regression line to the plot, along with these goodness-of-fit stats?

In general: It's best practice to provide a [reproducible example](http://adv-r.had.co.nz/Reproducibility.html), ready to copy, paste and run. By default, `plot` regards the 1st param as `x` and the 2nd as `y`. Try `plot(y=subdata2$PeakToGone, x=subdata2$NO3_AVG, ...` and look if the line is visible, now. You'll find an answer to your 2nd question here http://stackoverflow.com/questions/13114539/how-to-add-rmse-slope-intercept-r2-to-r-plot#answer-13119440 — lukeA, Oct 23 '14 at 23:17
@lukeA - that looks to be the drama, a simple mismatch in plotting or the specification of the `lm` formula. Add that as an answer as it is probably a common problem. — thelatemail, Oct 23 '14 at 23:19

score 1 · Accepted Answer · edited May 23 '17 at 11:57

1

By default, plot regards the 1st param as x and the 2nd as y. Try

plot(y = subdata2$PeakToGone, x = subdata2$NO3_AVG, xlim = c(0, 70))
abline(lm(PeakToGone~NO3_AVG, data = subdata2))

and look if the line is visible, now. You'll find an answer to your 2nd question here.

edited May 23 '17 at 11:57

Community

answered Oct 24 '14 at 10:24

lukeA

1 Answers1