Various probability density functions have been proposed to perform statistical analysis of lifetime data. The Weibull distribution is one of the most widely used distributions in the analysis of lifetimes data. It was introduced by the French Mathematicians Fréchet (1928)  . Indeed in the 1920s Fréchet developed a distribution to which he gave his name; Fréchet distribution, as an extreme value distribution. This distribution is in fact equal to the reciprocal of the Weibull distribution. Rosin and Rammler (1933)  applied Fréchet’s ditribution to describe the particle size distribution generated by grinding, milling and crushing operations of materials. This probability distribution has been widely used as a probabilistic model in studies on lifetimes. Mudholkar and Srivastava (1993)  introduced the exponentiated Weibull distribution to analyse bathtub failure rate data which cannot be handled well by the regular Weibull for monotonicity of its hazard rate. Also Zhang and Xie (2011)  worked on bathtub failure data using the truncated Weibull distribution. Soumaya and Soufiane (2014)  have given estimation of the parameters of the exponentiated Weibull distribution and the additive Weibull distribution, which are two specific generalizations of the Weibull distribution.
Cordeiro, et al. (2013)  introduced the exponentiated generalized class of distribution which is more general than the two classes of Lehmann’s (1953)  alternatives, it is a combination of the Lehmann type I and type II alternatives. Indeed, for any baseline (or parent) distribution it is possible to define the corresponding Exponentiated Generalized family of distribution. Cordeiro, et al. (2013) discussed four special models namely the Exponentiated Generalized Fréchet, the Exponentiated Generalized Normal distribution, the Exponentiated Generalized Gamma distribution and the Exponentiated Generalized Gumbel distribution. Oguntunde et al. (2015)  have discussed the special case of the Exponentiated Generalized Weibull distribution by using the Weibull distribution as baseline distribution. The proposed distribution has four parameters (three shape parameters and one scale parameter). The work of Oguntunde et al. is mainly focused on the mathematical properties of the distribution like the moments, the limiting behaviour of the functions (pdf and cdf), the reliability analysis, and the quantile function.
The probability density function and the cumulative distribution function of the Exponentiated Generalized Weibull are respectively given by:
where , , , , .
The Exponentiated Generalized Weibull generalizes the following distributions:
For , Generalized Weibull;
For , Exponentiated Weibull;
For , Weibull distribution;
For , Exponential distribution.
The survival function and the hazard function have respectively the following expressions:
2. Parameters Estimation
2.1. The Model
Let us assume that we have n independent variables in a trial, and the ordered m failures are observed under the progressive type-II censoring plan , where for and . Let the observed and censored data be respectively and , where for . Now consider to be the complete data (observed and censored data together). Then the joint probability that the complete sample (the complete data likelihood) is observed is given by
(Ng et al 2002)  .
From which we get the following log-likelihood by substituting in (5) the pdf by expression (1)
2.2. EM Algorithm
In oder to tackle the E-step the conditional expectation of the log-likelihood given the observed sample is computed. Let us denote it by where is the vector of parameters. That is,
The conditional expectation of the above log-likelihood becomes
Thus, to facilitate the E-step, the conditional distribution of Z for given Y and the current value of the parameters, needs to be determined.
The conditional distribution of Z for given Y is given by
see Ng et al. (2002).
Let us set
Using (8)-(11) the expressions for , and become
is the upper incomplete gamma function.
We therefore obtain an expression for the conditional expectation of the log- likelihood as
where the functions A, B, and C are respectively defined in (12)-(14).
In the M-step on the p-th iteration of the EM algorithm, the value of which maximizes will be used as the next estimate of . Where is the vector of parameters at the p-th iteration , and is the initial value of the vector of parameters.
Therefore, if at the p-th stage the estimate of is , then can be obtained by maximizing
Then is solution of the following system of equations
which is equivalent to
From the second equation in the above system we can express for known , , as:
The expressions for , and are not available in closed form.
The solution to the M-step does not exists in closed form. For this case Dempster et al. (1977)  defined what is called the generalized EM algorithm (GEM algorithm) for which the M-step requires to be chosen such that
Since we need only to increase the likelihood, we may replace the M-step with a single iteration of the Newton-Raphson (N-R) algorithm.
For the simulation the M-step is replaced by a single iteration of the Newton- Raphson algorithm.
For the values of and and progressively Type- II censored sample was generated from the Exponentiated Generalized Weibull distribution using the algorithm in Balakrishnan and Sandhu (1995)  .
The algorithm is defined as follows:
• Generate m independent Uniform (0, 1) observations
• Set for
• Set for . Then is the required progressive Type-II censored sample from the Uniform (0, 1) distribution.
• Finally, we set for , where is the inverse cdf of the Exponentiated Generalized Weibull distribution. Then is the required progressive Type-II censored sample from the Exponentiated Generalized Weibull distribution.
with censoring scheme .
The generated sample is:
0.0138 0.0230 0.0447 0.2401 0.3091 0.3264 0.4597 0.5448 0.5841 0.7274 0.9875 1.1164 1.2090 1.3519 1.4896 1.5041 1.6224 2.9952 3.4537 3.6385
Via the EM algorithm discussed in Section 2, the computed MLEs of the parameters become:
, , and
In Table 1 a Monte Carlo simulation for N = 500 was used to compute the RMSE and the mean estimates for different value of n, m and . The following formula was used to compute the RMSE
where is the i-th estimates of the parameter
• For fixed sample size n and by increasing m, we get smaller RMSE’s.
• By increasing the sample size n, we get smaller RMSE’s.
• The largest values of m in each case represent the complete sample case.
Table 1. RMSE of the estimators.
The parameters of the Exponentiated Generalized Weibull distribution were estimated using maximum likelihood estimation method via Expectation Maximization (EM) algorithm. The Root Mean Square Error were computed at different values of the sample size n and failures (observed data) m. It was observed that the RMSEs were smaller for fixed sample size n and increasing the size m of the observed data, and also for the increasing sample size n.
I would like to thank my supervisors Professor Leo Odongo and Doctor Ibrahim Ly for accompanying me through this work. Sincere thanks to the African Union for giving me the opportunity to do scientific reseach.