Next Article in Journal
Approximation Formula for Option Prices under Rough Heston Model and Short-Time Implied Volatility Behavior
Previous Article in Journal
Modeling of Extreme Values via Exponential Normalization Compared with Linear and Power Normalization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Slashed Lomax Distribution and Regression Model

1
Department of Applied Mathematics, Xi’an University of Technology, Xi’an 710048, China
2
Department of Mathematical Sciences, Eastern New Mexico University, Portales, NM 88130, USA
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(11), 1877; https://doi.org/10.3390/sym12111877
Submission received: 4 October 2020 / Revised: 22 October 2020 / Accepted: 12 November 2020 / Published: 14 November 2020

Abstract

:
In this article, the slashed Lomax distribution is introduced, which is an asymmetric distribution and can be used for fitting thick-tailed datasets. Various properties are explored, such as the density function, hazard rate function, Renyi entropy, r-th moments, and the coefficients of the skewness and kurtosis. Some useful characterizations of this distribution are obtained. Furthermore, we study a slashed Lomax regression model and the expectation conditional maximization (ECM) algorithm to estimate the model parameters. Simulation studies are conducted to evaluate the performances of the proposed method. Finally, two sets of data are applied to verify the importance of the slashed Lomax distribution.

1. Introduction

The Lomax distribution, which was introduced by Lomax [1], has been regarded as the mixed distribution of the exponential distribution and gamma distribution. It has a heavy-tailed probability distribution, often used in business, economics, and actuarial modeling. Let X be a non-negative random variable with a Lomax distribution, then its probability density function (pdf) is given by,
f x ; α , β = α β α β + x α + 1 ,
where α > 0 and β > 0 , and it is denoted by X L o m a x α , β . The Lomax distribution contains the monotone decreasing failure rate and the monotone increasing failure rate, which have been regarded as a quantitative life distribution. It has been applied broadly in practical production and real life. For example: Myhre and Saunders [2] applied the Lomax distribution to right censored data; Balakrishnan and Ahsanullah [3] discussed some important statistical properties of the Lomax distribution; Childs et al. [4] studied the properties of right-truncated Lomax distributions and discussed some practical applications; and Howlader and Hossain [5] considered estimating the survival function of the Lomax distribution with the Bayesian method.
The slashed distribution was proposed to model the bell shaped data with a heavier tail, by Rogers and Tukey [6], which is thicker than the tail of the normal distribution and is itself a symmetric distribution. It has a stochastic representation as Z U 1 q , where q > 0 , Z N 0 , 1 is independent of U U 0 , 1 , and is called the canonical slashed distribution when q = 1 . After that, statisticians have done in-depth research on and promoted this model; see Wang and Genton [7], Gomez et al. [8], Arslan and Genc [9], Reyes et al. [10], and Tian et al. [11,12]. Recently, there has been a new extension of the slashed distribution: Gui [13] introduced a three-parameter extension model called the Lindley slashed distribution; Iriarte et al. [14,15] proposed the slashed Rayleigh distribution and modified slashed Rayleigh distribution, which have a more flexible kurtosis than the Rayleigh distribution; Reyes et al. [16] discussed the modified slashed Birnbaum–Saunders distribution and concluded that it has greater kurtosis values than the usual BSdistribution. Similar to this methodology, the slashed Lindley–Weibull distribution was introduced by Reyes et al. [17]; the slashed power Lindley distributions was studied by Iriarte et al. [18]; and the modified slashed generalized exponential distribution was introduced by Astorga et al. [19].
Regression models are undoubtedly the most widely used, but the normality assumption of the residual errors is more restrictive. Recently, the error term subject to a more flexible distribution has been studied: Gómez [20] analyzed the regression model of the slashed half-normal distribution; Jamal [21] studied the properties of the Topp–Leone–Weibull Lomax distribution and its regression model; and Hamedani et al. [22] analyzed the regression model of the Burr XII distribution.
Based on the previous research, we propose the slashed Lomax distribution, which has a thicker tail than the Lomax distribution and has more flexibility in kurtosis. Furthermore, the construction of the slashed Lomax distribution makes it possible to estimate the parameters in a variety of ways. Therefore, the proposed new distribution can be used not only for fitting thick-tailed datasets, but also can analyze some phenomena in real life. The rest of the paper is organized as follows. In Section 2, we introduce the slashed Lomax distribution and obtain some of its properties. The ECM algorithm for the parameter estimation and simulation studies is proposed in Section 3. The slashed Lomax regression model is studied in Section 4. Two applications to real data are investigated in Section 5. Some conclusions are offered in Section 6.

2. Slashed Lomax Distribution

In this section, some basic properties of the slashed Lomax distribution are described.
Definition 1.
A random variable Y follows a slashed Lomax distribution, denoted as Y S l o m a x α , β , λ , if it can be represented as:
Y = V 1 2 X ,
where X L o m a x α , β and V b e t a λ , 1 are independent random variables and α , β , λ > 0 .
Proposition 1.
Let Y S l o m a x α , β , λ , then its pdf can be written as:
f Y y ; α , β , λ = 2 λ α β 2 λ y 2 λ + 1 0 y β w 2 λ 1 + w α + 1 d w .
Proof. 
Using the stochastic representation in Equation (1), we have the joint density function of ( Y , V ) as:
f Y , V y , v = λ α v λ 1 2 β ( 1 + v 1 2 y β ) α + 1 ,
and the pdf of Y is obtained by marginalizing the distribution with respect to V. □
Corollary 1.
Let Y S l o m a x α , β , λ .
(i)
For y β < 1 ,
f Y y ; α , β , λ = 2 λ α 2 λ + 1 β F α + 1 , 2 λ + 1 , 2 λ + 2 ; y β ,
where F α , β , γ ; x = 1 B β , γ β 0 1 1 t γ β 1 t β 1 1 x t α d t is the hyper geometric function.
(ii)
2 λ and α 2 λ are positive integers, and 0 < y β < 1 ,
f Y y ; α , β , λ = 2 λ α 2 λ + 1 β B y β ; 2 λ + 1 , α 2 λ ,
where B x , y = 0 t x 1 1 + t x + y d t is the incomplete beta function.
Remark 1.
For λ = 1 , Y S l o m a x α , β , λ is reduced to the canonical slashed Lomax distribution, and the pdf is:
f Y y ; α , β , λ = 0 1 α v 1 2 β 1 + v 1 2 y β α + 1 d v .
For an illustration of the new family of slashed Lomax distributions, we draw the density curves, with different values of α , β , and λ , as follows. After fixing two of the three parameters, we find that the slashed Lomax distributions have a heavier tail as β increases (and the same for λ , but the opposite for α ).
Proposition 2.
Let Y S l o m a x α , β , λ .
(i)
The reliability (survival) function of Y is given by:
R Y y = 1 2 λ α β 2 λ 0 y A u β ; α , β , λ u 2 λ + 1 d u ,
where A y β ; α , β , λ = 0 y β w 2 λ 1 + w α + 1 d w .
(ii)
The hazard rate function of Y is given by:
h Y y ; α , β , λ = 2 λ α β 2 λ y 2 λ + 1 A y β ; α , β , λ 1 2 λ α β 2 λ 0 y A u β ; α , β , λ u 2 λ + 1 d u 1 .
(iii)
The reversed hazard rate of Y is given by:
r Y ( y ) = f ( y ) F ( y ) = A y β ; α , β , λ y 2 λ + 1 0 y A u β ; α , β , λ u 2 λ + 1 d u .
Figure 1 shows probability density function curve of the slashed Lomax distribution with several different parameters. Figure 2 shows the curves of hazard rate functions under different parameters. It can be seen that the values of the hazard rate function based on λ , with α and β fixed, have small differences as y increases in these figures.
Next, the Renyi entropy and distributional moments of the slashed Lomax distribution are derived. In addition, the coefficients of skewness and kurtosis are calculated.
Proposition 3.
Let Y S l o m a x α , β , λ ; the Renyi entropy of order δ for Y is given by,
R δ ( y ) = 1 1 δ δ l o g α λ β + l o g 0 0 1 v λ 1 2 1 + v 1 2 y β α + 1 δ d y ,
where δ > 0 and δ 1 .
Proof. 
Renyi entropy is an important diversity index in ecology and statistics, which is defined as R δ ( y ) = 1 1 δ l o g 0 f ( y ) δ d y . Therefore, we have:
R δ ( y ) = 1 1 δ l o g α λ β δ 0 0 1 v λ 1 2 1 + v 1 2 y β α + 1 δ d y = 1 1 δ δ l o g α λ β + l o g 0 0 1 v λ 1 2 1 + v 1 2 y β α + 1 δ d y .
 □
Proposition 4.
Let Y Slomax α , β , λ ; the r-th moment of Y, for r = 1 , 2 and r < 2 λ , is given by:
μ r = 2 λ β Γ α r Γ r + 1 ( 2 λ r ) Γ α ,
where Γ ( · ) is the gamma function.
Proof. 
Using Equation (1), we have,
μ r = E Y r = E X r E V r 2 .
According to the properties of the beta distribution and the Lomax distribution, we have that E V r = λ λ r , for r < λ , and E X r = β Γ α r Γ r + 1 Γ α . Thus, the result is obtained. □
Corollary 2.
Let Y Slomax α , β , λ ; the mean and variance of Y are given by:
E Y = 2 λ β 2 λ 1 α 1 ,
V a r Y = 2 λ 1 2 α 1 λ β 2 λ 2 α 2 .
Corollary 3.
Let Y Slomax α , β , λ ; the coefficients of skewness and kurtosis for Y, γ 1 and γ 2 , are given by:
γ 1 = λ β 5 2 2 λ 2 1 2 α 2 1 2 Δ 1 α 1 9 2 α 3 2 λ 1 6 2 λ 3 ,
γ 2 = 16 λ 3 β 3 2 λ 2 α 2 Δ 2 α 1 6 α 3 α 4 2 λ 1 8 2 λ 3 2 λ 4 ,
where Δ 1 = 12 ( α 1 ) 2 ( 2 λ 1 ) 3 ( 2 λ 2 ) 24 λ β ( α 1 ) ( α 3 ) ( 2 λ 1 ) 2 ( 2 λ 3 ) + 16 β 2 λ 2 ( α 2 ) ( α 3 ) ( 2 λ 2 ) ( 2 λ 3 ) and Δ 2 = 3 ( α 1 ) 3 ( 2 λ 1 ) 4 ( 2 λ 2 ) ( 2 λ 3 ) 3 λ β ( α 1 ) 2 ( α 4 ) ( 2 λ 1 ) 3 ( 2 λ 2 ) ( 2 λ 4 ) + 4 λ 2 β 2 ( α 1 ) ( α 3 ) ( α 4 ) ( 2 λ 1 ) 2 ( 2 λ 3 ) ( 2 λ 4 ) 3 λ 3 β 3 ( α 2 ) ( α 3 ) ( α 4 ) ( 2 λ 2 ) ( 2 λ 3 ) ( 2 λ 4 ) .
Proof. 
From Proposition 4, γ 1 = μ 3 3 μ 2 μ 1 + 2 μ 1 3 μ 2 μ 1 2 3 2 and γ 2 = μ 4 4 μ 1 μ 3 + 6 μ 2 μ 2 3 μ 1 4 μ 2 μ 1 2 2 ; thus, the results are obtained directly. □
In the following propositions, we consider the linear transformation of the slashed Lomax distribution and discuss its scale mixture representation.
Proposition 5.
Let Y S l o m a x ( α , β , λ ) and the scalar a > 0 , then a Y S l o m a x ( α , a β , λ ) .
Proof. 
Let X L o m a x ( α , β ) , and a > 0 is a scalar; we have a X L o m a x ( α , a β ) . Therefore, the result is obtained by Equation (1). □
Proposition 6.
If Y | V = v L o m a x α , v 1 2 β and V∼ beta ( λ , 1 ), then Y∼ Slomax( α , β , λ ).
Proof. 
The result can be obtained by Equation (1) and Proposition 1. □

3. ECM Algorithm for Parameter Estimation

In this section, we consider the maximum likelihood (ML) estimation for the parameters θ = α , β , λ of the slashed Lomax distribution. It is well known that the EM algorithm is an important tool for ML estimation when no data or potential variables are observed. The E-step is used to find the expectation of the incomplete data based on the observed values, and the M-step works for the maximization. Usually, the M-step calculation is difficult when the maximum likelihood of the complete data is complex. Meng and Rubin [23] proposed an ECM algorithm that was an extension of the EM algorithm and satisfied all the properties of the EM algorithm. The basic idea of the ECM algorithm is to decompose the M-step of the EM algorithm into k times conditional maximization.
In the following, an ECM algorithm is proposed for obtaining the estimator of θ from the slashed Lomax distribution. Let y 1 , , y n be a random sample from S l o m a x α , β , λ . According to Equation (1), we have,
y i | v i = v i L o m a x α , v i 1 2 β , v i b e t a λ , 1 .
Let y = y 1 , . . . , y n and v = v 1 , . . . , v n be observed data, then ( y i , v i ), i = 1 , , n are the complete dataset. Therefore, the log-likelihood function for the complete data ( y i , v i ) , i = 1 , , n , is given by:
l θ = n l o g λ + n l o g α n l o g β + α + 1 n l o g β + λ 1 2 i = 1 n l o g v i α + 1 i = 1 n l o g β + v i 1 2 y i .
E-Step: Compute the conditional expectation of the log-likelihood function, Q θ , θ ^ ( k 1 ) = E [ l θ | y , θ ^ ( k 1 ) ] , by using the following equation,
Q θ , θ ^ ( k 1 ) = n l o g λ ( k 1 ) + n l o g α ( k 1 ) + α ( k 1 ) n l o g β ( k 1 ) + λ ( k 1 ) 1 2 i = 1 n δ 1 i ( k 1 ) α + 1 i = 1 n δ 2 i ( k 1 ) ,
where δ 1 i = E l o g v i | y i , θ ^ , δ 2 i = E l o g ( β + v i 1 2 y i ) | y i , θ ^ .
CM-Steps: Maximize Q θ , θ ^ ( k 1 ) with respect to θ to obtain ( θ ^ ) ( k ) .
Denote W ( k 1 ) = 0 y β ^ ( k 1 ) w 2 λ ^ ( k 1 ) 1 + w α ^ ( k 1 ) + 1 d w ; the following equations are obtained for the maximization.
(i)
Set up Q θ , θ ^ ( k 1 ) α = n α ^ ( k ) + n l o g β ^ ( k 1 ) i = 1 n δ 2 i ( k 1 ) = 0 , then:
α ^ ( k ) = n i = 1 n δ 2 i ( k 1 ) n l o g β ^ ( k 1 ) ,
where δ 2 i ( k 1 ) = y i 2 λ ^ ( k 1 ) + 1 2 W ( k 1 ) β ^ ( k 1 ) 2 λ ^ ( k 1 ) + 1 0 1 l o g β ^ ( k 1 ) + v 1 2 y i v λ ^ ( k 1 ) 1 2 ( 1 + y i β ^ ( k 1 ) ) α ^ ( k 1 ) + 1 d v .
(ii)
Set up Q θ , θ ^ ( k 1 ) β = α ^ ( k ) n β ^ ( k ) ( α ^ ( k ) + 1 ) i = 1 n δ 3 i ( k 1 ) = 0 , then:
β ^ ( k ) = n α ^ ( k ) α ^ ( k ) + 1 i = 1 n δ 3 i ( k 1 ) ,
where δ 3 i ( k 1 ) = y i 2 λ ^ ( k 1 ) + 1 2 W ( k 1 ) β ^ ( k 1 ) 2 λ ^ ( k 1 ) + 1 0 1 β ^ ( k 1 ) + v 1 2 y i 1 v λ ^ ( k 1 ) 1 2 ( 1 + y i β ^ ( k 1 ) ) α ^ ( k 1 ) + 1 d v .
(iii)
Set up Q θ , θ ^ ( k 1 ) λ = n λ ^ ( k ) + i = 1 n δ 1 i ( k 1 ) = 0 , then:
λ ^ ( k ) = n i = 1 n δ 1 i ( k 1 ) ,
where δ 1 i ( k 1 ) = y i 2 λ ^ ( k 1 ) + 1 2 W β ^ ( k 1 ) 2 λ ^ ( k 1 ) + 1 0 1 l o g ( v ) v λ ^ ( k 1 ) 1 2 1 + y i β ^ ( k 1 ) α ^ ( k 1 ) + 1 d v .
The iterations are repeated until a suitable convergence rule is satisfied, say | | θ ^ ( k ) θ ^ ( k 1 ) | | < ϵ , where | | · | | is the Euclidean norm and ϵ is sufficiently small.
In the following, we conduct simulation studies to test the efficiency of the estimation procedure discussed above. The random values following the slashed Lomax distribution can be generated by Equation (1). The parameters are set up as: θ = (2,1,0.5), (2,1,1), (2,1,1.5), (1,0.5,1), (1,1,1), (1,2,1), (0.5,2,2), (2,2,1), (4,2,1). The sample sizes n = 50, 100, and 200. The following procedure is for generating a random number with size n from S l o m a x ( α , β , λ ) ,
  • set α , β , λ , and n;
  • simulate U U n i f o r m ( 0 , 1 ) ;
  • simulate V b e t a ( λ , 1 ) ;
  • compute Z = ( 1 U ) 1 α 1 β ;
  • compute Y = V 1 2 Z .
For each scenario, we repeat the process N = 1000 times. The ECM algorithm is applied, and the results are computed using the software R. The procedures are put in the same conditions (same initial values, and ϵ = 10 4 ). The estimators are obtained by applying the nleqslvfunction in Equations (4)–(6). The mean values of the parameters and the corresponding standard deviation (SD) are shown in Table 1. The empirical estimated mean value and SDs based on the 1000 replicates are calculated by:
M e a n n ( ξ ) = 1 N i = 1 n ξ ^ i , S D n ( ξ ) = 1 N i = 1 N ( ξ ^ i ξ ) 2 ,
where ξ is the interesting parameter, α , β , λ , respectively.
From Table 1, it can be seen that as the sample size increases, the mean value of estimators comes closer to the true values, and the SD decreases in all case.

4. Slashed Lomax Regression Model

In this section, we study the slashed Lomax regression model and propose the ECM algorithm to estimate the parameters.
Definition 2.
The slashed Lomax regression model is defined as:
y = x γ + ε ,
where y , γ p is a p-dimensional vector of regression parameters, x = ( x 1 , x 2 , . . . , x n ) is a known p × n matrix, with x i p for i = 1 , , n , and ε S l o m a x ( α , β , λ ) .
According to Equation (2), the density function of y is given by:
f ( y ) = 0 1 α λ β α v λ 1 2 β + ( y x γ ) v 1 2 α + 1 d v .
Let y = ( y 1 , y 2 , . . . y n ) come from the slashed Lomax regression model with X the given design matrix x. The log-likelihood function of δ = ( α , β , λ , γ ) can be written as:
l ( δ | y , X ) = n l o g ( λ ) + n l o g ( α ) + α n l o g ( β ) + i = 1 n l o g M ( y i ) ,
where M ( y i ) = 0 1 v λ 1 2 β + ( y i x i γ ) v 1 2 α + 1 d v .
Now, given v = ( v 1 , . . . v n ) , we can get the complete log-likelihood function:
l ( δ | y , v ) = n l o g ( α ) + n l o g ( λ ) + α n l o g ( β ) + λ 1 2 i = 1 n l o g ( v i ) ( α + 1 ) i = 1 n l o g β + ( y i x i γ ) v 1 2 ,
and the conditional expectation of l ( δ | y , v ) ,
Q ( δ ) = E l ( δ ) | y , δ ^ = n l o g ( λ ) + n l o g ( α ) + α n l o g ( β ) + λ 1 2 i = 1 n δ 1 i ( α + 1 ) i = 1 n δ 2 i ,
where δ 1 i = E l o g v i | y i , δ ^ , δ 2 i = E l o g ( β + v i 1 2 ( y i x i γ ) ) | y i , δ ^ .
In the following, we illustrate the steps for the ECM algorithm for estimating the parameters in the regression model.
E-step: Given δ ^ ( k ) = α ^ ( k ) , β ^ ( k ) , λ ^ ( k ) , γ ^ ( k ) , find the expectation of the condition.
CM-step I: α ^ ( k ) must be updated as:
α ^ ( k + 1 ) = n i = 1 n δ 2 i n l o g ( β ^ ( k ) ) .
CM-step I I : β ^ ( k ) must be updated as:
β ^ ( k + 1 ) = α ^ ( k + 1 ) n ( α ^ ( k + 1 ) + 1 ) i = 1 n δ 3 i .
CM-step I I I : λ ^ ( k ) must be updated as:
λ ^ ( k + 1 ) = n i = 1 n δ 1 i .
CM-step I V : γ ^ ( k ) must be updated as:
γ ^ ( k + 1 ) = arg min δ i = 1 n δ 2 i .
The iterations are repeated until a suitable convergence rule is satisfied, say | δ ^ ( k + 1 ) δ ^ ( k ) | sufficiently small. On the basis of this theory, a simulation experiment is carried out. We take the two-dimensional regression model to carry out the numerical simulation. The value of the parameters for ( α , β , λ ) are chosen as (2,2,1), (1,1,1), (1,4,0.5),(0.5,3,0.5), and the regression parameters ( γ 1 , γ 2 ) are set to (1,4), (−2,3), and (−3,−5); x i N ( 0 , 1 ) , for i = 1 , 2 , are uncorrelated. For each scenario, we generate N = 1000 samples of size n = 50, 100, and 200, respectively. The empirical estimated mean value and SDs for the corresponding parameters based on the 1000 replicates are shown in Table 2, Table 3 and Table 4.
From the simulation results, we can see that, as the sample size increases, the standard deviations between the estimated parameters decrease, and the estimated values are close to the actual values.

5. Application

Next, we use the real data to verify the practicability of the slashed Lomax distribution. We compare the performance of the slashed Lomax distribution to that of the beta Lomax (Blomax) distribution (Lemonte and Cordeiro [24]), the exponentiated standard Lomax (ESlomax) distribution, and the Gumbel Lomax (Gulomax) distribution (Tahir et al. [25]). The measures of the goodness of fit including the Akaike information criterion (AIC) and Bayesian information criterion (BIC) values are computed to compare the fitted model and regression model.
Data 1: This dataset has 214 observations and describes the successive failure of the air conditioning systems in a fleet of 13 Boeing 720 jet airplanes, which was studied by Kus [26], Tahir et al. [25], and many others.
The parameter estimates, AIC, and BIC for all fitted distributions are shown in Table 5. Both criteria provide evidence in favor of the slashed Lomax distribution for this dataset, corroborating that the slashed Lomax distribution can be seen as a competitive distribution of practical interest in the real world.
Figure 3 displays the fitted models for the dataset. The left panel of Figure 3 shows the fitted densities to the dataset histogram and some estimated distributions, and the right panel displays the empirical distribution function for the dataset and the estimated distributions. Both figures reveal that the slashed Lomax distribution provides a qualified fit for the dataset.
Data 2: This dataset is about the lifespan of patients with heart failure and contains 96 patients with heart failure, the detailed information of which was provided by Tanvir Ahmad [27].
The numerical variables considered here are ejection fraction ( x 1 ) , serum creatinine ( x 2 ) , serum sodium ( x 3 ) , and age ( x 4 ) as potential variables that explain cardiovascular mortality ( y ) . We apply the slashed Lomax regression model developed in Section 4 to this dataset. The resulting estimates and other regression models are given in Table 6.
As can be noted from Table 6, the slashed Lomax regression model has the lowest AIC and BIC values among those of the other regression models. The values of these statistics indicate that the slashed Lomax regression model provides the best fit to the data. In addition, the increase of the content of serum creatinine will shorten the lifespan of patients and increase the risk of death, which are consistent with the conclusions by Tanvir Ahmad [27].

6. Discussion

In this article, the slashed Lomax distribution and its corresponding statistical properties are introduced. Maximum likelihood estimators through the ECM algorithm and simulation studies are discussed. In addition, the slashed Lomax regression model is studied. Applications to real data demonstrate the importance of the proposed distribution. In the future, we will consider extending this proposed model to the multivariate case and study its extremal properties, aiming to explore more practical values.

Author Contributions

W.T.: Conceptualization, Methodology, Validation, Investigation, Resources, Supervision, Project Administration, Visualization, Writing–review and editing; H.L.: Software, Formal analysis, Data curation, Writing–original draft preparation, Visualization. Both the authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the Editor and two anonymous referees for their careful reading of this article and for their constructive suggestions, which considerably improved this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lomax, K.S. Business failures: Another example of the analysis of failure data. J. Am. Stat. Assoc. 1954, 49, 847–852. [Google Scholar] [CrossRef]
  2. Myhre, J.; Saunders, S. Screen testing and conditional probability of survival. Lect. Notes-Monogr. Ser. 1982, 2, 166–178. [Google Scholar]
  3. Balakrishnan, N.; Ahsanullah, M. Relations for single and product moments of record values from Lomax distribution. Sankhy Indian J. Stat. Ser. B 1994, 56, 140–146. [Google Scholar]
  4. Childs, A.; Balakrishnan, N.; Moshref, M. Order statistics from non-identical right-truncated Lomax random variables with applications. Stat. Pap. 2001, 42, 187–206. [Google Scholar] [CrossRef]
  5. Howlader, H.A.; Hossain, A.M. Bayesian survival estimation of Pareto distribution of the second kind based on failure-censored data. Comput. Stat. Data Anal. 2002, 38, 301–314. [Google Scholar] [CrossRef]
  6. Rogers, W.H.; Tukey, J.W. Understanding some long-tailed symmetrical distributions. Stat. Neerl. 1972, 26, 211–226. [Google Scholar] [CrossRef]
  7. Wang, J.; Genton, M.G. The multivariate skew slash distribution. J. Stat. Plan. Inference 2006, 136, 209–220. [Google Scholar] [CrossRef]
  8. Gómez, H.W.; Quintana, F.A.; Torres, F.J. A new family of slash-distributions with elliptical contours. Stat. Probab. Lett. 2007, 77, 717–725. [Google Scholar] [CrossRef]
  9. Arslan, O.; Genc, A.I. A generalization of the multivariate slash distribution. J. Stat. Plan. Inference 2009, 139, 1164–1170. [Google Scholar] [CrossRef]
  10. Reyes, J.; Gomez, H.W.; Bolfarine, H. Modified slash distribution. Statistics 2013, 47, 929–941. [Google Scholar] [CrossRef]
  11. Tian, W.; Wang, T.; Gupta, A.K. A new family of multivariate skew slash distribution. Commun. Stat. Theory Methods 2018, 47, 5812–5824. [Google Scholar] [CrossRef]
  12. Tian, W.; Han, G.; Wang, T.; Pipitpojanakarn, V. EM estimation for multivariate skew slash distribution. In Robustness in Econometrics; Springer: Cham, Switzerland, 2017; pp. 235–248. [Google Scholar]
  13. Gui, W. Statistical properties and applications of the Lindley slash distribution. J. Appl. Stat. Sci. 2012, 20, 283. [Google Scholar]
  14. Iriarte, Y.A.; Gomez, H.W.; Varela, H.; Bolfarine, H. Slashed rayleigh distribution. Rev. Colomb. Estadstica 2015, 38, 31–44. [Google Scholar] [CrossRef]
  15. Iriarte, Y.A.; Castillo, N.O.; Bolfarine, H.; Gomez, H.W. Modified slashed-Rayleigh distribution. Commun. Stat. Theory Methods 2018, 47, 3220–3233. [Google Scholar] [CrossRef]
  16. Reyes, J.; Vilca, F.; Gallardo, D.I.; Gomez, H.W. Modified slash Birnbaum–Saunders distribution. Hacet. J. Math. Stat. 2017, 46, 969–984. [Google Scholar] [CrossRef]
  17. Reyes, J.; Iriarte, Y.A.; Jodr, P.; Gomez, H.W. The Slash Lindley–Weibull Distribution. Methodol. Comput. Appl. Probab. 2019, 21, 235–251. [Google Scholar] [CrossRef]
  18. Iriarte, Y.A.; Rojas, M.A. Slashed power-Lindley distribution. Commun. Stat. Theory Methods 2019, 48, 1709–1720. [Google Scholar] [CrossRef]
  19. Astorga, J.M.; Iriarte, Y.A.; Gomez, H.W.; Bolfarine, H. Modified slashed generalized exponential distribution. Commun. Stat. Theory Methods 2019, 49, 4603–4617. [Google Scholar] [CrossRef]
  20. Gómez, Y.M.; Gallardo, D.I.; de Castro, M. A regression model for positive data based on the slashed half-normal distribution. Available online: https://www.ine.pt/revstat/pdf/Aregressionmodelforpositivedata.pdf (accessed on 1 January 2020).
  21. Jamal, F.; Reyad, H.M.; Nasir, M.A.; Chesneau, C.; Shah, M.A.A.; Ahmed, S.O. Topp-Leone Weibull-Lomax distribution: Properties, Regression Model and Applications. 2019. Available online: https://hal.archives-ouvertes.fr/hal-02270561/ (accessed on 1 January 2020).
  22. Hamedani, G.G.; Rasekhi, M.; Najibi, S.M.; Yousof, H.M.; Alizadeh, M. Type II general exponential class of distributions. Pak. J. Stat. Oper. Res. 2019, 15, 503–523. [Google Scholar] [CrossRef]
  23. Meng, X.L.; Rubin, D.B. Maximum likelihood estimation via the ECM algorithm: A general framework. Biometrika 1993, 80, 267–278. [Google Scholar] [CrossRef]
  24. Lemonte, A.J.; Cordeiro, G.M. An extended Lomax distribution. Statistics 2013, 47, 800–816. [Google Scholar] [CrossRef]
  25. Tahir, M.H.; Hussain, M.A.; Cordeiro, G.M.; Hamedani, G.G.; Mansoor, M.; Zubair, M. The Gumbel-Lomax distribution: Properties and applications. J. Stat. Theory Appl. 2016, 15, 61–79. [Google Scholar] [CrossRef] [Green Version]
  26. Kus, C. A new lifetime distribution. Comput. Stat. Data Anal. 2007, 51, 4497–4509. [Google Scholar] [CrossRef]
  27. Ahmad, T.; Munir, A.; Bhatti, S.H.; Aftab, M.; Raza, M.A. Survival analysis of heart failure patients: A case study. PLoS ONE 2017, 12, e0181001. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The probability density function plots of Slomax ( α , β , λ ).
Figure 1. The probability density function plots of Slomax ( α , β , λ ).
Symmetry 12 01877 g001
Figure 2. Hazard rate function plots of Slomax ( α , β , λ ).
Figure 2. Hazard rate function plots of Slomax ( α , β , λ ).
Symmetry 12 01877 g002
Figure 3. Fitted curves of the Gulomax, Blomax, ESlomax, and Slomax distribution.
Figure 3. Fitted curves of the Gulomax, Blomax, ESlomax, and Slomax distribution.
Symmetry 12 01877 g003
Table 1. Estimates for parameters α , β , λ of the slashed Lomax distribution.
Table 1. Estimates for parameters α , β , λ of the slashed Lomax distribution.
Sample α β λ α ^ (SD) β ^ (SD) λ ^ (SD)
n = 50210.51.8883 (0.0848)1.3074 (0.3761)0.8684 (0.3455)
2111.9468 (0.0509)1.0003 (0.2849)1.1409 (0.5941)
211.51.9620 (0.0439)0.9147 (0.2761)1.3130 (0.7770)
10.510.9431 (0.0640)0.4590 (0.1624)1.1234 (0.6727)
1110.9334 (0.0660)0.8869 (0.3077)1.0288 (0.5105)
1210.9305 (0.0676)1.7732 (0.6035)0.9922 (0.4897)
0.5210.4480 (0.0408)1.6718 (0.6684)0.9473 (0.3461)
2211.9607 (0.0403)1.9989 (0.5392)1.1505 (0.6040)
4214.0289 (0.0886)1.9676 (0.5494)1.0652 (0.3331)
n = 100210.51.8869 (0.5755)1.3423 (0.2870)0.7882 (0.1667)
2111.9626 (0.4163)1.0579 (0.2244)1.2117 (0.5973)
211.51.9745 (0.0328)0.9262 (0.1931)1.4674 (0.6219)
10.510.9431 (0.0489)0.4598 (0.1270)1.1133 (0.6548)
1110.9556 (0.0464)0.9610 (0.2665)1.1394 (0.5322)
1210.9572 (0.0447)1.9481 (0.5402)1.1264 (0.5208)
0.5210.4677 (0.0304)1.8369 (0.5869)1.0834 (0.4289)
2211.9673 (0.3313)2.0052 (0.4470)1.1695 (0.5673)
4213.9852 (0.0452)2.1742 (0.4372)1.1036 (0.1969)
n = 200210.51.8961 (0.0425)1.3209 (0.1952)0.7660 (0.0732)
2111.9681 (0.0368)1.0377 (0.1805)1.1338 (0.5696)
211.51.9841 (0.0236)0.9704 (0.1689)1.5644 (0.5242)
10.510.9632 (0.0351)0.5031 (0.1172)1.1189 (0.4276)
1110.9651 (0.0353)0.9893 (0.2253)1.1067 (0.4699)
1210.9667 (0.0353)1.9817 (0.4290)1.1444 (0.4659)
0.5210.4755 (0.0235)1.9270 (0.4552)1.1307 (0.4125)
2211.9755 (0.0243)2.0940 (0.3621)1.1379 (0.5646)
4214.0136 (0.0758)1.9683 (0.5087)1.0997 (0.1847)
Table 2. Estimates for the parameters of the Slomax regression model (n = 50).
Table 2. Estimates for the parameters of the Slomax regression model (n = 50).
Actual ValuesEstimated Values (Standard Deviations)
α β λ γ 1 γ 2 α ^ β ^ λ ^ γ 1 ^ γ 2 ^
221142.31642.40011.38821.02993.9981
(1.0269)(1.0390)(0.5322)(0.4223)(0.4285)
111141.05781.19371.16430.94454.0047
(0.3028)(0.6507)(0.1921)(0.4653)(0.4501)
140.5140.85515.36151.11270.84413.9263
(0.2407)(3.3273)(0.1690)(1.9800)(2.2611)
0.530.5140.48434.62781.07620.84083.3788
0.09832.97160.13683.19682.2944
221−232.17302.30551.2889−2.01052.9500
(0.7806)(0.9017)(0.3916)(0.3225)(0.3497)
111−231.05841.21411.1524−1.95422.9697
(0.3037)(0.6631)(0.1387)(0.3959)(0.4510)
140.5−230.83385.26901.1026−1.93392.5350
(0.2292)(2.7569)(0.1260)(2.5031)(2.5080)
0.530.5−230.73873.25420.5982−1.89972.7328
(0.3441)(0.5349)(0.2066)(0.3553)(1.4824)
221−3−52.25372.34811.1121−3.2708−5.1619
(0.3472)(0.3213)(0.2769)(0.7314)(0.5327)
111−3−51.04101.12621.0871−2.8100−4.9051
(0.4803)(0.3203)(0.1356)(0.3583)(0.3044)
140.5−3−50.73343.25080.3213−2.5547−4.3723
(0.3510)(1.4735)(0.3225)(1.6271)(1.5240)
0.530.5−3−50.42872.85130.8711−3.4246−4.7502
(0.1308)(0.2210)(0.0805)(1.3251)(1.1469)
Table 3. Estimates for the parameters of the Slomax regression model (n = 100).
Table 3. Estimates for the parameters of the Slomax regression model (n = 100).
Actual ValuesEstimated Values (SD)
α β λ γ 1 γ 2 α ^ β ^ λ ^ γ 1 ^ γ 2 ^
221141.99192.12691.26820.99064.0009
(0.5305)(0.4709)(0.2391)(0.2311)(0.1808)
111141.05711.14891.14731.02134.0196
(0.2517)(0.4530)(0.1452)(0.2380)(0.2790)
140.5141.15835.09871.07340.89923.9050
(0.2011)(2.4680)(0.1343)(1.5287)(1.7334)
0.530.5140.46394.18231.09461.12633.7404
(0.0775)(1.6672)(0.1658)(1.4425)(1.3554)
221−232.04662.16321.2692−1.97913.0303
(0.5222)(0.4559)(0.2689)(0.1964)(0.2036)
111−231.03461.13971.1739−2.02672.9309
(0.2461)(0.5273)(0.1933)(0.2701)(0.2840)
140.5−230.83235.05061.0024−1.95082.5885
(0.1907)(2.0433)(0.1102)(1.5322)(1.4601)
0.530.5−230.7210(3.18150.5733−1.94632.7575)
(0.3066)(0.4831)(0.2032)(0.3026)(0.8780)
221−3−52.23342.20871.0897−2.8324−5.1021
(0.2031)(0.1894)(0.2003)(0.5485)(0.4097)
111−3−51.03271.18920.9358−2.9350−5.1317
(0.2566)(0.2387)(0.1348)(0.3369)(0.2447)
140.5−3−50.25133.70440.4024−2.7307−5.2522
(0.3310)(0.8761)(0.2550)(0.5103)(0.7322)
0.530.5−3−50.45532.90320.5041−3.3221−4.8077
(0.1003)(0.1757)(0.0533)(0.8258)(0.6045)
Table 4. Estimates for the parameters of the Slomax regression model (n = 300).
Table 4. Estimates for the parameters of the Slomax regression model (n = 300).
Actual ValuesEstimated Values (SD)
α β λ γ 1 γ 2 α ^ β ^ λ ^ γ 1 ^ γ 2 ^
221141.89192.00561.20361.01053.9960
(0.2624)(0.2031)(0.1497)(0.0966)(0.0889)
111141.00821.07881.17370.97643.9674
(0.1296)(0.2563)(0.1687)(0.1455)(0.1293)
140.5141.08584.73330.88260.94554.1321
(0.1582)(0.9974)(0.1125)(0.9275)(0.7030)
0.530.5140.48703.73240.87071.09733.8105
(0.0532)(0.7099)(0.0984)(1.1261)(0.8640)
221−231.89852.01021.2102−2.00193.0049
(0.3199)(0.3041)(0.1904)(0.0987)(0.0969)
111−230.99451.07231.1833−1.99642.9891
(0.1063)(0.1932)(0.2049)(0.1352)(0.1334)
140.5−230.69044.05330.7347−1.97612.7032
(0.0997)(0.8376)(0.0823)(0.8573)(0.7067)
0.530.5−230.61473.09040.4930−2.07032.0988
(0.1441)(0.2095)(0.2011)(0.2887)(0.5330)
221−3−52.10602.03131.0049−3.1371−5.0980
(0.1355)(0.1208)(0.1003)(0.2487)(0.1072)
111−3−51.03131.05110.9554−3.0508−5.0944
(0.1089)(0.1344)(0.0732)(0.1457)(0.2330)
140.5−3−51.14174.15300.5224−3.1811−5.2004
(0.2113)(0.3499)(0.2057)(0.3734)(0.2810)
0.530.5−3−50.51063.08330.4988−2.8740−5.1311
(0.0773)(0.0927)(0.0471)(0.3273)(0.2474)
Table 5. Parameter estimates and log-likelihood values for different models. Gulomax, Gumbel Lomax; Blomax, beta Lomax; ESlomax, exponentiated standard Lomax.
Table 5. Parameter estimates and log-likelihood values for different models. Gulomax, Gumbel Lomax; Blomax, beta Lomax; ESlomax, exponentiated standard Lomax.
Distribution α β λ σ a η loglikAICBIC
Slomax4.412251.7322.602 979.8261965.6531975.215
Gulomax5.80432.4263.4083.310 981.0741970.1471982.897
Blomax102.216319.152 1.6940.046979.0421966.0851978.834
ESlomax0.788 12.283 1001.5311999.0621992.688
Table 6. Parameter estimates and log-likelihood values for different models.
Table 6. Parameter estimates and log-likelihood values for different models.
Model
ParameterRSlomaxRGulomaxRBlomaxRESlomax
α 0.53405.23350.68190.6235
β 7.08807.683510.3537-
λ 4.57218.6872--
σ -3.8718--
a--7.28675.8261
η --2.2743-
γ 0 −297.3180−7.64121.14145.2382
γ 1 0.11640.01420.03270.0688
γ 2 −0.71480.2185−0.42581.1429
γ 3 2.16720.0843−0.02590.0532
γ 4 0.1286−0.0820−0.0674−0.2133
-loglik504.2635508.7053507.0078526.2180
AIC1024.52701035.41101032.01601066.4360
BIC1045.04201058.49001055.09501084.3860
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Tian, W. Slashed Lomax Distribution and Regression Model. Symmetry 2020, 12, 1877. https://doi.org/10.3390/sym12111877

AMA Style

Li H, Tian W. Slashed Lomax Distribution and Regression Model. Symmetry. 2020; 12(11):1877. https://doi.org/10.3390/sym12111877

Chicago/Turabian Style

Li, Huihui, and Weizhong Tian. 2020. "Slashed Lomax Distribution and Regression Model" Symmetry 12, no. 11: 1877. https://doi.org/10.3390/sym12111877

APA Style

Li, H., & Tian, W. (2020). Slashed Lomax Distribution and Regression Model. Symmetry, 12(11), 1877. https://doi.org/10.3390/sym12111877

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