How to verify that uncertainty (estimated from pymc3) is accurate?

Agreed, but I can’t figure out where there calculation of r differs from the calculation U_i above except for a factor of L.