Initialization for bounded random variable

You can set theinitval of a distribution so it starts at a valid point. But even better is to pass an interval transform to the distribution so that the sampler never proposes invalid values: pymc.distributions.transforms.Interval — PyMC 4.3.0+0.gea721e4a.dirty documentation