adding plot images and R code to reproduce images