Next Article in Journal
Math Learning in a Science Museum—Proposal for a Workshop Design Based on STEAM Strategy to Learn Mathematics. The Case of the Cryptography Workshop
Next Article in Special Issue
Power Families of Bivariate Proportional Hazard Models
Previous Article in Journal
Interactive Impacts of Overconfidence and Competitive Preference on Closed-Loop Supply Chain Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Alternative to the Log-Skew-Normal Distribution: Properties, Inference, and an Application to Air Pollutant Concentrations

by
Jaime Arrué
1,
Reinaldo Boris Arellano-Valle
2,
Osvaldo Venegas
3,*,
Heleno Bolfarine
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 Estadística, Facultad de Matemáticas, Pontificia Universidad Católica de Chile, Santiago 7820436, 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ítica, IME, Universidade de São Paulo, São Paulo 05508-090, Brazil
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(22), 4336; https://doi.org/10.3390/math10224336
Submission received: 28 September 2022 / Revised: 27 October 2022 / Accepted: 10 November 2022 / Published: 18 November 2022
(This article belongs to the Special Issue Probability, Statistics & Symmetry)

Abstract

:
In this study, we consider an alternative to the log-skew-normal distribution. It is called the modified log-skew-normal distribution and introduces greater flexibility in the skewness and kurtosis parameters. We first study several of the main probabilistic properties of the new distribution, such as the computation of its moments and the non-existence of the moment-generating function. Parameter estimation by the maximum likelihood approach is also studied. This approach presents an overestimation problem in the shape parameter, which in some cases, can even be infinite. However, as we demonstrate, this problem is solved by adapting bias reduction using Firth’s approach. We also show that the modified log-skew-normal model likewise inherits the non-singularity of the Fisher information matrix of the untransformed model, when the shape parameter is null. Finally, we apply the modified log-skew-normal model to a real example related to pollution data.

1. Introduction

The log-normal distribution is commonly used to model the behavior of data with positive asymmetry, in which most of the observations are concentrated near the minimum value. Some applications of the log-normal model are in species abundance patterns, environmental concentrations, stock prices, the distribution of the molecular weights of polymers, the production of copper nano-particles, etc.
The log-normal (LN) distribution arises from the transformation Y = e X , where X N ( μ , σ 2 ) is the normal distribution with mean μ and variance σ 2 and has the density given by
f Y ( y ) = 1 σ y ϕ log y μ σ , y > 0 ,
where ϕ ( · ) is the standard normal density. The notation Y L N ( μ , σ 2 ) is typically used for the log-normal distribution with parameters μ R and σ 2 > 0 . More simply, we can say that a random variable Y L N ( μ , σ 2 ) if and only if log Y N ( μ , σ 2 ) . As noted by Jones and Arnold [1], this distribution is both log-symmetric about its mean/median E ( Y ) = e μ and R-symmetric about its mode δ = e μ + σ 2 , a condition called double symmetry. As defined by Mudholkar and Wang [2], in density terms, a non-negative random variable Y is said to be R-symmetric about δ if f Y ( δ y ) = f Y ( δ / y ) for some δ > 0 and all y > 0 , which is equivalent to δ Y = d δ / Y , where “ = d ” means equal in distribution. When f Y is unimodal, then δ is the (unique) mode of Y, which is less than the mean of Y. Thus, R-symmetric distributions are always positively skewed. Despite these good properties, there are data that are not adequately modeled by the log-normal distribution, since they have symmetry and kurtosis indices that are outside their natural range. A model that presents this characteristic is the log-skew-normal (LSN) distribution, introduced and studied by Azzalini et al. [3]. It is a version with positive support of the skew-normal (SN) distribution, defined as follows: for a skew-normal random variable Z S N ( λ ) , where λ R is the skewness parameter, we say that the random variable Y = e Z follows a log-skew-normal distribution with skewness parameter λ , denoted by Y L S N ( λ ) , if its probability density function (pdf) is given by
f Y ( y ; λ ) = 2 y ϕ ( log ( y ) ) Φ ( λ log ( y ) ) , y > 0 ,
where Φ ( · ) denotes the cumulative distribution function (cdf) of the standard N ( 0 , 1 ) .
Applications of this model to real datasets are reported in Azzalini et al. [3], Marchenko and Genton [4], and Bolfarine et al. [5], where a bimodal version of the model was used to fit a real (bimodal) dataset. Chai and Bailey [6] extended the SN model to a situation of continuous datasets with a discrete component at zero. On the other hand, the modified skew-normal (MSN) distribution is a particular case of the generalized skew-normal (GSN) distribution introduced by Arellano-Valle et al. [7], for λ 2 = 1 . Later, Arellano-Valle et al. [8] investigated the Fisher information matrix (FIM) for the location-scale version of the GSN model. Thus, the MSN model is a fair competitor to the SN model since for both control asymmetry with a single scalar parameter, say λ R , such that, if λ = 0 , then the ordinary normal model results. Moreover, one advantage of the MSN model over the SN model in location-scale situations is that its FIM is nonsingular at λ = 0 , which is not the case with the SN model (see Arellano-Valle et al. [8] and Arrué et al. [9]). The present paper focuses on investigating the possibility of developing a more flexible distribution for positive data, such that the regularity conditions for inference in large samples remain valid when a maximum likelihood approach is used. We considered it natural to use the log version of the MSN for this situation, calling it the modified log-skew-normal(MLSN) model.
This paper is organized as follows. Section 2 presents some probabilistic properties of the MSN model. It includes the derivation of the first few moments of the distribution and the moment-generating function. Plots and ranges for the asymmetry and kurtosis coefficients are reported. The score functions are derived, and the observed information matrix is presented. Section 3 is devoted to the definition of the new distribution, termed the MLSN distribution. The survival and the associated risk functions are derived. A general expression for the moments is obtained, and the non-existence of the moment-generating function is proven. Parameter estimation is conducted by the maximum likelihood (ML) approach. Observed and expected (Fisher) information matrices are derived. It is shown that the FIM for a location-scale extension of the model is non-singular so that large sample properties of the ML estimators (MLE) are satisfied. In Section 4, we present a brief introduction to Firth’s approach (see Firth [10]) for bias reduction and some tables related to a small-scale simulation study, illustrating the amount of bias reduction in the asymmetry parameter. Finally, in Section 5, a real data illustration is presented indicating that the new model outperforms its most direct competitors. The closing section summarizes the main contributions of the paper.

2. Preliminaries

We say that a random variable Z follows a standard MSN distribution with parameter λ , denoted Z M S N ( λ ) , if its pdf is given by
f Z ( z ; λ ) = 2 ϕ ( z ) Φ λ u ( z ) , z R ,
where u ( z ) = z / 1 + z 2 and λ R . If λ = 0 , then the MSN pdf in (2) reduces to the pdf of the standard normal distribution. Non-null values of the parameter λ directly affect the model’s asymmetry so that, in the limit, as λ , the MSN model becomes the half-normal (HN) distribution. All that is required to obtain the location-scale version of the model is to make the transformation X = μ + σ Z , where μ R and σ > 0 are the location and scale parameters, respectively. We use the notation X M S N ( μ , σ , λ ) . As shown in Arellano-Valle et al. [8] the FIM for the location-scale version is nonsingular at λ = 0 . Thus, the ordinary properties of the MLE (consistency and asymptotic normality) remain valid for the MSN model.
Figure 1 shows plots of the pdf for the standard MSN model for several values of λ . The subtle difference between the MSN and SN densities for the same parameter value can be observed, say for the case λ = 1 (dotted line and segmented line).
The following properties follow from Arellano-Valle et al. [7].

2.1. Properties

Let Z MSN ( λ ) ; then:
  • Z MSN ( λ ) Z MSN ( λ ) ;
  • If Z MSN ( λ ) , then | Z | HN ( 0 , 1 ) ;
  • If Z | S = s SN ( s ) and S N ( λ , 1 ) , then Z MSN ( λ ) ;
  • MSN ( 0 ) = d N ( 0 , 1 ) , and MSN ( λ ) = d MGSN ( λ , 1 ) = d GSN ( λ , 1 ) .
According to Property 3, the MSN distribution can be represented as a mixture of the asymmetry parameter between the skew-normal and standard normal distributions. For the location-scale extension, that is X = μ + σ Z , so that X MSN ( μ , σ , λ ) , with pdf
f X ( x ; μ , σ , λ ) = 2 σ ϕ x μ σ Φ λ u x μ σ .

2.2. Moments of the MSN Model

According to Property 3, the moments of the MSN distribution can be obtained by using the fact that it is a mixture between the SN and the N ( 0 , 1 ) so that we can write E ( Z k ) = E ( E ( Z k | S ) ) , where E ( Z k | S ) , k = 0 , 1 , . . . , are the moments for the SN model with asymmetry parameter s. Thus, for even values of k, the moments of the SN ( s ) distribution are constant, coinciding with the moments for the MSN distribution. For the odd moments, making use of the stochastic representation for the skew-normal distribution in Henze [11], we can write
E ( Z 2 k + 1 ) = E ( E ( Z 2 k + 1 | S ) ) = E b ( 2 k + 1 ) ! 2 k ν = 0 k ν ! 2 2 ν ( 2 ν + 1 ) ! ( k ν ) ! S 2 ν + 1 ( 1 + S 2 ) ( 2 k + 1 ) / 2 = b ( 2 k + 1 ) ! 2 k ν = 0 k ν ! 2 2 ν ( 2 ν + 1 ) ! ( k ν ) ! E S 2 ν + 1 ( 1 + S 2 ) ( 2 k + 1 ) / 2 ,
where k = 0 , 1 , 2 , , b = 2 π . Considering
ψ k , ν = ψ k , ν ( λ ) = E S 2 ν + 1 ( 1 + S 2 ) ( 2 k + 1 ) / 2 = s 2 ν + 1 ( 1 + s 2 ) ( 2 k + 1 ) / 2 ϕ ( s λ ) d s
and ψ k = ψ k , k , we have that
E ( Z 2 k + 1 ) = b ( 2 k + 1 ) ! 2 k ν = 0 k ν ! 2 2 ν ( 2 ν + 1 ) ! ( k ν ) ! ψ k , ν .
Therefore, the first four moments of the standard MSN distribution are given by
E ( Z ) = b ψ 0 , E ( Z 2 ) = 1 , E ( Z 3 ) = b ( 3 ψ 0 ψ 1 ) a n d E ( Z 4 ) = 3 .
The odd moments can also be obtained from Arellano-Valle et al. [7]:
E ( Z 2 k + 1 ) = 2 c k 2 k k ! b ,
for k = 0 , 1 , 2 , where c k : = c k ( λ ) = 0 x k ϕ ( x ) Φ ( λ u ( x ) ) d x .

2.3. Moment-Generating Function

As mentioned previously, if Z | S = s S N ( s ) and S N ( λ , 1 ) , then Z MSN ( λ ) , so that we can write
M Z ( t ) = E ( E ( e Z t | S ) ) = E ( M Z | S ( t ) ) = 2 e t 2 / 2 E Φ S 1 + S 2 t .
Therefore, all the moments for the random variable Z are defined as follows:
E ( Z n ) = n M Z ( t ) n t t = 0 = E n M Z | S ( t ) n t t = 0 .

2.4. Observed Information Matrix for MSN Model

Consider a random sample x 1 , x 2 , . . , x n from the MSN( θ ) distribution, with θ = ( μ , σ , λ ) , so that the corresponding log-likelihood function is given by
l ( θ ) = n 2 log 2 π n log ( σ ) 1 2 i = 1 n z i 2 + i = 1 n log ( Φ ( λ u ( z i ) ) ,
where z i = ( x i μ ) / σ , u ( x ) = x / 1 + x 2 , μ R , σ > 0 , λ R , and x R . The entries of the observed information matrix for the MSN model have the same structure as the corresponding moments for the MLSN model considered in Section 4.3, with the appropriate transformation.

3. The Modified Log-Skew-Normal Distribution

If Z MSN ( λ ) , then we have that Y = e Z is distributed according to the standard MLSNdistribution with parameter λ , denoted by Y MLSN ( λ ) , if its pdf is given by
f Y ( y ; λ ) = 2 y ϕ ( log ( y ) ) Φ ( λ u ( log ( y ) ) ) ,
where u ( x ) = x / 1 + x 2 , y R + , and λ R .
Figure 2 depicts plots of the pdf for the MLSN( λ ) model for several values of λ . If λ = 0 , then it coincides with the log-normal distribution (solid line), illustrating the fact that the MLSN model is an extension of the LN model. Concerning the location-scale situation, that is Y MLSN ( μ , σ , λ ) , where μ R , σ > 0 , and λ R , its density is given by
f Y ( y ; μ , σ , λ ) = 2 y σ ϕ ( z ) Φ ( λ u ( z ) ) ,
where z = ( log ( y ) μ ) / σ , u ( x ) = x / 1 + x 2 , and y R + . In a survival analysis scenario, it is important to study the following functions: the survival function S ( t ) = 1 F ( t ) and the risk function r ( t ) = f ( t ) / S ( t ) , which for the model under study can be shown to be given by
S ( t ) = 1 0 t 2 w ϕ ( log ( w ) ) Φ ( λ u ( log ( w ) ) ) d w
and
r ( t ) = 2 t ϕ ( log ( t ) ) Φ ( λ u ( log ( t ) ) ) 1 0 t 2 w ϕ ( log ( w ) ) Φ ( λ u ( log ( w ) ) ) d w .
Clearly, using L’Hopital´s rule, it follows that r ( t ) 0 as t , as can also be appreciated graphically.
Figure 3 illustrates the behavior of the risk function r ( t ) for some values of λ . The maximum values of the risk function for each λ decrease for λ ( , 1.4484 ) and increase otherwise. Moreover, λ = , r ( t ) tends to a strictly increasing function defined in the interval ( 0 , 1 ] , otherwise being zero. On the other hand, λ = , r ( t ) is a strictly increasing function defined in the interval [ 1 , ) , coinciding with the risk function of the log-normal model in this interval and taking the value zero in the interval ( 0 , 1 ) .

3.1. Moments for the MLSN Model

The r-th moment is given by
μ r = 2 e r 2 / 2 E Φ S 1 + S 2 r ,
for r = 0 , 1 , 2 . . and S N ( λ , 1 ) . This expression is obtained directly from the moment-generating function of the MSN model, given in (5), since it is valid for all t > 0 , particularly for t = r . Alternatively, the moments can be obtained from
μ r = E ( Y r ) = 0 y r 2 y ϕ ( log ( y ) ) Φ ( λ u ( log ( y ) ) ) d y .

3.1.1. Non-Existence of the Moment-Generating Function for the MLSN Distribution

Proposition 1. 
For all λ R , variable Y M L S N ( λ ) has no moment-generating function.
Proof. 
Part of the proof parallels that of Lin and Stoyanov [12] for a similar proof. Thus, for each t > 0 ,
E ( e t Y ) = 0 h λ ( y ) d y ,
where
h λ ( y ) = e t y y ϕ ( log ( y ) ) Φ ( λ u ( log ( y ) ) ) > 0 , y > 0 ,
so that, for
λ 0 , lim y Φ ( λ u ( log ( y ) ) ) 1 2
and
λ < 0 , lim y Φ ( λ u ( log ( y ) ) ) Φ ( λ ) ,
in both cases, h λ ( y ) as y . Therefore, given t > 0 , E ( e t Y ) = , for any λ . □

3.1.2. Skewness and Kurtosis Coefficients

To obtain expressions for the kurtosis and skewness coefficients, we have to compute the central moments using the following relationships:
μ 2 = μ 2 μ 1 2 , μ 3 = μ 3 3 μ 1 μ 2 + 2 μ 1 3 a n d μ 4 = μ 4 4 μ 1 μ 3 + 6 μ 1 2 μ 2 3 μ 1 4 .
This allows us to compute the variance, standard deviation (SD), coefficient of variation (CV), asymmetry ( γ 1 ), and kurtosis ( γ 2 ) coefficients, respectively, given by:
μ 2 = V a r ( Y ) , C V ( Y ) = S D ( Y ) μ 2 , γ 1 = μ 3 ( μ 2 ) 3 / 2 and γ 2 = μ 4 ( μ 2 ) 2 .
Table 1 reveals that the variation ranges for the variance and CV are relatively short, while the variation range for the asymmetry and kurtosis coefficients are relatively long if compared with the ordinary LSN model.
Figure 4 shows the behavior of the variance, CV, skewness, and kurtosis as a function of the lambda parameter. It can be observed that the minimum values correspond to the asymptotes of the left tail. On the other hand, for small values of λ , say λ > 4 , the plots of the right tail stabilize around the horizontal asymptotes with values of 6.74, 0.936, 5.83, and 97.93, respectively, for the corresponding indices.

3.2. MLE for the MLSN Model

Given a random sample y 1 , y 2 , . . , y n from a random variable Y M L S N ( θ ) , with θ = ( μ , σ , λ ) , so that the corresponding log-likelihood is given by
l ( θ ) = n 2 log 2 π n log ( σ ) i = 1 n log ( y i ) 1 2 i = 1 n z i 2 + i = 1 n log ( Φ ( λ u ( z i ) ) ,
where z i = ( log ( y i ) μ ) / σ , u ( x ) = x / 1 + x 2 , μ R , σ > 0 , λ R , and y R + . Using the following notation:
ρ n m i = z i n ( 1 + z i 2 ) m / 2 ζ ( λ u ( z i ) ) , η n m i = z i n ( 1 + z i 2 ) m ζ 2 ( λ u ( z i ) ) , η ¯ n m = 1 n i = 1 n ρ n m i a n d ζ ( x ) = ϕ ( x ) Φ ( x ) ,
we have that the associated scoring vector is given by
S μ = n σ z ¯ λ ρ ¯ 03 , S σ = n σ 1 + z 2 ¯ λ ρ ¯ 13 a n d S λ = n ρ ¯ 11 .
Equating the scoring functions to zero, it follows that the likelihood equations are given as
z ¯ = λ ρ ¯ 03 z 2 ¯ 1 = λ ρ ¯ 13 and ρ ¯ 11 = 0 .
Solving this system of equations, which require numerical procedures, leads to the MLE for μ , σ , and λ .

4. Bias Prevention of the MLE

4.1. Abstract for Firth’s Approach

One problem with the ML approach is that the likelihood function can be unbounded for the parameter λ . In the standard case (just λ unknown), the MLE of λ is infinite given that all observations are positive. In the location-scale situation, the same happens because all observations are greater than μ , the location parameter. Therefore, the MLE overestimates parameter λ , which is always of concern in a statistical analysis. The percentage of overestimation of the shape parameter λ depends on the sample size n and on the true parameter value. It is of interest to reduce the amount of bias in a statistical problem. One approach considered previously in the literature was proposed by Firth [10] and had previously been used by Sartori [13] for the skew-normal and skew-t (ST) distributions. The reduced bias estimator is of order O ( n 2 ) (see also Cox and Snell [14]). To arrive at the reducedbias estimator, the following notation is used: Let l ( θ ) be the log-likelihood function for a parametric regular family of distributions, and denote by U ( θ ) = l ( θ ) the score function, j ( θ ) = l ( θ ) the second derivative of the log-likelihood function, and l its first derivative. We also consider the expectation of those functions that are of order O ( n ) . Therefore, we have
i ( θ ) = E θ { j ( θ ) } , ν θ , θ , θ ( θ ) = E θ { l ( θ ) 3 } , ν θ , θ θ ( θ ) = E θ { l ( θ ) l ( θ ) } .
Briefly speaking, Firth’s approach consists of modifying the score function keeping the score unchanged, that is
U M ( θ ) = U ( θ ) + M ( θ ) ,
where M ( θ ) is O ( 1 ) as n . The solution to the modified likelihood equation U M ( θ ) = 0 produces the modified MLE, say θ ^ M with bias of order O ( n 2 ) . The modification factor M ( θ ) is given by
M ( θ ) = i ( θ ) b ( θ ) = 1 2 i ( θ ) 1 ( ν θ , θ , θ + ν θ , θ θ ) .
From the modified score function, the modified quasi-log-likelihood function is defined as
l M ( θ ) = c θ U M ( t ) d t = l ( θ ) l ( c ) + c θ M ( t ) d t ,
where c is an arbitrary constant. This penalized log-likelihood function has a penalty of order O ( 1 ) . It is also possible to define a modified log-likelihood ratio statistic:
W M ( θ ) = 2 { l M ( θ ^ M ) l M ( θ ) } ,
which is asymptotically distributed as a χ 1 2 distribution, which can be used for testing hypotheses about θ , as well as for confidence interval construction.

4.2. Simulation Study

Below, we present a simulation study based on 5000 iterations for μ = 0 , σ = 1 , λ = 5 , 10 , and several sample sizes. Table 2 shows the bias for the modified and ordinary MLE, the empirical coverage of the confidence intervals based on W M P ( λ ) , and the percentage of cases in which the un-modified MLE for λ is finite. It can be observed that the MLEs μ ^ and σ ^ present satisfactory behavior concerning bias, and they seem to be unaffected by the sometimes erratic behavior of estimator λ ^ , so that Firth’s method is only applied for the parameter λ .
The notation W M P ( λ ) = 2 { l M P ( λ ^ M ) l M P ( λ ) } defines the likelihood ratio statistic based on the profile log-likelihood function l P ( λ ) = l ( μ ^ ( λ ) , σ ^ ( λ ) , λ ) , where μ ^ and σ ^ are the MLEs for μ and σ , respectively, for a fixed value of λ .

4.3. Observed Information Matrix for the MLSN Model

The entries for the observed information matrix of the MLSN distribution, corresponding to the log-likelihood given in (7), are J θ i θ j = 2 l ( θ ) θ i θ j , for i , j = 1 , 2 , 3 (see Appendix A):
J μ μ = n σ 2 1 + 3 λ ρ ¯ 15 + λ 3 ρ ¯ 17 + λ 2 η ¯ 03 , J μ σ = n σ 2 2 z ¯ + λ ρ ¯ 05 2 λ ρ ¯ 25 λ 3 ^ ρ ¯ 27 λ 2 η ¯ 13 ,
J μ λ = n σ ρ ¯ 03 λ 2 ρ ¯ 25 λ η ¯ 12 , J σ σ = n σ 2 1 3 z 2 ¯ + λ ρ ¯ 13 + λ ρ ¯ 15 2 λ ρ ¯ 35 λ 3 ρ ¯ 37 λ 2 η ¯ 23 ,
J σ λ = n σ ρ ¯ 13 + λ 2 ρ ¯ 35 + λ η ¯ 22 , J λ λ = n ( λ ρ ¯ 33 + η ¯ 21 ) .

4.4. Fisher Information

For a random sample y 1 , y 2 , . . , y n from the MLSN ( θ ) distribution, with θ = ( μ , σ , λ ) , the FIM associated with (7) is given by the following entries (see Appendix A):
I μ μ = 1 σ 2 1 + λ 2 η 03 , I σ μ = 1 σ 2 2 E ( z ) + λ ρ 05 2 λ ρ 25 λ 3 ρ 27 λ 2 η 13 ,
I μ λ = 1 σ ρ 03 λ 2 ρ 25 λ η 12 , I σ σ = 1 σ 2 2 + λ 2 η 23 , I λ σ = λ σ η 22 , I λ λ = η 21 ,
where
ρ n m = E z n ( 1 + z 2 ) m / 2 ζ , ρ n m = 0 i f n = odd   number η n m = E z n ( 1 + z 2 ) m ζ 2 , with ζ : = ζ ( λ u ( z ) ) = ϕ ( λ u ( z ) ) Φ ( λ u ( z ) ) .
For the case λ = 0 , we have that
FIM = 1 σ 2 0 d 1 σ 0 2 σ 2 0 d 1 σ 0 d 2 ,
where d 1 = 2 ( 2 / π ) 1 / 2 0 ϕ ( x ) ( 1 + x 2 ) 3 / 2 d x and d 2 = 2 π ( 1 ( 2 π ) 1 / 2 e 1 / 2 Φ ( 1 ) ) have to be computed numerically. It is clearly seen that the above FIM is non-singular at λ = 0 .

5. An Application

The dataset analyzed in this section was previously studied in Nadarajah [15] and Leiva et al. [16]. It consists of daily measurements of ozone concentration (in p p b = p p m × 1000 ) in New York city between May and September 1973. The data were supposed to be independent, without the presence of tendencies or cyclical components (see Gokhale and Khare [17]). Table 3 presents the summary statistics, in particular the asymmetry and kurtosis coefficients, which are represented by the sample quantities ( b 1 ) and ( b 2 ), respectively.
Table 4 shows the MLE for the three parameters for the MLSN, LN, and LSN distributions, respectively, where values in parentheses correspond to standard errors. Note also from the table that the AIC for the LSN model is slightly smaller than that for the MLSN model; it is our opinion that the latter model should be preferred because the LSN model has a singular covariance matrix, and so, the likelihood ratio statistics for testing H 0 : λ = 0 is not distributed according to a central chi-squared distribution in large samples. On the other hand, under the MLSN model, for which the FIM is non-singular, we have that the hypothesis H 0 : λ = 0 is rejected at the 5 % level, using the likelihood ratio statistics; this is given by
Λ = L L N ( μ ^ , σ ^ ) L M L S N ( μ ^ , σ ^ , λ ^ ) .
Replacing the MLE from Table 4, we have that 2 l o g Λ = 2 ( 543.883 + 542.105 ) = 9.516 , greater than the critical value χ 1 2 = 3.84 . Furthermore, from the table, a large sample 95% confidence interval for λ based on the MLSN model does not contain λ = 0 .
Figure 5 presents the data histogram with the corresponding fitted fdp for the MLSN ( μ ^ , σ ^ , λ ^ ) (solid line) and LN ( μ ^ , σ ^ , λ ^ ) distributions (dotted line) and the fit cdf for the MLSN and LN models, jointly with the empirical cdf. It can be shown from the figure that the MLSN model seems to provide a (graphically) satisfactory fit.
Graphical corroboration of the better fit of the MLSN model than the LN model is also illustrated with the QQ-plots in Figure 6.
Table 5 presents the MLEs μ ^ , σ ^ , and λ ^ and the modified λ ^ M with the corresponding standard errors (in parentheses), obtained by using the estimated FIM for the MLSN model, bearing in mind that the asymptotic distribution of θ ^ is N 3 ( θ , I ( θ ) 1 / n ) , where θ = ( μ , σ , λ ) or θ = ( μ , σ , λ M ) . The table indicates that the modified MLE λ ^ M is greater than the ordinary MLE λ ^ and is expected to have smaller bias.
Table 6 presents confidence intervals for λ ^ and λ ^ M for several confidence coefficients. Comparing the lengths of the intervals for the two estimators, there is strong evidence that the modified estimator λ ^ M presents shorter intervals.

6. Concluding Remarks

This paper focused on a transformation of the MSN model (Arellano-Valle et al. [8]), which led to a more flexible distribution (wider ranges for asymmetry and kurtosis). This model, which we call the MLSN distribution, is suitable for positive data, its main competitors being the LN and LSN distributions. One interesting aspect of the new model is that its FIM is non-singular so that the large sample theory for the MLE remains valid. This is not the case with the LSN model, for which the FIM is singular in the location-scale version. Thus, in particular for testing H 0 : L N against H 0 : M L S N , under the MLSN model, the likelihood ratio statistics in large samples are distributed as in the chi-squared distribution. Large sample confidence intervals could also be constructed and used for testing the hypothesis H 0 : λ = 0 , where λ is the skewness parameter. Rejection of H 0 indicates that the MLSN model should be preferred. It is also noticed with the simulation study that the MLE of λ can overestimate λ (which could be infinite for some samples). Thus, the bias-reducing approach of Firth [10] was used to derive a less-biased estimator. Estimations for the location and scale parameters remain stable, however, and do not need to be corrected. An application to a real dataset revealed that the new model can be a valuable alternative for modeling positive data.

Author Contributions

Conceptualization, J.A., R.B.A.-V., and H.W.G.; methodology, R.B.A.-V., H.B., and H.W.G.; software, J.A. and O.V.; validation, R.B.A.-V., H.B., and H.W.G.; formal analysis, R.B.A.-V. and H.W.G.; investigation, J.A., O.V., and H.W.G.; writing—original draft preparation, J.A. and R.B.A.-V.; writing—review and editing, J.A. and O.V.; funding acquisition, H.W.G. and O.V. All authors have read and agreed to the published version of the manuscript.

Funding

The research of Jaime Arrué and Héctor W. Gómez was supported by SEMILLERO UA-2022 (Chile).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset is available in [15,16].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Considering the notation:
θ = ( μ , σ ) , z = x μ σ , u ( x ) = x 1 + x 2 , ζ = ζ ( λ u ( z ) ) = ϕ ( λ u ( z ) ) Φ ( λ u ( z ) ) ,
and the following derivatives:
z μ = 1 σ , z σ = z σ , u ( z ) θ = 1 ( 1 + z 2 ) 3 / 2 z θ ,
ζ θ = λ 2 z ( 1 + z 2 ) 2 ζ λ ( 1 + z 2 ) 3 / 2 ζ 2 z θ , θ z ( 1 + z 2 ) 3 / 2 = 1 2 z 2 ( 1 + z 2 ) 5 / 2 z θ , ζ λ = λ z 2 ( 1 + z 2 ) ζ z ( 1 + z 2 ) 1 / 2 ζ 2 ,
we have that the score functions score for μ , σ , and λ for the MSN ( λ ) model are given by
S μ = 1 σ z λ ( 1 + z 2 ) 3 / 2 ζ , S σ = 1 σ 1 + z 2 λ z ( 1 + z 2 ) 3 / 2 ζ , S λ = z ( 1 + z 2 ) 1 / 2 ζ ,
so that the entries of the FIM are given by:
I μ μ = E ( S μ μ ) = 1 σ 2 E 1 3 λ z ( 1 + z 2 ) 5 / 2 ζ λ 3 z ( 1 + z 2 ) 7 / 2 ζ λ 2 ( 1 + z 2 ) 3 ζ 2 = 1 σ 2 1 + λ 2 η 03 , I μ σ = E ( S σ μ ) = 1 σ 2 E 2 z + λ ( 1 2 z 2 ) ( 1 + z 2 ) 5 / 2 ζ λ 3 z 2 ( 1 + z 2 ) 7 / 2 ζ λ 2 z ( 1 + z 2 ) 3 ζ 2 = 1 σ 2 2 E ( z ) + λ ρ 05 2 λ ρ 25 λ 3 ρ 27 λ 2 η 13 , I μ λ = E ( S μ λ ) = 1 σ E 1 ( 1 + z 2 ) 3 / 2 ζ λ 2 z 2 ( 1 + z 2 ) ζ λ z ( 1 + z 2 ) 1 / 2 ζ 2 = 1 σ ρ 03 λ 2 ρ 25 λ η 12 , I σ σ = E ( S σ σ ) = 1 σ 2 E 1 3 z 2 + λ z ( 1 + z 2 ) 3 / 2 ζ + λ ( z 2 z 3 ) ( 1 + z 2 ) 5 / 2 ζ λ 3 z 3 ( 1 + z 2 ) 7 / 2 ζ λ 2 z 2 ( 1 + z 2 ) 3 ζ 2 = 1 σ 2 2 + λ 2 η 23 , I σ λ = E ( S λ σ ) = 1 σ E z ( 1 + z 2 ) 3 / 2 ζ + λ 2 z 3 ( 1 + z 2 ) 5 / 2 ζ + λ z 2 ( 1 + z 2 ) 2 ζ 2 = λ σ η 22 , I λ λ = E ( S λ λ ) = E λ z 3 ( 1 + z 2 ) 3 / 2 ζ z 2 1 + z 2 ζ 2 = η 21 .

References

  1. Jones, M.C.; Arnold, B.C. Distributions that are both log-symmetric and R-symmetric. Electron. J. Stat. 2008, 2, 1300–1308. [Google Scholar] [CrossRef] [Green Version]
  2. Mudholkar, G.S.; Wang, H. IG-symmetry and R-symmetry: Interrelations and applications to the inverse Gaussian theory. J. Stat. Plan. Inference 2007, 137, 3655–3671. [Google Scholar] [CrossRef]
  3. Azzalini, A.; Cappello, D.; Kotz, S. Log-skew-normal and log-skew-t distributions as models for family income data. J. Income Distrib. 2003, 11, 12–20. [Google Scholar] [CrossRef]
  4. Marchenko, Y.V.; Genton, M.G. Multivariate log-skew-elliptical distributions with applications to precipitation data. Environmetrics 2010, 21, 318–340. [Google Scholar] [CrossRef]
  5. Bolfarine, H.; Gómez, H.W.; Rivas, L. The log-bimodal-skew-normal model. A geochemical application. J. Chemom. 2011, 25, 329–332. [Google Scholar] [CrossRef]
  6. Chai, H.; Bailey, K. Use of log-normal distribution in analysis of continuous data with a discrete component at zero. Stat. Med. 2008, 27, 3643–3655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Arellano-Valle, R.B.; Gómez, H.W.; Quintana, F.A. A new class of skew-normal distributions. Commun. Stat.-Theory Methods 2004, 33, 1465–1480. [Google Scholar] [CrossRef]
  8. Arellano-Valle, R.B.; Gómez, H.W.; Salinas, H.S. A note on the Fisher information matrix for the skew-generalized-normal model. SORT 2013, 37, 19–28. [Google Scholar]
  9. Arrué, J.; Arellano-Valle, R.B.; Gómez, H.W. Bias reduction of maximum likelihood estimates for a modified skew normal distribution. J. Stat. Comput. Simul. 2016, 86, 2967–2984. [Google Scholar] [CrossRef]
  10. Firth, D. Bias reduction of maximum likelihood estimates. Biometrika 1993, 80, 27–38, Amended in Biometrika 1995, 82, 667. [Google Scholar] [CrossRef]
  11. Henze, N. A probabilistic representation of the skew-normal distribution. Scand. J. Stat. 1986, 13, 271–275. [Google Scholar]
  12. Lin, G.D.; Stoyanov, J. The logarithmic skew-normal distributions are moment-indeterminate. J. Appl. Prob. 2009, 46, 909–916. [Google Scholar] [CrossRef]
  13. Sartori, N. Bias prevention of maximum likelihood estimates for scalar skew normal and skew t distributions. J. Stat. Plan. Inference 2006, 136, 4259–4275. [Google Scholar] [CrossRef]
  14. Cox, D.R.; Snell, E.J. A general definition of residuals. J. R. Stat. Soc. Ser. B 1968, 30, 248–275. [Google Scholar] [CrossRef]
  15. Nadarajah, S. A truncated inverted beta distribution with application to air pollution data. Stoch. Environ. Res. Risk Assess. 2008, 22, 285–289. [Google Scholar] [CrossRef]
  16. Leiva, V.; Vilca-Labra, F.; Balakrishnan, N.; Sanhueza, A. A Skewed Sinh-Normal Distribution and its Properties and Application to Air Pollution. Commun. Stat.-Theory Methods 2010, 39, 426–443. [Google Scholar] [CrossRef]
  17. Gokhale, S.; Khare, M. Statistical behavior of carbon monoxide from vehicular exhausts in urban environments. Environ. Model. Softw. 2007, 22, 526–535. [Google Scholar] [CrossRef]
Figure 1. Plots of the MSN ( λ ) model for different values of λ : MSN ( 0 ) (solid line), MSN ( 1 ) (dotted line), SN ( 1 ) (dashed line), and MSN ( 3 ) (thick solid line) models.
Figure 1. Plots of the MSN ( λ ) model for different values of λ : MSN ( 0 ) (solid line), MSN ( 1 ) (dotted line), SN ( 1 ) (dashed line), and MSN ( 3 ) (thick solid line) models.
Mathematics 10 04336 g001
Figure 2. Plots of the MLSN ( λ ) model for different values of λ .
Figure 2. Plots of the MLSN ( λ ) model for different values of λ .
Mathematics 10 04336 g002
Figure 3. Plots of the risk function of the MLSN model for different values of λ .
Figure 3. Plots of the risk function of the MLSN model for different values of λ .
Mathematics 10 04336 g003
Figure 4. Plots of the variance, CV, skewness, and kurtosis.
Figure 4. Plots of the variance, CV, skewness, and kurtosis.
Mathematics 10 04336 g004
Figure 5. Left panel: Kernel density plot of the ozone concentration level data fit by the LN and MLSN distributions. Right panel: Empirical cdf of the LN (dotted line) and MLSN (solid line) models, whose parameters are estimated by the ML method for ozone concentration level data.
Figure 5. Left panel: Kernel density plot of the ozone concentration level data fit by the LN and MLSN distributions. Right panel: Empirical cdf of the LN (dotted line) and MLSN (solid line) models, whose parameters are estimated by the ML method for ozone concentration level data.
Mathematics 10 04336 g005
Figure 6. Left panel: QQ-plot for the MLSN model. Right panel: QQ-plot for the LN model.
Figure 6. Left panel: QQ-plot for the MLSN model. Right panel: QQ-plot for the LN model.
Mathematics 10 04336 g006
Table 1. Range of values for statistics corresponding to variance, CV, skewness, and kurtosis.
Table 1. Range of values for statistics corresponding to variance, CV, skewness, and kurtosis.
StatisticsVarianceCV γ 1 γ 2
Minimum0.0620.4770.1101.952
Maximum6.9861.39510.985524.351
Table 2. Results are based on 5000 samples simulated from the MSN ( 0 , 1 , λ ) model.
Table 2. Results are based on 5000 samples simulated from the MSN ( 0 , 1 , λ ) model.
λ n μ ^ ( sd ) σ ^ ( sd ) λ ^ a ( sd ) λ ^ M ( sd ) W M P ( λ ) IC % ( λ ^ < )
550−0.003510.998537.001034.3595496.1686.48
(0.10073)(0.12670)(4.63146)(1.18831)
100−0.002080.998736.441164.8031594.1897.84
(0.07114)(0.08973)(2.91069)(0.93203)
2000.000380.998385.536224.8974494.0499.98
(0.05080)(0.06355)(1.43447)(0.42190)
10500.015990.9802611.360526.5616087.9866.00
(0.07556)(0.11462)(9.14682)(2.39887)
1000.003590.9931913.526248.7176192.0887.74
(0.04890)(0.08054)(8.63025)(2.92047)
2000.000950.9974812.950049.5494493.3898.42
(0.03409)(0.05698)(6.00251)(1.74399)
μ ^ , σ ^ , λ ^ , and   λ ^ M are the estimates of the true values of the parameters, and sd denotes the respective standard deviations; the empirical coverage of the 0.95 confidence interval is based on W M P ( λ ) and the empirical (theoretical) percentage for the situation that λ ^ exists. a Computed as λ ^ < .
Table 3. Descriptive statistics of ozone concentration level measurement data.
Table 3. Descriptive statistics of ozone concentration level measurement data.
DatanMean sd b 1 b 2
Ozone116 42.129 32.987 1.209 4.112
Table 4. MLEs for the MLSN, LN, and LSN models.
Table 4. MLEs for the MLSN, LN, and LSN models.
MLEsMLSNLNLSN
μ ^ 4.532 (0.105)3.418 (0.080)4.372 (0.234)
σ ^ 1.411 (0.121)0.861 (0.056)1.285 (0.193)
λ ^ −4.095 (1.057)-−2.411 (1.396)
log-likelihood−542.105−543.883−541.655
AIC1090.2111091.7661089.31
Table 5. MLEs of the μ , σ , λ , and λ M parameters.
Table 5. MLEs of the μ , σ , λ , and λ M parameters.
μ ^ σ ^ λ ^ λ ^ M l ( μ ^ , σ ^ , λ ^ ) l ( μ ^ , σ ^ , λ ^ M )
4.532 ( 0.105 ) 1.411 ( 0.121 ) 4.095 ( 1.057 ) 542.105
4.532 ( 0.123 ) 1.411 ( 0.130 ) 3.342 ( 0.791 ) 542.663
Table 6. Confidence intervals for the estimated values of λ ^ and λ ^ M .
Table 6. Confidence intervals for the estimated values of λ ^ and λ ^ M .
MLE 95 % 98 % 99 %
λ ^ ( 5.834 , 2.356 ) ( 6.266 , 1.924 ) ( 6.554 , 1.636 )
λ ^ M ( 4.644 , 2.040 ) ( 4.968 , 1.716 ) ( 5.183 , 1.500 )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arrué, J.; Arellano-Valle, R.B.; Venegas, O.; Bolfarine, H.; Gómez, H.W. An Alternative to the Log-Skew-Normal Distribution: Properties, Inference, and an Application to Air Pollutant Concentrations. Mathematics 2022, 10, 4336. https://doi.org/10.3390/math10224336

AMA Style

Arrué J, Arellano-Valle RB, Venegas O, Bolfarine H, Gómez HW. An Alternative to the Log-Skew-Normal Distribution: Properties, Inference, and an Application to Air Pollutant Concentrations. Mathematics. 2022; 10(22):4336. https://doi.org/10.3390/math10224336

Chicago/Turabian Style

Arrué, Jaime, Reinaldo Boris Arellano-Valle, Osvaldo Venegas, Heleno Bolfarine, and Héctor W. Gómez. 2022. "An Alternative to the Log-Skew-Normal Distribution: Properties, Inference, and an Application to Air Pollutant Concentrations" Mathematics 10, no. 22: 4336. https://doi.org/10.3390/math10224336

APA Style

Arrué, J., Arellano-Valle, R. B., Venegas, O., Bolfarine, H., & Gómez, H. W. (2022). An Alternative to the Log-Skew-Normal Distribution: Properties, Inference, and an Application to Air Pollutant Concentrations. Mathematics, 10(22), 4336. https://doi.org/10.3390/math10224336

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