add four plots, and R code needed to generate them