R code from A Primer