Hierarchical model-based clustering of large datasets through fractionation and refractionation