Analytical likelihood for Gaussian hierarchical distribution
If performing a hierarchical analysis and using the BoundedGaussianDistribution
as the distribution being inferred, there is potential to perform the integral for each pulsar analytically. This would required estimating the Q_{22}
distributions as a Gaussian mixture model and then analytically integrating each combination of products of Gaussian distributions in the sum (taking into account the different weights).