Percent point function (inverse of `CDF`) in pymc

Ahh that’s even better, thanks!

I just added the shift so it lines up with the (0,1) scale for a standard Uniform variable:
image

That other version you mention is the same calculation but the transformation is inside the erfinv function, whereas I just transformed first.