Initial work: plot 1 - 3 R code