Added the ability to fit using different likelihood functions within the
Maximum Likelihood estimator. These slow the fitting down but enable
slightly better precision than the pure Poisson likelihood model. Note
that they cannot be used with the gradient optimiser (BFGS).