Next Article in Journal
Relative Tail Pressure with Subadditive Potentials
Previous Article in Journal
Fractal Plasmons on Cantor Set Thin Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Power Topp–Leone Generated Family of Distributions with Applications

by
Rashad A. R. Bantan
1,2,
Farrukh Jamal
3,
Christophe Chesneau
4,* and
Mohammed Elgarhy
5
1
Department of Marine Geology, Faculty of Marine Sience, King AbdulAziz University, Jeddah 21551, Saudi Arabia
2
Deanship of Scientific Research, King AbdulAziz University, Jeddah 21551, Saudi Arabia
3
Department of Statistics, Govt. S.A Postgraduate College Dera Nawab Sahib, Bahawalpur, Punjab 63100, Pakistan
4
Department of Mathematics, Université de Caen, LMNO, Campus II, Science 3, 14032 Caen, France
5
Valley High Institute for Management Finance and Information Systems, Obour, Qaliubia 11828, Egypt
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(12), 1177; https://doi.org/10.3390/e21121177
Submission received: 12 November 2019 / Revised: 24 November 2019 / Accepted: 28 November 2019 / Published: 29 November 2019

Abstract

:
In this paper, we introduce a new general family of distributions obtained by a subtle combination of two well-established families of distributions: the so-called power Topp–Leone-G and inverse exponential-G families. Its definition is centered around an original cumulative distribution function involving exponential and polynomial functions. Some desirable theoretical properties of the new family are discussed in full generality, with comprehensive results on stochastic ordering, quantile function and related measures, general moments and related measures, and the Shannon entropy. Then, a statistical parametric model is constructed from a special member of the family, defined with the use of the inverse Lomax distribution as the baseline distribution. The maximum likelihood method was applied to estimate the unknown model parameters. From the general theory of this method, the asymptotic confidence intervals of these parameters were deduced. A simulation study was conducted to evaluate the numerical behavior of the estimates we obtained. Finally, in order to highlight the practical perspectives of the new family, two real-life data sets were analyzed. All the measures considered are favorable to the new model in comparison to four serious competitors.

1. Introduction

Owing to the growing amount of data from various applied fields and unstoppable computer progress, there is increasing motivation on developing efficient and flexible statistical models. Such models can be derived from general families of distributions having desirable properties, such as those constructed from a generator distribution. The main idea of this construction is to add shape parameter(s) to a baseline distribution with the aim to upgrade its flexibility level. Among the well-known examples of such families, there are the beta-G [1], Kumaraswamy-G [2], Weibull-GG [3], Garhy-G [4], type II half logistic-G [5], Transmuted Topp–Leone G [6], generalized odd log-logistic-G [7], odd Fréchet-G [8], power Lindley-G [9], Fréchet Topp–Leone-G [10], exponentiated generalized Topp–Leone-G [11], and truncated inverted Kumaraswamy-G [12]. We also refer to the exhaustive survey in [13]. Recently, several researchers used the Topp–Leone (TL) distribution as generator distribution to develop new general families, reaching the aims of simplicity and flexibility. Among them, Ref. [14] proposed the Topp–Leone-G (TL-G) family, Ref. [15] introduced the power TL-G (PTL-G) family, Ref. [16] introduced the generalized TL-G family, Ref. [17] studied the type II TL-G family, and [18] proposed the type II generalized TL-G family.
For the purposes of this paper, let us describe in detail the PTL-G family from [15]. The PTL-G family is defined by the following cumulative distribution function (cdf):
F ( x ; α , β , ξ ) = G ( x ; ξ ) α β [ 2 G ( x ; ξ ) β ] α , x R ,
with α , β > 0 , where G ( x ; ξ ) is a cdf of a baseline continuous distribution which may depend on a vector parameter ξ ; i.e., ξ = ( ξ 1 , ξ 2 , ) . As indicated by the name, the construction of the family uses the so-called power Topp–Leone distribution as the generator distribution. In comparison to the (power one) TL-G family, Ref. [15] demonstrated the significant impact of the parameter β on the shapes of the probability density and hazard rate functions, providing desirable modeling properties. This is particularly flagrant with the consideration of the gamma distribution as the baseline distribution, as illustrated by the graphics and applications of [15].
In a parallel work, beyond the TL distribution and its extensions, Ref. [19] introduced the inverse exponential-G (IE-G) family, based on the inverse exponential distribution as the generator distribution, and defined by the following cdf:
F ( x ; ξ ) = e 1 1 G ( x ; ξ ) , x R .
The main features of this family are being simple, with no new, additional parameters, and having a completely different nature of the former baseline cdf G ( x ; ξ ) owing to the combination the exponential (implicit) odd functions. An immediate remark illustrating this claim is the following: it has a fastest rate of decay to 0 when G ( x ; ξ ) 0 . By the consideration of a practical data set and the exponential distribution as baseline distribution, Ref. [19] shows that the corresponding model is better than the Lindley and exponential models (all having the same number of parameters). The nice results behind the IE-G family have been the driver for more investigations, with extended or modified versions of this family. We refer the reader to [8,20] for the odd Fréchet-G family, Ref. [21] for the extended odd Fréchet-G family, and [22] for the modified odd Fréchet-G family.
In this paper, in view of the previously mentioned literature, we introduce a new family of distributions by combining, in some senses, the PTL-G and IE-G families. It is defined by composition of their respective cdfs, i.e., by the cdf given by
F ( x ; α , β , ξ ) = e α β 1 1 G ( x ; ξ ) 2 e β 1 1 G ( x ; ξ ) α , x R ,
with α , β > 0 . Thus, this cdf can be view as a polyno-exponential transformation of the baseline cdf G ( x ; ξ ) . The new family is called the new power TL-G (NPTL-G) family. Thus, by construction, we aim to combine the benefits of the PTL-G and IE-G families, and thus, create new statistical perspectives of various kinds. The key motivations behind the NPTL-G family are the following.
  • To provide very simple models and create new simple distributions.
  • To improve the flexibility of existing distributions on various aspects (such as mode, median, skewness, and kurtosis…).
  • To provide better fits than competing modified models having the same of higher number of parameters.
We support these claims both in full generality and by putting the light on the special member of the NPTL-G family defined with the inverse Lomax (ILx) distribution as the baseline distribution (the reason of this choice will be explained later). The resulting distribution, called the new power Topp–Leone inverse Lomax (NPTLILx) distribution, offers a new three-parameter lifetime distribution, with a high potential of applicability. We illustrate that by the means of two practical data sets with different features: the first one is from [23] and is about active repair times for airborne communication transceiver, and the second one is from [24] and is about actual tax revenue in Egypt. Favorable results were obtained for the proposed model in comparison to serious competitors, motivating its use wider statistical uses.
The contents of this paper are organized as follows. In Section 2, the basics of the NPTL-G family are presented, as is the NPTLILx distribution. Various mathematical properties of the family are discussed in Section 3. Section 4 is devoted to the estimation of the unknown parameters from the NPTLILx model, with a comprehensive simulation study. The data analyses are shown in Section 5 with numerical and graphical illustrations. A conclusion and perspectives are formulated in Section 6.

2. Basics of the NPTL-G Family

The basics of the NPTL-G family are presented in this section, with a focus on the main functions of interest.

2.1. Probability Density Function

Upon differentiation of F ( x ; α , β , ξ ) according to x, owing to (1), the probability density function (pdf) of the NPTL-G family is given by
f ( x ; α , β , ξ ) = 2 α β g ( x ; ξ ) G ( x ; ξ ) 2 e α β 1 1 G ( x ; ξ ) 1 e β 1 1 G ( x ; ξ ) 2 e β 1 1 G ( x ; ξ ) α 1 , x R ,
where g ( x ; ξ ) is the probability density function corresponding to G ( x ; ξ ) . From this expression, some asymptotic results on f ( x ; α , β , ξ ) can be derived. When G ( x ; ξ ) 0 , we have
f ( x ; α , β , ξ ) 2 α β g ( x ; ξ ) G ( x ; ξ ) 2 e α β 1 1 G ( x ; ξ ) .
Furthermore, when G ( x ; ξ ) 1 , we have
f ( x ; α , β , ξ ) 2 α β 2 g ( x ; ξ ) 1 G ( x ; ξ ) .
The variations of f ( x ; α , β , ξ ) can be studied in a standard manner, starting with the critical point(s) given by the solution of the non-linear equation according to x: ln [ f ( x ; α , β , ξ ) ] = 0 , with
ln [ f ( x ; α , β , ξ ) ] = g ( x ; ξ ) g ( x ; ξ ) 2 g ( x ; ξ ) G ( x ; ξ ) + α β g ( x ; ξ ) G ( x ; ξ ) 2 β g ( x ; ξ ) e β 1 1 G ( x ; ξ ) G ( x ; ξ ) 2 1 e β 1 1 G ( x ; ξ ) β ( α 1 ) g ( x ; ξ ) e β 1 1 G ( x ; ξ ) G ( x ; ξ ) 2 2 e β 1 1 G ( x ; ξ ) .
Then, for a critical point x c , the sign of ln [ f ( x ; α , β , ξ ) ] x = x c is informative on its nature (minimum, maximum, or inflection point).

2.2. Hazard Rate Function

The hazard rate function (hrf) of the NPTL-G family is given by
h ( x ; α , β , ξ ) = f ( x ; α , β , ξ ) 1 F ( x ; α , β , ξ ) = 2 α β g ( x ; ξ ) e α β 1 1 G ( x ; ξ ) 1 e β 1 1 G ( x ; ξ ) 2 e β 1 1 G ( x ; ξ ) α 1 G ( x ; ξ ) 2 1 e α β 1 1 G ( x ; ξ ) 2 e β 1 1 G ( x ; ξ ) α , x R .
Some asymptotic results on h ( x ; α , β , ξ ) are presented below. When G ( x ; ξ ) 0 , we have
h ( x ; α , β , ξ ) 2 α β g ( x ; ξ ) G ( x ; ξ ) 2 e α β 1 1 G ( x ; ξ ) .
Additionally, when G ( x ; ξ ) 1 , we have
h ( x ; α , β , ξ ) 2 g ( x ; ξ ) 1 G ( x ; ξ ) 1 .
Thus, the parameters α and β have a significant effect on the asymptotes when G ( x ; ξ ) 0 , but no effect when G ( x ; ξ ) 1 . The variations of h ( x ; α , β , ξ ) can be studied in similar manner to f ( x ; α , β , ξ ) by using the relation ln [ h ( x ; α , β , ξ ) ] = ln [ f ( x ; α , β , ξ ) ] + h ( x ; α , β , ξ ) .

2.3. A Special Member: The NPTLILx Distribution

The NPTL-G family contains distributions of various natures, depending on the choice of the baseline distribution. In this study, as evoked in the introduction, we chose the inverse Lomax distribution with shape parameter θ > 0 as the baseline distribution to define the NPTLILx distribution. Thus, it is defined by the following cdf:
G ( x ; θ ) = 1 + x 1 θ , x > 0 ,
(another parameter of the former definition of the inverse Lomax distribution has been reduced to 1 for the purposes of the paper). Let us now briefly motivate this choice. As suggested by its name, the inverse Lomax distribution is the distribution of the random variable Y = 1 / X , where X denotes a random variable following the standard Lomax distribution (with parameters θ and 1). The corresponding pdf and hrf are, respectively, given by
g ( x ; θ ) = θ x 2 1 + x 1 θ 1
and
h ( x ; θ ) = θ x 2 1 + x 1 θ 1 1 1 + x 1 θ , x > 0 .
In addition to being simple, it has been proven to be a very flexible to model data having a subjacent non-monotonic hrf. Further details and applications can be found in [25,26,27].
Thus, the NPTLILx distribution is defined by the following cdf:
F ( x ; α , β , θ ) = e α β 1 1 + x 1 θ 2 e β 1 1 + x 1 θ α , x > 0 ,
with α , β , θ > 0 . The corresponding pdf and hrf are given by, respectively,
f ( x ; α , β , θ ) = 2 α β θ x 2 1 + x 1 θ 1 e α β 1 1 + x 1 θ 1 e β 1 1 + x 1 θ 2 e β 1 1 + x 1 θ α 1
and
h ( x ; α , β , θ ) = 2 α β θ x 2 1 + x 1 θ 1 e α β 1 1 + x 1 θ 1 e β 1 1 + x 1 θ 2 e β 1 1 + x 1 θ α 1 1 e α β 1 1 + x 1 θ 2 e β 1 1 + x 1 θ α , x > 0 .
Possible shapes of the pdf and hrf of the NPTLILx distribution are illustrated in Figure 1 and Figure 2, respectively. In particular, from Figure 1, we see that the pdf can be right skewed and reversed-J shaped. From Figure 2, we see that the hrf can be increasing, decreasing, upside down, and bathtub shaped. All these curvature properties are known to be desirable to create flexible statistical models.

3. Some Mathematical Properties

The section presents some important mathematical properties of the NPTL-G family.

3.1. On a Stochastic Ordering

The following result shows some inequalities involving F ( x ; α , β , ξ ) .
Proposition 1.
For any x R such that G ( x ; ξ ) > 0 , the following inequalities hold:
e α β 1 1 G ( x ; ξ ) 2 G ( x ; ξ ) β α F ( x ; α , β , ξ ) 2 α e α β 1 1 G ( x ; ξ ) .
Proof. 
The bracket term in the definition of F ( x ; α , β , ξ ) given by (1) is central. Since e β 1 1 G ( x ; ξ ) ( 0 , 1 ) , we have 2 e β 1 1 G ( x ; ξ ) 2 , implying the second inequality. For the first inequality, the following well-known logarithmic inequality: for y > 1 , y 1 + y ln ( 1 + y ) gives e 1 1 1 + y 1 + y , implying that e 1 1 G ( x ; ξ ) G ( x ; ξ ) by taking y = G ( x ; ξ ) 1 . Therefore, we have e β 1 1 G ( x ; ξ ) G ( x ; ξ ) β , and a fortiori, 2 G ( x ; ξ ) β 2 e β 1 1 G ( x ; ξ ) . The first inequality follows. This ends the proof of Proposition 1. □
An immediate consequence of Proposition 1 is the following stochastic ordering result:
F * ( x ; α , β , ξ ) F ( x ; α , β , ξ ) ,
where F * ( x ; α , β , ξ ) = e α β 1 1 G ( x ; ξ ) is the cdf of the exponentiated IE-G family (with power parameter  α β ).
Another stochastic ordering result comes from the following remark: the function F o ( x ; α , β , ξ ) given by
F o ( x ; α , β , ξ ) = e α β 1 1 G ( x ; ξ ) 2 G ( x ; ξ ) β α , x R ,
has the properties of a cdf, with the corresponding pdf given by
f o ( x ; α , β , ξ ) = α β g ( x ; ξ ) G ( x ; ξ ) 2 e α β 1 1 G ( x ; ξ ) 2 [ 1 + G ( x ; ξ ) ] G ( x ; ξ ) β 2 G ( x ; ξ ) β α 1 , x R .
To the best of our knowledge, it is new in the literature (and out the scope of this paper).

3.2. Quantile Function with Some Related Measures and Functions

The quantile function (qf) of the NPTL-G family is expressed in the following result.
Proposition 2.
The qf of the NPTL-G family is given by
Q ( u ; α , β , ξ ) = Q G 1 1 β ln 1 1 u 1 α 1 ; ξ , u ( 0 , 1 ) ,
where Q G ( u ; ξ ) is the qf corresponding to G ( x ; ξ ) .
Proof. 
For the sake of simplicity, let us set x u = Q ( u ; α , β , ξ ) for u ( 0 , 1 ) . Then, by the definition of a qf, x u satisfies the non-linear equation: u = F ( x u ; α , β , ξ ) , implying that u = e α β 1 1 G ( x u ; ξ ) 2 e β 1 1 G ( x u ; ξ ) α ; hence, u 1 α = e β 1 1 G ( x u ; ξ ) 2 e β 1 1 G ( x u ; ξ ) , which is equivalent to solving the polynomial equation according to y: y 2 2 y + u 1 α = 0 , with y = e β 1 1 G ( x u ; ξ ) . By determining the two roots of this polynomial, keeping only the one in the unit interval (since y ( 0 , 1 ) ), we get y = 1 1 u 1 α . After some algebra, we get G ( x u ; ξ ) = 1 1 β ln 1 1 u 1 α 1 . The desired result follows by compounding with Q G ( u ; ξ ) , ending the proof of Proposition 2. □
From the qf, we can define several quantities of importance, providing distributional properties of the family. Some of them are presented below.
The three quartiles of the NPTL-G family are defined by Q 1 = Q ( 1 / 4 ; α , β , ξ ) , Q 2 = Q ( 1 / 2 ; α , β , ξ ) , and Q 3 = Q ( 3 / 4 ; α , β , ξ ) . In particular, the median of the NPTL-G family is given by
M = Q 2 = Q G 1 1 β ln 1 1 1 2 1 α 1 ; ξ .
Additionally, the inter-quartile range is given by I Q R = Q 3 Q 1 , allowing one to define the Galton coefficient of skewness and the Moors coefficient of kurtosis, given by, respectively,
S = Q 3 + Q 1 2 M I Q R
and
K = Q ( 7 / 8 ; α , β , ξ ) Q ( 5 / 8 ; α , β , ξ ) + Q ( 3 / 8 ; α , β , ξ ) Q ( 1 / 8 ; α , β , ξ ) I Q R .
See [28,29] for more details on these coefficients, respectively.
On the other hand, upon differentiation of Q ( u ; α , β , ξ ) according to u, the corresponding quantile density function is given by
q ( u ; α , β , ξ ) = u 1 α 1 2 α β 1 u 1 / α 1 1 u 1 α 1 1 β ln 1 1 u 1 α 2 × q G 1 1 β ln 1 1 u 1 α 1 ; ξ , u ( 0 , 1 ) ,
where q G ( u ; ξ ) is the quantile density function corresponding to G ( x ; ξ ) . Also, the hazard quantile function is defined by
H ( u ; α , β , ξ ) = 1 ( 1 u ) q ( u ; α , β , ξ ) = 2 α β 1 u 1 / α 1 1 u 1 α 1 1 β ln 1 1 u 1 α 2 ( 1 u ) u 1 α 1 × q G 1 1 β ln 1 1 u 1 α 1 ; ξ 1 , u ( 0 , 1 ) .
These functions have central roles in reliability. Further details can be found in [30].
Last but not least, the qf allows us to generate values from members of the NPTL-G family. This property will be used in Section 4.2 in the context of the NPTLILx distribution; i.e., with the qf given by Q G ( u ; θ ) = ( u 1 θ 1 ) 1 , u ( 0 , 1 ) , so
Q ( u ; α , β , θ ) = 1 1 β ln 1 1 u 1 α 1 θ 1 1 , u ( 0 , 1 ) .
As a numerical illustration, Table 1 shows the values of Q 1 , M, Q 3 , S, and K of the NPTLILx distribution for some parameter values.
We see in Table 1 that the effects of α , β , and θ on the quartiles are significant (we always have S > 0 so the distribution is right-skewed and moderate variations for K).

3.3. Series Expansion

The exp-G family of distributions, introduced by [31], is defined by the following cdf: G γ ( x ; ξ ) = G ( x ; ξ ) γ , x R , with γ > 0 . The corresponding pdf is given by
g γ ( x ; ξ ) = γ g ( x ; ξ ) G γ 1 ( x ; ξ ) , x R .
The interesting part of the exp-G family is to have well-known properties for a lot of baseline cdfs G ( x ; ξ ) . For instance, the member of the exp-G family defined with the inverse Lomax distribution as baseline with shape parameter θ becomes the inverse Lomax distribution with shape parameter γ θ .
The following result concerns a series expansion for the pdf of the NPTL-G family in terms of pdfs of the exp-G family.
Proposition 3.
We have the following series expansion:
f ( x ; α , β , ξ ) = k , , m = 0 + q = 1 m + ω k , , m , q g q ( x ; ξ ) ,
where
ω k , , m , q = α k m m + q 2 α k 1 ! ( α + k ) β ( 1 ) k + + m + q ,
with the notation: b a = b ( b 1 ) ( b a + 1 ) / a ! .
Proof. 
We first investigate a series expansion of F ( x ; α , β , ξ ) based on the Equation (1). Since e β 1 1 G ( x ; ξ ) / 2 ( 0 , 1 ) , the generalized binomial formula gives
2 e β 1 1 G ( x ; ξ ) α = k = 0 + α k 2 α k ( 1 ) k e k β 1 1 G ( x ; ξ ) .
On the other hand, thanks to the power series of the exponential function, we get
e ( α + k ) β 1 1 G ( x ; ξ ) = = 0 + 1 ! ( α + k ) β 1 1 G ( x ; ξ ) .
Now, it follows from the generalized and standard binomial formulas that
1 1 G ( x ; ξ ) = m = 0 + m ( 1 ) + m [ 1 G ( x ; ξ ) ] m + = m = 0 + q = 0 m + m m + q ( 1 ) + m + q G q ( x ; ξ ) .
By combining all the above equalities together, we obtain
F ( x ; α , β , ξ ) = k , , m = 0 + q = 0 m + ω k , , m , q G q ( x ; ξ ) .
Upon differentiation of F ( x ; α , β , ξ ) according to x, we get the desired result, by removing the term in q = 0 , which vanished. Proposition 3 is proven. □

3.4. General Moments with Some Related Measures and Functions

Let X be a random variable having the cdf given by (1) (defined on a probability space ( Ω , A , P ) , with an expectation denoted by E). Then, for any function ϕ ( x ) (such that all the following introduced quantities exist or converge), we have
E ϕ ( X ) = + ϕ ( x ) f ( x ; α , β , ξ ) d x = + ϕ ( x ) 2 α β g ( x ; ξ ) G ( x ; ξ ) 2 e α β 1 1 G ( x ; ξ ) 1 e β 1 1 G ( x ; ξ ) 2 e β 1 1 G ( x ; ξ ) α 1 d x .
Two equivalent expressions involving already introduced qfs are as follows:
E ϕ ( X ) = 0 1 ϕ [ Q G ( u ; ξ ) ] 2 α β u 2 e α β 1 1 u 1 e β 1 1 u 2 e β 1 1 u α 1 d u
and
E ϕ ( X ) = 0 1 ϕ Q G 1 1 β ln 1 1 u 1 α 1 ; ξ d u .
Numerical solutions exist to evaluate them for given G ( x ; ξ ) , ϕ ( x ) and α , β , and θ . Alternatively, we can consider Proposition 3, which implies that
E ϕ ( X ) = k , , m = 0 + q = 1 m + ω k , , m , q U q ( ϕ , G ) ,
where
U q ( ϕ , G ) = + ϕ ( x ) g q ( x ; ξ ) d x = q 0 1 u q 1 ϕ [ Q G ( u ; ξ ) ] d u .
In some circumstances, truncated sums can be considered for practical purposes; for a large integer K, the following approximation reveals to be tractable and efficient:
E ϕ ( X ) k , , m = 0 K q = 1 m + ω k , , m , q U q ( ϕ , G ) .
Some specific choices for ϕ ( x ) are of particular interest. Some of them are discussed below.
  • By taking ϕ ( x ) = x s , we get the s-th moment of X—i.e., μ s = E ( X s ) , including the mean of X, i.e., μ = μ 1 = E ( X ) —and allow the expression the variance of X; i.e., σ 2 = μ 2 ( μ 1 ) 2 .
  • By taking ϕ ( x ) = ( x μ ) s , we get the s-th central moment of X, i.e., μ s = E [ ( X μ ) s ] , allowing one to calculate the s-th general coefficient of X given by C s = μ s / σ s , among others. This coefficient is useful to investigate the skewness and kurtosis properties of X.
  • By taking ϕ ( x ) = e t x , we get the moment generation function of X according to the variable t; i.e., M ( t ) = E ( e t X ) . It is well-known that μ s = M ( t ) ( s ) t = 0 .
  • By taking ϕ ( x ) = e i t x , we get the characteristic function of X according to the variable t; i.e., φ ( t ) = E ( e i t X ) . In a same title of the cdf, the characteristic function entirely determines the NPTL-G family.
  • By taking ϕ ( x ) = ϕ y ( x ) = x s 1 { x y } , which is equal to x s if x y and 0 otherwise, we get the s-th incomplete moment of X according to the variable y; i.e., μ s ( y ) = E ( X s 1 { X y } ) . This function is useful to define mean deviations of X, the corresponding residual life function, Bonferroni and Lorenz curves, and others.
In the case of the NPTLILx distribution, since f ( x ; α , β , θ ) 2 α β 2 θ 2 x 3 when x + , the mean exists but the variance does not exist, nor do moments of order greater to 2 (there is no problem when x 0 ). However, all the incomplete moments exist for any fixed y > 0 . In this regard, Table 2 provides the four first incomplete moments for X with y = 1000 .

3.5. Shannon Entropy

Here, we study the Shannon entropy of the NPTL-G family as defined by [32]. We recall that the Shannon entropy of a random variable measures the amount of uncertainty for the outcome of this variable. A high entropy reveals a high degree of uncertainty.
Now, let X be a random variable having the cdf given by (1). Then, the Shannon entropy of X is defined by
η = E { ln [ f ( X ; α , β , ξ ) } = + ln [ f ( x ; α , β , ξ ) ] f ( x ; α , β , ξ ) d x .
By the use of any mathematical software, for a given baseline cdf G ( x ; ξ ) , ϕ ( x ) and α , β , and θ , we can determine this integral. Another approach consists of developing η by the use of the pdf given by (2):
η = ln ( 2 ) ln ( α ) ln ( β ) α β E ln [ g ( X ; ξ ) ] + 2 E ln [ G ( X ; ξ ) ] + α β E 1 G ( X ; ξ ) E ln 1 e β 1 1 G ( X ; ξ ) ( α 1 ) E ln 2 e β 1 1 G ( X ; ξ ) .
Some expectation terms can be expressed by using (5) with an appropriate function ϕ ( x ) as soon as U q ( ϕ , G ) exists and the sums converge.
In the context of the NPTLILx distribution, some values of η are collected in Table 3 for some parameter values.
In Table 3, the values belongs to the wide interval [ 12.8 , 3.55 ] , meaning that α , β , and θ have an important impact on the amount of information quantified by η .

4. Estimation with Numerical Results

In this section, we investigate the NPTLILx model characterized by the cdf given by (3). Thanks to its attractive theoretical and practical properties, the maximum likelihood method is used to estimate the parameters α , β , and θ . Numerical results attest to the efficiency of the estimates obtained.
Hereafter, we consider a random variable X following the NPTLILx distribution with parameters α , β , and θ .

4.1. Maximum Likelihood Estimation

Let x 1 , , x n be a random sample of size n of X. Then, by using the pdf given by (4), the likelihood and log-likelihood functions are, respectively, given by
L ( α , β , θ ) = i = 1 n f ( x i ; α , β , θ ) = ( 2 α β θ ) n i = 1 n x i 2 1 + x i 1 θ 1 e α β 1 1 + x i 1 θ 1 e β 1 1 + x i 1 θ 2 e β 1 1 + x i 1 θ α 1
and
( α , β , θ ) = ln [ L ( α , β , θ ) ] = n ln ( 2 ) + n ln ( α ) + n ln ( β ) + n ln ( θ ) 2 i = 1 n ln ( x i ) + ( θ 1 ) i = 1 n ln ( 1 + x i 1 ) + α β i = 1 n 1 ( 1 + x i 1 ) θ + i = 1 n ln 1 e β 1 ( 1 + x i 1 ) θ + ( α 1 ) i = 1 n ln 2 e β 1 ( 1 + x i 1 ) θ .
The maximum likelihood estimates (MLEs) of α , β , and θ , say α ^ , β ^ , and θ ^ , respectively, are defined such that L ( α ^ , β ^ , θ ^ ) = max ( α , β , θ ) ( 0 , + ) 3 L ( α , β , θ ) or ( α ^ , β ^ , θ ^ ) = max ( α , β , θ ) ( 0 , + ) 3 ( α , β , θ ) . Let us work with the function ( α , β , θ ) for the sake of simplicity. Since ( α , β , θ ) is differentiable with respect to α , β , and θ , the MLEs can obtained by solving the non-linear equations defined by the first partial derivatives of ( α , β , θ ) with respect to α , β , and θ equal to 0, with
( α , β , θ ) α = n α + β i = 1 n 1 ( 1 + x i 1 ) θ + i = 1 n ln 2 e β 1 ( 1 + x i 1 ) θ ,
( α , β , θ ) β = n β + α i = 1 n 1 ( 1 + x i 1 ) θ i = 1 n 1 ( 1 + x i 1 ) θ e β 1 ( 1 + x i 1 ) θ 1 e β 1 ( 1 + x i 1 ) θ ( α 1 ) i = 1 n 1 ( 1 + x i 1 ) θ e β 1 ( 1 + x i 1 ) θ 2 e β 1 ( 1 + x i 1 ) θ
and
( α , β , θ ) θ = n θ + i = 1 n ln ( 1 + x i 1 ) α β i = 1 n ( 1 + x i 1 ) θ ln ( 1 + x i 1 ) + β i = 1 n ( 1 + x i 1 ) θ ln ( 1 + x i 1 ) e β 1 ( 1 + x i 1 ) θ 1 e β 1 ( 1 + x i 1 ) θ + β ( α 1 ) i = 1 n ( 1 + x i 1 ) θ ln ( 1 + x i 1 ) e β 1 ( 1 + x i 1 ) θ 2 e β 1 ( 1 + x i 1 ) θ .
The complexity of these expressions do not allow us to provide closed-forms for the MLEs. However, several numerical solutions exist to maximize ( α , β , θ ) based on Newton–Raphson algorithms, one of which is employed in this study.
The corresponding Fisher information matrix we observed is given by
J ( α , β , θ ) = 2 ( α , β , θ ) α 2 2 ( α , β , θ ) α β 2 ( α , β , θ ) α θ 2 ( α , β , θ ) β α 2 ( α , β , θ ) β 2 2 ( α , β , θ ) β θ 2 ( α , β , θ ) θ α 2 ( α , β , θ ) θ β 2 ( α , β , θ ) θ 2 .
(the elements of J ( α , β , θ ) are upon request from the authors). When n is large, the distribution of the subjacent random vector behind ( α ^ , β ^ , θ ^ ) can be approximated by a three dimensional normal distribution with mean vector ( α , β , θ ) and covariance matrix J ( α ^ , β ^ , θ ^ ) 1 . By denoting v α ^ , v β ^ and v θ ^ , the diagonal elements of this matrix, we are able to construct asymptotic confidence intervals for α , β , and θ . Indeed, with the adopted notations, the asymptotic (equitailed) confidence intervals (CIs) of α , β , and θ at the level 100 ( 1 γ ) % are given by, respectively,
C I α = α ^ z γ / 2 v α ^ , α ^ + z γ / 2 v α ^ , C I β = β ^ z γ / 2 v β ^ , β ^ + z γ / 2 v β ^
and
C I θ = θ ^ z γ / 2 v θ ^ , θ ^ + z γ / 2 v θ ^ ,
where z γ / 2 is the upper γ / 2 -th percentile of the normal distribution N ( 0 , 1 ) . For practical purposes, if lower bounds of these intervals are negative, we can put it at 0, since all the parameters are supposed to be positive. All the technical details can be found in [33].

4.2. Numerical Results

Here, we provide a simulation study to show the nice behavior of the MLEs for the NPTLILx model presented in the subsection above. First of all, let us mention that a random sample from X can be obtained by the use of the qf: for any random sample of size n from the uniform distribution U ( 0 , 1 ) , say u 1 , , u n , the corresponding random sample of size n of X is given by x 1 , , x n with x i = Q ( u i ; α , β , θ ) .
From N random samples of X, let ϵ be either α , β , or θ and ϵ ^ i be the MLE of ϵ constructed from the i-th sample. Then, we define the (mean) MLE, bias, and mean square error (MSE) by, respectively,
M L E ^ ϵ ( n ) = 1 N i = 1 N ϵ ^ i , B i a s ^ ϵ ( n ) = M L E ^ ϵ ( n ) ϵ , M S E ^ ϵ ( n ) = 1 N i = 1 N ( ϵ ^ i ϵ ) 2 .
Additionally, the asymptotic (mean) confidence intervals of α , β , and θ at the level 100 ( 1 γ ) % can be determined. We define the (mean) lower bounds (LBs), (mean) upper bounds (UBs), and (mean) average length (ALs) by, respectively,
L B ^ ϵ ( n ) = M L E ^ ϵ ( n ) z γ / 2 V ^ ϵ ( n ) , U P ^ ϵ ( n ) = M L E ^ ϵ ( n ) + z γ / 2 V ^ ϵ ( n ) , A L ^ ϵ ( n ) = 2 z γ / 2 V ^ ϵ ( n ) ,
where V ^ ϵ ( n ) = ( 1 / N ) i = 1 N v ϵ ^ i . For the purposes of this study, we consider the levels 90 % and 95 % , so z 0.05 = 1.644854 and z 0.025 = 1.959964 , respectively. The software Mathematica 9 was employed.
Our simulation study was based on the the following plan.
  • N = 1000 random samples of size n = 100 , 200, 300, and 1000 are to be generated from X.
  • Values of the true parameters ( α , β , θ ) are taken as, in order, ( 0.5 , 0.1 , 0.5 ) , ( 1.5 , 0.5 , 0.5 ) , and ( 1.8 , 0.4 , 1.2 ) .
  • The MLEs, MSEs, biases, LBs, UBs, and ALs for the selected values of the parameters are to be calculated.
Numerical outcomes are listed in Table 4, Table 5 and Table 6.
From Table 4, Table 5 and Table 6, we can see that, when n increases, biases, MSEs, and ALs decrease. This observation is consistent with the well-known convergence properties of the MLEs.

5. Data Analysis

In this section, we prove the flexibility of the NPTLILx model by analyzing two practical datasets. The fits of the NPTLILx model are compared to the competitive models listed in Table 7. The common point of all of them is the use the inverse Lomax distribution as the baseline distribution.
Except the former inverse Lomax distribution, the considered models possess three or four parameters. The comparison of these models was performed by using the following well-known statistical benchmarks: CVM (Cramér–von Mises); AD (Anderson–Darling); KS (Kolmogorov–Smirnov) statistic with the corresponding p-value, minus log-likelihood ( ^ ) ; AIC (Akaike information criterion); CAIC (corrected Akaike information criterion); BIC (Bayesian information criterion); and HQIC (Hannan–Quinn information criterion). For the CVM, AD, KS, ( ^ ) , AIC, CAIC, BIC, and HQIC, the smaller the value is, the better the fit to the data. Additionally, the higher the p-values of the KS test are, the better the fit to the data. All these measures were computed by using the R software.
Dataset I: The first data refer to [23]. It consists of 40 observations of the active repair times for airborne communication transceiver. The unit is the hour. The data are: 0.50, 0.60, 0.60, 0.70, 0.70, 0.70, 0.80, 0.80, 1.00, 1.00, 1.00, 1.00, 1.10, 1.30, 1.50, 1.50, 1.50, 1.50, 2.00, 2.00, 2.20, 2.50, 2.70, 3.00, 3.00, 3.30, 4.00, 4.00, 4.50, 4.70, 5.00, 5.40, 5.40, 7.00, 7.50, 8.80, 9.00, 10.20, 22.00, 24.50.
A basic statistical description of the data gives: n = 40 , mean = 4.01 , standard deviation = 5.17 , median = 2.1 , skewness = 2.62 , and kurtosis = 7.02 . One can notice that the data are skewed to the right with a high kurtosis.
Dataset II: Next, we use the actual taxes dataset as described in [24]. The data consist of the monthly actual taxes revenue in Egypt from January 2006 to November 2010. The unit is the 1000 million Egyptian pounds. The data are: 5.9, 20.4, 14.9, 16.2, 17.2, 7.8, 6.1, 9.2, 10.2, 9.6, 13.3, 8.5, 21.6, 18.5,5.1,6.7, 17, 8.6, 9.7, 39.2, 35.7, 15.7, 9.7, 10, 4.1, 36, 8.5, 8, 9.2, 26.2, 21.9, 16.7, 21.3, 35.4, 14.3, 8.5, 10.6, 19.1, 20.5, 7.1, 7.7, 18.1, 16.5, 11.9, 7, 8.6, 12.5, 10.3, 11.2, 6.1, 8.4, 11, 11.6, 11.9, 5.2, 6.8, 8.9, 7.1, 10.8.
A basic statistical description of the data gives: n = 59 , mean = 13.49 , standard deviation = 8.05 , median = 10.6 , skewness = 1.57 , and kurtosis = 2.08 . Thus, these data are skewed to the right with a moderate kurtosis.
The graphical and numerical analyses of these two datasets are as follows. Figure 3 presents the total test time (TTT) plots of the two datasets. The first plot shows a convex curve, indicating that a decreasing hrf for the fitting model is appropriate for Data set I, whereas the second plot shows a concave curve, indicating that an increasing hrf for the fitting model is appropriate for Data set II. These cases are covered by the NPTLILx model, as shown in Figure 2.
Table 8 and Table 9 present the CVM, AD, KS, and the related p-value, and the MLEs of the models’ parameters for Datasets I and II, respectively. The obtained p-values indicate that the NPTLILx model is the best. Table 10 and Table 11 communicate the ^ , AIC, BIC, CAIC, BIC, and HQIC of the models for Datasets I and II, respectively. Since the smallest values are obtained for the NPTLILx model, it can be considered the best with these criteria. The estimated pdfs and cdfs for the considered models are displayed in Figure 4 and Figure 5 for Datasets I and II, respectively. The plots of the estimated pdfs are visually refined via an individual treatment in Figure 6 and Figure 7. In order to give another point of view, we illustrate the adequateness of the models via the use of probability–probability (PP) plots in Figure 8 and Figure 9, for Datasets I and II, respectively. In particular, for Dataset II, in view of the perfect adjustment of the scatter plot by the PP line, it is clear that the NPTLILx model provides a better fit in comparison to the other models. To resume, the NPTLILx model reveals itself to be the more appropriate model for the two datasets, illustrating its applicability in a concrete setting.
We end this section by providing some additional graphical and numerical elements on the NPTLILx model, related to the quantities presented in Section 4.1. To illustrate the uniqueness of the MLEs of α , β and θ , the profiles of the log-likelihood function are proposed in Figure 10 and Figure 11 for Datasets I and II, respectively. The Fisher information matrices of the NPTLILx model taken at the MLEs for Datasets I and II are, respectively, given by
J I = 0.0080 0.0354 0.0018 0.0354 5.8404 0.0216 0.0018 0.0216 0.4128 , J I I = 1.0770 1.2877 0.0627 1.2877 1.2327 0.0654 0.0627 0.0654 0.0590 .
Then, the confidence intervals for α , β , and θ at the levels 90 % and 95 % are provided in Table 12.

6. Conclusion and Perspectives

In this paper, we introduced and studied a new general family of distributions, called the NPTL-G family, based on the so-called power Topp–Leone-G and inverse exponential-G families. Various mathematical properties were presented, including stochastic ordering, quantile function and related measures, general moments and related measures, and the Shannon entropy, with discussions. Then, we payed special attention to a member of the family defined with the inverse Lomax distribution, called the NPTLILx distribution. The estimation of the unknown model parameters was done with the maximum likelihood method, with numerical guarantees on their behavior via a simulation study. The applicability of the NPTLILx model was then illustrated by the consideration of two practical datasets. It was then proven that the NPTLILx model is a serious alternative to other models, also using the inverse Lomax distribution as the baseline. Future work will include the constructions of various regression models, Bayesian estimation of the parameters, and analyses of new datasets. Thanks to its numerous qualities, we believe that the NPTL-G family can be helpful for the practitioner, for statistical analyses beyond the scope of this paper.
Among the interesting perspectives of work, one could investigate the confidence bounds and supersaturation properties of the cdfs of the members of the NPTL-G family, which are useful for choosing an appropriate model for given data, following the spirit of [38,39,40,41,42]. All these aspects need further investigations that we leave for future works.

Author Contributions

R.A.R.B., F.J., C.C., and M.E. contributed equally to this work.

Funding

This work was funded by the Deanship of Scientific Research (DSR), King AbdulAziz University, Jeddah, under grant number DF-277-305-1441.

Acknowledgments

We would like to thank the two reviewers for their thoughtful efforts towards improving our manuscript. The authors gratefully acknowledge the DSR for technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eugene, N.; Lee, C.; Famoye, F. The beta-normal distribution and its applications. Commun. Stat. Theory Methods 2002, 31, 497–512. [Google Scholar] [CrossRef]
  2. Cordeiro, G.M.; de Castro, M. A new family of generalized distribution. J. Stat. Comput. Simul. 2002, 81, 883–898. [Google Scholar] [CrossRef]
  3. Bourguignon, M.; Silva, R.B.; Cordeiro, G.M. The Weibull-G family of probability distributions. J. Data Sci. 2014, 12, 53–68. [Google Scholar]
  4. Elgarhy, M.; Hassan, A.S.; Rashed, M. Garhy-Generated Family of Distributions with Application. Math. Theory Model. 2016, 6, 1–15. [Google Scholar]
  5. Hassan, A.S.; Elgarhy, M.; Shakil, M. Type II half logistic family of distributions with applications. Pak. J. Stats. Oper. Res. 2017, 13, 245–264. [Google Scholar]
  6. Yousof, H.M.; Alizadeh, M.; Jahanshahi, S.M.A.; Ramires, T.G.; Ghosh, I.; Hamedani, G.G. The Transmuted Topp-Leone G family of distributions: Theory, characterizations and applications. J. Data Sci. 2017, 15, 723–740. [Google Scholar]
  7. Cordeiro, G.M.; Alizadeh, M.; Ozel, G.; Hosseini, B.; Ortega, E.M.M.; Altun, E. The generalized odd log-logistic family of distributions: Properties, regression models and applications. J. Stat. Comput. Simul. 2017, 87, 908–932. [Google Scholar] [CrossRef]
  8. Haq, M.A.; Elgarhy, M. The odd Fréchet-G family of probability distributions. J. Stat. Appl. Probab. 2018, 7, 189–203. [Google Scholar] [CrossRef]
  9. Hassan, A.S.; Nassar, S.G. Power Lindley-G family. Ann. Data Sci. 2018, 6, 189–210. [Google Scholar]
  10. Reyad, H.; Korkmaz, M.C.; Afify, A.Z.; Hamedani, G.G.; Othman, S. The Fréchet Topp-Leone-G family of distributions: Properties, characterizations and applications. Ann. Data Sci. 2019. [Google Scholar] [CrossRef]
  11. Reyad, H.M.; Alizadeh, M.; Jamal, F.; Othman, S.; Hamedani, G.G. The exponentiated generalized Topp Leone-G family of distributions: Properties and applications. Pak. J. Stats. Oper. Res. 2019, 15, 1–24. [Google Scholar] [CrossRef]
  12. Bantan, R.A.; Jamal, F.; Chesneau, C.; Elgarhy, M. Truncated inverted Kumaraswamy generated family of distributions with applications. Entropy 2019, 21, 1089. [Google Scholar] [CrossRef]
  13. Tahir, M.H.; Cordeiro, G.M. Compounding of distributions: A survey and new generalized classes. J. Stat. Distribut. Appl. 2016, 3, 1–35. [Google Scholar] [CrossRef]
  14. Al-Shomrani, A.; Arif, O.; Shawky, A.; Hanif, S.; Shahbaz, M.Q. Topp-Leone family of distributions: Some properties and application. Pak. J. Stat. Oper. Res. 2016, 12, 443–451. [Google Scholar] [CrossRef]
  15. Rezaei, S.; Sadr, B.B.; Alizadeh, M.; Nadarajah, S. Topp-Leone generated family of distributions: Properties and applications. Commun. Stat. Theory Methods 2016, 46, 2893–2909. [Google Scholar] [CrossRef]
  16. Mahdavi, A. Generalized Topp-Leone family of distributions. J. Biostat. Epidemiol. 2017, 3, 65–75. [Google Scholar]
  17. Elgarhy, M.; Nasir, M.A.; Jamal, F.; Ozel, G. The type II Topp-Leone generated family of distributions: Properties and applications. J. Stat. Manag. Syst. 2018, 21, 1529–1551. [Google Scholar] [CrossRef]
  18. Hassan, A.S.; Elgarhy, M.; Ahmad, Z. Type II generalized Topp-Leone family of distributions: Properties and applications. J. Data Sci. 2019, 17, 638–659. [Google Scholar]
  19. Kumar, D.; Singh, S.K.; Singh, U. Life time distributions: Derived from some minimum guarantee distribution. Sohag J. Math. 2017, 4, 7–11. [Google Scholar] [CrossRef]
  20. Alrajhi, S. The odd Fréchet inverse exponential distribution with application. J. Nonlinear Sci. Appl. 2019, 12, 535–542. [Google Scholar] [CrossRef]
  21. Nasiru, S. Extended odd Fréchet-G family of distributions. J. Probab. Stat. 2018, 1, 1–12. [Google Scholar] [CrossRef]
  22. Chesneau, C.; Djibrila, S. The generalized odd inverted exponential-G family of distributions: properties and applications. Eurasian Bullet. Math. 2019, in press. [Google Scholar]
  23. Jorgensen, B. Statistical Properties of the Generalized Inverse Gaussian Distribution; Springer: New York, NY, USA, 1982. [Google Scholar]
  24. Mead, M.E. A new generalization of Burr XII distribution. J. Stat. 2014, 12, 53–73. [Google Scholar]
  25. Kleiber, C. Lorenz ordering of order statistics from log-logistic and related distributions. J. Stat. 2004, 120, 13–19. [Google Scholar] [CrossRef]
  26. Kleiber, C.; Kotz, S. Statistical Size Distributions in Economics and Actuarial Sciences; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2003. [Google Scholar]
  27. Singh, S.K.; Singh, U.; Kumar, D. Bayes estimators of the reliability function and parameter of inverted exponential distribution using informative and non-informative priors. J. Stat. Comput. Simul. 2013, 83, 2258–2269. [Google Scholar] [CrossRef]
  28. Galton, F. Inquiries into Human Faculty and Its Development; Macmillan and Company: London, UK, 1883. [Google Scholar]
  29. Moors, J.J.A. A quantile alternative for kurtosis. J. R. Stat. Soc. 1988, 37, 25–32. [Google Scholar] [CrossRef]
  30. Nair, N.U.; Sankaran, P.G. Quantile-based reliability analysis. Commun. Stat. Theory Methods 2009, 38, 222–232. [Google Scholar] [CrossRef]
  31. Gupta, R.D.; Kundu, D. Exponentiated exponential family: An alternative to gamma and Weibull distributions. Biometric. J. 2001, 43, 117–130. [Google Scholar] [CrossRef]
  32. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  33. Casella, G.; Berger, R.L. Statistical Inference; Brooks/Cole Publishing Company: Bel Air, CA, USA, 1990. [Google Scholar]
  34. Rahman, J.; Aslam, M.; Ali, S. Estimation and prediction of inverse Lomax model via Bayesian approach. Caspian J. Appl. Sci. Res. 2013, 2, 43–56. [Google Scholar]
  35. Hassan, A.S.; Abd-Allah, M. On the Inverse Power Lomax Distribution. Ann. Data Sci. 2019, 6, 259–278. [Google Scholar] [CrossRef]
  36. Hassan, A.S.; Ismail, D.M. Parameter Estimation of Topp-Leone Inverse Lomax Distribution. J. Modern Appl. Stat. Methods 2019, in press. [Google Scholar]
  37. Hassan, A.S.; Mohamed, R.E. Weibull Inverse Lomax Distribution. Pak. J. Stat. Oper. Res. 2019, 15, 587–603. [Google Scholar]
  38. Sendov, B. Hausdorff Approximations; Wolters Kluwer: Alphen aan den Rijn, The Netherlands, 1990. [Google Scholar]
  39. Iliev, A.; Rahnev, A.; Kyurkchiev, N.; Markov, S. A study on the unit-logistic, unit-Weibull and Topp-Leone cumulative sigmoids. Biomath. Commun. 2019, 6, 1–15. [Google Scholar] [CrossRef] [Green Version]
  40. Kyurkchiev, N. Uniform approximation of the generalized cut function by Erlang cumulative distribution function. Application in applied insurance mathematics. Int. J. Theor. Appl. Math. 2016, 2, 40–44. [Google Scholar]
  41. Kyurkchiev, N. Mathematical Concepts in Insurance and Reinsurance: Some Moduli in Programming Environment MATHEMATICA; LAP LAMBERT Academic Publishing: Riga, Latvia, 2016; ISBN 978-3-659-96906-5. [Google Scholar]
  42. Kyurkchiev, N.; Iliev, A.; Markov, S. Some Techniques for Recurrence Generating of Activation Functions: Some Modeling and Approximation Aspects; LAP LAMBERT Academic Publishing: Riga, Latvia, 2017; ISBN 978-3-330-33143-3. [Google Scholar]
Figure 1. Plots of some probability density functions (pdfs) of the new power Topp–Leone inverse Lomax (NPTLILx) distribution.
Figure 1. Plots of some probability density functions (pdfs) of the new power Topp–Leone inverse Lomax (NPTLILx) distribution.
Entropy 21 01177 g001
Figure 2. Plots of some hazard rate functions (hrfs) of the NPTLILx distribution.
Figure 2. Plots of some hazard rate functions (hrfs) of the NPTLILx distribution.
Entropy 21 01177 g002
Figure 3. Total test time (TTT) plots for Datasets I and II, respectively.
Figure 3. Total test time (TTT) plots for Datasets I and II, respectively.
Entropy 21 01177 g003
Figure 4. Plots for the estimated pdfs and cdfs for Dataset I.
Figure 4. Plots for the estimated pdfs and cdfs for Dataset I.
Entropy 21 01177 g004
Figure 5. Plots for the estimated pdfs and cdfs for Dataset II.
Figure 5. Plots for the estimated pdfs and cdfs for Dataset II.
Entropy 21 01177 g005
Figure 6. Plots of the pdfs estimated for Dataset I.
Figure 6. Plots of the pdfs estimated for Dataset I.
Entropy 21 01177 g006
Figure 7. Plots of the pdfs estimated for Dataset II.
Figure 7. Plots of the pdfs estimated for Dataset II.
Entropy 21 01177 g007
Figure 8. Probability–probability (PP) plots of considered models for Dataset I.
Figure 8. Probability–probability (PP) plots of considered models for Dataset I.
Entropy 21 01177 g008
Figure 9. PP plots of considered models for Dataset II.
Figure 9. PP plots of considered models for Dataset II.
Entropy 21 01177 g009
Figure 10. Profiles of the log-likelihood function of the NPTLILx model for Dataset I.
Figure 10. Profiles of the log-likelihood function of the NPTLILx model for Dataset I.
Entropy 21 01177 g010
Figure 11. Profiles of the log-likelihood function of the NPTLILx model for Dataset II.
Figure 11. Profiles of the log-likelihood function of the NPTLILx model for Dataset II.
Entropy 21 01177 g011
Table 1. The values of Q 1 , M, Q 3 , S, and K of the NPTLILx distribution for some parameter values.
Table 1. The values of Q 1 , M, Q 3 , S, and K of the NPTLILx distribution for some parameter values.
( α , β , θ ) Q 1 M Q 3 SK
( 0.5 , 0.5 , 0.5 ) 0.01620.04130.11080.47031.8976
( 1.5 , 0.5 , 0.5 ) 0.06670.13770.30000.39171.7017
( 2.5 , 0.5 , 0.5 ) 0.11500.22000.44640.36641.6542
( 3.5 , 0.5 , 0.5 ) 0.15940.29220.57010.35331.6323
( 5.0 , 0.5 , 0.5 ) 0.22000.38770.72980.34211.6151
( 5.0 , 1.0 , 0.5 ) 0.55130.91681.63300.32441.5894
( 5.0 , 2.0 , 0.5 ) 1.26242.01743.47240.31671.5807
( 5.0 , 3.0 , 0.5 ) 1.98963.13085.32070.31481.5788
( 5.0 , 5.0 , 0.5 ) 3.45635.36669.02310.31371.5778
( 5.0 , 5.0 , 1.0 ) 7.380811.211818.53310.31301.5772
( 5.0 , 5.0 , 2.0 ) 15.245822.912837.55970.31281.5771
( 5.0 , 5.0 , 3.0 ) 23.114334.616356.58770.31281.5771
( 5.0 , 5.0 , 5.0 ) 38.853458.024694.64450.31281.5771
Table 2. The values for the first four incomplete moments of the NPTLILx distribution; i.e., μ s ( y ) with s = 1 , 2 , 3 , 4 , with y = 1000 , for some parameter values.
Table 2. The values for the first four incomplete moments of the NPTLILx distribution; i.e., μ s ( y ) with s = 1 , 2 , 3 , 4 , with y = 1000 , for some parameter values.
( α , β , θ ) μ 1 ( 1000 ) μ 2 ( 1000 ) μ 3 ( 1000 ) μ 4 ( 1000 )
( 0.5 , 0.5 , 0.5 ) 0.13300.411462.015031180.1500
( 1.0 , 0.5 , 0.5 ) 0.22890.8177124.026962360.2800
( 3.0 , 0.5 , 0.5 ) 0.50372.4083372.0441187080.6
( 0.5 , 1.0 , 0.5 ) 0.32161.5544247.4997124628
( 0.5 , 3.0 , 0.5 ) 1.164812.25362209.28201118342
( 0.5 , 0.5 , 1.0 ) 0.44991.6965248.6065124813.2
( 0.5 , 0.5 , 3.0 ) 3.637641.68066308.09303138854
( 2.0 , 2.0 , 0.5 ) 1.882622.15193941.26401991066
( 2.0 , 2.0 , 2.0 ) 8.7637288.420162161.9131710832
( 2.0 , 2.0 , 15 ) 65.79219473.51931490441707182402
( 2.0 , 2.0 , 30 ) 125.144228393.86113310026455199311
( 2.0 , 5.0 , 20 ) 180.989755135.342560472615505315569
( 5.0 , 5.0 , 20 ) 268.5182106564.85609250735760970261
( 10 , 10 , 10 ) 337.2701160356.59326339262533636278
Table 3. The values of the Shannon entropy of the NPTLILx distribution for some parameter values.
Table 3. The values of the Shannon entropy of the NPTLILx distribution for some parameter values.
α β θ η
0.50.50.5 1.3743
1.00.50.5 0.7246
2.00.50.5 0.1757
3.00.50.50.1094
5.00.50.50.4405
100.50.50.8528
101.00.51.5895
102.00.52.2907
105.00.53.1614
1080.53.5497
0.50.50.1 12.8618
5.00.50.1 3.9876
5.05.00.11.1590
Table 4. Maximum likelihoods (MLEs), biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 0.5 , β = 0.1 , θ = 0.5 ) .
Table 4. Maximum likelihoods (MLEs), biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 0.5 , β = 0.1 , θ = 0.5 ) .
nPar.MLBiasMSE90%95%
LBUBALLBUBAL
100 α 0.290 0.210 0.105 0.528 1.1091.637 0.685 1.2661.951
β 0.1630.0630.016 0.090 0.4170.508 0.139 0.4660.605
θ 0.5720.0720.0120.3640.7800.4160.3240.8200.496
200 α 0.304 0.196 0.103 0.329 0.9371.266 0.450 1.0581.508
β 0.1460.0460.008 0.031 0.3240.354 0.065 0.3580.422
θ 0.5670.0670.0110.3860.7490.3630.3510.7840.433
300 α 0.340 0.160 0.084 0.264 0.9441.207 0.379 1.0591.438
β 0.1370.0370.0060.0210.2540.233 0.001 0.2760.277
θ 0.5480.0480.0070.4140.6830.2690.3880.7090.321
1000 α 0.342 0.158 0.074 0.064 0.7470.811 0.141 0.8250.966
β 0.1240.0240.0020.0470.2000.1530.0330.2150.182
θ 0.5450.0450.0050.4350.6550.2200.4140.6760.263
Table 5. MLEs, biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 1.5 , β = 0.5 , θ = 0.5 ) .
Table 5. MLEs, biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 1.5 , β = 0.5 , θ = 0.5 ) .
nPar.MLBiasMSE90%95%
LBUBALLBUBAL
100 α 2.0420.5425.123 211.872 215.955427.827 252.834 256.918509.752
β 0.373 0.127 0.055 8.964 9.71118.674 10.752 11.49922.250
θ 0.6180.1180.060 30.579 31.81462.393 36.553 37.78874.341
200 α 1.8430.3431.449 81.290 84.975166.265 97.209 100.894198.103
β 0.442 0.058 0.015 4.239 5.1249.363 5.136 6.02111.156
θ 0.5470.0470.019 10.112 11.20621.317 12.153 13.24725.400
300 α 1.6670.1670.874 45.858 49.19295.050 54.959 58.293113.252
β 0.458 0.042 0.012 2.019 2.9354.954 2.493 3.4105.903
θ 0.5420.0420.013 5.733 6.81812.550 6.934 8.01914.953
1000 α 1.358 0.142 0.460 2.292 5.0087.300 2.991 5.7078.698
β 0.5340.0340.0020.1840.8840.7000.1170.9510.834
θ 0.5340.0340.0100.2790.7880.5090.2300.8370.607
Table 6. MLEs, biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 1.8 , β = 0.4 , θ = 1.2 ) .
Table 6. MLEs, biases, MSEs, LBs, UBs, and ALs of the NPTLILx model for ( α = 1.8 , β = 0.4 , θ = 1.2 ) .
nPar.MLBiasMSE90%95%
LBUBALLBUBAL
100 α 2.8431.0431.889 1324.960 1330.6402655.600 1579.220 1584.9003164.120
β 0.6390.2390.675 94.183 95.462189.645 112.341 113.619225.960
θ 1.5730.3730.660 252.661 255.807508.468 301.344 304.490605.834
200 α 1.095 0.705 0.523 359.428 361.618721.046 428.464 430.654859.119
β 0.331−0.0690.014 27.830 28.49356.324 33.223 33.88667.109
θ 1.5710.3710.183 109.805 112.947222.752 131.132 134.274265.406
300 α 1.150 0.650 0.519 306.474 309.685616.159 365.468 368.679734.147
β 0.4210.0210.004 18.962 19.80438.766 22.674 23.51546.189
θ 1.3900.1900.165 92.650 95.370188.020 110.652 113.372224.024
1000 α 1.260 0.540 0.496 138.760 141.281280.040 165.572 168.093333.665
β 0.387 0.013 0.001 9.334 10.00719.341 11.186 11.85923.044
θ 1.3610.1610.150 45.112 48.29593.407 54.055 57.238111.293
Table 7. The competitive models considered.
Table 7. The competitive models considered.
DistributionReference
Inverse Lomax (ILx)[34]
Inverse Power Lomax (PILx)[35]
Topp–Leone Inverse Lomax (TILx)[36]
Weibull Inverse Lomax (WILx)[37]
Table 8. Goodness-of-fit measures, MLEs, and SEs for Dataset I.
Table 8. Goodness-of-fit measures, MLEs, and SEs for Dataset I.
ModelCVMADKSp-ValueMLEs with SEs (in Parentheses)
NPTLILx0.05500.34620.09430.8683 α β θ
0.068211.99021.5651
(0.0897)(2.4167)(0.6425)
WILx0.15221.08520.17840.1566ab λ β
0.00260.88670.01850.2581
(0.0006)(0.1938)(0.0285)(0.7650)
TILx0.06850.45780.11080.7100 α β λ
37.83242.98790.1676
(9.7598)(1.9941)(0.2682)
PILx0.10790.65820.12720.5369 α β λ
0.11306.85940.0571
(0.1242)(6.4470)(0.2480)
ILx0.06320.40650.09810.8355 α β
0.20038.2426
(0.1372)(5.1671)
Table 9. Goodness-of-fit measures, MLEs, and SEs for Dataset II.
Table 9. Goodness-of-fit measures, MLEs, and SEs for Dataset II.
ModelCVMADKSp-ValueMLEs with SEs (in Parentheses)
NPTLILx0.03570.26980.06150.9786 α β θ
14.43610.43785.0301
(1.0378)(1.1103)(0.2431)
WILx0.23631.48290.3248 7.79 × 10 6 ab λ β
0.00211.14040.01723.9985
(0.0002)(0.1202)(0.0064)(3.1371)
TILx0.03980.27010.09980.5988 α β λ
50.45790.090815.8617
(6.7245)(0.2094)(3.9080)
PILx0.11330.64400.14470.1689 α β λ
1.15452.3262300.7315
(0.4472)(0.2880)(121.3061)
ILx0.05290.30750.2928 8.06 × 10 5 α β
0.146471.1473
(0.2569)(23.8010)
Table 10. The values of ^ , AIC, and KS with its p-value for Dataset I.
Table 10. The values of ^ , AIC, and KS with its p-value for Dataset I.
Model ^ AICCAICBICHQIC
NPTLILx88.8229183.6459184.3125188.7125185.4778
WILx98.2937204.5874205.7303211.3429207.0300
TILx90.5459187.0919187.7586192.1586188.9239
PILx90.5908187.1817187.8484192.2483189.0136
ILx91.3612186.7226187.0469190.1003187.9439
Table 11. The values of ^ , AIC, and KS with its p-value for Dataset II.
Table 11. The values of ^ , AIC, and KS with its p-value for Dataset II.
Model ^ AICCAICBICHQIC
NPTLILx189.2811384.5622384.9985390.7948386.9951
WILx219.6212447.2424447.9832455.5526450.4864
TILx190.3769386.7538387.1902392.9864389.1868
PILx195.1056396.2113396.6476402.4439398.6442
ILx211.6436427.2872427.5015431.4422428.9091
Table 12. Confidence intervals for the parameters of the NPTLILx model for Datasets I and II, respectively.
Table 12. Confidence intervals for the parameters of the NPTLILx model for Datasets I and II, respectively.
CI α β θ
90 % [0 0.2157][8.0147 15.9656][0.5081 2.6220]
95 % [0 0.2440][7.2534 16.7269][0.3058 2.8244]
CI α β θ
90 % [12.7289 16.1432][0 2.2642][4.6302 5.4300]
95 % [12.4020 16.4701][0 2.6139][4.5536 5.5065]

Share and Cite

MDPI and ACS Style

Bantan, R.A.R.; Jamal, F.; Chesneau, C.; Elgarhy, M. A New Power Topp–Leone Generated Family of Distributions with Applications. Entropy 2019, 21, 1177. https://doi.org/10.3390/e21121177

AMA Style

Bantan RAR, Jamal F, Chesneau C, Elgarhy M. A New Power Topp–Leone Generated Family of Distributions with Applications. Entropy. 2019; 21(12):1177. https://doi.org/10.3390/e21121177

Chicago/Turabian Style

Bantan, Rashad A. R., Farrukh Jamal, Christophe Chesneau, and Mohammed Elgarhy. 2019. "A New Power Topp–Leone Generated Family of Distributions with Applications" Entropy 21, no. 12: 1177. https://doi.org/10.3390/e21121177

APA Style

Bantan, R. A. R., Jamal, F., Chesneau, C., & Elgarhy, M. (2019). A New Power Topp–Leone Generated Family of Distributions with Applications. Entropy, 21(12), 1177. https://doi.org/10.3390/e21121177

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