Optimal model-based clustering with multilevel data