load("bwt.Rdata")
load("sleep.Rdata")

L1 <- lm(log(sleep$BrainWt) ~ log(sleep$BodyWt))

sigmasq <- sum(L1$resid^2) / 56

stdresid <- L1$residuals /sqrt(sigmasq)

qqplot(c(1:58) / 59, pnorm(stdresid), xlab = "Theoretical probabilities",
       ylab = "Sample probabilities", pch=19)
abline(a = 0, b = 1, col = 2)

qqnorm(stdresid)
abline(a = 0, b = 1, col = 2)