Classification of Pulsars using Extreme Deconvolution

We carry out a classification of the observed pulsar dataset into distinct clusters, based on the $P-\dot{P}$ diagram, using Extreme Deconvolution based Gaussian Mixture Model. We then use the Bayesian Information Criterion to select the optimum number of clusters. We find in accord with previous works, that the pulsar dataset can be optimally classified into six clusters, with two for the millisecond pulsar population, and four for the ordinary pulsar population. Beyond that, however we do not glean any additional insight into the pulsar population based on this classification. Using numerical experiments, we confirm that Extreme Deconvolution-based classification is less sensitive to variations in the dataset compared to ordinary Gaussian Mixture Models. All our analysis codes used for this work have been made publicly available.