Single-suture craniosynostosis is a condition of the sutures of the infant's skull that causes major craniofacial deformities and is associated with an increased risk of cognitive deficits and learning/language disabilities. In this paper we adapt to classification of synostostic head shapes a Bayesian methodology that overcomes the limitations of our previously published shape representation and classification techniques. We evaluate our approach in a series of large-scale experiments and show performance superior to those of standard approaches such as Fourier descriptors, cranial spectrum, and Euclidian-distance-based analyses.