Better commented plot4.R