Added plots and R code