Add plot4.R and output