Research Article  Open Access
The Exponentiated HalfLogistic Family of Distributions: Properties and Applications
Abstract
We study some mathematical properties of a new generator of continuous distributions with two extra parameters called the exponentiated halflogistic family. We present some special models. We investigate the shapes of the density and hazard rate function. We derive explicit expressions for the ordinary and incomplete moments, quantile and generating functions, probability weighted moments, Bonferroni and Lorenz curves, Shannon and Rényi entropies, and order statistics, which hold for any baseline model. We introduce two bivariate extensions of this family. We discuss the estimation of the model parameters by maximum likelihood and demonstrate the potentiality of the new family by means of two real data sets.
1. Introduction
The use of new generators of continuous distributions from classic distributions has become very common in recent years. One example is the betagenerated family of distributions proposed by Eugene et al. [4]. Another example is the gammagenerated family of distributions defined by Zografos and Balakrishnan [5]. Based on a baseline continuous distribution with survival function and density , their families are defined by the cumulative distribution function (cdf) and probability density function (pdf) (for ): respectively, where is the gamma function.
Based on Zografos and Balakrishnan’s [5] paper, we replace the gamma distribution by the exponentiated halflogistic (“EHL” for short) distribution to define a new family of continuous distributions by the cdf: where is the baseline cdf depending on a parameter vector and and are two additional shape parameters. For any continuous distribution, the EHL distribution is defined by the cdf (2). Equation (2) is a wider family of continuous distributions and includes some special models as those listed in Table 1.

The density function corresponding to (2) is given by where is the baseline pdf. Equation (3) will be most tractable when and have simple analytic expressions. Hereafter, a random variable with density function (3) is denoted by . Further, we can omit sometimes the dependence on the vector of the parameters and simply write .
A physical interpretation of the EHL distribution can be given as follows. Consider a system formed by independent components having the halflogistic (“HL”) cdf given by Suppose that the system fails if all of the components fail and let denote the lifetime of the entire system. Then, the pdf of is given by (3).
The hazard rate function (hrf) of becomes The EHL family of distributions is easily simulated by inverting (2) as follows: if has a uniform distribution, the solution of the nonlinear equation has the density function (3).
This paper is organized as follows. In Section 2, some special cases of the EHL family of distributions are defined. In Section 3, the shapes of the density and hazard rate functions are described analytically. A useful expansion for the new density family is obtained and we derive a power series for the EHL quantile function in Section 4. General explicit expressions for some special EHL moments are obtained in Section 5.
In Section 6, we derive the generating function, the incomplete moments are investigated, we obtain the mean deviations and the reliability and provide expressions for the Rényi and Shannon entropies, and the order statistics and their moments are determined. We introduce two bivariate extensions of the new family in Section 7. Estimation of the model parameters by maximum likelihood is performed in Section 8. Applications to two real data sets illustrate the performance of the new family in Section 9. The paper is concluded in Section 10.
2. Special EHL Models
Here, we introduce only three of the many distributions which can arise as EHL special models, where and are positive shape parameters of the new generator. We consider three baseline distributions, namely, Fréchet, loglogistic, and generalized halfnormal distributions, although we can generate as many new distributions as desirable.
2.1. Exponentiated HalfLogisticFréchet (EHLF) Model
The Fréchet (or type II extreme value) distribution has been useful for modeling of marketreturns which are often heavytailed in applications to finance [6]. Now, we introduce a new fourparameter distribution called the EHLF distribution. Taking to be the Fréchet distribution with scale parameter and shape parameter , where , the EHLF density function (for ) is given by
The cdf and hrf corresponding to (7) are given by respectively. A characteristic of the EHLF distribution is that its hrf can be monotonically increasing or decreasing and upsidedown bathtub depending basically on the parameter values. Plots of its density function and hrf for some parameter values are displayed in Figures 1 and 2, respectively.
(a)
(b)
(c)
(a)
(b)
2.2. Exponentiated HalfLogisticLogLogistic (EHLLL) Model
The loglogistic (LL) distribution is widely used in practice and it is an alternative to the lognormal distribution since it presents a failure rate function that increases, reaches a peak after some finite period, and then declines gradually. The properties of the LL distribution make it an attractive alternative to the lognormal and Weibull distributions in the analysis of survival data [7]. This distribution can exhibit a monotonically decreasing failure rate function for some parameter values. For , let be the LL cdf, where is the shape parameter and is the scale parameter, where . The EHLLL density function becomes
In Figure 3, we display some possible shapes of the EHLLL density function. The corresponding cdf and hrf are given by respectively. Plots of the EHLLL hrf for some parameter values are displayed in Figure 4.
(a)
(b)
(c)
(a)
(b)
2.3. Exponentiated HalfLogistic Generalized HalfNormal (EHLGHN) Model
The most popular models used to describe the lifetime process under fatigue are the halfnormal (HN) and BirnbaumSaunders (BS) distributions. When modeling monotone hazard rates, the HN and BS distributions may be an initial choice because of their negatively and positively skewed density shapes. Consider to be the generalized halfnormal (GHN) distribution [8] with scale parameter and shape parameter , where , given by , where is the error function. Note that Then, the fourparameter EHLGHN density (for ) can be expressed as
If , the EHLGHN distribution model reduces to the exponentiated halflogistic halfnormal (EHLHN) distribution. The cdf and hrf corresponding to (12) are respectively. A characteristic of the EHLGHN distribution is that its hrf can be bathtub shaped, monotonically increasing or decreasing, and upsidedown bathtub depending basically on the parameter values. Plots of the EHLGHN density function and hrf for some parameter values are displayed in Figures 5 and 6, respectively.
(a)
(b)
(c)
(a)
(b)
3. Shapes
The shapes of the density and hazard rate functions can be described analytically. The critical points of the EHL density function are the roots of the equation: There may be more than one root to (14). Let . We have If is a root of (14), then it corresponds to a local maximum if for all and for all . It corresponds to a local minimum if for all and for all . It refers to a point of inflexion if either for all or for all .
The critical point of the hrf of , say , is obtained from the following equation:
There may be more than one root to (16). Let . We have If is a root of (16), then it refers to a local maximum if for all and for all . It corresponds to a local minimum if for all and for all . It gives an inflexion point if either for all or for all .
4. A Useful Expansion and Quantile Power Series
We can demonstrate that the cdf of given by (2) admits the following expansion: where denotes the exponentiated (“exp”) cumulative distribution with power parameter , Some structural properties of the exp distributions are investigated by Mudholkar et al. [9], Gupta and Kundu [10], and Nadarajah and Kotz [11], among others.
The density function of can be expressed as an infinite linear combination of exp density functions: where denotes the density function of the exp random variable with power parameter . Equation (20) reveals that the EHL density function is a linear combination of exp density functions. Thus, some mathematical properties of the new family can be obtained directly from those properties of the exp distribution.
Here, we derive a power series expansion for the quantile function of by expanding (6). If the quantile function, say , does not have a closedform expression, it can usually be expressed in terms of a power series where the coefficients are suitably chosen real numbers which depend on the parameters of the distribution. For several important distributions, such as the normal, the Student , and gamma and beta distributions, does not have explicit expressions but it can be expanded as in (21). As a simple example, for the normal distribution, for and , , , and .
We use throughout the paper a result of Gradshteyn and Ryzhik ([12], Section 0.314) for a power series raised to a positive integer (for ): where the coefficients (for ) are easily obtained from the recurrence equation (with ): Clearly, can be determined from and then from the quantities .
Next, we derive an expansion for the argument of in (6):
Using the generalized binomial expansion four times since , we can write and then where and for , , and Then, the quantile function of can be expressed from (6) as where for and For any baseline distribution, we can combine (21) with (28) to obtain and then using (22) and (23), we have where , , and, for , Equation (32) is the main result of this section since it allows to obtain various mathematical quantities for the EHL family as shown in the next sections.
The formulae derived throughout the paper can be easily handled in most symbolic computation software platforms such as Maple, Mathematica, and MATLAB. These platforms currently have the ability to deal with analytic expressions of formidable size and complexity. Established explicit expressions to calculate statistical measures can be more efficient than computing them directly by numerical integration. The infinity limit in these sums can be substituted by a large positive integer such as 20 or 30 for most practical purposes.
5. Moments
Hereafter, we will assume that is the cdf of a random variable and that is the cdf of the random variable having density function (3). The moments of can be obtained from the th probability weighted moments (PWMs) of given by
An alternative expression for can be determined using (22) and (23): The PWMs for several distributions can be calculated from (34) and (35).
We can write from (20) Thus, the moments of any EHL distribution can be expressed as an infinite weighted linear combination of the baseline PWMs. Equations (34)–(36) are the main results of this section.
Further, the central moments () and cumulants () of can be calculated as respectively, where . Then, , , , and so forth. The skewness and kurtosis quantities follow from the second, third, and fourth cumulants.
5.1. EHLF Model
Consider the Fréchet baseline cdf for and corresponding pdf discussed in Section 2.2. The EHLF density function can be written from (20) as where . This equation reveals that the EHLF density function can be expressed as an infinite mixture of Fréchet densities.
The th PWM of the Fréchet distribution becomes Setting , reduces to The integral converges absolutely for and then
Plots of the skewness and kurtosis for some choices of as functions of , for and , are displayed in Figure 7.
(a)
(b)
5.2. Exponentiated HalfLogistic Logistic (EHLLo) Model
For the EHLLo distribution, the baseline cumulative function is . Using a result from Prudnikov et al. ([13], Section 2.6.13, equation (4)), we can write from (34) (for ) the following: where is the beta function. The th moment of the EHLLo distribution comes from (36) as
5.3. Exponentiated HalfLogistic Gamma (EHLGa) Model
Using the power series expansion for the gamma cdf we obtain from (20) the following series expansion:
The EHLGa moments follow from (36) and the expression for given by
5.4. Exponentiated HalfLogistic Normal (EHLN) Model
The moments of can be obtained from the moments of using , and then we can work with the standard normal distribution. We can expand the EHLN cumulative function (18) (with and ) as From the series expansion for the error function we obtain a series expansion from (20) (with and ) given by The EHLN moments can be obtained from (36) and the PWMs given by Cordeiro and Nadarajah [14]. Plots of the skewness and kurtosis for some choices of as functions of , for , and , , are displayed for the EHLLL and EHLHGN distributions in Figures 8 and 9, respectively. These plots show that the skewness and kurtosis are very flexible.
(a)
(b)
(a)
(b)
6. Other Measures
In this section, we calculate the following measures: generating function, incomplete moments, mean deviations, reliability, entropies, and order statistics for the EHL family.
6.1. Generating Function
Here, we provide two formulae for the moment generating function (mgf) of . A first formula for comes from (20) as where is the generating function of the exp distribution with power parameter . Hence, can be determined from the exp generating function.
A second formula for can be derived from (20) as where
We can derive the mgf’s of several EHL distributions directly from (50)(51). For example, the mgf’s of the exponentiated halflogistic exponential (EHLE) (with parameter and ) and EHLLo (with ) distributions are given by respectively.
Clearly, two representations for the characteristic function (chf) of can be derived from (50)–(52) by , where .
6.2. Incomplete Moments
Incomplete moments of the income distribution form natural building blocks for measuring inequality. For example, the Lorenz and Bonferroni curves depend upon the incomplete moments of the income distribution. The th incomplete moment of is defined as . Here, we provide two formulae to calculate the incomplete moments of the EHL family. First, the th incomplete moment of can be expressed as
The integral in (54) can be computed at least numerically for most baseline distributions.
A second formula follows from (54) using (22) and (23). We can write where is given by (23).
The first incomplete moment can be used to obtain Bonferroni and Lorenz curves defined for a given probability by and , respectively, where is immediately calculated from the parent quantile function.
6.3. Mean Deviations
The mean deviations about the mean () and about the median () of can be expressed as respectively, where is the median of , and come from (2) and (36), respectively, and is the first incomplete moment.
Now, we provide two alternative ways to compute and . A general equation for can be derived from (20) as where Equation (58) is the basic quantity to compute the mean deviations for the EGL distributions.
A second general formula for can be derived by setting in (57): where
Equations (55)–(59) are the main results of this section.
6.4. Reliability
Here, we derive the reliability when and are independent random variables with a positive support. It has many applications especially in engineering concepts. Let denote the pdf of and let denote the cdf of . By expanding the binomial terms in and , we obtain where If , we obtain Further, if and , then .
6.5. Entropies
An entropy is a measure of variation or uncertainty of a random variable . Two popular entropy measures are the Rényi and Shannon entropies. The Rényi entropy of a random variable with pdf is defined (for and ) as The Shannon entropy of a random variable is given by , which is the special case of the Rényi entropy when . Direct calculation gives
After some algebraic manipulations, we obtain the following.
Proposition 1. Let be a random variable with pdf given by (3). Then,
The simplest formula for the entropy of becomes After some algebraic developments, we obtain an alternative expression for : where .
6.6. Order Statistics
Order statistics make their appearance in many areas of statistical theory and practice. Suppose that is a random sample from the EHL distribution. Let denote the th order statistic. From (18) and (20), the pdf of is given by where . Using (22) and (23), we can write where and Hence, where .
Equation (72) is the main result of this section. It reveals that the pdf of the EHL order statistics is a linear combination of exp density functions. So, several structural quantities of the EHL order statistics like ordinary, incomplete moments, generating function, mean deviations, and several others can be obtained from the corresponding quantities of exp distributions.
7. Bivariate Extensions
In this section, we introduce two extensions of the proposed model. The first extension is based on the idea of [15]. Let , , and be independent random variables. Further, we define and . Then, the pdf of the bivariate random variable is given by where . The marginal cdf’s are Clearly, if we consider and , the pdf of is given by