Script plot1.R