This paper combines two estimation procedures: Iterative Generalized Least Squares as used in the software MLwiN; Gibbs Sampling as employed in the software BUGS to produce a modelling strategy that respects the hierarchical nature of the Teaching Styles data and also allows for the endogeneity problems encountered when examining pupil progress.