added four plots and R code