Add R scripts to generate 4 plots