We used reweighted iterative least squares and penalised quasi-likelihood method with outputs fed into a Markov chain Monte-Carlo model with 10 000 iterations which resulted in convergence for all analyses (Browne WJ. MCMC estimation in MLwiN v2.31. Centre for Multilevel Modelling, University of Bristol).