Added R code to create plots