A new extended gumbel distribution: Properties and application.

A robust generalisation of the Gumbel distribution is proposed in this article. This family of distributions is based on the T-X paradigm. From a list of special distributions that have evolved as a result of this family, three separate models are also mentioned in this article. A linear combination of generalised exponential distributions can be used to characterise the density of a new family, which is critical in assessing some of the family's properties. The statistical features of this family are determined, including exact formulations for the quantile function, ordinary and incomplete moments, generating function, and order statistics. The model parameters are estimated using the maximum likelihood method. Further, one of the unique models has been systematically studied. Along with conventional skewness measures, MacGillivray skewness is also used to quantify the skewness measure. The new probability distribution also enables us to determine certain critical risk indicators, both numerically and graphically. We use a simulated assessment of the suggested distribution, as well as apply three real-world data sets in modelling the proposed model, in order to ensure its authenticity and superiority.

1 Introduction

The employment of traditional probability models to anticipate real-life occurrences is causing increasing dissatisfaction among applied practitioners. Tail characteristics and goodness of fit metrics may have a constraining tendency, which could be one of the reasons. As a response, in recent years, there has been a substantial rise in the generalisation of well-known probability distributions. The challenge is to find such versatile families that can fit both skew and symmetric data. It’s important to realize that the majority of generalised distributions described in the literature are constructed using the generalised classes approach (G-class) and the compounding principle. In [1], the authors provided a concise review of literature regarding generalization of distributions and transformation through versatile parameter induction techniques. We offer to the readers a few, but not exhaustive, lists in chronological sequence: [2-14].

According to [15], Emil J. Gumbel originated the use of the Gumbel distribution (GuD) on data bearing extreme values. By “extreme data” we mean the behaviour of a random variable that occurs at the sample threshold level and is seen using intense data and insights. In the reference [16], the authors defined the GuD, alternatively known as extreme value distribution-type I (γ = 0), as the predominant model for quantifying extreme occurrences such as flood frequency analysis, droughts, earthquakes, sea currents and wind speed in order to understand the trajectory, magnitude, and pattern of complex phenomena. Environmental sciences, geology, accelerated life testing, meteorology, risk assessment and epidemiology are just a few of the a set of diverse fields where it might well be utilized. The authors in [17] showed that the score statistics of global sequence alignment follows a Gumbel distribution. In the reference [18], a comprehensive list of real life scenarios to which GuD can be applicable is provided by the authors. To learn more, see references [19-24].

Earliest generalizations of GuD was reported by [25] by introducing a shape parameter to Gumbel distribution. [26] provided a generalization of GuD based on the asymptotic distribution of the mth extreme, tracing back to [27]. [28] emphasized on a trivial choice of distribution since the GuD with only location and scale parameter yields narrower confidence intervals and has also the danger of underestimating the level of return. In the reference [29], a unique modification of GuD was proposed. It is based on the logarithmic transformation of an odd Weibull variable and is defined as where Δ = e(, −∞ ≤ x ≤ + ∞, −∞ ≤ μ ≤ + ∞ and 0 ≤ βσ ≤ + ∞.

Since then, researchers adopted a more formalistic approach to generalize GuD. Some notable generalizations include [30] to define Beta-Gumbel (BGu), [31] to propose Kumaraswamy-gumbel (KumGu), [32] to define exponentiated-Gumbel (EGu), [33] to define exponentiated-generalized Gumbel (EGGu), just to mention some.

[14] proposed a simplified approach to generalize any continuous distribution viz. a viz. the transformed-transformer (T-X) family, which has become an indispensable part of modern distribution theory. Let z(t) be the probability density function (pdf) and Z(t) be the cumulative distribution function (cdf) of a random variable (av) T such that tϵ(a1, a2) with support −∞ < a1 < a2 < ∞. Let W[Z(x)] act as generator function of the cdf Z(x) of any baseline av such that K[Z(x)] is differentiable and increasing, lies in the defined range, i.e. a1 ≤ K[Z(x)] ≤ a2 such that when x → −∞ as K[Z(x)] → a1 and x → + ∞ as K[Z(x)] → a2.

The pdf corresponding to Eq (2) is

To generalize any continuous distribution, the methodology defined by the cdf in Eq (2) has become indispensable part of modern distribution theory. In the same vein, Al-Aqtash et al. [34] introduced the Gumbel-X family of distributions. Let z(t) be the density function (pdf) and Z(t) be the distribution function (cdf) of an arbitrary variable (av) T such that tϵ(a1, a2) with support −∞ < a1 < a2 < ∞. Let K[Z(x; φ)] = log [(Z(x; φ))/1−(Z(x; φ))] act as the generator function of the cdf of any baseline av such that K[Z(x; φ)] fulfills the defined criterion. For μ = 0, the cdf of Gumbel-X family is given as where Δ = Z(x)/(1 − Z(x)).

This study introduces a new class of distributions following the T-X methodology, viz. a viz. the exponentiated Gumbel-G (EGuG) family of distributions. This is achieved by replacing the link function K[Z(x; φ)] = log[−log(1 − Z(x; φ))] in T-X family. It is worthy to remember that the link function log[−log(1 − Z(x; φ))] = log[−log(Z(x; φ))] and either of the mentioned form can be employ. EGuG family has thus far not been reported in the literature. We choose EGu distribution to define new family due to its superiority over the ordinary GuD because of presence of shape parameter θ that entails the improvements in tail of the distribution. Moreover, to the best of our knowledge, majority of the extreme value theory literature is supported by data from meteorology, geology, seismology, and hydrology (see references [15, 18–22, 24, 30–34]). The health implications of climate-related shifts in extreme event exposure, on the other hand, have not been explored. This study’s theoretical investigation will presumably fill this void in existing literature. Additionally, log[−log(1 − Z(x; φ))] function involves double log transformation and cannot be employed on GuD, which somehow makes the link function redundant. Following the success of the proposed generator to generalize Logistic and Normal distributions, we use this generator to define EGuG distribution. We study some of its mathematical properties and provide general properties and application of one of its specific model.

This article is outlined as follows: In Section 2, we define the EGuG family and present some of its special models. In order to optimise the structure of the generalisation being proposed, we provide the linear representation of EGuG density along with some of the mathematical properties of the family such as shapes of density and hazard rate function, moments and generating function, order statistics and estimation of model parameters. In Section 3, we choose Nadarajah-Haghighi (NH) distribution as baseline model to form EGuNH distribution whose mathematical properties as well as some risk measures are established. A simulation study is also conducted for some parametric combinations. Section 4 comprises of the numerical illustrations based on three life data sets. In Section 5, the article’s concluding thoughts are summed up.

2 The EGuG family

Let T follows the EGu av with μ = 0 and shape parameters θ ≥ 0 and σ ≥ 0, say EGu (θ, σ), then its cdf is given by The corresponding pdf to Eq (5) is given as For any baseline distribution with cdf Z(x; Φ) and pdf z(x; Φ) = dZ(x; Φ)/dx, the cdf of EGuG family is given as where θ, σ are shape parameters and φ is the vector of baseline parameter.

The pdf corresponding to Eq (7) is given by where Z(x; φ) is the baseline cdf and z(x; φ) is the baseline pdf. Furthermore, the dependence on the vector φ of the parameters might be omitted at times and simply write Z(x) = Z(x; φ) and z(x) = z(x; φ). Henceforth, X ∼ EGu−G(θ, σ; φ) denotes an av having density Eq (8).

The survival function (sf), hazard rate function (hrf) and cumulative hazard rate function (chrf) of this new family are, respectively, given by and

Simulating the EGuG family is simply done by inverting Eq (7) as follows: If U has a uniform U(0, 1) distribution, then has the density function Eq (8), where Q(.) = Z−1(.) is the baseline quantile function (qf).

2.1 Special models

Eq (7) can be useful in modelling real life survival data with different shapes of hrf. Table (1) lists −log[Z(x; φ)] and the associated parameters for some special distributions.

Table 1

Distributions and corresponding −log[Z(x; φ)] functions.

Distribution	−log[Z(x; φ)]	φ
Burr XII (x > 0)	−log[1 − (1 + xa)−b]	(a, b)
Weibull (x > 0)	−log[1−e^−ax^b]	(a, b)
Normal (−∞ < x < ∞)	−log[ϕ(x−μσ)]	(μ, σ)
Nadarajah Haghighi (x > 0)	−log[1−e^1−(1+λx)^α]	(α, λ)
Rayleigh (x > 0)	−log[1−e^−ax^2]	(a)
Exponential (x > 0)	−log[1−e^−αx]	(α)
Power function (0 < x < a)	−blog[xa]	(a, b)
Fréchet (x > 0)	(λx)^σ	(λ, σ)
Inverted Rayleigh (x > 0)	(λx^2)^σ	(λ, σ)
Burr III (x > 0)	z log[1 + x−c]	(c, z)
Pareto (δ < x < ∞)	−log[1 − (δ/x)λ]	(δ, λ)

Here three special models of EGuG family of distribution are defined.

2.1.1 EGu-Weibull(EGuW)

The EGu-W model is defined from Eq (7) by taking Z(x; φ) = 1 − exp{−ax}, z(x; φ) = abx exp{−ax}, as cdf and pdf of the baseline Weibull distribution with a, b > 0, respectively.

The cdf and pdf of EGu-W distribution are, respectively, given by and where θ, σ and b are shape parameters while a is scale parameter.

2.1.2 EGu-BurrXII(EGuBXII)

Let us consider the parent distribution as BXII with power parameters a, b > 0 by taking Z(x; φ) = 1−(1 + x)−, z(x; φ) = abx(1 + x)− be the cdf and pdf of the distribution.

The cdf and pdf of EGu-BXII distribution are, respectively, given by and where θ, σ, a and b are shape parameters.

2.1.3 EGu-Nadarajah Haghighi(EGuNH)

Consider to take Nadaraah Haghighi (NH) as baseline distribution with cdf as Z(x; φ) = [1 − exp{1−(1 + λx)}] and pdf as z(x; φ) = λβ(1 + λx)exp{1 − (1 + λx)}. Then, the cdf and pdf of EGuNH reduces to and where θ, σ and α are shape parameters while λ is scale parameter.

2.2 Useful expansion for the EGuG cdf

We provide a useful expansion for Eq (7) in terms of linear combinations of exp-G distribution. For a random baseline cdf Z(x), an av is said to have the exp-G distribution having parameter ℓ such that ℓ > 0, say Y ∼ exp-G (ℓ), if its pdf and cdf are given as respectively. Thus, several properties of the proposed model can be derived from those properties of the exp-G distribution studied by the authors in [3-8], to mention few.

By expanding Eq (7) using binomial and power series expansion, the resultant expression is given

Using Mathematica software, it can be verified that we can start the limit of integers (i, j) from 1 instead of 0 in above equation. Further, we can write [−log Z(x; φ)]− as since .

Now consider, for any real parameter c and zϵ(0, 1), the following formula holds: where P0(c) = 1/2; P1(c) = c(3c + 5)/24; P2(c) = c(c2 + 5c + 6)/48 etc is the stirling’s polynomial. Then, the cdf F(x) in Eq (13) can be expressed (using Eq (14) and generalized binomial expansion) as where .

By differentiating Eq (15), we obtain where h(x) = mZ(x; Φ) z(x; Φ) is the exp-G density function with power parameter m.

2.3 Shapes of density and hazard function

Analytical descriptions of density and hrf forms are conceivable. The roots of the equation represent the EGuG density’s critical points: The equation is used to find the EGuG hrf’s crucial points. Any numerical software can be used to examine Eqs (17) and (18) to determine the local maximum and minimum and inflexion points.

2.4 Moments

The first formula for the sth moment of X follows from Eq (16) as where Setting s = 1 in Eq (refrthmoment1pdfmix1) can provide explicit expression for the mean of several parent distributions.

A second alternative formula for is obtained from Eq (19) in terms of the baseline qf as where

The central moments (μ) and cumulants (κ) of X can follow from Eq (19) as and , respectively, where . The skewness and kurtosis can be calculated from the third and fourth standardized cumulants.

The sth incomplete moment of X can be determined from Eq (16) as The last integral can be computed for most G distributions.

A crucial applicability of the first incomplete moment m1(⋅) has to do with the Bonferroni and Lorenz curves, which are extremely beneficial in a variety of fields. For a given probability π, they are given by and , respectively, where m1(y) comes from Eq (refincompletepdfmix1) with s = 1 and q = Q(π) follows from Eq (9). The Lorenz and Bonferroni curve for EGuNH are displayed graphically (Figs 13 and 14, subsequently).

The totality of excursions from the mean and median is used to estimate the degree of scatter in a population and is defined by and , respectively, where is the mean and M = Q(0.5) is the median. These measures can be expressed as and , where is given by Eq (refcdfEGuG).

The moment generating function (mgf) of X can be expressed as where M(t) is the mgf of Y. Hence, M(t) can be determined from the exp-G generating functions.

2.5 Order statistics

Order statistics are used in a wide range of statistical theory and practise. Let X1, …, X is a random sample from the EGu-G distribution and X denote the ith order statistic. Then, pdf of X can be written as Inserting Eq (refcdfEGuG) and Eq (refpdfEGuG) in the last equation, and expanding it as in section (3.1), we get where and

2.6 Estimation

The three alternate approaches for inference are point estimation, interval estimation, and hypothesis tests. Several approaches for parameter point estimation have been published in the literature, the most extensively utilised of which is the maximum likelihood method. MLEs (maximum likelihood estimates) have properties that can be used to construct confidence ranges for model parameters. Large sample theory provides simple approximations that work well in repeated sampling for these estimations. The normal approximation for MLEs can be tackled analytically or computationally in distribution theory.

We use the optimum likelihood method to estimate the unknown parameters of the new distribution. Let x1, ⋯, x be n observations from the EGu-G family given by Eq (8) with parameter vector Θ = (θ, β; Φ)⊤. The log-likelihood ℑ = ℑ(Θ) for Θ is given by Eq (refmleegug) can be maximized either directly by using the R (optim function), SAS (NLMixed procedure) or Ox (MaxBFGS function), or then by solving the nonlinear likelihood equations by differentiating it. The components of the score vector U(Θ) are Setting these equations to zero and solving them simultaneously yields the MLEs of the family parameters.

The observed information matrix for the parameter vector Θ = (θ, σ, Φ)⊤ is given by whose elements can be determined by using any mathematical software. Under normal conditions of regularity, the multivariate normal distribution, where is the observed information analysed at , can be used to estimate confidence ranges for model parameters. Furthermore, we may use likelihood ratio (LR) statistics to assess the EGuG model to any of its specific models.

3 Properties of EGuNH

In comparison to Gamma, Weibull, and exponentiated exponential distributions, NH distribution (also known as extended exponential distribution) is the preferred option for zero inflated data. The cdf and pdf of NH distribution has already been defined in Section (3.1.3). For λ = 1, we define the cdf and pdf of EGuNH distribution as and where C = 1 − (1 + x) and θ, σ and α are shape parameters.

Henceforth, we denote by X a av having density (25). The sf and hrf of X has the form and

3.1 Shapes of density and hazard rate function of EGuNH distribution

The crucial points of the pdf of X are obtained from the equation: Similarly, the critical points of the hrf of X are obtained from the equation:

Some plots of the density of EGuNH for selected parameter values are presented in Figs 1–4 while plots of the hrf of EGuNH for random parameter values are presented (Figs 5–8). It is apparent that the density of EGuNH can be reversed-J, unimodal, and symmetrical. Similarly, EGuNH hazard rate shapes may tend to be increasing, decreasing, bathtub, or upside-down bathtub. The new model is much superior at fitting data sets in a variety of risk evaluation scenarios.

Fig 1

Plots of EGuNH density for a variety of parameter combinations.

Fig 4

Plots of EGuNH density for a variety of parameter combinations.

Fig 5

Plots of EGuNH hazard rate for a variety of parameter combinations.

Fig 8

Plots of EGuNH hazard rate for a variety of parameter combinations.

3.2 Central properties of EGuNH distribution

In this section, some useful expressions for the linear expansion, moments and incomplete moments of EGuNH distribution have been deduced using the Eq (16).

Proposition 1. where

. Recalling the result defined in Eq (15) as A straightforward differentiation of the above result yields density by The result in (27) is the linear expansion of NH densities. Hence, we shall derive several core properties of EGuNH using the major result of Eq 27.

Proposition 2.

Let W be a av with density ψ(x;m, α). Then, several properties of W can follow from those of X. The sth ordinary moment of X can be written as where and Utilizing the results derived in Eq (19), the sth moments are defined in (28).

Proposition 3.

The sth incomplete moment expression can be written as where Following the results defined in Eq (21), the sth incomplete moments are defined in (29). The skewness and kurtosis of X can be calculated from the third and fourth standardized cumulants. The classical skewness (Fig 9) and kurtosis plots (Fig 10) of the EGuNH distribution are displayed. Additionally, we provide the graphical illustration of MacGvillary skewness (MGs), which is based on quantile approach, in Figs 11 and 12. These plots reveal that the parameters θ and σ play a decisive role in modeling the skewness and kurtosis behaviors of X.

Fig 9

Bowley skewness of EGuNH.

Fig 10

Moors kurtosis of EGuNH.

Fig 11

MacGillivray skewness of EGuNH for a variety of parameter combinations.

Fig 12

MacGillivray skewness of EGuNH for a variety of parameter combinations.

MacGillivray [35] proposed another method to evaluate the skewness measure based on the qf and is defined as where uϵ(0, 1), Q(.) is the qf defined in (32).

The MG skewness plots are very sensitive for extremely small values of parameter θ and σ which certainly signifies longer tails of EGuNH. Likewise, plots of the Lorenz (Fig 13) and Bonferroni (Fig 14) curves of EGuNH distribution for some random values are displayed. These plots reveal how the distribution parameters affect inequality measures which can be used to establish some orderings, an essential feature for applied statisticians. Some descriptive statistics related to EGuNH are presented in Tables 2 and 3, respectively.

Fig 13

Plots of the Lorenz curves of EGuNH distribution.

Fig 14

Plots of the Bonferroni curves of EGuNH distribution.

Table 2

Descriptive measures of EGuNH for some parameter values.

Parameter values	Descriptives
(θ, σ, α)	Q1	Q2	Q3	B	M
(2.1, 0.1, 1.5)	0.275	0.294	0.311	−0.060	1.248
(2.1, 0.29, 1.5)	0.255	0.310	0.361	−0.041	1.231
(85.5, 0.979, 0.85)	1.990	2.182	2.395	0.053	1.241
(6.5, 0.33, 1.35)	0.39	0.436	0.482	−0.003	1.234
(1.5, 1.9, 2.15)	0.012	0.189	0.451	0.192	0.866
(0.91, 2.1, 1.5)	0	0.048	0.497	0.805	1.314
(1.7, 9.1, 1.5)	0	0.604	2.536	0.524	0.757
(0.7, 5.1, 1.3)	0	1.394	1.996	0.667	1.757
(1.7, 0.1, 0.85)	0	0.314	2.267	0.609	0.888

Table 3

Moments and moment ratios of EGuNH for some parameter combinations.

Parameter values	Moments and moments ratio
(θ, σ, α)	E(x)	E(x2)	E(x3)	E(x4)	V(x)	σ(x)	CV	CS	CK
(0.83, 2.1, 1.5)	0.225	0.233	0.241	0.299	0.141	0.401	1.322	1.589	4.253
(0.91, 2.1, 1.5)	0.295	0.263	0.291	0.365	0.176	0.420	1.422	1.680	4.944
(1.7, 9.1, 1.5)	1.343	4.258	15.684	63.478	2.456	1.567	1.167	1.680	4.944
(2.1, 0.1, 1.5)	0.292	0.086	0.026	0.008	0.001	0.027	0.094	2.505	6.661
(2.1, 0.29, 1.5)	0.306	0.100	0.034	0.012	0.006	0.078	0.254	-0.285	2.957
(1.5, 1.9, 2.15)	0.597	0.767	1.219	2.205	0.410	0.640	1.072	88.455	47.815
(6.5, 0.33, 1.35)	0.436	0.194	0.089	0.041	0.005	0.068	0.157	-0.044	3.206
(85.5, 0.979, 0.85)	2.163	4.764	10.965	24.462	0.088	0.297	0.137	0.437	3.289

3.3 Acturial measures EGuNH: Value at risk

The theory of finance is based upon risk evaluation. Investors are particularly interested to invest in entities in which there is minimum risk (specified with high degree of confidence) of losing money. In finance, value at risk (VaR) is the most extensively used metric for assessing liability. It is also known as quantile risk measure or quantile premium principle of the distribution of aggregate losses. It is characterised by a level of assurance “q” (usually at 95% and 99%). To a layman, VaR answers a simple question that “What is the worst case scenario that can happen in a particular investment?”

If X has pdf (25), then VaR is the qth quantile of its cdf (24), defined as

3.4 Acturial measures EGuNH: Expected shortfall

Despite of the popularity of VaR measures, there are many shortcomings (see [36]). To counter inherent problems in VaR, Artzner et al. [37, 38] proposed the use of expected shortfall (ES). Expected shortfall quantifies the average loss in states beyond the VaR level. ES has a number of aliases such as “conditional VaR”, “mean excess loss” or “tail VaR”. We define the ES as follows For a combination of various parameter values, plots of VaRs (Fig 15) and ESs (Fig 16) are displayed respectively.

Fig 15

Plots of the VaR of EGuNH distribution for some random parameter values.

Fig 16

Plots of the ES of EGuNH distribution for some random parameter values.

3.5 Parameter estimation of EGuNH

The log-likelihood function ℑ for the vector of parameters Θ = (θ, σ, α)⊤ for the model defined in (25) is given by The components of the score vector U(Θ) are The MLE of Θ can also be obtained by solving the nonlinear equations U = 0, U = 0 and U = 0. Because these equations cannot be solved analytically, the estimates can be calculated numerically using statistical software.

3.6 Simulation study of EGuNH distribution

The qf of the EGuNH distribution has an explicit form as follows Here, we use Monte Carlo simulations to demonstrate the performance and correctness of maximum likelihood estimations of the EGuNH parameters by inverting Eq (24) to generate a sample data from the model. The simulation study is perform for sample sizes n = 50, 100, 200, 500, and parameter combinations: I: θ = 0.2, σ = 0.75 and α = 0.5, II: θ = 2.2, σ = 0.45 and α = 0.5, III: θ = 3.4, σ = 0.75 and α = 1.5 and IV: θ = 3.4, σ = 1.35 and α = 1.5. This study is carried out for N = 2000 times, each with given n and computed the average estimates (AEs) as well as their average biases (Bias), mean squared errors (MSEs) and coverage probabilities (CPs) of the MLEs. The AEs, Bias, MSEs and CPs for the parameters θ, σ and α are given in Tables (4)–(7). The empirical findings suggests that the bias and MSEs decreases as sample size increases. Further, the empirical CPs are quite close to the nominal level of 95%. As a result, MLEs and their approximate findings can be used to evaluate and build approximated confidence intervals of the EGuNH distribution parameters θ, σ and α.

Table 4

AEs, Biases, MSEs and CPs for combination-I.

		n = 50			n = 100
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	0.279	0.081	0.447	0.99	0.274	0.074	0.436	0.98
σ	0.927	-1.661	0.691	1.00	0.895	-1.651	0.683	0.97
α	0.925	-0.241	0.494	0.98	0.817	-0.203	0.441	0.92
		n = 200			n = 500
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	0.272	0.072	0.429	0.96	0.272	0.051	0.297	0.94
σ	0.784	-1.350	0.669	0.95	0.757	-1.115	0.662	0.95
α	0.629	-0.220	0.359	0.95	0.548	-0.204	0.320	0.96

Table 7

AEs, Biases, MSEs and CPs for combination-IV.

		n = 50			n = 100
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	3.877	-0.281	0.097	0.97	3.595	-0.223	0.019	0.92
σ	0.889	0.021	0.025	1.00	0.827	0.010	0.020	0.94
α	2.541	0.058	0.033	0.97	2.230	0.030	0.028	0.99
		n = 200			n = 500
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	3.495	-0.019	0.010	0.96	3.410	-0.013	0.009	0.95
σ	0.765	0.005	0.016	0.97	0.751	0.003	0.008	0.96
α	1.915	0.015	0.014	0.95	1.507	0.011	0.010	0.95

4 Applications of the EGuNH distribution

Statistical methods that fail to account for all of the uncertainties in the model are prone to produce an overly optimistic assessment of future extremes, are frequently contradicted by observed extreme events in a variety of scientific fields. The current literature regarding extreme value theory is full of such models in which data sets are meteorology data such as earthquakes, floods, rains, droughts, hurricanes etc. (see [15-32]). On the contrary, health hazards is an area of extreme value theory which should be explored. Death, damage, or disease; exacerbation of underlying medical disorders; and negative effects on mental health are some of the health hazards of climate-related increases in exposure to extreme occurrences.

In this section, we provide some applications of the EGuNH model on three real life phenomenons, two of which related to health hazards in extreme value theory. We estimate the unknown parameters of the distributions by the principal of maximum likelihood. We compute the log-likelihood function evaluated at the MLEs () using the method of a limited-memory quasi-Newton code for bound-constrained maximization (L-BFGS-B). In order to select the best probability model, a variety of criteria for evaluating information (ICs) can be considered. We considered the following well-known ICs: the maximized log-likelihood (), Akaike Information criterion (AIC), Anderson-Darling (A⋆), Cramér-von Mises (W⋆) and Kolmogorov-Smirnov measures (D⋆; P-value (p⋆)), where lower values of all these statistics except higher p⋆ values of K-S, indicate good fits. The required computations are carried out using the R script AdequacyModel which is freely accessible from http://cran.r-project.org/web/packages/AdequacyModel/AdequacyModel.pdf.

The fits of the EGuNH distribution is compared with other competitive models which are given in Table 8. The parameters are all positive real numbers of these densities.

Table 8

The comparative fitted models.

Distribution	Author(s)
GaNH	Cordiero et al., (2015) [31]
LxNH	Ramirez et al., (2020) [39]
TLNH	Yuwadee Sangsanit and Winai Bodhisuwan, (2016) [40]
ENH	Lemonte et al., (2013) [41]
MONH	Lemonte et al., (2016) [42]
NH	Nadarajah and Haghighi, (2011) [43]

4.1 Meteorology data

Meteorological phenomena are weather events that most individuals are affected by, due to changes in extreme weather and climatic events, such as earthquakes, heat waves, floods, hurricanes, droughts etc. The present data is taken from [44], denoted by D1, gives the time in days between successive serious earthquakes world-wide. An earthquake is included if its magnitude was at least 7.5 on the Richter scale, or if over 1000 people were killed. There were 63 earthquakes recorded altogether, and so 62 recorded waiting times. The data are: 840, 157, 145, 44, 33, 121, 150, 280, 434, 736, 584, 887, 263, 1901, 695, 294, 562, 721,40, 1336, 335, 1354, 454,139, 780, 203, 436, 30, 246, 1617, 638, 937, 735,76, 710, 36, 667,384, 129, 46, 402, 194, 40, 556, 99, 9, 209, 599, 38, 365, 92, 82, 220, 759, 304, 83, 319, 375, 832, 460, 567, 328.

4.2 Cancer data

According to [45], extreme events have the potential to disrupt the delivery of cancer care. For example, some deadly carcinogens may be released into communities as a result of hurricanes and wild fires; industry shutdowns may result in a shortage of life-saving medical equipment in hospitals, causing shortages in cancer facilities across the country; and infrastructure collapse may limit access to patients undergoing cancer therapies. The following two data, denoted by cancer 1 (D2) and cancer 2 (D3) are related to cancer patients.

Cancer data 1. The survival times, in weeks, of 33 patients who succumbed to Acute Myelogenous Leukemia are the subject of D2. This data was recently studied by the authors in [46] The data are: 65, 156, 100, 134, 16, 108, 121, 4,39, 143, 56, 26, 22, 1, 1, 5, 65, 56, 65, 17, 7, 16, 22, 3, 4, 2, 8, 4, 3, 30, 4,43.

Cancer data 2. D3 signifies the number of patients suffering from blood cancer. The Saudi Cancer Registry (SCR) provides such information, covering the time period from 1994 to the present day. The data is extracted from a report [47] which concerns an overview of cancer incidence statistics for Saudi Arabia in 2012. The data are: 1277, 1290, 1357, 1369, 1408, 1455, 1478, 1549, 115, 181, 255, 418, 441, 461, 516, 739, 743, 789, 807, 865, 924, 983, 1024, 1062, 1063, 1165, 1191, 1222, 1222, 1251, 1578, 1578, 1599, 1603, 1605, 1696, 1735, 1799, 1815, 1852. The descriptive statistics for each of the three data sets are given in Table 9.

Table 9

The descriptive statistics related to D1, D2 and D3.

Data	Sample Size	Arithmetic Mean	Standard Deviation	Lowest	Highest	Skewness	Kurtosis
1	62	437.21	399.93	9	1901	1.50	2.52
2	32	42.07	46.95	1	156	1.12	0.03
3	40	1137	481.60	115	1852	-0.49	-0.73

The empirical findings of all the three data are suggestive of the heavy tailed data. The TTT plots Figs (17)–(19) for the data sets are given. In particular, the TTT plots show bathtub, increasing and decreasing hrf, allowing us to fit EGuNH model on these data sets. The approximated hrf Figs (20)–(22) for each data point correlates to the TTT graphs. Table 10 summarizes the results of the MLEs and their related standard errors (in parentheses) of the model parameters for the proposed model while the ICs are listed in Table 11 for the D1, D2 and D3, respectively. It is customary to supplement the analytical result defined in Tables 9 and 10, by displaying it graphically. Hence, the estimated pdfs Figs (23)–(25), PP–plots Figs (26)–(28), estimated cdfs Figs (29)–(31) and estimated sfs Figs (32)–(34) for the three data sets are given. On the given data sets, the numerical values authenticates that the EGuNH model provides the best fit as compared to the other models.

Fig 17

TTT plot for D1.

Fig 19

TTT plot for D3.

Fig 20

Estimated hrf for D1.

Fig 22

Estimated hrf plots for D3.

Table 10

MLEs with their respective SEs (in parenthesis) for D1, D2 and D3.

	Data 1	Data 2	Data 3
Distribution	MLEs	MLEs	MLEs
EGuNH	3.655,11.787,0.275	0.071,0.031,0.076	22.753,20.069,0.469
(θ,σ,α)	(2.217),(2.371),(0.073)	(0.015),(0.003),(0.001)	(5.093),(2.193),(0.011)
ENH	1.285,0.764,0.004	1.256,1.716,0.011	4.515,0.971,0.002
(θ,α,λ)	(0.229),(0.117),(0.001)	(0.178),(0.556),(0.005)	(0.987),(0.066),(0.002)
GaNH	0.295,7.164,3.874	1.6,1.256,0.013	1.146,3.802,0.002
(θ,aα,λ)	(1.018),(6.192),(2.987)	(0.542),(0.174),(0.007)	(0.094),(0.615),(0.518)
LxNH	8.918,0.080,18.577	1.470,3.660,0.006	13.499,0.111,0.456
(θ,α,λ)	(3.211),(0.027),(58.410)	(0.126),(0.930),(0.002)	(3.672),(0.025),(0.643)
TLNH	1.499,0.600,0.003	1.374,1.5490.007	11.215,0.371,0.009
(θ,α,λ)	(0.262),(0.084),(0.987)	(0.173),(0.405),(0.003)	(6.699),(0.072),(0.017)
MONH	5.162,0.886,0.003	3.067,0.962,0.038	42.021,0.827,0.005
(θ,α,λ)	(0.559),(0.131),(0.001)	(0.159),(0.560),(0.055)	(2.127),(0.081),(0.311)
NH	1.054,0.212	2.616,0.006	0.067,158.079
(α,λ)	(0.104),(0.298)	(0.728),(0.012)	(0.095),(89.711)

Table 11

The statistics , AIC, BIC, A⋆, W⋆, D⋆ and p⋆ for D1, D2 and D3.

Distribution	−ℓ^	AIC	BIC	A ⋆	W ⋆	D ⋆	p ⋆
Data set 1
EGuNH	441.19	889.90	896.33	0.41	0.05	0.08	0.79
GaNH	444.78	895.56	901.94	1.05	0.17	0.11	0.54
LxNH	448.86	903.71	910.14	1.58	0.25	0.18	0.45
TLNH	442.73	891.24	897.67	0.61	0.09	0.10	0.68
ENH	442.67	890.98	897.74	0.46	0.07	0.09	0.77
MONH	441.90	890.14	897.01	0.34	0.05	0.09	0.78
NH	442.27	891.75	897.15	0.31	0.06	0.08	0.78
Data set 2
EGuNH	149.027	304.054	308.451	0.567	0.08	0.111	0.817
GaNH	150.203	306.406	310.803	0.588	0.09	0.122	0.729
LxNH	152.672	311.345	315.741	0.860	0.14	0.126	0.687
TLNH	150.104	306.205	310.602	0.565	0.08	0.127	0.681
ENH	150.101	306.202	310.599	0.581	0.08	0.126	0.692
MONH	150.036	306.071	310.468	0.593	0.09	0.143	0.534
NH	152.875	310.751	312.005	0.636	0.10	0.159	0.487
Data set 3
EGuNH	307.586	621.173	626.240	1.401	0.226	0.151	0.317
ENH	310.775	627.566	632.616	1.912	0.319	0.199	0.084
GaNH	308.249	622.499	627.566	1.494	0.243	0.155	0.284
LxNH	313.300	632.601	637.668	2.176	0.369	0.172	0.183
TLNH	316.845	639.691	644.757	2.868	0.500	0.184	0.131
MONH	307.698	622.019	627.107	1.626	0.235	0.161	0.289
NH	400.885	805.770	809.184	2.518	0.433	0.607	0.000

Fig 23

Estimated density for D1.

Fig 25

Estimated density plot for D3.

Fig 26

PP plot for D1.

Fig 28

PP plot for D3.

Fig 29

Estimated cdf plot for D1.

Fig 31

Estimated cdf plot for D3.

Fig 32

Estimated sf plot for D1.

Fig 34

Estimated sf plot for D3.

The variance-covariance matrices of the MLEs of the EGuNH distribution for D1 is The variance-covariance matrices of the MLEs of the EGuNH distribution for D2 is The variance-covariance matrices of the MLEs of the EGuNH distribution for D3 is

4.2.1 Numerical calculations of VaRs and ESs

We were able to further investigate EGuNH’s application to these risk measures thanks to the results reported in Section 4. To quantify the volatility associated with these measures, we take the values of MLEs of D1, D2 and D3., respectively, from Table 11. Higher risk measures indicate heavier tails, while lower risk measures indicate a model with a much lighter tail tendency. It’s pertinent to mention that the EGuNH model yielded significantly more impressive results than others, implying that the model has a longer tail. The numerical findings of VaRs and ESs for data 1, data 2, and data 3 of the proposed model at respective level of significance (LoS) are shown in Table 12. The summarized output of these risk measures (VaRs in Figs 35–37 and ESs in Figs 38–40), graphically, for the reader’s expedience.

Table 12

Numerical measures of VaRs and ESs of EGuNH for D1, D2 and D3.

LoS	Data 1		Data 2		Data 3
	VaRs	ESs	VaRs	ESs	VaRs	ESs
0.55	327.4119	146.4568	37.86083	12.29339	1125.354	727.9202
0.60	373.4939	163.4189	45.77802	14.74769	1201.044	764.1592
0.65	425.7544	181.5435	54.61479	17.46849	1282.559	800.8571
0.70	486.2076	201.0896	64.41790	20.46624	1372.164	838.4086
0.75	557.9849	222.4155	75.25685	23.75171	1473.216	877.2859
0.80	646.3563	246.0443	87.25890	27.33904	1591.174	918.1102
0.85	761.2607	272.8042	100.69662	31.25081	1736.442	961.7964
0.90	925.2607	304.1764	116.24233	35.52842	1927.397	1009.9029
0.95	1210.5273	343.4825	136.04737	40.26890	2210.153	1065.6140

Fig 35

Estimated VaRs for D1.

Fig 37

Estimated VaRs for D3.

Fig 38

Estimated ES for D1.

Fig 40

Estimated ES for D3.

5 Concluding remarks

We propose and study the EGuG model and obtain a wide range of mathematical and statistical modelling methods to characterise the model’s structural and dynamic aspects including properties such as quantile function, ordinary and incomplete moments, mean deviations, bonferroni and lorenz curves, generating function and order statistics. The parameters of the family are estimated by the method of maximum likelihood. An extended exponential distribution is taken as baseline model to propose EGuNH distribution. Some simulations are performed to check the asymptotic properties of the estimates. Three applications to real data set are presented to illustrate the potentiality of the proposed models. For future research, the proposed model can further be extended using compounding. We expect that the modification may facilitate in estimating analytically tractable Bayesian estimates of the reliability function under different priors.

9 Feb 2022

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. Please include in your Methods section further details of how the dataset was obtained for the study, and provide the data sources.

3. Please ensure that you refer to Figures 4 and 6 in your text as, if accepted, production will need this reference to link the reader to the figure.

4. We note you have included a table to which you do not refer in the text of your manuscript. Please ensure that you refer to Table 7 in your text; if accepted, production will need this reference to link the reader to the Table.

Additional Editor Comments:

The manuscript were reviewed by three experts, and their comments are included. Please carefully address all the comments, especially the comment on the originality of my paper.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: No

Reviewer #3: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: No

Reviewer #3: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: This paper proposed a new generalization of the Gumbel distribution as well as three separate models evolved in the same family of distribution. As the author claimed, this should be the first time that the Nadarajah-Haghighi(NH) distribution is chosen as the baseline model to apply to the EGuG family. They formed the EGuNH distribution, conducted completed simulation study and applied it to lift data sets.

In all, this is a well-organized and written work, the idea and method are innovative and thoughtful. The reviewer would suggest the following modifications:

1. When claiming the start of limit of integers could be changed in formula (11), the author need to clarify why it could be verified by Mathematica, maybe list the result given by Mathematica, to make the point more convincing.

2. In section 3.6, there are 4 groups of parameter combinations but only 3 out of them are tested, the author should explain why they chose not to test all of them and explain how to decide which one to be tested.

3. The author should make some explanation or do some analysis for figure 5 instead of just leaving it there

4. At the beginning of section 3, the citation “The cdf and pdf…” should be Section 2.1.3 instead of 1.1.3

5. There are many typesetting errors and some punctuation errors, the author need carefully check the paper.

Reviewer #2: Thank you for the opportunity to review the manuscript. This work builds on prior work on Gumbel distribution. The authors studied the exponentiated Gumbel-G (EGuG) model and proposed a new model (called the EGuNH model) by extending the EGuG model. The new model takes the Nadarajah-Haghighi distribution as the baseline model. The simulations and real data are used to examine the model properties. The proofs in the paper seem to be mathematically correct. However, I have some concerns about the paper I would like to discuss with the authors.

Major concerns relating to methods or significance:

This work is very similar to A. A. Ogunde et.al (https://www.hindawi.com/journals/jam/2020/2798327/ ) and Hormatollah Pourreza et.al (https://journals.sagepub.com/doi/full/10.1177/09622802211009262 ). Did authors read these papers to compare the result? In addition, the manuscript is not clearly presented overall. It seems to me that this paper did not add much to the scientific value in terms of combined level of methodological and practical innovation.

Minor concerns related to clarity:

1. In the Abstract, the authors used the term “greatest likelihood method”. Is it the maximum likelihood method? If so, I think the common term should be the maximum likelihood method.

2. In Section 1, the authors did not explain what T-X methodology is. It would be good if the authors can provide a citation.

3. In Section 2, it would be good to use the hazard rate function to define “hrf” rather than failure rate function.

4. In-Page 13, the authors did not cite figure # for “Plots of MGs for some parameter values”.

5. The tables should be self-explanatory. Please explain the acronyms in the tables. For example, what is LoS in Table 12?

6. The authors seem did not explain most of figures and tables. For example, what is the message for Figure 13?

7. In the simulations, it would be good if the authors could report the coverage probability of the 95% confidence intervals as another evaluation metric for MLE estimation of model parameters.

8. In real data examples, the authors did not show how to interpret model parameters and their practical meanings in real data.

Reviewer #3: In the manuscript “A New Extended Gumbel distribution: Properties and Application”, Fayomi et al. introduced a new generalization of the Gumbel distribution. The author presented the distribution through the data simulation as well as three real-world applications.

The paper is very well presented, and the author provided a good review of the literature. The equations are well explained and clearly presented; I am glad to see that the author used the real-world applications to show the superiority of the distribution.

I would have like to see how the extended Gumbel distribution behaves for global sequence alignment application. For example, Sardiu et al. 2005, showed that the score statistics of sequence alignment follows a Gumbel distribution. Is the extended Gumbel distribution presented in this paper superior to the classical Tracy-Widom distribution for example for global alignments? I suggest the author cite this paper since the application of the Gumbel distribution in biological applications were used since 2005. The author does not talk in the manuscript about the Tracy-Widom distributions. This needs to be included in the introduction.

I would also suggest that the author comment more on the limitation of this distribution. How the distribution change when the size is large or too small, for example. Is it a critical point where the distribution does not apply anymore? These are important points which need to be addressed to prove superiority and utility of the extended Gumbel distribution.

Overall, the manuscript is good, however I suggest that the author address these concerns.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes

Reviewer #2: No

Reviewer #3: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

Submitted filename: Plos_review_2022.docx

Click here for additional data file.

12 Mar 2022

It was the evident wish and will of the Reviewers that we make the changes to our manuscript on their insightful observations. Without losing the originality, we have tried to accommodate their suggestions in best possible manner. We have also improved the overall presentation of the paper.

We thank the Editor and the three referees for the constructive comments and hope that the revision is now appropriate for publishing.

Submitted filename: Response to Reviewers.pdf

Click here for additional data file.

4 Apr 2022

A New Extended Gumbel distribution: Properties

and Application

PONE-D-22-00530R1

Dear Dr. Khan,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

The revised manuscript was reviewed by two experts. Both of them believe that all comments have been addressed. I agree with them and suggest acceptance.

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #3: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #3: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #3: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #3: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #3: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: (No Response)

Reviewer #3: All concerns have been addressed. The author addressed my questions. No need for further work. I suggest this paper for acceptance.

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #3: No

18 Apr 2022

PONE-D-22-00530R1

A New Extended Gumbel Distribution: Properties And Application

Dear Dr. Khan:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Jiangtao Gou

Academic Editor

PLOS ONE

Table 5

AEs, Biases, MSEs and CPs for combination-II.

		n = 50			n = 100
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	2.307	-0.031	0.027	1.00	2.197	-0.023	0.012	0.92
σ	0.499	-0.021	0.011	0.91	0.495	-0.005	0.002	0.98
α	0.511	0.018	0.018	0.99	0.510	0.010	0.005	0.93
		n = 200			n = 500
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	2.196	-0.014	0.002	0.97	2.198	-0.011	0.001	0.95
σ	0.501	0.001	0.001	0.95	0.498	0.001	0.001	0.94
α	0.508	0.008	0.003	0.94	0.501	0.000	0.002	0.95

Table 6

AEs, Biases, MSEs and CPs for combination-IV.

		n = 50			n = 100
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	3.507	-0.201	0.027	0.98	3.405	-0.123	0.019	0.93
σ	0.769	0.021	0.015	0.97	0.761	0.010	0.013	0.99
α	2.541	0.058	0.023	0.90	2.130	0.030	0.018	0.97
		n = 200			n = 500
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
	AEs	Bias	MSEs	CPs	AEs	Bias	MSEs	CPs
θ	3.405	-0.014	0.010	0.96	3.401	-0.011	0.007	0.96
σ	0.755	0.005	0.006	0.94	0.753	0.003	0.004	0.95
α	1.615	0.015	0.004	0.96	1.511	0.011	0.002	0.94