R code from my Linear Models class