optimal_phyloregion function

Determine optimal number of clusters