BICompare function

Identify modalities in a phylogeny