Mirrored and bounded HalfNormal Distribution

See https://docs.pymc.io/api/bounds.html ; but you have to use keyword arguments

with pm.Model() as model:
    BHN = pm.Bound(pm.Normal, lower=0., upper=10.)
    u = BHN('u', mu=0., sd=2.)
    x = pm.Deterministic('x', 10. - u)
    trace = pm.sample_prior_predictive(25000)
1 Like