Rixing problems with diversity_stats.R, now generating correct taxonomy matrix and figures