Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Local and Global Stability of Certain Mixed Monotone Fractional Second Order Difference Equation with Quadratic Terms
Next Article in Special Issue
Spatial Statistical Models: An Overview under the Bayesian Approach
Previous Article in Journal
Foundations of Engineering Mathematics Applied for Fluid Flows
Previous Article in Special Issue
Is Football/Soccer Purely Stochastic, Made Out of Luck, or Maybe Predictable? How Does Bayesian Reasoning Assess Sports?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Discrete Weibull Marshall–Olkin Family of Distributions: Properties, Characterizations, and Applications

by
Jiju Gillariose
1,
Oluwafemi Samson Balogun
2,*,
Ehab M. Almetwally
3,
Rehan Ahmad Khan Sherwani
4,
Farrukh Jamal
5,* and
Joshin Joseph
6
1
Department of Statistics, CHRIST (Deemed to Be University), Bangalore 560029, India
2
School of Computing, University of Eastern Finland, 70211 Kuopio, Finland
3
Department of Statistics, Delta University of Science and Technology, Belkas 35631, Egypt
4
College of Statistical and Actuarial Sciences, University of the Punjab, Lahore 54000, Pakistan
5
Department of Statistics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
6
Department of Finance, KMML a Company under Government of Kerala, Sankaramangalam, Kollam 691583, India
*
Authors to whom correspondence should be addressed.
Axioms 2021, 10(4), 287; https://doi.org/10.3390/axioms10040287
Submission received: 16 September 2021 / Revised: 27 October 2021 / Accepted: 28 October 2021 / Published: 30 October 2021

Abstract

:
In this article, we introduce a new flexible discrete family of distributions, which accommodates wide collection of monotone failure rates. A sub-model of geometric distribution or a discrete generalization of the exponential model is proposed as a special case of the derived family. Besides, we point out a comprehensive record of some of its mathematical properties. Two distinct estimation methods for parameters estimation and two different methods for constructing confidence intervals are explored for the proposed distribution. In addition, three extensive Monte Carlo simulations studies are conducted to assess the advantages between estimation methods. Finally, the utility of the new model is embellished by dint of two real datasets.

1. Introduction

The Marshall–Olkin family of distributions was introduced by [1]. This is an interesting endeavour for associating additional parameter to an existing baseline. The resulting new probability distributions provide more flexibility to model various kinds of data. The early studies of the this scheme by [2,3,4,5,6,7,8], etc. offered its features. For a comprehensive discussion, one can refer a write-up by [9]. Keeping in view its applications, researchers focused on generalizations and extensions of the Marshall–Olkinn family of distribution. For example, exponentiated Marshall–Olkin extended family [10], Harris extended family [11], truncated negative-binomial generalized Marshall–Olkin family [12], beta Marshall–Olkin family [13], Kumaraswamy Marshall–Olkin family [14], truncated discrete Mittag–Leffler family [15], Marshall–Olkin generalized family [16], truncated discrete Linnik family [17], etc.
Recently, ref. [18] proposed a new family of distributions called the Weibull Marshall–Olkin generalized family. Its cumulative distribution function (cdf) with arbitrary baseline F ( x ; η ) , is given by
G ( x ; α , b , η ) = 1 exp log α F ¯ ( x , η ) 1 α ¯ F ¯ ( x , η ) b x > 0 , b > 0 , α > 0
and α ¯ = 1 α . This family is derived based on T-X family introduced by [19].
In many practical situations, discretizing statistical models from existing continuous distributions brings a crucial role. Consequently, some discrete distributions derived based on popular continuous models, for example, see, refs. [20,21,22,23,24,25], etc. Furthermore, several discrete distributions appeared in the literature by using continuous Marshall–Olkin distributions [26,27,28].
The most prominent motivating factor that pushed us towards the introduction of this new discrete distribution is the that while comparing with the volume of literature in continuous case, only limited research works have been inscribed about the discrete version of continuous family of distributions. The properties of the proposed model itself is another source of motivation. That is, the proposed new discrete model has positively skewed, decreasing and symmetric probability mass function (pmf). Moreover, the new distribution has decreasing, increasing, unimodal and bathtub shaped hazard rate functions (hrf). In addition to this, we offered a comparison of the estimation methods and real-life data application to explain how well the proposed model give consistently better fit than the other well-known discrete distributions.
The rest of the paper is outlined as follows: In Section 2, we propose a new discrete family of distributions and a special case of the derived family together with its probabilistic properties. Two distinct estimation procedures for parameters estimation are illustrated in Section 3. In Section 4, we give demonstration of the two different methods for constructing confidence intervals for the proposed distribution. Section 5 inculcates the Monte Carlo simulations study. Applications to the two real data sets to illustrate the performance of the new family are explained in Section 6. Section 7 offers some concluding remarks.

2. Discrete Weibull Marshall–Olkin Exponential Distribution

In this section, we propose a new discrete family of distributions namely discrete Weibull Marshall–Olkin family of distributions. Different methods and constructions of discrete analogs of continuous distributions are given in [29]. If the underlying continuous life time X has the survival function G ¯ ( x ) = P ( X > x ) , the pmf of the discrete random variable associated with that continuous distribution can be written as
P ( X = x ) = G ¯ ( x ) G ¯ ( x + 1 ) ; x = 0 , 1 , 2 , . .
The new family is generated by discretizing the continuous survival in Equation (1) using Equation (2), we obtain new family of discrete distribution with the pmf g ( x ) . The pmf can be written as
g ( x ) = exp log α F ¯ ( x ) 1 α ¯ F ¯ ( x ) b exp log α F ¯ ( x + 1 ) 1 α ¯ F ¯ ( x + 1 ) b ; x = 0 , 1 , 2 . . .
The survival function of the discrete random variable having the pmf (3) is given by
G ¯ ( x ) = exp log α F ¯ ( x + 1 ) 1 α ¯ F ¯ ( x + 1 ) b ; x = 0 , 1 , 2 . . .
We study a special case of this family, namely, discrete Weibull Marshall–Olkin exponential (DWMOE) distribution. The exponential distribution was chosen because it is the most tractable model. In practice for modeling and analyzing of random phenomenons, other distributions could be used.
Let the parent distribution be exponential with parameter λ and survival function F ¯ ( x ) = e λ x . We set ρ = e λ and 0 < ρ < 1 . Then, the pmf of the new model using (3) is given by
g ( x ) = exp log α ρ x 1 α ¯ ρ x b exp log α ρ ( x + 1 ) 1 α ¯ ρ ( x + 1 ) b ; x = 0 , 1 , 2 . . .
We called this new distribution the DWMOE distribution with parameters α , b and ρ . Note that, when b = 1 , the distribution with pmf (5) reduces to discrete Marshall–Olkin distributions, when b = 1 and α = 1 , (5) becomes geometric distribution. The corresponding cumulative distribution function (cdf) corresponding to (5) can be written as
G ¯ ( x ) = exp log α ρ ( x + 1 ) 1 α ¯ ρ ( x + 1 ) b ; x = 0 , 1 , 2 . . .
The corresponding hrf is given by
h ( x ) = exp log α ρ x 1 α ¯ ρ x b exp log α ρ ( x + 1 ) 1 α ¯ ρ ( x + 1 ) b exp log α ρ ( x + 1 ) 1 α ¯ ρ ( x + 1 ) b
Figure 1 shows the plots of the pmf of the DWMOE distribution for various values of α , b and ρ . The pmf can be decreasing and upside-down bathtub shaped. Furthermore, Figure 2 illustrates some of the possible shapes of the hrf of DWMOE distribution for selected values of the parameters. From figures, we can see that the hrf can be increasing, decreasing, bathtub, and upside-down bathtub shaped. Therefore, the DWMOE distribution is very efficient for modeling various data sets.

2.1. General Characteristics of the DWMOE Distribution

In this section, we study some general properties of the DWMOE distribution.

2.1.1. Quantiles, Probability Generating Function, Mean and Variance

Suppose X follows the DWMOE with cdf G ( x ) = 1 G ¯ ( x ) . The quantile function X = Q ( u ) , 0 < u < 1 of DWMOE is given as
X = 1 log ρ log exp log 1 1 u 1 b α + α ¯ exp log 1 1 u 1 b 1 ,
where u follows uniform distribution with support (0,1). Equation (7) can be used to simulate the DWMOE random variable. First, simulate a random variable u and compute the value of X in (7), which is not necessarily an integer. The DWMOE random variate X is the largest integer X , denoted as [ X ] .
In particular, the Median M is given as
M = 1 log ρ log exp log 2 1 b α + α ¯ exp log 2 1 b 1
The probability generating function (pgf) of the DWMOE( α , b, ρ ) distribution is given by
P ( s ) = 1 + ( s 1 ) x = 1 s x 1 exp log α ρ x 1 α ¯ ρ x b
The mean and variance do not have compact forms, but still the solutions of following equation gives mean and variance of random variable having DWMOE( α , b, ρ ) model.
E ( X ) = x = 1 exp log α ρ x 1 α ¯ ρ x b
and
V ( X ) = x = 1 ( 2 x 1 ) exp log α ρ x 1 α ¯ ρ x b x = 1 exp log α ρ x 1 α ¯ ρ x b 2
Using statistical software, the mean and variance of the DWMOE( α , b, ρ ) distribution for different values of α , b and ρ are calculated in Table 1. From this, we can say that the mean increases with ρ and α for different values of b. Moreover, depending on the values of α and ρ , the mean of the distribution can be smaller or greater than its variance. Therefore, the parameters of the DWMOE( α , b, ρ ) distribution can be used to model different data sets.

2.1.2. Shannon Entropy

In probabilistic context, Shannon entropy is a measure of variation of the uncertainty, with higher entropy corresponding to less information. For a discrete random variable X with pmf g(x), the Shannon entropy is defined as
S ( x ) = x g ( x ) log g ( x )
Combining (8) and (5) gives
S ( x ) = x exp log α F ¯ ( x ) 1 α ¯ F ¯ ( x ) b exp log α F ¯ ( x + 1 ) 1 α ¯ F ¯ ( x + 1 ) b log exp log α F ¯ ( x ) 1 α ¯ F ¯ ( x ) b exp log α F ¯ ( x + 1 ) 1 α ¯ F ¯ ( x + 1 ) b
Consider the another representation of pmf of DWMOE, that is after rearranging the pmf of DWMOE, we get
g ( x ) = exp log α ρ x 1 α ¯ ρ x b 1 exp log α ρ x + 1 1 α ¯ ρ x + 1 b log α ρ x 1 α ¯ ρ x b
Note that when b and 0 < α < 1 then S ( x ) 0 . This indicates that smaller values of b increase the uncertainty in the distribution.

3. Estimation Methods

To acquire the unknown model parameters, we will look at two estimation methods in this section. Maximum likelihood estimation (MLE) and the Bayesian estimation method are the two methods.

3.1. MLE Method

Let X = ( X 1 , , X n ) a random sample of size n from the DWMOE distribution and x = ( x 1 , , x n ) its observed values. The log-likelihood equation of the vector Ω = ( α , b , ρ ) are given by
l ( Ω ) = i = 1 n log exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b .
The non-linear likelihood equations are obtained by differentiating Equation (10) with respect to the parameters α , b , and ρ , respectively:
l ( Ω ) α = b α i = 1 n 1 ρ x i 1 α ¯ ρ x i log α ρ x i 1 α ¯ ρ x i b 1 exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b b α i = 1 n 1 ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b 1 exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
l ( Ω ) b = i = 1 n log log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b i = 1 n log log α ρ x i 1 α ¯ ρ x i log α ρ x i 1 α ¯ ρ x i b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
and
l ( Ω ) α = b ρ i = 1 n x i 1 α ¯ ρ x i log α ρ x i 1 α ¯ ρ x i b 1 exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b b ρ i = 1 n x i + 1 1 α ¯ ρ ( x i + 1 ) log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b 1 exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
Equating (11)–(13) with zeros and solving simultaneously by using iterative methods such as Newton–Raphson, we obtain the MLEs of α , b , and ρ . The asymptotic inference for the parameter vector ξ can be based on the normal approximation of the MLE of ξ . Indeed, under some regular conditions stated in [30], we have ξ ^ is approximately normally distributed with mean ξ and asymptotic variance-covariance matrix Σ ξ .

3.2. Bayesian Estimation Model

We explore the Square error type of loss function in Bayesian estimation. Several authors discussed this loss function based on various models [31,32]. Moreover, several authors discussed Bayesian estimation for parameters of the discrete distribution. For example, see [33,34,35], etc. For the Bayesian analysis, the prior distributions for each parameter of DWMOE distribution are essential. When prior information about the parameters is lacking, the Bayesian study might use the non-informative prior. Most previous studies of Weibull and Marshall–Olkin family distributions have Gamma prior distribution. Therefore, we supposed the informative priors for parameters α and b. The independent gamma distribution is our prior distribution of choice since it plays a significant role in estimating the parameters. While the parameter ρ , we used Jeffreys’ rule.
The independent joint prior density function of α , b can be written as follows:
π ( α , b ) α s 1 1 b s 2 1 e ( h 1 α + h 2 b )
The prior distribution of ρ may be considered to be
π ( ρ ) 1 ρ ( 1 ρ )
The joint prior density function of Ω can be written as follows:
π ( Ω ) α s 1 1 b s 2 1 ρ ( 1 ρ ) e ( h 1 α + h 2 b )
From the likelihood function and joint prior function (16), the joint posterior density function of Ω is obtained. The joint posterior of the DWMOE distribution can be written as
Π ( Ω | x ) α s 1 1 b s 2 1 ρ ( 1 ρ ) e ( h 1 α + h 2 b ) i = 1 n exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
Using the most common loss function (symmetric), which is a squared error. Bayesian estimators of Ω ˜ based on the squared error loss function are defined by Markov Chain Monte Carlo (MCMC).
S ( Ω ˜ ) = E Ω ˜ Ω 2 0 0 0 Ω ˜ Ω 2 Π ( Ω | x ) d α d p d ρ .
As the integrals given by (18) cannot be estimated directly. Because of this, we used the MCMC to find an approximate value of integrals. A significant sub-class of the MCMC techniques is the Gibbs sampling and more general Metropolis within Gibbs samplers. The Metropolis-Hastings (MH) bound with the Gibbs sampling were the two most common instances of the MCMC method. Because of the similarity of the MH algorithm with acceptance-rejection sampling, each iteration of the algorithm produces a corresponding candidate value from the proposal distribution. We used the MH within the Gibbs sampling steps to produce random samples of conditional posterior densities from the DWMOE distribution family:
Π ( α | p , ρ , x ) α s 1 1 e h 1 α i = 1 n exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
Π ( b | α , ρ , x ) b s 2 1 e h 2 b i = 1 n exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b ,
and
Π ( ρ | α , b , x ) 1 ρ ( 1 ρ ) i = 1 n exp log α ρ x i 1 α ¯ ρ x i b exp log α ρ ( x i + 1 ) 1 α ¯ ρ ( x i + 1 ) b .

4. Confidence Intervals

In this section, we propose two different methods to construct confidence intervals (CI) for the unknown parameters of the DWMOE distribution, which are asymptotic confidence interval (ACI) in MLE and credible confidence interval in MCMC of α , b , and ρ .

4.1. Asymptotic Confidence Intervals

The most common method to set confidence bounds for the parameters is to use the asymptotic normal distribution of the MLE. In relation to the asymptotic variance-covariance matrix of the MLE of the parameters, Fisher information matrix I ( Ω ) , where it is composed of the negative second derivatives of the natural logarithm of the likelihood function evaluated at Ω ^ = ( α ^ , b ^ , ρ ^ ) . Suppose the asymptotic variance-covariance matrix of the parameter vector Ω is
I ( Ω ^ ) = I α ^ α ^ I b ^ α ^ I b ^ b ^ I ρ ^ α ^ I ρ ^ b ^ I ρ ^ ρ ^
where V ( Ω ^ ) = I 1 ( Ω ^ ) .
A 100 ( 1 γ ) % confidence interval for parameter Ω can be constructed based on the asymptotic normality of the MLE.
α ^ ± Z γ 2 I α ^ α ^ , b ^ ± Z γ 2 I b ^ b ^ , and ρ ^ ± Z γ 2 I ρ ^ ρ ^
where Z γ 2 is the percentile of the standard normal distribution with right tail probability γ 2 .

4.2. Highest Posterior Density

The Highest Posterior Density (HPD) Intervals: [36] discussed this technique to generate the HPD intervals of unknown parameters of the benefit distribution. In this study, samples drawn with the proposed MH algorithm should be used to generate time-lapse estimates. From the percentile tail points, for instance, a 1 γ % HPD interval can be obtained with two points for 3 t h parameters of DWMOE distribution from the MCMC sampling outputs. It is sometimes useful to present the posterior median to informally check on possible asymmetry in the posterior density of a parameter.
According to [36], the BCIs of the parameters of DWMOE distribution α , b , ρ can be obtained through the following steps:
1.
Arrange α ˜ , b ˜ , ρ ˜ as α ˜ [ 1 ] α ˜ [ 2 ] α ˜ [ M ] , b ˜ [ 1 ] b ˜ [ 2 ] b ˜ [ M ] and ( ρ ˜ [ 1 ] ρ ˜ [ 2 ] ρ ˜ [ M ] ) , where M denotes the length of the generated of MH algorithm.
2.
The 100 ( 1 γ ) % symmetric credible intervals of α , b , ρ are obtained as: α ˜ [ M γ 2 ] , α ˜ M 1 γ 2 , b ˜ [ M γ 2 ] , b ˜ M 1 γ 2 and ρ ˜ [ M γ 2 ] , ρ ˜ M 1 γ 2 .

5. Simulation

In order to compare the classical estimation methods, the Monte-Carlo simulation procedure is carried out: MLE, and Bayesian estimation method under square error loss function based on MCMC, for estimation of DWMOE lifetime distribution parameters by R program. Monte-Carlo experiments are performed based on t data generated on 10,000 random DWMOE distribution samples using Equation (7), where x has DWMOE lifetime for various parameter actual values as:
  • Case 1: α = 2 ; b = 0.75 with different ρ = 0.25 , 0.5 , 0.85 .
  • Case 2: α = 2 ; b = 2 with different ρ = 0.25 , 0.5 , 0.85 .
  • Case 3: α = 0.65 ; b = 2 with different ρ = 0.25 , 0.5 , 0.85 .
and different sample sizes n : ( 25 , 80 , and 150 ) . We could describe the best methods of estimators as minimizing estimators’ Bias, mean squared error (MSE) and confidence interval when γ = 0.05 .
The simulation results of the methods presented in this paper for point estimation are summarized in the Table 2, Table 3 and Table 4. We consider the Bias, MSE, lower and upper of confidence interval values in order to perform the required comparison between various point estimation methods. The following remarks can be noted from these tables:
1.
For fixed actual parameters of DWMOE distribution, the Bias, and MSE decrease as n increases.
2.
When ρ increases, the Bias and MSE for all parameters decrease.
3.
When b increases, the Bias and MSE for α and ρ parameters decrease.
4.
Bayesian estimation is the best estimation method according Bias and MSE.
5.
The confidence interval of Bayesian estimation is shorter confidence interval than MLE ACI for parameters of DWMOE distribution.

6. Applications

In this section, we illustrate the usefulness of the newly proposed DWMOE distribution. We fit the DWMOE distribution to two data sets and compare the results with those of the fitted Kumaraswamy-geometric (KG), discrete Lindley (DL), and geometric models. For comparing the these distributions, we estimated the values of unknown parameters by the maximum likelihood method with standard errors (SE). Moreover, the fitted distributions are compared using the −log likelihood (−logL), AIC (Akaike Information Criterion), BIC (Bayesian Information Criterion), Kolmogorov–Smirnov (K-S) statistic, the corresponding p-values, the values of Anderson-Darling ( A * ) and Cramér von Mises ( W * ). The considered data sets are presented below.
Data set 1: The first dataset represents COVID-19 data belong to Australia of 32 days, which are recorded from 3 September to 4 October 2020, see, ref. [37]. The data set is:
6 15 59 11 5 9 8 11 7 9 6 7 6 0 8 8 5 7 5 2 35 2 8 1 2 3 7 4 2 2 3
Data set 2: The second data set refers the integer part of the lifetime of fifty devices in weeks is given by [38]. The data set is:
0 0 1 1 1 1 1 2 3 6 7 11 12 18 18 18 18 18 21 32 36 40 45 46 47 50 55
60 63 63 67 67 67 67 72 75 79 82 82 83 84 84 84 85 85 85 85 85 86 86
Table 5 and Table 6 represent the values of descriptive study for the fitted DWMOE, KG, DL and geometric models for Data set 1 and Data set 2, respectively. From the study, the smallest −logL, AIC, BIC, K-S statistic, A * , W * and the highest p-values are obtained for the DWMOE distribution. From these results, the new DWMOE distribution is the adequate model than others. Therefore, it can be used for fitting these data sets. Figure 3 and Figure 4 give the estimated pdfs for Data set 1 and Data set 2, respectively. These plots also indeed affirm that DWMOE distribution is clearly a competitive model for considered datasets. Thus, the new model may be a best alternative to the considered discrete models for modeling these real data sets.

7. Conclusions

In this paper, we proposed a new discrete Weibull Marshall–Olkin family of distributions with two shape parameters. A special case of this family, discrete Weibull Marshall–Olkin exponential (DWMOE) distribution is studied in detail. The model parameters have been estimated by the maximum likelihood estimation and Bayesian methods. Simulation study is carried out to assess the performance of the different estimation methods. Finally, two real data sets are analyzed to show the importance and flexibility of specified distribution. We hope that the proposed model will serve as a substitute to various discrete distributions exist in Statistical literature.

Author Contributions

Conceptualization, J.G.; Data curation, E.M.A.; Funding acquisition, O.S.B.; Investigation, J.J.; Methodology, J.G.; Project administration, O.S.B.; Resources, J.J.; Software, F.J.; Validation, J.J.; Writing—original draft, J.G.; Writing—review & editing, R.A.K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This manuscript is supported by Digiteknologian TKI-ymparisto project A74338(ERDF, Regional Council of Pohjois-Savo).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marshall, A.W.; Olkin, I. A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrica 1997, 84, 641–652. [Google Scholar] [CrossRef]
  2. Castellares, F.; Lemonte, A.J. On the Marshall-Olkin extended distributions. Commun. Stat.-Theory Methods 2016, 45, 4537–4555. [Google Scholar] [CrossRef]
  3. Cordeiro, G.M.; Lemonte, A.J.; Ortega, E.M.M. The Marshall-Olkin family of distributions: Mathematical properties and new models. J. Stat. Theory Pract. 2014, 8, 343–366. [Google Scholar] [CrossRef]
  4. Economou, P.; Caroni, C. Parametric proportional odds frailty models. Commun. Stat.-Simul. Comput. 2007, 36, 579–592. [Google Scholar] [CrossRef]
  5. Lam, K.F.; Leung, T.L. Marginal likelihood estimation for proportional odds models with right censored data. Lifetime Data Anal. 2001, 7, 39–54. [Google Scholar] [CrossRef]
  6. Nanda, A.K.; Das, S. Stochastic orders of the Marshall-Olkin extended distribution. Stat. Probab. Lett. 2012, 82, 295–302. [Google Scholar] [CrossRef]
  7. Rao, G.S.; Ghitany, M.E.; Kantam, R.R.L. Reliability test plans for Marshall-Olkin extended exponential distribution. Appl. Math. Sci. 2009, 55, 2745–2755. [Google Scholar]
  8. Rubio, F.J.; Steel, M.F.J. On the Marshall-Olkin transformation as a skewing mechanism. Comput. Stat. Data Anal. 2012, 56, 2251–2257. [Google Scholar] [CrossRef] [Green Version]
  9. Tahir, M.H.; Nadarajah, S. Parameter induction in continuous univariate distributions:Well established G families. Ann. Braz. Acad. Sci. 2015, 87, 539–568. [Google Scholar] [CrossRef]
  10. Jayakumar, K.; Thomas, M. On a generalization of Marshall-Olkin scheme and its application to Burr type XII distribution. Stat. Pap. 2008, 49, 421–439. [Google Scholar] [CrossRef]
  11. Aly, E.; Benkherouf, E. A new family of distributions based on probability generating functions. Sankhya B-Appl. Interdiscip. Stat. 2011, 73, 70–80. [Google Scholar] [CrossRef]
  12. Nadarajah, S.; Jayakumar, K.; Ristić, M.M. A new family of lifetime models. J. Stat. Comput. Simul. 2012, 83, 1–16. [Google Scholar] [CrossRef]
  13. Alizadeh, M.; Cordeiro, G.; Brito, E.; Demétrio, C. The beta Marshall-Olkin family of distributions. J. Stat. Distrib. Appl. 2015, 23, 546–557. [Google Scholar] [CrossRef] [Green Version]
  14. Alizadeh, M.; Tahir, H.M.; Cordeiro, G.M.; Zubair, M.; Hamedani, G.G. The Kumaraswamy Marshal-Olkin family of distributions. J. Egypt. Math. Soc. 2015, 23, 546–557. [Google Scholar] [CrossRef] [Green Version]
  15. Sankaran, K.K.; Jayakumar, K. A new extended uniform distribution. International J. Stat. Distrib. Appl. 2016, 2, 35–41. [Google Scholar] [CrossRef] [Green Version]
  16. Yousof, H.M.; Afify, A.; Nadarajah, S.; Hamedani, G.; Aryal, G. The Marshall-Olkin Generalized-G Family of Distributions with Applications. Statistica 2018, 78, 273–295. [Google Scholar]
  17. Jayakumar, K.; Sankaran, K.K. Discrete Linnik Weibull distribution. Commun. Stat.-Simul. Comput. 2019, 48, 3092–3117. [Google Scholar] [CrossRef]
  18. Korkmaz, M.C.; Cordeiro, G.M.; Yousof, H.M.; Pescim, R.R.; Afify, A.Z.; Nadarajah, S. The Weibull Marshall-Olkin family: Regression model and applications to censored data. Commun. Stat.-Theory Methods 2019, 48, 4171–4194. [Google Scholar] [CrossRef]
  19. Alzaatreh, A.; Lee, C.; Famoye, F. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef] [Green Version]
  20. Chakraborty, S.; Chakravarty, D. A new discrete probability distribution with integer support on (−,). Commun. Stat.-Theory Methods 2016, 45, 492–505. [Google Scholar] [CrossRef]
  21. Krishna, H.; Pundir, P.S. Discrete Burr and discrete Pareto distributions. Stat. Methodol. 2009, 6, 177–188. [Google Scholar] [CrossRef]
  22. Lisman, J.H.C.; van Zuylen, M.C.A. Note on the generation of the most probable frequency distribution. Stat. Neerl. 1972, 26, 19–23. [Google Scholar] [CrossRef]
  23. Nakagawa, T.; Osaki, S. The discrete Weibull distribution. IEEE Trans. Reliab. 1975, 24, 300–301. [Google Scholar] [CrossRef]
  24. Sato, H.; Ikota, M.; Aritoshi, S.; Masuda, H. A new defect distribution meteorology with a consistent discrete exponential formula and its applications. IEEE Trans. Semicond. Manuf. 1999, 12, 409–418. [Google Scholar] [CrossRef]
  25. Stein, W.E.; Dattero, R. A new discrete Weibull distribution. IEEE Trans. Reliab. 1984, 33, 196–197. [Google Scholar] [CrossRef]
  26. Jayakumar, K.; Sankaran, K.K. A Discrete Generalization of Marshall-Olkin Scheme and its Application to Geometric Distribution. J. Kerala Stat. Assoc. 2017, 28, 1–21. [Google Scholar]
  27. Jayakumar, K.; Sankaran, K.K. A generalization of discrete Weibull distribution. Commun. Stat.-Simul. Comput. 2017, 47, 6064–6078. [Google Scholar] [CrossRef]
  28. Sandhya, E.; Prasanth, C.B. Marshall-Olkin discrete uniform distribution. J. Probab. 2014, 2014, 979312. [Google Scholar] [CrossRef]
  29. Chakraborty, S. A new discrete distribution related to generalized gamma distribution and its properties. Commun. Stat.-Theory Methods 2015, 44, 1691–1705. [Google Scholar] [CrossRef]
  30. Cox, D.R.; Hinkley, D.V. Theoretical Statistics; Chapman & Hall: London, UK, 1974. [Google Scholar]
  31. Bantan, R.; Hassan, A.S.; Almetwally, E.; Elgarhy, M.; Jamal, F.; Chesneau, C.; Elsehetry, M. Bayesian Analysis in Partially Accelerated Life Tests for Weighted Lomax Distribution. CMC-Comput. Mater. Contin. 2021, 68, 2859–2875. [Google Scholar] [CrossRef]
  32. Sabry, M.A.; Almetwally, E.M.; Alamri, O.A.; Yusuf, M.; Almongy, H.M.; Eldeeb, A.S. Inference of fuzzy reliability model for inverse Rayleigh distribution. Aims Math. 2021, 6, 9770–9785. [Google Scholar] [CrossRef]
  33. Archer, E.; Park, I.M.; Pillow, J.W. Bayesian entropy estimation for countable discrete distributions. J. Mach. Learn. Res. 2014, 15, 2833–2868. [Google Scholar]
  34. Ashour, S.K.; Muiftah, M.S.A. Bayesian Estimation of the Parameters of Discrete Weibull Type (I) Distribution. J. Mod. Appl. Stat. Methods 2020, 18, 12. [Google Scholar] [CrossRef]
  35. Irony, T.Z. Bayesian estimation for discrete distributions. J. Appl. Stat. 1992, 19, 533–549. [Google Scholar] [CrossRef]
  36. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  37. Almetwally, E.M.; Ibrahim, G.M. Discrete Alpha Power Inverse Lomax Distribution with Application of COVID-19 Data. Int. J. Appl. Math. Stat. Sci. 2020, 9, 11–22. [Google Scholar]
  38. Aarset, M.V. How to identify a bathtub hazard rate. IEEE Trans. Reliab. 1987, 36, 106–108. [Google Scholar] [CrossRef]
Figure 1. pmfs of the DWMOE ( α , b, ρ ) distribution for some parameter values.
Figure 1. pmfs of the DWMOE ( α , b, ρ ) distribution for some parameter values.
Axioms 10 00287 g001
Figure 2. hrfs of the DWMOE( α , b, ρ ) distribution for some parameter values.
Figure 2. hrfs of the DWMOE( α , b, ρ ) distribution for some parameter values.
Axioms 10 00287 g002
Figure 3. Plots of estimated pdfs of distributions for Data set 1.
Figure 3. Plots of estimated pdfs of distributions for Data set 1.
Axioms 10 00287 g003
Figure 4. Plots of estimated pdfs of distributions for Data set 2.
Figure 4. Plots of estimated pdfs of distributions for Data set 2.
Axioms 10 00287 g004
Table 1. The mean (variance) of DWMOE( α , b, ρ ) for different values of parameters.
Table 1. The mean (variance) of DWMOE( α , b, ρ ) for different values of parameters.
ρ 0.20.50.7
α
1.223.8886 (6880.8534)29.2683 (11,009.7047)36.9219 (14,310.5902)
b = 0.2 0.923.8283 (6879.0820)29.1251 (11,001.4607)36.6520 (14,291.7716)
1.20.2296 (0.2294)1.0254 ( 1.4259)2.4205 ( 5.4258)
b = 1.2 0.90.1680 (0.1741)0.8224 (1.1803)2.0072 (4.5549)
1.20.1166 (0.1032)0.9268 (0.5908)2.2734 (2.0174)
b = 2 0.90.0566 (0.0534)0.6978 (0.4811)1.8268 (1.6304)
1.20.0011 (0.0011)0.9671 (0.1231)2.3879 (0.4845)
b = 5 0.98.5316 ×  10 7 (8.5316× 10 7 )0.7933 (0.1662)1.9016 (0.3986)
1.21.1009 ×  10 20 ( 1.1009 ×  10 20 )0.9934 (0.0067)2.5263 (0.2673)
b = 10 0.91.55 ×  10 85 ( 1.55 ×  10 85 )0.9471 (0.0500)1.9548 (0.0917)
Table 2. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 1 with different value of ρ .
Table 2. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 1 with different value of ρ .
MLEBayesian
ρ n BiasMSELowerUpperBiasMSELowerUpper
α 0.24950.48250.97823.52080.11080.17551.43692.9437
25b0.36040.19410.61371.60720.13060.02720.68551.0716
ρ 0.14880.02760.25380.54390.13390.02190.26830.5074
α 0.05220.18491.21522.88930.02390.02731.70412.3423
0.2580b0.27320.09780.72511.32130.11670.02150.69691.0366
ρ 0.14140.02410.26640.51650.12060.01630.30000.4651
α 0.19580.15361.53012.86150.02310.02011.78412.3271
150b0.24620.07110.79561.19670.13110.02030.74211.0302
ρ 0.13800.02060.31150.46450.12220.01600.30770.4321
α −0.04020.20961.04592.87380.12750.04081.90282.3497
25b0.31610.17070.53281.59940.11940.01830.74610.9896
ρ 0.06240.01120.39060.73410.07230.00600.42380.6467
α −0.01440.12351.13612.73510.02490.02411.72442.3009
0.580b0.20030.05270.73061.17000.09060.01440.70171.0009
ρ 0.05920.00910.48850.69570.06450.00580.51750.6249
α 0.02370.03371.68432.38940.02400.02241.74172.3112
150b0.17860.03750.78141.07580.08560.01350.68220.9720
ρ 0.08380.00790.52700.64060.02300.00480.48480.6360
α −0.07060.12261.25672.60210.07940.06461.68482.5829
25b0.14130.04060.60961.17300.03970.00720.65690.9542
ρ 0.01020.00110.79930.9211−0.01010.00100.81690.8953
α −0.09490.07491.40132.40900.05270.06271.56942.5802
0.8580b0.09630.01460.70280.98970.05470.00670.69010.9238
ρ 0.01420.00050.82980.89870.00420.00030.82250.8907
α −0.05380.02611.64752.24500.00800.02381.74252.3330
150b0.08710.01050.73110.94310.05420.00580.69580.9046
ρ 0.01270.00030.84110.88440.00670.00020.83010.8810
Table 3. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 2 with different value of ρ .
Table 3. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 2 with different value of ρ .
MLEBayesian
ρ n BiasMSELowerUpperBiasMSELowerUpper
α 0.31390.13661.93122.69660.13650.07511.68672.5906
25b0.95321.21371.86994.03640.12570.02481.93022.3134
ρ 0.11610.01450.30310.42920.10260.01410.29320.4037
α 0.29350.09542.10482.48210.04930.02501.72612.3206
0.2580b0.90930.91202.33683.48180.09480.02091.86942.2937
ρ 0.12040.01410.34320.39760.11230.01380.34180.4161
α 0.27420.07872.15722.39120.05050.01751.78382.2669
150b0.91000.87282.49523.32480.07620.00881.96582.1777
ρ 0.11220.01350.30510.39930.10270.01240.34830.4075
α 0.13880.04161.84602.43170.11100.04071.86592.4055
25b0.52840.40351.83713.21980.05900.00701.95592.1861
ρ 0.06000.00420.51160.60840.05910.00400.52320.6014
α 0.14040.02611.98332.29740.03620.02411.77582.3729
0.580b0.52100.31192.12662.91540.06700.00602.18682.2654
ρ 0.06270.00410.53580.58960.06240.00400.53270.5860
α 0.10700.01241.78522.25100.02460.01212.04712.1669
150b0.36930.14651.90602.25380.05890.00582.19722.2566
ρ 0.06260.00400.53670.59150.06300.00390.54400.5812
α 0.01170.00281.91012.11330.00950.00211.95482.1258
25b0.05840.02171.79302.32380.01590.00551.86392.1513
ρ 0.00480.000130.83410.87550.00430.000120.82790.8804
α 0.02420.00231.94232.10610.00380.00201.97352.2802
0.8580b0.05630.02081.85652.35620.02250.00191.82792.2248
ρ 0.00610.000080.84390.86830.00600.000070.84210.8707
α 0.00820.000151.99082.02560.00360.000101.99782.0210
150b0.03240.00201.97082.09400.02110.00181.98552.0919
ρ 0.00620.000060.84820.86420.00610.000050.84560.8675
Table 4. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 3 with different value of ρ .
Table 4. Bias, MSE, Lower, and Upper of CIs for parameters of WMOE distribution: Case 3 with different value of ρ .
MLEBayesian
ρ n BiasMSELowerUpperBiasMSELowerUpper
α 0.62150.42910.86541.67760.24480.11640.54501.3491
25b0.80580.73662.22573.38580.02580.03041.94062.1294
ρ 0.17940.07610.25830.40040.17640.05340.30450.5271
α 0.39940.18250.75211.34660.09770.02840.50571.0224
0.2580b0.48320.28352.04452.92190.02870.01381.81912.2592
ρ 0.12360.01570.33410.41320.10870.01370.35120.5246
α 0.24880.14970.69301.36470.11780.02310.54190.9682
150b0.36100.19642.30102.90820.04280.01181.85862.2454
ρ 0.10680.01180.31590.39760.08360.01040.36420.5011
α 0.52230.33380.68801.65660.38910.24320.54031.5798
25b0.77860.74652.04403.51310.06410.07021.95072.1617
ρ 0.01220.00300.40720.61720.04530.00300.42420.6474
α 0.39750.18860.70111.39390.11130.02910.49440.9766
0.580b0.66000.46712.30843.01170.05720.01471.89942.2347
ρ 0.03820.00290.46340.61290.03590.00280.52040.6423
α 0.34280.18490.86471.29030.12060.02610.57500.9880
150b0.65970.45532.38122.93820.04900.01211.94152.3053
ρ 0.02890.00140.48140.57640.02180.00130.53060.6450
α 0.13510.04770.44571.12450.12190.04110.53291.3685
25b0.27230.12361.83262.71210.07460.01511.94712.3112
ρ 0.00210.00040.81050.8938−0.00460.00040.81830.8853
α 0.13860.03430.53821.03900.04050.01220.49710.9237
0.8580b0.20660.07561.83662.57670.05040.01031.86742.2138
ρ −0.00170.00020.81910.87760.00160.00020.82150.8798
α 0.12740.02120.63850.91630.05910.01200.53350.8672
150b0.19260.07272.09802.43040.02080.00561.90202.1739
ρ 0.00650.00020.78320.82990.00410.00020.78580.8482
Table 5. Estimated values, −logL, AIC, BIC, K-S statistics, p-value, A * and W * for data set 1.
Table 5. Estimated values, −logL, AIC, BIC, K-S statistics, p-value, A * and W * for data set 1.
DistributionEstimates(SE)−logLAICBICK-Sp-Value A * W *
DWMOE α ^  = 1657.191491.11431188.2286192.62580.12090.7370.540660.0831
(4233.4787)
b ^  = 0.421(0.1050)
ρ ^  = 0.3446 (0.5445)
KG a ^  = 3.1744(0.0171)94.9914195.9828200.380.193790.18071.890.2395
b ^  = 0.2943(0.0526)
p ^  = 0.5717(0.0062)
DL a ^  = 3.1744(0.0171)97.12205196.2441197.70980.16200.37031.38170.1894
geometric a ^  = 0.1194 (0.0198)98.01666198.0333199.4990.254640.031531.38170.1894
Table 6. Estimated values, −logL, AIC, BIC, K-S statistics, p-value, A * and W * for data set 2.
Table 6. Estimated values, −logL, AIC, BIC, K-S statistics, p-value, A * and W * for data set 2.
DistributionEstimates(SE)−logLAICBICK-Sp-Value A * W *
DWMOE α ^  = 11.5621(12.5867)238.0759482.1519487.88790.178920.08142.3640.3069
b ^  = 0.6771(0.1510)
ρ ^  = 0.9488(0.0138)
KG a ^  =0.4987(0.3082)240.1928486.3855492.12160.18550.06412.8760.4761
b ^  = 0.1129(0.1299)
p ^  = 0.8367(0.1758)
DL a ^  = 0.0424(0.0042)249.9767501.9534503.86540.201470.03456.29770.6529
geometric p ^  = 0.02142(0.0030)241.6264485.2527487.16470.193100.0480.19313.0769
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gillariose, J.; Balogun, O.S.; Almetwally, E.M.; Sherwani, R.A.K.; Jamal, F.; Joseph, J. On the Discrete Weibull Marshall–Olkin Family of Distributions: Properties, Characterizations, and Applications. Axioms 2021, 10, 287. https://doi.org/10.3390/axioms10040287

AMA Style

Gillariose J, Balogun OS, Almetwally EM, Sherwani RAK, Jamal F, Joseph J. On the Discrete Weibull Marshall–Olkin Family of Distributions: Properties, Characterizations, and Applications. Axioms. 2021; 10(4):287. https://doi.org/10.3390/axioms10040287

Chicago/Turabian Style

Gillariose, Jiju, Oluwafemi Samson Balogun, Ehab M. Almetwally, Rehan Ahmad Khan Sherwani, Farrukh Jamal, and Joshin Joseph. 2021. "On the Discrete Weibull Marshall–Olkin Family of Distributions: Properties, Characterizations, and Applications" Axioms 10, no. 4: 287. https://doi.org/10.3390/axioms10040287

APA Style

Gillariose, J., Balogun, O. S., Almetwally, E. M., Sherwani, R. A. K., Jamal, F., & Joseph, J. (2021). On the Discrete Weibull Marshall–Olkin Family of Distributions: Properties, Characterizations, and Applications. Axioms, 10(4), 287. https://doi.org/10.3390/axioms10040287

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop