Added R script for plot1