Next Article in Journal
A Complex Matrix Private Key to Enhance the Security Level of Image Cryptography
Previous Article in Journal
Fourth Hankel Determinant Problem Based on Certain Analytic Functions
Previous Article in Special Issue
An Asymmetric Bimodal Double Regression Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Bimodal Model Based on Truncation Positive Normal with Application to Height Data

1
Departamento de Ciencias Matemáticas y Físicas, Facultad de Ingeniería, Universidad Católica de Temuco, Temuco 4780000, Chile
2
Departamento de Matemática, Facultad de Ingeniería, Universidad de Atacama, Copiapó 1530000, Chile
3
Department of Statistics, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora 36036-900, Brazil
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(4), 665; https://doi.org/10.3390/sym14040665
Submission received: 13 February 2022 / Revised: 17 March 2022 / Accepted: 17 March 2022 / Published: 24 March 2022
(This article belongs to the Special Issue Symmetric and Asymmetric Bimodal Distributions with Applications)

Abstract

:
In this work, we propose a new bimodal distribution with support in the real line. We obtain some properties of the model, such as moments, quantiles, and mode, among others. The computational implementation of the model is presented in the tpn package of the software R. We perform a simulation study in order to assess the properties of the maximum likelihood estimators in finite samples. Finally, we present an application to a bimodal data set, where our proposal is compared with other models in the literature.

1. Introduction

Describing a phenomenon by a probability distribution is very useful because of the properties associated with it: expectation, shape, range, etc. However, this description can be difficult when a phenomenon (in practice, an observed dataset) is bimodal, which occurs commonly in areas like astrophysics, ecology and genetics; see [1,2,3], respectively. The first approach to fit a bimodal data is using a mixture of two unimodal distributions, for instance, a mixture of gaussian distributions; see [4]. The main disadvantage of this procedure is the non-identifiability of the proposed mixture model. The second and the most workable practical approach is to use distributions which already have bimodal properties. Because of these properties, there is an increasing interest to derive bimodal distributions in the literature: refs. [5,6] presented extensions of the skew-normal, ref. [7] proposed a generalization of the Burr type X distribution and [8] derived an extension of the sinh Cauchy distribution. In this paper, we will discuss an extension of the half normal distribution proposed by [9], the truncated positive normal (tpn) model. The probability density function (pdf) for the tpn model is given by
f ( x ; σ , λ ) = 1 σ Φ ( λ ) ϕ x σ λ , x , σ + , λ ,
where σ and λ are the scale and shape parameters, respectively, and ϕ ( · ) and Φ ( · ) are the pdf and cumulative distribution function (cdf), respectively, of the standard normal distribution. The corresponding cdf of the tpn model is
F X ( x ; σ , λ ) = Φ ( x σ λ ) + Φ ( λ ) 1 Φ ( λ ) .
Note that the cdf above has a closed-form expression, which is useful for generating random data besides defining quantiles. For more properties of the tpn model, see [9]. The restriction to positive values is a limitation of the tpn model. To overcome this limitation, the chief goal of this paper is to derive an extension of the tpn model which has support in the real line. We describe in detail the model, studying its main properties and related functions. Moreover, we show analytically the regions in which the model is unimodal and bimodal, and such regions depend only on one parameter.
The paper is organized as follows. In Section 2, we derive an extension of the tpn with support in the real line and study some properties of the distribution. The inference for parameter estimation in the proposed model and computational aspects are presented in Section 3. In Section 4, we perform a simulation study to evaluate the parameter estimation in finite samples. An application to real data is discussed in Section 5. Finally, conclusions are given in Section 6.

2. A Bimodal Truncation Positive Normal Distribution

In this section, we present the stochastic representation for the bimodal truncation positive normal (btpn) distribution and some properties, such as its pdf and its cdf. We also discuss some particular cases of the model.

2.1. Stochastic Representation, pdf and cdf

Let T be a discrete random variable such as
T = { ( 1 + ϵ ) , with   probability   ( 1 + ϵ ) / 2 1 ϵ , with   probability   ( 1 ϵ ) / 2 ,
where ϵ ( 1 , 1 ) . If Z tpn ( σ , λ ) , independent from T, then we define a new random variable given by X = Z T . We say that X follows a btpn distribution.
Proposition 1.
The pdf for the btpn distribution is given by
f ( x ; σ , λ , ϵ ) = { ϕ ( x σ ( 1 + ϵ ) + λ ) 2 σ Φ ( λ ) , i f   x < 0 ϕ ( x σ ( 1 ϵ ) λ ) 2 σ Φ ( λ ) , i f   x 0 ,
where σ > 0 , λ R and ϵ ( 1 , 1 ) .
Proof. 
If x < 0 , then the cdf for X is
F X ( x ) = P ( X x ) = ( 1 + ϵ ) 2 P ( z x 1 + ϵ ) = ( 1 + ϵ ) 2 [ 1 P ( z x 1 + ϵ ) ] = ( 1 + ϵ ) 2 [ 1 F Z ( x 1 + ϵ ) ] .
Deriving the last expression in relation to x, we have
f X ( x ) = ( 1 + ϵ ) 2 [ f Z ( x 1 + ϵ ) 1 ( 1 + ϵ ) ] = 1 2 [ 1 σ Φ ( λ ) ϕ ( x ( 1 + ϵ ) σ λ ) ] = ϕ ( x σ ( 1 + ϵ ) + λ ) 2 σ Φ ( λ ) .
A similar routine calculation shows that for x 0 , we have that
f X ( x ) = ϕ ( x σ ( 1 ϵ ) λ ) 2 σ Φ ( λ ) ,
completing the proof. □
Figure 1 shows the pdf function for the btpn model with different combination of parameters. Note that the model can assume different shapes, including unimodal, bimodal, symmetric and asymmetric.
Proposition 2.
The cdf of X b t p n ( σ , λ , ϵ ) is given by
F X ( x ) = { ( 1 + ϵ ) 2 Φ ( λ ) Φ ( x σ ( 1 + ϵ ) + λ ) , i f   x 0 ( 1 ϵ ) 2 Φ ( λ ) [ Φ ( x σ ( 1 ϵ ) λ ) + Φ ( λ ) 1 ] , i f   x 0
Proof. 
It is immediate from the last proof. □
Proposition 3.
Let there be X b t p n ( σ , λ , ϵ ) . Its quantile function is given by
Q X ( p ) = F X 1 ( p ) = { σ ( 1 + ϵ ) [ Φ 1 ( 2 p Φ ( λ ) 1 + ϵ ) λ ] , i f   0 < p 1 + ϵ 2 σ ( 1 ϵ ) [ Φ 1 ( 2 p Φ ( λ ) 1 ϵ Φ ( λ ) + 1 ) + λ ] , i f   1 + ϵ 2 < p < 1
Proof. 
It is immediate from inverting the cdf for the btpn distribution given in Proposition 2. □
Corollary 1.
The median for X b t p n ( σ , λ , ϵ ) is given by
M e ( X ) = { σ ( 1 + ϵ ) [ Φ 1 ( Φ ( λ ) 1 + ϵ ) λ ] , i f   ϵ 0 σ ( 1 ϵ ) [ Φ 1 ( ϵ Φ ( λ ) ( 1 ϵ ) + 1 ) + λ ] , i f   ϵ < 0 .
Corollary 2.
The median for X b t p n ( σ , λ , ϵ ) is < 0 , = 0 and > 0 if ϵ is > 0 , = 0 and < 0 , respectively.

2.2. Moments and Moment-Generating Function

The following proposition presents the central moments of the btpn distribution.
Proposition 4.
Let X b t p n ( σ , λ , ϵ ) . The r-th central moment of X is given by
E ( X r ) = σ r 2 2 π Φ ( λ ) [ ( 1 ) r ( 1 + ϵ ) r + 1 + ( 1 ϵ ) r + 1 ] k = 0 r ( r k ) λ r k ( 2 ) ( k 1 ) / 2 Γ ( k + 1 ) / 2 , λ 2 / 2 ) ,
where Γ ( a , b ) = b + t a 1 e t d t is the upper incomplete gamma function.
Proof. 
Note that E ( X r ) = E 1 ( X r ) + E 2 ( X r ) , where E 1 ( X r ) = 0 x r f ( x ) d x and E 2 ( X r ) = 0 + x r f ( x ) d x . For the first term, we perform the change of variable u = x σ ( 1 + ϵ ) λ . With this,
E 1 ( X r ) = 0 x r 2 σ ( 1 + ϵ ) ϕ ( x σ ( 1 + ϵ ) λ ) d x = ( σ ) r ( 1 + ϵ ) ( r + 1 ) 2 Φ ( λ ) λ ( u + λ ) r ϕ ( u ) d u .
Using the binomial theorem and the change of variable t = u 2 / 2 in the last expression, we obtain
E 1 ( X r ) = ( σ ) r ( 1 + ϵ ) ( r + 1 ) 2 2 π Φ ( λ ) k = 0 r ( r k ) λ r k λ u k e u 2 / 2 d u , = ( σ ) r ( 1 + ϵ ) ( r + 1 ) 2 2 π Φ ( λ ) k = 0 r ( r k ) λ r k 2 ( k 1 ) / 2 λ 2 / 2 t ( k 1 ) / 2 e t d t ,
Note that the last integral corresponds to Γ ( ( k + 1 ) / 2 , λ 2 / 2 ) .
On the other hand, for E 2 ( X r ) , we perform the change of variable u = x σ ( 1 ϵ ) λ , obtaining
E 2 ( X r ) = 0 x r 2 σ ( 1 ϵ ) ϕ ( x σ ( 1 ϵ ) λ ) d x = σ r ( 1 ϵ ) ( r + 1 ) 2 Φ ( λ ) λ ( u + λ ) r ϕ ( u ) d u .
Using the same routine calculation, we obtain
E 2 ( X r ) = σ r ( 1 ϵ ) ( r + 1 ) 2 2 π Φ ( λ ) k = 0 r ( r k ) λ r k 2 ( k 1 ) / 2 λ 2 / 2 t ( k 1 ) / 2 e t d t ,
where again the last integral corresponds to Γ ( ( k + 1 ) / 2 , λ 2 / 2 ) . The final result is obtained by summing E 1 ( X r ) and E 2 ( X r ) . □
Proposition 5.
Let X b t p n ( σ , λ , ϵ ) . The moment-generating function (mgf) for X is given by
M x ( t ) = ( 1 + ϵ ) exp { ( 1 + ϵ ) 2 2 t 2 σ 2 t σ λ ( 1 + ϵ ) } Φ ( λ t σ ( 1 + ϵ ) ) 2 Φ ( λ ) , + ( 1 ϵ ) exp { ( 1 ϵ ) 2 2 t 2 σ 2 + t σ λ ( 1 ϵ ) } Φ ( λ + t σ ( 1 ϵ ) ) 2 Φ ( λ ) .
Proof. 
Note that M X ( t ) = M X 1 ( t ) + M X 2 ( t ) , where M X 1 ( t ) = 0 e t x f ( x ) d x and M X 2 ( t ) = 0 + e t x f ( x ) d x . For the first integral and using the change of variable u = x / [ σ ( 1 + ϵ ) ] + λ , we obtain
M X 1 ( t ) = ( 1 + ϵ ) 2 Φ ( λ ) λ e t σ ( 1 + ϵ ) ( μ λ ) ϕ ( u ) d u .
Completing the square of a binomial in the last term of the exponential and using the change of variable z = μ t σ ( 1 + ϵ ) , we have
M X 1 ( t ) = ( 1 + ϵ ) exp { t 2 σ 2 ( 1 + ϵ ) 2 2 t σ λ ( 1 + ϵ ) } 2 Φ ( λ ) λ t σ ( 1 + ϵ ) e z 2 / 2 2 π d z , = ( 1 + ϵ ) exp { t 2 σ 2 ( 1 + ϵ ) 2 2 t σ λ ( 1 + ϵ ) } 2 Φ ( λ ) Φ λ t σ ( 1 + ϵ ) .
For X 0 , and similarly to the previous development, we use the change of variable u = x / [ σ ( 1 ϵ ) ] λ , obtaining that
M X 2 ( t ) = ( 1 ϵ ) 2 Φ ( λ ) λ e t σ ( 1 ϵ ) ( μ + λ ) ϕ ( u ) d u .
Again, completing the square and using z = μ t σ ( 1 ϵ ) , we obtain
M X 2 ( t ) = ( 1 ϵ ) exp { t 2 σ 2 ( 1 ϵ ) 2 2 + t σ λ ( 1 ϵ ) } 2 Φ ( λ ) λ t σ ( 1 ϵ ) e z 2 / 2 2 π d z , = ( 1 ϵ ) exp { t 2 σ 2 ( 1 ϵ ) 2 2 + t σ λ ( 1 ϵ ) } 2 Φ ( λ ) Φ ( λ + t σ ( 1 ϵ ) ) .
Finally, the result is obtained by summing M X 1 ( t ) and M X 2 ( t ) . □
Corollary 3.
Using properties of the mgf, the first four moments of X b t p n ( σ , λ , ϵ ) can be obtained from the expression μ r = ( r M X ( t ) / t r ) | t = 0 .
  • μ 1 = E [ X 1 ] = 2 ϵ σ [ λ + Ω ( λ ) ]
  • μ 2 = E [ X 2 ] = σ 2 ( 1 + 3 ϵ 2 ) [ λ 2 + λ Ω ( λ ) + 1 ]
  • μ 3 = E [ X 3 ] = 4 ϵ σ 3 ( 1 + ϵ 2 ) [ λ 3 + λ 2 Ω ( λ ) + λ + 2 Ω ( λ ) ]
  • μ 4 = E [ X 4 ] = σ 4 [ 1 + 5 ϵ 2 ( 2 + ϵ ) ] [ λ 4 + λ 3 Ω ( λ ) + 6 λ 2 + 5 λ Ω ( λ ) + 3 ]
where Ω ( λ ) = ϕ ( λ ) / Φ ( λ ) is the reciprocal of the Mill’s ratio for the standard normal distribution.
Corollary 4.
The variance, coefficients of skewness and kurtosis for X b t p n ( σ , λ , ϵ ) are given by
V a r ( X ) = σ 2 [ λ 2 ( 1 ϵ 2 ) + λ Ω ( λ ) ( 1 5 ϵ 2 ) + ϵ 2 ( 3 4 Ω ( λ ) ) + 1 ] , b 1 = 4 ϵ ( 1 + ϵ 2 ) [ λ 3 + λ 2 Ω ( λ ) + 2 Ω ( λ ) ] [ ( 1 + 3 ϵ ) ( λ 2 + λ Ω ( λ ) + 1 ) ] 3 / 2 , a n d b 2 = [ 1 + 5 ϵ 2 ( 2 + ϵ ) ] [ λ 4 + λ 3 Ω ( λ ) + 6 λ 2 + 5 λ Ω ( λ ) + 3 ] [ ( 1 + 3 ϵ 2 ) ( λ 2 + λ Ω ( λ ) + 1 ) ] 2 ,
respectively.
Figure 2 shows the plots for asymmetry and kurtosis coefficients. Note that a more right-skewed distribution is obtained when ϵ 1 and λ , whereas a more left-skewed model is obtained when ϵ 1 and λ . On the other hand, a greater kurtosis is obtained when ϵ 1 and λ , whereas a lower kurtosis is obtained when | ϵ | 1 and λ . Note that this pattern is consistent with the pdf for different parameters presented in Figure 1.

2.3. Mode and Unimodality and Bimodality Regions

The next proposition presents the unimodality and bimodality property of the btpn distribution.
Proposition 6.
Let X b t p n ( σ , λ , ϵ ) . For λ 0 , the model is unimodal, and for λ > 0 , the model is bimodal. Moreover, for the unimodal case, the mode of the model is 0, and for the bimodal case, the two modes are σ λ ( 1 + ϵ ) and σ λ ( 1 ϵ ) , respectively.
Proof. 
By definition, the mode is the value that maximizes the pdf or, equivalently, the logarithm of the pdf. For X b t p n ( σ , λ , ϵ ) , it is straighforward to show that
log f ( x ; σ , λ , ϵ ) x = { 1 σ ( 1 + ϵ ) ( x σ ( 1 + ϵ ) + λ ) , if   x < 0 1 σ ( 1 ϵ ) ( x σ ( 1 ϵ ) λ ) , if   x 0           and 2 log f ( x ; σ , λ , ϵ ) x 2 = { 1 σ 2 ( 1 + ϵ ) 2 , if   x < 0 1 σ 2 ( 1 ϵ ) 2 , if   x 0 .
Therefore, solving the equation log f ( x ; σ , λ , ϵ ) / x = 0 , we obtain x 1 = σ λ ( 1 + ϵ ) and x 2 = σ λ ( 1 ϵ ) as the potential mode for each branch, because the second derivative is negative for each respective case. However, this is valid if and only if x 1 < 0 and x 2 > 0 , respectively. In other words, if σ λ ( 1 + ϵ ) < 0 , then x 1 is a mode and if σ λ ( 1 ϵ ) > 0 , then x 2 also is a mode. This is equivalent to
{ λ > 0 1 + ϵ > 0 } { λ < 0 1 + ϵ < 0 } x 1   is   a   mode   and { λ > 0 1 ϵ > 0 } { λ < 0 1 ϵ < 0 } x 2   is   a   mode .
where 1 + ϵ > 0 and 1 ϵ > 0 , ϵ ( 1 , 1 ) . On the other hand, 1 + ϵ 0 and 1 ϵ 0 . For this reason, it is immediate that for λ > 0 , the btpn distribution have two modes, and such modes are x 1 and x 2 . Finally, for λ 0 , it is immediate that log f ( x ; σ , λ , ϵ ) / x > 0 , for x < 0 and log f ( x ; σ , λ , ϵ ) / x 0 for x 0 . In other words, the pdf for the btpn distribution is an increasing function in ( , 0 ) and a decreasing function in ( 0 , ) , where we can deduce that the model is unimodal and the respective mode is attached in zero. □
Figure 3 shows the regions of unimodality and bimodality for the btpn depending on the parameters λ and ϵ .

2.4. Particular Cases

By construction, the following models are particular cases for the btpn distribution:
  • btpn ( σ = σ / 2 , λ , ϵ = 1 ) tpn ( σ , λ ) ;
  • btpn ( σ , λ = 0 , ϵ = 0 ) N ( 0 , σ 2 ) , i.e., the normal distribution with mean 0 and variance σ 2 ;
  • btpn ( σ = 1 , λ = 0 , ϵ ) esn ( σ , ϵ ) , i.e., the epsilon skew-normal distribution (Mudholkar and Hutson [10]).
Figure 4 summarizes the relationships among the btpn and its particular cases.

3. Inference

In this section, we discuss the maximum likelihood (ML) method for parameter estimation for the btpn model. We also provide details about the computational aspects.

3.1. Maximum Likelihood Function

Hereafter, and to simplify the estimation procedure, we consider the reparameterization η = ϵ 1 ϵ 2 R . Therefore, henceforth, we denote X b t p n ( σ , λ , η ) , with σ , λ and η scale, shape and asymmetry parameters, respectively, if its pdf is given by
f ( x ; σ , λ , η ) = { 1 2 σ Φ ( λ ) ϕ ( x 1 + η 2 σ ( 1 + η 2 + η ) λ ) , if   x < 0 1 2 σ Φ ( λ ) ϕ ( x 1 + η 2 σ ( 1 + η 2 η ) λ ) , if   x 0
Given z 1 , , z n , a random sample from the btpn ( σ , λ , η ) distribution, the log-likelihood function for θ = ( σ , λ , η ) is given by
( θ ) = n [ log ( 2 σ Φ ( λ ) ) + 1 2 log ( 2 π ) + λ 2 2 ] + 1 ( θ ) + 2 ( θ ) ,
where
1 ( θ ) = 1 2 i : z i 0 { ( z i 1 + η 2 σ ( 1 + η 2 η ) ) 2 2 z i λ 1 + η 2 σ ( 1 + η 2 η ) } ,   and 2 ( θ ) = 1 2 i : z i > 0 { ( z i 1 + η 2 σ ( 1 + η 2 + η ) ) 2 + 2 z i λ 1 + η 2 σ ( 1 + η 2 + η ) } .
To find the ML estimator of θ , say θ ^ , we need to maximize ( θ ) in (1) in relation to θ . However, no closed-form expressions for the ML estimates are possible. Therefore, we must use an iterative method for nonlinear optimization. For instance, we solve this problem using the Broyden-–Fletcher–-Goldfarb–-Shanno (BFGS) quasi-Newton method; see [11] (p. 199).

3.2. Computational Aspects

The ML estimators for the btpn model and the obtaining of their standard errors are included in the tpn package [12] from the R [13] software. The following function can be used to obtain these results:
  •     est.btpn(y)
where y is the sample. The function returns a list with the estimates, the iterations used for the maximization algorithm, the log-likelihood function evaluated in the parameter estimations and the corresponding Akaike information criterion (AIC, see [14]) and the Bayesian information criterion (BIC, see [15]). Models with lower AIC and/or BIC are preferable. The package also includes the functions to drawn values to evaluate the pdf and the cdf for the btpn model named rbtpn, dbptn and pbtpn, respectively.

4. Simulation Study

In this section, we present a simulation study in order to evaluate the behaviour of the ML estimators in finite samples. The study was conducted using the tpn package [12]. Specifically, random samples were generated using the rbtpn function, and the estimation was performed using the est.btpn function. We considered 5000 Monte Carlo replicates for 3 sample sizes: 50, 100 and 200. We also considered 2 combinations for the scale parameter σ : 2 and 10; 3 values for λ : 0.75 , 1 and 3; and 2 values for η : 0.5 and 0.75 . This setting provides 36 combinations of the parameters σ , λ and η and the sample size. Table 1 and Table 2 summarize the empirical bias, the standard errors of the MLE (SE), the root-mean-squared error (RMSE) and the 95% coverage probability (CP) based on the asymptotic distribution of the MLE. In general terms, the bias and RMSE terms are reduced when the sample size is increased, suggesting the consistency of the MLE. Note also that the SE and RMSE terms are closer when the sample size is increased, suggesting that the standard errors of the estimators are also well estimated. Additionally, the CP terms converge reasonably to the nominal value used to their construction (95%), suggesting that the normality is reasonable as an asymptotic distribution to the ML estimators in the btpn model, even in reasonable sample sizes.

5. Application

In this section, we present an application to a real data set in order to illustrate the btpn model. We consider the height data set, which consists of the height of 126 students from the University of Pennsylvania (Cruz-Medina [16]). We compare our proposal with other bimodal proposals, such as the epsilon skew inverted gamma (esig, see Abdulah et al. [17]) and the alpha skew-normal (asn, Elal-Olivero [18]). The pdf for the esig model is given by:
f ( x ; σ , λ , η ) = λ σ 2 Γ ( σ ) { ( x 1 η ) ( σ + 1 ) e λ ( 1 η ) x x 0 ( x 1 + η ) ( σ + 1 ) e λ ( 1 + η ) x x < 0 ,
where λ > 0 , σ > 0 and | η | < 1 are the scale, shape and skewness parameters, respectively.
The pdf for the asn model is given by:
f ( x ; η , λ , σ ) = ( [ 1 η ( x λ σ ) ] 2 + 1 σ ( 2 + η 2 ) ) ϕ ( x λ σ ) ,
where η , λ R and σ > 0 are the shape, location and scale parameters, respectively.
Table 3 summarizes some descriptive statistics for the sample, where we highlight the symmetrical behaviour of the data ( b 1 = 0.05 ).
Table 4 presents the estimatives, standard errors, AIC and BIC criteria for the mentioned models. Note that, based on both criteria, btpn presents a better fit than the rest of the distributions. Figure 5 shows the histogram for the height data and the pdf for the three considered distributions, where the better performance for the btpn in this data set is demonstrated. Moreover, as discussed in Proposition 6, λ ^ = 0.496 > 0 implies a bimodal model, and such modes are equal to x 1 = 0.402 and x 2 = 0.404 . In addition, the distribution of height is very close to symmetry.
We also compute the randomized quantile residuals [19] for the three fitted models. If the model was correctly specified, these residuals should be a random sample from the standard normal distribution. Figure 6 shows the qqplot for such residuals, also suggesting that the btpn is a more appropriated model for this data set.

6. Conclusions

The importance of fitting an observable dataset by a probability distribution is well-known, since it will be covered with convenient properties. Difficulty arises when the data is bimodal, because there are not traditional distributions with this property. This gap is being filled by an increasing movement in the statistical literature to develop probability distributions which already have a bimodality feature. In this paper, we made our contribution with the bimodal positive truncation normal distribution. The btpn distribution has the following advantages: support in the real line, closed-form cdf and moments, and the ability to generalize the standard normal and treatable maximum likelihood estimators. The ML procedure works very reasonably, i.e, as the sample size increases, the bias and the SE decrease. Since there are models for which the estimation procedure does not work even for large samples, the btpn distribution also has this strength. We ended the advantages of our proposed distribution with an application where btpn was the best choice of fitting. As suggestions for future work, we can mention two possibilities: the first entails the improvement of the asymptotic properties of the ML estimation through bias and variance corrections (see [20,21], respectively), and the second involves the addition of a regression structure. A closed-form cdf allows even a quantile regression structure, see [22,23], for instance, as [24] did for the gamma–sinh Cauchy distribution. For the applicability and possibilities of future works, we think the bimodal positive truncation normal distribution is useful for practitioners and researchers of many different areas.

Author Contributions

Conceptualization, H.J.G., Y.M.G. and D.I.G.; Data curation, M.C.; Formal analysis, W.E.C., Y.M.G. and M.C.; Investigation, W.E.C. and T.M.M.; Methodology, H.J.G., T.M.M. and D.I.G.; Software, H.J.G., Y.M.G. and D.I.G.; Supervision, Y.M.G. and D.I.G.; Validation, M.C.; Visualization, W.E.C.; Writing—review editing, T.M.M. and D.I.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research of H.J.G. was supported by the Proyecto de Investigación de Facultad de Ingeniería. Universidad Católica de Temuco. UCT-FDI032020. W.E.C., Y.M.G., M.C. and D.I.G. also acknowledge the support of Proyecto Gidi: “La estadística como respuesta a problemas de otras áreas” supported by the University of Atacama.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ashman, K.M.; Bird, C.M.; Zepf, S.E. Detecting Bimodality in Astronomical Datasets. Astron. J. 1994, 108, 2348–2361. [Google Scholar] [CrossRef] [Green Version]
  2. Michele, C.; Accatino, F. Tree cover bimodality in savannas and forests emerging from the switching between two fire dynamics. PLoS ONE 2014, 9, e91195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, J.; Wen, S.; Symmans, W.F.; Pusztai, L.; Coombes, K.R. The bimodality index: A criterion for discovering and ranking bimodal signatures from cancer gene expression profiling data. Cancer Inform. 2009, 7, 199–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Robertson, C.A.; Fryer, J.G. Some descriptive properties of normal mixtures. Scand. Actuar. J. 1969, 3–4, 137–146. [Google Scholar] [CrossRef]
  5. Gómez, H.W.; Elal-Olivero, D.; Salinas, H.S.; Bolfarine, H. Bimodal extension based on the skew-normal distribution with application to pollen data. Environmetrics 2011, 22, 50–62. [Google Scholar] [CrossRef]
  6. Venegas, O.; Salinas, H.S.; Gallardo, D.I.; Bolfarine, H.; Gómez, H.W. Bimodality based on the generalized skew-normal distribution. J. Stat. Comput. Simul. 2018, 88, 156–181. [Google Scholar] [CrossRef]
  7. Butt, N.S.; Khalil, M.G. A New Bimodal Distribution for Modeling Asymmetric Bimodal Heavy-Tail Real Lifetime Data. Symmetry 2020, 12, 2058. [Google Scholar] [CrossRef]
  8. Gómez, Y.M.; Gómez-Déniz, E.; Venegas, O.; Gallardo, D.I.; Gómez, H.W. An Asymmetric Bimodal Distribution with Application to Quantile Regression. Symmetry 2019, 11, 899. [Google Scholar] [CrossRef] [Green Version]
  9. Gómez, H.J.; Olmos, N.M.; Varela, H.; Bolfarine, H. Inference for a truncated positive normal distribution. Appl. Math. J. Chin. Univ. 2018, 33, 163–176. [Google Scholar] [CrossRef]
  10. Mudholkar, G.S.; Hutson, A.D. The epsilon–skew–normal distribution for analyzing near-normal data. J. Stat. Plan. Inference 2000, 83, 291–309. [Google Scholar] [CrossRef]
  11. Mittelhammer, R.C.; Judge, G.G.; Miller, D.J. Econometric Foundations; Cambridge University Press: New York, NY, USA, 2000. [Google Scholar]
  12. Gallardo, D.I.; Gómez, H.J.; Gómez, Y.M. tpn: Truncated Positive Normal Model and Extensions. R Package Version 1.1. 2021. Available online: https://cran.r-project.org/web/packages/tpn/index.html (accessed on 25 January 2022).
  13. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 25 January 2022).
  14. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  15. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  16. Cruz-Medina, I.R.; Olmos, N.M. Almost Nonparametric and Nonparametric Estimation in Mixture Model. Ph.D. Thesis, Pennsylvania State University, State College, PA, USA, 2001. [Google Scholar]
  17. Abdulah, E.K.; Elsalloukh, H. Bimodal Class based on the Inverted Symmetrized Gamma Distribution with Applications. J. Stat. Appl. Probab. 2014, 3, 1–7. [Google Scholar] [CrossRef]
  18. Elal-Olivero, D. Alpha-skew-normal distribution. Proyecciones 2010, 29, 224–240. [Google Scholar] [CrossRef] [Green Version]
  19. Dunn, P.K.; Smyth, G.K. Randomized quantile residuals. J. Comput. Graph. Stat. 1996, 5, 236–244. [Google Scholar]
  20. Magalhães, T.M.; Gómez, Y.M.; Gallardo, D.I.; Venegas, O. Bias reduction for the Marshall-Olkin extended family of distributions with application to an airplane’s air conditioning system and precipitation data. Symmetry 2020, 12, 851. [Google Scholar] [CrossRef]
  21. Magalhães, T.M.; Botter, D.A.; Sandoval, M.C. A general expression for second-order covariance matrices—An application to dispersion models. Braz. J. Probab. Stat. 2021, 35, 37–49. [Google Scholar] [CrossRef]
  22. Cade, B.S.; Noon, B.R. A gentle introduction to quantile regression for ecologists. Front. Ecol. Environ. 2003, 1, 412–420. [Google Scholar] [CrossRef]
  23. Alencar, A.P.; Santos, B.R. Association of pollution with quantiles and expectations of the hospitalization rate of elderly people by respiratory diseases in the city of São Paulo, Brazil. Environmetrics 2014, 25, 165–171. [Google Scholar] [CrossRef]
  24. Gómez, Y.M.; Gallardo, D.I.; Venegas, O.; Magalhães, T.M. An asymmetric bimodal double regression model. Symmetry 2021, 13, 2279. [Google Scholar] [CrossRef]
Figure 1. Pdf for btpn ( σ = 1 , λ , ϵ ) with different fixed values for λ and varying ϵ : (a) λ = 2 ; (b) λ = 0.2 ; (c) λ = 0.4 and; (d) λ = 1.5 .
Figure 1. Pdf for btpn ( σ = 1 , λ , ϵ ) with different fixed values for λ and varying ϵ : (a) λ = 2 ; (b) λ = 0.2 ; (c) λ = 0.4 and; (d) λ = 1.5 .
Symmetry 14 00665 g001
Figure 2. (a) Asymmetry coefficient and (b) kurtosis coefficient for btpn ( σ = 1 , λ , ϵ ) distribution.
Figure 2. (a) Asymmetry coefficient and (b) kurtosis coefficient for btpn ( σ = 1 , λ , ϵ ) distribution.
Symmetry 14 00665 g002
Figure 3. Regions of unimodality and bimodality for the btpn model in terms of λ and η .
Figure 3. Regions of unimodality and bimodality for the btpn model in terms of λ and η .
Symmetry 14 00665 g003
Figure 4. Particular cases for the btpn distribution.
Figure 4. Particular cases for the btpn distribution.
Symmetry 14 00665 g004
Figure 5. Histogram for the h e i g h t data set and the estimated pdf for the btpn, esig and asn models.
Figure 5. Histogram for the h e i g h t data set and the estimated pdf for the btpn, esig and asn models.
Symmetry 14 00665 g005
Figure 6. Quantile residuals for fitted models: (a) asn, (b) esig, and (c) btpn.
Figure 6. Quantile residuals for fitted models: (a) asn, (b) esig, and (c) btpn.
Symmetry 14 00665 g006
Table 1. Empirical bias, SE, RMSE and 95% CP for the ML estimators of σ , λ and η in the btpn distribution with different combinations of parameters (case true σ = 2 ).
Table 1. Empirical bias, SE, RMSE and 95% CP for the ML estimators of σ , λ and η in the btpn distribution with different combinations of parameters (case true σ = 2 ).
True Value n = 50 n = 100 n = 200
λ η par.biasSERMSECPbiasSERMSECPbiasSERMSECP
−0.75−0.5 σ 0.8943.2718.4030.8140.1701.0893.1650.8560.0350.5140.5840.898
λ −0.8313.4058.5610.862−0.1321.2373.0500.891−0.0280.6380.6910.917
η −0.0180.1120.1210.940−0.0080.0780.0810.945−0.0030.0550.0570.947
0.75 σ 0.7362.4907.3810.8020.2331.2183.4970.8620.0410.5380.6300.901
λ −0.6632.6637.3250.855−0.2051.3683.4530.897−0.0290.6600.7280.924
η 0.0300.1420.1670.9360.0140.0990.1060.9450.0090.0700.0730.947
1−0.5 σ −0.0570.3500.3600.887−0.0260.2480.2500.914−0.0120.1750.1760.936
λ 0.0730.4060.4160.9380.0350.2850.2890.9460.0140.2010.2010.949
η −0.0140.0860.0950.932−0.0080.0610.0640.940−0.0030.0430.0450.944
0.75 σ −0.0550.3500.3550.884−0.0240.2480.2510.919−0.0150.1740.1760.931
λ 0.0720.4060.4120.9350.0350.2850.2880.9430.0200.2000.2030.942
η 0.0280.1090.1270.9400.0120.0760.0820.9420.0060.0540.0550.943
3−0.5 σ −0.0490.2030.2140.919−0.0220.1450.1520.930−0.0130.1030.1060.936
λ 0.1030.3560.3830.9480.0470.2480.2620.9470.0280.1740.1800.948
η −0.0070.0510.0540.937−0.0030.0360.0360.949−0.0010.0250.0260.950
0.75 σ −0.0480.2030.2140.917−0.0230.1450.1480.932−0.0120.1030.1030.946
λ 0.1050.3560.3850.9460.0460.2480.2530.9520.0240.1740.1750.953
η 0.0110.0640.0720.9330.0050.0450.0480.9390.0020.0320.0330.942
Table 2. Empirical bias, SE, RMSE and 95% CP for the ML estimators of σ , λ and η in the btpn distribution with different combinations of parameters (case true σ = 10 ).
Table 2. Empirical bias, SE, RMSE and 95% CP for the ML estimators of σ , λ and η in the btpn distribution with different combinations of parameters (case true σ = 10 ).
True Value n = 50 n = 100 n = 200
λ η par.biasSERMSECPbiasSERMSECPbiasSERMSECP
−0.75−0.5 σ 2.27910.59719.5710.8020.9615.25411.6950.8640.2212.5383.4110.900
λ −0.4002.3364.0230.857−0.1611.2092.3280.896−0.0340.6320.7800.924
η −0.0200.1120.1260.938−0.0090.0790.0830.949−0.0040.0550.0560.950
0.75 σ 1.7489.14717.6940.7970.6085.0637.3120.8590.1572.4963.2250.891
λ −0.2722.0463.4980.856−0.0931.1721.5480.895−0.0210.6230.7340.920
η 0.0380.1430.1770.9330.0130.0990.1030.9510.0070.0690.0720.945
1−0.5 σ −0.2251.7721.8290.887−0.1261.2371.2600.916−0.0570.8760.8910.932
λ 0.0650.4080.4150.9380.0370.2850.2890.9400.0150.2010.2050.943
η −0.0140.0860.0950.936−0.0050.0600.0630.943−0.0020.0430.0430.952
0.75 σ −0.2841.7441.7560.886−0.1361.2361.2550.913−0.0590.8750.8690.933
λ 0.0730.4050.4050.9410.0370.2850.2880.9430.0160.2010.2010.946
η 0.0280.1090.1310.9260.0120.0760.0810.9450.0080.0540.0560.946
3−0.5 σ 1.9852.78430.1650.9131.6991.82426.9820.9281.1130.79922.8870.932
λ 0.0230.4201.3200.941−0.0320.2871.1780.947−0.0250.1840.9550.941
η −0.0040.0510.0540.934−0.0020.0360.0370.944−0.0010.0260.0260.948
0.75 σ 0.3864.61813.3990.9180.5961.81614.1210.9240.3550.92910.5610.943
λ 0.0710.4740.6920.9540.0160.2850.6800.9410.0030.1880.5030.952
η 0.0190.0650.0760.9250.0150.0460.0510.9350.0080.0330.0340.941
Table 3. Descriptive statistics for the h e i g h t data set.
Table 3. Descriptive statistics for the h e i g h t data set.
Data Setn X ¯ S 2 b 1 b 2
weight measured12601 0.05 3.053
Table 4. Estimated parameters and their standard errors (in parentheses) for the btpn, esig and asn models for the h e i g h t data set. The AIC and BIC criteria are also presented.
Table 4. Estimated parameters and their standard errors (in parentheses) for the btpn, esig and asn models for the h e i g h t data set. The AIC and BIC criteria are also presented.
Estimatedbtpnesigasn
σ 0.813 (0.113)1.304 (0.148)0.996 (0.063)
λ 0.496 (0.316)0.527 (0.073)0.014 (3.422)
η −0.002 (0.048)0.095 (0.059)0.014 (3.409)
AIC360.76415.79362.67
BIC369.27424.30375.91
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gómez, H.J.; Caimanque, W.E.; Gómez, Y.M.; Magalhães, T.M.; Concha, M.; Gallardo, D.I. A Bimodal Model Based on Truncation Positive Normal with Application to Height Data. Symmetry 2022, 14, 665. https://doi.org/10.3390/sym14040665

AMA Style

Gómez HJ, Caimanque WE, Gómez YM, Magalhães TM, Concha M, Gallardo DI. A Bimodal Model Based on Truncation Positive Normal with Application to Height Data. Symmetry. 2022; 14(4):665. https://doi.org/10.3390/sym14040665

Chicago/Turabian Style

Gómez, Héctor J., Wilson E. Caimanque, Yolanda M. Gómez, Tiago M. Magalhães, Miguel Concha, and Diego I. Gallardo. 2022. "A Bimodal Model Based on Truncation Positive Normal with Application to Height Data" Symmetry 14, no. 4: 665. https://doi.org/10.3390/sym14040665

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