Categorical model with continuous dependent variable

[junpenglao]

@pdeboer

So I want to try modeling the probability as a beta function with a mean that depends on b

Do you mean beta function or beta distribution? A log-beta function is available as pm.distributions.dist_math.betaln. If you are trying to get a beta distribution you should do p1 = pm.Beta('p', balance_drift, balance_sd)