Perform model fit evaluation in Bayesian way when sampling from the custom distribution is not known

If someone would like to help with adding an InverseGaussian to pymc-experimental, check out: Implement InverseGaussian distribution · Issue #253 · pymc-devs/pymc-experimental · GitHub