Create diversity_stats.Rscript