Hierarchical mixtures of Gaussians for combined dimensionality reduction and clustering