Creating diversity plotter R script 'calculate_diversity.R'