Abstract
In this article, we propose the discrete version of the binomial exponential II distribution for modelling count data. Some of its statistical properties including hazard rate function, mode, moments, skewness, kurtosis, and index of dispersion are derived. The shape of the failure rate function is increasing. Moreover, the proposed model is appropriate for modelling equi-, over- and under-dispersed data. The parameter estimation through the classical point of view has been done using the method of maximum likelihood, whereas, in the Bayesian framework, assuming independent beta priors of model parameters, the Metropolis-Hastings algorithm within Gibbs sampler is used to obtain sample-based Bayes estimates of the unknown parameters of the proposed model. A detailed simulation study is carried out to examine the outcomes of maximum likelihood and Bayesian estimators. Finally, two distinctive real data sets are analyzed using the proposed model. These applications showed the flexibility of the new distribution.