Initial R code plot 4