Hi, Junpeng – I was just rereading this, and I realized that PyMC3 is itself confusing your two kinds of transformation, or at least I think so.
So the bounding transformation fits your model, because it transforms to an unbounded underlying RV (IIUC).
But there are also transforms like ordered which I do not believe fits your model (or does it?)
Is there a good source of information on what the transformation kwarg should and should not be used for?