The Mixed Weibull Distribution

Introduction
Besides the Weibull, exponential, normal and lognormal, there are other distributions that are used to model reliability and life data. However, these four represent the most prominent distributions in Weibull++. In the following chapters, we will discuss other distributions that are used under special circumstances: the mixed Weibull, the generalized gamma, the Gumbel, the logistic and the loglogistic distributions. In this chapter, we will introduce mixed Weibull distribution. Other distributions will be discussed in later chapters.

The mixed Weibull distribution (also known as a multimodal Weibull) is used to model data that do not fall on a straight line on a Weibull probability plot. Data of this type, particularly if the data points follow an S-shape on the probability plot, may be indicative of more than one failure mode at work in the population of failure times. Field data from a given mixed population may frequently represent multiple failure modes. The necessity of determining the life regions where these failure modes occur is apparent when it is realized that the times-to-failure for each mode may follow a distinct Weibull distribution, thus requiring individual mathematical treatment. Another reason is that each failure mode may require a different design change to improve the component's reliability [19].

A decreasing failure rate is usually encountered during the early life period of components when the substandard components fail and are removed from the population. The failure rate continues to decrease until all such substandard components fail and are removed. This corresponds to a decreasing failure rate. The Weibull distribution having $$\beta <1$$  is often used to depict this life characteristic.

A second type of failure prevails when the components fail by chance alone and their failure rate is nearly constant. This can be caused by sudden, unpredictable stress applications that have a stress level above those to which the product is designed. Such failures tend to occur throughout the life of a component. The distributions most often used to describe this failure rate characteristic are the exponential distribution and the Weibull distribution with $$\beta \approx 1$$.

A third type of failure is characterized by a failure rate that increases as operating hours are accumulated. Usually, wear has started to set in and this brings the component's performance out of specification. As age increases further, this wear-out process removes more and more components until all components fail. The normal distribution and the Weibull distribution with a $$\beta >1$$  have been successfully used to model the times-to-failure distribution during the wear-out period.

Several different failure modes may occur during the various life periods. A methodology is needed to identify these failure modes and determine their failure distributions and reliabilities. This section presents a procedure whereby the proportion of units failing in each mode is determined and their contribution to the reliability of the component is quantified. From this reliability expression, the remaining major reliability functions, the probability density, the failure rate and the conditional-reliability functions are calculated to complete the reliability analysis of such mixed populations.

Statistical Background
Consider a life test of identical components. The components were placed in a test at age  $$t=0$$  and were tested to failure, with their times-to-failure recorded. Further assume that the test covered the entire lifespan of the units, and different failure modes were observed over each region of life, namely early life (early failure mode), chance life (chance failure mode), and wear-out life (wear-out failure mode). Also, as items failed during the test, they were removed from the test, inspected and segregated into lots according to their failure mode. At the conclusion of the test, there will be $$n$$  subpopulations of  $${{N}_{1}},{{N}_{2}},{{N}_{3}},...,{{N}_{n}}$$  failed components. If the events of the test are now reconstructed, it may be theorized that at age  $$t=0$$  there were actually  $$n$$  separate subpopulations in the test, each with a different times-to-failure distribution and failure mode, even though at   $$t=0$$  the subpopulations were not physically distinguishable. The mixed Weibull methodology accomplishes this segregation based on the results of the life test.

If $$N$$  identical components from a mixed population undertake a mission of  $$t$$  duration, starting the mission at age zero, then the number of components surviving this mission can be found from the following definition of reliability:


 * $${{R}_{1,2,...,n}}(t)=\frac{{{N}_{1,2,3,..,{{n}_{S}}}}(t)}{N}$$

Then:


 * $$\begin{align}

{{N}_{1,2,...,{{n}_{S}}}}(t)= & N[{{R}_{1,2,...,n}}(t)] \\ \\  {{N}_}(t)=& {{N}_{1}}{{R}_{1}}(t);{{N}_}(t)={{N}_{2}}{{R}_{2}}(t) \\ {{N}_}(t)=& {{N}_{3}}{{R}_{3}}(t);...;{{N}_}={{N}_{n}}{{R}_{n}}(t) \end{align}$$

The total number surviving by age $$t$$  in the mixed population is the sum of the number surviving in all subpopulations or:


 * $${{N}_{1,2,...,{{n}_{S}}}}(t)={{N}_}(t)+{{N}_}(t)+{{N}_}(t)+\cdots +{{N}_}(t)$$

Substituting into the reliability equation yields:


 * $${{R}_{1,2,...,n}}(t)=\frac{1}{N}[{{N}_{1}}{{R}_{1}}(t)+{{N}_{2}}{{R}_{2}}(t)+{{N}_{3}}{{R}_{3}}(t)+\cdots +{{N}_{n}}{{R}_{n}}(t)]$$

or:


 * $${{R}_{1,2,...,n}}(t)=\frac{N}{{R}_{1}}(t)+\frac{N}{{R}_{2}}(t)+\frac{N}{{R}_{3}}(t)+\cdots +\frac{N}{{R}_{n}}(t)$$

This expression can also be derived by applying Bayes's theorem [20], which says that the reliability of a component drawn at random from a mixed population composed of $$n$$  types of failure subpopulations is its reliability,  $${{R}_{1}}(t)$$, given that the component is from subpopulation 1, or  $$\tfrac{N}$$  plus its reliability,  $${{R}_{2}}(t)$$ , given that the component is from subpopulation 2, or  $$\tfrac{N}$$  plus its reliability,  $${{R}_{3}}(t)$$ , given that the component is from subpopulation 3, or  $$\tfrac{N}$$ , and so on, plus its reliability,  $${{R}_{n}}(t)$$ , given that the component is from subpopulation  $$n$$ , or  $$\tfrac{N}$$ , and:


 * $$\underset{i=1}{\overset{n}{\mathop \sum }}\,\frac{N}=1$$

This may be written mathematically as:


 * $${{R}_{1,2,...,n}}(t)=\frac{N}{{R}_{1}}(t)+\frac{N}{{R}_{2}}(t)+\frac{N}{{R}_{3}}(t)+\cdots +\frac{N}{{R}_{n}}(t)$$

Other functions of reliability engineering interest are found by applying the fundamentals to the above reliability equation. For example, the probability density function can be found from:


 * $$\begin{align}

{{f}_{1,2,...,n}}(t)= & -\frac{d}{dT}[{{R}_{1,2,...,n}}(t)] \\ {{f}_{1,2,...,n}}(t)= & \frac{N}\left( -\frac{d}{dT}[{{R}_{1}}(t)] \right)+\frac{N}\left( -\frac{d}{dT}[{{R}_{2}}(t)] \right) \\ & +\ \ \frac{N}\left( -\frac{d}{dT}[{{R}_{3}}(t)] \right)+\cdots +\frac{N}\left( -\frac{d}{dT}[{{R}_{n}}(t)] \right) \\ {{f}_{1,2,...,n}}(t)= & \frac{N}{{f}_{1}}(t)+\frac{N}{{f}_{2}}(t) \\ & +\ \ \frac{N}{{f}_{3}}(t)+\cdots +\frac{N}{{f}_{n}}(t) \end{align}$$

Also, the failure rate function of a population is given by:


 * $$\begin{align}

{{\lambda }_{1,2,...,n}}(t)= & \frac{{{f}_{1,2,...,n}}(t)}{{{R}_{1,2,...,n}}(t)}, \\ {{\lambda }_{1,2,...,n}}(t)= & \frac{\tfrac{N}{{f}_{1}}(t)+\tfrac{N}{{f}_{2}}(t)+\tfrac{N}{{f}_{3}}(t)+\cdots +\tfrac{N}{{f}_{n}}(t)}{\tfrac{N}{{R}_{1}}(t)+\tfrac{N}{{R}_{2}}(t)+\tfrac{N}{{R}_{3}}(t)+\cdots +\tfrac{N}{{R}_{n}}(t)}. \end{align}$$

The conditional reliability for a new mission of duration $$t$$, starting this mission at age  $$T$$ , or after having already operated a total of  $$T$$  hours, is given by:


 * $$\begin{align}

{{R}_{1,2,...,n}}(T,t)= & \frac{{{R}_{1,2,...,n}}(T+t)}{{{R}_{1,2,...,n}}(T)} \\ {{R}_{1,2,...,n}}(T,t)= & \frac{\tfrac{N}{{R}_{1}}(T+t)+\tfrac{N}{{R}_{2}}(T+t)+\cdots +\tfrac{N}{{R}_{n}}(T+t)}{\tfrac{N}{{R}_{1}}(T)+\tfrac{N}{{R}_{2}}(T)+\cdots +\tfrac{N}{{R}_{n}}(T)} \end{align}$$

The Mixed Weibull Equations
Depending on the number of subpopulations chosen, Weibull++ uses the following equations for the reliability and probability density functions:


 * $${{R}_{1,...,S}}(t)=\underset{i=1}{\overset{S}{\mathop \sum }}\,\frac{N}{{e}^{-{{\left( \tfrac{t}{{{\eta }_{i}}} \right)}^}}}$$

and:


 * $${{f}_{1,...,S}}(t)=\underset{i=1}{\overset{S}{\mathop \sum }}\,\frac{N{{\eta }_{i}}}{{\left( \frac{t}{{{\eta }_{i}}} \right)}^{{{\beta }_{i}}-1}}{{e}^{-{{(\tfrac{t}{{{\eta }_{i}}})}^}}}$$

where $$S=2$$,  $$S=3$$ , and  $$S=4$$  for 2, 3 and 4 subpopulations respectively. Weibull++ uses a non-linear regression method or direct maximum likelihood methods to estimate the parameters.

Regression Solution
Weibull++ utilizes a modified Levenberg-Marquardt algorithm (non-linear regression) when performing regression analysis on a mixed Weibull distribution. The procedure is rather involved and is beyond the scope of this reference. It is sufficient to say that the algorithm fits a curved line of the form:


 * $${{R}_{1,...,S}}(t)=\underset{i=1}{\overset{S}{\mathop \sum }}\,{{\rho }_{i}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{i}}} \right)}^}}}$$

where:


 * $$\underset{i=1}{\overset{S}{\mathop \sum }}\,{{\rho }_{i}}=1$$

to the parameters $$\widehat$$   $$\widehat,$$   $$\widehat,$$   $$\widehat\widehat,$$   $$\widehat,...,$$   $$\widehat{{{\rho }_{S,}}\text{ }}\widehat,$$   $$\widehat,$$  utilizing the times-to-failure and their respective plotting positions. It is important to note that in the case of regression analysis, using a mixed Weibull model, the choice of regression axis, (i.e., $$RRX$$  or  $$RRY,$$  is of no consequence since non-linear regression is utilized).

MLE
The same space of parameters, namely $$\widehat$$   $$\widehat,$$   $$\widehat,$$   $$\widehat\widehat,$$   $$\widehat,...,$$   $$\widehat{{{\rho }_{S,}}\text{ }}\widehat,$$   $$\widehat,$$  is also used under the MLE method, using the likelihood function as given in Appendix C of this reference. Weibull++ uses the EM algorithm, short for Expectation-Maximization algorithm, for the MLE analysis. Details on the numerical procedure are beyond the scope of this reference.

Mixed Weibull Confidence Bounds
In Weibull++, two methods are available for estimating the confidence bounds for the mixed Weibull distribution. The first method is the beta binomial, described in Confidence Bounds. The second method is the Fisher matrix confidence bounds. For the Fisher matrix bounds, the methodology is the same as described in Confidence Bounds. The variance/covariance matrix for the mixed Weibull is a $$(3\cdot S-1)\times (3\cdot S-1)$$  matrix, where  $$S$$  is the number of subpopulations. Bounds on the parameters, reliability and time are estimated using the same transformations and methods that were used for the The Weibull Distribution. Note, however, that in addition to the Weibull parameters, the bounds on the subpopulation portions are obtained as well. The bounds on the portions are estimated by:


 * $$\begin{align}

& {{\rho }_{U}}= & \frac{\hat{\rho }+(1-\hat{\rho }){{e}^{-\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\rho })}}{\hat{\rho }(1-\hat{\rho })}}}} \\ & &  \\  & {{\rho }_{L}}= & \frac{\hat{\rho }+(1-\hat{\rho }){{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\rho })}}{\hat{\rho }(1-\hat{\rho })}}}} \end{align}$$

where $$Var(\widehat{\rho })$$  is obtained from the variance/covariance matrix. When using the Fisher matrix bounds method, problems can occur on the transition points of the distribution, and in particular on the Type 1 confidence bounds (bounds on time). The problems (i.e.. the departure from the expected monotonic behavior) occur when the transition region between two subpopulations becomes a "saddle" (i.e., the probability line is almost parallel to the time axis on a probability plot). In this case, the bounds on time approach infinity. This behavior is more frequently encountered with smaller sample sizes. The physical interpretation is that there is insufficient data to support any inferences when in this region.

This is graphically illustrated in the following figure. In this plot it can be seen that there are no data points between the last point of the first subpopulation and the first point of the second subpopulation, thus the uncertainty is high, as described by the mathematical model.



Beta binomial bounds can be used instead in these cases, especially when estimations are to be obtained close to these regions.

Reliability Bathtub Curves
A reliability bathtub curve is nothing more than the graph of the failure rate versus time, over the life of the product. In general, the life stages of the product consist of early, chance and wear-out. Weibull++ allows you to plot this by simply selecting the failure rate plot, as shown next.



Determination of the Burn-in Period
If the failure rate goal is known, then the burn-in period can be found from the failure rate plot by drawing a horizontal line at the failure rate goal level and then finding the intersection with the failure rate curve. Next, drop vertically at the intersection, and read off the burn-in time from the time axis. This burn-in time helps insure that the population will have a failure rate that is at least equal to or lower than the goal after the burn-in period. The same could also be obtained using the Function Wizard and generating different failure rates based on time increments. Using these generated times and the corresponding failure rates, one can decide on the optimum burn-in time versus the corresponding desired failure rate.

Using the Mixed Weibull Distribution in Weibull++
To use the mixed Weibull distribution, simply select the Mixed Weibull Distribution from the Distribution drop down list. You can choose Mixed Weibull distribution with 2, 3, or 4 subpopulations?



Viewing the Calculated Parameters

When using the Mixed Weibull option, the parameters given in the result area apply to different subpopulations. To view the results for a particular subpopulation, select the subpopulation, as shown next.





About the Calculated Parameters

Weibull++ uses the numbers 1, 2, 3 and 4 (or first, second, third and fourth subpopulation) to identify each subpopulation. These are just designations for each subpopulation, and they are ordered based on the value of the scale parameter, $$\eta $$. Since the equation used is additive or:


 * $${{R}_{1,..,S}}(T)=\underset{i=1}{\overset{S}{\mathop \sum }}\,\frac{N}{{e}^{-{{\left( \tfrac{t}{{{\eta }_{i}}} \right)}^}}}$$

the order of the subpopulations which are given the designation 1, 2, 3, or 4 is of no consequence. For consistency, the application will always return the order of the results based on the magnitude of the scale parameter.

Example 1: