Add R script to generate boxplots