How to implement Zero-Inflated Multinomial?

Ok thanks Junpeng, I’ll try the loop then :wink:
I think a large blockwise covariance matrix is what I’m looking for, but I don’t know how to do that :sweat_smile: