A Unimodal/Bimodal Skew/Symmetric Distribution Generated from Lambert’s Transformation

: The generalized bimodal distribution is especially efﬁcient in modeling univariate data exhibiting symmetry and bimodality. However, its performance is poor when the data show important levels of skewness. This article introduces a new unimodal/bimodal distribution capable of modeling different skewness levels. The proposal arises from the recently introduced Lambert transformation when considering a generalized bimodal baseline distribution. The bimodal-normal and generalized bimodal distributions can be derived as special cases of the new distribution. The main structural properties are derived and the parameter estimation is carried out under the maximum likelihood method. The behavior of the estimators is assessed through simulation experiments. Finally, two applications are presented in order to illustrate the utility of the proposed distribution in data modeling in different real settings.


Introduction
Analysts often have to deal with data that exhibit bimodality; for example, when observing the size of worker ants in weaver ant colonies [1], the duration of volcanic eruptions [2], the amount of excretion of mercury in urine [3], the grain size in sintered Zirconia [4], or the amount of tropospheric water vapor in the tropics [5].
Two-component mixture distributions are usually used to model data that exhibit bimodality. These distributions have a very flexible density function (unimodal/bimodal), a highly valued feature when trying to model data in different real settings. However, one difficulty in working with them is that it is necessary to deal with the problem of the non-identifiability of their parameters, see McLachlan et al. [6] for more detail.
It is possible to find in the literature guidelines to deal with the problem of nonidentifiability in a mixture distribution; for example, as suggested by Aitkin and Rubin [7], it is convenient to assume a certain restriction on the mixing parameter to evaluate the behavior of the maximum likelihood estimates. An alternative that has been considered in various studies is to impose restrictions on the components of the mixture distribution; for example, assuming that the components have the same variance. In this way, subfamilies of mixture distributions with a simpler parametric structure (where the parameters are identifiable) are defined.
Taking into account certain restrictions for the components of a mixture distribution implies working with a sub-family that exhibits a less flexible density function than the unrestricted case. However, it is possible to find in the literature sub-families of twocomponent mixture distributions that have been shown to be useful in data modeling in various real settings. In this context, it is possible to find the generalized bimodal (GB) distribution originally proposed by Rao [8] and later studied as a special case of the bimodal symmetric distribution proposed in Sarma et al. [9].
A random variable X has a generalized bimodal distribution, denoted as X ∼ GB(γ), if its probability density function (pdf) is given by and its cumulative distribution function (cdf) is given by where γ is a shape parameter that controls bimodality and φ(·) and Φ(·) denote the pdf and the cdf of the standard normal distribution, respectively. It is easy to verify that Equation (1) corresponds to the pdf of a two-component mixture distribution with mixing parameter 1/(1 + γ), for γ ∈ [0, 2), one component having a standard normal distribution, and the other a standard bimodal-normal (BN) distribution. If γ = 0, the Lambert generalized bimodal (LGB) distribution reduces to the BN distribution. Details on the structural properties of the BN distribution can be found in Hassan and Hijazi [10] and Elal-Olivero [11].
The parametric space bounded to the interval [0,2) for γ can be explained by the fact that the pdf given in Equation (1) is bimodal for any value of γ in such interval. However, we observe that the range of γ can be extended so that γ assumes values in the interval [0, ∞). In such case, it is possible to verify the following properties: The GB pdf is symmetric-bimodal when γ ∈ [0, 2) and symmetric-unimodal when γ ∈ [2, ∞). The GB distribution reduces to the bimodal-normal distribution when γ = 0 and tends to the standard normal distribution as γ → ∞.
The symmetry characteristic of the GB pdf can be considered a desirable characteristic in the analysis of certain data sets, but a limitation in the analysis of others. In the literature, it is possible to find different construction methods that allow generating an asymmetric distribution from a symmetric baseline distribution, see Azzalini [12], Eugene et al. [13], Cordeiro and de Castro [14], Ferreira and Steel [15], Goerg [16], and Alzaatreh et al. [17], among others.
Recently, Iriarte et al. [18] introduce the Lambert-F distribution generator defined as where α is an extra shape parameter and F(·) is the cdf of an arbitrary baseline distribution. The transformation given in Equation (3) defines a new family of distributions more flexible in terms of skewness than the baseline distribution. Iriarte et al. [18] study two special cases of Equation (3), extending the classical exponential and Rayleigh distributions and showing that the hazard rate functions induced by this transformation can be understood as modifications in the early times of the baseline hazard rate functions.
In this article, we introduce a new unimodal/bimodal distribution that generalizes the GB distribution and is capable of modeling different levels of skewness. The proposal, called the Lambert generalized bimodal (LGB) distribution, arises from Equation (3) when considering a baseline GB distribution. The result is a new distribution that generalizes to the GB and BN distributions and that can serve as an alternative to other asymmetric bimodal distributions in the literature.
The article is organized as follows. In Section 2, we define the LGB random variable and derive the density and distribution functions. In Section 3, we describe the characteristics of unimodality, bimodality, asymmetry, and kurtosis. In addition, alias distributions are discussed. In Section 4, we consider the problem of parameter estimation using the maximum likelihood (ML) method. In Section 5, the behavior of the ML estimators and the utility of the proposed distribution are evaluated through simulation experiments. Section 6 presents two application examples illustrating the usefulness of the LGB distribution in real settings. Finally, the main conclusions are reported in Section 7.

The LGB Distribution
In this section, we introduce the LGB distribution and study some of its main properties.

Remark 1.
The quantile function of the GB distribution does not have a closed analytical form. However, it can be calculated numerically from the cdf given in Equation (2).
Then, the cdf of X is given by Proof. From the Equation (4), for α = 1, we have that is the inverse function of Q GB (·; γ), that is, the standard GB cdf. Then, by definition of the Lambert W function, it follows that and the result is obtained considering that P(U ≤ u) = u, as U ∼ uniform(0, 1). Finally, note that the analytical expression obtained for the cdf of X is valid for α = 1 once . The pdf of X can be obtained in a straightforward way from the cdf given in Equation (5). Corollary 1. Let X ∼ LGB(µ, σ, η, α). Then, the pdf of X is given by In accordance with Definition 1, note that Equations (5) and (6) reduce, respectively, to the cdf and pdf of the GB distribution when α = 1. Consequently, the LGB distribution can be understood as an extension with one extra parameter of the GB distribution. An interesting property of the LGB distribution is that its pdf corresponds to a modification in a multiplicative fashion of the GB pdf. If α ∈ (0, 1) or α ∈ (1, e), the GB pdf is modified in a multiplicative fashion by the expression α }, allowing asymmetric shapes for the LGB pdf. Figure 1 shows some pdf curves for the LGB distribution considering different values of γ and α. Here, it can be seen that the LGB pdf is bimodal symmetric when α = 1. Note that the parameter α has an effect on the shape of the LGB pdf allowing unimodal or bimodal asymmetric shapes. This will be discussed in more detail in Section 3.  The quantile function (qf) of the LGB distribution can be easily derived by inverting Equation (5), considering steps very similar to those of the proof of Proposition 1. The resulting analytical expression for this function is given by As the Lambert W function is included in different statistical software, Equation (7) can be easily computed. A code in the R programming language [19] is provided in Appendix A.

Related Distributions
In the previous section, it was shown that the LGB and GB distributions are nested distributions, with the GB distribution being a special case of the LGB distribution. Now, taking into account that the BN distribution is a special case of the GB distribution, we observe that the LGB distribution reduces to the BN distribution when γ = 0 and α = 1. Consequently, a new special case can be highlighted. If γ = 0, the pdf given in Equation (6) reduces to The pdf given in Equation (8) is bimodal symmetric or bimodal asymmetric depending on the value of α. If α = 1, then Equation (8) reduces to the BN pdf. We refer to Equation (8) as the Lambert-bimodal normal (LBN) pdf. It is important to note that, unlike the LGB distribution, the pdf given in Equation (8) valued in the antimode is equal to 0 for all µ ∈ R, σ > 0, and α ∈ (0, e).
We observe that the parameter α induced by the Lambert-F transformation plays a determining role in the skewness of the LBN distribution. We also see that the LBN pdf can be understood as a modification in a multiplicative fashion of the BN pdf.

Shapes and Aliases
In this section, we discuss the forms of LGB distribution and analyze the possible existence of alias distributions for members of this LGB family.

Shapes
Let f (s) X (x) = ∂ s f X (x; µ, σ, γ, α)/∂x s be the s-th derivative of the LGB pdf with respect to x, s = 1, 2, we have Thus, if X ∼ LGB(µ, σ, γ, α), and if (x * , f X (x * ; γ, α)) and (x * * , f X (x * * ; γ, α)) are, respectively, critical points and inflection points of the pdf of X, then x * is a root of the equation and x * * a root of the equation In the case α = 1, Equations (9) and (10) lead to establish that the LGB pdf is bimodal with modes given by µ ± σ √ 2 − γ, antimode given by µ, and abscissa of inflection points given by µ ± σ √ w j , with j = 1, 2, where w 1 = [(5 − γ) + (5 − γ) 2 − 4(2 − γ)]/2 and In the case where α = 1, it is not possible to obtain closed expressions for the modes, antimode, and abscissa of inflection points. Therefore, these values must be obtained by solving Equations (9) and (10) by numerical procedures. Considering Equations (9) and (10) as functions of x, say w(x) and v(x), respectively, we observe that which means that Equation (9) has at least one root, associated with the unimodal case. Thus, the above also implies that the Equation (10) has at least two roots and, taking into account that it follows that Equation (10) has at least one negative minimum value at x = x 0 , where x 0 satisfies the equation Figure 2 shows a good summary of the shapes that the LGB pdf can take. In this figure, profiles are presented for the equations of critical points and inflection points, Equations (9) and (10), together with the corresponding LGB pdf. As in Figures 1 and 2, it is observed that the LGB pdf can be unimodal or bimodal depending on the values of γ and α, but here it is also observed that a unimodal LGB pdf can have two or four inflection points. In this sense, Equation (9) has one or three roots (depending on whether the pdf is unimodal or bimodal) and Equation (10) has two or four roots (four roots when the pdf is bimodal and two or four roots when the pdf is unimodal). In Figure 3, the regions of unimodality and bimodality established in the plane defined by the ranges of γ and α are presented. This figure was drawn by solving Equation (9) using the uniroot.all function [20] in the R language [19].

Skewness and Kurtosis
Next, the skewness behavior of the LGB distribution is described. The behavior of the kurtosis is also described in the case in which the LGB distribution is unimodal. In this case, first the raw moments of the LGB distribution are derived and later the Fisher's skewness and kurtosis coefficients are analyzed. Proposition 2. Let Z ∼ LGB(0, 1, γ, α) and X ∼ LGB(µ, σ, γ, α). Then, for r = 1, 2, . . ., the r-th raw moment of Z and X are given by E(Z r ) = a r (γ, α) and where Proof. Under the change of variable u = Φ(z) − z γ+1 φ(z), the r-th moment of Z is given by which is the result given in Equation (12). Now, as X can be represented as

Corollary 2.
The mean and the variance of X are given by E( Corollary 3. The skewness (β 1 (γ, α)) and kurtosis (β 2 (γ, α)) coefficients of X are given by Note that the r-th raw moment of the LGB distribution should be computed using numerical integration. Table 1 presents some values for the first four raw moments of the LGB distribution. Figure 4 (left panel) presents some curves for the skewness coefficient of the LGB distribution considering different values of γ and α. In the right panel of the same figure, some curves for the kurtosis coefficient are presented when the LGB pdf is unimodal.
In Table 1, we observe that the odd raw moments are 0 when α = 1, an expected result because in this case the LGB distribution is symmetric around µ. In the Figure 4, we observe that 1.
The LGB pdf is symmetric when α = 1, regardless of the value assumed by the parameter γ.

2.
The LGB distribution can be skewed (positively or negatively) depending on the value assumed by α. If α ∈ (0, 1) or α ∈ (1, e), then the LGB pdf is skewed and the skewness is also controlled by the parameter γ. However, the effect of γ on skewness is important when it assumes small values.

3.
If the LGB pdf is unimodal, then it is asymmetric. 4.
In the unimodal case, the excess kurtosis β 2 (γ, α) − 3 is less than 0; that is, the LGB distribution is a platykurtic distribution.    In the case where E(X) > Median(X), we observe from Equations (7) and (11) that the skewness of the LGB distribution is positive if γ and α satisfy the inequality

Parameters Moments
On the other hand, values of γ and α (α = 1) that do not satisfy the previous inequality are related to a negative skewness.
Based on these last results, together with the results of Section 3.1, we observe that the LGB distribution can be considered as a feasible alternative for modeling data exhibiting asymmetry and uni/bimodality. In this way, the LGB distribution can be used in the analysis of random phenomena raised in various areas of knowledge. Below, we describe some situations in which the LGB distribution could be used: (1) In entomology, volcanology, medicine, materials science, and climatology, according to the scenarios exposed at the beginning of Section 1. (2) In meteorology; for example, when analyzing wind speed data in certain geographic regions that tend to have two speed spikes per day. (3) In astronomy; for example, when analyzing data associated with certain solar wind parameters. Studies based on annual measurements have shown that the distribution of solar wind speed (as well as those of other parameters such as proton density, temperature and magnetic field) can exhibit bimodality. (4) In population dynamics. It is well known that the interactions between the characteristics of the individuals that make up a certain population can produce changes in the distribution of the size of the individuals, including changes that lead to a bimodal distribution.

Alias Distributions
Previously, it was noted that the original range for the parameter γ of the GB distribution can be extended so that γ assumes values in the interval [0,∞). In this way, the normal distribution can be obtained as a limiting case of the GB distribution. Despite this interesting property, we have considered the range bounded to the interval [0,2) for the parameter γ of the LGB distribution. The rationale for this is to avoid alias distributions for some members of the LGB distribution, that is, avoid distributions that are virtually identical in relation to minimizing the Kullback-Leibler divergence [21]. The Kullback-Leibler divergence measures the degree of divergence between the distributions of two random variables. For two LGB, random variables consider the following proposition.
Proof. From the definition of the Kullback-Leibler divergence, we obtain where the expectation is taken with respect to X 1 . Thus, the result is obtained by considering the change of variable u = z 1 = (x − µ 1 )/σ 1 , once z 2 = (x − µ 2 )/σ 2 can be written as Considering θ 1 known, minimizing Equation (13) with respect to θ 2 is equivalent to solving a nonlinear system of four equations, see Appendix B.
We observe that the existence of alias distributions in the LGB distributions family is related to its skewness behavior when considering that γ can also assume values greater than or equal to 2. In the case γ ≥ 2, it can be verified that the LGB family is unimodal regardless of the value of α. In Figure 4 (left panel), it is observed that the skewness of the LGB distribution is positive for those values of α ∈ (0, 1). However, we observe that the skewness can become negative if α assumes a value close to 0, once γ has assumed a value greater than 2 that is large enough. This favors the existence of alias distributions for certain unimodal members of the family, as it suggests that a negative skewness value may be associated with two different values of α. Under this scenario, it must be taken into account that this is very serious when both members of the family are unimodal, as it would be enough to make an adequate consideration of the locations and scales so that the members are virtually identical.

Maximum Likelihood Estimator
In this section, we deal with the problem of parameter estimation in the LGB distribution under the maximum likelihood (ML) method.
If X ∼ LGB(µ, σ, γ, α), then (for θ = (µ, σ, γ, α)') the log-likelihood is given by Thus, the elements of the score vector are For a random sample X 1 , . . . , X n from X ∼ LGB(µ, σ, γ, α), we observe that the ML estimator θ ML = ( µ ML , σ ML , γ ML , α ML )' for θ = (µ, σ, γ, α)' cannot be expressed in closed form. The solution of the likelihood equations gives rise to a system of four nonlinear equations (See Appendix D) that must be solved with the help of some computational routine in search of ML estimates.
In this case, as the ML estimators do not have a closed form, a good alternative to obtain ML estimates is to solve the following optimization problem, where (θ; x) is given in Equation (14). We solved (15) using the function optim of the R language [19] and, specifically, the L-BFGS-B algorithm [22] was applied. This algorithm requires the declaration of a feasible starting point in the parametric space to start the iterative process. Considering that the bimodal-normal distribution is a special case of the LGB distribution, we verify through simulation experiments that (x, s x , 0, 1), where x is the mean of the observations and s x the corresponding standard deviation, is a good starting point. Under regularity conditions, the asymptotic distribution of ( θ ML − θ) is N 4 (0, K(θ) −1 ), where K(θ) is the expected information matrix. As the function ∑ n i=1 (θ; x i ) is not simple, it is not easy to obtain the analytical expression of this matrix. However, we obtain an approximation from the observed information matrix, whose elements are computed as minus the second partial derivatives of the log-likelihood function with respect to all the parameters (evaluated at the ML estimates). Thus, for a random sample X 1 , . . . , X n from X ∼ LGB(µ, σ, γ, α), the observed information matrix is given by with θ 1 = µ, θ 2 = σ, θ 3 = γ and θ 4 = α, where the analytical expressions of the second partial derivatives are presented in Appendix E.

Simulation Studies
In the analysis of data exhibiting bimodality, it is common to use a two-component mixture normal (MN) distribution. In this section, we initially carry out a simulation study to evaluate the behavior of the ML estimators of the LGB distribution parameters. Subsequently, we conducted a second simulation study in order to evaluate the usefulness of the LGB distribution in a context where the MN distribution performs well.

First Simulation Study
In this study, 1000 random samples from the LGB distribution were generated considering the sample sizes n = 100, 200, 300, 500, and 1000 in the following two scenarios:

2.
Scenario B: µ = −5, σ = 4, γ = 0.75 and α = 1.5. Simulated random samples were generated using the qf given in Equation (7). The Lam-bertW package [23] in the R language was used to compute the principal branch of the Lambert W function. A code in the R language is provided in Appendix A.
For each simulated sample, we obtain the ML estimates by solving (15) under the considerations mentioned in Section 4. Table 3 reports the average estimate (AE), the empirical standard deviation (SD), and the root of the mean square error (RMSE) for the 1000 estimates obtained in each scenario and sample size considered. Table 4 reports the average of the asymptotic standard error (SE) for the ML estimates along with the coverage probability (CP) of the 95% asymptotic confidence intervals. Table 3 indicates that the AEs tend to be close to the true values of the parameters as the sample size increases. The SDs and RMSEs are close and decrease towards 0 as the sample size increases, as expected in the standard asymptotic theory. Table 4 indicates that the SEs are close to the SDs and RMSEs given in Table 3. As expected, the SEs decrease towards 0 and the CPs converge to the nominal values used to construct the confidence intervals as the sample size increases.
Finally, considering a procedure analogous to the one described above, we simulate random samples from the LGB distribution and calculate the hit and non-rejection rates for the MN distribution. The scenarios considered here are the following: Scenario A, µ = −2.5, σ = 2.8, γ = 0.2, and α = 2.5. Scenario B, µ = 15, σ = 22, γ = 0.2, and α = 0.5. Scenario C, µ = 14, σ = 22.4, γ = 1.0, and α = 0.3. Scenario D, µ = −20, σ = 30, γ = 1.5, and α = 1.8. Tables 5 and 6 report the hit and non-rejection rates for the LGB and MN distributions, respectively. In Table 5, we observe that the non-rejection rates are high, which means that a considerable proportion of samples generated from the MN distribution can be appropriately fitted with the LGB distribution. On the other hand, we observe that the hit rates are high, exceeding the value 0.5, even in moderate sample sizes, n = 300. Note that the non-rejection rates decrease considerably in scenario D as the sample size increases and that the hit rates are lower than the other scenarios. This is because the samples are generated from a MN population where the scales σ 1 and σ 2 are considerably different. This shows that the LGB distribution can perform well in settings where the MN distribution is used and where the estimates for σ 1 and σ 2 in this distribution are similar. Table 5. Non-rejection rate based on modified statistics W * and A * and hit rates based on the AIC, CAIC, and BIC for the LGB distribution, when the data are simulated from the MN distribution. In Table 6, we observe that the non-rejection rates are very high, which was expected as the MN distribution has one more parameter than the LGB distribution. This means that a very considerable proportion of samples generated from a LGB population can be appropriately fitted with the MN distribution. However, due to having one more parameter, the hit rates in the different scenarios are small, as the AIC, CAIC, and BIC values depend on the number of parameters of the distribution. Therefore, it can be expected that in a possible real setting where both distributions appropriately fit a certain dataset, the information criteria AIC, CAIC, and BIC will provide favorable indications for the use of the LGB distribution due to the fact of having to estimate a smaller amount of parameters.

Data Analysis
In this section, two applications are presented in order to illustrate the usefulness of the LGB distribution and its special cases in data modeling in different real settings. Other symmetric/asymmetric unimodal/bimodal distributions are also considered to illustrate that the LGB distribution or some of its special cases may have a better fit than other distributions in the literature. Specifically, the odd log-logistic skew-normal (OLLSN) [28] and gamma sinh-Cauchy (GSC) [29] distributions are considered. Like the LGB distribution, these distributions have four parameters: two shape parameters (which together control skewness and bimodality), a location parameter, and a scale parameter. A mixture distribution of two normal components (MN) [6] is also included into the analysis as it is a commonly used distribution for analyzing data exhibiting bimodality.

1.
The first dataset corresponds to 188 observations on the inflation rate (in %) registered quarterly between the years 1950 and 1996 in Canada. This dataset can be found with the name Tbrate in the R language [30].

2.
The second dataset refers to 128 observations on the electrical resistance (in ohms) of nectarine fruits. This data can be found with the name fruitohms in the R language [31].
For the datasets described above, we test hypothesis H 0 : the data have exactly one mode versus the alternative hypothesis H 1 : the data have at least two modes. For this, we consider the excess mass test [32] using the modetest function in R language [33]. For the inflation rate data, the observed statistic was 0.05 with a p-value equal to 0.05. For the electrical resistance data, these values were 0.074 and 0.01, respectively. Thus, at a significance level equal to 0.05, in both datasets H 0 is rejected; that is, the distributions of the inflation rate and electrical resistance data are at least bimodal.
We compared the distributions fitted by the ML method using the information criteria AIC, CAIC, and BIC. We also calculate the statistics W * and A * to test the hypothesis H 0 : X 1 , . . . , X n is a random sample from a continuous distribution F(x; θ), where F(·; ·) is known but θ is unknown. In these tests, H 0 is rejected at a significance level equal to 0.05 if W * > 0.126 and A * > 0.752. Table 7 reports the ML estimates with the corresponding standard errors for each distribution fitted to the inflation rate and electrical resistance data. In addition, the values associated with the statistics A * and W * and with the information criteria are reported. In Table 7, with respect to the inflation rate data, based on the values of the statistics W * and A * , it can be seen that the hypothesis that the data correspond to an observed random sample of the GSC, GB, LBN, or BN distributions is rejected at a significance level equal to 0.05. In addition, it can be seen that the LGB distribution is the one with the lowest AIC, CAIC, and BIC values among the fitted distributions, indicating that this distribution should be selected over the others for the modeling of these data.
With respect to the electrical resistance data, it is observed that the hypothesis that the data correspond to an observed random sample of the OLLSN, GB, BN, or GSC distributions is rejected at a significance level equal to 0.05. Note that the AIC, CAIC, and BIC values for the LGB, MN, and LBN distributions are close, the values associated with the LBN distribution being slightly lower. Thus, the LBN distribution (which has a smaller parametric dimension than the LGB and MN distributions) is capable of fitting the electrical resistance data as well as the LGB and MN distributions.
In the left panels of Figures 7 and 8, the histograms of inflation rate and electrical resistance are displayed along with the fitted densities. In the right panels of the same figures, the empirical cdf and the fitted cdf s are compared. In these plots, we see that the LGB distribution fits the inflation rate data appropriately, while the LBN distribution fits the electrical resistance data appropriately. Note that the LGB and LBN distributions have one and two fewer parameters, respectively, than the MN distribution, and that despite this fact they are capable of presenting good fits to the analyzed data.   Figure 8. Left panels: Histogram for electrical resistance data and the fitted pdf curves via the ML method. Right panels: Empirical cdf for the electrical resistance data and the fitted cdf curves.

Final Comments
This article introduces a new symmetric/asymmetric unimodal/bimodal distribution called the Lambert generalized bimodal (LGB) distribution. Some special cases of the LGB distribution are discussed. One of the special cases, the Lambert-bimodal normal (LBN) distribution, can be considered as an alternative to other symmetric/asymmetric bimodal distributions, including the LGB distribution. The LGB distribution arises using the Lambert-F transformation when the generalized bimodal distribution is considered as baseline distribution. We study the main structural properties of the LGB distribution, such as the pdf, cdf, qf, and raw moments that are used for a description of the skewness and kurtosis characteristics. Parameter estimation of the LGB distribution is discussed using the ML method. Through simulation experiments, we observe that the ML method provide acceptable estimates of the parameters of the LGB distribution. Furthermore, through simulation experiments, we observed that the LGB distribution can adequately fit datasets generated from the mixture normal distribution, despite having one less parameter. Finally, two applications that illustrate that the LGB distribution and the LBN especial case can present a better fit of data in real settings than other symmetric/asymmetric unimodal/bimodal distributions such as the odd log-logistic skew-normal distribution (OLLSN), gamma sinh-Cauchy (GSC), and mixture normal (MN) distributions.
As a final consideration, we leave the question open, does the LGB distribution have an intuitive stochastic representation? In the literature, it is possible to find distribution families that have an attractive intuitive generation mechanism, such as the skew-elliptical (SE) distributions [34] and the closed-skew-normal (CSN) distribution [35]. According to Loperfido et al. [36], any linear combination of the largest and smaller component of a bivariate, exchange elliptical random vector has a skew-elliptical distribution. According to Loperfido [37], any order statistic from a random vector with exchangeable normal distribution has a closed-skew-normal distribution. As far as we know, no similar property is known for the Lambert-F distributions class or for some of its special cases.