Reliability Data - Logistic Model

This example appears in the Reliability Growth and Repairable System Analysis Reference book.

Using the reliability growth data given in the table below, do the following:


 * 1)	Find a Gompertz curve that represents the data and plot it with the raw data.


 * 2)	Find a Logistic reliability growth curve that represents the data and plot it with the raw data.

Solution


 * 1) The figure below shows the entered data and the estimated parameters using the standard Gompertz model.




 * Therefore:


 * $$\begin{align}

& \widehat{a}= & 0.9999 \\ & \widehat{b}= & 0.2485 \\ & \widehat{c}= & 0.6858 \end{align}\,\!$$


 * $$R=(0.9999){{(0.2485)}^}\,\!$$


 * The values of the predicted reliabilities are plotted in the figure below.




 * Notice how the standard Gompertz model is not really capable of handling the S-shaped characteristics of this data.


 * 2. The least squares estimators of the Logistic growth curve parameters are given by Crow [9]:


 * where:


 * $$\begin{align}

{{{\hat{b}}}_{1}}= & \frac{\underset{i=0}{\overset{N-1}{\mathop{\sum }}}\,{{T}_{i}}{{Y}_{i}}-N\cdot \bar{T}\cdot \bar{Y}}{\underset{i=0}{\overset{N-1}{\mathop{\sum }}}\,T_{i}^{2}-N\cdot {{{\bar{T}}}^{2}}} \\ \\ {{{\hat{b}}}_{0}}= & \bar{Y}-{{{\hat{b}}}_{1}}\bar{T} \\ \\ {{Y}_{i}}= & \ln \left( \frac{1}-1 \right) \\ \\ \bar{Y}= & \frac{1}{N}\underset{i=0}{\overset{N-1}{\mathop \sum }}\,{{Y}_{i}} \end{align}\,\!$$

In this example $$N=9\,\!$$, which gives:


 * $$\begin{align}

\overline{Y}&=\frac{1}{9}\sum_{i=0}^{8}ln\left (\frac{1}{R_{i}}-1 \right ) \\ &= -1.7355 \\ \\ \overline{T}&=\frac{1}{9}\sum_{i=0}^{8}T_{i} = 4 \\ \sum_{i=0}^{8}T_{i}^{2} &= 204 \\ \sum_{i=0}^{8}T_{i}Y_{i} &= -106.8630 \end{align}\,\!$$

From the equations for $$b_{i}\,\!$$ and $$\hat{b_{0}}\,\!$$:


 * $$\begin{align}

\hat{b_{1}} &= \frac{-106.8630 - 9(4)(-1.7355)}{204-9(4)_{2}} \\ & = 0.7398 \\ \hat{b_{0}} &= -1.7355 - (-0.7398)(4)\\ &= 1.2235 \end{align}\,\!$$

And from the least squares estimators for $$\hat{b}\,\!$$ and $$\hat{k}\,\!$$:


 * $$\begin{align}

\widehat{b}= & {{e}^{1.2235}} \\ = & 3.3991 \\  \widehat{k}= & -(-0.7398) \\ = & 0.7398 \end{align}\,\!$$

Therefore, the Logistic reliability growth curve that represents this data set is given by:


 * $$R=\frac{1}{1+3.3991\,{{e}^{-0.7398\,T}}}\,\!$$

The following figure shows the Reliability vs. Time plot. The plot shows that the observed data set is estimated well by the Logistic reliability growth curve, except in the region closely surrounding the inflection point of the observed reliability. This problem can be overcome by using the modified Gompertz model.