In statistical analysis, especially for the lifetime data the model with monotone risk functions such as gamma and Weibull distributions is preferred. Recently modified extensions of the exponential distributions have monotone risk functions to be suitable for life time data analysis. Gupta and Kundu [1, 2] proposed and discussed the generalized exponential distribution [GED] under probability density function [PDF] given by (1)
Recently, Nadarajah and Haghighi  proposed also another extension of the exponential distribution with pdf (2)
Recently, Gemez et al.  proposed another extension of the exponential distribution call extension exponential distribution [EED], which is adopted to studying the coefficient of variation. More recently, Makhdoom and Nasiri  estimated the exponential parameter under the type-II censoring scheme when the data observed is a fuzzy sample. The generalized extended exponential distribution is presented by Thiago et. al. . Also - an extended exponential geometric distribution applied for medical data see, Louzada .
Let random variable X has an EED with parameters α and β. The PDF and CDF of EED is given respectively by (3) and (4)
Also the reliability and hazard rate functions given respectively by (5) and (6)
The EED satisfies the following properties see (Gemez et al. ).
The population mean μ of X (7)
The population standard deviation σ of X (8)
The population CV is defined as the ratio of the standard deviation to the population mean. The CV used to measure the stability or uncertainty, also the relative dispersion of population data to the population mean can be presented. Different authors addressed the problems of CV (Pang et al. ). For applications of CV in, medical sciences, engineering, physics, chemistry and telecommunications see, Miller and Karson , Hamer et al. , Reh and Scheffler , Ahn  and Gong and Li . Several authors studied CV in normal phenomena but the works in non normal phenomena are rare, see Soliman et al.  and Abd-Elmougod et al. .
The main aim of this paper iss to study CV of EED when the data obtained from type-II censoring experiment and the joint likelihood function as the form (9) where we will execute maximum likelihood estimation with the an extension of the exponential distribution. Also the parametric bootstrap confidence intervals are calculated. The posterior empirical distribution of CV with the help of MCMC method are obtained. The obtaining results are discussed with real data set.
This article is organized as follows. In Section 2, a maximum likelihood estimation is presented. In Section 3, parametric bootstrap confidence intervals are discussed. In Section 4, the point and interval estimation of CV under posterior empirical distribution are discussed. A real life data set applied to illustrate our proposed method is presented in Section 5. Conclusion can be found in Section 6.
2 Maximum likelihood estimation
Without any constant the log L of the likelihood function (10) can be presented by (11)
The likelihood equations obtained by taking the first derivatives of (11) with respect to parameters α and β, hence equating each of them to zero, we obtain (12) and (13)
The equations (12) and (13) cannot be solved analytically for , these are/can be solved numerical with the Newton’s method with the initial parameters values obtained by the method of moments as (14) and (15)
3 Bootstrap confidence intervals
To obtain interval estimations of the parameters or any function of the parameters, we can use the Bootstrap technique which builds confidence intervals or widely improves estimators other than asymptotic results. Parametric and nonparametric bootstrap methods, are discussed in  and . More details about parametric bootstrap percentile technique which is used in this section to build a bootstrap confidence interval of CV are also discussed in . To obtain for estimation of CV of EED we consider the following algorithm under type-II bootstrap sample.
The estimate values of are used for generation a bootstrap data x⋆ under the fixed values of n, and r.
Based on x⋆ as given in step 1, we built the bootstrap samples estimation of and CVML say are computed.
Steps 2 and 3 are repeated s times, the s bootstrap samples of based on s different data x⋆.
The ordering bootstrap samples are obtained by arrange all in an ascending order.
The bootstrap empirical distribution of say Ψ(t) = for given t. Then the corresponding bootstrap 100(1 − 2 γ)% confidence interval of given by (17)
4 Bayesian estimation
Bayes estimation of CV of EED under squared error loss [SEL] function can be achieved under independent gamma prior for each parameters α and β. So that the prior density of α and β can be represented as (18) and (19)
The joint posterior density function [JPDF] h(α,β |x) respected to parameters α and β for given data represented as (20) under SEL, the estimation of any function φ (α, β) say [α, β or CV], given by (21)
In several cases the ratio in (21) of two integrals does not have a closed form. So numerical integration or Lindely approximation can be used for this problem. Another important method, can be used for generating samples from the posterior distributions called MCMC methods, which present estimation of φ(α, β) under the SEL.
From equations (20) and (21), it is not possible to compute analytically form in this case. Different approaches can be used here, the important and more accurate method is Markov Chain Monte Carlo (MCMC) which can be used not only to approximate the point estimation of parameters and any function of the parameters such as a CV. It has the advantage of interval estimation of parameters and any function of the parameters such as a CV. From equations (10), (18), (19) and (20) the JPDF of α and β given data can be obtained as (22)
The full conditional distribution of α and β written by (23) and (24)
A more variety of MCMC schemes are existing. The important one of MCMC methods is Gibbs sampling. The property of two distributions of (23) and (24) shows that each are similar to normal distribution, so Metropolis-within-Gibbs samplers are employed. We use the MH method Metropolis et al.  with normal proposal distribution as follows.
MH under Gibbs sampling algorithm
Put the indicator I = 1 and ML estimates as the intial values of the parameters α and β.
Compute φ(i) (α, β) = φ(α(I), β(I))
Put I = I + 1 and.
Repeat steps from 2 and 3, N⋆ times.
The Bayes estimation of φ(α, β) under the MCMC methods as (25) where M⋆ is the number of iterations, need to the stationary distribution.
1. Put the values φ(1), φ(2), …, φ(N⋆−M⋆) in ordered forms as φ(1), φ(2), …, φ(N⋆−M⋆) so the symmetric credible interval 100(1 − 2γ)% given by (26)
5 Data analysis
Let us consider the real data set presented by, which fits to EED. The data presented the life of fatigue fracture of Kevlar 373/epoxy that are subject to constant pressure at the 90% stress level until all had failed, for more surveys about this data see Andrews and Herzberg  and Barlow et al. . Let us consider that type-II censoring sample with censoring schemes n = 76 and r = 40. The obtained type-II censoring data is presented in Table 1. The classical point estimation of parameters α and β and CV using each of maximum likelihood and bootstrap methods as will as the approximate percentile bootstrap confidence intervals are given in Table 2. In Bayesian approach based on the MCMC samples of size 11000 with 1000 as ‘burn-in’, the Bayes point estimates and 95% credible intervals for α, β and CV are computed and presented in Table 3. Figures (1-6) show the simulation number of α, β and CV generated by MCMC method and the corresponding histogram. The figures (1-6) plots the first 11000 steps for the Gibbs sampler and shown the convergence in the empirical posterior results. The corresponding histogram given by figures (2), (4) and (6) show that the empirical posterior distribution of each parameter and CV. Figures (2) and (4) show that the empirical posterior distribution of parameter α and β is similar to normal distribution, so the posterior mean is used to get a better approximation for the posterior estimation of parameters α and β. Also the histogram in figure (6) it is quite skewed for a skewed distribution, so the posterior mode is adopted as the best estimation for the CV more details on the CV see Upadhyay and Peshwani .
In practice, the theoretical sampling distribution interpreted in CV is more difficult to be represented mathematically. Therefore, the statistical inference about the CV in a several applications need to obtain an interval estimation for the CV. In this article we proposed Bayesian and non-Bayesian approach for estimation results of CV. The proposed method for obtaining a point estimate as well as an interval estimate for the CV. The results show an effective MCMC method therefore our method used in this article can be easily extended for another extension exponential distribution discussed in this article.
This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant no.(G-541-558-1437). The authors, therefore, acknowledge with thanks DSR for technical and financial support.
Nadarajah S., Haghighi F., An extension of the exponential distribution, A Journal of Theoretical and Applied Statistics, 2011, 45, 543. Google Scholar
Yolanda M.G., Heleno B., Gomez H.W., On the performance evaluation of different measures of association, Revista Colombiana de Estadstica Junio, 2014, 37, 25. Google Scholar
Makhdoom I., Nasiri P., Maximum likelihood estimation of exponential distribution under type-II censoring from imprecise data, Journal of Fundamental and Applied Sciences, 2016, 8, 697. CrossrefGoogle Scholar
Thiago A.N., Bourguignon A.N. M., Cordeiro G.M., The exponentiated generalized extended exponential distribution, Journal of Data Science, 2016, 14, 393. Google Scholar
Louzada F., The exponentiated exponential?geometric distribution: a distribution with decreasing, increasing and unimodal failure rate, Revista Colombiana de Estadstica Junio, 2014, 37, 1. Google Scholar
Pang W.K., Leung P.K., Huangand W.K., Wei L., On interval estimation of the coefficient of variation for the three-parameter Weibull, lognormal and gamma distribution: A simulation bas ed approach, European Journal of Operational Research, 2005, 164, 367. CrossrefGoogle Scholar
Miller E.G., Karson M.J., Testing Equality of Two Coefficients of Variation, American statistical Association Posceedings of the Business and economics Section, 1977, 1, 278. Google Scholar
Hamer A.J., Strachan J.R., Black M.M., Ibbotson C., Elson R.A., A new method of comparative bone strength measurement, Journal of Medical Engineering and Technology. 1995, 19, 1. CrossrefGoogle Scholar
Gong J., Li Y., Relationship between the estimated Weibull modulus and the coefficient of variation of the measured strength for ceramics, Journal of the American Ceramics Society, 1999, 82, 449. Google Scholar
Soliman A.A., Abd Ellah A. H., Abou-Elheggag N.A., Abd-Elmougod G.A., A simulation-based approach to the study of coefficient of variation of Gompertz distribution under progressive first-failure censoring, Indian Journal of Pure and Applied Mathematics. 2011, 42, 335. CrossrefWeb of ScienceGoogle Scholar
Abd-Elmougod G.A., El-Sayed M.A., Abdel-Rahman E.O., Estimation for coefficient of variation of Topp-Leone distribution under adaptive Type-II progressive censoring scheme: Bayesian and non-Bayesian approaches, Journal of Computational and Theoretical Nanoscience, 2015, 12, 4028. CrossrefGoogle Scholar
Davison A.C., Hinkley D.V., Bootstrap methods and their applications, Cambridge University Press, Cambridge United Kingdom, 1997. Google Scholar
Efron B., Tibshirani R.J., An introduction to the bootstrap, New York Chapman and Hall, 1993. Google Scholar
Efron B., The jackknife, the bootstrap and other resampling plans, SIAM, Phiadelphia, PA, 1982, 38, 86. Google Scholar
Andrews D.F., Herzberg A. M., Data: a collection of problems from many fields for the student and research worker, Springer Series in Statistics, New York, 1985. Google Scholar
Barlow R.E., Toland R.H., Freeman T., A Bayesian analysis of stress rupture life of kevlar 49/epoxy spherical pressure vessels, in ‘Proc, In Proc. Conference on Applications of Statistics’, Marcel Dekker, New York, 1984. Google Scholar
Upadhyay S.K., Peshwani M., Full posterior analysis of three parameter lognormal distribution using Gibbs sampler, Journal of Statistical Computation and Simulation, 2001, 71, 215. CrossrefGoogle Scholar
About the article
Published Online: 2017-08-10