We describe an empirical Bayesian approach to determine the most likely size of an astronomical population of sources of which only a small subset are observed above some limiting flux density threshold. The method is most naturally applied to astronomical source populations at a common distance (e.g.,stellar populations in globular clusters), and can be applied even to populations where a survey detects no objects. The model allows for the inclusion of physical parameters of the stellar population and the detection process. As an example, we apply this method to the current sample of radio pulsars in Galactic globular clusters. Using the sample of flux density limits on pulsar surveys in 94 globular clusters published by Boyles et al., we examine a large number of population models with different dependencies. We find that models which include the globular cluster two-body encounter rate, $\Gamma$, are strongly favoured over models in which this is not a factor. The optimal model is one in which the mean number of pulsars is proportional to $\exp(1.5 \log \Gamma)$. This model agrees well with earlier work by Hui et al. and provides strong support to the idea that the two-body encounter rate directly impacts the number of neutron stars in a cluster. Our model predicts that the total number of potentially observable globular cluster pulsars in the Boyles et al. sample is 1070$^{+1280}_{-700}$, where the uncertainties signify the 95% confidence interval. Scaling this result to all Galactic globular clusters, and to account for radio pulsar beaming, we estimate the total population to be 2280$^{+2720}_{-1490}$.

Currently three isolated radio pulsars and one binary radio pulsar with no evidence of any previous recycling are known in 97 surveyed Galactic globular clusters. As pointed out by Lyne et al., the presence of these pulsars cannot be explained by core-collapse supernovae, as is commonly assumed for their counterparts in the Galactic disk. We apply a Bayesian analysis to the results from surveys for radio pulsars in globular clusters and find the number of potentially observable non-recycled radio pulsars present in all clusters to be < 3600. Accounting for beaming and retention considerations, the implied birth rate for any formation scenario for all 97 clusters is < 0.25 pulsars per century assuming a Maxwellian distribution of velocities with a dispersion of 10 km s^{-1}. The implied birth rates for higher velocity dispersions are substantially higher than inferred for such pulsars in the Galactic disk. This suggests that the velocity dispersion of young pulsars in globular clusters is significantly lower than those of disk pulsars. These numbers may be substantial overestimates due to the fact that the currently known sample of young pulsars is observed only in metal-rich clusters. We propose that young pulsars may only be formed in globular clusters with metallicities with log[Fe/H] > -0.6. In this case, the potentially observable population of such young pulsars is 447^{+1420}_{-399} (the error bars give the 95% confidence interval) and their birth rate is 0.012^{+0.037}_{-0.010} pulsars per century. The mostly likely creation scenario to explain these pulsars is the electron capture supernova of a OMgNe white dwarf.

We place limits on the population of non-recycled pulsars originating in globular clusters through Monte Carlo simulations and frequentist statistical techniques. We set upper limits on the birth rates of non-recycled cluster pulsars and predict how many may remain in the clusters, and how many may escape the cluster potentials and enter the field of the Galaxy.