add plot4.R