added plot4.R