Plot 3 script and graphic. setwd commented in plot1.R