plot1

R code for plot1