R code for plot1.png