adding R code for plot