Univariate and multivariate versions of the negative binomial-inverse Gaussian distributions with applications