This is a clustering approach based on a probability model. The key idea is to fit alternative models to the data and select the best-fitting model using indices of model fit. The question of taxonicity can be decided by comparing single cluster and multiple cluster models.