Probability of a certain value in a discrete distribution

Questions
4

You will have to declare a model and put the rv in context like this:

>>> with pm.Model() as model:
...  x = pm.Bound(pm.Geometric, upper=100)('x', 0.02685)
...

Then you can get logp at a sample point using:

>>> lp = model.logp({'x': 34})
>>> np.exp(lp)
0.010936472378101412

This is one of the FAQ’s ==> here’s the thread

1 Like