added and ran plot4.R