Defining a custom multivariate distribution

you mean NegativeBinomial only accept vector mu? I think we will need to modify its logp as well. What is the final shape of mu?