qqnorm_ggplot.R

# For plotting qq
qqplot.data <- function(vec){
  y <- quantile(vec[!is.na(vec)], c(0.25, 0.75))
  x <- qnorm(c(0.25, 0.75)) 
  slope <- diff(y)/diff(x)
  int <- y[1L] - slope * x[1L]
  d <- data.frame(resids = vec)
  ggplot(d, aes(sample = resids)) + 
    stat_qq() + 
    geom_abline(slope = slope, intercept = int)
}
rstd.1.1 <- rstandard(model.1.1)
p <- qqplot.data(df$rstd.1.1)
print(p)