In this paper, we discuss the mixture model of two extreme lower bound distributions. First, some properties we obtain of the model with hazard function are discussed. In addition, the estimates of the unknown parameters via the EM algorithm are obtained. The performance of the findings in the paper is showed by demonstrating some numerical illustrations through Monte Carlo simulation.
Mixture Extreme Lower Bound Distribution Reliability Estimation EM Algorithm Monte Carlo Simulation1. Introduction
Recently, the extreme value distribution is becoming increasingly important in engineering statistics as a suitable model to represent phenomena with usually large maximum observations. In engineering circles, this distribution is often called the extreme lower bound model. It is one of the pioneers of extreme value statistics. The extreme lower bound distribution is one of the probability distributions used to model extreme events. The generalization of the standard extreme lower bound distribution has been introduced by Nadarajah and Kotz [1] and Abd-Elfattah [2] . There are over fifty applications ranging from accelerated life testing through to earthquakes, floods, rain fall, queues in supermarkets, sea currents, wind speeds and track race records, see Kotz and Nadarajah [3] . Mixture models play an important role in many practical applications. For example, direct applications of finite mixture models are in fisheries research, economics, medicine, psychology, palaeoanthropology, botany, agriculture, zoology, life testing and reliability. Direct applications include outliers, Gaussian sums, cluster analysis, latent structure models, modeling prior densities, empirical Bayes method and nonparametric density estimation. In many applications, the available data can be considered as data coming from a mixture population of two or more distributions. This data enables us to mix statistical distributions to get a new distribution which has the properties of its components. For an excellent survey of estimation techniques, discussion and applications, mixture distribution have been considered extensively by many authors, see Titterington [4] , Maclachlan and Basford [5] , Lindsay [6] , Maclachlan and Krishnan [7] and Maclachlan and Peel [8] . Recently, there are many authors [4] [9] [10] who discuss the mixture models, Mohie El-Din et al. [11] [12] [13] . In this paper, we discuss some important measures of two extreme lower bound distributions. Also, we estimate the vector of unknown parameters of a mixture model via the EM algorithm proposed by Dempster et al. [9] . Further, we carry out some simulated illustrations using Monte Carlo method.
2. Description of the Model
The mixture of two extreme lower bound distributions has its pdf as
where and, the density func- tion of component, is given by
The cdf of the mixture of two extreme lower bound distributions is given by
where, the cdf component, is given by
Such that, We study this case, when and are the parameters unknown and is the parameter known.
3. Properties
In this section we obtain some properties for two extreme lower bound distri- bution by extending the corresponding results of the two parameters extreme lower bound distribution where (2.3) is known as follow.
3.1. The Expected Value and Variance
The expected value of the pdf of the two extreme lower bound distribution obtain in (2.1) and (2.3) is
and the variance is given by
3.2. Mode and Median
The mode of the mixture of two extreme lower bound distribution is obtained by solving the following nonlinear equation with respect to x
By using (2.4) and (3.5), the median of the mixture of two extreme lower bound distribution is obtained by solving the following nonlinear equation with respect to x
From Table 1, we obtain the median and the mode of the mixture two extreme lower bound distribution based on different choices of the parameters and for each from this table we observe that the mode is slightly affected by the variation in the values of the mixing proportion, while one mode is stable in the bimodal case. In addition, for unimodal case, the median increases when increases. From the bimodal case, we observe that the median decreases when icreases.
3.3. Reliability and Failure Rate Function
The reliability function of the mixture two extreme lower bound distribution is given by
The median and the mode of the mixture of two extreme lower bound distribution
Bimodel case
(0.2, 2.5, 2, 1, 2.9)
(0.4, 2.5, 2, 1, 2.9)
(0.6, 2.5, 2, 1, 2.9)
Median
1.0286
0.9201
0.7852
Mode
0.326, 0.885
0.326, 0.863
0.326, 0.798
Unimodel case
(0.2, 1, 2, 2, 3)
(0.4, 1, 2, 2, 3)
(0.6, 1, 2, 2, 3)
Median
0.635
0.733
0.882
Mode
0.459
0.469
0.882
By using (2.4) and (3.5) it can be seen the failure rate function (hazard rate function) of the mixture two extreme lower bound distribution is given by
Which can be written as
where
and
The failure rate function of the mixture two extreme lower bound distribution given in (3. 10) satisfies the following limits
4. Estimation via EM Algorithm
The EM algorithm provides a simple computational method for fitting mixture models. We use the EM algorithm to estimate the parameters of the pdf of the mixture two extreme lower bound distribution which given in (2.1) and (2.3). We focus in this section, the Maximum likelihood fitting of two extreme lower bound distributions mixture via the EM algorithm. Maclachlan and Peel [9] , the essential nature of the algorithm is the alternation of expectation and maximization steps.
then, Concerning the E-step on the iteration, the updated estimate of the mixing proportion is given by
From (4.13) we obtain the M-step of the iteration, the updated estimates and for each are obtained, respectively, by solving the following systems of equations
and
where
The estimates of and are obtained by solving (4.15), (4.17) and (4.18). Equations (4.15) and (4.17) are written explicitly but Equation (4.18) has to be solved numerically with random choices of the initial values.
5. Numerical Illustration
In order to calculate the estimates of the five parameters and where and are known that appear in the pdf of the mixture two extreme lower bound distribution given in (2.1) and (2.3) by using EM algorithm in a Monte Carlo simulation as follows:
Generate random sample of size and 100 from the mixture two extreme lower bound distribution distribution with for each choice of the parameters and. Some of choices caver the unimodal case and other caver the bimodal case.
The random samples of the mixtures are generated with respect to two uniform variables and. If, then use to generate a random variable x from the mixture two extreme lower bound distribution by using (3.5) as, but if, then.
The bias and the mean square errors of the estimates are calculated based on 10,000 Monte Carlo simulation and the results are illustrated in Table 2 and Table 3. We see that in most of the considered cases, the mean square errors of the estimated parameters decrease as n increase.
6. Conclusion
In this paper, the behaviors of the mode and median of the mixture two extreme
Bais of the estimate of based on EM algorithm
Bimodal
(0.2, 2.5, 2, 1, 2.9)
50
−0.116
−1.204
−0.313
0.299
−1.203
100
−0.056
−1.189
−0.333
0.313
−1.242
(0.4, 2.5, 2, 1, 2.9)
50
−0.068
−0.770
−0.411
0.729
−1.309
100
−0.001
−0.757
−0.421
0.743
−1.333
(0.6, 2.5, 2, 1, 2.9)
50
0.006
−0.614
−0.398
0.894
−1.270
100
0.015
−0.611
−0.379
0.886
−1.295
Unimodal
(0.2, 1, 2, 2, 3)
50
−0.042
0.864
0.487
−0.136
−0.514
100
−0.025
0.850
0.423
−0.131
−0.542
(0.4, 1, 2, 2, 3)
50
0.011
0.565
−0.019
−0.483
−1.025
100
0.005
0.569
−0.015
−0.438
−0.101
(0.6, 1, 2, 2, 3)
50
0.024
0.465
−0.089
−0.613
−1.170
100
0.006
0.494
−0.074
−0.573
−1.132
MSE of based on EM algorithm
Unimodal
(0.2, 2.5, 2, 1, 2.9)
50
0.025
1.0002
0.1003
0.098
1.003
100
0.006
1.0001
0.120
0.095
1.0001
(0.4, 2.5, 2, 1, 2.9)
50
0.202
1.127
0.1003
0.200
2.005
100
0.309
1.0007
0.1001
0.199
2.002
(0.6, 2.5, 2, 1, 2.9)
50
0.002
0.377
0.108
0.798
2.002
100
0.0005
0.369
0.1008
0.789
2.000
Unimodal
(0.2, 1, 2, 2, 3)
50
0.0021
0.7002
0.2384
0.0189
0.2005
100
0.0005
0.7003
0.1798
0.0183
0.2004
(0.4, 1, 2, 2, 3)
50
0.0084
0.5987
0.0725
0.0482
0.5000
100
0.0024
0.5169
0.0248
0.0430
0.5023
(0.6, 1, 2, 2, 3)
50
0.0020
0.2145
0.0083
0.3001
1.000
100
0.0002
0.2003
0.0061
0.3000
0.9998
lower bound distribution are investigated, based on different choices of the parameters. Also, the behaviors of the failure rate function are discussed through some different graphs. In addition, the estimation of the unknown parameters is obtained using the EM algorithm. Finally, a Monte Carlo simulation based on 10,000 runs is carried out.
Cite this paper
Riad, F.A. (2017) Inference and Properties of Mixture Two Extreme Lower Bound Distributions. Engineering, 9, 517-523. https://doi.org/10.4236/eng.2017.96032
ReferencesAbd-Elfattah, A.M. and Omima, A.M. (2009) Estimation of the Unknown Parameters of the Generalized Frechet Distribution. Journal of Applied Sciences Research, 5, 1398-1408.Dempster A. P., Laird N. M., Rubin D. B. ,et al. (1977)Maximum-Likelihood from Incomplete Data via the EM Algorithm. Journal of the Royal Statistical Society 39, 1-38.Everitt, B.S. and Hand, D.J. (1981) Finite Mixture Distribution. Chapman and Hall, London. https://doi.org/10.1007/978-94-009-5897-5Jiang, R., Zuo, M.J. and Li, H. (1999) Weibull and Inverse Weibull Mixture Models Allowing Negative Weights. Reliability Engineering & System Safety, 66, 227-234.Kotz, S. and Nadarajah, S. (2000) Extreme Value Distributions: Theory and Applications. Imperial College Press, London. https://doi.org/10.1142/p191Landsay, B.G. (1995) Mixture Models: Theory, Geometry and Applications. The Institute of Mathematical Statistics, Hayward, CA.Maclachlan, G.J. and Basford, K.E. (1988) Mixture Models: Applications to Clustering. Marcel Dekker, New York.Maclachlan, G.J. and Krishnan, T. (1997) The EM Algorithm and Extensions. Wiley, New York.Maclachlan, G.J. and Peel, D. (2000) Finite Mixture Models. Wiley, New York.
https://doi.org/10.1002/0471721182Mohie El-Din, M.M., Abdel-Aty, Y. and Riad, F.H. (2010) Estimations of the Parameters of the Extreme Lower Bound Distribution with Progressively Censored Data. Far East Journal of Theoretical Statistics, 32, 25-34.Mohie El-Din, M.M., Riad, F.H. and El-Sayed, M.A. (2014) Statistical Inference and Prediction for the Inverse Weibull Distribution Based on Record Data. Journal of Statistics Applications and Probability, 3, 171-177.
https://doi.org/10.12785/jsap/030207Mohie El-Din, M.M., Riad, F.H. and El-Sayed, M.A. (2014) Confidence Intervals for Parameters of IWD Based on MLE and Bootstrap. International Journal of Advanced Statistics and Probability, 2, 28-33.Nadarajah, S. and Kotz, S. (2003) The Exponentiated Frechet Distribution. Interstat Lectronic Journal.