added plots and R code to produce them