created break point analysis function
moved data import into a modular script, added reference in diversity
analysis script
created a diversity metrics function and added source command in the
analysis script
created df with best model values and output to global.env