Bayesian and likelihood methods for fitting multilevel models with complex level-1 variation