Playing around with D. Hogg et al's Mixture of Gaussians code