I included the four plots and the four R codes