Grouped per Configuration Data - Logistic Model

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

Some equipment underwent testing in different stages. The testing may have been performed in subsequent days, weeks or months with an unequal number of units tested every day. Each group was tested and several failures occurred. The data set is given in columns 1 and 2 of the following table. Find the Logistic model that best fits the data, and plot it along with the reliability observed from the raw data.

Solution

The observed reliability is $$1-\tfrac{\#\text{ of failures}}{\#\text{ of units}}\,\!$$ and the last column of the table above shows the values for each group. With $$N=9\,\!$$, the least square estimator $$\overline{Y}$$ becomes:


 * $$\begin{align}

\bar{Y}= & \frac{1}{9}\underset{i=0}{\overset{8}{\mathop \sum }}\,\ln \left( \frac{1}-1 \right) \\ = & -1.4036 \end{align}\,\!$$

and:


 * $$\begin{align}

\bar{T}= & \frac{1}{9}\underset{i=0}{\overset{8}{\mathop \sum }}\,{{T}_{i}} \\ = & 4 \\  \underset{i=0}{\overset{8}{\mathop \sum }}\,T_{i}^{2}= & 204 \\ \underset{i=0}{\overset{8}{\mathop \sum }}\,{{T}_{i}}{{Y}_{i}}= & -68.33 \end{align}\,\!$$

Now from the least squares estimators, $$\hat{b_{i}}\,\!$$ and $$\hat{b_{0}}\,\!$$, we have:


 * $$\begin{align}

{{{\hat{b}}}_{1}}= & \frac{\underset{i=0}{\overset{8}{\mathop{\sum }}}\,{{T}_{i}}{{Y}_{i}}-N\cdot \bar{T}\cdot \bar{Y}}{\underset{i=0}{\overset{8}{\mathop{\sum }}}\,T_{i}^{2}-N\cdot {{{\bar{T}}}^{2}}} \\ = & \frac{-68.33-9\cdot 4\cdot \left( -1.4036 \right)}{204-9\cdot {{4}^{2}}} \\ = & -0.2967 \\   &  \\   {{{\hat{b}}}_{0}}= & \bar{Y}-{{{\hat{b}}}_{1}}\bar{T} \\ = & \left( -1.4036 \right)-\left( -0.2967 \right)\cdot 4.0 \\ = & -0.2168 \end{align}\,\!$$

Therefore:


 * $$\begin{align}

\widehat{b}= & {{e}^{-0.2168}} \\ = & 0.8051 \\  \widehat{k}= & -(-0.2967) \\ = & 0.2967 \end{align}\,\!$$

The Logistic reliability model that best fits the data is given by:


 * $$R=\frac{1}{1+0.8051\cdot \ \,{{e}^{-0.2967T}}}\,\!$$

The figure below shows the Reliability vs. Time plot.