Added plots 1 through 4 plus R code