Next Article in Journal
Federated Learning Incentive Mechanism with Supervised Fuzzy Shapley Value
Next Article in Special Issue
Transfer Learning for Logistic Regression with Differential Privacy
Previous Article in Journal
Vector-Valued Fuzzy Metric Spaces and Fixed Point Theorems
Previous Article in Special Issue
Randomly Stopped Minimum, Maximum, Minimum of Sums and Maximum of Sums with Generalized Subexponential Distributions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Extension of the Fréchet Distribution and Applications

by
Yolanda M. Gómez
1,*,
Inmaculada Barranco-Chamorro
2,*,
Jaime S. Castillo
3 and
Héctor W. Gómez
3
1
Departamento de Estadística, Facultad de Ciencias, Universidad del Bío-Bío, Concepción 4081112, Chile
2
Departamento de Estadística e Investigación Operativa, Facultad de Matemáticas, Universidad de Sevilla, 41012 Sevilla, Spain
3
Departamento de Estadística y Ciencias de Datos, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta 1240000, Chile
*
Authors to whom correspondence should be addressed.
Axioms 2024, 13(4), 253; https://doi.org/10.3390/axioms13040253
Submission received: 21 February 2024 / Revised: 28 March 2024 / Accepted: 9 April 2024 / Published: 11 April 2024
(This article belongs to the Special Issue Probability, Statistics and Estimation)

Abstract

:
This paper presents the Slash-Exponential-Fréchet distribution, which is an expanded version of the Fréchet distribution. Through its stochastic representation, probability distribution function, moments and other relevant features are obtained. Evidence supports that the updated model displays a lighter right tail than the Fréchet model and is more flexible as for skewness and kurtosis. Results on maximum likelihood estimators are given. Our proposition’s applicability is demonstrated through a simulation study and the evaluation of two real-world datasets.

1. Introduction

The French mathematician Maurice René Fréchet introduced the Fréchet (Fr) distribution in the 1920s as a maximum value distribution, see Fréchet [1] or Fisher and Tippet [2]. It is also recognized as the inverse Weibull distribution and this is a specific instance of what is commonly referred to as Generalized Extreme Value (GEV) distributions, (see Gumbel [3] Embrechts et al. [4], Resnick [5], or Haan and Ferreira [6]). Since its origin, this model has been applied in different areas, such as reliability, life testing, extreme events, rainfall, wind speeds, among others. Information regarding the theoretical characteristics and uses of Fréchet distribution can be found in the works of Kotz and Nadarajah [7], Gupta et al. [8], Coles [9], and Ramos et al. [10]. Results using the term inverse Weibull distribution can be seen in Calabria and Pulcini [11], Maswadah [12], and Salman [13], and from a Bayesian point of view in Abbas and Tang [14]. These papers along with the references appearing therein, show the interest of Fréchet model. In this paper a new extension of the Fréchet distribution is proposed, which is based on the Fréchet and the Slash-Exponential distribution recently introduced in Punathumparambath [15]. Next, properties of these models are given, which will allow us to reach our end. So, a continuous random variable (rv) X follows a Fréchet distribution if its probability density function (pdf) is
f X ( x ; λ ) = λ x ( 1 + λ ) exp x λ , x > 0 ,
where λ > 0 is a shape parameter. Equation (1) is denoted X F r ( λ ) . Some properties of interest are listed in next lemma.
Lemma 1.
(a) 
Let X F r ( λ ) . Then its cumulative distribution function (cdf) is
F X ( x ; λ ) = exp x λ , x > 0 .
(b) 
Let V E x p ( 1 ) . Then V 1 λ F r ( λ ) .
(c) 
If X F r ( λ ) , then the rth-moment of X exists for r < λ and
E ( X r ) = Γ 1 r λ ,
where Γ ( · ) denotes the gamma function.
(d) 
Let X F r ( λ ) . Then the pdf of X is unimodal with mode m = λ 1 + λ 1 / λ or decreasing (for λ close to zero), the hazard function is always unimodal, see [10].
The gamma distribution will be used throughout this paper. In Remark 1, its pdf and cdf are given.
Remark 1.
If a rv T adheres to a gamma distribution, represented as T G a ( a , b ) , it implies that its pdf is as such
g ( t ; a , b ) = b a Γ ( a ) t a 1 e b t , t > 0 ,
with a , b > 0 . The cdf of T is
G ( z ; a , b ) = 0 z g ( t ; a , b ) d t = γ ( a , b z ) Γ ( a ) ,
where γ ( a , b z ) = 0 b z u a 1 e u d u is the (lower) incomplete gamma function, (Abramowitz and Stegun [16]).
Next, the basis of our proposal is introduced, that is, the Slash-Exponential distribution. As origin of the Slash-Exponential model, the Generalized Exponential distribution can be cited (Gupta and Kundu [17,18]). Astorga et al. [19] proposed an extension called Slashed Generalized Exponential, which contains as a submodel, the Slash Exponential. Details can be seen in Punathumparambath [15].
Definition 1.
A rv Y follows a Slash-Exponential distribution with shape parameter α > 0 , Y S E ( α ) , if its pdf is presented by
f Y ( y ; α ) = α 2 y ( 1 + α ) Γ ( α ) G ( y ; 1 + α , 1 ) , y > 0 ,
where G ( · ) denotes the cdf of the gamma distribution given in (5).
The following lemma presents some properties of Slash-Exponential distribution.
Lemma 2.
(a) 
Let Y S E ( α ) . Then the cdf is
F Y ( y ; α ) = 1 α y α Γ ( α ) G ( y ; 1 + α , 1 ) exp ( y ) , y > 0 .
(b) 
If Y S E ( α ) then the survival and hazard rate functions are as follows
S Y ( y ; α ) = α y α Γ ( α ) G ( y ; 1 + α , 1 ) + exp ( y ) , y > 0 ,
h Y ( y ; α ) = α 2 y ( 1 + α ) Γ ( α ) G ( y ; 1 + α , 1 ) α y α Γ ( α ) G ( y ; 1 + α , 1 ) + exp ( y ) , y > 0 .
(c) 
Let Y S E ( α ) . Then the rth-moment of Y exists for α > r and
E ( Y r ) = α α r r ! .
The objective of this article is to present a novel expansion of the Fréchet distribution taking as starting point the Slash-Exponential distribution defined in (6). Specifically, by proceeding similarly to property b given in Lemma 1, an updated extension of the Fréchet distribution has been achieved. This is called Slash-Exponential-Fréchet (SEFr) model. The SEFr distribution can be used in different fields where the Fréchet distribution is of common use, with the advantage of being a model more flexible as for skewness and kurtosis than the Fréchet model. On the other hand, the SFEr model can be used as an alternative to the Slash Fréchet distribution proposed in Castillo et al. [20] (model with two shape parameters), since the proposed model incorporates a scale parameter, resulting a distribution with greater flexibility. As novelty of this paper, we highlight that, the new model has a lighter right tail than the baseline Fréchet distribution, this is a new feature in the field of applications of slash methodology. The Slash distribution is statistically characterized as the ratio formed by two independent rv’s, with one conforming to a standard normal distribution and the other representing a power of a uniform distribution. Therefore, it can be mentioned that Y possesses a slash distribution if it can be represented as such:
Y = X 1 X 2
where X 1 N ( 0 , 1 ) and X 2 B e t a ( q , 1 ) , X 1 is independent of X 2 and q > 0 . In light of the groundbreaking studies by Rogers and Tukey [21], Andrews et al. [22] on slash distribution; multivariate slash models proposed in Gómez et al. [23] or Arslan and Genc [24], slash methodology has been proven to be useful to increase the weight of tails of a baseline distribution. See for instance Reyes et al. [25], Del Castillo [26], Zörnig [27], Olmos et al. [28] for the half-normal and generalized half-normal, Astorga et al. [19] for the generalized exponential distribution, Barranco-Chamorro et al. [29] for the Rayleigh distribution, Barrios et al. [30] for the power half-normal, Gui [31] for the Lindley, and Castillo et al. [32] for weighted Lindley, among others. However, in this paper, the slash methodology is applied in such a way, that we get a class of distributions more flexible than Fréchet baseline distribution as for skewness and kurtosis.
The structure of this paper is as follows. In Section 2, the stochastic representation and relevant properties are given. These are: pdf, mode, cdf, survival and hazard rate function, right tail behaviour, approximations to Fréchet model, moments, Shannon entropy, and other properties of interest in reliability. Section 3 is devoted to inferential results. Maximum Likelihood (ML) estimation method is studied in detail. A simulation study is conducted to evaluate the effective performance of ML estimators. In Section 4, a pair of practical applications are presented where our suggestion is benchmarked against other rival models, demonstrating its superior efficacy. The concluding remarks can be found in Section 5.

2. Results for the SEFr Distribution

First, the stochastic representation of the Slash-Exponential-Fréchet distribution is given as
X = σ Y 1 / λ ,
where Y S E ( α ) with σ , λ , α > 0 . Equation (7) is denoted by X S E F r ( σ , λ , α ) .
In Equation (7), σ > 0 is a scale parameter, whereas λ > 0 and α > 0 are shape parameters.
Next proposition provides the pdf of the SEFr model.
Proposition 1.
Let X S E F r ( σ , λ , α ) . Then the pdf of X is given by
f X ( x ; σ , λ , α ) = λ α 2 Γ ( α ) σ λ α x λ α 1 G σ x λ ; 1 + α , 1 , x > 0 ,
where σ , λ , α > 0 and G ( z ; a , b ) denotes the cdf of the gamma distribution introduced in (5).
Proof. 
By employing the representation outlined in (7) and the rv transformation method, the result is obtained.    □
Given that σ is a scale parameter, we can assume σ = 1 without any loss of generality. Figure 1 provides illustrative plots for varying values of α and λ . These plots suggest that, like the Fréchet model, the SEFr distribution is unimodal or decreasing. These appreciations are proved in Proposition 2.
Proposition 2.
Let X S E F r ( 1 , λ , α ) . Then the mode of X is given as the solution of
e x λ x λ ( α + 1 ) γ ( α + 1 , x λ ) = λ λ α 1 , p r o v i d e d t h a t λ α > 1 .
For λ α 1 , the pdf of X is strictly decreasing.
Proof. 
Straightforward, studying the first derivative of (8) with respect to x.    □

2.1. Properties

Proposition 3.
Let X S E F r ( σ , λ , α ) . Then, the cdf of X is provided by
F X ( x ; σ , λ , α ) = α Γ ( α ) x σ λ α G σ x λ ; 1 + α , 1 + exp σ x λ , x > 0 .
Proof. 
By definition
F X ( x ; σ , λ , α ) = λ α 2 Γ ( α ) σ λ α 0 x t λ α 1 G σ t λ , 1 + α , 1 d t .
By employing the methods of integration by parts, it can be shown that
u = G σ t λ ; 1 + α , 1 d u = λ σ λ t λ + 1 g σ t λ , 1 + α , 1 d t = λ σ λ t λ + 1 Γ ( 1 + α ) σ t λ α exp σ t λ d t , d v = t λ α 1 v = t λ α λ α ,
where g ( · ) denotes the pdf of the gamma distribution. Then,
F X ( x ; σ , λ , α ) = λ α 2 Γ ( α ) σ λ α × G σ t λ ; 1 + α , 1 x λ α λ α + σ λ + λ α α Γ ( 1 + α ) 0 x t λ 1 exp σ t λ d t ,
and making the substitution w = σ t λ , the proposed result is obtained.    □
Corollary 1.
Let X S E F r ( σ , λ , α ) . Then the survival function and hazard rate function of X are provided by
S ( x ) = 1 α Γ ( α ) x σ λ α G σ x λ ; 1 + α , 1 exp σ x λ , h ( x ) = λ α 2 Γ ( α ) x λ α 1 G σ x λ ; 1 + α , 1 σ λ α 1 α Γ ( α ) x σ λ α G σ x λ ; 1 + α , 1 exp σ x λ .
Proof. 
It follows by applying S ( x ) = 1 F X ( x ) and h ( x ) = f X ( x ) / S ( x ) .    □
As an illustration, plots for the cdf, survival function, and hazard rate function in the SEFr ( 1 , λ , 0.5 ) model are given in Figure 2.
Next we focus on studying the right tail of SEFr model. Specifically, in Proposition 4, it is proven that the right tail of SEFr model is lighter than the right tail of baseline Fréchet model. Moreover, in Proposition 5, it is proven that the SEFr model tends to a Fréchet model when α .
Proposition 4.
Let X S E F r ( 1 , λ , α ) and F F r ( λ ) with survival functions S X and S F , respectively. Then
S X ( t ) S F ( t ) , t > 0 .
That is, for all fixed α > 0 the right tail of S E F r ( 1 , λ , α ) model is lighter than the right tail in the F r ( λ ) distribution.
Proof. 
From Corollary 1 and (10), note that t > 0
S X ( t ) = 1 exp t λ α Γ ( α ) t λ α G ( t λ ; 1 + α , 1 ) = S F ( t ) α Γ ( α ) t λ α G ( t λ ; 1 + α , 1 ) .
Since α > 0 and G ( · ) is the cdf of a gamma distribution, we have that S X ( t ) S F ( t ) , t > 0 .    □
For λ = 5 and increasing values of α , the property given in Proposition 4 is illustrated in Table 1.
Moreover we have that the cdf of a S E F r ( 1 , λ , α ) approaches to the F r ( λ ) distribution when α . This appreciation is formalised in next proposition.
Proposition 5.
Let X α S E F r ( 1 , λ , α ) . If  α , then X α converges in distribution to a rv F F r ( λ ) .
Proof. 
Let X α S E F r ( 1 , λ , α ) . By using the stochastic representation given in (7) and the definition of the Slash Exponential proposed in [15], we can write X α = U 1 1 λ U 1 α λ where U 1 E x p ( 1 ) and U U ( 0 , 1 ) independent. Whereas U 1 α λ a . s 1 if α where a . s denotes for almost surely convergence, then we have,
U 1 α λ P 1 , α ,
where P signifies the convergence in probability. By applying Slutsky’s lemma, see for instance Lehmann [33], it follows that
U 1 1 λ U 1 α λ D U 1 1 λ = F F r ( λ ) ,
where D denotes convergence in distribution.    □
Note that Proposition 5 establishes that the for large α , the SEFr model can be approximated by a Fréchet distribution.

2.2. Moments

Proposition 6.
Let X S E F r ( σ , λ , α ) . Then, for r a positive integer, E X r exists, if and only if r < λ , and in this case,
E X r = σ r α λ α λ + r Γ 1 r λ .
Proof. 
Through the utilization of the stochastic representation provided in (7)
E X r = σ r E U 1 r / λ E U r α λ
where U 1 E x p ( 1 ) and U U ( 0 , 1 ) independent. Therefore, we have the proposed result.    □
The subsequent corollary is a direct result of Proposition 6, as indicated
Corollary 2.
Let X S E F r ( σ , λ , α ) . Then
1. 
E [ X ] = σ α λ α λ + 1 Γ 1 1 λ , provided that λ > 1 .
2. 
V a r [ X ] = σ 2 α λ 1 α λ + 2 Γ 1 2 λ α λ ( α λ + 1 ) 2 Γ 2 1 1 λ , provided that λ > 2 .
3. 
Let μ j = E [ X j ] . Then, the skewness, β 1 , and kurtosis, β 2 , coefficients can be obtained by using
β 1 = μ 3 3 μ 1 μ 2 + 2 μ 1 3 ( μ 2 μ 1 2 ) 3 2 , λ > 3 , β 2 = μ 4 4 μ 1 μ 3 + 6 μ 2 μ 1 2 3 μ 1 4 ( μ 2 μ 1 2 ) 2 , λ > 4 .
In Figure 3 and Figure 4, plots and contour plots for the skewness and kurtosis coefficients in the SEFr ( 1 , λ , α ) model are given. These plots suggest that for fixed α and an increasing value of λ the skewness and kurtosis coefficients decrease quickly. On the other hand, for  λ fixed and decreasing values of α , the skewness and kurtosis coefficients decrease quite slower.
  • Shannon Entropy.
Recall that the Shannon entropy of a continuous rv X, H ( X ) , is defined as
H ( X ) = E log f X ( X )
where log denotes neperian logarithm and f X is the pdf of X. A useful result which can be used to get the entropy of a rv is given in Lemma 3, c.f. Jones [34] (p. 137) or Awad [35].
Lemma 3.
If X is a continuous rv and Z = g ( X ) is one-to-one transformation, then the relationship between the Shannon entropies of X and Z is
H ( Z ) = H ( X ) + E log d Z d X .
Proposition 7.
Let X S E F r ( σ , λ , α ) . Then the Shannon entropy of X, H ( X ) , can be obtained as
H ( X ) = log σ λ α + α 1 λ E [ log Y ] E log γ ( 1 + α , Y ) ,
where Y S E ( α ) .
Proof. 
Taking into account the definition of the SEFr model given in (7) and (14)
H ( X ) = H ( Y ) + log σ λ 1 + λ λ E [ log Y ]
with Y S E ( α ) . Finally, by applying (13) to Y S E ( α ) , (15) follows.    □

2.3. Other Properties of Interest in Reliability

In this subsection additional properties of interest in reliability are given. Specifically, these are: the reverse hazard rate function, the stress-strength parameter and the pdf of the order statistics.
  • The reverse hazard rate function.
The reverse hazard rate function, r ( x ) , has received increasing interest in reliability field. This function is defined as the quotient of the pdf to its cdf, r ( x ) = f ( x ) / F ( x ) . Roughly speaking, r ( x ) gives the probability of a recent failure occurring given that a failure has already taken place, additional details can be seen in Block et al. [36]. Its expression for the S E F r ( σ , λ , α ) model is given next.
Corollary 3.
Let X S E F r ( σ , λ , α ) . Then the reverse hazard rate function of X is provided by
r ( x ) = λ α 2 Γ ( α ) x λ α 1 G σ x λ , 1 + α , 1 σ λ α α Γ ( α ) x σ λ α G σ x λ , 1 + α , 1 + exp σ x λ .
Proof. 
It follows by applying (8) and (10).    □
  • The stress-strength parameter.
Proposition 8.
Let X and Z independent Slash-Exponential-Fréchet rv’s, X S E F r ( σ x , λ x , α x ) , Z S E F r ( σ z , λ z , α z ) , and R, the stress-strength parameter. Then
R = λ x α x σ x 0 s σ x λ x α x 1 γ 1 + α x , σ x s λ x s σ z λ z α z γ 1 + α z , σ z s λ z + e x p σ z s λ z d s .
Proof. 
It follows from the fact that R = P [ Z < X ] , see [37], by using (8) and (10).    □
Corollary 4.
Let X and Z independent Slash-Exponential-Fréchet rv’s with the same scale parameter, X S E F r ( σ , λ x , α x ) , Z S E F r ( σ , λ z , α z ) . Then
R = λ x α x 0 t λ x α x 1 γ 1 + α x , t λ x t λ z α z γ 1 + α z , t λ z + exp ( t λ z ) d t .
Proof. 
It is immediate making the change of variable t = s / σ in the expression of R given in Proposition 8.    □
  • Order Statistics.
Let X 1 , , X n be independent identically distributed rv’s, X i S E F r ( σ , λ , α ) . Let us consider the order statistics from these n rv’s, denoted by X r : n , r = 1 , , n . Then the pdf of X r : n is
f r : n ( x ) = n ! ( r 1 ) ! ( n r ) ! f ( x ) { F ( x ) } r 1 { 1 F ( x ) } n r , r = 1 , , n ,
where f ( · ) and F ( · ) are the pdf and cdf of the S E F r ( σ , λ , α ) model given in (8) and (10).
For r = 1 the pdf of the minimum is obtained and for r = n the maximum. We highlight that the distribution of the maximum is not SEFr.

3. Inference

3.1. ML Estimators

Given X 1 , X 2 , , X n a random sample of size n from S E F r ( σ , λ , α ) , then from (8), the log-likelihood function is given by
l ( θ ) n log λ + n log α n λ α log ( σ ) + ( λ α 1 ) i = 1 n log ( x i ) + i = 1 n log γ α + 1 , x i σ λ ,
where θ = ( σ , λ , α ) , ∝ means proportional to.
The components of the score vector are procured by performing partial derivatives with respect to σ , λ , and  α , S ( θ ) = l σ , l λ , l α , these are
l σ = n λ α σ + λ σ i = 1 n x i σ λ ( α + 1 ) exp x i σ ) λ γ 1 + α , x i σ λ , l λ = n λ n α log ( σ ) + α i = 1 n log ( x i ) + i = 1 n log σ x i x i σ λ ( α + 1 ) exp x i σ ) λ γ 1 + α , x i σ λ , l α = n α n λ log ( σ ) + λ i = 1 n log ( x i ) + i = 1 n I α + 1 , x i σ λ γ 1 + α , x i σ λ ,
where I ( a , v ) = 0 v t a 1 log ( t ) e t d t , a > 0 , and  v > 0 . It can be seen in Milgram [38] that I ( a , v ) is related to the generalized integro-exponential function when v = .
The MLE of θ , represented as θ ^ , can be secured by finding the solution to the equation S ( θ ) = 0 through numerical methodologies, for instance, the Newton-Rapson algorithm. In another way, the MLEs can be directly secured by optimizing the log-likelihood function as provided in (17) and using the “BFGS” method of the “optim” subroutine in R software [39]. The “BFGS” method is a limited-memory quasi-Newton method for approximating the Hessian matrix of the target distribution. It is worth mentioning that the parameter vector θ = ( σ , λ , α ) can be easily obtained, thanks to the properties of the pdf f. The smooth and continuous nature of the function f, along with the existence and finiteness of its first and second derivatives, ensure that the equation S ( θ ) = 0 has roots. These roots correspond to the MLEs of the vector θ . By employing relevant calculus techniques, it is possible to verify that the solutions correspond to a maximum.

3.2. Observed Fisher Information Matrix

We can estimate the asymptotic variance of the MLEs, denoted as θ ^ = ( σ ^ , λ ^ , α ^ ) , using the Fisher information matrix. The Fisher information matrix, denoted as I ( θ ) , is calculated as the negative expectation of the second derivative of the log-likelihood function, l ( θ ) , with respect to θ . Under regularity conditions, the MLEs are asymptotically normal. In other words, as the sample size n approaches infinity, the distribution of I ( θ ) 1 / 2 θ ^ θ converges to a standard trivariate normal distribution, denoted as N 3 ( 0 3 , I 3 ) . To obtain I ( θ ) , we calculate the elements of the matrix 2 l ( θ ) / θ θ . Specifically, the elements are given by I σ σ = 2 l ( θ ) / σ 2 , I σ λ = 2 l ( θ ) / σ λ , and so on.
Typically, it is difficult to obtain an exact form for the expected value of previous expressions. Consequently, an estimate of the covariance matrix of MLEs, denoted as I ( θ ) 1 , can be obtained by evaluating the previous elements at the MLE θ ^ , this is called the observed information matrix, and is denoted as I ( θ ^ ) .
I ( θ ^ ) = 2 l ( θ ) / θ θ | θ = θ ^ .
The estimation of the asymptotic variances for σ ^ , λ ^ , and  α ^ is done by calculating the diagonal elements of I ( θ ^ ) 1 . The standard errors are then obtained by taking the square root of these asymptotic variances. For more information on the theoretical results used in this subsection, please refer to [40,41].
It is of interest to recall that the normality of maximum likelihood estimators is not satisfied for the case of small sample size. In this case, confidence interval estimates of parameters can be obtained by applying generalized estimation methods such as those described in Wang et al. [42] and Luo et al. [43].

3.3. Simulation Study

A simulation study has been carried out to assess the performance of ML estimators of σ , λ and α in the SEFr model. 1000 samples were generated with sample sizes n = 50 , 100 and 200. The different possibilities for the parameters σ , λ and α can be seen in Table 2. Next Algorithm 1 to generate values of X S E F r ( σ , λ , α ) is proposed.
Algorithm 1: To simulate values from the X S E F r ( σ , λ , α ) .
Step 1: Generate X 1 E x p ( 1 ) and X 2 U n i f o r m ( 0 , 1 ) .
Step 2: Compute Y = X 1 X 2 1 / α .
Step 3: Compute X = σ Y 1 / λ .
In Table 2, the bias, Standard Error (SE), Root of Mean Squared Error (RMSE), and the empirical Coverage Probability (CP) to 95% for the asymptotic intervals based on ML estimators are given. Note that the bias, SE, and RMSE decrease when the sample size n increases. These facts suggest that the MLE are consistent. Moreover, the empirical CP’s also approach to the nominal 95% when n increases.

4. Applications

In this section, we provide applications to two actual datasets. In each setting the fit provided by the SEFr model is compared to the Fréchet (Fr), Slashed Quasi-Gamma (SQG), and Slash Fréchet (SFr). The criteria for comparison are: the Akaike Information Criterion (AIC), given by Akaike [44], and the Bayesian Information Criterion (BIC) proposed by Schwarz [45]. Next the pdf’s of SQG and SFr are given.
  • Slashed Quasi-Gamma, X S Q G ( β , θ , q ) , introduced in [46]. Its pdf is:
    f X ( x ; β , θ , q ) = q β q x ( q + 1 ) Γ 1 10 Γ q θ + 1 10 F x β θ , q θ + 1 10 , 1 , x > 0 ,
    where β > 0 , θ > 0 , and q > 0 .
  • Slash Fréchet, X S F r ( λ , q ) , introduced in [20]. Its pdf is:
    f X ( x ; λ , q ) = q x q + 1 Γ 1 q λ , x λ , x > 0 ,
    where λ > q > 0 and Γ ( a , t ) = t w a 1 e w d w is the upper incomplete gamma function.

4.1. Application 1 (Patients with Bladder Cancer)

The remission times (in months) of 128 patients with bladder cancer are considered as first application. This dataset was first studied by Lee and Wang [47]. Table 3 provides the descriptive summaries. These are: sample mean, sample standard deviation (S), sample skewness coefficient ( b 1 ) and sample kurtosis coefficient ( b 2 ). Also, the boxplot is given in Figure 5. We highlight the high values of the skewness and kurtosis coefficients, and the existence of possible outliers in this dataset.
The results for the different models under consideration are listed in Table 4, that is, the estimates of parameters in the Fr, SFr, SQG, and SEFr models, along with their standard errors, log-likelihood, AIC, and BIC. It can be seen that the SEFr model achieves the lower values of AIC and BIC, and therefore it provides the best fit to this dataset.
The histogram and the fitted pdf’s for the previously mentioned distributions are given in Figure 6. In Figure 7, the corresponding QQ-plots can be seen. Both graphical summaries support our conclusions about the good fit of SEFr model to this dataset.

4.2. Application 2 (Air Conditioning System Failures)

The second dataset was first studied in Proschan [48]. It consists of the time interval (in hours) between successive failures of the air-conditioning system of Boeing 720 number 7912. Table 5 provides the descriptive summary and Figure 8 the corresponding boxplot. Note the presence of outliers and that a model with a long right tail seems appropriate. These summaries suggest that the proposed distributions can be good candidates for modelling this dataset.
The results for the fitted models are given in Table 6. According to AIC and BIC, the SEFr distribution provides a better fit than the other ones, since its AIC and BIC are the smallest ones.
The histogram for this dataset is given in Figure 9 along with the fitted pdf’s for Fr, SFr, SQG, and SEFr models. The QQ-plots are given in Figure 10. it can also be verified that the SEFr model outperforms other rival models in terms of fitting better. Hence, based on all these summaries and plot interpretations, it can be inferred that the SEFr model outperforms the Fr, SFr, and SQG distributions in fitting this dataset.

5. Conclusions

The Slash-Exponential-Fréchet distribution has been introduced. This model is an extension of the Fréchet, from which it inherits its main features, and moreover, it is more flexible as for its skewness and kurtosis. The following points have been considered:
  • The stochastic representation of the new model in terms of the Slash-Exponential is given. In this way, an additional shape parameter is added to Fréchet model.
  • Closed expressions for the pdf and cdf are given, therefore also for the survival and hazard rate function.
  • It is shown that the new model is unimodal or decreasing. It is proven that if the new shape parameter tends to infinity then the SEFr approaches to Fréchet model.
  • Closed expressions are given for the moments, with particular interest on skewness and kurtosis coefficients.
  • We highlight that the new model presents less kurtosis than the basal Fréchet distribution. For the best of our understanding, it is the first time in literature that, as result of applying slash methodology, the new model exhibits a lighter right tail and less kurtosis compared to basal model.
  • Maximum likelihood method has been proposed to estimate the parameters in the model. Score equations and the observed Fisher information matrix are studied.
  • A simulation study has been carried out. There, bias, standard error, RMSE and empirical coverage probability for MLEs have been obtained for increasing sample size. The good asymptotic properties of MLEs can be seen.
  • Two real applications are included where the SEFr model is compared to Fr, Slashed Quasi-Gamma and Slash-Fréchet. By using AIC and BIC, it has been seen that the new model provides a better fit compared to others.

Author Contributions

Conceptualization, Y.M.G., I.B.-C. and J.S.C.; software, J.S.C.; formal analysis, Y.M.G., H.W.G. and I.B.-C.; investigation, Y.M.G., I.B.-C. and J.S.C.; writing—original draft preparation, J.S.C., I.B.-C. and H.W.G.; writing—review and editing, Y.M.G., I.B.-C. and J.S.C.; supervision, Y.M.G. and H.W.G.; funding acquisition, H.W.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research of H.W. Gómez was supported by Semillero UA-2024. The research of I. Barranco-Chamorro was supported by IOAP of University of Seville, Spain.

Data Availability Statement

The datasets are available in the references given in Section 4.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fréchet, M. Sur la loi de probabilité de l’écart maximum. Ann. Soc. Polon. Math. 1927, 6, 93. [Google Scholar]
  2. Fisher, R.A.; Tippett, L.H.C. Limiting forms of the frequency distribution of the largest and smallest member of a sample. Proc. Camb. Philos. Soc. 1928, 24, 180–190. [Google Scholar] [CrossRef]
  3. Gumbel, E.J. Statistics of Extremes; Columbia University Press: New York, NY, USA, 1958. [Google Scholar]
  4. Embrechts, P.; Klüppelberg, C.; Mikosch, T. Modelling Extremal Events for Insurance and Finance; Springer: Berlin, Germany, 1997. [Google Scholar]
  5. Resnick, S.I. Extreme Values, Regular Variation and Point Processes; Springer: Berlin, Germany, 1987. [Google Scholar]
  6. Haan, L.; Ferreira, A. Extreme Value Theory: An Introduction; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  7. Kotz, S.; Nadarajah, S. Extreme Value Distributions: Theory and Applications; Imperial College Press: London, UK, 2000. [Google Scholar]
  8. Gupta, R.C.; Gupta, P.L.; Gupta, R.D. Modeling failure time data by Lehman alternatives. Commun. Stat. Methods 1998, 27, 887–904. [Google Scholar] [CrossRef]
  9. Coles, S. An Introduction to Statistical Modeling of Extreme Values; Springer: London, UK, 2001. [Google Scholar]
  10. Ramos, P.L.; Louzada, F.; Ramos, E.; Dey, S. The Fréchet distribution: Estimation and application—An overview. J. Stat. Manag. Syst. 2020, 23, 549–578. [Google Scholar] [CrossRef]
  11. Calabria, R.; Pulcini, G. Confidence limits for reliability and tolerance limits in the inverse Weibull distribution. Reliab. Eng. Syst. Saf. 1989, 24, 77–85. [Google Scholar] [CrossRef]
  12. Maswadah, M. Conditional confidence interval estimation for the inverse Weibull distribution based on censored generalized order statistics. J. Stat. Comput. Simul. 2003, 73, 887–898. [Google Scholar] [CrossRef]
  13. Salman, M.A.W.M.K.; Amer, S.S.M. Order statistics from inverse Weibull distribution and characterizations. Metron 2003, 61, 389–401. [Google Scholar]
  14. Abbas, K.; Tang, Y. Analysis of Fréchet distribution using reference priors. Commun. Stat. Theory Methods 2015, 44, 2945–2956. [Google Scholar] [CrossRef]
  15. Punathumparambath, B. Slash Exponential Distribution: Theory and Applications. Bull. Math. Stat. Res. 2020, 8, 38–49. [Google Scholar]
  16. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. In National Bureau of Standards Applied Mathematics Series 55; United State Department of Commerce: Washington, DC, USA, 1968. [Google Scholar]
  17. Gupta, R.D.; Kundu, D. Generalized Exponential Distributions. Aust. New Zealand J. Stat. 1999, 41, 173–188. [Google Scholar] [CrossRef]
  18. Gupta, R.D.; Kundu, D. Generalized Exponential Distribution: Different Methods of Estimations. J. Stat. Comput. Simul. 2001, 69, 315–337. [Google Scholar] [CrossRef]
  19. Astorga, J.M.; Gómez, H.W.; Bolfarine, H. Slashed generalized exponential distribution. Commun. Stat. Theory Methods 2017, 46, 2091–2102. [Google Scholar] [CrossRef]
  20. Castillo, J.S.; Rojas, M.A.; Reyes, J. A More Flexible Extension of the Fréchet Distribution Based on the Incomplete Gamma Function and Applications. Symmetry 2023, 15, 1608. [Google Scholar] [CrossRef]
  21. Rogers, W.H.; Tukey, J.W. Understanding some long-tailed symmetrical distributions. Stat. Neerl. 1972, 26, 211–226. [Google Scholar] [CrossRef]
  22. Andrews, D.F.; Bickel, P.J.; Hampel, F.R.; Huber, P.J.; Rogers, W.H.; Tukey, J.W. Robust Estimates of Location: Survey and Advances; Princeton University Press: Princeton, NJ, USA, 1972. [Google Scholar]
  23. Gómez, H.W.; Quintana, F.A.; Torres, F.J. A new Family of Slash-Distributions with Elliptical Contours. Stat. Probab. Lett. 2007, 77, 717–725. [Google Scholar] [CrossRef]
  24. Arslan, O.; Genc, A. A generalization of the multivariate slash distribution. J. Stat. Plan. Inference 2009, 139, 1164–1170. [Google Scholar] [CrossRef]
  25. Reyes, J.; Barranco-Chamorro, I.; Gómez, H.W. Generalized modified slash distribution with applications. Commun. Stat. Theory Methods 2020, 49, 2025–2048. [Google Scholar] [CrossRef]
  26. Del Castillo, J.M. Slash distributions of the sum of independent logistic random variables. Stat. Probab. Lett. 2016, 110, 111–118. [Google Scholar] [CrossRef]
  27. Zörnig, P. On Generalized Slash Distributions: Representation by Hypergeometric Functions. Stats 2019, 2, 371–387. [Google Scholar] [CrossRef]
  28. Olmos, N.M.; Varela, H.; Bolfarine, H.; Gómez, H.W. An extension of the generalized half-normal distribution. Stat. Pap. 2014, 55, 967–981. [Google Scholar] [CrossRef]
  29. Barranco-Chamorro, I.; Iriarte, Y.A.; Gómez, Y.M.; Astorga, J.M.; Gómez, H.W. A generalized Rayleigh family of distributions based on the modified slash model. Symmetry 2021, 13, 1226. [Google Scholar] [CrossRef]
  30. Barrios, L.; Gómez, Y.M.; Venegas, O.; Barranco-Chamorro, I.; Gómez, H.W. The Slashed Power Half-Normal Distribution with Applications. Mathematics 2022, 10, 1528. [Google Scholar] [CrossRef]
  31. Gui, W. Statistical properties and applications of the Lindley slash distribution. J. Appl. Stat. Sci. 2012, 20, 283–298. [Google Scholar]
  32. Castillo, J.S.; Barranco-Chamorro, I.; Venegas, O.; Gómez, H.W. Slash-Weighted Lindley Distribution: Properties, Inference, and Applications. Mathematics 2023, 11, 3980. [Google Scholar] [CrossRef]
  33. Lehman, L.E. Elements of Large-Sample Theory; Springer: New York, NY, USA, 1999. [Google Scholar]
  34. Jones, D. Elementary Information Theory; Clarendon Press: Oxford, UK, 1979. [Google Scholar]
  35. Awad, A.M. The Shannon entropy of generalized gamma and related distributions. In Proceedings of the First Jordamian Mathematics Conference, Amman, Jordan, 2–4 September 1991; pp. 13–27. [Google Scholar]
  36. Block, H.; Savits, T.; Singh, V. The reversed hazard rate function. Probab. Eng. Inf. Sci. 1998, 12, 69–70. [Google Scholar] [CrossRef]
  37. Johnson, R.A. Stress-Strength Models for Reliability. In Handbook of Statistics; North-Holland Press: Amsterdam, The Netherlands, 1988; Volume 7. [Google Scholar]
  38. Milgram, M. The generalized integro-exponential function. Math. Comput. 1985, 44, 443–458. [Google Scholar] [CrossRef]
  39. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2023; Available online: https://www.R-project.org/ (accessed on 25 March 2024).
  40. Rohatgi, V.K.; Saleh, A.K. An Introduction to Probability and Statistics, 3rd ed.; John Wiley& Sons: New York, NY, USA, 2001. [Google Scholar]
  41. Lehman, L.E.; Casella, G. Theory of Point Estimation, 2nd ed.; Springer: New York, NY, USA, 1998. [Google Scholar]
  42. Wang, W.; Cui, Z.; Chen, R.; Wang, Y.; Zhao, X. Regression analysis of clustered panel count data with additive mean models. Stat. Pap. 2023. [Google Scholar] [CrossRef]
  43. Luo, C.; Shen, L.; Xu, A. Modelling and estimation of system reliability under dynamic operating environments and lifetime ordering constraints. Reliab. Eng. Syst. Saf. 2022, 218 Pt A, 108136. [Google Scholar] [CrossRef]
  44. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  45. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  46. Iriarte, Y.A.; Varela, H.; Gómez, H.J.; Gómez, H.W. A Gamma-Type Distribution with Applications. Symmetry 2020, 12, 870. [Google Scholar] [CrossRef]
  47. Lee, E.T.; Wang, J. Statistical Methods for Survival Data Analysis; John Wiley & Sons: New York, NY, USA, 2003; Volume 476. [Google Scholar]
  48. Proschan, F. Theoretical explanation of observed decreasing failure rate. Technometrics 1963, 5, 375–383. [Google Scholar] [CrossRef]
Figure 1. pdf of S E F r ( 1 , λ , α ) : (a) for α = 0.5 and λ { 2 , 4 , 8 } , (b) for λ = 5 and α { 0.5 , 1 , 8 } .
Figure 1. pdf of S E F r ( 1 , λ , α ) : (a) for α = 0.5 and λ { 2 , 4 , 8 } , (b) for λ = 5 and α { 0.5 , 1 , 8 } .
Axioms 13 00253 g001
Figure 2. The cdf, survival function and hazard rate function for the SEFr ( 1 , λ , 0.5 ) model can be visualized through their respective plots.
Figure 2. The cdf, survival function and hazard rate function for the SEFr ( 1 , λ , 0.5 ) model can be visualized through their respective plots.
Axioms 13 00253 g002
Figure 3. Plot and contour plot of the skewness coefficient in the SEFr ( 1 , λ , α ) model.
Figure 3. Plot and contour plot of the skewness coefficient in the SEFr ( 1 , λ , α ) model.
Axioms 13 00253 g003
Figure 4. Plot and contour plot of the kurtosis coefficient in the SEFr ( 1 , λ , α ) model.
Figure 4. Plot and contour plot of the kurtosis coefficient in the SEFr ( 1 , λ , α ) model.
Axioms 13 00253 g004
Figure 5. Boxplot for remission times of patients with bladder cancer.
Figure 5. Boxplot for remission times of patients with bladder cancer.
Axioms 13 00253 g005
Figure 6. Fitted pdf for the remission times of patients with bladder cancer dataset in the Fr and SFr, SQG, and SEFr distributions.
Figure 6. Fitted pdf for the remission times of patients with bladder cancer dataset in the Fr and SFr, SQG, and SEFr distributions.
Axioms 13 00253 g006
Figure 7. QQ plots for the remission times of patients with bladder cancer dataset: (a) SEFr Model; (b) SQG model; (c) SFr model; (d) Fr model.
Figure 7. QQ plots for the remission times of patients with bladder cancer dataset: (a) SEFr Model; (b) SQG model; (c) SFr model; (d) Fr model.
Axioms 13 00253 g007
Figure 8. Boxplot of the Air Conditioning System Failures dataset.
Figure 8. Boxplot of the Air Conditioning System Failures dataset.
Axioms 13 00253 g008
Figure 9. Fitted Density for Air Conditioning System Failures dataset in the Fr and SFr, SQG, and SEFr distributions.
Figure 9. Fitted Density for Air Conditioning System Failures dataset in the Fr and SFr, SQG, and SEFr distributions.
Axioms 13 00253 g009
Figure 10. QQ plots for Air Conditioning System Failures dataset: (a) SEFr Model; (b) SQG model; (c) SFr model; (d) Fr model.
Figure 10. QQ plots for Air Conditioning System Failures dataset: (a) SEFr Model; (b) SQG model; (c) SFr model; (d) Fr model.
Axioms 13 00253 g010
Table 1. Some probabilities on right tails of SEFr and Fr distributions.
Table 1. Some probabilities on right tails of SEFr and Fr distributions.
α P [ X > 1 ] P [ X > 1.5 ] P [ X > 2 ] P [ X > 2.5 ]
S E F r ( 1 , 5 , 0.5 ) 0.2530.0420.0100.003
S E F r ( 1 , 5 , 1 ) 0.3680.0630.0160.005
S E F r ( 1 , 5 , 8 ) 0.5870.1100.0270.009
F r ( 5 ) 0.6320.1230.0310.010
Table 2. Bias, SE, RMSE, and empirical CP for the simulation results in SEFr model.
Table 2. Bias, SE, RMSE, and empirical CP for the simulation results in SEFr model.
True Value n = 50n = 100n = 200
σ λ α EstimatorBiasSERMSECPBiasSERMSECPBiasSERMSECP
1.520.5 σ ^ 0.0290.3050.34790.50.0090.2080.23090.80.0020.1490.15393.8
λ ^ 0.3160.6822.06794.00.0990.3560.40795.10.0370.2410.26293.3
α ^ 0.0730.3740.72690.50.0150.1330.14493.70.0070.0910.09693.5
1.2 σ ^ −0.0370.2530.29190.6−0.0250.1770.19693.8−0.0060.1220.12494.9
λ ^ 0.1240.4040.59892.80.0270.2590.27894.10.0240.1810.18595.1
α ^ 1.5237.6294.55891.20.4271.4601.79594.80.0770.3280.43693.8
40.5 σ ^ −0.0100.1490.17689.1−0.0040.1060.10993.3−0.0030.0740.07793.9
λ ^ 0.5121.2223.10092.90.1870.7220.85494.80.0840.4830.52194.7
α ^ 0.1040.4430.71589.90.0190.1370.15593.50.0080.0920.09494.1
1.2 σ ^ −0.0190.1320.14092.4−0.0140.0890.09795.3−0.0040.0610.06294.8
λ ^ 0.2560.8101.02393.40.0960.5270.58393.10.0490.3620.36895.2
α ^ 1.2276.4904.11491.10.3971.4522.01392.40.0660.3170.41094.5
230.7 σ ^ −0.0180.2420.26290.9−0.0050.1700.18392.9−0.0020.1200.12594.3
λ ^ 0.2590.7330.94694.60.1080.4720.51895.30.0650.3230.34395.6
α ^ 0.1800.8421.15392.00.0520.2520.43494.00.0130.1330.14493.4
1 σ ^ −0.0280.2340.27291.2−0.0090.1610.17295.2−0.0020.1110.11195.7
λ ^ 0.1800.6310.82092.20.0680.4160.43995.30.0400.2860.29794.8
α ^ 0.8894.2613.23391.30.1640.5950.91293.50.0360.2220.26694.6
50.7 σ ^ −0.0160.1470.16691.70.0030.1030.10994.40.0010.0720.07494.5
λ ^ 0.4421.2381.67993.50.2180.7970.88595.70.1050.5380.56196.1
α ^ 0.1970.8981.11091.80.0240.2040.23493.10.0090.1320.14593.8
1 σ ^ −0.0250.1470.16592.9−0.0080.0960.10193.6−0.0050.0670.06794.7
λ ^ 0.2671.0451.30591.70.1170.6900.73794.60.0810.4800.51394.0
α ^ 0.9746.5293.37092.40.1520.5610.85893.90.0460.2220.25495.2
31.50.3 σ ^ 0.1700.9571.18488.10.0540.6840.73392.10.0200.4750.49594.5
λ ^ 0.5350.9663.14793.90.1120.3500.42394.40.0550.2270.25195.0
α ^ 0.0130.1230.15288.90.0050.0830.09892.10.0020.0570.05994.6
0.9 σ ^ −0.0330.6890.76690.6−0.0140.4870.52692.8−0.0060.3400.33895.0
λ ^ 0.1040.3270.39094.40.0390.2150.24993.90.0230.1480.15594.7
α ^ 0.6222.7552.59891.70.1420.5440.80793.40.0340.1890.20995.2
3.50.3 σ ^ 0.0220.3940.46189.90.0070.2870.30692.50.0010.2040.20295.3
λ ^ 0.9231.6494.37893.20.2890.8321.18295.20.0930.5240.54795.5
α ^ 0.0150.1230.16488.30.0040.0810.08992.10.0030.0570.05794.5
0.9 σ ^ −0.0210.3010.35289.7−0.0050.2090.23293.2−0.0030.1450.15093.3
λ ^ 0.2800.7941.34493.70.1230.5090.54995.50.0640.3460.37095.1
α ^ 0.6813.2112.88590.60.1000.4380.72894.90.0220.1840.20293.8
Table 3. Descriptive Summary for Patients with Bladder Cancer dataset.
Table 3. Descriptive Summary for Patients with Bladder Cancer dataset.
n x ¯ S b 1 b 2
1289.36610.5083.28718.483
Table 4. Estimates, SE in parentheses, log-likelihood, AIC, and BIC values for the remission times of patients with bladder cancer dataset.
Table 4. Estimates, SE in parentheses, log-likelihood, AIC, and BIC values for the remission times of patients with bladder cancer dataset.
ParametersFr (SE)SFr (SE)SQG (SE)SEFr (SE)
σ ^ ---9.9436 (1.3592)
λ ^ 0.6726 (0.0479)0.9242 (0.0688)-1.8586 (0.2660)
α ^ ---0.6329 (0.1558)
β ^ --7.7993 (0.9893)-
θ ^ --10.8627 (1.3211)-
q ^ -0.9623 (0.1302)1.5211 (0.2282)-
log-likelihood−481.0559−448.1104−411.7342−410.0634
AIC964.1118900.2208829.4683826.1268
BIC966.9638905.9249838.0244834.6829
Table 5. Descriptive statistics for the Air Conditioning System Failures dataset.
Table 5. Descriptive statistics for the Air Conditioning System Failures dataset.
n x ¯ s b 1 b 2
3059.633371.89961.69144.9595
Table 6. Estimates of parameters, SE in parentheses, log-likelihood, AIC and BIC values for Air Conditioning System Failures.
Table 6. Estimates of parameters, SE in parentheses, log-likelihood, AIC and BIC values for Air Conditioning System Failures.
ParametersFr (SE)SFr (SE)SQG (SE)SEFr (SE)
σ ^ ---38.5732 (13.9807)
λ ^ 0.3924 (0.0601)0.9508 (0.3089)-1.0968 (0.2273)
α ^ ---1.1344 (0.5565)
β ^ --14.7397 (3.1532)-
θ ^ --14.3648 (4.2492)-
q ^ -0.4067 (0.0960)0.7165 (0.1559)-
log-likelihood−177.5930−163.9272−153.0741−152.3953
AIC357.1859331.8543312.1481310.7905
BIC358.5871334.6567316.3517314.9941
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gómez, Y.M.; Barranco-Chamorro, I.; Castillo, J.S.; Gómez, H.W. An Extension of the Fréchet Distribution and Applications. Axioms 2024, 13, 253. https://doi.org/10.3390/axioms13040253

AMA Style

Gómez YM, Barranco-Chamorro I, Castillo JS, Gómez HW. An Extension of the Fréchet Distribution and Applications. Axioms. 2024; 13(4):253. https://doi.org/10.3390/axioms13040253

Chicago/Turabian Style

Gómez, Yolanda M., Inmaculada Barranco-Chamorro, Jaime S. Castillo, and Héctor W. Gómez. 2024. "An Extension of the Fréchet Distribution and Applications" Axioms 13, no. 4: 253. https://doi.org/10.3390/axioms13040253

APA Style

Gómez, Y. M., Barranco-Chamorro, I., Castillo, J. S., & Gómez, H. W. (2024). An Extension of the Fréchet Distribution and Applications. Axioms, 13(4), 253. https://doi.org/10.3390/axioms13040253

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop