thankyou so much, I think I got the idea, I have my last question. How can I assign a ‘tag’ to that dist_bound to identify it for plotting ? Because if I put
dist_bounded = pm.Bound(pm.Geometric, upper=100).dist(**'names_customers'**,p=0.02685)
the program give me an error. Thankyou so much another time.