Combine two distributions to new one & scale parameter modifying probability depending on category

Did you perhaps intend:
p = pm.math.exp(pm.logcdf(combined, x))?

Logcdf gives you a value in the range (-inf, 0], whereas p should be in the range of [0, 1]

1 Like