Next Article in Journal
Double-Diffusive Convection in Bidispersive Porous Medium with Coriolis Effect
Next Article in Special Issue
Prony Method for Two-Generator Sparse Expansion Problem
Previous Article in Journal
Learning Motion Primitives Automata for Autonomous Driving Applications
Previous Article in Special Issue
A Note on Gerber–Shiu Function with Delayed Claim Reporting under Constant Force of Interest
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Generalized Odd Linear Exponential Family of Distributions with Applications to Reliability Theory

1
Department of Statistics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
2
Department of Statistics, Darrang College, Gauhati University, Tezpur 784001, India
3
Department of Mathematical Statistics, Institute of Statistical Studies and Research, Cairo University, Giza 12613, Egypt
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2022, 27(4), 55; https://doi.org/10.3390/mca27040055
Submission received: 20 May 2022 / Revised: 18 June 2022 / Accepted: 21 June 2022 / Published: 23 June 2022
(This article belongs to the Special Issue Computational Mathematics and Applied Statistics)

Abstract

:
A new family of continuous distributions called the generalized odd linear exponential family is proposed. The probability density and cumulative distribution function are expressed as infinite linear mixtures of exponentiated-F distribution. Important statistical properties such as quantile function, moment generating function, distribution of order statistics, moments, mean deviations, asymptotes and the stress–strength model of the proposed family are investigated. The maximum likelihood estimation of the parameters is presented. Simulation is carried out for two of the mentioned sub-models to check the asymptotic behavior of the maximum likelihood estimates. Two real-life data sets are used to establish the credibility of the proposed model. This is achieved by conducting data fitting of two of its sub-models and then comparing the results with suitable competitive lifetime models to generate conclusive evidence.

1. Introduction

Analysis of lifetime data is an important subject in many fields, including reliability, social sciences, biomedical, engineering and other fields. In practice, it has been observed that many phenomena do not follow any of the classical distributions; for this reason, many efforts have been made in the last few decades to introduce new generators or families of distributions that extend these classical distributions to provide considerable flexibility in modeling data in diverse spectrums. Many authors have suggested new generators or families in the literature, for example, and not exclusively: Marshall and Olkin (1997) [1] introduced the Marshall–Olkin family, Gupta et al. (1998) [2] introduced the exponentiated-G family, Eugene et al. (2002) [3] proposed the beta-G family, Cordeiro and Castro (2011) [4] suggested the Kumaraswamy-G family, Alexander et al. (2012) [5] presented the McDonald-G family, Alzaatreh et al. (2013) [6] proposed the transformed-transformer (T-X) family, Bourguignon et al. (2014) [7] presented the Weibull-G family, Tahir et al. (2015) [8] studied the odd generalized exponential-G family, Cordeiro et al. (2016) [9] discussed the Zografos Balakrishnan odd log-logistic family, Gomes-Silva et al. (2017) [10] presented the odd Lindley-G family, Alizadeh et al. (2017) [11] provided the Gompertz-G family and Jamal et al. (2017) [12] defined the odd Burr-III family, among others. For a clearer understanding of the odds ratio to define new G-classes, we motivate the readers to Khan et al. (2021) [13], in which the authors adopted a unique odd function to propose an alternate generalized odd generalized exponential-G family.
The linear exponential or (linear failure rate) distribution is the distribution of the minimum of two independent random variables Z1 and Z2 having exponential (a) and Rayleigh (b) (Sen and Bhattacharyya, 1995 [14]). Therefore, the variables have exponential and Rayleigh distributions as special cases, which are well-known distributions for modeling lifetime data in reliability and medical studies. The linear exponential distribution is used to model phenomena with linearly increasing failure rates, but it does not provide a reasonable fit for modeling phenomena with decreasing, non-linear increasing, or non-monotonic failure rates, which include the bathtub and upside-down bathtub, among others. These phenomena are common in reliability and biological studies. This motivated us to introduce generalizations of linear, exponential distribution so that their goodness of fit measures may improve the tail properties. Our motivations and the main goals of this paper are to propose a random variable that follows the linear exponential distribution as a new generator to introduce new models which can yield all types of the hazard rate functions with improved goodness of fit properties for real-life data.

2. The Generalized Odd Linear Exponential (GOLE-F) Family

Suppose the random variable Z has a linear exponential distribution with parameters a , b 0 where a + b > 0 , then its cumulative distribution function (CDF) and probability density function (PDF) are, respectively,
R ( z ) = 1 e ( a z + b 2 z 2 ) ,   z 0  
r ( z ) = ( a + b z ) e ( a z + b 2 z 2 ) .   z > 0 .
Adopting the T-X framework defined by the authors in [6], for any power parameter c > 0 , we define the CDF of a new wider family called the generalized odd linear exponential (“GOLE-F” for short) family by
G ( x ; a , b , c , ϕ ) = 0 F ( x ; ϕ ) c 1 F ( x ; ϕ ) c ( a + b z ) e ( a z + b 2 z 2 ) d z = 1 exp [ ( a F ( x ; ϕ ) c 1 F ( x ; ϕ ) c + b 2 ( F ( x ; ϕ ) c 1 F ( x ; ϕ ) c ) 2 ) ] ,  
where W [ F ( x ) ] = F ( x ; ϕ ) c 1 F ( x ; ϕ ) c is the link function with F ( x ; ϕ ) as the baseline CDF of an absolutely continuous distribution with parameter vector ϕ and pdf f ( x ; ϕ ) .
The PDF of GOLE-F corresponding to the CDF in Equation (3) is provided by
g ( x ; a , b , c , ϕ ) = [ c f ( x ; ϕ ) F ( x ; ϕ ) c 1 ( a + ( b a ) F ( x ; ϕ ) c ) ( 1 F ( x ; ϕ ) c ) 3 ] × exp [ ( a F ( x ; ϕ ) c 1 F ( x ; ϕ ) c + b 2 ( F ( x ; ϕ ) c 1 F ( x ; ϕ ) c ) 2 ) ] .
Henceforth, for any parent model, we will simply write F ( x ) = F ( x ; ϕ ) as the distribution function and f ( x ) = f ( x ; ϕ ) as the density function. Further, any random variable X with density function (4) is denoted by X ~ G O L E F   ( a , b , c , ϕ ) .
The hazard rate function (HRF) and reversed hazard rate function (RHRF) of the random variable X are, respectively,
h ( x ; a , b , c , ϕ ) = c f ( x ) F ( x ) c 1 ( a + ( b a ) F ( x ) c ) ( 1 F ( x ) c ) 3 ,  
and
τ ( x ; a , b , c , ϕ ) = { c f ( x ) F ( x ) c 1 ( a + ( b a ) F ( x ) c ) } e ( a F ( x ) c 1 F ( x ) c + b 2 ( F ( x ) c 1 F ( x ) c ) 2 ) ( 1 F ( x ) c ) 3 { 1 e ( a F ( x ) c 1 F ( x ) c + b 2 ( F ( x ) c 1 F ( x ) c ) 2 ) } .  
The quantile function of the random variable X can be obtained by inverting Equation (3), and hence the GOLE-F distribution can be simulated easily from the following Equation.
X = Q ( U ) = F 1 ( [ a + a 2 2 b log ( 1 U ) b a + a 2 2 b log ( 1 U ) ] 1 / c ) ,  
where U has a uniform distribution over the interval (0,1), in particular, if u = 1 / 2 we obtain the median of the random variable X as follows:
M = Q ( 1 2 ) = F 1 ( [ a + a 2 2 b log ( 1 1 / 2 ) b a + a 2 2 b log ( 1 1 / 2 ) ] 1 / c ) .  

3. Special Model of the GOLE-F Family

In this section, we provide two extended distributions as special models of the GOLE-F family and display their plots of density and hazard rate functions.

3.1. The Generalized Odd Linear Exponential-Weibull (GOLE-W) Distribution

Consider the Weibull distribution with density and distribution functions f ( x ; λ , β ) = λ β x β 1 e λ x β and F ( x ; λ , β ) = 1 e λ x β , respectively, where λ , β > 0 and x 0 . Then, the GOLE-W distribution has (PDF) provided by
g ( x ; a , b , c , λ , β ) = [ c λ β x β 1 e λ x β ( 1 e λ x β ) c 1 ( a + ( b a ) ( 1 e λ x β ) c ) ( 1 ( 1 e λ x β ) c ) 3 ] × exp [ ( a ( 1 e λ x β ) c 1 ( 1 e λ x β ) c + b 2 ( ( 1 e λ x β ) c 1 ( 1 e λ x β ) c ) 2 ) ] .
Figure 1a show a wealth of possible shapes of the distribution once different choices of the parameters are made. For example, the shape can be U and inverted-U, right-skewed, reversed-J shape or symmetrical. Additionally, Figure 1b reveal that the HRF of the GOLE-W distribution can be increasing–constant, constant–monotone–increasing or monotone–increasing shapes.

3.2. The Generalized Odd Linear Exponential-Exponential (GOLE-E) Distribution

Consider the Exponential distribution with density and distribution functions f ( x ; λ ) = λ e λ x and F ( x ; λ ) = 1 e λ x , respectively, where λ > 0 and x 0 . Then, the GOLE-E distribution has (PDF) provided by
g ( x ; a , b , c , λ ) = [ c λ e λ x ( 1 e λ x ) c 1 ( a + ( b a ) ( 1 e λ x ) c ) ( 1 ( 1 e λ x ) c ) 3 ] × exp [ ( a ( 1 e λ x ) c 1 ( 1 e λ x ) c + b 2 ( ( 1 e λ x ) c 1 ( 1 e λ x ) c ) 2 ) ] .
Figure 2a show possible shapes of the GOLE-E distribution for different choices of the parameters. The shapes of pdf can be right-skewed, or symmetrical. Further, Figure 2b reveal that the HRF of the GOLE-E distribution can be decreasing–constant, monotone–increasing or bathtub shape. The PDF and HRF of the GOLE-W and GOLE-E distributions for some selected values of the parameters indicate the flexibility of the new family.

4. Mathematical Properties of the GOLE-F Family

In this section, some mathematical properties of the GOLE-F family are obtained.

4.1. Asymptotic Behavior of GOLE-F Family

First of all, for the statements of the following results, we recall that F ( x ) is the CDF of an absolutely continuous distribution with pdf f ( x ) .
Proposition 1.
The asymptotes corresponding to Equations (3)–(5) when x are provided by
G ( x ) ~ a   F ( x ) c ,
g ( x ) ~ c   a   f ( x ) F ( x ) c 1 ,
h ( x ) ~ c   a   f ( x ) F ( x ) c 1 .  
Proposition 2.
The asymptotes corresponding to Equations (3)–(5) when x are provided by
1 G ( x )   ~   1 e ( a 1 F ( x ) c   + b 2 { 1 1 F ( x ) c   } 2 ) ,
g ( x ) ~ b c f ( x ) ( 1 F ( x ) c ) 3     e ( a 1 F ( x ) c   + b 2 { 1 1 F ( x ) c   } 2 )  
h ( x ) ~ b c f ( x ) ( 1 F ( x ) c ) 3 .  
For detail seeAppendix A.

4.2. Useful Expansions for CDF and PDF of the New Family

Using the power series for the exponential function and the generalized binomial expansion
e z = i = 0 ( 1 ) i z i i ! ,
and
( 1 v ) n = i = 0 ( 1 ) i ( n i ) v i ,
respectively, where | v | < 1 and n is any real number, we can rewrite the CDF of the GOLE-F family as follows:
G ( x ; a , b , c , ϕ ) = 1 i = 0 j = 0 i k = 0 ( 1 ) i b j a i ( i j ) ( i + j + k 1 k ) i ! 2 j F ( x ) c ( i + j + k ) .
Again, based on the binomial expansion, we find
F ( x ) c ( i + j + k ) = ( 1 ( 1 F ( x ) ) c ( i + j + k ) = m = 0 l = m ( 1 ) l + m ( l m ) ( c ( i + j + k ) l ) F ( x ) m .  
From (15) and (16), we obtain
G ( x ; a , b , c , ϕ ) = 1 m = 0 ω m F ( x ) m ,  
where
ω m = l = m i = 0 j = 0 i k = 0 ρ i , j , k , l ( a , b , c ) ,
and
ρ i , j , k , l ( a , b , c ) = ( 1 ) i + l + m b j a i ( i j ) ( l m ) ( i + j + k 1 k ) ( c ( i + j + k ) l ) i ! 2 j .
Now, we can write the CDF of the GOLE-F family in Equation (17), as
G ( x ; a , b , c , ϕ ) = m = 0 δ m F ( x ) m ,
where δ 0 = 1 ω 0 , and δ m = ω m for m = 1 , 2 , By differentiating Equation (18), we obtain the expansion of the density function of the GOLE-F family as an infinite linear mixture of exp-F densities in the following form
g ( x ; a , b , c , ϕ ) = m = 0 δ m + 1 π m + 1 ( x ) ,
where π m + 1 ( x ) = ( m + 1 ) f ( x ) F ( x ) m is the exp-F density function with power parameter ( m + 1 ) . Now, if the random variable Y m + 1 has the density function π m + 1 ( x ) , then many mathematical properties of the random variable X , including the ordinary and incomplete moments and moment generating function can easily be obtained based on the exp-F distribution.

4.3. Moments

Suppose that the random variable Y m + 1 has the density function π m + 1 ( x ) in (19), then the nth moment of the random variable X can be obtained from
μ n = E ( X n ) = m = 0 δ m + 1 E ( Y m + 1 n ) .
A second alternative formula for μ n   in terms of the baseline qf. Q F ( u ) can be obtained as
μ n = m = 0 δ m + 1 ( m + 1 ) 0 1 Q F ( u ) n u m d u ,  
where Q F ( u ) = F 1 ( u ) is the qf of the parent distribution and u ( 0 , 1 ) .
The incomplete moments have an important role in measuring inequality, for example, income quantiles, the mean deviations and Lorenz and Bonferroni curves. The nth incomplete moment of X is provided by
η n ( z ) = m = 0 δ m + 1 ( m + 1 ) 0 F ( z ) Q F ( u ) n u m d u .
The last integral can be computed analytically or numerically for most baseline distributions. Bonferroni and Lorenz curves have applications in many different areas such as economics to study income and poverty, reliability, demography, insurance and medicine. For a random variable X , the Bonferroni and Lorenz curves are defined by B ( p ) = η 1 ( q ) / p E ( X ) and L ( p ) = η 1 ( q ) / E ( X ) , respectively, where p is a given probability, q = Q ( p ) and η 1 ( q ) is the first incomplete moment that can be calculated from the above Equation with r = 1 at q . Table 1 display the mean, variance, skewness and kurtosis of the GOLE-E distribution for some choices values of the parameters. We note from Table 1 that the skewness of the GOLE-E distribution is always positive, whereas the kurtosis of the GOLE-E distribution varies only in the interval (1.0571, 2.6112).

4.4. Generating Function

Here, we provide three formulae for the mgf M ( t ) = E ( e t X ) of the random variable X . The first one is provided by
M ( t ) = n = 0 t n n ! μ n ,  
where μ n   is the nth moment of the random variable X . A second formula for M ( t ) comes from (19) as
M ( t ) = m = 0 δ m + 1 M m + 1 ( t ) ,  
where M m + 1 ( t )   is the mgf of the random variable Y m + 1 ~ exp-F ( m + 1 ). A third formula for M ( t ) can also be derived based on (19) in terms of the baseline qf. Q F ( u ) as
M ( t ) = m = 0 δ m + 1 ( m + 1 ) 0 1 exp ( t Q F ( u ) ) u m d u ,  
where Q F ( u ) = F 1 ( u ) is the qf of the baseline distribution and u ( 0 , 1 ) .

4.5. Mean Deviations

The amount of scattering in a population is evidently measured to some extent by the totality of deviations from the mean and median. These are known as the mean deviation about the mean and the mean deviation about the median. These measures can be calculated using the following relationships:
δ 1 ( X ) = 2 μ G ( μ ) 2 μ x g ( x ) d x and δ 2 ( X ) = μ 2 M x g ( x ) d x , respectively, where μ = E ( X ) and M = Q ( 1 2 ) .

4.6. Order Statistics

Let X 1 , X 2 , , X n   be a random sample from the GOLE-F family with CDF and PDF defined in Equations (3) and (4), respectively. Suppose X 1 : n , X 2 : n , , X n : n denote the order statistics obtained from this sample and X r : n is the i th order statistic, then the density function of the r th order statistic is provided by
g r : n ( x ) = n ! ( r 1 ) ! ( n r ) ! s = 0 n r ( 1 ) s ( n r s ) g ( x ) G ( x ) r + s 1 .  
From (17), we determine
G ( x ) r + s 1 = [ m = 0 δ m F ( x ) m ] r + s 1 = m = 0 d r + s 1 , m F ( x ) m ,  
where
d r + s 1 , 0 = δ 0 r + s 1 and d r + s 1 , m = ( m δ 0 ) 1 q = 1 m [ q ( r + s ) m ] δ q d r + s 1 , m q .
By replacing t   instead of m in Equation (19), we obtain
g ( x ) = t = 0 δ t + 1 ( t + 1 ) f ( x ) F ( x ) t .  
By substituting (26) in (27) and (28), we determine the PDF of the r th order statistic X r : n as
g r : n ( x ) = t , m = 0 π t , m h t + m + 1 ( x ) ,  
where h t + m + 1 ( x ) denotes the PDF of exp-F distribution with power parameter ( t + m + 1 ) , and
π t , m = s = 0 n r ( 1 ) s ( n r s ) n ! δ t + 1 ( t + 1 ) d r + s 1 , m ( r 1 ) ! ( n r ) ! ( t + m + 1 )   .
Based on Equation (29), several mathematical properties of these order statistics such as ordinary and incomplete moments, factorial moments, moment generating function, mean deviations and several others, can be obtained.

4.7. Stochastic Orderings

Stochastic orders and inequalities are used in many different areas of probability and statistics. Such areas include reliability theory, survival analysis, economics, insurance, actuarial science, queuing theory, biology, operations research, management science, etc. For more detail regarding stochastic ordering, see (Shaked et al., 1994 [15]). Given two random variables X and Y , we say that X is smaller than Y in the:
  • usual stochastic order, denoted by X s t Y , if G X ( x ) G Y ( x ) , for all x ;
  • hazard rate order, denoted by X h r Y , if h X ( x ) h Y ( x ) , for all x ;
  • reversed hazard rate order, denoted by X r h Y , if G X ( x ) / G Y ( x ) , is decreases in x ;
  • mean residual life order, denoted by X m r l Y , if m X ( x ) m Y ( x ) , for all x ;
  • likelihood ratio order, denoted by X l r Y , if g X ( x ) / g Y ( x ) , is decreases in x .
For all the previous orders, we determine the following chains of implications:
X l r Y X h r Y X s t Y
and
X l r Y X r h Y X s t Y ,
also
X h r Y X m r l Y .
For the proposed GOLE-F family, the following theorem provides the stochastic comparison results with respect to the above orderings.
Theorem 1.
Let X ~ G O L E ( a 1 , b 1 ,   c 1 , ϕ ) and Y ~ G O L E ( a 2 , b 2 ,   c 2 , ϕ ) . If a 1 a 2 . and b 1 b 2   and c 1 c 2 , then X s t Y .
Proof. 
If c 1 c 2 , then
F ( x )   c 1 1 F ( x )   c 1 F ( x )   c 2 1 F ( x )   c 2 .  
Hence, if a 1 a 2 and b 1 b 2 , then
a 1 F ( x )   c 1 1 F ( x )   c 1 + b 1 2 ( F ( x )   c 1 1 F ( x )   c 1 ) 2 a 2 F ( x )   c 2 1 F ( x )   c 2 + b 2 2 ( F ( x )   c 2 1 F ( x )   c 2 ) 2 .  
Therefore,
[ ( a 1 F ( x )   c 1 1 F ( x )   c 1 + b 1 2 ( F ( x )   c 1 1 F ( x )   c 1 ) 2 ) ] [ ( a 2 F ( x )   c 2 1 F ( x )   c 2 + b 2 2 ( F ( x )   c 2 1 F ( x )   c 2 ) 2 ) ] .  
Thus,
1 exp [ ( a 1 F ( x )   c 1 1 F ( x )   c 1 + b 1 2 ( F ( x )   c 1 1 F ( x )   c 1 ) 2 ) ]      1 exp [ ( a 2 F ( x )   c 2 1 F ( x )   c 2 + b 2 2 ( F ( x )   c 2 1 F ( x )   c 2 ) 2 ) ] .
That means G X ( x ) G Y ( x ) and X s t Y . □
Theorem 2.
Let X ~ G O L E ( a 1 , b 1 ,   c , ϕ ) and Y ~ G O L E ( a 2 , b 2 ,   c , ϕ ) . If a 1 > a 2 and b 1 = b 2 , then X l r Y .
Proof. 
We determine
g X ( x ) g Y ( x ) = [ ( a 1 + ( b 1 a 1 ) F ( x ) c ) ] × exp [ ( a 1 F ( x ) c 1 F ( x ) c + b 1 2 ( F ( x ) c 1 F ( x ) c ) 2 ) ] [ ( a 2 + ( b 2 a 2 ) F ( x ) c ) ] × exp [ ( a 2 F ( x ) c 1 F ( x ) c + b 2 2 ( F ( x ) c 1 F ( x ) c ) 2 ) ] .
Thus,
log ( g X ( x ) g Y ( x ) ) = log [ ( a 1 + ( b 1 a 1 ) F ( x ) c ) ] ( a 1 F ( x ) c 1 F ( x ) c + b 1 2 ( F ( x ) c 1 F ( x ) c ) 2 ) log [ ( a 2 + ( b 2 a 2 ) F ( x ) c ) ] + ( a 2 F ( x ) c 1 F ( x ) c + b 2 2 ( F ( x ) c 1 F ( x ) c ) 2 ) .
By differentiating the last Equation and after some simplifications, we obtain
d d x ( log ( g X ( x ) g Y ( x ) ) ) = ( a 2 b 1 a 1 b 2 ) c f ( x ) F ( x ) c 1 ( a 1 + ( b 1 a 1 ) F ( x ) c ) ( a 2 + ( b 2 a 2 ) F ( x ) c ) + ( a 2 a 1 ) c f ( x ) F ( x ) c 1 ( 1 F ( x ) c ) 2 + ( b 2 b 1 ) c f ( x ) F ( x ) 2 c 1 ( 1 F ( x ) c ) 3 .
Now, if a 1 > a 2 and b 1 = b 2 , then d d x [ log ( g X ( x ) g Y ( x ) ) ] < 0 , and hence g X ( x ) / g Y ( x ) is decreases in x . This implies that   X l r Y . □

4.8. Stress-Strength Model

The stress–strength model defines the life of an element which has a random strength Y that is subjected to an accidental stress X . The component fails at the instant that the stress applied to it exceeds the strength, and the component will function suitably whenever X < Y . Hence, R = P ( X < Y ) is a measure of component reliability (Kotz et al., 2003 [16]). It has many applications, especially in reliability engineering. We derive the reliability R when Y and X are two independent continuous random variables from the GOLE-F ( a 1 , b 1 , c 1 , ϕ 1 ) and GOLE-F ( a 2 , b 2 , c 2 , ϕ 2 ) distributions, respectively. The reliability is defined by
R = 0 g Y ( x ) G X ( x ) d x .  
Using the PDF in (19) and the CDF in (18), we obtain
R = m , t = 0 δ m + 1 δ t R m + 1 , t   ,  
where
R m + 1 , t = 0 π m + 1 ( x , ϕ 1 ) Π t ( x , ϕ 2 ) d x , and
π m + 1 ( x , ϕ 1 ) = ( m + 1 ) f ( x , ϕ 1 ) F ( x , ϕ 1 ) m , Π t ( x , ϕ 2 ) = F ( x , ϕ 2 ) t .
The constants δ t , δ m + 1 are defined as:
δ t = l = t i = 0 j = 0 i k = 0 ( 1 ) i + l + t + 1 b 2 j a 2 i ( i j ) ( l t ) ( i + j + k 1 k ) ( c 2 ( i + j + k ) l ) i ! 2 j ,
for t 1 . For t = 0 , then
δ 0 = 1 l = 0 i = 0 j = 0 i k = 0 ( 1 ) i + l + m b 2 j a 2 i ( i j ) ( i + j + k 1 k ) ( c 2 ( i + j + k ) l ) i ! 2 j ,
and
δ m + 1 = l = m + 1 i = 0 j = 0 i k = 0 ( 1 ) i + l + m + 2 b 1 j a 1 i ( i j ) ( l m + 1 ) ( i + j + k 1 k ) ( c 1 ( i + j + k ) l ) i ! 2 j ,
for m 0 .
If ϕ 1 = ϕ 2 , then the model reduces to
R = m , t = 0 δ m + 1 δ t ( m + 1 ) m + t + 1 .  

5. Estimation and Simulation

5.1. Estimation of the Parameters

Here, we find the maximum likelihood estimates (MLEs) of the parameters of the new family of distributions from complete samples only. Let x 1 ,   x 2 , ,   x n be observed values from the GOLE F family with parameters a , b , c   and ϕ . Let ξ = ( a , b , c , ϕ ) T be the parameters vector. The total log-likelihood function for ξ is obtained by
l ( ξ ) = n   log c + i = 1 n log f ( x i ; ϕ ) + ( c 1 ) i = 1 n log F ( x i ; ϕ ) + i = 1 n log ( a + ( b a ) F ( x i ; ϕ ) c ) 3 i = 1 n log ( 1 F ( x i ; ϕ ) c ) a i = 1 n H ( x i , ϕ ) b 2 i = 1 n H (   x i , ϕ ) 2 ,
where H (   x i , ϕ ) = F ( x i ; ϕ ) c 1 F ( x i ; ϕ ) c and H (   x i , ϕ ) 2 = ( F ( x i ; ϕ ) c 1 F ( x i ; ϕ ) c ) 2 .
The components of the score vector U ( ξ )   are obtained by
U a = i = 1 n 1 F ( x i ; ϕ ) c a + ( b a ) F ( x i ; ϕ ) c i = 1 n H ( x i , ϕ ) ,
U b = i = 1 n F ( x i ; ϕ ) c a + ( b a ) F ( x i ; ϕ ) c 1 2 i = 1 n H (   x i , ϕ ) 2 ,
U c = n c + i = 1 n log F ( x i ; ϕ ) + i = 1 n ( b a ) F ( x i ; ϕ ) c log F ( x i ; ϕ ) a + ( b a ) F ( x i ; ϕ ) c + 3 i = 1 n F ( x i ; ϕ ) c log F ( x i ; ϕ ) 1 F ( x i ; ϕ ) c a i = 1 n F ( x i ; ϕ ) c log F ( x i ; ϕ ) ( 1 F ( x i ; ϕ ) c ) 2 b i = 1 n F ( x i ; ϕ ) 2 c log F ( x i ; ϕ ) ( 1 F ( x i ; ϕ ) c ) 3 ,
and
U ϕ k = i = 1 n f ( x i ; ϕ ) ϕ k f ( x i ; ϕ ) + ( c 1 ) i = 1 n F ( x i ; ϕ ) ϕ k F ( x i ; ϕ ) + i = 1 n   c ( b a ) F ( x i ; ϕ ) c 1 F ( x i ; ϕ ) ϕ k a + ( b a ) F ( x i ; ϕ ) c + 3 i = 1 n c F ( x i ; ϕ ) c 1 F ( x i ; ϕ ) ϕ k 1 F ( x i ; ϕ ) c a i = 1 n H (   x i , ϕ ) ϕ k b i = 1 n H (   x i , ϕ ) H (   x i , ϕ ) ϕ k
Setting U a ,   U b ,   U c and U ϕ equal to zero, and solving the equations simultaneously, yields the MLE   ξ ^ = ( a ^ , b ^ , c ^ , ϕ ^ ) T of ξ = ( a , b , c , ϕ ) T . These equations cannot be solved analytically, and statistical software can be used to solve them numerically using iterative methods such as the Newton–Raphson type algorithms.

5.2. Simulation Study

In this section, a graphical Monte Carlo simulation study is conducted to compare the performance of the different estimators of the unknown parameters for the GOLE-E ( a , b , c , λ ) distribution. All the computations in this section are conducted using the R program. We generate N = 1000 samples of size n = 20, 25, …, 500 from the GOLE-W and GOLE-E distributions. The true parameter values for GOLE-W ( λ = 1 ) are a = 1.8 , b = 0.5 , c = 1.7 and β = 2.8 , and those for GOLE-E are a = 2 , b = 1.5 , c = 2 and λ = 2.5 , respectively. We also calculate the bias and mean square error (MSE) of the MLEs empirically. The bias and MSE are computed by
B i a s ^ h = 1 N i = 1 N ( h ^ i h ) ,   M S E ^ h = 1 N i = 1 N ( h ^ i h ) 2 .
For h = a , b , c , λ , respectively.
We provide the results of this simulation study in Figure 3, Figure 4, Figure 5 and Figure 6. From these figures, we can perceive that when the sample size increases, the empirical biases and MSEs approach zero in all cases for the two models.

6. Applications on Real-Life Data Sets

In this section, we illustrate the suitability of the proposed family by fitting two real data sets on the special models viz-a-viz GOLE W ( a , b , c , λ , β ) and GOLE E   ( a , b , c , λ ) , arising due to this family with PDF mentioned in Section 3.1 and Section 3.2, respectively. The comparison is conducted with some of the existing models via numerical maximizations of log-likelihood functions using the method of a limited memory quasi-Newton code for bound–constrained maximization (L-BFGS-B). We determine the log-likelihood function adjudicated at the MLEs by estimating the parameters.
Data I: The first data set is related to the measurements of nicotine levels in 346 cigarettes. [https://arxiv.org/ftp/arxiv/papers/1509/1509.08108.pdf, accessed on 19 May 2022]. Data II: The second data set consists of 74 observations of gauge lengths of 20 mm of single carbon fibers pertaining to failure stresses. (Kundu and Raqab, 2009 [17]). The descriptive statistics related to this data sets are given in Table 2.
The total time on test (TTT) plot proposed by Aarset (1987) [18] is a technique to extract information about the shape of the hazard function. This is drawn by plotting T ( i / n ) = { ( r = 1 i y ( r ) ) + ( n i ) y ( i ) } / r = 1 n y ( r ) , where i = 1 , 2 , , n and y ( r ) where r = 1 , 2 , , n is the order statistics of the sample against ( i / n ) . The constant hazard plot is a straight diagonal, while for decreasing (increasing) hazards, it is convex (concave), respectively. The TTT plots for the data sets in Figure 7 indicate that the data sets have an increasing hazard rate.
The best model is chosen on the basis of information criteria such as AIC (Akaike Information Criterion), CAIC (Consistent Akaike Information Criterion), BIC (Bayesian Information Criterion) and HQIC (Hannan–Quinn Information Criterion) with the goodness of fit measures as A* (Anderson–Darling criterion), W* (Cramér–von Mises criterion) and Kolmogorov–Smirnov (K-S) tests with p-values. The model with minimum values for these statistics could be chosen as the best model to fit the data except for the KS p-value, whose maximum value is the desired outcome. Asymptotic standard errors and 95% confidence intervals of the MLEs of the parameters for each competing model are also computed. For visual comparison, the fitted PDFs and the fitted CDFs are plotted with the corresponding observed histograms and ogives.

6.1. Application of GOLE-E

The GOLE-E   ( a , b , c , λ ) distribution is compared with some models, namely exponential (E), moment exponential (ME) (Dara and Ahmad, 2012 [19]), exponentiated moment exponential (EM-E) (Hasnain et al., 2015 [20]), exponentiated exponential (E-E) (Gupta and Kundu, 2001 [21]), beta exponential (B-E) (Nadarajah and Kotz, 2006 [22]) and Kumaraswamy exponential (Kw-E) (Cordeiro and de Castro, 2011 [4]) distributions for all data sets.
In Table 3, Table 4, Table 5 and Table 6, the MLEs, standard errors (SEs) and confidence interval (in parentheses) of the parameters from all the fitted distributions along with the AIC, BIC, CAIC and HQIC for the two data sets are presented. From Table 3, Table 4, Table 5 and Table 6, it is evident that for the data sets, the GOLE-E distribution is the best model with the lowest values of the AIC, BIC, CAIC, HQIC, A*, W* and highest p-value of the K-S statistics. Hence, it is a better model than some recently introduced models, namely exponential (E), moment exponential (ME), exponentiated moment exponential (EM-E), exponentiated exponential (E-E), beta exponential (B-E) and Kumaraswamy exponential (Kw-E) distribution, for the two data sets. More information is provided for a visual comparison in the form of histograms, ogives or cumulative frequency curves of the observed data with the fitted densities and fitted cdfs displayed in Figure 8 and Figure 9. These plots show that the proposed distributions provide the closest fit to all the observed data sets.

6.2. Application of GOLE-W

The GOLE-W ( a , b , c , λ , β ) distribution with ( λ = 1 ) is compared with some models, namely Weibull (W), moment exponential (ME), exponentiated Weibull (EW) (Mudholker and Srivastava, 1993 [23]), generalized Weibull (GW) (Lai 2014 [24]), beta Weibull (B-W) (Lee et al., 2007 [25]) and Kumaraswamy Weibull (Kw-W) (Cordeiro et al. 2010 [26]) distributions for all data sets.
Likewise, in Table 7, Table 8, Table 9 and Table 10, the MLEs, standard errors (in parentheses) and confidence interval [in brackets] of the parameters from all the competitive models along with AIC, CAIC, BIC and HQIC for the two data sets are presented. From these tables, it is quite obvious that for the two data sets, GOLE-W distribution is the best model with the lowest values of AIC, BIC, CAIC, HQIC, A*, W* and highest p-value of the K-S statistics. Hence, it is worth emphasizing that the proposed GOLE-F provides a more useful generalization (with exponential and Weibull as special models) than the competitive models for both of the datasets. A much more useful depiction is presented in the form of a visual comparison in Figure 10 and Figure 11, where the densities and distribution function of observed data are compared against the fitted models, respectively. These plots reveal that the proposed distributions provide the closest fit to all the observed data sets.

7. Conclusions

Through this paper, we provide a new general family of distributions to generalize any continuous baseline distribution. The main properties of the new family and other properties associated with the area of reliability are discussed. It was noted that the distributions generated by the new family are highly flexible in data modeling where we used one member to fit two real data to illustrate the importance of this family. This member provided consistently better fits than the other comparative distributions.

Author Contributions

Data curation, L.H.; Formal analysis, L.H.; Investigation, S.S.; Project administration, A.H.N.A.; Resources, S.K.; Supervision, W.M.; Writing—original draft, F.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are thankful to the Editor-in-Chief and the anonymous referees for their meticulous and thorough reading, which significantly enhanced the readability of this paper, and special thanks to Christophe Chesneau, Department of Mathematics, University of Caen-Normandie, LMNO, France, for their appreciable directions regarding Propositions 1 and 2.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Recalling Equations (3) and (4) and assigned new numbers as (A1) and (A2), respectively, as follows
G ( x ; a , b , c , ϕ ) = 1 exp [ ( a F ( x ; ϕ ) c 1 F ( x ; ϕ ) c + b 2 ( F ( x ; ϕ ) c 1 F ( x ; ϕ ) c ) 2 ) ] .  
g ( x ; a , b , c , ϕ ) = [ c f ( x ; ϕ ) F ( x ; ϕ ) c 1 ( a + ( b a ) F ( x ; ϕ ) c ) ( 1 F ( x ; ϕ ) c ) 3 ] × exp [ ( a F ( x ; ϕ ) c 1 F ( x ; ϕ ) c + b 2 ( F ( x ; ϕ ) c 1 F ( x ; ϕ ) c ) 2 ) ] .
Proposition A1.
Given x = F ( x ) , by using the equivalence: e y ~ 1 + y when y 0 since lim x F ( x ) c 0 . Then, by the properties of the CDF in Equation (A1), we arrive at
G ( x )   ~   a F ( x ) c 1 F ( x ) c + b 2 ( F ( x ) c 1 F ( x ) c ) 2 ,
and, by asymptotic dominance, we obtain
G ( x ) ~ a   F ( x ) c .
Using the same arguments, we obtain
g ( x ) ~ c   a   f ( x ) F ( x ) c 1 .
In addition, the survival function is close to one; thus, the denominator in the hazard function is close to one. Then, using Equations (A3) and (A4), we obtain
h ( x )   ~   c   a   f ( x ) F ( x ) c 1 .
Proposition A2.
Similarly, using the same arguments when lim x + F ( x ) c 1 , we can prove that the survival function can be approximately reduced as follows
1 G ( x ) ~ 1 e ( a 1 F ( x ) c   + b 2 { 1 1 F ( x ) c   } 2 ) .
Using the same arguments, we obtain
g ( x ) ~ b c f ( x ) ( 1 F ( x ) c ) 3     e ( a 1 F ( x ) c   + b 2 { 1 1 F ( x ) c   } 2 ) .
Using Equations (A6) and (A7), we obtain
h ( x ) ~ b c f ( x ) ( 1 F ( x ) c ) 3 .
This completes the proof.

References

  1. Marshall, A.W. 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]
  2. Gupta, R.C.; Gupta, P.L.; Gupta, R.D. Modeling failure time data by Lehman alternatives. Commun. Stat.-Theory Methods 1998, 27, 887–904. [Google Scholar] [CrossRef]
  3. Eugene, N.; Lee, C.; Famoye, F. Beta-normal distribution and its applications. Commun. Stat. Theory Methods 2002, 31, 497–512. [Google Scholar] [CrossRef]
  4. Cordeiro, G.M.; De Castro, M. A new family of generalized distributions. J. Stat. Comput. Simul. 2011, 81, 883–898. [Google Scholar] [CrossRef]
  5. Alexander, C.; Cordeiro, G.M.; Ortega, E.M.; Sarabia, J.M. Generalized beta-generated distributions. Comput. Stat. Data Anal. 2012, 56, 1880–1897. Available online: https://EconPapers.repec.org/RePEc:eee:csdana:v:56:y:2012:i:6 (accessed on 19 May 2022). [CrossRef]
  6. 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]
  7. Bourguignon, M.; Silva, R.B.; Cordeiro, G.M. The Weibull-G Family of Probability Distributions. J. Data Sci. 2014, 12, 53–68. Available online: http://www.jds-online.com/volume-12-number-1-january-2014 (accessed on 19 May 2022). [CrossRef]
  8. Tahir, M.H.; Cordeiro, G.M.; Alizadeh, M.; Mansoor, M.; Zubair, M.; Hamedani, G.G. The odd generalized exponential family of distributions with applications. J. Stat. Distrib. Appl. 2015, 2, 1. [Google Scholar] [CrossRef] [Green Version]
  9. Cordeiro, G.M.; Alizadeh, M.; Ortega, E.M.; Serrano, L.H.V. The Zografos-Balakrishnan odd log-logistic family of distributions: Properties and Applications. Hacet. J. Math. Stat. 2015, 46, 11781–11803. Available online: https://dergipark.org.tr/hujms/issue/43489/524407 (accessed on 19 May 2022). [CrossRef]
  10. Gomes-Silva, F.; Percontini, A.; De Brito, E.; Ramos, M.W.; Venâncio, R.; Cordeiro, G.M. The Odd Lindley-G Family of Distributions. Austrian J. Stat. 2017, 46, 65–87. [Google Scholar] [CrossRef] [Green Version]
  11. Alizadeh, M.; Cordeiro, G.M.; Pinho, L.G.B.; Ghosh, I. The Gompertz-G family of distributions. J. Stat. Theory Pract. 2016, 11, 179–207. [Google Scholar] [CrossRef]
  12. Jamal, F.; Nasir, M.A.; Tahir, M.H.; Montazeri, N.H. The odd Burr-III family of distributions. J. Stat. Appl. Probab. 2017, 6, 105–122. Available online: http://www.naturalspublishing.com/files/published/4nk6g57u5e512l.pdf (accessed on 19 May 2022). [CrossRef]
  13. Khan, S.; Balogun, O.S.; Tahir, M.H.; Almutiry, W.; Alahmadi, A.A. An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data. Symmetry 2021, 13, 2064. [Google Scholar] [CrossRef]
  14. Sen, A.; Bhattacharyya, G.K. Inference procedures for the linear failure rate model. J. Stat. Plan. Inference 1995, 46, 59–76. [Google Scholar] [CrossRef]
  15. Shaked, M.; Shanthikumar, J.G. Stochastic Orders and Their Applications; Academic Press: San Diego, CA, USA, 2014. [Google Scholar] [CrossRef]
  16. Kotz, S.; Pensky, M. The Stress-Strength Model and Its Generalizations: Theory and Applications; World Scientific: Singapore, 2003. [Google Scholar] [CrossRef]
  17. Kundu, D.; Raqab, M.Z. Estimation of R = P (Y < X) for three-parameter Weibull distribution. Stat. Probab. Lett. 2009, 79, 1839–1846. [Google Scholar] [CrossRef]
  18. Aarset, M.V. How to Identify a Bathtub Hazard Rate. IEEE Trans. Reliab. 1987, 36, 106–108. [Google Scholar] [CrossRef]
  19. Dara, S.T.; Ahmad, M. Recent Advances in Moment Distribution and Their Hazard Rates; Lap Lambert Academic Publishing: Chisinau, Republic of Moldova, 2012. [Google Scholar]
  20. Hasnain, S.A.; Iqbal, Z.; Ahmad, M. On exponentiated moment exponential distribution. Pak. J. Stat. 2015, 31, 267–280. Available online: https://www.statindex.org/journals/1313/31/2 (accessed on 19 May 2022).
  21. Gupta, R.D.; Kundu, D. Exponentiated Exponential Family: An Alternative to Gamma and Weibull Distributions. Biom. J. 2001, 43, 117–130. [Google Scholar] [CrossRef]
  22. Nadarajah, S.; Kotz, S. The beta exponential distribution. Reliab. Eng. Syst. Saf. 2006, 91, 689–697. [Google Scholar] [CrossRef]
  23. Mudholkar, G.; Srivastava, D. Exponentiated Weibull family for analyzing bathtub failure-rate data. IEEE Trans. Reliab. 1993, 42, 299–302. [Google Scholar] [CrossRef]
  24. Lai, C.D. Generalized Weibull Distributions. In Generalized Weibull Distributions; Springer Briefs in Statistics; Springer: Berlin/Heidelber, Germany, 2014. [Google Scholar] [CrossRef]
  25. Lee, C.; Famoye, F.; Olumolade, O. Beta-Weibull Distribution: Some Properties and Applications to Censored Data. J. Mod. Appl. Stat. Methods 2007, 6, 173–186. [Google Scholar] [CrossRef]
  26. Cordeiro, G.M.; Ortega, E.M.; Nadarajah, S. The Kumaraswamy Weibull distribution with application to failure data. J. Frankl. Inst. 2010, 347, 1399–1429. [Google Scholar] [CrossRef]
Figure 1. (a) Density function and (b) hazard rate plots of the GOLE-W distribution for different parameter values.
Figure 1. (a) Density function and (b) hazard rate plots of the GOLE-W distribution for different parameter values.
Mca 27 00055 g001
Figure 2. Plots of (a) density function and (b) hazard rate of the GOLE-E distribution for different parameter values.
Figure 2. Plots of (a) density function and (b) hazard rate of the GOLE-E distribution for different parameter values.
Mca 27 00055 g002
Figure 3. The biases of the estimates of parameters of the GOLE-W distribution.
Figure 3. The biases of the estimates of parameters of the GOLE-W distribution.
Mca 27 00055 g003
Figure 4. The MSEs of the estimates of parameters of the GOLE-W distribution.
Figure 4. The MSEs of the estimates of parameters of the GOLE-W distribution.
Mca 27 00055 g004
Figure 5. The biases of the estimates of parameters of the GOLE-E distribution.
Figure 5. The biases of the estimates of parameters of the GOLE-E distribution.
Mca 27 00055 g005
Figure 6. The MSE of the estimates of parameters of the GOLE-E distribution.
Figure 6. The MSE of the estimates of parameters of the GOLE-E distribution.
Mca 27 00055 g006
Figure 7. TTT plots of the data set I and II.
Figure 7. TTT plots of the data set I and II.
Mca 27 00055 g007
Figure 8. Plots of (a) the fitted PDF and (b) estimated CDF for the GOLE-E distribution for data set I.
Figure 8. Plots of (a) the fitted PDF and (b) estimated CDF for the GOLE-E distribution for data set I.
Mca 27 00055 g008
Figure 9. Plots of (a) the fitted PDF and (b) estimated CDF for the GOLE-E distribution for data set II.
Figure 9. Plots of (a) the fitted PDF and (b) estimated CDF for the GOLE-E distribution for data set II.
Mca 27 00055 g009
Figure 10. Plots of (a) the fitted PDF for the GOLE-W distribution and (b) estimated CDF for the GOLE-W distribution for data set I.
Figure 10. Plots of (a) the fitted PDF for the GOLE-W distribution and (b) estimated CDF for the GOLE-W distribution for data set I.
Mca 27 00055 g010
Figure 11. Plots of (a) the fitted PDF for the GOLE-W distribution and (b) estimated CDF for the GOLE-W distribution for data set II.
Figure 11. Plots of (a) the fitted PDF for the GOLE-W distribution and (b) estimated CDF for the GOLE-W distribution for data set II.
Mca 27 00055 g011
Table 1. Mean, variance, skewness and kurtosis of the GOLE-E distribution with different values of a , b , c and λ = 1 .
Table 1. Mean, variance, skewness and kurtosis of the GOLE-E distribution with different values of a , b , c and λ = 1 .
abcMeanVarianceSkewnessKurtosis
0.50.510.67040.13301.29791.8368
21.16670.21861.17761.4755
51.95170.30041.09601.2512
102.59870.33511.06351.1656
203.26830.35411.04401.1146
504.17040.36601.02881.0753
1004.85870.37011.02181.0571
1110.46140.08241.38692.1392
20.89950.15871.21801.5987
51.64110.24081.11031.2923
102.27210.27761.06991.1839
202.93340.29821.04671.1226
503.83030.31141.02941.0774
1004.51690.31591.02181.0574
21.510.30910.04881.51192.6112
20.68600.11301.27031.7688
51.37900.19191.12691.3424
101.99140.22971.07681.2044
202.64300.25141.04941.1308
503.53390.26541.02991.0793
1004.21850.27021.02171.0574
Table 2. Descriptive Statistics for the data set I and data set II.
Table 2. Descriptive Statistics for the data set I and data set II.
Data SetsMin.MeanMedianS.D.SkewnessKurtosis1st Q.3rd Q.Max.
I0.100.850.900.330.170.290.601.102.00
II1.3122.4772.5130.487−0.151−0.1272.1502.8163.5
Table 3. MLEs, standard error (in parentheses), confidence interval values [in brackets] for the data set I.
Table 3. MLEs, standard error (in parentheses), confidence interval values [in brackets] for the data set I.
Models a ^ b ^ c ^ λ ^
GOLE-E ( a , b , c , λ )
0.2188.9491.3450.554
(0.315)(3.246)(0.237)(0.083)
[0, 0.84][2.58, 15.31][0.88, 1.81][0.39, 0.72]
Kw-E ( a , b , λ )
3.020105.575-0.252
(0.163)(38.348)(0.045)
[2.70, 3.34][30.41, 180.73][0.160.34]
B-E ( a , b , λ )
4.92217.433-0.298
(0.364)(8.216)(0.128)
[4.21, 5.64][1.32, 33.54][0.05, 0.55]
E-E ( b , λ )
-5.526-2.726
(0.514)(0.128)
[4.52, 6.53][2.475, 2.98]
EM-E ( a , b )
2.5740.284--
(0.229)(0.012)
[2.13, 3.02][0.26, 0.31]
M-E ( b )
-0.406--
(0.016)
[0.37, 0.44]
E ( λ )
---1.173
(0.063)
[1.04, 1.29]
Table 4. The AIC, BIC, CAIC, HQIC, A*, W* and KS (p-value) values for data set I.
Table 4. The AIC, BIC, CAIC, HQIC, A*, W* and KS (p-value) values for data set I.
ModelsAICBICCAICHQICA*W*KS
(p-Value)
GOLE-E ( a , b , c , λ )
232.14247.54232.28238.302.670.470.25
(0.29)
Kw-E ( a , b , λ )
236.92248.46236.99241.513.370.580.12
(0.03)
B-E ( a , b , λ )
276.04287.59276.11280.666.481.090.24
(0.16)
E-E ( b , λ )
302.44310.14302.47305.529.421.590.22
(0.23)
EM-E ( a , b )
290.62298.32290.65293.708.481.460.24
(0.20)
M-E ( b )
388.70392.55388.71390.246.491.090.23
(0.008)
E ( λ )
583.66587.51583.67585.206.541.110.34
(0.002)
Table 5. MLEs, standard error (in parentheses) and confidence interval values [in brackets] for data set II.
Table 5. MLEs, standard error (in parentheses) and confidence interval values [in brackets] for data set II.
Models a ^ b ^ c ^
λ ^
GOLE-E ( a , b , c , λ )
0.3651.2994.0912.748
(0.160)(0.657)(1.248)(0.531)
[0.05, 0.68][0.01, 2.59][1.64, 6.54][1.71, 3.78]
Kw-E ( a , b , λ )
12.47324.773-0.559
(3.939)(23.936)(0.194)
[4.75, 20.19][0, 71.68][0.17, 0.93]
B-E ( a , b , λ )
26.25914.354-0.421
(5.838)(17.832)(0.376)
[14.81, 37.70][0, 49.30][0, 1.16]
E-E ( b , λ )
-89.394-2.018
(32.458)(0.171)
[25.77, 153.01][1.68, 2.35]
EM-E ( a , b )
32.3190.418--
(10.705)(0.032)
[11.33, 53.30][0.35, 0.48]
M-E ( b )
-1.238--
(0.101)
[1.04, 1.44]
E ( λ )
---0.403
(0.046)
[0.31, 0.49]
Table 6. The AIC, BIC, CAIC, HQIC, A*, W* and KS (p-value) values for data set II.
Table 6. The AIC, BIC, CAIC, HQIC, A*, W* and KS (p-value) values for data set II.
ModelsAICBICCAICHQICA*W*KS
(p-Value)
GOLE-E ( a , b , c , λ )
107.90116.20108.48111.580.430.040.06
(0.83)
Kw-E ( a , b , λ )
112.66119.56113.00115.360.520.060.07
(0.79)
B-E ( a , b , λ )
116.82123.72117.16119.520.620.090.08
(0.71)
E-E ( b , λ )
121.60126.20121.76123.401.040.160.01
(0.44)
EM-E ( a , b )
119.90124.50120.07121.700.630.100.09
(0.52)
M-E ( b )
230.16232.46230.22231.060.580.080.35
(0.002)
E ( λ )
284.24286.54284.29285.140.570.090.44
(0.01)
Table 7. MLEs, standard errors (in parentheses) and confidence interval [in brackets] values for data set I.
Table 7. MLEs, standard errors (in parentheses) and confidence interval [in brackets] values for data set I.
Models a ^ b ^ c ^
λ ^
β ^
GOLE-W ( a , b , c , β )
2.3893114.66534.8673-0.506777
(2.1340)(50.1098)(2.1515)(0.1582)
[0, 6.57][16.45, 212.88][0.65, 9.08][0.20, 0.82]
Kw-W ( a , b , λ , β )
0.71030.2623-3.04643.8368
(0.0233) (0.0263)(0.0174)
[0.66, 0.76][0.23,0.29][2.99, 3.10][3.80, 3.87]
B-W ( a , b , λ , β )
0.77300.2276-3.02014.3742
(0.0673)(0.0137)(0.0042)(0.0042)
[0.64, 0.90][0.20, 0.25][3.01, 3.02][4.37, 4.38]
E-W ( b , λ , β )
-0.8090-3.0689220.9440
(0.1515)(0.3541)(0.1732)
(4.52, 6.53)(2.475, 2.98)[0.60, 1.28]
G-W ( a , λ , β )
0.5597--2.71902.0240
(11.2701)(0.1140)(0.7523)
[0, 22.65][2.50, 2.94][0.55, 3.50]
M-E ( b )
-0.406---
(0.016)
[0.37, 0.44]
W ( λ , β )
---1.1329262.71898
(0.0623)(0.1140)
[1.01, 1.26][2.50, 2.94]
Table 8. The AIC, CAIC, BIC, HQIC, A*, W* and KS (p-value) values for data set I.
Table 8. The AIC, CAIC, BIC, HQIC, A*, W* and KS (p-value) values for data set I.
ModelsAICCAICBICHQICA*W*KS
(p-Value)
GOLE-W ( a , b , c , β )
230.11230.23245.49236.232.510.440.10
(0.17)
Kw-W ( a , b , λ , β )
231.86231.96247.23237.972.570.450.11
(0.03)
B-W ( a , b , λ , β )
232.84232.95248.22238.962.670.460.12
(0.01)
E-W ( b , λ , β )
232.42232.39243.86236.912.810.480.12
(0.000)
G-W ( a , λ , β )
233.56233.63245.09238.152.970.510.24
(0.000)
M-E ( b )
388.70392.55388.71390.246.491.090.23
(0.008)
W ( λ , β )
231.56231.88239.25234.622.970.510.14
(0.000)
Table 9. MLEs, standard errors (in parentheses), confidence interval values [in brackets] for data set II.
Table 9. MLEs, standard errors (in parentheses), confidence interval values [in brackets] for data set II.
Models a ^ b ^
c ^
λ ^
β ^
GOLE-W ( a , b , c , β )
6.9553114.665320.1203-0.8448
(6.6794)(50.1098)(6.6638)(0.31433)
[0,20.04][0, 27.65][7.06, 33.18][0.23, 1.46]
Kw-W ( a , b , λ , β )
1.66461.0950-4.36750.0187
(0.8438)(0.5829)(2.8871)(0.0192)
[0.01, 3.32][0.23, 0.29][0, 10.0262][0, 0.0564]
B-W ( a , b , λ , β )
1.74010.9961-4.29870.0222
(1.3064)(0.9249)(2.0778)(0.0253)
[0, 4.30][0, 2.81][0.23, 8.37][0, 0.07]
E-W ( b , λ , β )
-1.7298-4.30830.9440
(0.7208)(0.9066)(0.1732)
[0.32, 3.14][2.53, 6.09][0, 0.07]
G-W ( a , λ , β )
1.5460--5.38160.0033
(0.9021)(0.4906)(0.0008)
[0, 3.3142][4.42, 6.34][0.002, 0.005]
M-E ( b )
-0.406---
(0.016)
[0.37, 0.44]
W ( λ , β )
---0.00365.7342
(0.0009)(0.2428)
[0.0002, 0.0053][5.26, 6.21]
Table 10. The AIC, CAIC, BIC, HQIC, A*, W* and KS (p-value) values for data set II.
Table 10. The AIC, CAIC, BIC, HQIC, A*, W* and KS (p-value) values for data set II.
ModelsAICCAICBICHQICA*W*KS
(p-Value)
GOLE-W ( a , b , c , β )
110.57111.15119.79114.250.250.0310.06
(0.93)
Kw-W ( a , b , λ , β )
111.06112.83120.08114.982.270.0370.08
(0.91)
B-W ( a , b , λ , β )
111.32112.90120.13115.990.260.0380.07
(0.92)
E-W ( b , λ , β )
118.33118.72122.89120.680.310.0750.098
(0.89)
G-W ( a , λ , β )
113.84113.13123.02119.150.300.0520.08
(0.88)
M-E ( b )
388.70392.55388.71390.246.491.090.23
(0.008)
W ( λ , β )
117.45116.61121.77117.910.290.0370.09
(0.87)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jamal, F.; Handique, L.; Ahmed, A.H.N.; Khan, S.; Shafiq, S.; Marzouk, W. The Generalized Odd Linear Exponential Family of Distributions with Applications to Reliability Theory. Math. Comput. Appl. 2022, 27, 55. https://doi.org/10.3390/mca27040055

AMA Style

Jamal F, Handique L, Ahmed AHN, Khan S, Shafiq S, Marzouk W. The Generalized Odd Linear Exponential Family of Distributions with Applications to Reliability Theory. Mathematical and Computational Applications. 2022; 27(4):55. https://doi.org/10.3390/mca27040055

Chicago/Turabian Style

Jamal, Farrukh, Laba Handique, Abdul Hadi N. Ahmed, Sadaf Khan, Shakaiba Shafiq, and Waleed Marzouk. 2022. "The Generalized Odd Linear Exponential Family of Distributions with Applications to Reliability Theory" Mathematical and Computational Applications 27, no. 4: 55. https://doi.org/10.3390/mca27040055

APA Style

Jamal, F., Handique, L., Ahmed, A. H. N., Khan, S., Shafiq, S., & Marzouk, W. (2022). The Generalized Odd Linear Exponential Family of Distributions with Applications to Reliability Theory. Mathematical and Computational Applications, 27(4), 55. https://doi.org/10.3390/mca27040055

Article Metrics

Back to TopTop