Next Article in Journal
MHD Laminar Boundary Layer Flow of a Jeffrey Fluid Past a Vertical Plate Influenced by Viscous Dissipation and a Heat Source/Sink
Previous Article in Journal
A Hybrid Approach of VIKOR and Bi-Objective Decision Model for Emergency Shelter Location–Allocation to Respond to Earthquakes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Extension of the Truncated-Exponential Skew- Normal Distribution

by
Pilar A. Rivera
1,
Diego I. Gallardo
2,
Osvaldo Venegas
3,*,
Marcelo Bourguignon
4 and
Héctor W. Gómez
1
1
Departamento de Matemáticas, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta 1240000, Chile
2
Departamento de Matemática, Facultad de Ingeniería, Universidad de Atacama, Copiapó 1530000, Chile
3
Departamento de Ciencias Matemáticas y Físicas, Facultad de Ingeniería, Universidad Católica de Temuco, Temuco 4780000, Chile
4
Departamento de Estatística, Universidade Federal do Rio Grande do Norte, Natal 59078-970, Brazil
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(16), 1894; https://doi.org/10.3390/math9161894
Submission received: 24 June 2021 / Revised: 15 July 2021 / Accepted: 27 July 2021 / Published: 9 August 2021
(This article belongs to the Section Probability and Statistics)

Abstract

:
In the paper, we present an extension of the truncated-exponential skew-normal (TESN) distribution. This distribution is defined as the quotient of two independent random variables whose distributions are the TESN distribution and the beta distribution with shape parameters q and 1, respectively. The resulting distribution has a more flexible coefficient of kurtosis. We studied the general probability density function (pdf) of this distribution, its survival and hazard functions, some of its properties, moments and inference by the maximum likelihood method. We carried out a simulation and applied the methodology to a real dataset.

1. Introduction

The Slash (S) distribution is a generalization of the normal model. Its stochastic representation is given by
S = Z U 1 q ,
where Z N ( 0 , 1 ) is independent of U U ( 0 , 1 ) and q > 0 .
q = 1 represents the canonical Slash model and the standard normal model is obtained for q . The pdf of the canonical S distribution is
p ( x ) = ϕ ( 0 ) ϕ ( x ) x 2 x 0 1 2 ϕ ( 0 ) x = 0
where ϕ ( · ) represents the pdf of the standard normal model (see Johnson et al. [1]). This distribution is characterized by having heavier tails than normal distribution, i.e., it has greater kurtosis.
Properties of the S distribution are discussed by Rogers and Tukey [2] and Mosteller and Tukey [3]. The ML parameters for location and scale in the S model are discussed in Kadafar [4]. Wang et al. [5] studied a multivariate version of the S distribution and a multivariate skew version. Gómez et al. [6] extended the S distribution using the family of univariate and multivariate elliptical distributions also was extended by using the S model in Gómez et al. [7].
Nadarajah et al. [8] proposed the idea of constructing biased distributions, motivated by Azzalini [9], including asymmetry in these. A unified approach for the construction of models of this kind is given in Ferreira and Steel [10]. If X is a symmetrical random variable around zero with pdf f X ( · ) and cumulative distribution function (cdf) F X ( · ) , the new random variable Y is defined with pdf given by:
f Y ( y ) = f X ( y ) ω ( F X ( y ) ) , y R ,
with ω ( · ) denoting a pdf in the interval ( 0 , 1 ) . Then Y is a skew version of the variable X. The most commonly-used versions of (2) are the skew distributions proposed by Azzalini [9] in the form:
f Y ( y ) = 2 f X ( y ) F X ( λ y ) , y R , λ R .
The skew-normal (SN) model is obtained as a particular case of (3) considering f X ( · ) = ϕ ( · ) and F X ( · ) = Φ ( · ) , where Φ ( · ) denotes the cdf of the standard normal model.
In the present paper, we extend the TESN model introduced by Nadarajah et al. [8], based on the Slash methodology. The pdf of the TESN distribution is given by:
f ( x ; λ ) = λ 1 e λ ϕ ( x ) e { λ Φ ( x ) } ,
where x , λ R . Hereafter, we use the notation X T E S N ( λ ) to indicate that X is a random variable following a TESN distribution. According to Barreto-Souza and Simas [11], the distribution presents different behavior for a large | λ | , suggesting that this is a rich class of distributions. Furthermore, the parameter λ can be interpreted as a concentration parameter. Figure 1 shows the graph of the TESN pdf function with variations of the parameter λ .
The extension of this model is based on the quotient of two independent random variables, with TESN distribution and a power of the uniform distribution ( 0 , 1 ) , respectively, obtaining a distribution with a more flexible coefficient of kurtosis and so generating an appropriate model for fitting data. In practical terms, this generalization is based on the search for distributions that are more flexible, which may provide a “better fit” than the TESN distribution. For example, see the works by Gomes et al. [12], Maurya and Nadarajah [13] and the references therein.
The article is organized as follows. In Section 2 we study the representation, pdf, properties, and graphs. In Section 3, we present a Monte Carlo simulation experiment to evaluate the maximum likelihood estimates of the model parameters in Section 2. In Section 4, we provide an application of the proposed distribution. In Section 5, we conclude with some final comments.

2. Incorporating Kurtosis

In this section, we introduce a new extension of the TESN distribution. We studied its pdf, survival and hazard functions, moments, location and scale parameters, and their log-likelihood equations.

2.1. Representation

Following the representation of the Slash distribution, the representation of this new distribution is given by the following definition:
Definition 1.
A random variable (r.v.) Z has a Slashed Truncated Exponential Skew-Normal (STESN) distribution, denoted by Z S T E S N ( λ , q ) if it is represented by:
Z = X Y 1 / q ,
where X T E S N ( λ ) and Y U ( 0 , 1 ) are independent variables, λ R , q > 0 .

2.2. Probability Density Function

The following proposition shows the pdf function for the STESN distribution, generated using the stochastic representation given in (5) and using the Jacobian method for transforming the r.v.:
Proposition 1.
If Z S T E S N ( λ , q ) the pdf of Z is given by:
f Z ( z ; λ , q ) = λ · q 1 e λ R ( z ; λ , q ) ,
where z R , λ R , q > 0 and R ( z ; λ , q ) = 0 1 w q ϕ ( z w ) e { λ Φ ( z w ) } d w .
Proof. 
The pdf is generated using the representation given in (5). Using the Jacobian method for transforming the r.v. we obtain X = Z · W ; W = Y 1 / q , calculating the Jacobian, we obtain:
J ( z , w ) =   w z 0 q w q 1   = q w q ,
replacing the joint pdf f Z , W :
f Z , W ( z , w ) = q w q f X ( z w ) f Y ( w q ) = q w q λ 1 e λ ϕ ( z w ) e { λ Φ ( z w ) } ,
where z R , 0 < w < 1 , and λ R . Hence, marginalizing with respect to variable W, we obtain:
f Z ( z ) = λ · q 1 e λ R ( z ; λ , q ) ,
where z R , λ R , q > 0 and:
R ( z ; λ , q ) = 0 1 w q ϕ ( z w ) e { λ Φ ( z w ) } d w .
Proposition 2.
Let Z S T E S N ( λ , q ) . If q the r.v. Z converges in law to the r.v. X T E S N ( λ ) .
Proof. 
Let Z S T E S N ( λ , q ) and Z = X Y 1 / q given by (5), where Y U ( 0 , 1 ) T = Y 1 / q B e t a ( q , 1 ) . We obtain:
E [ ( T 1 ) 2 ] = 2 ( q + 1 ) 2 ( q + 2 ) ,
where if q E [ ( T 1 ) 2 ] 0 T P 1 (see Lehmann [14]) here P denote converges in probability. Applying Slutsky’s Lemma to Z = X V , we have:
Z L X T E S N ( λ ) , q ,
where L denote converges in law, i.e., for increasing values of q, the r.v Z converges in law to a T E S N ( λ ) distribution. □
Remark 1.
The above proposition implies that if q then the pdf of the STESN distribution approaches the pdf of a TESN distribution.
Figure 2a,b show the graphs for the pdf of this model for some parameter values.

2.3. Reliability Analysis

For the study of failure times, we need to consider a time y, where we have y 0 . Therefore, in our model we study the case of non-negative variables Y = exp ( Z ) where Z S T E S N ( λ , q ) , thus, the model must be transformed to obtain the following pdf:
f Y ( y ) = q · λ y · ( 1 e λ ) · R log ( y ) ; λ , q ,
where y > 0 , q > 0 , λ R , and R ( log ( y ) ; λ , q ) = 0 1 w q ϕ ( w log ( y ) ) e { λ Φ ( w log ( y ) ) } d w . Figure 2c shows the graphs of the pdf for different parameter values. Once the transformation is complete, we obtain the survival and hazard functions. The survival function is defined as the probability that a subject does not experience the event of interest before a moment t, and in our model is given by:
S ( y ) = 1 1 1 e λ 0 y κ ( t ) d t ,
where κ ( t ) = q · λ t · R ( log ( t ) ; λ , q ) . It also gives the hazard function defined as the probability of failure during a time interval given in our model by:
λ ( y ) = κ ( y ) 1 e λ 0 y κ ( t ) d t ,
where κ ( a ) = q · λ a · R ( log ( a ) ; λ , q ) . Figure 3a,b show the graphs of survival and hazard functions respectively, for different parameter values.

2.4. Moments

Let Z be a r.v. where Z S T E S N ( λ , q ) , the r-th moment for the variable is given by the following proposition.
Proposition 3.
Using the representation (5), the r-th moment of the r.v. Z is:
E [ Z r ] = λ · q ( q r ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 r ] ,
where q > r and X k : n the k-th order statistic of a random variable with distribution N ( 0 , 1 ) .
Proof. 
The r-th moments of Z can be calculated as:
E [ Z r ] = E X Y 1 q r   = E X r · Y r q   = E [ X r ] · E Y r q ,
where E [ X r ] is the r-th moment for the model proposed by Nadarajah et al. [8] given by:
E [ X ] = λ 1 e λ k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 r ] ,
where X k : n the k-th order statistic of a random variable with distribution N ( 0 , 1 ) and E Y r q the r-th moment of the r.v. Y 1 q .
E [ Y r / q ] = q q r .
Therefore the r-th moment for the variable is given by:
E [ Z r ] = λ · q ( q r ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 r ] .
Using this proposition, the first four moments of the r.v. Z are given in the following corollary.
Corollary 1.
From the r -th moment of the r.v. Z S T E S N ( λ , q ) represented by (10), we obtain the first four moments of the variable, given by:
E [ Z ] = λ · q ( q 1 ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 ] , q > 1 ,
E [ Z 2 ] = λ · q ( q 2 ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 2 ] , q > 2 ,
E [ Z 3 ] = λ · q ( q 3 ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 3 ] , q > 3 ,
E [ Z 4 ] = λ · q ( q 4 ) ( 1 e λ ) k = 0 ( λ ) k ( k + 1 ) ! E [ X k + 1 : k + 1 4 ] , q > 4 .

2.5. Incorporation of Parameters

To produce a more flexible distribution, we will extend this model to location μ and scale σ parameters as X = μ + σ Z , where Z S T E S N ( λ , q ) obtaining the following proposition.
Proposition 4.
If X S T E S N ( λ , q , μ , σ ) the pdf of X is given by:
f X ( x ) = q · λ σ ( 1 e λ ) R x μ σ ; λ , q , μ , σ ,
where x R ; λ R ; q > 0 ; μ R ; σ > 0 ; and R ( · ) is presented in (6).

2.6. Log Likelihood Equations

Let x 1 , , x n be a random sample of the r.v. X with S T E S N ( λ , q , μ , σ ) distribution, the log-likelihood function can be written as:
l ( θ ; X ) = n · log ( q ) + n · log ( λ ) n · log ( σ ) n · log ( 1 e λ ) + i = 1 n log ρ ( x i ) ,
where θ = ( λ , q , μ , σ ) and ρ ( x i ) = R x i μ σ ; λ , q , μ , σ .
For each parameter we have the following likelihood equations:
n λ n · e λ 1 e λ + i = 1 n ρ 1 ( x i ) ρ ( x i ) = 0 n q + i = 1 n ρ 2 ( x i ) ρ ( x i ) = 0 1 σ i = 1 n ρ 3 ( x i ) ρ ( x i ) = 0 n σ 1 σ 2 i = 1 n ( x i μ ) · ρ 4 ( x i ) ρ ( x i ) = 0 ,
where ρ ( x i ) = R x i μ σ ; λ , q , μ , σ . Furthermore, ρ 1 ( x i ) = d ρ ( x i ) d λ , ρ 2 ( x i ) = d ρ ( x i ) d q , ρ 3 ( x i ) = d ρ ( x i ) d μ , and ρ 4 ( x i ) = d ρ ( x i ) d σ .
As can be observed, this system can only be resolved by iterative procedures such as Newton–Raphson. As an alternative, it is also possible is to use the optim routine implemented in R software [15]. Standard errors for parameters can be estimated using the hessian matrix of the log-likelihood function, which can be estimated, for instance, using the pracma package (see Borchers [16]).

2.7. STESN or TESN Model?

In order to decide between the STESN and TESN models, we can use the traditional Akaike (AIC, Akaike [17]) and Schwarz (BIC, Schwartz [18]) criteria. As the TESN model corresponds to the STESN model with q + , we can use the likelihood ratio test (LRT) to decide between the two models considering H 0 : q = + (TESN model) versus H 1 : q < (STENS model). This is a problem where the null hypothesis is exactly on the boundary of the parameter space. This kind of problem was first discussed in Chernoff [19]. This problem is also presented, for instance, in a random effects model when we are interested in testing if the variance of such random effects is zero (see Stram and Lee [20] and Gallardo et al. [21]) or in a cure rate model when we are interested in testing the presence of cured individuals in the population (see Maller and Zhou [22]). In this particular case, the statistic for the LRT, say d n , does not converge asymptotically to the usual χ ( 1 ) 2 distribution, i.e., the chi-squared distribution with 1 degree of freedom, but converges to 1 2 χ ( 0 ) 2 + 1 2 χ ( 1 ) 2 , i.e., a 50-50 mixture between a point mass χ ( 0 ) 2 in 0 and χ ( 1 ) 2 distribution.

3. Simulation Study

In this section, we will study the behavior of ML estimators in finite samples, verifying empirically whether these estimators have desirable properties (unbiased, asymptotically efficient, verification of the normal asymptotic distribution of ML estimators).
The random variables of the TESN distribution and the Beta distribution were generated to obtain our new variable with pdf shown in (15). The initial values used for optimization were obtained by a sequence of values which maximized this function. In this sequence, λ takes values between 3 and 3, q between 1 and 5, μ between 2 and 2, and σ between 2 and 10. This process was repeated 5000 times with sample size n = 25 , n = 50 , n = 100 and 200 for different combinations of parameters. Table 1 presents the empirical bias, the standard errors (SE), root of the mean squared error (RMSE), and 95% coverage probabilities (CP) for the estimators of the parameters of the STESN distribution with different combinations of parameters and sample sizes. From those tables, notice that the biases SE and RMSE decrease as the sample size increases, suggesting that estimators are consistent. Furthermore, the asymptotic confidence intervals have an empirical CP differing from the nominal values, especially when the sample size is small. However, we observe that the asymptotic confidence intervals converge to the nominal values when the sample size is increased. Figure 4 shows the estimated pdf for the ML estimators of μ , σ , λ , and q for two combinations of parameters, showing graphically that the skewness of the estimators disappears progressively when the sample sizes increases. We also note that the distribution of the estimators for λ and q are more asymmetric than the distribution of the estimators for μ and σ , especially in small sample sizes.

4. Application to a Data Set

In this section, we will present a real data application to illustrate the STESN model compared with other models discussed in the literature. These data were presented by Barlow et al. [23] and represent the fatigue fracture life of Kevlar 373/epoxi subjected to a constant pressure of 90 % stress until they all fail. To obtain the parameter estimations, the optim command was used and its estimation errors were calculated by the Hessian matrix, both in R software. Codes are available as supplementary material.
Table 2 shows a summary of the dataset, including the sample size n, the mean X ¯ , the standard deviation S, the asymmetry coefficient b 1 , the kurtosis coefficient b 2 , the minimum min ( X ) , and the maximum max ( X ) . A high kurtosis value is observed.
Table 3 shows the results of the fit; the TESN distribution was compared with the STESN distribution by AIC and BIC criteria. It is concluded that the distribution which achieves the best fit for this dataset is the STESN distribution, since it presents a lower value in the criteria. Furthermore, Table 3 provides the Kolmogorov–Smirnov statistic (KSS), a formal goodness-of-fit test to verify which distribution gives a better fit for these data. Small values of this statistic suggests a better fit. Thus, according to the Kolmogorov–Smirnov test, the STESN model fits the current data better than the TESN model.
Figure 5a,b present a histogram of the data with the densities fitted for the data set and Figure 6a,b present the QQ-plot of the densities fitted for the dataset, showing the good fit given by the new distribution.
In our problem, the observed statistic d n for the LRT to decide between the TESN and STESN models, discussed in Section 2.7, is d n = 23.59 with an associated p-value < 0.001 . Therefore, the H 0 is rejected under any usual level of significance and the STESN model is preferred over the TESN model.

5. Final Comments

In this paper, we introduced an extension of the TESN distribution from which we obtained a distribution that showed greater flexibility in the coefficient of kurtosis. Some mathematical properties of the new distribution were studied. Note that the formulae derived easily implemented in different softwares. Inference was implemented based on the ML approach, and its performance was assessed by Monte Carlo simulations using R software. An application to a real dataset showed that the new model produced a better fit than the TESN model. This application demonstrated the practical importance of the new model, and also showed the advantage of STESN over TESN. We hope this new distribution may attract wider applications.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/math9161894/s1.

Author Contributions

Conceptualization, P.A.R. and D.I.G.; methodology, H.W.G.; software, D.I.G.; validation, P.A.R., D.I.G., and H.W.G.; formal analysis, O.V.; investigation, O.V. and M.B.; writing—original draft preparation, P.A.R.; writing—review and editing, O.V. and M.B.; funding acquisition, M.B. All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The research of Pilar A. Rivera and Héctor W. Gómez was supported by SEMILLERO UA-2021 (Chile). The research of O. Venegas was supported by Vicerrectoría de Investigación y Postgrado of the Universidad Católica de Temuco, Projecto interno FEQUIP 2019-INRN-03.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in Section 4 can be obtained from the corresponding reference.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Johnson, N.L.; Kotz, S.; Balakrishnan, N. Continuous Univariate Distributions, 2nd ed.; Wiley Series in Probability and Statistics; Wiley: New York, NY, USA, 1995. [Google Scholar]
  2. Rogers, W.H.; Tukey, J.W. Understanding some long-tailed symmetrical distributions. Stat. Neerl. 1972, 26, 211–226. [Google Scholar] [CrossRef]
  3. Mosteller, F.; Tukey, J.W. Data Analysis and Regression: A Second Course in Statistics; Addison-Wesley: Reading, MA, USA, 1977. [Google Scholar]
  4. Kadafar, K. A biweight approach to the one-sample problem. J. Am. Stat. Assoc. 1982, 77, 416–424. [Google Scholar] [CrossRef]
  5. Wang, J.; Boyer, J.; Genton, M.G. A skew-symmetric representation of multivariate distributions. Stat. Sin. 2004, 14, 1259–1270. [Google Scholar]
  6. Gómez, H.W.; Quintana, F.A.; Torres, F.J. New family of slash-distributions with elliptical contours. Stat. Probab. Lett. 2007, 77, 717–725. [Google Scholar] [CrossRef]
  7. Gómez, H.W.; Olivares-Pacheco, J.F.; Bolfarine, H. An extension of the generalized Birnbaum-Saunders distribution. Stat. Probab. Lett. 2009, 79, 331–338. [Google Scholar] [CrossRef]
  8. Nadarajah, S.; Nassiri, V.; Mohammadpour, A. Truncated-exponential skew-symmetric distributions. Statistics 2014, 48, 872–895. [Google Scholar] [CrossRef]
  9. Azzalini, A. A Class of Distributions Which Includes the Normal Ones. Scand. J. Stat. 1985, 12, 171–178. [Google Scholar]
  10. Ferreira, J.T.A.S.; Steel, M.F.J. A constructive representation of univariate skewed distributions. J. Am. Stat. Assoc. 2006, 101, 823–829. [Google Scholar] [CrossRef] [Green Version]
  11. Barreto-Souza, W.; Simas, A.B. The exp- G family of probability distributions. Braz. J. Probab. Stat. 2013, 27, 84–109. [Google Scholar] [CrossRef]
  12. Gomes, A.E.; Da-Silva, C.Q.; Cordeiro, G.M. The Exponentiated G Poisson Model. Commun. Stat.-Theory Methods 2015, 44, 4217–4240. [Google Scholar] [CrossRef]
  13. Maurya, S.K.; Nadarajah, S. Poisson Generated Family of Distributions: A Review. Sankhya B 2020, 1–57. [Google Scholar] [CrossRef]
  14. Lehmann, E.L. Elements of Large-Sample Theory; Springer: New York, NY, USA, 1999. [Google Scholar]
  15. R Development Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  16. Borchers, H.W. Pracma: Practical Numerical Math Functions, R package version 2.3.3; 2021. Available online: https://CRAN.R-project.org/package=pracma (accessed on 24 June 2021).
  17. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  18. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  19. Chernoff, H. On the distribution of the likelihood ratio. Ann. Stat. 1954, 54, 573–578. [Google Scholar] [CrossRef]
  20. Stram, D.O.; Lee, J.W. Variance components testing in the longitudinal mixed effects model. Biometrics 1994, 50, 1171–1177. [Google Scholar] [CrossRef] [PubMed]
  21. Gallardo, D.I.; Bolfarine, H.; Pedroso-de-Lima, A.C. A clustering cure rate model with application to a sealant study. J. Appl. Stat. 2017, 44, 2949–2962. [Google Scholar] [CrossRef]
  22. Maller, R.; Zhou, S. Testing for the Presence of Immune or Cured Individuals in Censored Survival Data. Biometrics 1995, 51, 1197–1205. [Google Scholar] [CrossRef] [PubMed]
  23. Barlow, R.E.; Toland, R.H.; Freeman, T. A Bayesian analysis of stress rupture life of Kevlar 49/epoxy spherical pressure vessels. In Procedings Conference on Applications of Statistics; Marcel Dekker: New York, NY, USA, 1984. [Google Scholar]
Figure 1. TESN pdf for different values of λ .
Figure 1. TESN pdf for different values of λ .
Mathematics 09 01894 g001
Figure 2. STESN for values of (a) λ = 2 ; (b) λ = 2 ; and (c) λ = 7 for different values of q.
Figure 2. STESN for values of (a) λ = 2 ; (b) λ = 2 ; and (c) λ = 7 for different values of q.
Mathematics 09 01894 g002
Figure 3. (a) Survival function and (b) hazard function for log-STESN model with different combinations of values for λ and q.
Figure 3. (a) Survival function and (b) hazard function for log-STESN model with different combinations of values for λ and q.
Mathematics 09 01894 g003
Figure 4. Estimated pdf for the ML estimators of μ , σ , λ , and q in the TESN distribution for: μ = 0 , σ = 1 , λ = 3 , q = 1 (upper panels), and μ = 10 , σ = 16 , λ = 1 , q = 3 (lower panels).
Figure 4. Estimated pdf for the ML estimators of μ , σ , λ , and q in the TESN distribution for: μ = 0 , σ = 1 , λ = 3 , q = 1 (upper panels), and μ = 10 , σ = 16 , λ = 1 , q = 3 (lower panels).
Mathematics 09 01894 g004
Figure 5. Estimated pdf for STESN and TESN for kevlar data set (a) and a zoom for the right tail (b).
Figure 5. Estimated pdf for STESN and TESN for kevlar data set (a) and a zoom for the right tail (b).
Mathematics 09 01894 g005
Figure 6. QQ-plot for the (a) TESN and (b) STESN distributions for the dataset.
Figure 6. QQ-plot for the (a) TESN and (b) STESN distributions for the dataset.
Mathematics 09 01894 g006
Table 1. Empirical bias, SE, RMSE, and 95% CP for the ML estimators of μ , σ , λ , and q with different combinations of parameters and sample sizes.
Table 1. Empirical bias, SE, RMSE, and 95% CP for the ML estimators of μ , σ , λ , and q with different combinations of parameters and sample sizes.
True Values n = 25 n = 50 n = 100 n = 200
μ σ λ q biasSERMSECPbiasSERMSECPbiasSERMSECPbiasSERMSECP
01−31 μ 0.0531.1401.1410.9030.0310.8820.8850.9090.0130.6620.6750.9120.0090.4490.4580.939
σ 0.0840.4890.4960.9720.0690.3670.3730.9650.0520.2780.2870.9610.0420.1970.2030.959
λ −0.1762.3362.3360.887−0.0931.9141.9360.892−0.0531.3661.3980.906−0.0420.8920.9100.921
q0.1950.5870.6180.9730.1190.3500.3700.9720.0820.2200.2350.9650.0570.1510.1610.961
01−11 μ −0.0991.4731.4730.918−0.0781.0451.0490.928−0.0510.7040.7110.935−0.0330.4730.4770.948
σ 0.1590.5320.5550.9700.1110.3720.3880.9650.0750.2520.2630.9580.0440.1790.1840.953
λ −0.1542.3982.3990.887−0.1171.8101.8220.889−0.0831.1831.1960.907−0.0600.7650.7700.915
q0.2140.5550.5940.9620.1630.4790.5060.9600.0800.2290.2430.9580.0410.1470.1520.953
0131 μ −0.1001.1661.1700.901−0.0910.9070.9130.928−0.0770.6150.6240.936−0.0470.4480.4580.945
σ 0.0940.5120.5140.9740.0760.3940.4010.9680.0560.2540.2610.9610.0530.1940.2010.957
λ 0.2512.3002.2990.8710.1921.8431.8750.8960.1671.2951.3220.9170.0940.8640.8870.925
q0.1680.6570.6780.9650.1150.3760.3930.9610.0710.2200.2310.9600.0450.1480.1540.959
01−22 μ −0.0611.1401.1390.905−0.0331.0971.0970.918−0.0310.9410.9500.929−0.0100.6640.6720.932
σ 0.0530.3860.3890.9860.0420.3430.3570.9790.0340.2970.3180.9710.0250.1890.2000.964
λ −0.1922.8752.8800.894−0.1512.7782.7830.901−0.1222.3182.3550.915−0.0971.6251.6520.921
q0.3621.3011.3500.9710.2781.2541.3410.9690.2261.0081.0940.9630.1800.5100.5410.959
0132 μ −0.2051.0041.0240.899−0.0790.9150.9180.910−0.0610.7580.7620.928−0.0440.5800.5810.936
σ 0.0790.3890.3890.9710.0530.3130.3140.9680.0410.2640.2710.9610.0300.2010.2050.959
λ 0.2782.6552.6680.8960.2122.5702.5690.9060.1442.1002.1270.9170.1001.5661.5780.929
q0.2131.2561.2730.9650.1921.0351.0750.9610.1190.8800.9260.9590.0990.5070.5290.958
0122 μ 0.1641.1711.1910.9040.1151.1141.1150.9120.0790.9010.9040.9390.0520.6890.6930.945
σ 0.0990.3690.3690.9860.0820.3350.3500.9820.0740.2800.2990.9710.0400.2040.2150.960
λ 0.3813.0813.1030.8660.2422.7722.7820.8990.1712.2522.2670.9190.1271.6681.6830.935
q0.2801.3181.3470.9760.2231.1371.2130.9730.1121.0231.1030.9670.0870.5830.6290.959
01−13 μ 0.0231.1021.1020.9190.0221.0591.0590.9320.0051.0131.0130.9380.0040.9130.9140.941
σ 0.1020.3210.3210.9790.0710.2460.2560.9750.0680.2220.2510.9710.0570.1910.2240.961
λ −0.1343.1893.1880.972−0.1163.0013.0020.965−0.0722.7262.7250.961−0.0492.4142.4170.955
q0.3861.8351.8370.9740.3201.5261.5820.9700.2021.4201.4840.9650.1701.2581.2440.961
0123 μ 0.1281.0141.0220.8990.0690.9300.9300.9120.0610.8300.8320.9430.0460.8130.8150.944
σ 0.1340.3240.3250.9840.0790.2740.2780.9770.0520.2460.2660.9620.0490.2090.2270.959
λ 0.1592.9392.9420.8880.1282.7392.7520.9110.1162.6432.6550.9270.0822.2532.2640.935
q0.3221.5561.5560.9810.2621.4111.4750.9720.1741.2711.2370.9610.0541.0271.0390.958
−54−22 μ −0.3544.4844.4950.900−0.2254.4054.4050.917−0.1373.5623.6000.927−0.0402.6412.6700.935
σ 0.3171.5501.5530.9700.2811.3111.3650.9640.2371.0701.1320.9580.1700.7930.8370.957
λ −0.2472.8622.8610.885−0.2232.8082.8170.906−0.1912.2402.2730.919−0.1301.6311.6580.929
q0.3091.2921.3280.9740.2451.1551.2370.9660.2220.8510.9090.9610.1570.5810.6170.955
101613 μ 0.67717.50717.5710.9120.59413.64313.9440.9140.43011.89111.8890.9260.1849.2229.2400.945
σ 1.2374.9484.9510.9040.9574.1594.2050.9270.6513.4543.8270.9370.4802.0302.1130.949
λ 0.1283.2063.2040.8980.1183.0843.0880.9150.1022.7262.7270.9220.0752.3392.3440.938
q0.4701.4011.4010.9810.4531.2731.2400.9790.3331.0871.0520.9610.1420.9270.9380.959
Table 2. Descriptive statistics for kevlar dataset.
Table 2. Descriptive statistics for kevlar dataset.
n X ¯ S b 1 b 2 min ( X ) max ( X )
76 1.959 1.573 1.979 8.161 0.025 9.096
Table 3. Estimated parameters and standard errors (in parentheses), log-likelihood, AIC and BIC values, and KSS with p-values for TESN and STESN models in kevlar dataset.
Table 3. Estimated parameters and standard errors (in parentheses), log-likelihood, AIC and BIC values, and KSS with p-values for TESN and STESN models in kevlar dataset.
EstimationsTESNSTESN
μ 3.374 ( 0.384 ) 0.870 ( 0.790 )
σ 1.668 ( 0.166 ) 1.077 ( 0.247 )
λ 3.939 ( 1.052 ) 15.562 ( 11.223 )
q 2.873 ( 0.896 )
log-likelihood 133.846 122.052
AIC 273.693 250.103
BIC 280.685 261.426
KSS 0.122 0.077
p-value 0.191 0.728
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rivera, P.A.; Gallardo, D.I.; Venegas, O.; Bourguignon, M.; Gómez, H.W. An Extension of the Truncated-Exponential Skew- Normal Distribution. Mathematics 2021, 9, 1894. https://doi.org/10.3390/math9161894

AMA Style

Rivera PA, Gallardo DI, Venegas O, Bourguignon M, Gómez HW. An Extension of the Truncated-Exponential Skew- Normal Distribution. Mathematics. 2021; 9(16):1894. https://doi.org/10.3390/math9161894

Chicago/Turabian Style

Rivera, Pilar A., Diego I. Gallardo, Osvaldo Venegas, Marcelo Bourguignon, and Héctor W. Gómez. 2021. "An Extension of the Truncated-Exponential Skew- Normal Distribution" Mathematics 9, no. 16: 1894. https://doi.org/10.3390/math9161894

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