Add beta negative binomial distribution for extreme over dispersion

[LHMarshall]

Need first four moments for mgcv which might not be easy to calculate for this distribution.

Might be available in gamlss but this package might not be that stable.