A very useful approach [288, 390], is to define a scaling factor as the ratio of the total Galactic volume weighted by pulsar density to the volume in which a pulsar is detectable:
In practice, is calculated for each pulsar separately using a Monte Carlo simulation to model the volume of the Galaxy probed by the major surveys . For a sample of observed pulsars above a minimum luminosity , the total number of pulsars in the Galaxy.
For small samples of observationally-selected objects, the detected sources are likely to be those with larger-than-average luminosities. The sum of the scale factors (5), therefore, will tend to underestimate the true size of the population. This “small-number bias” was first pointed out [177, 181] for the sample of double neutron star binaries where we know of only five systems relevant for calculations of the merging rate (see Section 3.4.1). Only when does the sum of the scale factors become a good indicator of the true population size.
Despite a limited sample size, it has been demonstrated  that rigorous confidence intervals of can be derived using Bayesian techniques. Monte Carlo simulations verify that the simulated number of detected objects closely follows a Poisson distribution and that , where is a constant. By varying the value of in the simulations, the mean of this Poisson distribution can be measured. The Bayesian analysis  finds, for a single object, the probability density function of the total population is
The “beaming fraction” in Equation (5) is the fraction of steradians swept out by a pulsar’s radio beam during one rotation. Thus is the probability that the beam cuts the line-of-sight of an arbitrarily positioned observer. A naïve estimate for of roughly 20% assumes a circular beam of width and a randomly distributed inclination angle between the spin and magnetic axes . Observational evidence summarised in Figure 16 suggests that shorter period pulsars have wider beams and therefore larger beaming fractions than their long-period counterparts [268, 246, 39, 353].
When most of these beaming models were originally proposed, the sample of millisecond pulsars was and hence their predictions about the beaming fractions of short-period pulsars relied largely on extrapolations from the normal pulsars. An analysis of a large sample of millisecond pulsar profiles  suggests that their beaming fraction lies between 50 and 100%. Independent constraints on for millisecond pulsars come from deep Chandra observations of the globular cluster 47 Tucanae  and radio pulsar surveys  which suggest that and likely close to unity . The large beaming fraction and narrow pulses often observed strongly suggests a fan beam model for millisecond pulsars .
This work is licensed under a Creative Commons License.