plot1 R code