A weighted transmuted exponential distribution with environmental applications

02/08/2020 ∙ by Christophe Chesneau, et al. ∙ 0

In this paper, we introduce a new three-parameter distribution based on the combination of re-parametrization of the so-called EGNB2 and transmuted exponential distributions. This combination aims to modify the transmuted exponential distribution via the incorporation of an additional parameter, mainly adding a high degree of flexibility on the mode and impacting the skewness and kurtosis of the tail. We explore some mathematical properties of this distribution including the hazard rate function, moments, the moment generating function, the quantile function, various entropy measures and (reversed) residual life functions. A statistical study investigates estimation of the parameters using the method of maximum likelihood. The distribution along with other existing distributions are fitted to two environmental data sets and its superior performance is assessed by using some goodness-of-fit tests. As a result, some environmental measures associated with these data are obtained such as the return level and mean deviation about this level.



There are no comments yet.


page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

The precise analysis of a wide variety of data sets is limited by the use of models based on the classical distributions (normal, exponential, logistic…). For instance, the analysis of environmental data sets collecting from observations of complex natural phenomena needs special treatments to reveal all the underlying informations. Over the last decades, numerous solutions have been provided by the statisticians, including the elaboration of several methods which aim to increase the flexibility of the former classical distributions. Among these methods, a popular one that aims to construct a generator of distributions by compounding continuous distributions with well-known discrete distributions. This compounding is always motivated by practical problems as those involving cdf of minimum or maximum of several independent and identically random variables. An exhaustive survey on the construction of such generators, with the presentation of new ones, can be found in

[22], and the references therein. Among the long list, let us briefly present the EGNB2 distribution introduced by [22, Remark 2 (ii)]

. Using a cumulative distribution function (cdf)

, the general form of the associated cdf is given by


The EGNB2 distribution can be viewed as an extension of the G-negative binomial families introduced by [9] and [17]. It enjoys remarkable theoretical and practical properties.

In this study, we consider a particular case of this EGNB2 distribution consisting in a re-parametrization for the parameters , and appearing in (1) as described below. Let , , and . That yields a cdf of the (simple) form:


Let us now explain the importance of this re-parametrization of (1), with some statistical features. One can observe that as the following integral form: , where denotes the pdf: . So it reveals to be a new particular case of the T-X family cdf introduced by [4]. Another remark is that, when , we have and when , we have . This transformation of cdf corresponds to the one proposed in [10]

. All the resulting distributions have demonstrated nice properties in terms of analysis of real life data sets. Furthermore, let us observe that the probability density function (pdf) associated to (

2) is given by

Note that we can also express it as a weighted pdf: , where is a weight function and is a normalizing constant. It thus belongs to the family of weighted distributions. Further details on such family of distributions can be found in [19]. On the other side, [5] introduced the transmuted exponential distribution defined by the following cdf: , , where denotes the cdf of the exponential distribution. Then, it is proved that the additional parameter can significantly increase the flexibility of the former exponential distribution, demonstrating a superiority in terms of fit in comparison to the former exponential distribution. We may refer the reader to [16], and the references therein.

In this paper, we introduce a new three-parameter distribution which combines the features of the distribution characterized by (2) and the transmuted exponential distribution. This combination aims to modify the former transmuted exponential distribution by incorporating the parameter

and takes benefit of the flexibility of the EGNB2 distribution. Its main role is to add a high degree of flexibility on the mode, and the skewness and kurtosis of the tail. We thus obtain a very flexible distribution, which opens new perspectives in terms of the construction of statistical models for data analysis. The theoretical and practical aspects are explored in an exhaustive way. The theoretical ones include expansions of the cdf, pdf, hazard rate function (hrf), quantile function, moments, moment generating function, various entropy measures, residual life functions, conditional moments, mean deviations and reversed residual life function. We investigate the estimation of its parameters via the maximum likelihood method. Two real-life data sets in environmental sciences are analyzed to show its superior performance in terms of fit in comparison to well-known distributions: The gamma distribution, the Marshal-Olkin exponential distribution

[11], the Nadarajah-Haghighi exponential distribution [15], the exponentiated exponential distribution [7], the transmuted Weibull distribution [5], the transmuted generalized exponential distribution [8], the transmuted linear exponential distribution [23] and the Kappa distribution [13]. The best performance of the proposed distribution recommends it as a hydrologic probability model, such as the most known distributions: Kappa and gamma distributions. This motivates to estimate important hydrologic parameters of those data sets by making use of the distribution.

The rest of this article is organized as follows. In Section 2, we present our main distribution. Some of its mathematical properties are studied in Section 3. Residual life functions are determined in Section 4. Estimations of the parameters are investigated in Section 5. Applications to two real-life data sets are provided in Section 6. Concluding remarks are addressed in Section 7.

2 A new weighted transmuted exponential distribution

In this section, we precise what is the considered cdf given by (2). [21] and [5] introduced the quadratic rank transmutation map (QRTM) to propose a new distribution based on the Weibull/exponential one with great flexibility and nice fit for real-life data. In the current studies, it remains a serious competitor in terms of precision in modelling (see [16]). For these reasons, we use it in our study. We consider the cdf:

where is considered to be the cdf of the exponential distribution of parameter :

Set the above expression into (2), we introduce a new cdf defined by

Another useful expression is the following one:


We will refer to the distribution given by (3) as the new weighted transmuted exponential and denote it by NWTE(, , ) with the considered parameters.

The corresponding pdf is given by


The associated hrf is given by


Let us now discuss the possible shapes of pdf (4) and hrf (5) as follows.

On the other side, we have

In order to visualize the wide variety of shapes, some plots of the pdf (4) and hrf (5) are given in Figures 1 and 2. We see that has a great impact on the mode of the NWTE distribution. Moreover, the hrf also exhibits sudden spikes at the end of upside-down bathtub shapes, which manages the model to analyze a non-stationary real-life data.

Figure 1: Plots of the NWTE pdf.
Figure 2: Plots of the NWTE hrf.

3 Structural properties of the NWTE distribution

3.1 Expansion for the associated functions

Expansion for the cdf function. First of all, set , , . Note that we have , so is increasing. Since and , we have for all . Since and , the generalized binomial expansion, we have



Therefore we can expand the cdf function as


Expansion for the pdf function. Similar mathematical arguments used for (6) give





On the survival function. Note that


Using (7), we have the following expansion


Expansion for the hrf function. Using (5), (8) and (10), an expansion of the hrf function is given by


Another expansion comes from the geometric series decomposition:

By (5) and similar mathematical arguments used for (6) give:


3.2 Quantile function

The quantile functions are in widespread use in general statistics to obtain mathematical properties of a distribution and often find representations in terms of lookup tables for key percentiles. For generating data from the NWTE model, let . Then, by inverting the cdf (3) and after some algebra, we get the quantile function


The analysis of the variability of the skewness and kurtosis of X can be investigated based on quantile measures. The Bowley skewness is given by

and the Moors’ kurtosis by

where is given by (12).

These measures are less sensitive to outliers and they exist even for distributions without moments. Figure

3 displays plots of S and K as functions of and , which show their variability in terms of the shape parameters.

Figure 3: Plots of the skewness and kurtosis of the NWTE distribution for .

3.3 Moments and moment generating function

Moments. Using equation (8) and the gamma function , the -th moments about the origin is given by


The moment generating function. Similarly the moment generating function associated to the NWTE distribution is given by, for ,


3.4 Entropies

An entropy can be considered as a measure of uncertainty of probability distribution of a random variable. Therefore, we obtain three entropies for the NWTE distribution with investigating a numerical study among them.

Entropy 1. Let us consider the Shannon entropy [20]: . One can observe that


Let us now expand the two integrals by using the logarithmic expansion: , . Since , we have


denotes the moment generating function defined by (3.3).

For the second integral in (3.4), since , we have


Entropy 2. Let us now focus our attention on the Rényi entropy [18]: , with and . Similar mathematical arguments used for (6) give :


On the other side, observing that , similar mathematical arguments used for (6) give :


Hence can be expanded as




Entropy 3. We now focus our attention on the entropy introduced by [12]: , with and . Proceeding as for with instead of , we obtain



Some numerical values for the three entropies are given in Table 1. It can be observed that these entropies decrease with increasing the parameter values. Moreover, one can see that has the smallest values comparing with the other entropies considered here.

0.1 0.94207 1.32902 0.63437
0.4 0.91244 1.30828 0.61079
0.8 0.88415 1.28871 0.58774
1.2 0.86349 1.27456 0.57056
1.5 0.85122 1.26622 0.56021
1.8 0.84092 1.25926 0.55144
2.0 0.83492 1.25522 0.54629
-0.9 1.40499 1.66797 0.94970
-0.5 1.33563 1.61265 0.91027
-0.2 1.24327 1.54839 0.84921
0.1 1.11731 1.46168 0.75998
0.4 0.95418 1.34297 0.64014
0.6 0.82110 1.23549 0.54159
0.8 0.66351 1.08602 0.42654
Table 1: Entropy for several arbitrary parameter values with .

3.5 Conditional moments and mean deviations

Here, we introduce an important lemma which will be used in the next sections.

Lemma 1.

Let and be the lower incomplete gamma function. Then we have


Using the equation (8), we have

The -th conditional moments of the NWTE distribution is given by


It can be expressed using (5), (3.3) and Lemma 1. The same remark holds for the -th reversed moments of the NWTE distribution given by

The mean deviations of about the mean can be expressed as and the mean deviations of about the median has the form .

4 (Reversed) Residual life functions

4.1 Residual lifetime function

The residual life is described by the conditional random variable , . Using (10), the survival function of the residual lifetime for the NWTE distribution is given by

The associated cdf is given by

The corresponding pdf is given by

The associated hrf is given by

The mean residual life is defined as

where is given by (4), is mentioned in (9), is given by (3.3) and is stated in Lemma 1.

Further, the variance residual life is given by

where is given by (3.3) and is given by Lemma 1. Some numerical values for the mean residual life are displayed in Table 2 for various choices of the parameters and at the time points It can be seen that, the mean residual life increases with increasing the time points t, also decreases with increasing and .

1.0 3.0 5.0 7.0 10
0.1 0.918095 0.982144 0.997423 0.999648 0.999982
0.7 0.900668 0.980089 0.997152 0.999612 0.999981
1.1 0.894316 0.979369 0.997057 0.999599 0.999980
1.6 0.889012 0.978778 0.996979 0.999588 0.999979
2.0 0.885996 0.978447 0.996936 0.999582 0.999979
1.0 3.0 5.0 7.0 10
-0.9 1.219460 1.027861 1.003735 1.000504 1.000025
-0.5 1.162075 1.020915 1.002810 1.000380 1.000018
-0.1 1.088250 1.011465 1.001542 1.000208 1.000010
0.1 1.040826 1.004807 1.000638 1.000086 1.000004
0.5 0.905030 0.980593 0.997218 0.999621 0.999981
1.0 0.511566 0.500207 0.500004 0.500000 0.500001
Table 2: Mean residual life function for arbitrary parameter values with .

4.2 Reversed residual life function

The reverse residual life is described by the conditional random variable , . Using (3), the survival function of the reversed residual lifetime for the NWTE distribution is given by