We summarize our efforts at understanding the galactic neutron star (NS) population using likelihood and bayesian analyses. These include determination of the velocity distribution of young, high-magnetic field objects; the spatial, velocity, period and luminosity distributions of millisecond pulsars; and a full analysis of high-field objects including modeling of pulsars’ radio beams with the implicit assumption that such objects, for fixed period and period derivative, are standard candles.