R scripts for plots 1-4 were added