Since the theta in Figure 8 is generated from a Gamma distribution, it is already positive. In PyMC3 you just need to specify the right shape, and then tt.sort it along the right axis.
You can also try to use an order transformation (by @aseyboldt), it’s not in our code base yet so if it doesnt work ask @aseyboldt
