Next Article in Journal
On Quantum Duality of Group Amenability
Next Article in Special Issue
Generation of Julia and Mandelbrot Sets via Fixed Points
Previous Article in Journal
On Generalized Hardy–Rogers Type α-Admissible Mappings in Cone b-Metric Spaces over Banach Algebras
Previous Article in Special Issue
Parametric Jensen-Shannon Statistical Complexity and Its Applications on Full-Scale Compartment Fire Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Asymmetric Alpha-Power Skew-t Distribution

by
Roger Tovar-Falón
1,*,
Heleno Bolfarine
2 and
Guillermo Martínez-Flórez
1
1
Departamento de Matemáticas y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Montería 230027, Colombia
2
Departamento de Estatística, IME, Universidade de São Paulo, São Paulo 1010, Brazil
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(1), 82; https://doi.org/10.3390/sym12010082
Submission received: 13 November 2019 / Revised: 24 December 2019 / Accepted: 25 December 2019 / Published: 2 January 2020
(This article belongs to the Special Issue Symmetry in Applied Mathematics)

Abstract

:
In this paper, we propose a new asymmetric and heavy-tail model that generalizes both the skew-t and power-t models. Properties of the model are studied in detail. The score functions and the elements of the observed information matrix are given. The process to estimate the parameters in model is discussed by using the maximum likelihood approach. Also, the observed information matrix is shown to be non-singular at the whole parametric space. Two applications to real data sets are reported to demonstrate the usefulness of this new model.

1. Introduction

In recent years, there has been considerable interest in the statistical literature related to flexible families of distributions able of modeling data that present high degree of asymmetry, with kurtosis index greater or smaller than the captured by normal model. In this context, two proposals that have shown a promising behavior in this type of situations are the skew-normal (SN) distribution of Azzalini [1] and the power-normal (PN) distribution of Durrans [2]. The SN distribution has been widely studied by many authors, and its main drawback is that it presents singular Fisher information matrix, implying the inference is useless from the theory of large samples using the maximum likelihood (ML) approach. Although the PN model has a shorter asymmetry range than SN distribution, it presents non-singular information matrix and can easily be extended to censored scenarios, as it has a simple distribution function, see, for example, in Martínez-Flórez et al. [3].
The PN model is part of a wide family of distributions known as alpha-power, which has been widely studied by many authors. In addition to the normal distribution, the Birnbaum–Saunders (BS) distribution [4] has also been considered, see, for example, in Martínez-Flórez et al. [5], who propose an extension of the BS distribution based on the asymmetric alpha-power family of distributions to illustrate the applicability of the new proposal with a data set is related to the lifetimes in cycles × 10 3 n = 101 aluminum 6061 T 6 pieces cut in parallel angle to the rotation direction of rolling at the rate of 18 cycles per second and maximum stress of 21.000 psi. More details of the PN distribution can be found in Gupta and Gupta [6] and Pewsey et al. [7].
An alternative propose for modeling asymmetric data that unifies the two previous approaches was introduced by Martínez-Flórez et al. [8]. The proposed model, which is called alpha-power skew-normal (APSN), has non-singular Fisher information matrix, and it can fit data with much more asymmetry than PN models it can handle. In addition, symmetry can be tested by using the likelihood ratio statistic, as the properties of large samples are satisfied for the ML estimator.
Another set of distributions with non-singular information matrices, useful for modeling asymmetric and heavy-tailed data, are based on generalizations of the Student-t distribution, see, for example, in [9,10,11,12,13]. Azzalini and Capitanio [9] for example, introduced a skew-t (ST) distribution as an extension of the SN model for modeling asymmetric and heavy-tailed data as follows; The random variable X is said to have the ST distribution with parameter λ and degrees of freedom ν , if X has the probability density function (PDF) given by
f S T ( x ; λ , ν ) = 2 f T x ; ν F T λ ν + 1 x 2 + ν x ; ν + 1 , x R
where λ R is a parameter that controls the skewness of the distribution, and f T · ; ν and F T · ; ν denote the PDF and the cumulative distribution function (CDF) of a standard Student-t distribution with ν degree of freedom, respectively. The ST distribution, like an extension of the SN model, inherits the problem of the singularity of the information matrix and before this inconvenience Zhao and Kim [14] proposed the power Student-t (PT) distribution, whose information matrix is non-singular and for a given degree of freedom, the kurtosis range surpasses the kurtosis range of the skew-t model at all times. The PT distribution is defined as follows. The random variable X is said to have the PT distribution with parameter α , and degrees of freedom ν , if X has PDF given by
f P T ( x ; α , ν ) = α f T x ; ν F T x ; ν α 1 , x R
where α > 0 is a parameter that controls the form of the distribution, and, again, f T · ; ν and F T · ; ν denote the PDF and the CDF of a standard Student-t distribution, respectively.
Based on the properties of the ST model, to fit data with high degree of asymmetry and the characteristic of the PN model to capture kurtosis larger than the normal model, in this paper, we introduce a new distribution for modeling asymmetric and heavy-tailed data. The proposed model possess non-singular information matrix, and it is able to fit data with far more asymmetry than ST and PT models can handle and with large sample properties satisfied for the ML estimator. The model introduced in this paper is named as alpha-power skew-t (APST) model and it extends both, ST and PT models. The APSN model by Martínez-Flórez et al. [8] is also a particular case when ν tends to infinite. Note that symmetry can be tested using the likelihood ratio statistics with its large sample chi-square distribution.
The rest of this paper is organized as follows. Section 2 introduces the APST model and some of its properties like moments are studied. In particular, skewness and kurtosis indices are computed showing that their ranges surpass those of the ST and PT models. Section 3 deals with the ML estimation for the location-scale situation and its observed information matrix is derived. The extension to censored data is also presented. Finally, two applications are shown in Section 4, revealing that the model proposed can present much improvement over competitors.

2. The Alpha-Power Skew-t Distribution

Definition 1.
The random variable X is said to have an alpha-power skew-t (APST) distribution, if X has PDF given by
f A P S T ( x ; λ , α , ν ) = α f S T ( x ; λ , ν ) F S T ( x ; λ , ν ) α 1 ,
for x R , λ R , and α , ν R + . Functions f S T ( · ) and F S T ( · ) denote the PDF and the CDF of the standard ST distribution. A random variable having f A P S T ( x ; λ , α , ν ) distribution is denoted shortly by X APST ( λ , α , ν ) .
Figure 1 displays the form of the APST distribution for some selected values of the parameters λ and α for ν = 6 . Note from the figure that the asymmetry and kurtosis of the APST distribution are affected by the parameters α and λ ; therefore, the APST model is more flexible to model data that can be highly skewed, as well as heavier tails than ST and PT models.
The following result provides some special cases of the model (3), which occur for different values of λ , α , and ν .
Proposition 1.
Let X APST ( λ , α , ν ) ,
(i) 
if λ = 0 , then X PT ( α , ν ) ,
(ii) 
if α = 1 , then X ST ( λ , ν ) ,
(iii) 
if λ = 0 and α = 1 , then X T ( ν ) , where T ( ν ) denotes the Student-t disribution with ν degree of freedom.
(iv) 
if ν + , then X APSN ( λ , α ) ,
(v) 
if λ = 0 and ν + , then X PN ( α ) ,
(vi) 
if α = 1 and ν + , then X SN ( λ ) ,
(vii) 
if λ = 0 , α = 1 and ν + , then X N ( 0 , 1 ) ,
Proof. 
The proof of (i)–(vii) is immediate from the definition of APST distribution. □

2.1. Moments

The following proposition presents an expression to compute the k-th moment of a random variable APST ( λ , α , ν ) .
Proposition 2.
Let X APST ( λ , α , ν ) , then
E X k = E F S T 1 ( Y ; λ , ν ) k
where Y follows a Beta ( α , 1 ) distribution and F S T 1 ( · ; λ , ν ) is the inverse of the function F S T ( · ; λ , ν ) .
Proof. 
We have by definition that
E X k = R x k α f S T ( x ) F S T ( x ; λ , ν ) α 1 d x
thus, letting y = F S T ( x ; λ , ν ) , then x = F S T 1 ( y ; λ , ν ) , it follows that
E X k = 0 1 α F S T 1 ( y ; λ , ν ) k y α 1 d y
which is the expected value of the function F S T 1 ( Y ; λ , ν ) k , where Y follows a beta distribution with parameters α and 1. □
The indices of skewness ( β 1 ) and kurtosis ( β 2 ) of APST distribution can be calculated by using the moments (4) as follows,
β 1 = μ 3 3 μ 1 μ 2 + 2 μ 1 3 ( μ 2 μ 1 2 ) 3 / 2 and β 2 = μ 4 4 μ 1 μ 3 + 6 μ 2 μ 1 2 3 μ 1 4 ( μ 2 μ 1 2 ) 2
where μ k = E [ X k ] for k = 1 , , 4 . Table 1 presents the ranges of possible values for the indices of asymmetry and kurtosis for ST ( λ , ν ) , PT ( α , ν ) , and APST ( λ , α , ν ) distributions, for values of λ between −40 and 40, values of α between 0.5 and 50, and for values of ν = 2 , 3 , 4 , 5 , 6 , 7 . It can seen from Table 1 that the length of the admissible intervals for the skewness and the kurtosis parameters of the APST distribution are larger than the corresponding intervals of the ST and PT distributions. This is an indicator that the APST model is more flexible in terms of asymmetry and kurtosis than the ST and PT models.

2.2. Distribution Function

Proposition 3.
Let X APST ( λ , α , ν ) , then the CDF of X, namely, F A P S T ( x ; λ , α , ν ) is
F A P S T ( x ; λ , α , ν ) = F S T ( x ; λ , ν ) α , x R .
Proof. 
The proof is immediate and it follows from results of Durrans [2]. □
The inversion method can be used to generate a random variable with APST distribution. Thus, taking λ R , α , ν R + and a random variable with uniform distribution, namely, U U ( 0 , 1 ) , random variable X with APST ( λ , α , ν ) distribution is generated by taking
X = F S T 1 U 1 / α ; λ , ν .
Remark 1.
We consider a truncated APST ( λ , α ) distribution to obtain a new and useful lifetime distribution. A random variable T has a truncated alpha-power skew-t distribution (at zero), denoted by TAPST ( λ , α , ν ) , if its PDF is given by
f ( t ) = α f S T ( t , λ , ν ) F S T ( t , λ , ν ) α 1 1 F S T ( 0 , λ , ν ) α ; t > 0
The survival and hazard rate functions of a random variable T following a TAPST ( λ , α , ν ) distribution are given by
S T ( t ) = P ( T > t ) = 1 F S T ( 0 , λ , ν ) α F S T ( t , λ , ν ) α 1 F S T ( 0 , λ , ν ) α ; t > 0
and
h T ( t ) = α f S T ( t , λ , ν ) f S T ( t , λ , ν ) α 1 1 F S T ( 0 , λ , ν ) α F S T ( t , λ , ν ) α ; t > 0
respectively.

2.3. Location and Scale Extension

We can also consider a generalization of a APST distribution by adding location and scale parameters. The following definition gives a generalization of the APST model.
Definition 2.
Let X APST ( λ , α , ν ) . The APST density of location and scale is defined as the distribution of Y = μ + σ X , for μ R and σ > 0 . The corresponding PDF is given by
f A P S T ( y ; μ , σ , λ , α , ν ) = α σ f S T y μ σ ; λ , ν F S T y μ σ ; λ , ν α 1 , x R ,
for λ R and α , ν R + . A random variable following a APST distribution of location and scale is denoted by Y APST ( μ , σ , λ , α , ν ) .
The k-th moment of a random variable Y APST ( μ , σ , λ , α , ν ) can be obtained from the formula
E Y k = i = 0 k k i μ i σ k i E X k i ,
where X APST ( λ , α , ν ) .

3. Statistical Inference for APST Distribution

This section concerns likelihood inference about the parameter vector θ = ( μ , σ , λ , α , ν ) of the location-scale family defined in Equation (9). Let Y = ( Y 1 , , Y n ) be a random sample of the distribution APST ( μ , σ , λ , α , ν ) . The log-likelihood function for θ = ( μ , σ , λ , α , ν ) can be written as follows,
( θ ; Y ) n log α n log σ n 2 log ν   + n log Γ ν + 1 2 n log Γ ν 2 ν + 1 2 i = 1 n log 1 + z i 2 ν   + i = 1 n log F T λ z i ν + 1 z i 2 + ν ; ν + 1 + ( α 1 ) i = 1 n log F S T z i ; λ , ν
where z i = ( y i μ ) / σ . Thus, by differentiating the log-likelihood function, we obtain the following score equations,
( θ ; Y ) μ = ν + 1 σ ν i = 1 n z i 1 + z i 2 ν 1 λ σ i = 1 n w i 1 + z i 2 ν 1 f T λ z i w i ; ν + 1 F T λ z i w i ; ν + 1 α 1 σ i = 1 n f S T ( z i ; λ , ν ) F S T ( z i ; λ , ν ) = 0
( θ ; Y ) σ = n σ + ν + 1 σ ν i = 1 n z i 2 1 + z i 2 ν 1 λ σ i = 1 n z i w i 1 + z i 2 ν 1 f T λ z i w i ; ν + 1 F T λ z i w i ; ν + 1 α 1 σ i = 1 n z i f S T ( z i ; λ , ν ) F S T ( z i ; λ , ν ) = 0
( θ ; Y ) λ = i = 1 n z i w i f T λ z i w i ; ν + 1 F T λ z i w i ; ν + 1 α 1 π ( 1 + λ 2 ) i = 1 n 1 + ( 1 + λ 2 ) z i 2 / ν ν 2 F S T ( z i ; λ , ν ) = 0 ,
( θ ; Y ) α = n α + i = 1 n log F S T ( z i ; λ , ν ) = 0 ,
( θ ; Y ) ν = n α 2 ψ ν + 1 2 ψ ν 2 1 ν 1 2 i = 1 n log 1 + z i 2 ν + ν + 1 2 ν 2 i = 1 n z i 2 1 + z i 2 ν 1 + λ 2 ν ( ν + 1 ) i = 1 n z i 3 w i 1 + z i 2 ν 1 f T λ z i w i ; ν + 1 F T λ z i w i ; ν + 1 λ 2 ν ( ν + 1 ) i = 1 n z i w i 1 + z i 2 ν 1 f T λ z i w i ; ν + 1 F T λ z i w i ; ν + 1 ( α 1 ) 2 π ( ν + 1 ) λ ( 1 + λ 2 ) i = 1 n 1 + ( 1 + λ 2 ) z i 2 / ν ν 2 F S T ( z i ; λ , ν ) + α 1 2 i = 1 n g ( z i ; ν ) F S T z i ; λ , ν = 0
where ψ ( · ) is the digamma function, w i = ν + 1 x i 2 + ν for i = 1 , , n , and g ( x ; ν ) is the function defined by
g ( x ; ν ) = x ( ν + 1 ) ν 2 s 2 1 + s 2 ν 1 log 1 + s 2 ν f S T ( s ; λ , ν ) d s λ π ν x s 1 + s 2 ν 1 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s
Equations (11)–(15) include nonlinear functions; therefore, it is not possible to obtain explicit forms of the maximum likelihood estimators (MLEs), and they must be calculated by using numerical methods. In this work, we used the maxLik function of R Development Core Team [15] which uses the Newton–Raphson optimization method. The elements of the observed information matrix are easily obtained after calculating the second derivative of the log-likelihood function and multiplying by −1, that is,
j θ i θ k = ( θ ; Y ) θ i θ k , i , k = 1 , 2 , , 5
where θ = ( μ , σ , λ , α , ν ) . This elements are given in the Appendix A. To find the standard errors (EE) of the MLEs and calculate confidence intervals, the information matrix I (or Fisher information) must be calculated, which is defined as the expected value of the second derived from the log-likelihood function or less the expected value of the Hessian matrix; from this matrix, we calculate the EE as the diagonal elements of the inverse of this matrix. The elements of the I matrix are obtained as
I ( i , k ) = E ( θ ; Y ) θ i θ k , i , k = 1 , 2 , , 5
The role of the Fisher information in the asymptotic theory of maximum-likelihood estimation was emphasized by Ronald Fisher following some initial results by Francis Edgeworth, see Lehman and Casella [16] and Frieden [17] for more details. The Fisher-information matrix is used to calculate the covariance matrices associated with maximum-likelihood estimates, and it can also be used in the formulation of test statistics, such as the Wald test.
As the expected value under the APST distribution and the second-order derivatives are not direct, numerical methods must be used to obtain the explicit form of the information matrix I. Therefore, we use the observed information matrix to calculate the standard errors in the rest of the document.
When ν tends to infinite the ST distribution converges to the SN distribution and we recall that the information matrix of a random variable X SN ( μ , σ , λ ) which is denoted by I λ ( φ ) , where φ = ( μ , σ , λ ) , is singular for λ = 0 . Therefore, it is convenient to use a centered parameterization of the ST distribution proposed by Arellano-Valle and Azzalini [18].
The centered parameterization of the SN distribution was proposed as an alternative to the problem of singularity of the information matrix of the SN when λ = 0 . Arellano-Valle and Azzalini [19] proposed a second representation of the SN by defining a new random variable X as
X = μ + σ Z E [ Z ] Var [ Z ] ,
where μ R and σ > 0 are parameters of the random variable X and Z SN ( λ ) . This representation is called centered parameterization, as E [ X ] = μ and Var [ X ] = σ 2 and it is denoted by CSN ( μ , σ , γ 1 ) , where 0 . 9953 < γ 1 < 0 . 9953 . Under the centered parameterization model, μ , σ , and γ 1 = β 1 represent the mean, the standard deviation and the skewness index of X, respectively. If Z SN ( λ ) then E [ Z ] = b δ and Var [ Z ] = 1 ( b δ ) 2 , where b = 2 / π and δ = λ / 1 + λ 2 ; it has that the random variable X can be written as X = μ + σ Z which has SN ( λ 1 , λ 2 , λ ) distribution, where
λ 1 = μ c σ γ 1 1 / 3 , λ 2 = σ 1 + c 2 γ 1 2 / 3 , λ = c γ 1 1 / 3 b 2 + c 2 ( b 2 1 ) γ 1 2 / 3
with c = { 2 / ( 4 π ) } 1 / 3 . Under this denomination, the information matrix can be written as I γ 1 = D I λ D , where D is a matrix that represents the derivative of the parameters of the standard representation ( λ 1 , λ 2 and λ ) regarding to the new parameters ( μ , σ and γ 1 ). It also follows that the information matrix converges to a diagonal matrix Σ c 1 = diag ( σ 2 , σ 2 / 2 , 6 ) when λ 0 . This guarantees the existence and uniqueness of the MLEs of λ 1 and λ 2 for each fixed value of λ .
Following this same line of thought, we suppose that Y follows the model (1) with location parameter μ R and scale parameter σ > 0 , that is,
f S T ( y ; μ , σ , λ , ν ) = 2 σ f T y μ σ ; ν F T λ ν + 1 Q y + ν y μ σ ; ν + 1 , y R
where λ R and Q y = ( ( y μ ) / σ ) 2 . This representation relates to the direct parameterization of the ST distribution with parameter vector ρ = ( μ , σ , λ , ν ) . It follows that Z T = ( Y μ ) / σ ST ( λ , ν ) , and by the stochastic representation of the ST distribution is given by Z T = Z / V , where Z SN ( λ ) and V χ ν 2 / ν . This entails to compute the first four cumulants of Z T denoted by μ 1 ( δ , ν ) , μ 2 ( δ , ν ) , μ 3 ( δ , ν ) and μ 4 ( δ , ν ) , see [18]. The centered parameterization of the ST distribution of a random variable Y comes by defining
μ t = E [ Y ] = μ + σ μ 1 ( δ , ν ) = μ + σ b ν δ σ t 2 = Var [ Y ] = σ 2 μ 2 ( δ , ν ) = η 2 ν ν 2 b ν 2 δ 2 ,
γ 1 t = μ 3 ( δ , ν ) μ 2 ( δ , ν ) 3 / 2 = b ν δ μ 2 ( δ , ν ) 3 / 2 ν ( 3 δ 2 ) ν 3 3 ν ν 2 + 2 b ν 2 δ 2 γ 2 t = μ 4 ( δ , ν ) μ 2 ( δ , ν ) 2 = 1 μ 2 ( δ , ν ) 2 3 ν 2 ( ν 2 ) ( ν 4 ) 4 b ν 2 δ 2 ν ( 3 δ 2 ) + 6 b ν 2 δ 2 ν ν 2 4 b ν 4 δ 4 ν 3 3 .
The new representation is defined as the centered skew-t distribution with parameter vector ρ ˜ = ( μ , σ 2 , γ 1 , γ 2 ) . According to Arellano-Valle and Azzalini [18], the information matrix of this representation can be written as
I ( ρ ˜ ) = B I ( ρ ) B ,
where B is a matrix representing the derivative of the parameter vector ρ with respect to the new vector ρ ˜ . It can shown that b ν b when ν , see [18]. Therefore, the parameters of the centered ST model converge to μ t μ , σ t 2 σ 2 , and γ 1 t γ 1 when ν , that is, the parameters of the CSN. As Z T SN ( λ ) when ν , it follows that the random variable Y converges to a distribution with information matrix
I ( μ , σ 2 , γ 1 , α ) = I θ 1 θ 1 I θ 1 , α I θ 1 , α I α , α ,
where the elements of the diagonal correspond to the information of the parameter vector θ 1 = ( μ , σ 2 , γ 1 ) and α , and I θ 1 , α is the joint information of θ 1 = ( μ , σ 2 , γ 1 ) and α . Now, when λ 0 and α = 1 , it can be shown that I θ 1 θ 1 diag ( σ 2 , σ 2 / 2 , 6 ) , with determinant equal to 0 . 3333 / σ 4 , and I θ 1 , α = ( 0 . 9031 / σ , 0 . 5956 / σ , 0 . 7206 ) ; therefore, the determinant | I ( μ , σ 2 , γ 1 , α ) | 0 , and it concludes that the random variable Y converges to a distribution with information matrix non-singular when ν tends to infinite.

3.1. Extension to Censored Data

Based on the goodness of the APST distribution to fit asymmetric and heavy-tailed data, in this section we introduce the censored APST model which we will be denote by CAPST.
Definition 3.
Suppose that the random variable Y follows APST distribution, and consider a random sample Y = ( Y 1 , Y 2 , , Y n ) where only the Y i values greater than a constant k are recorded. In addition, for values Y i k only the value of k is recorded. Therefore, for i = 1 , 2 , , n , the observed values Y i o can be written as
Y i o = k , i f Y i k , Y i , i f Y i > k .
The resulting sample is said to be a censored APST, and we say that Y is a censored random variable APST. We will use the notation Y CAPST ( θ ) , where θ = ( μ , σ , λ , α , ν ) .
From Definition 3 it follows that P ( Y i o = k ) = P ( Y i k ) = F S T ( k μ ) / σ α and for the observations Y i o = Y i , the distribution of Y i o is the same of Y i , i.e., Y i o APST ( θ ) . For convenience, we choose to work with the case of left-censored data; however, the followings results can be extended to other types of censorship.

3.2. Properties of the CAPST Model

Let Y CAPST ( μ , σ , λ , α , ν ) ,
  • If α = 1 , then Y CST ( μ , σ , λ , ν ) , where CST indicates the censored skew-t model.
  • If λ = 0 , then Y CPT ( μ , σ , α , ν ) , where CPT indicates the censored power-t model.
  • If α = 1 and λ = 0 , then Y CT ( μ , σ , ν ) , that is, the censored Student-t model follows.
  • If ν + , then Y CAPSN ( μ , σ , λ , α ) , where CAPSN indicates the censored alpha-power skew-normal model.
  • If α = 1 and ν + , then Y CSN ( μ , σ , λ ) , that is, the censored skew-normal model follows.
  • If λ = 0 and ν + , then Y CPN ( μ , σ , α ) , that is, the censored power-normal model follows.
  • If α = 1 , λ = 0 and ν + , then Y CN ( μ , σ 2 ) , that is, the censored normal model follows.
The estimates of the parameters of the model can be obtained via maximum likelihood method, where the log-likelihood function is given by
( θ ; Y ) α 0 log F S T k μ σ ; λ , ν + n 1 log α n 1 log σ n 1 2 log ν + n 1 log Γ ν + 1 2 n 1 log Γ ν 2 ν + 1 2 1 log 1 + x i 2 ν + 1 log F T λ x i ν + 1 x i 2 + ν ; ν + 1 + ( α 1 ) 1 log F S T x i ; λ , ν
where x i = ( y i μ ) / σ ; 1 and 0 are the sum over censored individuals and uncensored individuals, respectively; and n 1 is the number of uncensored individuals.

4. Real Data Applications

In this section, we illustrate the applicability of the proposed model in Section 2 by analyzing two data sets. We use the statistical software R [15], version 3.5.3 with the package maxLike for maximizing the corresponding likelihood functions. For comparing purposes of various models, the AIC Akaike [20], BIC Schwarz [21], and corrected AIC (CAIC) Bozdogan [22] information criteria were used.

4.1. Application 1: Volcano Heights Data

Consider the data set related to heights of 1520 volcanoes in the world which is available in website dx.doi.org/10.5479/si.GVP.VOTW4-2013 [23]. Table 2 presents the summary statistics for the data set. It can be noted that the asymmetry and kurtosis indices seem to indicate that the use of an asymmetric and heavy-tailed model is appropriate to analyze this data set. We analyzed these data by fitting the Student-t, ST, PT, and APST distributions.
Table 3 shows the parameter estimates, together with their corresponding standard errors (SE). Note that the values of the standard errors of the μ and σ estimates for the APST model are smaller than the corresponding standard errors of the respective parameters for the Student-t, ST, and PT models. Table 3 also presents some model selection criteria, together with the values of the log-likelihood. The AIC, BIC, and CAIC criteria indicate that the APST model seems to provide better fit to the volcanoes heights data than the T, ST, and PT models, supporting the asymmetry assertion of the volcano’s heights variable. Figure 2 shows the graphs QQplot of the fitted models. It can be clearly seen from the figure that the APST model fits the data better than the Student-t, ST, and PT models. In addition, we can use the likelihood ratio (LR) test statistic to conform our claim. To do this, we consider the following hypotheses,
H 0 : ( λ , α ) = ( 0 , 1 ) ( T ( μ , σ , ν ) ) v . s H 1 : ( λ , α ) ( 0 , 1 ) ( APST ( μ , σ , λ , α , ν ) ) ,
The value of the LR test statistic is 2 log ( Λ ) = 2 T ( θ ^ ) A P S T ( θ ^ ) = 134.823 and comparing this quantity with χ 2 2 = 5.9914, the null hypotheses is rejected. The APST model is also compared with the ST and PT models by considering the hypotheses
H 01 : α = 1 ( ST ( μ , σ , λ , ν ) ) v . s H 11 : α 1 ( APST ( μ , σ , λ , α , ν ) ) ,
and
H 02 : λ = 0 ( PT ( μ , σ , α , ν ) ) v . s H 12 : λ 0 ( APST ( μ , σ , λ , α , ν ) ) ,
respectively. The respective values of the LR test statistic are given by 2 log ( Λ 1 ) = 2 S T ( θ ^ ) A P S T ( θ ^ ) = 26.620 and 2 log ( Λ 2 ) = 2 P T ( θ ^ ) A P S T ( θ ^ ) = 45.660 and comparing these quantities with χ 1 2 = 3.8414, both null hypotheses are rejected. Finally, Figure 3left shows the histogram of the volcano heights variable, whereas Figure 3right presents the empirical CDF (solid line) together with the CDF of the fitted APST model (dotted line).

4.2. Application 2: Stellar Abundances Data

The second data set is related to measurements for 68 solar-type stars, which are available in the package astrodatR of the software R [24] under the name Stellar abundances. These data were previously analyzed Mattos et al. [25] by using the Scale Mixture of Skew Normal Censored Regression (SMSNCR) models. We take only the response variable: log N(Be), which represents the log of the abundance of beryllium scaled to Sun’s abundance (i.e., the Sun has log N ( B e ) = 0 . 0 )
In astronomical research, a previously identified sample of objects (stars, galaxies, quasars, X-ray sources, etc.) is observed at some new wavebands. According to Feigelson [24], due to limited sensitivities, some objects may be undetected, leading to upper limits in their derived luminosities. For this dataset we have 12 left-censored data points, i.e., 12 undetected beryllium measurement, that represents 19.35% of observations. Table 4 presents the ML estimates for the parameters of the censored Studen-t (CT), censored skew-t (CST), censored power-t (CPT), and censored alpha-power skew-t (CAPST) models, together with their corresponding standard errors. Table 4 also compares the fit of the four models using the model selection criteria (AIC, CAIC and BIC). Note that, again, the CAPST model with heavy tails have better fit than the CT, CST, and CPT models.
To identify atypical observations and/or model mispecification, we analyzed the transformation of the martingale residual, r M T i , proposed in Barros et al. [26]. These residuals are defined by
r M T i = sign ( r M i ) 2 [ r M i + δ i log ( δ i r M i ) ] , i = 1 , , n
where r M i = δ i + log S ( y i ; θ ^ ) is the martingal residual proposed by Ortega et al. [27], where δ i = 0 , 1 indicates whether the i-th observation is censored or not, respectively; sign ( r M i ) denotes the sign of r M i ; and S ( y i ; θ ^ ) = P θ ^ ( Y i > y i ) represents the survival function evaluated at y i , where θ ^ are the MLE for θ . The plots of r M T i with generated confidence envelopes are presented in Figure 4. From this figure, we can see clearly that the CST, CPT, and CAPST models fit better the data than the CT model, since, in that cases, there are not observations which lie outside the envelopes. The Figure 5 shows the graph of the densities of the different models fitted to the stellar abundances data. From the figure, the CAPST model seems to fit better the stellar abundances data than CT, CST and CPT models.

5. Conclusions

In this work, a new asymmetric model has been introduced. It is based on the combination of skew-t [1] and power-t [2] models. The new model presents greater ranges of asymmetry and kurtosis, which is very useful for modeling skewed and heavy-tailed data. The problem of estimating the parameters in the model is dealt by using the maximum likelihood approach which is also used for developing large sample properties for the estimators. The elements of the observed information matrix are analytically obtained. The likelihood ratio statistics can be used for testing the APST null hypothesis since the Student-t, ST, and PT models are special cases of the model entertained. Two applications to volcano heights data and stellar abundances data indicate that the proposed model can be a useful alternative to the ST and PT models.

Author Contributions

Individual contributions to this article: conceptualization, R.T.-F., H.B., and G.M.-F.; methodology, R.T.-F., H.B., and G.M.-F.; software, R.T.-F., H.B., and G.M.-F.; validation, R.T.-F., H.B., and G.M.-F.; formal analysis, R.T.-F., H.B., and G.M.-F.; investigation, R.T.-F., H.B., and G.M.-F.; resources, R.T.-F., H.B., and G.M.-F.; writing-original draft preparation, R.T.-F., H.B., and G.M.-F.; writing-review and editing, R.T.-F., H.B., and G.M.-F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The research of R. Tovar-Falón was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq/Brazil) under grant 140831/2014-2.

Acknowledgments

We thank the anonymous referees for helpful suggestions which improved the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this section, expressions for the elements of the observed information matrix of the alpha-power skew-t model are provided. Initially we suppose that Y APST ( μ , σ , λ , α , ν ) , and for i = 1 , , n we define z i = ( y i μ ) / σ , w i = ( ν + 1 ) / ( z i 2 + ν ) , r 1 ( z ; ν ) = f T ( z ; ν ) / F T ( z ; ν ) , r 2 ( z ; λ , ν ) = f S T ( z ; λ , ν ) / F S T ( z ; λ , ν ) , and r 3 ( z ; λ , ν ) = 1 + ( 1 + λ 2 ) z 2 ν ν 2 / F S T ( z ; λ , ν ) . Denoting the elements of the observed information matrix of the APST model by j μ μ , j μ σ , , j α α , and after some algebraic manipulations, we obtain
j μ μ = 1 σ 2 ν + 1 ν 2 i = 1 n z i 2 1 + z i 2 ν 2 + 1 σ 2 ν + 1 ν i = 1 n 1 + z i 2 ν 2 + λ σ 2 i = 1 n z i w i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + 2 σ 2 λ ν i = 1 n z i w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 3 σ 2 ν + 2 ν i = 1 n z i w i 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ 2 ν + 1 ν i = 1 n 1 + z i 2 ν 3 r 1 ( λ z i w i ; ν + 1 ) 2
λ σ 2 α 1 π i = 1 n 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 σ 2 ν + 1 ν i = 1 n z i 1 + z i 2 ν 1 r 2 ( z i ; λ , ν ) + α 1 σ 2 i = 1 n r 2 ( z i ; λ , ν ) 2 j μ σ = 2 σ 2 ν + 1 ν i = 1 n z i 1 + z i 2 ν 2 + λ 3 σ 2 ν + 2 ν i = 1 n z i 2 w i 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + 2 σ 2 λ ν i = 1 n z i 2 w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) λ σ 2 i = 1 n w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ 2 ν + 1 ν i = 1 n z i 1 + z i 2 ν 3 r 1 ( λ z i w i ; ν + 1 ) 2 λ σ 2 α 1 π i = 1 n z i 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) α 1 σ 2 i = 1 n r 2 ( z i ; λ , ν ) + α 1 σ 2 ν + 1 ν i = 1 n z i 2 1 + z i 2 ν 1 r 2 ( z i ; λ , ν ) + α 1 σ 2 i = 1 n z i r 2 ( z i ; λ , ν ) 2
j μ λ = λ 2 σ ν + 2 ν i = 1 n z i 2 w i 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + 1 σ i = 1 n w i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) λ σ ν + 1 ν i = 1 n z i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) 2 + α 1 π σ i = 1 n z i 2 1 + z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 π σ 1 1 + λ 2 i = 1 n r 2 ( z i ; λ , ν ) r 3 ( z i ; λ , ν ) j μ α = 1 σ i = 1 n r 2 ( z i ; λ , ν ) j μ ν = 1 σ ν i = 1 n z i 1 + z i 2 ν 1 + ν + 1 σ ν 2 i = 1 n z i 1 + z i 2 ν 2 + λ σ ν 2 i = 1 n z i 2 w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ 1 ν ( ν + 1 ) i = 1 n ( z i 2 1 ) w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 3 2 σ ν + 2 ν 2 ( ν + 1 ) i = 1 n z i 2 ( z i 2 1 ) w i 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) λ 2 2 σ 1 ν 2 i = 1 n z i ( z i 2 1 ) 1 + z i 2 ν 3 r 1 ( λ z i w i ; ν + 1 ) 2 + λ 2 π σ α 1 ν + 1 i = 1 n z i 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 2 π σ ( ν + 1 ) λ 1 + λ 2 i = 1 n r 2 ( z i ; λ , ν ) r 3 ( z i ; λ , ν ) + α 1 2 σ ν + 1 ν 2 i = 1 n z i 2 1 + z i 2 ν 1 r 2 ( z i ; λ , ν ) α 1 2 σ i = 1 n log 1 + z i 2 ν r 2 ( z i ; λ , ν ) α 1 2 σ i = 1 n g ( z i , ν ) F S T ( z i ; λ , ν ) r 2 ( z i ; λ , ν ) + λ 2 π σ α 1 ν i = 1 n z i 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν )
j σ σ = n σ 2 + 1 σ 2 ν + 1 ν i = 1 n z i 2 1 + z i 2 ν 1 + 2 σ 2 ν + 1 ν i = 1 n z i 2 1 + z i 2 ν 2 2 λ σ 2 i = 1 n z i w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) λ σ 2 ν i = 1 n z i 3 w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 3 σ 2 ν + 2 ν i = 1 n z i 3 w i 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ 2 ν + 1 ν i = 1 n z i 2 1 + z i 2 ν 3 r 1 ( λ z i w i ; ν + 1 ) λ σ 2 α 1 π i = 1 n z i 2 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) 2 ( α 1 ) σ 2 i = 1 n z i r 2 ( z i ; λ , ν ) α 1 σ 2 ν + 1 ν i = 1 n z i 3 1 + z i 2 ν 1 r 2 ( z i ; λ , ν ) + α 1 σ 2 i = 1 n z i 2 r 2 ( z i ; λ , ν )
j σ λ = 1 σ i = 1 n z i w i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ ν + 2 ν i = 1 n z i 3 w i 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ σ ν + 1 ν i = 1 n z i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) 2 α 1 π σ i = 1 n z i 2 1 + z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 π σ 1 1 + λ 2 i = 1 n r 2 ( z i ; λ , ν ) r 3 ( z i ; λ , ν ) j σ α = 1 σ i = 1 n z i r 2 ( z i ; λ , ν )
j σ ν = 1 σ ν i = 1 n z i 2 1 + z i 2 ν 1 + 1 σ ν + 1 ν 2 i = 1 n z i 2 1 + z i 2 ν 2 + λ σ ν 2 i = 1 n z i 3 w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) + λ 2 σ 1 ν ( ν + 1 ) i = 1 n z i ( z i 2 1 ) w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) λ 3 2 σ ν + 2 ν 2 ( ν + 1 ) i = 1 n z i 3 ( z i 2 1 ) w i 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) λ 2 2 σ ν 2 i = 1 n z i 2 ( z i 2 1 ) 1 + z i 2 ν 3 r 1 ( λ z i w i ; ν + 1 ) 2 + λ 2 π σ α 1 ν ( ν + 1 ) i = 1 n z i 2 ( z i 2 1 ) 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 2 π σ ( ν + 1 ) λ 1 + λ 2 i = 1 n z i r 2 ( z i ; λ , ν ) r 3 ( z i ; λ , ν ) + α 1 2 σ ν + 1 ν 2 i = 1 n z i 3 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) α 1 2 σ i = 1 n z i log 1 + z i 2 ν r 1 ( λ z i w i ; ν + 1 ) α 1 2 σ i = 1 n z i g ( z i , ν ) F S T ( z i , λ , ν ) r 2 ( z i ; λ , ν ) j λ λ = λ ( ν + 2 ) ν i = 1 n z i 3 w i 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + ν + 1 ν i = 1 n z i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) 2 2 ( α 1 ) π λ ( 1 + λ 2 ) 2 i = 1 n 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) α 1 π λ 1 + λ 2 ν + 2 ν i = 1 n z i 2 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) + α 1 π 2 ( 1 + λ 2 ) 2 i = 1 n r 3 ( z i ; λ , ν ) 2
j λ α = 1 π ( 1 + λ 2 ) i = 1 n r 3 ( z i ; λ , ν ) j λ ν = 1 2 ν ( ν + 1 ) i = 1 n z i ( z i 2 1 ) w i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 2 2 ν 2 ν + 2 ( ν + 1 ) 2 i = 1 n z i 3 ( z i 2 1 ) w i 1 + z i 2 ν 1 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 2 ν 2 i = 1 n z i 2 ( z i 2 1 ) 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) 2 + α 1 2 π ( ν + 1 ) 1 λ 2 ( 1 + λ 2 ) 2 i = 1 n r 3 ( z i ; λ , ν ) + α 1 2 π ( ν + 1 ) λ ( 1 + λ 2 ) 2 i = 1 n [ r 3 ( z i ; λ , ν ) ] 2 α 1 2 π ( ν + 1 ) λ 2 1 + λ 2 i = 1 n z i 2 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) α 1 2 π ( 1 + λ 2 ) i = 1 n g ( z i , ν ) F S T ( z i , λ , ν ) r 3 ( z i ; λ , ν ) α 1 2 π i = 1 n g 1 ( z i , ν ) F S T ( z i , λ , ν ) j α α = n α 2 j α ν = n 2 ψ ν + 1 2 + n 2 ψ ν 2 + n 2 ν + 1 2 π ( ν + 1 ) λ 1 + λ 2 i = 1 n r 3 ( z i ; λ , ν ) 1 2 i = 1 g ( z i , ν ) F S T ( z i , λ , ν )
j ν ν = n α 2 ν 2 n α 4 ψ 1 ν + 1 2 + n α 4 ψ 1 ν 2 ν 1 2 ν 3 i = 1 n z i 2 1 + z i 2 ν 1 + ν + 1 2 ν 3 i = 1 n z i 2 1 + z i 2 ν 2 + λ 4 ( ν + 1 ) 2 1 ν 2 i = 1 n z i ( z i 2 1 ) ( z i 2 + 4 ν + 3 ) w i 1 + z i 2 ν 2 r 1 ( λ z i w i ; ν + 1 ) λ 4 ν ( ν + 1 ) ψ ν + 2 2 ψ ν + 1 2 1 ν + 1 i = 1 n z i ( z i 2 1 ) w i 1 + z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) λ 3 4 ( ν + 1 ) ν + 2 ν 3 i = 1 n z i 3 ( z i 2 1 ) 1 + z i 2 ν 2 1 + 1 + λ 2 z i 2 ν 1 r 1 ( λ z i w i ; ν + 1 ) + λ 4 ν ( ν + 1 ) i = 1 n z i ( z i 2 1 ) log 1 + λ 2 z i 2 ν + z i 2 r 1 ( λ z i w i ; ν + 1 ) + λ 2 4 ν 3 ( ν + 1 ) i = 1 n z i 2 ( z i 2 1 ) 2 1 + z i 2 ν 3 [ r 1 ( λ z i w i ; ν + 1 ) ] 2 α 1 2 π ( ν + 1 ) 2 λ 1 + λ 2 i = 1 n r 3 ( z i ; λ , ν ) + α 1 4 π ( ν + 1 ) λ ν i = 1 n z i 2 1 + 1 + λ 2 z i 2 ν 1 r 3 ( z i ; λ , ν ) α 1 4 π ( ν + 1 ) λ 1 + λ 2 i = 1 n log 1 + 1 + λ 2 z i 2 ν r 3 ( z i ; λ , ν ) + α 1 4 π ( ν + 1 ) λ 1 + λ 2 ψ ν + 1 2 ψ ν 2 1 ν i = 1 n r 3 ( z i ; λ , ν ) α 1 2 π ( ν + 1 ) λ 1 + λ 2 i = 1 n g ( z i , ν ) F S T ( z i , λ , ν ) r 3 ( z i ; λ , ν ) + α 1 4 π 2 ( ν + 1 ) 2 λ 2 ( 1 + λ 2 ) 2 i = 1 n r 3 ( z i ; λ , ν ) 2 + α 1 4 i = 1 n g ( z i , ν ) F S T ( z i , λ , ν ) 2 α 1 2 i = 1 n g 2 ( z i , ν ) F S T ( z i , λ , ν )
where g ( z ; ν ) is given in Equation (16), and g 1 ( z ; ν ) and g 2 ( z ; ν ) are given in Equations (A1) and (A3), respectively.
g 1 ( x ; ν ) = x ( ν + 1 ) s 2 ν ( s 2 + ν ) log 1 + s 2 ν 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 s d s x s s 2 + ν 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s + x λ 2 ( ν + 2 ) s 3 ( s 2 + ν ) ( ν + ( 1 + λ 2 ) s 2 ) 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s
g 2 ( x ; ν ) = x s 2 ( s 2 ν 2 ν s 2 ) ν 2 ( s 2 + ν ) 2 + 1 2 ( ν + 1 ) s 2 ν ( s 2 + ν ) log 1 + s 2 ν 2 f S T ( s ; λ , ν ) d s + λ 2 π ( ν + 1 ) x s ( s 2 1 ) ( s 2 + ν ) ( ν + 1 ) s 2 ν ( s 2 + ν ) log 1 + s 2 ν × 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s + λ π x s ( s + ν ) 2 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s
+ λ 2 π ψ ν + 1 2 ψ ν 2 1 ν x s s 2 + ν 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s λ 2 π x s s 2 + ν ( ν + 2 ) ( 1 + λ 2 ) s 2 ν ( ν + ( 1 + λ 2 ) s 2 ) log 1 + ( 1 + λ 2 ) s 2 ν × 1 + ( 1 + λ 2 ) s 2 ν ν + 2 2 d s

References

  1. Azzalini, A. A class of distributions which includes the normal ones. Scand. J. Stat. 1985, 12, 171–178. [Google Scholar]
  2. Durrans, S.R. Distributions of fractional order statistics in hydrology. Water Resour. Res. 1992, 28, 1649–1655. [Google Scholar] [CrossRef]
  3. Martínez-Flórez, G.; Bolfarine, H.; Gómez, H.W. The alpha–power tobit model. Commun. Stat. Theory Methods 2013, 42, 633–643. [Google Scholar] [CrossRef]
  4. Birnbaum, Z.W.; Saunders, S.C. A new family of life distributions. J. Appl. Probab. 1969, 6, 319–327. [Google Scholar] [CrossRef]
  5. Martínez-Flórez, G.; Bolfarine, H.; Gómez, H.W. An alpha-power extension for the Birnbaum-Saunders distribution. Statistics 2014, 48, 896–912. [Google Scholar] [CrossRef]
  6. Gupta, R.D.; Gupta, R.C. Analyzing skewed data by power-normal model. Test 2008, 17, 197–210. [Google Scholar] [CrossRef]
  7. Pewsey, A.; Gómez, H.W.; Bolfarine, H. Likelihood–based inference for power distributions. Test 2012, 21, 775–789. [Google Scholar] [CrossRef]
  8. Martínez-Flórez, G.; Bolfarine, H.; Gómez, H.W. Skew-normal alpha-power model. Stat. J. Theor. Appl. Stat. 2014, 48, 1414–1428. [Google Scholar] [CrossRef]
  9. Azzalini, A.; Capitanio, A. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew-t distribution. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2003, 65, 367–389. [Google Scholar] [CrossRef]
  10. Branco, M.D.; Dey, D.K. General class of multivariate skew-elliptical distributions. J. Multivar. Anal. 2001, 79, 99–113. [Google Scholar] [CrossRef] [Green Version]
  11. Durrans, S.R. Multivariate skew t-distribution. Stat. J. Theor. Appl. Stat. 2003, 37, 359–363. [Google Scholar]
  12. Sahu, S.K.; Dey, D.K.; Branco, M.D. A new class of multivariate skew distributions with applications to Bayesian regression models. Can. J. Stat. 2003, 31, 129–150. [Google Scholar] [CrossRef] [Green Version]
  13. Jones, M.C.; Faddy, M.J. A skew extension of the t-distribution, with Applications. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2003, 65, 159–174. [Google Scholar] [CrossRef]
  14. Zhao, J.; Kim, H.M. Power t distributions. Commun. Stat. Appl. Methods 2016, 23, 321–334. [Google Scholar]
  15. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: http://www.R-project.org (accessed on 10 October 2019).
  16. Lehman, E.L.; Casella, G. Theory of Point Estimation, 2nd ed.; Springer: New York, NY, USA, 1998. [Google Scholar]
  17. Frieden, B.R. Science from Fisher Information: A Unification; Cambridge Univerisity Press: Cambridge, UK, 2004. [Google Scholar]
  18. Arellano-Valle, R.B.; Azzalini, A. The centered parameterization and related quantities of the skew–t distribution. J. Multivar. Anal. 2013, 113, 73–90. [Google Scholar] [CrossRef]
  19. Arellano-Valle, R.B.; Azzalini, A. The centered parametrization for the multivariate skew-normal distribution. J. Multivar. Anal. 2008, 99, 1362–1382. [Google Scholar] [CrossRef]
  20. Akaike, H. A new look at statistical model identification. IEEE Trans. Autom. Contr. 1974, 19, 716–722. [Google Scholar] [CrossRef]
  21. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  22. Bozdogan, H. Model selection and akaike’s information criterion (AIC): The general theory and its analytical extensions. Psychometrika 2010, 52, 345–370. [Google Scholar] [CrossRef]
  23. Siebert, L.; Simkin, T.; Kimberly, P. Global Volcanism Program. In Volcanoes of the World; v. 4.6.0.; Venzke, E., Ed.; Smithsonian Institution: Washington, DC, USA, 2013; Available online: https://doi.org/10.5479/si.GVP.VOTW4-2013 (accessed on 10 October 2019).
  24. Feigelson, E.D. astrodatR: Astronomical Data. R Package v. 0.1. Available online: http://CRAN.R-project.org/package=astrodatR (accessed on 10 October 2019).
  25. Mattos, T.; Garay, A.M.; Lachos, V.H. Likelihood-based inference for censored linear regression models with scale mixtures of skew-normal distributions. J. Appl. Stat. 2018, 45, 2019–2066. [Google Scholar] [CrossRef]
  26. Barros, M.; Galea, M.; González, M.; Leiva, V. Influence diagnostics in the tobit censored response model. Stat. Methods Appl. 2010, 19, 379–397. [Google Scholar] [CrossRef]
  27. Ortega, E.M.; Bolfarine, H.; Paula, G.A. Influence diagnostics in generalized log-gamma regression models. Comput. Stat. Data Anal. 2003, 42, 165–186. [Google Scholar] [CrossRef]
Figure 1. Probability density function of APST ( λ , α , 10 ) for some values of λ and α .
Figure 1. Probability density function of APST ( λ , α , 10 ) for some values of λ and α .
Symmetry 12 00082 g001
Figure 2. Volcano height data: QQplot for Student-t, ST, PT, and APST models.
Figure 2. Volcano height data: QQplot for Student-t, ST, PT, and APST models.
Symmetry 12 00082 g002
Figure 3. (Left) Graph of fitted densities to volcano height data. (Right) Empirical CDF and CDF of fitted APST model.
Figure 3. (Left) Graph of fitted densities to volcano height data. (Right) Empirical CDF and CDF of fitted APST model.
Symmetry 12 00082 g003
Figure 4. Stellar abundances data. Envelopes of transformed martingale residuals for CT, CST, CPT, and CAPST models.
Figure 4. Stellar abundances data. Envelopes of transformed martingale residuals for CT, CST, CPT, and CAPST models.
Symmetry 12 00082 g004
Figure 5. Graph of fitted densities to stellar abundances data.
Figure 5. Graph of fitted densities to stellar abundances data.
Symmetry 12 00082 g005
Table 1. Skewness and kurtosis for the models ST ( λ , ν ) , PT ( α , ν ) , and APST ( λ , α , ν ) , for λ ( 40 , 40 ) , α ( 0 . 5 , 50 ) and ν = 2 , 7 .
Table 1. Skewness and kurtosis for the models ST ( λ , ν ) , PT ( α , ν ) , and APST ( λ , α , ν ) , for λ ( 40 , 40 ) , α ( 0 . 5 , 50 ) and ν = 2 , 7 .
Skew t Power t Alpha—Power Skew t
ν SkewnessKurtosisSkewnessKurtosisSkewnessKurtosis
2 ( 0.963 , 0.963 ) ( 3.170 , 3.489 ) ( 0.119 , 3.040 ) ( 1.552 , 10.436 ) ( 2.452 , 14.314 ) ( 1.395 , 864.385 )
3 ( 0.950 , 0.950 ) ( 3.146 , 3.357 ) ( 0.086 , 1.362 ) ( 1.325 , 3.223 ) ( 2.130 , 4.902 ) ( 1.628 , 114.098 )
4 ( 1.853 , 1.853 ) ( 5.099 , 7.824 ) ( 0.530 , 1.178 ) ( 3.461 , 5.299 ) ( 1.898 , 3.215 ) ( 3.153 , 29.874 )
5 ( 0.947 , 0.947 ) ( 3.051 , 3.327 ) ( 0.475 , 0.271 ) ( 1.176 , 3.130 ) ( 1.968 , 3.046 ) ( 3.862 , 19.925 )
6 ( 1.681 , 1.681 ) ( 4.554 , 7.279 ) ( 0.533 , 1.118 ) ( 3.974 , 5.173 ) ( 1.681 , 2.145 ) ( 3.892 , 11.893 )
7 ( 0.944 , 0.944 ) ( 3.007 , 3.367 ) ( 0.710 , 0.243 ) ( 1.264 , 3.082 ) ( 1.535 , 2.536 ) ( 3.136 , 15.924 )
Table 2. Volcano heights data: Statistical summary.
Table 2. Volcano heights data: Statistical summary.
nMeanVariance b 1 b 2
152016.776015.66820.64614.3809
Table 3. Parameter estimates (SE) for the fitted models to the volcano height data.
Table 3. Parameter estimates (SE) for the fitted models to the volcano height data.
Distribution
EstimatesStudent-tSTPTAPST
μ ^ 14.7835(0.3615 )4.7469(0.6892)8.4027(0.7923)11.5509(0.1337)
σ ^ 11.0045(0.3975)14.1532(0.7237)11.8146(0.4707)22.6885(0.0792)
λ ^ 1.5673(0.1838)5.2347(0.2870)
α ^ 1.7912(0.1147)0.3205(0.0347)
ν ^ 3.4156(0.3601)3.4075(0.3454)2.7473(0.2566)12.8734(2.9729)
^ −6273.35−6219.25−6228.77−6205.94
AIC12,552.7012,446.4912,465.5312,421.87
BIC12,568.6812,467.7912,486.5312,448.50
CAIC12,571.6812,471.7912,490.8312,453.50
Table 4. Parameter estimates (SE) for the fitted models to the stellar abundances data.
Table 4. Parameter estimates (SE) for the fitted models to the stellar abundances data.
Distribution
EstimatesCTCSTCPTCAPST
μ ^ 1.0314(0.0010)1.2306(0.0018)1.2098(0.0052)1.1761(0.0054)
σ ^ 0.1596(0.0012)0.2712(0.0058)0.0818(0.0008)0.0905(0.0020)
λ ^ −3.5655(3.7748)0.6580(0.5031)
α ^ 0.1705(0.0208)0.1518(0.0251)
ν ^ 0.9974(0.0884)1.2501(0.1774)6.0927(0.7501)6.0999(0.7326)
^ −29.50743−18.87016−17.67113−14.80241
AIC65.0148745.7403343.3422739.60482
BIC71.6733954.6183652.2203050.70236
CAIC59.3898738.3752535.9771930.57256

Share and Cite

MDPI and ACS Style

Tovar-Falón, R.; Bolfarine, H.; Martínez-Flórez, G. The Asymmetric Alpha-Power Skew-t Distribution. Symmetry 2020, 12, 82. https://doi.org/10.3390/sym12010082

AMA Style

Tovar-Falón R, Bolfarine H, Martínez-Flórez G. The Asymmetric Alpha-Power Skew-t Distribution. Symmetry. 2020; 12(1):82. https://doi.org/10.3390/sym12010082

Chicago/Turabian Style

Tovar-Falón, Roger, Heleno Bolfarine, and Guillermo Martínez-Flórez. 2020. "The Asymmetric Alpha-Power Skew-t Distribution" Symmetry 12, no. 1: 82. https://doi.org/10.3390/sym12010082

APA Style

Tovar-Falón, R., Bolfarine, H., & Martínez-Flórez, G. (2020). The Asymmetric Alpha-Power Skew-t Distribution. Symmetry, 12(1), 82. https://doi.org/10.3390/sym12010082

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