We employ a flexible multilevel model using a Dirichlet process mixture.
The standard parametric assumption in the random effects model is relaxed.
The expected number of latent components in data is estimated as part of the analysis.
Latent subpopulations and outliers are identified among different Canadian regions.