Hi! I’m working on a dataset where observations are best described as zero-inflated, Poisson distributed, with an upper bound (that’s a relatively low number). I was hoping that I could simply apply pm.Truncated on pm.ZeroInflatedPoisson, along the lines of:

However, it seems that truncation is not implemented for this distribution (I get an error NotImplementedError: Truncation not implemented for SymbolicRandomVariable MarginalMixtureRV{inline=True}). What would be the best way to go around this limitation? Any hints would be much appreciated!