Beta distribution with a different scale

Do you mean:

with pm.Model():
    p = pm.Beta('p', 1., 1.,)
    p_ = a + p * (b-a) # p_ is now in the range of (a, b)

? (All I did was remove the underscore from the equation involving p_.)