Zero-Inflated Poisson Example

Yes, custom log likelihood could be implemented using pm.DensityDist