Ah so I see now that PyMC automatically transforms the bounded variables to be unbounded.
So in all my formulas i would have to apply the backwards transform then the forwards to get the right value back.
So what I’m looking for is, how can I figure out which transform and name mangling have been applied to a bounded variable?
edit: perusing the source code I’ve come up with model.k01_interval__.distribution.transform_used to get the transform. Is splitting the name at the 3rd _ starting from the right sufficiently general to retrieve the original name?