I am working on phylogenetic beta diversity in microbial assemblages across biomes, and interested in calculating betaNTI through the null distribution of betaMNTD by randomizing OTUs across the phylogeny and recalculating bMNTD 999 times. Does anyone have a script for this?