KIBNB function

K-inflated Beta Negative Binomial distributions for fitting a GAMLSS model