How to implement Zero-Inflated Multinomial?

Otherwise another possible solution is to construct a large blockwise covariance matrix.