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]
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]