How can I plot the residuals of lm() with ggplot?

r ggplot2 lm
Lanza picture Lanza · Apr 20, 2016 · Viewed 22.6k times · Source

I would like to have a nice plot about residuals I got from an lm() model. Currently I use plot(model$residuals), but I want to have something nicer. If I try to plot it with ggplot, I get the error message:

ggplot2 doesn't know how to deal with data of class numeric

Answer

Gopala picture Gopala · Apr 20, 2016

Fortify is no longer recommended and might be deprecated according to Hadley.

You can use the broom package to do something similar (better):

library(broom)
y <-rnorm(10)
x <-1:10
mod <- lm(y ~ x)
df <- augment(mod)
ggplot(df, aes(x = .fitted, y = .resid)) + geom_point()

Related questions

Plot two graphs in same plot in R

I would like to plot y1 and y2 in the same plot. x <- seq(-2, 2, 0.05) y1 <- pnorm(x) y2 <- pnorm(x, 1, 1) plot(x, y1, type = "l", col = "red") plot(x, y2, type = "l", col = "…

Read More
Rotating and spacing axis labels in ggplot2

I have a plot where the x-axis is a factor whose labels are long. While probably not an ideal visualization, for now I'd like to simply rotate these labels to be vertical. I've figured this part out with the code …

Read More
How to set limits for axes in ggplot2 R plots?

I plot the following: library(ggplot2) carrots <- data.frame(length = rnorm(500000, 10000, 10000)) cukes <- data.frame(length = rnorm(50000, 10000, 20000)) carrots$veg <- 'carrot' cukes$veg <- 'cuke' vegLengths <- rbind(carrots, cukes) ggplot(vegLengths, aes(length, …

Read More