Next Article in Journal
Material Symmetries in Homogenized Hexagonal-Shaped Composites as Cosserat Continua
Next Article in Special Issue
A Multichannel Data Fusion Method Based on Multiple Deep Belief Networks for Intelligent Fault Diagnosis of Main Reducer
Previous Article in Journal
A Low-Rank Tensor Factorization Using Implicit Similarity in Trust Relationships
Previous Article in Special Issue
Evaluation of the Mass Diffusion Coefficient and Mass Biot Number Using a Nondominated Sorting Genetic Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Flexible Three-Parameter Model: Properties, Clayton Copula, and Modeling Real Data

by
Abdulhakim A. Al-babtain
1,
I. Elbatal
2,3 and
Haitham M. Yousof
4,*
1
Department of Statistics and Operations Research, King Saud University, Riyadh 11362, Saudi Arabia
2
Department of Mathematics and Statistics, College of Science Imam Muhammad ibn Saud Islamic University, Riyadh 11432, Saudi Arabia
3
Department of Mathematical Statistics, Faculty of Graduate Studies for Statistical Research, Cairo University, Benha 13513, Egypt
4
Department of Statistics, Mathematics and Insurance, Benha University, Benha 13513, Egypt
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(3), 440; https://doi.org/10.3390/sym12030440
Submission received: 5 February 2020 / Revised: 21 February 2020 / Accepted: 27 February 2020 / Published: 9 March 2020

Abstract

:
In this article, we introduced a new extension of the binomial-exponential 2 distribution. We discussed some of its structural mathematical properties. A simple type Copula-based construction is also presented to construct the bivariate- and multivariate-type distributions. We estimated the model parameters via the maximum likelihood method. Finally, we illustrated the importance of the new model by the study of two real data applications to show the flexibility and potentiality of the new model in modeling skewed and symmetric data sets.

1. Introduction and Motivation

The monotonicity of the hazard (failure) rate function (HRF) of a life model plays an important role in modeling failure time data. Probability distributions with an increasing failure rate (IFR) have various applications in pricing and supply chain contracting studies. The IFR property is a well-known and useful concept in reliability theory, dynamic programming, and other areas of applied probability and statistics (See [1,2]). Recently, [3] introduced a new two-parameter lifetime model with IFR. The model of [3] is named the binomial-exponential 2 (BE-2) model, which is constructed as the distribution of the random sum (RSum) of independent exponential random variables (IID RVs) when the sample size ( n ) has a zero truncated binomial (ZTB) model. The BE-2 distribution can be used as an alternative to the standard Weibull (W), standard gamma (Ga), exponentiated exponential (EE), and weighted exponential (WhE) distributions. The cumulative distribution function (CDF) of BE-2 distribution is given by:
G ϕ _ ( x ) | ( ϕ _ = ( α , β ) ) = 1 e α x ( 1 + β α x 2 β ) | ( x > 0 ) ,
where α > 0 is the scale parameter, and β is the shape parameter, where 0 β 1 . The probability density function (PDF) corresponding to (1) can be expressed as:
g ϕ _ ( x ) = α [ 1 + ( α x 1 ) β 2 β ] e α x = α 2 β e α x [ 2 ( 1 β ) + α x β ] .
The PDF in (1) can be written as:
g ϕ _ ( x ) =   P α e α x + ( 1 P ) α 2 x e α x ,
where P = 2 ( 1 β ) 2 β . The BE-2 distribution is a mixture of the exponential (E) distribution (with scale parameter α ), and the Ga distribution (with shape parameter 2 and scale parameter α ), with mixing proportion P . We notice that when β = 0 we get the standard exponential distribution, and when β = 1 the BE-2 distribution reduces to the gamma distribution with shape parameter 2 and scale parameter α . The BE-2 distribution has a PDF whose shape is like those of Ga, W, WhE, and EE distributions.
Recently, [4] proposed a general family of distributions called the Marshall–Olkin (MO-G) family of distributions. The MO-G family of distributions is also known as the proportional odds (PO) family. The CDF of the MO-G family is defined by:
F γ , ϕ _ ( x ) = G ϕ _ ( x ) 1 γ ¯ G ¯ ϕ _ ( x ) | ( x > 0 , γ > 0 ) .
The survival function (SF) F ¯ γ ( x ) is given by:
F ¯ γ , ϕ _ ( x ) = γ G ¯ ϕ _ ( x ) 1 γ ¯ G ¯ ϕ _ ( x ) ,
where γ ¯ = ( 1 γ ) , for γ = 1 , we get the baseline model, i.e.; F ¯ ( x ) = G ¯ ( x ) , where the shape parameter γ is called the tilt parameter. The PDF corresponding to (4) can be expressed as:
f γ , ϕ _ ( x ) = γ g ϕ _ ( x ) [ 1 γ ¯ G ¯ ϕ _ ( x ) ] 2 ,
and the HRF is given by:
h γ , ϕ _ ( x ) = f γ , ϕ _ ( x ) F ¯ γ , ϕ _ ( x ) = γ g ϕ _ ( x ) G ¯ ϕ _ ( x ) [ 1 γ ¯ G ¯ ϕ _ ( x ) ] .
The new PDF of the proposed lifetime model distribution can be right-skewed, symmetric, and left-skewed with many different useful shapes (see Figure 1), and this means that the new model will be suitable for modeling different real data sets, and the HRF of the new model exhibits many important HRF shapes such as the “increasing-constant”, “decreasing”, “increasing”, “constant”, and “bathtub” shapes (see Figure 2).
Practically, the proposed lifetime model is much better than many competitive versions of the exponential model, such as the odd Lindley exponential, the Marshall–Olkin exponential, moment exponential, the logarithmic Burr–Hatke exponential, the generalized Marshall–Olkin exponential, beta exponential, the Marshall–Olkin–Kumaraswamy exponential, the Kumaraswamy exponential, and the Kumaraswamy–Marshall–Olkin exponential, so the new lifetime model may be a good alternative to these models in modeling relief times and survival times data sets.

2. Genesis of the New Model

In this section, we introduce the three parameters of the MOBE-2 distribution. Using (1) and (2) in Equations (4)−(6), we obtain the CDF, SF, and PDF of the MOBE-2 distribution, (for x > 0 ) with vector of parameters Ψ _ = ( α , β , γ ) . The CDF and SF can be written as:
F Ψ _ ( x ) = 1 ( 1 + β α x 2 β ) e α x 1 γ ¯ ( 1 + β α x 2 β ) e α x ,
F ¯ Ψ _ ( x ) = γ ( 1 + β α x 2 β ) e α x 1 γ ¯ ( 1 + β α x 2 β ) e α x ,
respectively. The corresponding PDF can be derived as:
f Ψ _ ( x ) = γ α [ 1 + ( α x 1 ) β 2 β ] e α x [ 1 γ ¯ ( 1 + β α x 2 β ) e α x ] 2 .
Henceforth, let X MOBE-2 ( Ψ _ ) , with PDF (10). For the MOBE-2 distribution, the HRF can be written as:
h Ψ _ ( x ) = γ g ( x ) G ¯ ( x ) [ 1 γ ¯ G ¯ ( x ) ] = γ α [ 1 + ( α x 1 ) β 2 β ] ( 1 + β α x 2 β ) [ 1 γ ¯ ( 1 + β α x 2 β ) e α x ] .
The MOBE-2 distribution is a very flexible model that approaches different distributions when its parameters are changed. For = 0 , the MOBE-2 distribution reduces to the Marshall–Olkin extended exponential (MOEE) distribution. For = 1 , we get MOEGa distribution with shape parameter 2 and scale parameter α . For = 1 , we get BE-2 distribution (see Bakouch et al. (2014)). For β = 0 and = 1 , we get the exponential (E) distribution. For β = γ = 1 , the MOBE-2 distribution reduces to the Ga model with shape parameter 2 and scale parameter α . A useful representation for the new PDF is given in Appendix A. Figure 1 below gives some plots of the new PDF based on some selected parameters values. Based on Figure 1, we note that the new MOBE-2 distribution PDF can be right-skewed and left-skewed with many different useful shapes.
Figure 2 below gives some plots of the new HRF based on some selected parameters values. From Figure 2 we note that the HRF of the new model exhibits many important HRF shapes, such as the increasing-constant ( α = 1 , β = 1 , γ = 1 ) , decreasing ( α = 1 , β = 0.5 , γ = 0.1 ) , increasing ( α = 0.005 , β = 1 , γ = 1 ) , constant ( α = 0.001 , β = 0.001 , γ = 1 ) , and bathtub ( α = 0.05 , β = 0.05 , γ = 0.65 ) shapes.
The solution of the following relationship is used to find the quantile function (QF) of the MOBE-2 distribution, as follows:
( 1 + β α x q 2 β ) e α x q 1 q 1 γ ¯ q = 0 .
Since the uniform RVs are easily generated numerically in most statistical packages, the above scheme in (12) is very useful to generate MOBE-2 RVs and therefore can be easily implemented. It facilitates ready quantile-based statistical modeling. In particular, the median of X is Q ( 1 2 ) , given by setting q = 1 2 in (12). Also using (12), we can determine the Bowley’s skewness and the Moors’ kurtosis. The Bowley’s skewness is based on quartiles. Figure 3 indicates that both measures depend very much on the shape parameters β .

3. Properties

3.1. Moments

Theorem 1.
If X MOBE-2 ( Ψ _ ) , the r ( t h ) moment of X is given by:
μ r ( x ) = κ = 0 τ = 0 κ C τ , κ Γ ( r + τ ) ,
where:
C τ , κ = w τ , κ [ α ( 1 + κ ) ] r + τ + 2 [ 2 α ( 1 β ) ( 1 + κ ) + β α ( 1 + τ + r ) ] .
Proof of Theorem 1.
Let X be a RV following the MOBE-2 distribution. The r ( t h ) ordinary moment can be obtained using the well-known formula:
μ r ( x ) = 0 x r f Ψ _ ( x ) d x .
Then:
μ r ( x ) = γ α κ = 0 τ = 0 κ w τ , κ 0 x r + τ [ 2 ( 1 β ) + β α x ] e α ( 1 + κ ) x d x ,
setting = α ( 1 + κ ) x , and after some algebra, the r ( t h ) ordinary moment is given by:
μ r ( x ) = κ = 0 τ = 0 κ C τ , κ Γ ( r + τ ) ,
which complete the proof, where Γ ( n ) = 0 x n 1 e x d x refers to the gamma function. □

3.2. Moment Generating Function (MGF)

Theorem 2.
If X has the MOBE-2 ( Ψ _ ) , then the MGF of X is given as follows:
M X ( t ) = κ = 0 τ = 0 κ C τ , κ * Γ ( τ + 1 ) ,
where:
C τ , κ * = w τ , κ [ α ( 1 + κ ) t ] τ + 2 [ 2 α ( 1 β ) [ α ( 1 + κ ) t ] + β α ( τ + 1 ) ] .
Proof of Theorem 2.
The MGF can be derived from:
M X ( t ) = E ( e t X ) = 0 e t x f ( x ) d x ,
then we have:
M X ( t ) = κ = 0 τ = 0 κ w τ , κ 0 x τ [ 2 ( 1 β ) + β α x ] e [ α ( 1 + κ t ] x d x ,
which can be written as:
M X ( t ) = κ = 0 τ = 0 κ C τ , κ * Γ ( τ + 1 ) ,
which completes the proof. □

3.3. Conditional Moments

For any lifetime model, the s ( t h ) lower ( δ s ( t ) ) and upper ( π s ( t ) ) IM of X is defined by:
δ s ( t ) = E ( X s | ( X < t ) ) = 0 t x s f Ψ _ ( x ) d x ,
and
π s ( t ) = E ( X s | ( X > t ) ) = t x s f Ψ _ ( x ) d x ,
respectively, for any real s > 0 . The s ( t h ) lower incomplete moment of MOBE-2 distribution is:
δ s ( t ) = 0 t x s f ( x ) d x = κ = 0 τ = 0 κ w τ , κ 0 t x s + τ [ 2 ( 1 β ) + β α x ] e α ( 1 + κ ) x d x = κ = 0 τ = 0 κ ( { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( s + τ + 1 , α ( 1 + κ ) t ] } + β α [ Γ ( s + τ + 2 , α ( 1 + κ ) t ] ) ,
where:
ϒ τ , κ = w τ , κ [ α ( 1 + κ ) ] s + τ + 2 ,
and:
Γ ( s , t ) = 0 t x s 1 e x d x ,
is the lower incomplete gamma function, where:
Γ ( s , q ) | ( s 0 , 1 , 2 , ) = 0 q t s 1 e x p ( t ) d t = q s s { 1 F 1 [ s ; s q + 1 ; q ] } = κ = 0 ( 1 ) κ κ ! ( s + κ ) q s + κ ,
the function 1 F 1 [ , , ] is a called the confluent hypergeometric function. The first incomplete moment of X , denoted by, δ 1 ( t ) , is computed using Equation (24) by setting s = 1 as:
δ 1 ( t ) = κ = 0 τ = 0 κ V τ , κ ( { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( τ + 2 , α ( 1 + κ ) t ] } + β α [ Γ ( τ + 3 , α ( 1 + κ ) t ] ) ,
where:
V τ , κ = γ α ( 1 + κ ) γ ¯ κ ( β α ) τ ( 2 β ) τ + 1 [ α ( 1 + κ ) ] τ + 3 .
Similarly, the s ( t h ) upper incomplete moment of MOBE-2 distribution is:
π s ( t ) = t x s f ( x ) d x = κ = 0 τ = 0 κ w τ , κ t x s + τ [ 2 ( 1 β ) + β α x ] e α ( 1 + κ ) x d x = κ = 0 τ = 0 κ ϒ τ , κ ( { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( s + τ + 1 , α ( 1 + κ ) t ] } + β α [ Γ ( s + τ + 2 , α ( 1 + κ ) t ] ) ,
where ϒ τ , κ = w τ , κ [ α ( 1 + κ ) ] s + τ + 2 , and Γ ( s , t ) = t x s 1 e x d x , is the upper incomplete gamma function. The MRL is given by m X ( t ) = E ( X | ( X > t ) ) = π 1 ( t ) F ¯ ( t ) t , where π 1 ( t ) is the first incomplete moment of X and by setting s = 1 in Equation (37), we get:
m X ( t ) = 1 F ¯ ( t ) κ = 0 τ = 0 κ K τ , κ ( { 2 α ( 1 β ) ( 1 + κ ) Γ ( τ + 21 , α ( 1 + κ ) t } + β α Γ ( τ + 3 , α ( 1 + κ ) t ) ,
where K τ , κ = w τ , κ [ α ( 1 + κ ) ] τ + 3 . The mean inactivity time (MIT) of X is defined (for t > 0 ) by:
M X ( t ) = E ( X | X < t ) = t δ 1 ( t ) F ( t ) .
Then, we have:
M X ( t ) = t 1 F ( t ) κ = 0 τ = 0 κ K τ , κ { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( τ + 2 , α ( 1 + κ ) t ] + β α [ Γ ( τ + 3 , α ( 1 + κ ) t ] } .

3.4. Residual Life and Reversed Failure Rate Function

The n t h order moment of the residual life is given by the general formula (see [5]):
μ n ( t ) = E ( ( X t ) n | ( X > t ) ) = 1 F ¯ ( t ) t ( x t ) n f Ψ _ ( x ) d x   | ( r 1 ) .
Applying the binomial expansion of ( x t ) n and substituting f Ψ _ ( x ) given by (10) into the above formula gives:
μ n ( t ) = 1 F ¯ ( t ) h = 0 n ( t ) n h ( n h ) t x n f Ψ _ ( x ) d x = 1 F ¯ ( t ) κ = 0 τ = 0 κ h = 0 n ζ τ , κ { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( 1 + n + τ , α ( 1 + κ ) t ] + β α [ Γ ( 2 + n + τ , α ( 1 + κ ) t ] } ,
where ζ τ , κ = ( t ) n h [ α ( 1 + κ ) ] 2 + n + τ w τ , κ ( n h ) . The MRL of the MOBE-2 distribution is obtained by setting n = 1 in μ n ( t ) . The variance of the residual life of the MOBE-2 distribution can be obtained easily by using μ 2 ( t ) and μ(t). The r- t h moment of the reversed residual life (MRRL) can be obtained by the well-known formula:
m n ( t ) = E ( ( t X ) n | ( X t ) ) = 1 F ( t ) 0 t ( t x ) n f Ψ _ ( x ) d x   | ( n 1 ) .
Applying the binomial expansion of ( t x ) n and substituting f Ψ _ ( x ) , given before, into the above formula gives:
m n ( t ) = 1 F ( t ) h = 0 n ( t ) n h ( n h ) 0 t x n f Ψ _ ( x ) d x = 1 F ( t ) κ = 0 τ = 0 κ h = 0 n ζ τ , κ { 2 α ( 1 β ) ( 1 + κ ) [ Γ ( 1 + n + τ , α ( 1 + κ ) t ] + β α [ Γ ( 2 + n + τ , α ( 1 + κ ) t ] } .
The mean waiting time (MWT) of the MOBE-2 distribution can be obtained by setting n = 1 in m n ( t ) . Using m ( t ) and m 2 ( t ) , one can obtain the variance and the coefficient of variation of the reversed residual life of the MOBE-2 distribution (for more details see [2]).

4. Simple Type Copula-Based Construction

4.1. The Bivariate MOBE-2 Using the Morgenstern Family

First, we start with CDF for the Morgenstern family of two RVs ( X 1 , X 2 ) , which has the following form:
F λ ( x 1 , x 2 ) | ( | λ | 1 ) = F 1 ( x 1 ) F 2 ( x 2 ) { 1 + λ [ 1 F 1 ( x 1 ) ] [ 1 F 2 ( x 2 ) ] } ,
setting:
F γ 1 , α 1 , β 1 ( x 1 ) = 1 ( 1 + α 1 β 1 x 1 β 1 + 2 ) e α 1 x 1 1 γ ¯ 1 ( 1 + α 1 β 1 x 1 β 1 + 2 ) e α 1 x 1 ,
and:
F γ 2 , α 2 , β 2 ( x 2 ) = 1 ( 1 + α 2 β 2 x 2 β 2 + 2 ) e α 2 x 2 1 γ ¯ 2 ( 1 + α 2 β 2 x 2 β 2 + 2 ) e α 2 x 2 ,
then we have a seven-dimension parameter model.

4.2. Via Clayton Copula

4.2.1. The Bivariate MOBE-2 Model

The bivariate extension via Clayton copula can be considered as a weighted version of the Clayton Copula, which is of the form
C ( u , v ) = [ u ( δ 1 + δ 2 ) + v ( δ 1 + δ 2 ) 1 ] 1 δ 1 + δ 2 .
This is indeed a valid copula. Next, let us assume that X MOBE-2 ( γ 1 , α 1 , β 1 ) and Y MOBE-2 ( γ 2 , α 2 , β 2 ) . Then, setting:
u = u ( x ) = 1 ( 1 + α 1 β 1 x β 1 + 2 ) e α 1 x 1 γ ¯ 1 ( 1 + α 1 β 1 x β 1 + 2 ) e α 1 x ,
and:
v = v ( y ) = 1 ( 1 + α 2 β 2 y β 2 + 2 ) e α 2 y 1 γ ¯ 2 ( 1 + α 2 β 2 y β 2 + 2 ) e α 2 y ,
the associated CDF bivariate MOBE-2 type distribution will be:
H ( x , y ) = = { [ 1 ( 1 + α 1 β 1 x β 1 + 2 ) e α 1 x 1 γ ¯ 1 ( 1 + α 1 β 1 x β 1 + 2 ) e α 1 x ] ( δ 1 + δ 2 ) + [ 1 ( 1 + α 2 β 2 y β 2 + 2 ) e α 2 y 1 γ ¯ 2 ( 1 + α 2 β 2 y β 2 + 2 ) e α 2 y ] ( δ 1 + δ 2 ) 1 } 1 δ 1 + δ 2 .
Note: depending on the specific baseline CDF, one may construct various bivariate MOBE-2 type models in which ( δ 1 + δ 2 ) 0 .

4.2.2. The Multivariate Extension

The d -dimensional version from the above will be:
H ( x 1 , x 2 , , x d ) = { i = 1 d [ 1 ( 1 + α i β i x i β 1 + 2 ) e α i x i 1 γ ¯ i ( 1 + α i 1 β i x i β i + 2 ) e α i x i ] ( δ 1 + δ 2 ) + 1 d } 1 / ( δ 1 + δ 2 ) .
Further future works could be allocated for studying the bivariate and the multivariate extensions of the MOBE-2 model.

5. Estimation and Inference

Let X 1 , X 2 , , X n be a random sample of size n from MOBE-2 ( Ψ _ ) where Ψ _ = ( α , β , γ ) . The log likelihood function ( l o g L ( Ψ _ ) ) can be written as:
l o g L ( Ψ _ ) = n l o g ( α ) + n l o g ( γ ) n l o g ( 2 β ) + i = 0 n l o g [ 2 2 β + β α x i ] α i = 0 n x i 2 i = 0 n l o g [ 1 γ ¯ ( 1 + β α x i 2 β ) e α x i ] .
The associated score function is given by U n ( Ψ _ ) = ( l o g L ( Ψ _ ) α , l o g L ( Ψ _ ) β , l o g L ( Ψ _ ) γ ) T . The log-likelihood can be maximized either directly or by solving the nonlinear likelihood equations obtained by differentiating l o g L ( Ψ _ ) . The components of the score vector are given by:
l o g L ( Ψ _ ) α = n α + i = 0 n β x i 2 2 β + β α x i i = 0 n x i + 2 γ ¯ i = 0 n β x i e α x i 2 β x i e α x i ( 1 + β α x i 2 β ) 1 γ ¯ ( 1 + β α x i 2 β ) e α x i ,
l o g L ( Ψ _ ) β = n 2 β + i = 0 n α x i 2 2 2 β + β α x i + 4 i = 0 n γ ¯ α x i ( 2 β ) 2 [ 1 γ ¯ ( 1 + β α x i 2 β ) e α x i ] ,
and:
l o g L ( Ψ _ ) γ = n γ 2 i = 0 n ( 1 + β α x i 2 β ) e α x i 1 γ ¯ ( 1 + β α x i 2 β ) e α x i .

6. Modeling

In this Section, two real data applications are presented for illustrating the importance and flexibility of the new model. The first data set (1.1, 1.4, 1.3, 1.7, 1.9, 1.8, 1.6, 2.2, 1.7, 2.7, 4.1, 1.8, 1.5, 1.2, 1.4, 3, 1.7, 2.3, 1.6, 2), called the failure time data, represents the lifetime data relating to relief times (in minutes) of patients receiving an analgesic (for more applications to this data see [6,7,8,9,10,11,12]).
The second data set (0.1, 0.33, 0.44, 0.56, 0.59, 0.72, 0.74, 0.77, 0.92, 0.93, 0.96, 1, 1, 1.02, 1.05, 1.07, 07, 1.08, 1.08, 1.08, 1.09, 1.12, 1.13, 1.15, 1.16, 1.2, 1.21, 1.22, 1.22, 1.24, 1.3, 1.34, 1.36, 1.39, 1.44, 1.46, 1.53, 1.59, 1.6, 1.63, 1.63, 1.68, 1.71, 1.72, 1.76, 1.83, 1.95, 1.96, 1.97, 2.02, 2.13, 2.15, 2.16, 2.22, 2.3, 2.31, 2.4, 2.45, 2.51, 2.53, 2.54, 2.54, 2.78, 2.93, 3.27, 3.42, 3.47, 3.61, 4.02, 4.32, 4.58, 5.55) represents the survival times (in days) of 72 guinea pigs infected with virulent tubercle bacilli, observed and reported by Bjerkedal (1960). Many other real data sets related to failure times can be found in [13,14,15,16,17,18,19,20].
Figure 3 gives the total time test (TTT) plots. Figure 4 and Figure 5 gives the estimated PDF (EPDF), estimated survival function (ESF), P-P plots, and estimated CDF (ECDF) for the two data sets, respectively. From Figure 3, we note that the empirical HRF is increasing for the two data sets.
We compared the fits of the MOBE-2 distribution with some competitive models, namely: exponential (E ( β ) ), odd Lindley exponential (OLiE), MO exponential (MOE ( α , β ) ), moment exponential (MomE ( β ) ), the logarithmic Burr–Hatke exponential (Log BrHE ( β ) ), generalized MO exponential (GMOE ( α , α , β ) ), beta exponential (BE ( a , b , β ) ), MO–Kumaraswamy exponential (MOKwE ( α , a , b , β ) ), Kumaraswamy exponential (KwE ( a , b , β ) ), and Kumaraswamy MO exponential (KwMOE ( α , a , b , β ) ). See the PDFs of the competitive moels in [21,22,23,24,25,26,27,28,29,30,31]. We considered the Cramér-Von Mises ( W * ), the Anderson–Darling ( A * ), and the Kolmogorov–Smirnov (KS) statistics. The W * and A * statistics are given by:
W * = ( 1 + 1 / 2 n ) [ [ 1 / ( 12 n ) ] + κ = 1 n ω κ ] ,
and:
A * = a ( n ) ( n + n 1 κ = 1 n a κ ) ,
where:
ω κ = [ z i ( 2 κ 1 ) / ( 2 n ) ] 2 ,   a ( n ) = 1 + 9 4 n 2 + 3 4 n 1 ,
and:
a κ = ( 2 κ 1 ) l o g [ z i ( 1 z n 1 + κ ) ] ,
where z i = F ( y κ ) and y κ ’s values are the ordered observations. Moreover, we considered some other goodness-of-fit measures, including the Akaike Information Criterion (AIC), Consistent Akaike Information Criterion (CAIC), Hannan–Quinn Information Criterion (HQIC), and Bayesian Information Criterion (BIC).
Table 1 gives the maximum likelihood estimations (MLE) and SE values for the relief times data. Table 2 give the AIC, BIC, CAIC, and HQIC for the relief times data. Table 3 gives A * , W * , KS, and p-value for the relief times data. Table 4 gives the MLE and SE values for the survival times data. Table 5 gives the AIC, BIC, CAIC, and HQIC for the survival times data. Table 6 gives the A * , W * , KS, and p-value for the survival times data.
Based on Table 2, Table 3, Table 5 and Table 6, the proposed lifetime MOBE-2 model is much better than many competitive models, such as the E, MomE, MOE, GMOE, KwE, BE, MOKE, and KMOE models, so the new lifetime model may be a good alternative to these models in modeling relief times and survival times data sets. From Figure 5 and Figure 6, we note that the MOBE-2 model gives an adequate fit with the two real data sets.

7. Concluding Remarks

In this paper, we introduced a new version of the BE-2 model. The new model is called MOBE-2 distribution. Some of its structural properties are also presented. A simple type Copula-based construction was also presented to construct the bivariate and the multivariate type distributions. We illustrated the importance of the new version by the study of two real data applications. The proposed lifetime model was much better than many competitive versions, such as the exponential, the odd Lindley exponential, the Marshall–Olkin exponential, moment exponential, the logarithmic Burr–Hatke exponential, the generalized Marshall–Olkin exponential, beta exponential, the Marshall–Olkin–Kumaraswamy exponential, the Kumaraswamy exponential, and the Kumaraswamy–Marshall–Olkin Exponential, so the new lifetime model may be a good alternative to these models in modeling relief times and survival times data sets.

Author Contributions

A.A.A.-b. and I.E.: writing, conceptualization, and methodology. H.M.Y.: writing—review, editing, and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by King Saud University, Deanship of Scientific Research, and College of Sciences Research Center.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this Appendix, we derive a useful expansion for the MOBE-2 density function as follows. Using the generalized binomial expression
( 1 z ) c = κ = 0 Γ ( c + κ ) Γ ( c ) κ ! z κ |   ( | z | < 1   and   c > 0   real   non - integer ) .
Then using (12), the PDF of the MOBE-2 can be written as
f Ψ _ ( x ) = γ α κ = 0 γ ¯ κ ( 1 + κ ) ( 1 + β α x 2 β ) κ [ 1 + ( α x 1 ) β 2 β ] e α ( 1 + κ ) x ,
applying the binomial expression for ( 1 + β α x 2 β ) κ , the new PDF becomes
f Ψ _ ( x ) = γ α κ = 0 γ ¯ κ ( 1 + κ ) τ = 0 κ ( β α 2 β ) τ x τ [ 1 + ( α x 1 ) β 2 β ] e α ( 1 + κ ) x ,
after some simplification, MOBE-2 density function can be expressed as
f Ψ _ ( x ) = κ = 0 τ = 0 κ w τ , κ   x τ [ 2 ( 1 β ) + β α x ] e α ( 1 + κ ) x ,
where
w τ , κ = γ α ( 1 + κ ) γ ¯ κ ( β α ) τ ( 2 β ) τ + 1 .

References

  1. Shaked, M.; Shanthikumar, J. Stochastic Orders; Springer: New York, NY, USA, 2007. [Google Scholar]
  2. Barlow, R.E.; Proschan, F. Statistical Theory of Reliability and Life Testing; Holt, Rinehart, and Winston: New York, NY, USA, 1975. [Google Scholar]
  3. Bakouch, S.; Jazi, M.; Nadarajah, S.; Dolati, A.; Roozegar, R. A lifetime model with increasing failure rate. Appl. Math. Model. 2014, 38, 5392–5406. [Google Scholar] [CrossRef]
  4. 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. Biometrika 1997, 84, 641–652. [Google Scholar] [CrossRef]
  5. Navarro, J.; Franco, M.; Ruiz, J.M. Characterization through moments of the residual life and conditional spacing. Indian J. Stat. Ser. A 1998, 60, 36–48. [Google Scholar]
  6. Alizadeh, M.; Ghosh, I.; Yousof, H.M.; Rasekhi, M.; Hamedani, G.G. The generalized odd generalized exponential family of distributions: Properties, characterizations and applications. J. Data Sci. 2017, 15, 443–466. [Google Scholar]
  7. Yousof, H.M.; Afify, A.Z.; Alizadeh, M.; Nadarajah, S.; Aryal, G.R.; Hamedani, G.G. The Marshall-Olkin generalized-G family of distributions with Applications. Statistica 2018, 78, 273–295. [Google Scholar]
  8. Yousof, H.M.; Altun, E.; Hamedani, G.G. A new extension of Frechet distribution with regression models, residual analysis and characterizations. J. Data Sci. 2018, 16, 743–770. [Google Scholar]
  9. Yousof, H.M.; Afify, A.Z.; Alizadeh, M.; Hamedani, G.G.; Jahanshahi, S.M.A.; Ghosh, I. The generalized transmuted Poisson-G family of Distributions. Pak. J. Stat. Oper. Res. 2018, 14, 759–779. [Google Scholar] [CrossRef]
  10. Yousof, H.M.; Altun, E.; Ramires, T.G.; Alizadeh, M.; Rasekhi, M. A new family of distributions with properties, regression models and applications. J. Stat. Manag. Syst. 2018, 21, 163–188. [Google Scholar] [CrossRef]
  11. Yousof, H.M.; Butt, N.S.; Alotaibi, R.; Rezk, H.; Alomani, G.A.; Ibrahim, M. A new compound Fréchet distribution for modeling breaking stress and strengths data. Pak. J. Stat. Oper. Res. 2019, 15, 1017–1035. [Google Scholar] [CrossRef]
  12. Yousof, H.M.; Rasekhi, M.; Altun, E.; Alizadeh, M. The extended odd Frechet family of distributions: properties, applications and regression modeling. Int. J. Appl. Math. Stat. 2018, 30, 1–30. [Google Scholar]
  13. Nascimento, A.D.C.; Silva, K.F.; Cordeiro, G.M.; Alizadeh, M.; Yousof, H.M. The odd Nadarajah-Haghighi family of distributions: Properties and applications. Stud. Scientiarum Math. Hungarica 2019, 56, 1–26. [Google Scholar] [CrossRef]
  14. Bantan, R.A.R.; Jamal, F.; Chesneau, C.; Elgarhy, M. Truncated inverted Kumaraswamy generated family of distributions with applications. Entropy 2019, 21, 1089. [Google Scholar] [CrossRef] [Green Version]
  15. Arshad, R.M.I.; Chesneau, C.; Jamal, F. The odd gamma Weibull-geometric model: Theory and applications. Mathematics 2019, 7, 399. [Google Scholar] [CrossRef] [Green Version]
  16. Korkmaz, M.C.; Altun, E.; Yousof, H.M.; Hamedani, G.G. The odd power Lindley generator of probability distributions: Properties, characterizations and regression modeling. Int. J. Stat. Probab. 2019, 8, 70–89. [Google Scholar] [CrossRef]
  17. Korkmaz, M.C.; Yousof, H.M.; Hamedani, G.G. The exponential Lindley odd log-logistic G family: properties, characterizations and applications. J. Stat. Theory Appl. 2018, 17, 554–571. [Google Scholar] [CrossRef] [Green Version]
  18. Elsayed, H.A.H.; Yousof, H.M. The Burr X Nadarajah Haghighi distribution: Statistical properties and application to the exceedances of flood peaks data. J. Math. Stat. 2019, 15, 146–157. [Google Scholar] [CrossRef]
  19. Merovci, F.; Alizadeh, M.; Yousof, H.M.; Hamedani, G.G. The exponentiated transmuted-G family of distributions: theory and applications. Commun. Stat.-Theory Methods 2017, 46, 10800–10822. [Google Scholar] [CrossRef]
  20. Korkmaz, M.C.; Yousof, H.M.; Ali, M.M. Some theoretical and computational aspects of the odd Lindley Fréchet distribution. J. Stat.: Stat. Actuarial Sci. 2017, 2, 129–140. [Google Scholar]
  21. Ibrahim, M.; Yadav, A.S.; Yousof, H.M.; Goual, H.; Hamedani, G.G. A new extension of Lindley distribution: modified validation test, characterizations and different methods of estimation. Commun. Stat. Appl. Methods 2019, 26, 473–495. [Google Scholar] [CrossRef]
  22. Idika, E.; Anthony, C.; Johnson, O. Marshall–Olkin generalized Erlange-truncated exponential distribution: Properties and Applications. Cogent Math. 2017, 4, 1–19. [Google Scholar]
  23. Hamedani, G.G.; Altun, E.; Korkmaz, M.C.; Yousof, H.M.; Butt, N.S. A new extended G family of continuous distributions with mathematical properties, characterizations and regression modeling. Pak. J. Stat. Oper. Res. 2018, 14, 737–758. [Google Scholar] [CrossRef] [Green Version]
  24. Alizadeh, M.; Rasekhi, M.; Yousof, H.M.; Hamedani, G.G. The transmuted Weibull G family of distributions. Hacet. J. Math. Stat. 2018, 47, 1–20. [Google Scholar] [CrossRef]
  25. Yadav, A.S.; Goual, H.; Alotaibi, R.M.; Rezk, H.; Ali, M.M.; Yousof, H.M. Validation of the Topp-Leone-Lomax model via a modified Nikulin-Rao-Robson goodness-of-fit test with different methods of estimation. Symmetry 2020, 12, 57. [Google Scholar] [CrossRef] [Green Version]
  26. Hamedani, G.G.; Rasekhi, M.; Najibi, S.M.; Yousof, H.M.; Alizadeh, M. Type II general exponential class of distributions. Pak. J. Stat. Oper. Res. 2019, 15, 503–523. [Google Scholar] [CrossRef]
  27. Hamedani, G.G.; Yousof, H.M.; Rasekhi, M.; Alizadeh, M.; Najibi, S.M. Type I general exponential class of distributions. Pak. J. Stat. Oper. Res. 2017, 14, 39–55. [Google Scholar] [CrossRef] [Green Version]
  28. Ibrahim, M. A new extended Fréchet distribution: Properties and estimation. Pak. J. Stat. Oper. Res. 2019, 15, 773–796. [Google Scholar]
  29. Alizadeh, M.; Yousof, H.M.; Rasekhi, M.; Altun, E. The odd log-logistic Poisson-G Family of distributions. J. Math. Ext. 2019, 12, 81–104. [Google Scholar]
  30. Ghitany, M.E.; Al-Awadhi, F.A.; Alkhalfan, L.A. Marshall-Olkin extended Lomax distribution and its application to censored data. Commun. Statist. Theor. Meth. 2007, 36, 1855–1866. [Google Scholar] [CrossRef]
  31. Aryal, G.R.; Yousof, H.M. The exponentiated generalized-G Poisson family of distributions. Econ. Qual. Control 2017, 32, 1–17. [Google Scholar] [CrossRef]
Figure 1. Plots of the new PDF based on some selected parameters values.
Figure 1. Plots of the new PDF based on some selected parameters values.
Symmetry 12 00440 g001
Figure 2. Plots of the new PDF based on some selected parameters values.
Figure 2. Plots of the new PDF based on some selected parameters values.
Symmetry 12 00440 g002
Figure 3. 3-D plot for skewness and kurtosis of the new model when γ = 2.25 .
Figure 3. 3-D plot for skewness and kurtosis of the new model when γ = 2.25 .
Symmetry 12 00440 g003
Figure 4. The total time test TTT plots.
Figure 4. The total time test TTT plots.
Symmetry 12 00440 g004
Figure 5. Estimated PDF (EPDF), estimated survival function (ESF), Probability-Probability (P-P) plot, and estimated cumulative distribution function (ECDF) for the relief times data.
Figure 5. Estimated PDF (EPDF), estimated survival function (ESF), Probability-Probability (P-P) plot, and estimated cumulative distribution function (ECDF) for the relief times data.
Symmetry 12 00440 g005
Figure 6. Estimated PDF (EPDF), estimated survival function (ESF), Probability-Probability (P-P) plot, and estimated cumulative distribution function (ECDF) for the survival times data.
Figure 6. Estimated PDF (EPDF), estimated survival function (ESF), Probability-Probability (P-P) plot, and estimated cumulative distribution function (ECDF) for the survival times data.
Symmetry 12 00440 g006
Table 1. MLE and SE values for the relief times data.
Table 1. MLE and SE values for the relief times data.
ModelsEstimates
E(β)0.526
(0.117)
OLiE(β)0.6044
(0.0535)
MomE(β)0.950
(0.150)
Log BrHE(β)0.5263
(0.118)
MOE (α, β)54.474, 2.316
(35.582), (0.374)
GMOE (α, α, β)0.519, 89.462, 3.169
(0.256), (66.278), (0.772)
KwE(a,b,β)83.756, 0.568, 3.330
(42.361), (0.326), (1.188)
BE(a,b,β)81.633, 0.542, 3.514
(120.41), (0.327), (1.410)
MOKwE (α, a,b,β)0.133, 33.232, 0.571, 1.669
(0.332), (57.837), (0.721), (1.814)
KwMOE (α, a,b,β)8.868, 34.826, 0.299, 4.899
(9.146), (22.312), (0.239), (3.176)
BrXE (α, β)1.1635, 0.3207
(0.33), (0.03)
MOBE2 (γ, α, β)1.83 × 103, 6.707 × 10−2, 6.096 × 10−3
(2.206 × 103), (4.991 × 10−3), (1.069 × 10−3)
Table 2. Akaike Information Criterion (AIC), Consistent Akaike Information Criterion (CAIC), Hannan–Quinn Information Criterion (HQIC), and Bayesian Information Criterion (BIC) for the relief times data.
Table 2. Akaike Information Criterion (AIC), Consistent Akaike Information Criterion (CAIC), Hannan–Quinn Information Criterion (HQIC), and Bayesian Information Criterion (BIC) for the relief times data.
ModelsAIC, BIC, CAIC, HQIC
MOBE240.1, 40.2, 40.3, 39.1
E67.67, 68.67, 67.89, 67.87
OLiE49.1, 50.1, 49.3, 49.3
MomE54.32, 55.31, 54.54, 54.50
Log BrHE67.67,68.67,67.89,67.87
MOE43.51, 45.51, 44.22, 43.90
GMOE42.75, 45.74, 44.25, 43.34
KwE41.78, 44.75, 43.28, 42.32
BE43.48, 46.45, 44.98, 44.02
MOKwE41.58, 45.54, 44.25, 42.30
KwMOExp42.8, 46.84, 45.55, 43.60
BrXE48.1, 50.1, 8.8, 48.5
Table 3. Cramér-Von Mises ( W * ), the Anderson–Darling ( A * ), and the Kolmogorov–Smirnov (KS) statistics, and p-value for the relief times data.
Table 3. Cramér-Von Mises ( W * ), the Anderson–Darling ( A * ), and the Kolmogorov–Smirnov (KS) statistics, and p-value for the relief times data.
Models A * ,   W * , KS and p-Value
MOBE20.33, 0.046, 0.12(0.95)
E4.60, 0.96, 0.44(0.004)
OLiE1.3, 0.22, 0.85(6.23 × e−13)
MomE2.76,0.53, 0.32(0.07)
Log BrHE0.62, 0.105, 0.44(0.0009)
MOE0.8, 0.14, 0.1(0.55)
GMOE0.51, 0.08, 0.15(0.78)
KwE0.45, 0.07, 0.14(0.86)
BE0.70, 0.12, 0.16(0.80)
MOKwE0.60, 0.11, 0.14(0.87)
KwMOExp1.08, 0.19, 0.15(0.86)
BrXE1.39, 0.24, 0.248(0.1705)
Table 4. MLEs and SEs values for the survival times data.
Table 4. MLEs and SEs values for the survival times data.
ModelsEstimates
E(β)0.540
(0.063)
OLiE(β)0.38145
(0.0209)
MomE(β)0.925
(0.077)
Log BrHE(β)0.54
(0.064)
MOE (α, β)8.778, 1.379
(3.555), (0.193)
GMOE (α, α, β)0.179, 47.635, 4.465
(0.070), (44.901), (1.327)
KwE (a, b, β)3.304, 1.100, 1.037
(1.106), (0.764), (0.614)
BE (a, b, β)0.807, 3.461, 1.331
(0.696), (1.003), (0.855)
MOKwE (α, a, b, β)0.008, 2.716, 1.986, 0.099
(0.002), 1.316), (0.784), (0.048)
KwMOE (α, a, b, β)0.373, 3.478, 3.306, 0.299
(0.136), (0.861), (0.779), (1.112)
BrXE (α, β)0.475, 0.2055
(0.06), (0.012)
MOBE2(γ, α, β)11.0365, 0.12054, 0.013601
(4.8066), (0.02246), (0.0077)
Table 5. AIC, BIC, CAIC, and HQIC for the survival times data.
Table 5. AIC, BIC, CAIC, and HQIC for the survival times data.
ModelsAIC, BIC, CAIC, HQIC
MOBE2207.3, 213.15, 206.6, 209.01
E234.63, 236.91, 234.68, 235.54
OLiE229.1, 231.4, 229.2, 230
MomE210.40, 212.68, 210.45, 211.30
Log BrHE234.63, 236.9, 234.7, 235.5
MOE210.36, 214.92, 210.53, 212.16
GMOE210.54, 217.38, 210.89, 213.24
KwE209.42, 216.24, 209.77, 212.12
BE207.38, 214.22, 207.73, 210.08
MOKwE209.44, 218.56, 210.04, 213.04
KwMOE207.82, 216.94, 208.42, 211.42
BrXE235.3, 239.9, 235.5, 237.1
Table 6. A * , W * , KS, and p-value for the survival times data.
Table 6. A * , W * , KS, and p-value for the survival times data.
Models A * ,   W * , KS, p-Value
MOBE20.68, 0.09, 0.089(0.64)
E6.53, 1.25, 0.27(0.06)
OLiE1.94, 0.33, 0.49(9.992 × e−16)
MomE1.52, 0.25, 0.14(0.13)
Log BrHE0.71, 0.115, 0.28(2.382 × e−5)
MOE1.18, 0.17, 0.1(0.43)
GMOE1.02, 0.16, 0.09(0.51)
KwE0.74, 0.11, 0.09(0.50)
BE0.98, 0.15, 0.11(0.34)
MOKwE0.79, 0.12, 0.10(0.44)
KwMOE0.61, 0.11, 0.09(0.53)
BrXE2.9, 0.52, 0.22(0.002)

Share and Cite

MDPI and ACS Style

A. Al-babtain, A.; Elbatal, I.; M. Yousof, H. A New Flexible Three-Parameter Model: Properties, Clayton Copula, and Modeling Real Data. Symmetry 2020, 12, 440. https://doi.org/10.3390/sym12030440

AMA Style

A. Al-babtain A, Elbatal I, M. Yousof H. A New Flexible Three-Parameter Model: Properties, Clayton Copula, and Modeling Real Data. Symmetry. 2020; 12(3):440. https://doi.org/10.3390/sym12030440

Chicago/Turabian Style

A. Al-babtain, Abdulhakim, I. Elbatal, and Haitham M. Yousof. 2020. "A New Flexible Three-Parameter Model: Properties, Clayton Copula, and Modeling Real Data" Symmetry 12, no. 3: 440. https://doi.org/10.3390/sym12030440

APA Style

A. Al-babtain, A., Elbatal, I., & M. Yousof, H. (2020). A New Flexible Three-Parameter Model: Properties, Clayton Copula, and Modeling Real Data. Symmetry, 12(3), 440. https://doi.org/10.3390/sym12030440

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