Many applicative studies deal with multinomial responses and hierarchical data. Performing clustering at the highest level of grouping, in multilevel regression, is also often interest. In this study we analyse Politecnico di Milano data aim profiling students, modelling their probabilities belonging to different categories considering nested structure within engineering degree programmes. part...