plot4.R script