Next Article in Journal
A Quantitative Approach to Evaluate the Application of the Extended Situational Teaching Model in Engineering Education
Previous Article in Journal
General Formulas for the Central and Non-Central Moments of the Multinomial Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kumaraswamy Generalized Power Lomax Distributionand Its Applications

by
Vasili B.V. Nagarjuna
1,
R. Vishnu Vardhan
1 and
Christophe Chesneau
2,*
1
Department of Statistics, Pondicherry University, Pondicherry 605 014, India
2
Department of Mathematics, LMNO, Université de Caen, Campus II, Science 3, 14032 Caen, France
*
Author to whom correspondence should be addressed.
Submission received: 24 December 2020 / Revised: 2 January 2021 / Accepted: 2 January 2021 / Published: 7 January 2021

Abstract

:
In this paper, a new five-parameter distribution is proposed using the functionalities of the Kumaraswamy generalized family of distributions and the features of the power Lomax distribution. It is named as Kumaraswamy generalized power Lomax distribution. In a first approach, we derive its main probability and reliability functions, with a visualization of its modeling behavior by considering different parameter combinations. As prime quality, the corresponding hazard rate function is very flexible; it possesses decreasing, increasing and inverted (upside-down) bathtub shapes. Also, decreasing-increasing-decreasing shapes are nicely observed. Some important characteristics of the Kumaraswamy generalized power Lomax distribution are derived, including moments, entropy measures and order statistics. The second approach is statistical. The maximum likelihood estimates of the parameters are described and a brief simulation study shows their effectiveness. Two real data sets are taken to show how the proposed distribution can be applied concretely; parameter estimates are obtained and fitting comparisons are performed with other well-established Lomax based distributions. The Kumaraswamy generalized power Lomax distribution turns out to be best by capturing fine details in the structure of the data considered.

1. Introduction

For several decades, researchers have been working to come up with several new distributions to meet certain practical requirements. The motivation is that, in concrete applications related to disciplines such as hydrology, econometrics and many others, the standard distributions have been observed to lack fit. For instance, for daily precipitation and daily vapor flow data, Kumaraswamy [1] showed that the beta distribution does not provide a suitable fit. Continuing this work, in References [2,3], it is also explored that distributions like Johnson, sinepower and extended sinepower distributions were satisfactory in fitting the above mentioned data type. However, in later years, an alternative distribution of finite range was suggested by Kumaraswamy [1] which is named later as Kw distribution. It reveals to fit such data appropriately. The cumulative distribution function (cdf) and probability density function (pdf) are given in (1) and (2), respectively:
F ( x ; a , b ) = 1 ( 1 x a ) b , x ( 0 , 1 ) ,
and
f ( x ; a , b ) = a b x a 1 ( 1 x a ) b 1 ,
where a > 0 and b > 0 are shape parameters, with the usual modifications for x ( 0 , 1 ) . The main advantage of the Kw distribution is that it has the shape parameter that addresses data that has an extended tail nature. Also, the Kw distribution has been well received by many researchers for fitting skewed types of data sets from hydrology and other engineering disciplines. Other works on Kw distribution has been planned by several researchers, pointing out that it is a special case of the three-parameter beta distribution [4]. Also, the similarities along with basic properties have been extensively studied by Jones [5]. In the recent past, a generalized version of the Kw distribution has been proposed by Cordeiro and de Castro [6] with the cdf and pdf given in (3) and (4), respectively:
F G ( x ; a , b ) = 1 [ 1 G ( x ) a ] b , x R ,
and
f G ( x ; a , b ) = a b g ( x ) G ( x ) a 1 [ 1 G ( x ) a ] b 1 .
Here, G ( x ) is the cdf of a chosen base distribution, with pdf given as g ( x ) . Clearly, if a = b = 1 , the forms in (3) and (4) reduce to the pdf and cdf from the base distribution. This generalized version of the Kw distribution is called the Kw-G family of distributions. Extensive work on the Kw-G family of distributions has been observed exponentially by proposing new distributions of various asymmetric natures. To cite a few, there are the Kumaraswamy-Weibull distribution [7], Kumaraswamy-Gumbel distribution [8], Kumaraswamy generalized gamma distribution [9] and Kumaraswamy-Burr XII (KBXII) distribution [10]. Applications of the Kw-G family can also be found in References [11,12,13], among others.
In general, the goal of providing new distributions is to create flexible mathematical models capable of handling non-normal data scenarios. This flexibility can be achieved in a simple way by adding additional parameters such as location, scale and shape. In similar lines of thought process, several distributions extending the famous Lomax distribution, such as the exponentiated-Lomax (EL) distribution [14], extended Lomax distribution [15], Kumaraswamy-generalized Lomax distribution [16], exponential-Lomax distribution [17], Weibull-Lomax (WL) distribution [18], Weibull Fréchet (WFr) distribution [19], power Lomax (PL) distribution [20], half-logistic Lomax distribution [21], inverse PL distribution [22], Topp-Leone Lomax (TLGL) distribution [23], type II Topp-Leone power Lomax (TIITLPL) distribution [24], Marshall-Olkin exponential Lomax distribution [25] and Marshall-Olkin length biased Lomax distribution [26] were proposed and developed.
In this present work, an attempt to propose a new distribution by compounding the PL distribution into the general Kw-G family of distributions is made. The proposed distribution is called the Kumaraswamy generalized PL distribution, KPL for short. The rationality of considering the PL distribution is that it equips the most famous extensions of the Lomax distribution [20] and it allows applications dealing with heavy tailed data. Because of its nature, we wish to exhibit mathematical flexibility by adding two additional shape parameters that are presented in the Kw-G family of distributions. To be more precise, a short retrospective on the PL distribution is necessary. First, let us mention that the cdf and pdf of the PL distribution have the forms in (5) and (6), respectively:
G ( x ; α , β , λ ) = 1 λ λ + x β α , x > 0
and
g ( x ; α , β , λ ) = α β x β 1 λ λ λ + x β α + 1 ,
where α > 0 is a shape parameter, and β > 0 and λ > 0 are scale parameters. Rady et al. [20] mentioned that the hazard rate function (hrf) of the PL distribution does not have an increasing curve, which remains a serious limitation for some modeling purposes. This issue is addressed in this work by making the use of the shape parameter and it is performed in the proposed KPL distribution. In particular, the various forms of the pdf of the KPL distribution show that, with increasing values of the new parameters that will be denoted “a” and “b”, it is unimodal and can attain the symmetric nature curve. Also, the corresponding hrf possesses decreasing, increasing and inverted (upside-down) bathtub shapes. In addition, decreasing-increasing-decreasing shapes are observed, which is a clear plus for various statistical purposes. The KPL distribution and its statistical properties such as quantiles, moments, information measures, order statistics and maximum likelihood (ML) estimation are detailed out in subsequent sections of the article. Using two famous data sets, namely turbo charger data set and radiation therapy data set, we demonstrate that the proposed distribution is better suited compared to different types of Lomax distribution, including the PL distribution.
The rest of the article is structured by the following sections. Section 2 completes the presentation of the KPL distribution. Section 3 is devoted to its moments analysis. Section 4 is about information measures of the KPL distribution. Section 5 discusses the related order statistics. The estimation of the parameters of the KPL distribution is afforded in Section 6, including a simulation study. Section 7 focuses on the applications by considering two practical data sets. A summary is given in Section 8.

2. The Kumaraswamy Generalized Power Lomax Distribution

By considering (5) and (6) in (3) and (4), we obtain the following cdf and pdf, respectively,
F K P L ( x ; ξ ) = 1 1 1 λ λ + x β α a b , x > 0 ,
where ξ = ( α , β , λ , a , b ) ( 0 , ) 5 , and
f K P L ( x ; ξ ) = a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 1 1 λ λ + x β α a b 1 .
The expressions (7) and (8) constitute the cdf and pdf of the KPL distribution, respectively. The KPL distribution contains several existing distribution, including the PL distribution for a = b = 1 and the TIITLPL distribution for a = 2 . It can be also viewed as a re-parametrized version of the KBXII distribution. As already mentioned, the roles of a and b will be major in the interests of the KPL distribution, reaching new levels of flexibility compared to those of the PL distribution, among others.
As preliminary properties, note that F K P L ( x ; ξ ) is a decreasing function with respect to a and λ , an increasing function with respect to b and α , and a non-monotonic function with respect to β . This implies various first-order stochastic dominance properties. For instance, for ξ 1 = ( α 1 , β , λ , a , b 2 ) and ξ 2 = ( α 2 , β , λ , a , b 2 ) with α 1 α 2 and b 1 b 2 , since F K P L ( x ; ξ ) is an increasing function with respect to b and α , we have F K P L ( x ; ξ 1 ) F K P L ( x ; ξ 2 ) . Similar inequalities can be presented by taking into account the other parameters. Also, it is worth mentioning that the KPL distribution is heavy-tailed; for all t > 0 , we can prove that 0 e t x f K P L ( x ; ξ ) d x = .
Considering different values of the parameters, variant density forms of the KPL distribution are obtained, and are shown in Figure 1.
From Figure 1, we observe that the pdf of the KPL distribution can be decreasing or unimodal, with very flexible skewness, peakness and plateness. One can show that the decreasing case corresponds to a β < 1 . When the pdf is unimodal, we see that it is mainly ‘almost symmetric’ or ‘right-skewed’, which is ideal for the modelling of diverse lifetime phenomena.
Since the KPL distribution belongs to the family of lifetime distributions, its hrf is of interest to deals with some of the statistical properties of the proposed distribution. These properties will help out to exhibit the practical applications and characterizations of real data phenomenon. The hrf of the KPL distribution is given by
h K P L ( x ; ξ ) = f K P L ( x ; ξ ) 1 F K P L ( x ; ξ ) ,
that is, by substituting the Equations (7) and (8) in h K P L ( x ; ξ ) ,
h K P L ( x ; ξ ) = a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 1 1 λ λ + x β α a 1 , x > 0 .
Considering various values of the parameters, the hrf of the KPL distribution contains different kinds of shapes as shown in Figure 2.
From Figure 2, we highlight a crucial difference between the PL and KPL distributions. Indeed, an immediate limitation of the PL distribution is that its hrf cannot be increasing (see Reference [20]). This limitation is overcome in this work by using shape parameters and it is shown that the hrf of the proposed KPL distribution may be increasing. At the same time, we are able to present a decreasing-increasing-decreasing hrf. It is another point to emphasize that the proposed distribution has a better way of expressing different natures of data.
We end this part by presenting the quantile function (qf) of the KPL distribution which also defines it in the mathematical sense. This qf is rigorously defined as the inverse function of F K P L ( x ; ξ ) . By solving the following nonlinear equation: F K P L [ Q K P L ( u ; ξ ) ; ξ ] = u with u ( 0 , 1 ) and respect to Q K P L ( u ; ξ ) , we obtain
Q K P L ( u ; ξ ) = λ 1 / β 1 1 ( 1 u ) 1 / b 1 / a 1 / α 1 1 / β , u ( 0 , 1 ) .
We immediately derive the median of the KPL distribution as
M = λ 1 / β 1 1 ( 0.5 ) 1 / b 1 / a 1 / α 1 1 / β .
Additional quantile analysis can be performed based on Q K P L ( u ; ξ ) . In this regard, one may refer to Reference [27].

3. Moments of the KPL Distribution

By definition, for any positive integer r, the r t h order moment about the origin of a random variable X following the KPL distribution is given as
μ r = E ( X r ) = 0 x r f K P L ( x ; ξ ) d x ,
where E denotes the expectation. Clearly, in view of the mathematical complexity of f K P L ( x ; ξ ) , a simple expression of this integral is not possible. Let us first study its existence according to the values of the parameters of the KPL distribution. At the neighborhood of x = 0 , we have f K P L ( x ; ξ ) ( a b α a β / λ a ) x β a 1 and, by the Riemann integrability criterion, since r + β a > 0 , the function x r + β a 1 is integrable over ( 0 , ϵ ) with ϵ > 0 . Now, at the neighborhood of x = , we have f K P L ( x ; ξ ) a b b α β λ α b x β α b 1 and, by the Riemann integrability criterion, the function x r β α b 1 is integrable over ( ϵ , ) with ϵ > 0 if and only if r < β α b . In summary, μ r exists if and only if r < β α b . In this case, one can approximate it via various numerical procedures.
For an analytical approach, one can derive series expansion for f K P L ( x ; ξ ) and plug-in into (9). With this in mind, the next result presents a series expansion of a power transformation of the pdf of the KPL distribution.
Proposition 1.
For any ρ > 0 , we have
f K P L ( x ; ξ ) ρ = i = 0 j = 0 υ i , j ( ξ , ρ ) d j ( x ; α , β , λ , ρ ) , x > 0 ,
where, by introducing generalized binomial coefficients,
υ i , j ( ξ , ρ ) = a b α β λ ρ ρ ( b 1 ) i a i + ρ ( a 1 ) j ( 1 ) i + j
and
d j ( x ; α , β , λ , ρ ) = x ρ ( β 1 ) λ λ + x β α j + ρ ( α + 1 ) .
Proof. 
We have
f K P L ( x ; ξ ) ρ = a b α β λ ρ x ρ ( β 1 ) λ λ + x β ρ ( α + 1 ) 1 λ λ + x β α ρ ( a 1 ) × 1 1 λ λ + x β α a ρ ( b 1 ) .
Now, the generalized binomial theorem states that ( 1 Z ) b = i = 0 b i ( 1 ) i Z i for any real numbers b and Z such that | Z | < 1 , condition that can be removed if b is a positive integer. Therefore, by the application of this theorem two times in a row, we obtain
f K P L ( x ; ξ ) ρ = a b α β λ ρ x ρ ( β 1 ) λ λ + x β ρ ( α + 1 ) i = 0 ρ ( b 1 ) i ( 1 ) i 1 λ λ + x β α a i + ρ ( a 1 ) = a b α β λ ρ x ρ ( β 1 ) i = 0 j = 0 ρ ( b 1 ) i a i + ρ ( a 1 ) j ( 1 ) i + j λ λ + x β α j + ρ ( α + 1 ) = i = 0 j = 0 υ i , j ( ξ , ρ ) d ( x ; α , β , λ , ρ ) .
The desired expansion is obtained. □
The particular case ρ = 1 in Proposition 1 gives
f K P L ( x ; ξ ) = i = 0 j = 0 υ i , j ( ξ , 1 ) d j ( x ; α , β , λ , 1 ) .
Hence, under the condition that r < β α inf ( b , 1 ) permitting to interchange the integral and sum signs, the r t h order moment about the origin of X is given as
μ r = 0 x r i = 0 j = 0 υ i , j ( ξ , 1 ) d j ( x ; α , β , λ , 1 ) d x = i = 0 j = 0 υ i , j ( ξ , 1 ) 0 x r d j ( x ; α , β , λ , 1 ) d x .
Let us now discuss a tractable expression for the integral term. We have
0 x r d j ( x ; α , β , λ , 1 ) d x = 0 x r + β 1 λ λ + x β α ( j + 1 ) + 1 d x .
Let us now set u = x β / λ , so x = ( u λ ) 1 / β and d u = ( β / λ ) x β 1 d x with no change at the boundaries, implying that
0 x r + β 1 λ λ + x β α ( j + 1 ) + 1 d x = λ β 0 ( u λ ) r / β ( 1 + u ) α ( j + 1 ) + 1 d u = λ r / β + 1 β 0 u ( r / β + 1 ) 1 ( 1 + u ) ( r / β + 1 ) + [ α ( j + 1 ) r / β ] d u = λ r / β + 1 β B r β + 1 , α ( j + 1 ) r β ,
where B ( a , b ) denotes the standard beta function defined by B ( a , b ) = 0 1 t a 1 ( 1 t ) b 1 d t with a > 0 and b > 0 . We finally get
μ r = i = 0 j = 0 υ i , j ( ξ , 1 ) λ r / β + 1 β B r β + 1 , α ( j + 1 ) r β .
This formula is exact, without approximation. It can serve to determine the exact numerical values of μ r and all the associated measures. The most basic of them are the mean of X defined by μ 1 and the variance given by V a r = μ 2 ( μ 1 ) 2 . A practical approximation of μ r is given by
μ r = i = 0 50 j = 0 50 υ i , j ( ξ , 1 ) λ r / β + 1 β B r β + 1 , α ( j + 1 ) r β ,
the bound 50 being an integer chosen arbitrary large.
As an illustration, with the use of the R software, Table 1 provides the moments about the origin and variance of X with different parameter values of the KPL distribution.
From Table 1, we see that the fourth moments about the origin and variance have notable numerical variations. This is particularly obvious for the last combination of parameters: b = 2.7 , α = 0.8 , λ = 0.5 , β = 2.2 with a = 4 , 11 , 17 , 22 and 36. The high numerical variability of these moment measures also testifies to the flexibility of the KPL distribution.
The incomplete versions of the moments about the origin can have a similar mathematical treatment. They allow us to define various deviation measures and diverse types of residual life function, as those managed in References [10,20,24].

4. Information Measures

In this section, some information measures of the KPL distribution are discussed, namely the Rényi entropy and β -entropy measures. Both measuring the variation or uncertainty of the considered distribution.

4.1. Rényi Entropy

Rényi [28] provided an useful extension of the Shannon entropy. The Rényi entropy of the KPL distribution can be defined as
I R ( ρ ) = 1 1 ρ log 0 f K P L ( x ; ξ ) ρ d x ,
with ρ > 0 and ρ 1 . Let us now study the existence of this entropy measure which depends on the existence of its integral term. At the neighborhood of x = 0 , we have f K P L ( x ; ξ ) ρ ( a b α a β / λ a ) ρ x ρ ( β a 1 ) and, by the Riemann integrability criterion, the function x ρ ( β a 1 ) is integrable over ( 0 , ϵ ) with ϵ > 0 if and only if ρ ( β a 1 ) > 1 . Now, at the neighborhood of x = , we have f K P L ( x ; ξ ) ρ [ a b b α β λ α b ] ρ x ρ ( β α b + 1 ) and, by the Riemann integrability criterion, the function x ρ ( β α b + 1 ) is integrable over ( ϵ , ) with ϵ > 0 if and only if ρ ( β α b + 1 ) > 1 . Hence, I R ( ρ ) exists if and only if ρ ( β a 1 ) > 1 and ρ ( β α b + 1 ) > 1 . If these conditions are satisfied, one can approximate I R ( ρ ) through the approximation of its integral term via various numerical procedures.
In a more analytical manner, one can use Proposition 1. Under the conditions above plus ρ ( β 1 ) > 1 and ρ ( β α + 1 ) > 1 , a direct application of this result gives
0 f K P L ( x ; ξ ) ρ d x = 0 i = 0 j = 0 υ i , j ( ξ , ρ ) d j ( x ; α , β , λ , ρ ) d x = i = 0 j = 0 υ i , j ( ξ , ρ ) 0 d j ( x ; α , β , λ , ρ ) d x .
A comprehensive expression of the integral term is developed below. By setting u = x β / λ , so x = ( u λ ) 1 / β and d u = ( β / λ ) x β 1 d x with no change at the boundaries, it comes
0 d j ( x ; α , β , λ , ρ ) d x = 0 x ρ ( β 1 ) λ λ + x β α j + ρ ( α + 1 ) d x = λ β 0 ( u λ ) ( ρ 1 ) ( 1 1 / β ) ( 1 + u ) α j + ρ ( α + 1 ) d u = λ ( ρ 1 ) ( 1 1 / β ) + 1 β 0 u [ ( ρ 1 ) ( 1 1 / β ) + 1 ] 1 ( 1 + u ) α j + ρ ( α + 1 ) [ ( ρ 1 ) ( 1 1 / β ) + 1 ] + [ ( ρ 1 ) ( 1 1 / β ) + 1 ] d u = λ ( ρ 1 ) ( 1 1 / β ) + 1 β B ( ρ 1 ) 1 1 β + 1 , α j + ρ ( α + 1 ) ( ρ 1 ) 1 1 β 1 .
Therefore, we obtain an expression for I R ( ρ ) as
I R ( ρ ) = 1 1 ρ log { i = 0 j = 0 υ i , j ( ξ , ρ ) λ ( ρ 1 ) ( 1 1 / β ) + 1 β × B ( ρ 1 ) 1 1 β + 1 , α j + ρ ( α + 1 ) ( ρ 1 ) 1 1 β 1 } .
This formula is exact, without approximation. The following simple approximation can be derived for practical purposes:
I R ( ρ ) 1 1 ρ log { i = 0 50 j = 0 50 υ i , j ( ξ , ρ ) λ ( ρ 1 ) ( 1 1 / β ) + 1 β × B ( ρ 1 ) 1 1 β + 1 , α j + ρ ( α + 1 ) ( ρ 1 ) 1 1 β 1 } .
Here again, the bound 50 must be viewed as a large integer arbitrarily chosen.

4.2. Tsallis Entropy

The Tsallis entropy or q-entropy was discovered by Havrda and Charvat [29]. Later, it was developed by Tsallis [30] in the context of Physics. The Tsallis entropy of the KPL distribution can be defined as
I H ( q ) = 1 q 1 1 0 f K P L ( x ; ξ ) q d x ,
where q > 0 and q 1 . Then, based on the previous work on the Rényi entropy, I H ( q ) exists if and only if q ( β a 1 ) > 1 and q ( β α b + 1 ) > 1 . In all cases, we can approximate it via numerical procedures. With the following additional assumptions: q ( β 1 ) > 1 and q ( β α + 1 ) > 1 , proceeding as for I R ( ρ ) , we can expand I H ( q ) as
I H ( q ) = 1 q 1 { 1 i = 0 j = 0 υ i , j ( ξ , q ) λ ( q 1 ) ( 1 1 / β ) + 1 β × B ( q 1 ) 1 1 β + 1 , α j + q ( α + 1 ) ( q 1 ) 1 1 β 1 ] .
Based on this formula, analytical approximation can be conducted.

5. Order Statistics

The modeling of certain random systems requires the concept of order statistic. Basically, for r = 1 , , n , the r t h order statistic of a statistical sample is equal to its r t h smallest value. In what follows, some immediate distributional properties of the order statistics of the KPL distribution are presented.
Let X ( 1 ) , X ( 2 ) , , X ( n ) be an ordered random sample distributed with the KPL distribution. Then, the pdf of X ( r ) is computed as
f r ( x ; ξ ) = 1 B ( r , n r + 1 ) F K P L ( x ; ξ ) r 1 [ 1 F K P L ( x ; ξ ) ] n r f K P L ( x ; ξ ) ,
that is, by substituting the Equations (7) and (8) in f r ( x ; ξ ) ,
f r ( x ; ξ ) = 1 B ( r , n r + 1 ) a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 × 1 1 λ λ + x β α a b ( n r + 1 ) 1 1 1 1 λ λ + x β α a b r 1 , x > 0 .
For r = 1 , we get the pdf of the first order statistics X ( 1 ) = min X 1 , X 2 , , X n as follows
f 1 ( x ; ξ ) = n a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 1 1 λ λ + x β α a b n 1 , x > 0 .
One can remark that f 1 ( x ; ξ ) = f K P L ( x ; α , β , λ , a , b n ) , meaning that the distribution of X ( 1 ) is also a KPL distribution.
Similarly, for r = n , we get the pdf of the n t h order statistics X ( n ) = max X 1 , X 2 , , X n as follows
f n ( x ; ξ ) = n a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 × 1 1 λ λ + x β α a b 1 1 1 1 λ λ + x β α a b n 1 , x > 0 .
The next result exhibits the linear relation existing between the pdf of X ( r ) and some pdfs of the KPL distribution.
Proposition 2.
The pdf of X ( r ) can be expressed as a linear combination of pdfs of the KPL distribution, and, more precisely,
f r ( x ; ξ ) = 1 B ( r , n r + 1 ) k = 0 r 1 r 1 k ( 1 ) k 1 n r + k + 1 f K P L ( x ; α , β , λ , a , b ( n r + k + 1 ) ) , x > 0 .
Proof. 
It follows from (10) and the (standard) binomial theorem that
f r ( x ; ξ ) = 1 B ( r , n r + 1 ) a b α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 × 1 1 λ λ + x β α a b ( n r + 1 ) 1 k = 0 r 1 r 1 k ( 1 ) k 1 1 λ λ + x β α a b k = 1 B ( r , n r + 1 ) k = 0 r 1 r 1 k ( 1 ) k 1 n r + k + 1 × a b ( n r + k + 1 ) α β λ x β 1 λ λ + x β α + 1 1 λ λ + x β α a 1 × 1 1 λ λ + x β α a b ( n r + k + 1 ) 1 = 1 B ( r , n r + 1 ) k = 0 r 1 r 1 k ( 1 ) k 1 n r + k + 1 f K P L ( x ; α , β , λ , a , b ( n r + k + 1 ) ) .
This ends the proof of Proposition 2. □
An immediate consequence of Proposition 2 is the determination of some properties for X ( r ) based on those of the KPL distribution. For example, the s t h order moment of X ( r ) about the origin can be written as
μ r , s = E ( X ( r ) s ) = 1 B ( r , n r + 1 ) k = 0 r 1 r 1 k ( 1 ) k 1 n r + k + 1 μ s ( k ) ,
where μ s ( k ) denotes the s t h order moment about the origin of a random variable with the KPL distribution with parameters α , β , λ , a and b ( n r + k + 1 ) .

6. Maximum Likelihood Estimates of the Parameters

The ML estimation method is used for estimating the unknown parameters of the distribution. Let x = ( x 1 , x 2 , x n ) be a random sample drawn from the KPL distribution. Then the likelihood function and log-likelihood function corresponding to the Equation (8) are, respectively, as follows
L ( x ; ξ ) = a b α β λ n i = 1 n x i β 1 λ λ + x i β α + 1 1 λ λ + x i β α a 1 1 1 λ λ + x i β α a b 1
and
log L ( x ; ξ ) = n log a b α β λ + ( β 1 ) i = 1 n log x i + ( α + 1 ) i = 1 n log λ λ + x i β + ( a 1 ) i = 1 n log ϑ i + ( b 1 ) i = 1 n log ( 1 ϑ i a ) ,
where it is set ϑ i = 1 λ / ( λ + x i β ) α . The ML estimates (MLEs) of the parameters a , b , α , β and λ , say a ^ , b ^ , α ^ , β ^ and λ ^ , are given by ξ ^ = ( a ^ , b ^ , α ^ , β ^ , λ ^ ) making L ( x ; ξ ^ ) or log L ( x ; ξ ^ ) maximal. These MLEs can be obtained by solving the following nonlinear equations:
a log L ( x ; ξ ) = n a + i = 1 n log ϑ i + ( 1 b ) i = 1 n log ϑ i ( ϑ i a 1 ) = 0 ,
b log L ( x ; ξ ) = n b + i = 1 n log ( 1 ϑ i a ) = 0 ,
α log L ( x ; ξ ) = n α + i = 1 n log λ λ + x i β + i = 1 n λ λ + x i β α log λ λ + x i β a ( b 1 ) ϑ i ( ϑ i a 1 ) ( a 1 ) ϑ i = 0 ,
β log L ( x ; ξ ) = n β + i = 1 n log x i ( α + 1 ) i = 1 n x i β λ + x i β λ λ + x i β α 1 λ x i β log x i ( λ + x i β ) 2 × α ( a 1 ) ϑ i a α ( b 1 ) ϑ i ( ϑ i a 1 ) log x i = 0
and
λ log L ( x ; ξ ) = n λ + α + 1 λ i = 1 n x i β λ + x i β + i = 1 n x i β λ λ + x i β α + 1 a α ( b 1 ) λ 2 ϑ i ( ϑ i a 1 ) α ( a 1 ) λ 2 ϑ i = 0 .
Based on data, the MLEs can be obtained numerically by the iterative procedure of Newton-Raphson method for a system of simultaneous nonlinear equations. As an example of use, Monte Carlo simulations are carried out to assess the finite sample behavior of the MLEs a ^ , b ^ , α ^ , β ^ and λ ^ . For a given sample size, 1000 random samples drawn from the KPL distribution with given parameters are generated by using the qf technique. In this setting, the MLEs of the five model parameters along with the respective bias and mean square error (MSE) for the sample sizes n = { 50 , 100 , 250 } are shown in Table 2.
As mentioned initially about the advantage of having additional shape parameters, the same is true between bias and MSE. From the results of Table 2, it is clearly evident that the estimates are quite stable and close to the true values of the parameters for these sample sizes. Additionally, as the sample size increases, the biases and MSEs of the MLEs decrease as expected.

7. Applications of the KPL Model

Two real data sets are used as applications of the proposed KPL distribution as heavy tailed distribution.
Data set 1: (Strength data [31]) The data represent 69 strength data for single carbon fibers (and impregnated 1000-carbon fiber tows). The measures in GPA by subtracting 1 are: 0.0312, 0.314, 0.479, 0.552, 0.700, 0.803, 0.861, 0.865, 0.944, 0.958, 0.966, 0.977, 1.006, 1.021, 1.027, 1.055, 1.063, 1.098, 1.140, 1.179, 1.224, 1.240, 1.253, 1.270, 1.272, 1.274, 1.301, 1.301, 1.359, 1.382, 1.382, 1.426, 1.434, 1.435, 1.478, 1.490, 1.511, 1.514, 1.535, 1.554, 1.566, 1.570, 1.586, 1.629, 1.633, 1.642, 1.648, 1.684, 1.697, 1.726, 1.770, 1.773, 1.800, 1.809, 1.818, 1.821, 1.848, 1.880, 1.954, 2.012, 2.067, 2.084, 2.090, 2.096, 2.128, 2.233, 2.433, 2.585, 2.585.
Data set 2: (Theft data [32]) The data represent the amounts of 120 theft claims made in a home insurance portfolio. The 120 theft claims data are: 3, 11, 27, 36, 47, 49, 54, 77, 78, 85, 104, 121, 130, 138, 139, 140, 143, 153, 193, 195, 205, 207, 216, 224, 233, 237, 254, 257, 259, 265, 273, 275, 278, 281, 396, 405, 412, 423, 436, 456, 473, 475, 503, 510, 534, 565, 656, 656, 716, 734, 743, 756, 784, 786, 819, 826, 841, 842, 853, 860, 877, 942, 942, 945, 998, 1029, 1066, 1101, 1128, 1167, 1194, 1209, 1223, 1283, 1288,1296, 1310, 1320, 1367, 1369, 1373, 1382, 1383, 1395, 1436, 1470, 1512, 1607, 1699, 1720, 1772, 1780, 1858, 1922, 2042, 2247, 2348, 2377, 2418, 2795, 2964, 3156, 3858, 3872, 4084, 4620, 4901, 5021, 5331, 5771, 6240, 6385, 7089, 7482, 8059, 8079, 8316, 11,453, 22,274, 32,043. After analyzing the histograms, data set 1 shows an offset deviation from the symmetrical pattern while data set 2 shows a decreasing histogram shape.
In order to show the best fit of the KPL distribution, some other distributions based on the Lomax distribution are considered and used for comparison. These competing distributions have already been mentioned in the introduction, and are the KBXII distribution [10], PL distribution [20], WL distribution [18], WFr distribution [19], TLGL distribution [23], EL distribution [14] and the basic Lomax distribution.
The ML estimation method is applied for all the distribution parameters. The MLEs are obtained by iterative procedures. The MLEs of the distribution parameters are given in Table 3.
Based on the notations of the KPL distribution, we now present the measures of adequacy that we use. Let x 1 , , x n represent the data and x ( 1 ) , , x ( n ) be their ordered values. First, we consider the Cramér-von Mises (W*), Anderson Darling (A*) and Kolmogorov-Smirnov (K-S) statistics ( D n ) defined by
W * = 1 12 n + i = 1 n F ( x ( i ) ; ξ ^ ) 2 i 1 n 2 ,
A * = n i = 1 n 2 i 1 n log ( F ( x ( i ) ; ξ ^ ) ) + log ( 1 F ( x i ; ξ ^ ) )
and
D n = max i = 1 , , n i n F ( x ( i ) ; ξ ^ ) , F ( x ( i ) ; ξ ^ ) i 1 n ,
respectively, where n is the sample size, ξ denotes the parameters of the distribution ( ξ = ( a , b , α , β , λ ) for the KPL distribution) and ξ ^ its vectorial MLE. The p-Value of the K-S test related to D n is also considered. The adequacy measures are widely used to know which distribution suits in a better manner. The distribution having the minimum value for the W* or A*, and maximum value for the p-Value, is chosen as the best one that is in adequacy to the data.
Also, we consider the Akaike information criterion (AIC), correct Akaike information criterion (CAIC), Bayesian information criterion (BIC) and Hannan-Quinn information criterion (HQIC), defined as
AIC = 2 log L ( x ; ξ ^ ) + 2 k , BIC = 2 log L ( x ; ξ ^ ) + k log ( n ) , CAIC = 2 log L ( x ; ξ ^ ) + 2 k n n k 1 , HQIC = 2 log L ( x ; ξ ^ ) + 2 k log [ log ( n ) ] ,
respectively, where k is the number of parameters ( k = 5 for the KPL distribution). As commonly accepted, the distribution having the minimum value for the AIC or CAIC or BIC or HQIC value is chosen as the best one that fits the data. For the considered data sets and distributions, the values of the measures above are computed and reported in Table 4.
From the results of Table 4, it is evident that best fit is observed with the proposed KPL distribution and other distributions based on the Lomax distribution attained worst information criterion values. This is also witnessed through the fits of the pdfs that are depicted in Figure 3.
Among all the listed distributions, the KPL distribution has a better fit for the data considered. We confirm this visual result by plotting the Probability-Probability (PP) plots of the estimated distributions in Figure 4 and Figure 5 for data sets 1 and 2, respectively.
From Figure 4 and Figure 5, for the two data sets, it is clear that the best fit of the black diagonal line is assigned to the red line, corresponding to the one of the KPL distribution.
From these analyses, it can be seen that having additional shape parameters in the distribution is an advantage when talking about the extended tail version of the data. In addition, the dispersion and trend of the data can be characterized using the scale and location parameters of the basic distribution. Overall, a distribution that is compared to measurements of location, scale, and shape parameters always has a better advantage in handling non-normal data structures. Here, this advantage is captured by having two additional shape parameters which helped to see the characterization and fit in the best way.

8. Summary

The work carried out in this paper is to address some limitations of the PL distribution and also to illustrate the usefulness of the shape parameters in handling non-normal data. An attempt is made to introduce a new distribution, namely the KPL distribution, which is obtained by inducing the cdf of the PL distribution into the functional form of the Kw-G family of distributions. It contains five parameters consisting of one location, one scale and three shape parameters. In the work of Rady et al. [20], it is learned that the hrf of the PL distribution fails to attain the increasing pattern and decreasing-increasing-increasing pattern. This is not the case for the hrf of the KPL distribution. Important properties are studied, including moments, entropy measures and order statistics. With the help of two practical data sets, the fit of the KPL distribution is done. Information criterion measures are compared between the KPL distribution and some distributions also based on the Lomax distribution, including the PL distribution. It is shown that the KPL distribution has a better fit. Hence, the KPL distribution can work in a better manner for some kind of non-normal data.

Author Contributions

V.B.V.N., R.V.V. and C.C. have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to the four anonymous referees for a careful checking of the details and for helpful comments that improved this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumaraswamy, P. A generalized probability density function for double-bounded random processes. J. Hydrol. 1980, 46, 79–88. [Google Scholar] [CrossRef]
  2. Kumaraswamy, P. Sinepower probability density function. J. Hydrol. 1976, 31, 181–184. [Google Scholar] [CrossRef]
  3. Kumaraswamy, P. Extended sinepower probability density function. J. Hydrol. 1978, 37, 81–89. [Google Scholar] [CrossRef]
  4. Nadarajah, S. On the distribution of Kumaraswamy. J. Hydrol. 2008, 348, 568–569. [Google Scholar] [CrossRef]
  5. Jones, M.C. Kumaraswamy’s distribution: A beta-type distribution with some tractability advantages. Stat. Methodol. 2009, 6, 70–81. [Google Scholar] [CrossRef]
  6. Cordeiro, G.M.; de Castro, M. A new family of generalized distributions. J. Stat. Comput. Simul. 2011, 81, 883–898. [Google Scholar] [CrossRef]
  7. 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]
  8. Cordeiro, G.M.; Nadarajah, S.; Ortega, E.M. The Kumaraswamy Gumbel distribution. Stat. Methods Appl. 2012, 21, 139–168. [Google Scholar] [CrossRef]
  9. De Pascoa, M.A.; Ortega, E.M.; Cordeiro, G.M. The Kumaraswamy generalized gamma distribution with application in survival analysis. Stat. Methodol. 2011, 8, 411–433. [Google Scholar] [CrossRef]
  10. Paranaíba, P.F.; Ortega, E.M.; Cordeiro, G.M.; Pascoa, M.A.d. The Kumaraswamy Burr XII distribution: Theory and practice. J. Stat. Comput. Simul. 2013, 83, 2117–2143. [Google Scholar] [CrossRef]
  11. Malinova, A.; Golev, A.; Rahneva, O.; Kyurkchiev, V. Some notes on the Kumaraswamy-Weibull-Exponential cumulative sigmoid. Int. J. Pure Appl. Math. 2018, 120, 521–529. [Google Scholar]
  12. Malinova, A.; Kyurkchiev, V.; Iliev, A.; Kyurkchiev, N. A note on the transmuted Kumaraswamy quasi Lindley cumulative distribution function. Int. J. Sci. Res. Dev. 2018, 6, 561–564. [Google Scholar]
  13. Angelova, E.; Arnaudova, V.; Terzieva, T.; Malinova, A. A note on the new Kumaraswamy alpha power inverted exponential family of C.D.F. Neural Parallel Sci. Comput. 2020, 28, 59–67. [Google Scholar]
  14. Abdul-Moniem, I.B. Recurrence relations for moments of lower generalized order statistics from exponentiated Lomax distribution and its characterization. J. Math. Comput. Sci. 2012, 2, 999–1011. [Google Scholar]
  15. Lemonte, A.J.; Cordeiro, G.M. An extended Lomax distribution. Statistics 2013, 47, 800–816. [Google Scholar] [CrossRef]
  16. Shams, T.M. The Kumaraswamy-Generalized Lomax Distribution. Middle-East J. Sci. Res. 2013, 17, 641–646. [Google Scholar] [CrossRef]
  17. El-Bassiouny, A.H.; Abdo, N.F.; Shahen, H.S. Exponential Lomax Distribution. Int. J. Comput. Appl. 2015, 121, 24–29. [Google Scholar] [CrossRef]
  18. Tahir, M.H.; Cordeiro, G.M.; Mansoor, M.; Zubair, M. The Weibull-Lomax distribution: Properties and applications. Hacet. J. Math. Stat. 2015, 44, 461–480. [Google Scholar] [CrossRef]
  19. Afify, A.Z.; Yousof, H.M.; Cordeiro, G.M.; Ortega, E.M.; Nofal, Z.M. The Weibull Fréchet distribution and its applications. J. Appl. Stat. 2016, 43, 2608–2626. [Google Scholar] [CrossRef]
  20. Rady, E.H.A.; Hassanein, W.A.; Elhaddad, T.A. The power Lomax distribution with an application to bladder cancer data. SpringerPlus 2016, 5, 1–22. [Google Scholar] [CrossRef] [Green Version]
  21. Anwar, M.; Zahoor, J. The Half-Logistic Lomax Distribution for Lifetime Modeling. J. Probab. Stat. 2018, 3152807. [Google Scholar] [CrossRef]
  22. Hassan, A.S.; Abd-Allah, M. On the Inverse Power Lomax Distribution. Ann. Data Sci. 2019, 6, 259–278. [Google Scholar] [CrossRef]
  23. Oguntunde, P.E.; Khaleel, M.A.; Okagbue, H.I.; Odetunmibi, O.A. The Topp–Leone Lomax (TLLo) Distribution with Applications to Airbone Communication Transceiver Dataset. Wirel. Pers. Commun. 2019, 109, 349–360. [Google Scholar] [CrossRef]
  24. Al-Marzouki, S.; Jamal, F.; Chesneau, C.; Elgarhy, M. Type II Topp Leone Power Lomax Distribution with Applications. Mathematics 2020, 8, 4. [Google Scholar] [CrossRef] [Green Version]
  25. Nagarjuna, V.B.V.; Vardhan, V. Marshall-Olkin Exponential Lomax Distribution: Properties and its Application. Stoch. Model. Appl. 2020, 24, 161–177. [Google Scholar]
  26. Mathew, J.; Chesneau, C. Some New Contributions on the Marshall–Olkin Length Biased Lomax Distribution: Theory, Modelling and Data Analysis. Math. Comput. Appl. 2020, 25, 79. [Google Scholar] [CrossRef]
  27. Nair, N.U.; Sankaran, P.; Balakrishnan, N. Quantile-Based Reliability Analysis; Birkhäuser: Basel, Switzerland, 2013. [Google Scholar]
  28. Rényi, A. On Measures of Entropy and Information; Technical Report; Hungarian Academy of Sciences: Budapest, Hungary, 1961. [Google Scholar]
  29. Havrda, J.; Charvát, F. Quantification method of classification processes. Concept of structural a-entropy. Kybernetika 1967, 3, 30–35. [Google Scholar]
  30. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  31. Bader, M.; Priest, A. Statistical aspects of fibre and bundle strength in hybrid composites. In Progress in Science and Engineering Composites; Hayashi, T., Kawata, K., Umekawa, S., Eds.; ICCM-IV: Tokyo, Japan, 1982; pp. 1129–1136. [Google Scholar]
  32. Louzada-Neto, F.; Mazucheli, J.; Achcar, J.A. Uma Introdução à análise de Sobrevivência e Confiabilidade; Sociedad Chilena de Estadística: Valparaíso, Chili, 2001. [Google Scholar]
Figure 1. Curves of the pdf of the KPL distribution at different parameter values.
Figure 1. Curves of the pdf of the KPL distribution at different parameter values.
Stats 04 00003 g001
Figure 2. Shapes of the hrf of the KPL distribution at different parameter values.
Figure 2. Shapes of the hrf of the KPL distribution at different parameter values.
Stats 04 00003 g002
Figure 3. Fitted pdf curves of the distributions to the histograms for (a) data set 1 and (b) data set 2.
Figure 3. Fitted pdf curves of the distributions to the histograms for (a) data set 1 and (b) data set 2.
Stats 04 00003 g003
Figure 4. Fitted PP plots of the distributions for data set 1.
Figure 4. Fitted PP plots of the distributions for data set 1.
Stats 04 00003 g004
Figure 5. Fitted PP plots of the distributions for data set 2.
Figure 5. Fitted PP plots of the distributions for data set 2.
Stats 04 00003 g005
Table 1. Moments about the origin and variance with different parameter values of the KPL distribution.
Table 1. Moments about the origin and variance with different parameter values of the KPL distribution.
Parametersa μ 1 μ 2 μ 3 μ 4 Var
b = 1    α = 2    λ = 3    β = 4 11.0963431.3603501.8989223.0000000.1583815
21.3133281.8704812.9077245.0000000.1456513
31.4341322.1999403.6411456.6000000.1432048
41.5181952.4483634.2336607.9714290.1434478
51.5828602.6500834.7384489.1904760.1446359
b = 2    α = 4    λ = 6    β = 8 10.91671150.86067790.82514980.80597740.020317974
21.02865841.06950941.12330401.19124260.011371281
31.08089451.17688141.29050141.42488880.008548480
41.11371161.24752271.40535241.59201040.007169093
51.13722471.29962381.49242691.72208340.006343810
b = 1.5    α = 1.5    λ = 2    β = 5 10.95153460.97711631.0759391.2676210.07169832
21.13987211.36090601.7037222.2423010.06159767
31.24285131.60396522.1544673.0225830.05928593
41.31417951.78596222.5167013.6913750.05889455
51.36900081.93334372.8246484.2862710.05918053
b = 1.5 >   α = 1.5 >   λ = 0.5 >   β = 2 30.92466641.1130721.9363377.6510630.2580643
71.33938902.2146265.04625224.9030520.4206629
101.54521392.9090267.44782540.6923160.5213404
141.75984843.73789810.69072864.5045850.6408318
171.89377664.30934413.14166884.0655160.7229542
b = 3.5    α = 0.8    λ = 2.5    β = 2 32.3818987.09844827.62178156.15131.425008
74.20973421.533759141.107211329.59623.811897
105.29148033.862448276.773663248.83815.862686
146.54754351.716382520.950117533.24748.846066
177.39818165.967625749.8002612233.855211.234540
b = 2.7    α = 0.8    λ = 0.5    β = 2.2 41.4201612.5423236.28713728.972990.5254662
112.5733738.20459535.761130289.640941.5823481
173.30104713.47325475.107768778.338052.5763424
223.82293018.059557116.5025751397.612703.4447658
365.05648931.580735269.3640024275.425756.0126525
Table 2. Bias in parenthesis and MSEs for different sample sizes in the context of the KPL distribution.
Table 2. Bias in parenthesis and MSEs for different sample sizes in the context of the KPL distribution.
( a , b , α , λ , β ) a ^ b ^ α ^ β ^ λ ^
( 10 , 4 , 5 , 0.5 , 20 )
n = 50 (8.651935)(6.67144)(4.366995)(−6.868542)(1.093022)
74.8559944.5081119.0706547.176871.194697
n = 100 (8.002526)(4.365969)(−2.73485)(0.207459)(−0.3513344)
64.0440619.084347.4799760.086664410.1234489
n = 250 (6.22343)(−0.751908)(−1.357911)(2.745664)(−0.3055339)
38.737970.56798161.8459817.5512920.09335099
( 4 , 2 , 5 , 0.5 , 1.5 )
n = 50 (−3.279042)(−1.111076)(−4.199035)(3.525116)(−0.4977061)
10.752121.2344917.631912.426450.2477113
n = 100 (−3.003313)(−0.9706962)(−4.235572)(0.7994587)(−0.4568123)
9.1242031.13406817.940190.6410230.2086831
n = 250 (−0.2668285)(−1.074196)(−2.093235)(0.7068528)(−0.4322879)
0.076454461.1875024.38660.49970610.1868756
Table 3. MLEs of the considered distribution parameters for the data sets.
Table 3. MLEs of the considered distribution parameters for the data sets.
Distribution Data Set 1 a ^ b ^ α ^ λ ^ β ^
KPL21.10826.54185.19570.44852.2026
KBXII0.43420.26397.57591.96394.4371
PL -3.457413.57593.1625
WL4.04052.7110.7077-1.7575
WFr6.46476.86520.7611-0.2297
TLGL4.8177-12.3340-16.0790
EL5.4551-21.9720-13.2469
Lomax--12.3594-17.9267
Distribution Data Set 2 a ^ b ^ α ^ λ ^ β ^
KPL8.99537.43962.05065.55810.2604
WFr0.08233.83035.0670-0.1522
WL0.00560.16454.3790-1.2422
KBXII8.35956.46010.05804.19351.5752
EL6.1659-0.6294-17.4538
TLGL9.5741-0.3008-7.2733
PL--0.918626.50700.5654
Lomax--0.2847-24.8517
Table 4. The goodness of fit tests and adequacy values for the data sets.
Table 4. The goodness of fit tests and adequacy values for the data sets.
Distribution
Data Set 1W*A* D n p-ValueAICCAICBICHQIC
KPL0.02320.22720.06040.9605125.6272126.5647136.8697130.0929
KBXII0.06450.51080.07470.8287129.838130.7755141.0804134.3036
PL0.13450.94500.07780.7909131.9917132.3553138.7372134.6711
WL0.18601.29950.10300.4476138.8186139.4340147.8126142.3911
WFr0.21861.49560.12050.2610141.3542141.9696150.3482144.9267
TLGL0.40602.52910.14430.1085152.7280153.0916159.4734155.4073
EL0.42652.64400.14400.1098153.4122153.7758160.1577156.0916
Lomax0.32132.05400.35544.18 × 10 08 204.3163204.4954208.8133206.1026
Data Set 2W*A* D n p-ValueAICCAICBICHQIC
KPL0.09360.47400.07260.55302035.7682036.2942049.7052041.428
WFr0.15230.87610.09570.22212039.4422039.7902050.5922043.970
WL0.19071.11220.11750.07302043.7632044.1112054.9132048.292
KBXII0.31341.76920.15580.00592060.2512060.7782074.1892065.912
EL0.57173.24530.13990.01822072.0462072.2522080.4082075.442
TLGL0.61553.52410.13720.02182077.8532078.0602086.2152081.249
PL0.14880.75850.24867.22 × 10 7 2117.3922117.5992125.7552120.788
Lomax0.41672.32680.31101.65 × 10 10 2157.8682157.9712163.4432160.132
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nagarjuna, V.B.V.; Vardhan, R.V.; Chesneau, C. Kumaraswamy Generalized Power Lomax Distributionand Its Applications. Stats 2021, 4, 28-45. https://0-doi-org.brum.beds.ac.uk/10.3390/stats4010003

AMA Style

Nagarjuna VBV, Vardhan RV, Chesneau C. Kumaraswamy Generalized Power Lomax Distributionand Its Applications. Stats. 2021; 4(1):28-45. https://0-doi-org.brum.beds.ac.uk/10.3390/stats4010003

Chicago/Turabian Style

Nagarjuna, Vasili B.V., R. Vishnu Vardhan, and Christophe Chesneau. 2021. "Kumaraswamy Generalized Power Lomax Distributionand Its Applications" Stats 4, no. 1: 28-45. https://0-doi-org.brum.beds.ac.uk/10.3390/stats4010003

Article Metrics

Back to TopTop