1 Introduction
The Weibull (W) distribution is one of the most popular and widely used distributions for failure time in reliability theory (see, Weibull (1951)). The cumulative distribution function (CDF) of W distribution is given by
(1) 
where is scale parameter and is shape parameter. Clearly, the exponential (E) distribution and the Rayleigh (R) distribution are special cases for and respectively. Unfortunately, the shape of the hazard rate function (HRF) of W distribution can be only increasing, decreasing or constant. So, more modifications and generalizations of W distribution are presented in the statistical literature to describe various phenomena in different fields, because in many applications, empirical hazard rate curves often exhibit nonmonotonic shapes such as a bathtub, unimodal and others. For example:

Keller et al. (1985) proposed inverse Weibull (IW) distribution. The CDF of IW distribution is given by
(2) 
Lai et al. (2003) introduced modified Weibull (MW) distribution. The CDF of MW distribution is given by
(3) where is an accelerating parameter. The exponentiated MW distribution was proposed by Jalmar et al. (2008).

Bebbington et al. (2007) proposed flexible Weibull (FxW) distribution. The CDF of FxW distribution is given by
(4) where is scale parameter. The exponentiated FxW distribution was presented by ElGohary et al. (2015a), the inverse flexible Weibull (IFxW) distribution was proposed by ElGohary et al. (2015b) and the exponentiated of it was studied by ElMorshedy et al. (2017).

Cordeiro et al. (2013) introduced exponentialWeibull (EW) distribution. The CDF of EW distribution is given by
(5) where is shape parameter.

Nadarajah et al. (2013) proposed exponentiated Weibull (EW) distribution. The CDF of EW distribution is given by
(6) where is shape parameter.

ElBassiouny et al. (2017) introduced exponentiated generalized WeibullGompertz (EGWGz) distribution. The CDF of EGWGz distribution is given by
(7) where is shape parameter. The mixure of 2EGWGz distribution was studied by ElBassiouny et al. (2016).
Moreover, some discrete versions of E, R, W distributions and its generalizations have been presented in the literature because in several cases, lifetimes need to be recorded on a discrete scale rather than on a continuous analogue. So, discretizing continuous distributions has received much attention in the literature. For example:

Toshio and Shunji (1975) introduced discrete Weibull (DsW) distribution. The probability mass function (PMF) of DsW distribution is given by
(8) where . Clearly, the discrete Rayleigh (DsR) distribution is a special case for which was presented by Dilip (2004).

GómezDéniz (2010) proposed generalization of geometric (GGo) distribution. The PMF of GGo distribution is given by
(9) The GGo distribution reduces to geometric or discrete exponential (DsE) distribution when

Jazi et al. (2010) introduced DsIW distribution. The PMF of DsIW distribution is given by
(10) 
Vahid et al. (2013) proposed discrete generalized exponential type II (DsGET2) distribution. The PMF of DsGET2 distribution is given by
(11) 
Vahid and Hamid (2015a) introduced exponentiated discrete Weibull (EDsW) distribution. The PMF of EDsW distribution is given by
(12) 
Vahid et al. (2015b) proposed discrete beta exponential (DsBE) distribution. The PMF of DsBE distribution is given by
(13) where .
On the other hand, in many practical situations, it is important to consider different bivariate continuous and discrete distributions that could be used to model bivariate lifetime data in many fields. So, several bivariate continuous and discrete distributions are available in the statistical literature. For example, Lee (1997), Karlis and Ntzoufras (2000), Wu and Yuen (2003), Yuen et al. (2006), Sarhan and Balakrishnan (2007), Kundu and Dey (2009), Morata (2009), Kundu and Gupta (2009), Ong and Ng (2013), Balakrishnan and Shiji (2014), Lee and Cha (2015), Rasool and Akbar (2016), Hiba (2016), ElBassiouny et al. (2016), ElGohary et al. (2016), Vahid and Kundu (2017), Mohamed et al. (2017), Kundu and Vahid (2018), ElMorshedy and Khalil (2018) among others. An excellent encyclopedic survey of various continuous and discrete bivariate distributions can be found in Balakrishnan and Lai (2009) and Johnson et al. (1997) respectively.
In this regard, we focus the aim of this paper on presenting a flexible discrete bivariate distribution called BDsIW distribution, which can be usefully applied not only by statisticians, but also by data analysis in many different disciplines, such as sports, engineering, and medical applications. The proposed discrete model can be obtained from 3independent DsIW distributions by using the maximization method as suggested by Lee and Cha (2015). The main reasons for introducing BDsIW distribution are:

The proposed model is a very flexible bivariate discrete distribution, and its joint PMF can take different shapes depending on the parameter values.

The generation from the proposed model is straight forward. So, the simulation experiments can be performed quite conveniently.

The marginals of the proposed model are DsIW distributions. Hence, the marginals are able to analyze the hazard rates in the discrete case.

The DsE and DsR distributions are special cases from the proposed model.
2 The BDsIW Distribution and Its Statistical Properties
2.1 Definition and interpretations
Suppose DsIW(), DsIW() and DsIW() and they are independently distributed. If
, then we can say that the bivariate vector
has a BDsIW distribution with the parameter vector = . We will denote this discrete bivariate distribution by BDsIW. If BDsIW, then the joint CDF of for and for is given by(14) 
The marginal CDFs for BDsIW distribution can be represented as follows
(15) 
The corresponding joint PMF of for is given by
(16) 
where
The expressions , and for can be easily obtained by using the relation
(17) 
Figure 1 shows the plots of the joint PMF of BDsIW distribution for different parameter values.
Figure 1. The scatter plots of the joint PMF for different parameter values.
From Figure 1, it is clear that the joint PMF can take different shapes depending on the model parameter values. Assume BDsIW,then the joint reliability function of can be expressed as
(18) 
where
and
Moreover, the bivariate hazard rate function (BHRF) of can be represented as
(19) 
where = Figure 2 shows the plots of the BHRF of BDsIW distribution for different parameter values.
Figure 2. The scatter plots of the BHRF for different parameter values.
From Figure 2, it is clear that the joint BHRF can take different shapes depending on the parameter values. Assume , then the HRF of the conditional distribution given is given by
(20) 
and the HRF of the conditional distribution given is given by
(21) 
Similarly, when , then
(22) 
and
(23) 
On the other hand, assume a parallel system contains 2component. Then, we can defined the BHRF as a vector which is useful to measure the total life span of a 2component as follows
(24) 
where gives the HRF of the system using the information that the 2component has survived beyond , gives the HRF span of the first component given that it has survived to an age and the other has failed at . Similar argument holds for , (see Cox (1972)). If BDsIW, then
and
The following shock model and maintenance model interpretations can be provided for BDsIW distribution.

Shock model: Consider a system has 2component, and it is assumed that the amount of shocks is measured in a discrete unit. Each component is subjected to individual shocks say and respectively. The system faces an overall shock , which is transmitted to both the component equally, independent of their individual shocks. So, the observed shocks at the 2 component are and respectively.

Maintenance model: Consider a system has 2component, and it is assumed that each component has been maintained independently and also there is an overall maintenance. Due to component maintenance, assume the lifetime of the individual component is increased by amount and because of the overall maintenance, the lifetime of each component is increased by amount. Here, , and are all measured in a discrete unit. So, the increased lifetimes of the 2component are and respectively.
2.2 Some statistical properties
Assume BDsIW, then and are positive quadrant dependent (PQD) where
(25) 
Furthermore, for every pair of increasing functions and , we get ; see for example Nelsen (2006). Let us recall that, the function , is said to have the total positivity of order two () property if satisfies
(26) 
for all . It is consider a very strong and an important property in lifetime testing, see for example Hu et al. (2003). Assume and from BDsIW, then the joint reliability function of satisfies the property where
(27) 
Similarly, when etc. Now, we present some statistical properties of the proposed model in a form of results.
Result 1. If the bivariate vector has the BDsIW, then

DsIW

The stressstrenght probability is given by
(28) 
The median of and is given by
(29) where has a uniform distribution.

The coefficient of median correlation between and is given by
(30) 
The conditional PMF of given , is given by
(31) where
and

The conditional CDF of given , is given by
(32)
Proof. The proofs are quite standard and the details are avoided.
Result 2. Assume BDsIW for and they are independently distributed. If
Proof. It is easy to proof that using the joint CDF.
Result 3. If the bivariate vector BDsIW, then the joint probability generating function (PGF) of and can be written as infinite mixtures,
(33) 
Proof. The proof can be easily obtained by using the fact that
Hence, different moments and product moments of BDsIW distribution can be obtained, as infinite series, using the joint PGF.
3 Statistical Inference
3.1 Maximum likelihood estimation (MLE)
In this section, we use the maximum likelihood method to estimate the unknown parameters and of BDsIW distribution. Suppose that, we have a sample of size , of the form from BDsIW distribution. We use the following notations: and Based on the observations, the likelihood function is given by
(34) 
The loglikelihood function becomes
(35) 
where The MLEs of the model parameters can be obtained by computing the first partial derivatives of Equation (35) with respect to and and then putting the results equal zeros. We get the likelihood equations as in the following form
(36) 
Comments
There are no comments yet.