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

:man_facepalming: If there’s a Wald then your life is much easier. Just jump to the step where I created the shifted wald and use pm.Wald directly.