Next Article in Journal
The Role of Risk Factors in System Performance: A Comprehensive Study with Adaptive Progressive Type-II Censoring
Previous Article in Journal
Matrix Pencil Optimal Iterative Algorithms and Restarted Versions for Linear Matrix Equation and Pseudoinverse
Previous Article in Special Issue
Exploring Complex Survival Data through Frailty Modeling and Regularization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Generalization of the Truncated Gumbel Distribution with Quantile Regression and Applications

1
Departamento de Ciencias Matemáticas y Físicas, Facultad de Ingeniería, Universidad Católica de Temuco, Temuco 4780000, Chile
2
Departamento de Estadística y Ciencia de Datos, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta 1240000, Chile
3
Departamento de Matemáticas, Facultad de Ciencias, Universidad Católica del Norte, Antofagasta 1240000, Chile
4
Facultad de Ciencias, Universidad Arturo Prat, Avenida Arturo Prat 2120, Iquique 1110939, Chile
5
Departamento de Estadística, Facultad de Ciencias, Universidad del Bío-Bío, Concepción 4081112, Chile
6
Department of Statistics, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora 36036-900, MG, Brazil
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(11), 1762; https://doi.org/10.3390/math12111762
Submission received: 8 April 2024 / Revised: 20 May 2024 / Accepted: 29 May 2024 / Published: 5 June 2024
(This article belongs to the Special Issue Mathematical and Computational Statistics and Applications)

Abstract

:
In this article, we introduce a new model with positive support. This model is an extension of the truncated Gumbel distribution, where a shape parameter is incorporated that provides greater flexibility to the new model. The model is parameterized in terms of the p-th quantile of the distribution to perform quantile regression in this model. An extensive simulation study demonstrates the good performance of the maximum likelihood estimators in finite samples. Finally, two applications to real datasets related to the level of beta-carotene and body mass index are presented.

1. Introduction

The Gumbel distribution, also known as the type-I generalized extreme value distribution, is commonly used to model data with extreme observations. This distribution and its extensions have a wide range of applications in several disciplines such as hydrology, economics, finance, climatology and seismology. The probability density function (pdf), the cumulative distribution function (cdf) and the quantile function of a random variable X that follows the Gumbel distribution are given by
g ( x ; μ , σ ) = 1 σ exp x μ σ exp x μ σ , x R , G ( x ; μ , σ ) = exp exp x μ σ , x R , Q ( p ) = μ σ log log ( p ) , p ( 0 , 1 ) ,
where μ R and σ > 0 . Applications of this model in different scenarios can be found in Bhaskaran et al. [1], Gurung et al. [2], Purohit et al. [3], Li et al. [4] and Kang et al. [5].
Several extensions of the Gumbel distribution have been recently proposed in the literature. Hossam et al. [6] presented a statistical model that combines the new alpha power transformation method and Gumbel distribution. Watthanawisut and Bodhisuwan [7] proposed a new extension of the so called Topp-Leone Gumbel distribution that is used to model minimum flow data. Fayomi et al. [8] presented the exponentiated Gumbel-G family of distributions and explored a special case called EGuNH. Nagode et al. [9] introduced a three-parameter Gumbel distribution, which was applied to rope failure data. Oseni and Okasha [10] derived the Gumbel-geometric distribution, which was applied to precipitation and maximum annual wind speed data. Note that all these extensions do not consider a regression framework, and their main objectives rely on the fit of univariate data.
It is evident that regression models have become relevant tools in the era of Data Science. Among them, the so-called quantile regression models (introduced by Koenker and Bassett [11]) are an alternative to the usual regression techniques where the mean response conditional to values of covariates (or explanatory variables) is estimated. The quantile regression models allow us to measure the effects of covariates at different quantiles of the response variable distribution. Thus, they provide an analysis across the entire conditional distribution, as can be seen in Cade et al. [12], Koenker [13] and Wei et al. [14]. The mean, as the only summary measure, is generally quite poor for assessing risk, as it is greatly affected by the presence of outlier observations. Outlier data can be quite strange, but at the same time, these can be enough to cause serious problems when analyzing the information obtained; see, for example, Gómez-Déniz et al. [15], who analyzed extreme values in insurance companies. To our knowledge, there are no studies on quantile regression models based on the Gumbel distribution. Thus, the objectives of this work were to introduce a new generalization of the truncated Gumbel distribution and then establish a quantile regression model based on this novel generalization. To do this, a reparametrization was obtained of the new truncated Gumbel generalization by incorporating a parameter that represents the quantile. We should note that the proposed generalization was achieved by considering the work of Neamah and Qasim [16] and the transformation provided by Cooray and Ananda [17]. The latter authors developed an extension of the half-normal (HN) distribution through the relation Y = β X 1 / α , where X HN ( 1 ) .
The rest of the paper is organized as follows. In Section 2, we introduce our proposal, the generalized truncated Gumbel (GTG) distribution, and several important properties of this new model are presented. In Section 3, inference is performed, including some initial points to obtain maximum likelihood (ML) estimators and present the observed Fisher information matrix for the proposed model. In Section 4, the reparametrized model in terms of a quantile is presented. In Section 5, the simulation study carried out to analyze the performances of the ML estimators in finite samples for the proposed model without and with covariates is discussed. In Section 6, two real-data applications are presented to illustrate the proposed models, without and with covariates. Finally, in Section 7, some concluding comments are presented.

2. Generalized Truncated Gumbel Distribution

Neamah and Qasim [16] derived a new model with positive support for the Gumbel distribution by truncating its pdf from the left. We will refer to the resulting model of these authors as the truncated Gumbel (TG) distribution, which is defined in the interval ( 0 , ) . In considering the reparametrization λ = μ / σ , the pdf of the TG distribution can be written as follows:
f ( y ; β , λ , α ) = 1 β ( 1 G ( λ ) ) g y β λ , y > 0 ,
where β > 0 is a scale parameter, λ R is a shape parameter, and g ( u ) = exp ( u exp ( u ) ) and G ( u ) = exp ( exp ( u ) ) are the pdf and cdf for the standard Gumbel distribution, respectively.
In this work, we considered the transformation developed by Cooray and Ananda [17] to extend the TG distribution. That is, we considered the transformation Z = β Y 1 / α , where Y TG ( 1 , λ ) . We will refer to this extension as the generalized truncated Gumbel (GTG) distribution. Important functions, such as the pdf, cdf, hazard and quantile functions of the GTG distribution are provided below.

2.1. Pdf, cdf and Hazard Function

Proposition 1. 
Let Z G T G ( β , λ , α ) . Then, the pdf of Z is given by
f ( z ; β , λ , α ) = α z α 1 β α ( 1 G ( λ ) ) g z β α λ , z > 0 ,
where β > 0 , λ R and α > 0 .
Proof. 
Considering the stochastic representation of Z, we have that z = h ( y ) = β y 1 / α . Then,
f Z ( z ) = f Y h 1 ( z ) d h 1 ( z ) d z = 1 1 G ( λ ) g ( h 1 ( z ) λ ) α z α 1 β α .
Therefore, the result is obtained by replacing h 1 ( z ) = z β α in g ( · ) .    □
Remark 1. 
We previously mentioned that Z G T G ( β , λ , α ) if Z = β Y 1 / α , where Y T G ( 1 , λ ) . Thus, when α = 1 , we obtain the TG distribution; that is, G T G ( β , λ , 1 ) T G ( β , λ ) .
Proposition 2. 
Let Z G T G ( β , λ , α ) . Then, the cdf and hazard function of Z are given by
F ( z ; β , λ , α ) = G z β α λ G ( λ ) 1 G ( λ )
and
h ( z ; β , λ , α ) = α z α 1 g z β α λ β α 1 G z β α λ ,
respectively, for all z > 0 .
Proof. 
Both functions are obtained immediately from their definitions.   □
Figure 1 shows the pdf, cdf and hazard function for the G T G ( 1 , λ , α ) model, considering some combinations for λ and α . We observe that the GTG model can have decreasing or unimodal shapes for the pdf, whereas for the hazard function, we can have decreasing or increasing shapes. Also, we observe that for some combinations of λ and α , the cdf rapidly increases, although all of them tend to 1 when z increases.

2.2. Mode

The shape of the pdf of Z G T G ( β , λ , α ) can be examined based on its inflection points. By computing the first derivative of log ( f ( z ) ) with respect to z, where f ( z ) is the pdf for the GTG model, we obtain that
log ( f ( z ) ) z = α 1 z α z α 1 β α 1 exp ( v ) ,
where v = z β α λ . By equating the previous expression to 0, we obtain that
β α ( α 1 ) α = z α 1 exp ( v ) ,
from which the mode of Z can be numerically obtained. The nature of the points are determined by 2 log ( f ( z ) ) / z 2 = u ( z ) , where u ( z ) is given by
u ( z ) = α 1 z 2 ( 1 exp ( v ) ) α ( α 1 ) z α 2 β α exp ( v ) α 2 z 2 ( α 1 ) β 2 α .
Depending on whether u ( z 0 ) < 0 or u ( z 0 ) > 0 , where z = z 0 is a solution of Equation (5), the inflection points can be local maxima or minima. Figure 2 shows the shape of u ( z ) for β = 1 and selected values of λ and α . From here, we observe that the pdf of the GTG distribution is zero when z , both for when λ takes a positive or negative value.

2.3. Quantiles

Proposition 3. 
If Z G T G ( β , λ , α ) , then the quantile function of Z is given by
Q ( p ) = β λ + G 1 p + ( 1 p ) G ( λ ) + G ( λ ) 1 α , 0 < p < 1 .
Proof. 
It follows from a direct computation, by applying the definition of the quantile function.    □
Corollary 1. 
The quartiles of the GTG distribution are as follows:
1. 
(First quartile) Q ( 0.25 ) = β λ + G 1 0.25 + 0.75 G ( λ ) 1 α .
2. 
(Median) Q ( 0.5 ) = β λ + G 1 0.5 ( 1 + G ( λ ) ) 1 α .
3. 
(Third quartile) Q ( 0.75 ) = β λ + G 1 0.75 + 0.25 G ( λ ) 1 α .
Proof. 
It is immediate from Proposition 3.    □

2.4. Moments

Proposition 4. 
Let Z G T G ( β , λ , α ) and n be a positive integer. Then, the n-th moment of Z is given by
μ n = E [ Z n ] = β n 1 G ( λ ) b n ( λ , α ) ,
where b n ( λ , α ) = k = 0 ( 1 ) k n / α k λ n / α k γ k , γ k = 0 e λ ( ln ( t ) ) k e t d t , and n / α k = 1 k ! Π i = 0 k 1 ( n / α i ) is the generalized binomial coefficient. When n / α N , the sum in b n ( α , λ ) stops at n / α .
Proof. 
Given the stochastic representation of Z, it is immediate that E [ Z n ] = β n E [ Y n / α ] , where Y T G ( 1 , λ ) . Then, the  ( n / α ) moment of Y can be computed by following the properties presented in Neamah and Qasim [16].    □
Corollary 2. 
If Z G T G ( β , λ , α ) , then the first four moments and the variance of Z are obtained as follows:
1. 
μ 1 = E ( Z ) = β 1 G ( λ ) b 1 ( λ , α ) ;
2. 
μ 2 = E ( Z 2 ) = β 2 1 G ( λ ) b 2 ( λ , α ) ;
3. 
μ 3 = E ( Z 3 ) = β 3 1 G ( λ ) b 3 ( λ , α ) ;
4. 
μ 4 = E ( Z 4 ) = β 4 1 G ( λ ) b 4 ( λ , α ) ;
5. 
V a r ( Z ) = β 2 ( 1 G ( λ ) ) 2 ( 1 G ( λ ) ) b 2 ( λ , α ) b 1 2 ( λ , α ) .
Proof. 
It is immediate from Proposition 4.    □
Corollary 3. 
Let Z G T G ( β , λ , α ) . Then, the skewness coefficient ( β 1 ) and the kurtosis coefficient ( β 2 ) are given by
β 1 = G ¯ 2 b 3 3 G ¯ b 1 b 2 + 2 b 1 3 ( G ¯ b 2 b 1 2 ) 3 / 2 a n d β 2 = G ¯ 3 b 4 4 G ¯ 2 b 1 b 3 + 6 G ¯ b 1 2 b 2 3 b 1 4 ( G ¯ b 2 b 1 2 ) 2 ,
where G ¯ k = ( 1 G ( λ ) ) k , and b i = b i ( λ , α ) , for  i = 1 , 2 , 3 , 4 .
Proof. 
The expressions above are obtained using the definitions of the skewness and kurtosis coefficients; that is,
β 1 = μ 3 3 μ 1 μ 2 + 2 μ 1 3 ( μ 2 μ 1 2 ) 3 / 2 a n d β 2 = μ 4 4 μ 1 μ 3 + 6 μ 1 2 μ 2 3 μ 1 4 ( μ 2 μ 1 2 ) 2 .
where μ n = E [ X n ] , for  n = 1 , 2 , 3 , 4 , are given in Corollary 2.    □
Remark 2. 
Proposition 4 shows that the moments of the G T G distribution basically depend on the moments of the T G ( 1 , λ ) model. Plots for the expected value, variance, skewness and kurtosis coefficients of the G T G ( 1 , λ , α ) model are given in Figure 3 for different values of the λ and α parameters. The bottom plots in Figure 3 reflect the effect of the α parameter: a lower value of α produces higher values of the skewness and kurtosis coefficients. This fact can also be appreciated in Table 1 and Table 2.

2.5. Bonferroni Curves

In different disciplines, such as socio-economics and public health sciences, there is a necessity to compare and analyze the inequality of non-negative distributions. Generally, Bonferroni curves are used as graphical methods to achieve the required comparison/analysis (see Bonferroni [18], and Arcagni and Porro [19] for a further discussion about these curves). The following result shows the expressions of these curves for the GTG model.
Proposition 5. 
If Z G T G ( β , λ , α ) , then the Bonferroni curves, say B ( p ) , are given by
B ( p ) = β μ 1 p [ 1 G ( λ ) ] k = 0 ( 1 ) k 1 / α k λ 1 / α k γ k , 0 < p < 1 ,
where μ 1 = E ( Z ) , γ k = v e λ ( ln ( t ) ) k e t d t , and v = exp ( q / β ) α + λ .
Proof. 
The expression above is obtained using the definition of the Bonferroni curves; that is,
B ( p ) = 1 μ p 0 p F 1 ( t ) d t = 1 μ p 0 q z f ( z ) d z , 0 < p < 1 ,
where μ is the expected value of the corresponding non-negative random variable, and q = F 1 ( p ) .    □
Figure 4 shows the Bonferroni curves for the G T G ( 1 , λ , α ) model, considering different values for λ and α .

3. Inference for the GTG Distribution

In this section, we discuss the maximum likelihood (ML) approach for parameter estimation in the GTG model.

3.1. Maximum Likelihood Estimators

Let z 1 , z 2 , , z n be a random sample of size n from the GTG ( β , λ , α ) model. Then, the log-likelihood function for θ = ( β , λ , α ) is given by
( θ ) = n log ( α ) + ( α 1 ) i = 1 n log ( z i ) n α log ( β ) n log ( 1 G ( λ ) ) i = 1 n z i β α + n λ i = 1 n exp v i ,
where v i = z i β α λ . Therefore, the score assumes the form S ( θ ) = ( S β ( θ ) , S λ ( θ ) , S α ( θ ) ) , where
S β ( θ ) = n α β + i = 1 n α z i α β α + 1 1 exp v i ,
S λ ( θ ) = n g ( λ ) 1 G ( λ ) + n i = 1 n exp v i ,
and
S α ( θ ) = n α + i = 1 n log ( z i ) n log ( β ) + i = 1 n z i β α log z i β exp ( v i ) 1 .
The ML estimators are then obtained by numerically solving the equation S ( θ ) = 0 3 , where 0 p denotes a vector of zeros with length p. Solutions for Equations (9)–(11) can be obtained using numerical procedures in R [20], such as the Newton–Raphson method. To initialize the numerical algorithm that solves S ( θ ) = 0 3 , in the next subsection, we propose an initial point for the vector θ .

3.2. Initial Points

In this subsection, we propose estimators based on the quantiles for the GTG distribution, and these estimators are an alternative to the moment estimators, which meets the objective of using them as initial values to calculate the maximum likelihood estimators of the GTG distribution.
Let q 1 , q 2 and q 3 be the sample quartiles that are based on z 1 , z 2 , , z n . Initial values for θ can be obtained by equating the sample quartiles with the theoretical quartiles. The resulting equations are given by
q 1 = β G 1 1 4 + 3 4 G ( λ ) + λ 1 α , q 2 = β G 1 1 2 + 1 2 G ( λ ) + λ 1 α ,
and
q 3 = β G 1 3 4 + 1 4 G ( λ ) + λ 1 α .
The solutions for β and α , say β ˜ and α ˜ , can be expressed in terms of λ ˜ (the solution for λ ) as follows:
β ˜ = β ˜ ( λ ˜ ) = q 1 G 1 1 4 + 3 4 G ( λ ˜ ) + λ ˜ 1 / α ˜ ( λ ˜ ) a n d α ˜ = α ˜ ( λ ˜ ) = log G 1 3 4 + 1 4 G ( λ ˜ ) + λ ˜ G 1 1 4 + 3 4 G ( λ ˜ ) + λ ˜ log q 3 q 1 ,
whereas λ ˜ is obtained from the non-linear equation
β ˜ ( λ ˜ ) G 1 1 2 + 1 2 G ( λ ˜ ) + λ ˜ 1 / α ˜ ( λ ˜ ) = q 2 .
Therefore, the initial point based on this method is given by θ quart = β ˜ , λ ˜ , α ˜ .

3.3. Observed Fisher Information Matrix

The asymptotic variance of the ML estimators, say θ ^ = ( β ^ , λ ^ , α ^ ) , can be estimated from the observed Fisher information matrix defined as I ( θ ) = E 2 ( θ ) / θ θ , with  ( θ ) given in Equation (8). Under regularity conditions,
I ( θ ) 1 / 2 θ ^ θ D N 3 ( 0 3 , I 3 ) , a s   n + ,
where D stands for convergence in distribution, and N 3 ( 0 3 , I 3 ) denotes the standard trivariate normal distribution (see Wang et al. [21]). Moreover, I ( θ ) can be estimated from the matrix 2 ( θ ) / θ θ , whose elements are given by I β β = 2 ( θ ) / β 2 , I β λ = 2 ( θ ) / β λ , and so on. Explicitly, we have that
I β β = n α β 2 + ( α + 1 ) α β 2 i = 1 n z i β α ( 1 exp ( v i ) ) + α 2 β 2 ( α + 1 ) i = 1 n z i 2 α exp ( v i ) , I β λ = α β α + 1 i = 1 n z i α exp ( v i ) , I β α = n β + i = 1 n z i α β α + 1 1 exp ( v i ) 1 + α log z i β α β α + 1 i = 1 n z i α z i β α log z i β exp ( v i ) I λ λ = n g 2 ( λ ) ( 1 G ( λ ) ) 2 n g ( λ ) ( 1 exp ( λ ) ) 1 G ( λ ) + i = 1 n exp ( v i ) , I λ α = i = 1 n z i β α log z i β exp ( v i ) , I α α = n α 2 i = 1 n log 2 z i β z i β α exp ( v i ) 1 + i = 1 n log 2 z i β z i β 2 α exp ( v i ) ,
where v i = z i β λ .

4. GTG Quantile Regression Model

For the GTG model, the mean has a complicated form, and then, it is not recommendable to consider a mean-parameterized version of the model. On the other hand, and thinking in a context of heterogeneous observations, quantile regression is a more appropriate tool for analyzing data in presence of covariates because they allow for a complete description of the distribution of the response variable (not just a particular measure as is the case when regression on the mean is conducted).
Specifically, for the GTG model and considering that τ represents the pth quantile of the distribution, we obtain the equation τ = Q ( p ; β , α , λ ) , τ ( 0 , ) . By solving such an equation, we obtain
β = τ ( h ( λ , p ) ) 1 / α ,
where h ( λ , p ) = λ + G 1 [ p ( 1 G ( λ ) ) + G ( λ ) ] .
Thus, we can reparameterize the pdf and cdf of the GTG model as
f ( z ; τ , λ , α ) = α z α 1 h ( λ , p ) τ α ( 1 G ( λ ) ) g z τ α h ( λ , p ) λ
and
F ( z ; τ , λ , α ) = G z τ α h ( λ , p ) λ G ( λ ) 1 G ( λ ) ,
respectively, where z > 0 , α > 0 , τ ( 0 , ) , and 0 < p < 1 is fixed. We refer to this model as the reparameterized GTG (RGTG) model.
The consideration of z i = ( z i 1 , z i 2 , , z i q ) as a set of q known covariates related to the p-th quantil of the i-th individual can be introduced in the model as follows:
ψ ( τ i ( p ) ) = z i β ( p ) ,
where β ( p ) = ( β 1 ( p ) , β 2 ( p ) , , β q ( p ) ) is a q-dimensional vector of unknown regression parameters ( q < n ), and ψ ( · ) is a link function, which is continuous, invertible and at least twice differentiable. A natural choice in this context is the logarithm link, i.e.,  ψ ( u ) = log ( u ) .
With this framework, the corresponding log-likelihood function for the RGTG quantile regression model is given by
( θ ( p ) ) = k n ( p ) + ( α ( p ) 1 ) i = 1 n log ( z i ) α ( p ) i = 1 n log ( τ i ( p ) ) h ( λ ( p ) , p ) i = 1 n z i τ i ( p ) α ( p ) + i = 1 n exp z i τ i ( p ) α ( p ) h ( λ ( p ) , p ) + λ ( p ) ,
where k n ( p ) = n log ( α ( p ) ) + h ( λ ( p ) , p ) log ( 1 G ( λ ( p ) ) ) + λ ( p ) . The estimation of the regression parameters is obtained by directly maximizing this function.

5. Simulation Study

In this section, we present two simulation studies related to assessing the performances of the ML estimators for the GTG model and the RGTG quantile regression model.

5.1. Without Covariates

In this study, we carried out a simulation study to evaluate the performances of the ML estimators given in Section 3.1. We generated random values from the GTG ( β , λ , α ) distribution with Algorithm 1.
Algorithm 1 Simulating values from the G T G ( β , λ , α ) distribution
1:
Simulate U U n i f o r m ( 0 , 1 )
2:
Compute Z = β λ + G 1 [ U ( 1 G ( λ ) ) + G ( λ ) ] 1 α
We used the following sequence to perform a simulation study to evaluate the behavior in finite samples of the MLEs of the GTG model. For β , we fixed three values: 1, 2 and 3; for λ , we fixed two values: 2 and 3; for α , we fixed two values: 1 and 2; and for the sample size n we fixed four values: 150, 300, 600 and 1000. For each combination of β , λ , α and n, we simulated 1000 replicates of that size and calculated ML estimators and their standard errors. Table 3 summarizes the mean of the estimated biases (bias), the mean of the estimated standard errors (SE), and the squared root of the estimated mean squared errors (RMSE), and each estimated coverage probability (CP) was obtained by taking into account the asymptotic distribution of the ML estimator using a 95% confidence level. Note that as the sample size increases, the bias, SE and RMSE decrease, which suggests that the ML estimators of the GTG model have an acceptable behavior even in finite samples. Moreover, the SE and RMSE terms tend to become closer as the sample size increases, suggesting that the variance of the estimators is well estimated. Finally, the CP terms come closer to the nominal value as n increases, which suggests that the asymptotic approach to the normal of the ML estimators of the GTG model is reasonable, even in finite samples.

5.2. With Covariates

Under this condition, we generated 3000 random samples with a sample size ranging from n { 100 , 150 , , 1000 } from the RGTG quantile regression model using three different values for p : { 0.1 , 0.5 , 0.9 } . The structure of the scenarios is defined as follows:
log ( ρ i ( p ) ) = 2.071 + 0.013 z 1 i + 0.034 z 2 i , λ ( p ) = 3.881 and log ( α ( p ) ) = 1.220 ,
where z 1 i p Logis ( 0 , 1 ) and z 2 i p Bern ( 0.6 ) . Here, Logis and Bern denote the standard logistic distribution and the Bernoulli distribution, respectively. The code for reproducing the simulations is available from the following repository https://github.com/isaaccortes1989/RGTG/tree/main/Simulation%20Study (accessed on 5 April 2024).
Figure A1, Figure A2, Figure A3, Figure A4 and Figure A5 in Appendix A.1 display the standard deviation (SD), mean of SE, RMSE, mean of the relative bias (RB) and CP of the 95% asymptotic confidence intervals of the estimates, under different sample sizes. It can be observed in Figure A1, Figure A2, Figure A3 and Figure A4 that the first four measures decrease as the sample size n increases, as expected in standard asymptotic theory. Finally, the CPs in Figure A5 indicate convergence to the nominal values.

6. Applications

In this section, we present two applications to real datasets to illustrate the performance of the GTG model and RGTG quantile regression model. The first application is related to the beta-carotene levels in a certain sample of persons. The second application involves explaining Mexican American individuals’ body mass index (BMI) in terms of their waist circumference and sex. The BMI is calculated by dividing a person’s weight in kilograms by the square of their height in meters, Kg / m 2 . Furthermore, it is a tool that health organizations use to monitor and plan public health programs.

6.1. Without Covariates

We considered the retinol plasma database available at http://lib.stat.cmu.edu/datasets/Plasma_Retinol (accessed on 5 April 2024), which focuses on understanding the determinants of the plasma levels of retinol, beta-carotene and other carotenoids. The main variable of interest is BETADIET, which represents the amount of beta-carotene consumed daily by each individual in micrograms (mcg). The importance of analyzing this variable lies in its direct relationship with the plasma levels of beta-carotene (BETAPLASMA), a key nutrient with antioxidant properties and a precursor of vitamin A. Table 4 presents a summary statistics of BETADIET, in addition to the skewness ( b 1 ) and kurtosis ( b 2 ).
It can be seen from Figure 5 that the GTG has a better fit compared to the GT, Weibull and Slash truncation positive normal (STPN) (see Gómez et al. [22]) models, in addition to a good behavior of the fitted GTG cdf compared to the empirical cdf. Based on the Akaike information criterion (AIC) (see Akaike [23]) and Bayesian information criterion (BIC) (see Schwarz [24]) given in Table 5, we also see that the GTG model is preferred (among the fitted models) for this dataset.

6.2. With Covariates

In this application, we fit the quantile regression model to a dataset provided by Cortés et al. [25]. The dataset comprises the body mass index (BMI) measured in Kg / m 2 , waist circumference (Waist) in centimeters, age in years and Sex (1 for female and 0 for male) of 1743 individuals who self-identified as Mexican American in the National Health and Nutrition Examination Survey (NHANES) conducted between 2017 and 2018.
Here, we assumed that the BMI follows an RGTG distribution, denoted as Z i R G T G ( p ; ρ i p , α p , λ p ) . Accordingly, we propose the following structure for modeling:
log ( ρ i ( p ) ) = β 1 ( p ) + β 2 ( p ) W a i s t i + β 3 ( p ) G e n d e r i and log ( α ( p ) ) = ν ( p ) ,
where β 1 ( p ) , β 2 ( p ) , β 3 ( p ) , λ ( p ) and ν ( p ) are the parameters used for the estimation, for i = 1 , , 1743 and p { 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 } .
The parameter estimates and their corresponding standard errors are presented in Table 6, where it is observed that all estimates are statistically significant at all quantiles p. In relation to the interpretations of the regression coefficients, we can consider the following interpretations. We will present the interpretations in relation to the median, but for this problem, they are also valid for the other quantiles considered:
  • The estimated median BMI for females, assuming a waist circumference of 0, is exp ( 2.076 ) 7.973 .
  • As exp ( β ^ 2 ( p ) ) = exp ( 0.013 ) 1.013 , then for a given gender the estimated median, the BMI is increased by 1.013% for each unit that increases the waist circumference.
  • As exp ( β ^ 3 ( p ) ) = exp ( 0.035 ) 1.036 , then for a given waist circumference, the estimated median BMI for females is increased by approximately 1.036 times in comparison with that for males.
  • λ ^ ( p ) = 3.504 and exp ( ν ^ ( p ) ) = exp ( 1.309 ) 3.702 are the estimates of the shape parameters associates with the median BMI.
Additionally, to assess the model’s adequacy, we calculated the normalized quantile residuals (NQRs) along with their respective envelopes; see Dunn and Smyth [26]. These can be observed graphically in Figure A6 of Appendix A.2, where observation #1267 is highlighted as atypical. Also, we obtained the measures of a generalized Cook’s distance and likelihood displacement for p = 0.9 ; see Figure 6. Here, we highlight that the potentially influential observations are #264, #486, #516, #1267 and #1299.
Table 7 presents a classification (https://www.nhlbi.nih.gov/health/educational/lose_wt/BMI/bmi_dis.htm (accessed on 5 April 2024)) of overweight and obesity by sex, BMI, waist circumference, and risk of diseases (type 2 diabetes, hypertension, and cardiovascular disease) for the highlighted observations in the residuals, generalized Cook’s distance and likelihood displacement, considering p = 0.9 (obese individuals). We observe from the table that observations #264 and #1267 indicate a very poor health status, emphasizing the importance of their detection.

7. Final Discussion

In this work, we proposed a new distribution called generalized truncated Gumbel. The model has quite a few interesting properties, mainly associated with having a cumulative distribution function and a quantile function in closed form. For this reason, an extension of the model was proposed to be used in the context of quantile regression. Future extensions of the model are directed in the context of random effects.

Author Contributions

Conceptualization, H.J.G., K.I.S., D.A. and I.E.C.; formal analysis, K.I.S., D.A., D.I.G. and T.M.M.; funding acquisition, H.J.G. and D.A.; investigation, H.J.G., K.I.S. and I.E.C.; methodology, H.J.G., K.I.S. and D.I.G.; software, I.E.C.; writing—original draft, H.J.G., K.I.S., D.A., D.I.G. and T.M.M.; writing—review and editing, D.A., D.I.G. and T.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data are available for Applications 1 and 2. (http://lib.stat.cmu.edu/datasets/Plasma_Retinol, https://www.nhlbi.nih.gov/health/educational/lose_wt/BMI/bmi_dis.htm) (accessed on 5 April 2024).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Appendix A.1. Behaviors of SD, SE, RMSE, RB and CP

Figure A1. SDs of the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Figure A1. SDs of the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Mathematics 12 01762 g0a1
Figure A2. Means of the asymptotic SE on the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Figure A2. Means of the asymptotic SE on the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Mathematics 12 01762 g0a2
Figure A3. RMSEs of the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Figure A3. RMSEs of the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Mathematics 12 01762 g0a3
Figure A4. Means of the RB on the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression under different sample sizes.
Figure A4. Means of the RB on the 3000 estimates of ρ q , λ q and α q obtained in the RGTG quantile regression under different sample sizes.
Mathematics 12 01762 g0a4
Figure A5. 95% CPs of the ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Figure A5. 95% CPs of the ρ q , λ q and α q obtained in the RGTG quantile regression model under different sample sizes.
Mathematics 12 01762 g0a5

Appendix A.2. QQ-Plots

Figure A6. QQ-plots with envelopes for the NQRs under the model from Equation (16) at different probabilities p.
Figure A6. QQ-plots with envelopes for the NQRs under the model from Equation (16) at different probabilities p.
Mathematics 12 01762 g0a6

References

  1. Bhaskaran, S.; Verma, A.S.; Goupee, A.J.; Bhattacharya, S.; Nejad, A.R.; Shi, W. Comparison of Extreme Wind and Waves Using Different Statistical Methods in 40 Offshore Wind Energy Lease Areas Worldwide. Energies 2023, 16, 6935. [Google Scholar] [CrossRef]
  2. Gurung, B.; Sarkar, K.P.; Singh, K.N.; Lama, A. Modelling annual maximum temperature of India: A distributional approach. Theor. Appl. Climatol. 2021, 145, 979–988. [Google Scholar] [CrossRef]
  3. Purohit, S.U.; Lalit, P.N. European option pricing using Gumbel distribution. Int. J. Financ. Eng. (IJFE) 2022, 9, 2141002. [Google Scholar] [CrossRef]
  4. Li, Y.; Qiao, H.; Yang, A. Reliability Study of Magnesium Oxychloride-Coated Reinforced Concrete Based on Gumbel Distribution. Materials 2023, 16, 2521. [Google Scholar] [CrossRef] [PubMed]
  5. Kang, D.; Ko, K.; Huh, J. Determination of extreme wind values using the Gumbel distribution. Energy 2015, 86, 51–58. [Google Scholar] [CrossRef]
  6. Hossam, E.; Abdulrahman, A.T.; Gemeay, A.M.; Alshammari, N.; Alshawarbeh, E.; Mashaqbah, N.K. A novel extension of Gumbel distribution: Statistical inference with COVID-19 application. Alex. Eng. J. 2022, 61, 8823–8842. [Google Scholar] [CrossRef]
  7. Watthanawisut, A.; Bodhisuwan, B. The Beta Topp-Leone Gumbel Distribution for Modeling the Minimum Flow Data. Thai J. Math. 2020, 165–174. [Google Scholar]
  8. Fayomi, A. A new extended Gumbel distribution: Properties and application. PLoS ONE 2022, 17, e0267142. [Google Scholar] [CrossRef] [PubMed]
  9. Nagode, M. Gumbel mixture modelling for multiple failure data. Reliab. Eng. Syst. Saf. 2023, 230, 108946. [Google Scholar] [CrossRef]
  10. Oseni, B.M.; Okasha, H.M. Gumbel-geometric Distribution: Properties and Applications. Gazi Univ. J. Sci. 2020, 33, 935–941. [Google Scholar] [CrossRef]
  11. Koenker, R.; Bassett, G. Regression Quantiles. Econometrica 1978, 46, 33–50. [Google Scholar] [CrossRef]
  12. Cade, B.S.; Terrell, J.W.; Schroeder, R.L. Estimating Effects of Limiting Factors with Regression Quantiles. Ecology 1999, 80, 311–323. [Google Scholar] [CrossRef]
  13. Koenker, R. Quantile Regression; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  14. Wei, Y.; Pere, A.; Koenker, R.; He, X. Quantile regression methods for reference growth charts. Stat. Med. 2006, 25, 1369–1382. [Google Scholar] [CrossRef] [PubMed]
  15. Gómez-Déniz, E.; Leiva, V.; Calderín-Ojeda, E.; Chesneau, C. A novel claim size distribution based on a Birnbaum-Saunders and gamma mixture capturing extreme values in insurance: Estimation, regression, and applications. Comput. Appl. Math. 2022, 41, 171. [Google Scholar] [CrossRef]
  16. Neamah, M.W.; Qasim, B.A. A New Left Truncated Gumbel Distribution: Properties and Estimation. J. Phys. Conf. Ser. 2021, 1897, 12015. [Google Scholar] [CrossRef]
  17. Cooray, K.; Ananda, M.M. A generalization of the half-normal distribution with applications to lifetime data. Commun. Stat.-Theory Methods 2008, 10, 195–224. [Google Scholar] [CrossRef]
  18. Bonferroni, C.E. Elementi di Statistica Generale; Libreria Seber: Firenze, Italy, 1930. [Google Scholar]
  19. Arcagni, A.; Porro, F. The Graphical Representation of Inequality. Rev. Colomb. Estadística 2014, 37, 419–436. [Google Scholar] [CrossRef]
  20. 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 12 October 2023).
  21. Wang, H.; Chen, W.; Li, B. Large sample properties of maximum likelihood estimator using moving extremes ranked set sampling. J. Korean Stat. Soc. 2024, 53, 398–415. [Google Scholar] [CrossRef]
  22. Gómez, H.J.; Gallardo, D.I.; Santoro, K.I. Slash Truncation Positive Normal Distribution and its Estimation Based on the EM Algorithm. Symmetry 2021, 13, 2164. [Google Scholar] [CrossRef]
  23. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar]
  24. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  25. Cortés, I.E.; De Castro, M.; Gallardo, D.I. A new family of quantile regression models applied to nutritional data. J. Appl. Stat. 2023, 1–21. [Google Scholar] [CrossRef] [PubMed]
  26. Dunn, P.K.; Smyth, G.K. Randomized quantile residuals. J. Comput. Graph. Stat. 1996, 5, 236–244. [Google Scholar] [CrossRef]
Figure 1. Pdf, cdf and hazard function for the G T G ( β = 1 , λ , α ) model with different value combinations for λ and α .
Figure 1. Pdf, cdf and hazard function for the G T G ( β = 1 , λ , α ) model with different value combinations for λ and α .
Mathematics 12 01762 g001
Figure 2. Shape of u ( z ) for β = 1 and some value selections of α and λ .
Figure 2. Shape of u ( z ) for β = 1 and some value selections of α and λ .
Mathematics 12 01762 g002
Figure 3. Plots of the (a) mean, (b) variance, (c) skewness and (d) kurtosis coefficients for GTG ( β = 1 , λ , α ) for α { 1 , 1.5 , 2 } as a function of λ R .
Figure 3. Plots of the (a) mean, (b) variance, (c) skewness and (d) kurtosis coefficients for GTG ( β = 1 , λ , α ) for α { 1 , 1.5 , 2 } as a function of λ R .
Mathematics 12 01762 g003
Figure 4. Bonferroni curves for the G T G ( β , λ , α ) model, with β = 1 and for selected values of λ and α .
Figure 4. Bonferroni curves for the G T G ( β , λ , α ) model, with β = 1 and for selected values of λ and α .
Mathematics 12 01762 g004
Figure 5. (a) GTG, TG, Weibull and STPN models adjusted using the maximum likelihood method for BETADIET. (b) Empirical (black) and fitted GTG (blue) cdf for the BETADIET dataset.
Figure 5. (a) GTG, TG, Weibull and STPN models adjusted using the maximum likelihood method for BETADIET. (b) Empirical (black) and fitted GTG (blue) cdf for the BETADIET dataset.
Mathematics 12 01762 g005
Figure 6. (a) Generalized Cook’s distance and (b) likelihood displacement for the NHANES dataset.
Figure 6. (a) Generalized Cook’s distance and (b) likelihood displacement for the NHANES dataset.
Mathematics 12 01762 g006
Table 1. Skewness coefficient of the G T G ( β = 1 , λ , α ) model for different values of λ and α .
Table 1. Skewness coefficient of the G T G ( β = 1 , λ , α ) model for different values of λ and α .
α
λ 0.511.522.5
−26.44871.95021.03770.1428−0.3964
−1.56.34331.91901.01610.1265−0.4114
05.55701.67820.8441−0.0092−0.5421
1.53.92151.20190.5286−0.2238−0.7331
23.46301.14570.5728−0.0296−0.3987
2.53.12081.13960.66100.2027−0.0336
Table 2. Kurtosis coefficient of the G T G ( β = 1 , λ , α ) model for different values of λ and α .
Table 2. Kurtosis coefficient of the G T G ( β = 1 , λ , α ) model for different values of λ and α .
α
λ 0.511.522.5
−283.42298.70634.29032.72263.0664
−1.580.81818.52704.22972.72013.0884
062.82587.27263.82312.75893.3450
1.533.31575.52173.58973.44254.6563
226.83415.40813.75743.45044.1342
2.522.46255.40013.89293.28713.3145
Table 3. Estimated bias, SE, RMSE and CP for ML estimators in finite samples from the GTG model.
Table 3. Estimated bias, SE, RMSE and CP for ML estimators in finite samples from the GTG model.
True Value n = 150 n = 300 n = 600 n = 1000
β λ α EstimatorbiasSERMSECPbiasSERMSECPbiasSERMSECPbiasSERMSECP
121 β ^ 0.00370.29190.33910.9400.00720.20090.22360.947−0.00060.13810.14510.940−0.00120.10620.10720.956
λ ^ −0.01990.413911.2130.9920.00270.24790.54460.9740.01230.16470.17460.9520.00650.12550.12720.962
α ^ 0.01600.16690.20520.9520.01100.11280.12900.9520.00210.07650.07960.9500.00030.05860.05870.950
2 β ^ −0.00730.14990.15870.972−0.00850.10130.10680.962−0.00200.06950.07200.952−0.00130.05320.05410.950
λ ^ −0.01910.437110.0330.9950.00890.26050.59560.9760.00950.16500.17160.9590.00570.12540.12830.955
α ^ 0.03440.33510.39040.9620.00490.22260.25000.9620.00700.15360.16220.9390.00350.11730.11840.950
31 β ^ −0.00120.42180.41010.928−0.02430.30050.28900.942−0.02940.21380.21700.935−0.00590.16590.17050.943
λ ^ 0.13980.73790.85530.9550.10230.48180.48650.9650.07650.32890.34610.9570.02960.24540.25770.946
α ^ 0.00890.20960.21480.969−0.00760.14360.13900.966−0.01150.10050.10140.947−0.00140.07760.07980.943
2 β ^ −0.03440.22940.22250.946−0.02210.15820.16510.942−0.01500.11030.10940.954−0.00190.08380.08620.931
λ ^ 0.18800.75450.77760.9580.10420.48590.53980.9550.05990.32580.32870.9650.01870.24330.25250.935
α ^ −0.00270.41860.39750.963−0.00970.28790.29130.951−0.01160.20150.19680.9540.00850.15540.16080.931
221 β ^ −0.02490.58930.71020.939−0.00430.40280.46030.9390.00010.27650.29150.948−0.00180.21150.20210.947
λ ^ −0.01040.512811.6820.9950.00180.25260.65240.9710.01310.16490.17310.9550.00610.12500.11950.963
α ^ 0.00680.16810.21910.9490.00540.11270.13460.9520.00320.07670.08070.9510.00160.05850.05500.959
2 β ^ −0.01670.30020.33780.965−0.01040.20120.20310.971−0.00960.13910.14370.954−0.00270.10660.10790.949
λ ^ −0.02770.543611.1010.9940.01840.24220.24410.9840.01550.16500.17080.9640.00660.12560.12900.949
α ^ 0.04610.33830.41600.9530.01130.22320.22780.9660.00070.15290.15760.9400.00260.11740.11770.947
31 β ^ −0.05130.84240.80940.920−0.05780.60300.56860.947−0.01140.42860.41230.958−0.02260.33200.32240.951
λ ^ 0.19620.75720.76380.9550.11020.48440.48060.9620.04110.32190.32600.9690.03510.24590.24210.963
α ^ −0.00310.21000.20150.964−0.00920.14400.13580.960−0.00070.10080.09620.962−0.00370.07750.07550.954
2 β ^ −0.08620.46120.46740.953−0.04500.31830.31480.949−0.03020.22060.22410.955−0.01990.16920.16930.950
λ ^ 0.21770.77210.83100.9620.10530.48600.49690.9560.05860.32570.33950.9590.03530.24620.24930.949
α ^ −0.01630.41750.41280.960−0.00900.28980.28280.949−0.01300.20110.20050.951−0.01040.15500.15520.955
321 β ^ 0.08080.893511.1770.9430.02660.60270.64920.9560.00620.41460.43510.9490.01070.31870.34270.944
λ ^ −0.06490.517911.6780.9920.00320.25560.40530.9840.01010.16460.17350.9540.00280.12550.13370.945
α ^ 0.02920.17090.22970.9600.01160.11280.12600.9660.00440.07670.08170.9420.00380.05880.06340.945
2 β ^ −0.03740.45710.51530.968−0.01520.30380.29850.970−0.01220.20940.21420.957−0.00180.15920.16120.957
λ ^ −0.04200.460712.5090.9910.01950.24310.24440.9850.01470.16550.16940.9660.00380.12520.12880.952
α ^ 0.02850.33970.42350.9550.00650.22330.21790.9650.00120.15350.15610.9540.00510.11710.11840.947
31 β ^ −0.065312.56612.4620.920−0.02680.90290.86710.948−0.02900.64150.63350.945−0.01720.49860.50850.934
λ ^ 0.18420.75090.75660.9600.08140.47590.47280.9710.04920.32300.33450.9620.02790.24550.25510.943
α ^ −0.00320.20850.20640.964−0.00020.14400.13800.959−0.00290.10050.09980.952−0.00150.07770.07870.941
2 β ^ −0.08380.68460.67060.942−0.09340.47910.48280.956−0.03050.32820.33200.951−0.01310.25230.24040.966
λ ^ 0.16530.74700.77410.9490.12900.49130.51760.9640.04350.32200.33130.9550.02290.24400.23620.962
α ^ 0.00900.41900.40550.973−0.02570.28790.28210.954−0.00260.20110.20180.9450.00120.15530.14680.967
Table 4. Descriptive statistics of the amount of beta-carotene consumed daily by each individual in micrograms (BETADIET).
Table 4. Descriptive statistics of the amount of beta-carotene consumed daily by each individual in micrograms (BETADIET).
Datasetn X ¯ S 2 b 1 b 2
BETADIET31521862,172,342 1.606 6.399
Table 5. Estimated parameters and their standard errors (in parentheses) for the GTG, TG, Weibull, and STPN models for the BETADIET dataset. The AIC and BIC criteria are also presented.
Table 5. Estimated parameters and their standard errors (in parentheses) for the GTG, TG, Weibull, and STPN models for the BETADIET dataset. The AIC and BIC criteria are also presented.
ParametersGTG (SE)TG (SE)WEI (SE)STPN (SE)
β ^ 290.83 (126.04)1041.09 (56.53)1.61 (0.07)660.57 (95.37)
λ ^ 2.61 (0.34)1.46 (0.10)2454.96 (91.19)1.97 (0.22)
α ^ 0.60 (0.08)--2.46 (0.35)
LLF−2672.896−2682.618−2686.675−2685.109
AIC5351.7925369.2365377.3515376.219
BIC5357.2975376.7425384.8565387.476
Table 6. ML estimates of the parameters and their corresponding standard errors (in parentheses) for the RGTG quantile regression model.
Table 6. ML estimates of the parameters and their corresponding standard errors (in parentheses) for the RGTG quantile regression model.
p β ^ 1 ( p ) β ^ 2 ( p ) β ^ 3 ( p ) λ ^ ( p ) ν ^ ( p )
0.11.975 (0.0082)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.22.009 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.32.034 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.42.055 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.52.076 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.62.096 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.72.118 (0.0081)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.82.145 (0.0082)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
0.92.183 (0.0084)0.013 (0.0001)0.035 (0.0039)3.504 (0.2444)1.309 (0.0675)
Table 7. Classification of overweight and obesity by BMI, waist circumference and associated disease risks.
Table 7. Classification of overweight and obesity by BMI, waist circumference and associated disease risks.
CaseSexBMIClassificationObesityWaistDisease
of OverweightClassRisk
264Female48.45Extreme obesityIII115.8Extremely high
486Female37.21ObesityII97.8Very high
516Female18.04Underweight-79.8-
1267Female71.72Extreme obesityIII132.3Extremely high
1299Male23.33Normal-101.8-
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, H.J.; Santoro, K.I.; Ayma, D.; Cortés, I.E.; Gallardo, D.I.; Magalhães, T.M. A New Generalization of the Truncated Gumbel Distribution with Quantile Regression and Applications. Mathematics 2024, 12, 1762. https://doi.org/10.3390/math12111762

AMA Style

Gómez HJ, Santoro KI, Ayma D, Cortés IE, Gallardo DI, Magalhães TM. A New Generalization of the Truncated Gumbel Distribution with Quantile Regression and Applications. Mathematics. 2024; 12(11):1762. https://doi.org/10.3390/math12111762

Chicago/Turabian Style

Gómez, Héctor J., Karol I. Santoro, Diego Ayma, Isaac E. Cortés, Diego I. Gallardo, and Tiago M. Magalhães. 2024. "A New Generalization of the Truncated Gumbel Distribution with Quantile Regression and Applications" Mathematics 12, no. 11: 1762. https://doi.org/10.3390/math12111762

APA Style

Gómez, H. J., Santoro, K. I., Ayma, D., Cortés, I. E., Gallardo, D. I., & Magalhães, T. M. (2024). A New Generalization of the Truncated Gumbel Distribution with Quantile Regression and Applications. Mathematics, 12(11), 1762. https://doi.org/10.3390/math12111762

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