R code and plots for graphs 1,2, and 3