Animals2 <- data(Animals2, package = "robustbase") mod.lm <- lm(log(Animals2$brain) ~ log(Animals2$body)) x <- rstudent(mod.lm) gg_qq(x) # See http://stackoverflow.com/questions/4357031/qqnorm-and-qqline-in-ggplot2/#27191036