Create plot_diversity_results.Rscript