The RELIABILITY Procedure

Nonparametric Confidence Intervals for Cumulative Failure Probabilities

The method used in the RELIABILITY procedure for computation of approximate pointwise and simultaneous confidence intervals for cumulative failure probabilities relies on the Kaplan-Meier estimator of the cumulative distribution function of failure time and approximate standard deviation of the Kaplan-Meier estimator. For the case of arbitrarily censored data, the Turnbull algorithm, discussed previously, provides an extension of the Kaplan-Meier estimator.

For multiply censored data, the Kaplan-Meier estimator of the cumulative distribution function at failure time $t_ i$ is $\hat{F}(t_ i) = 1 - \hat{S}(t_ i)$, where

\[ \hat{S}(t_ i) = \prod _{j=1}^{i}( 1 - \hat{p}_ j ), \]
\[ \hat{p}_ i = \frac{d_ i}{n_ i}, \]

$d_ i$ is the number of failures in the interval $(t_{i-1}, t_ i)$ , and $n_ i$ is the number of unfailed units at the beginning of the interval. This definition of the Kaplan-Meier estimator is equivalent to the one previously given.

An estimator of the variance $v_ i$ of the Kaplan-Meier estimator $\hat{F}(t_ i)$ is given by

\[ \hat{v}_ i = [\hat{S}(t_ i)]^2 \sum _{j=1}^ i\frac{\hat{p}_ j}{n_ j(1-\hat{p}_ j)} \]

An estimator of the standard deviation of $\hat{F}(t_ i)$ is $\mbox{se}_{\hat{F}} = \sqrt {\hat{v}_ i}$.

For arbitrarily censored data, the Kaplan-Meier estimator is replaced by the nonparametric maximum likelihood estimator computed with the Turnbull algorithm, and the approximate variance of the estimator of $F(t_ i)$ is computed from the inverse of the Fisher information matrix.

Pointwise Confidence Intervals

Approximate $(1-\alpha )100\% $ pointwise confidence intervals are computed as in Meeker and Escobar (1998, section 3.6) as

\[ [ F_ L, F_ U ] = \left[ \frac{\hat{F}}{\hat{F} + ( 1-\hat{F})w}, \; \; \frac{\hat{F}}{\hat{F} + ( 1-\hat{F})/w} \right] \]

where

\[ w = \exp \left[\frac{z_{1-\alpha /2}\mbox{se}_{\hat{F}}}{(\hat{F}(1-\hat{F}))} \right] \]

where $z_ p$ is the pth quantile of the standard normal distribution.

Simultaneous Confidence Intervals

Approximate $(1-\alpha )100\% $ simultaneous confidence bands valid over the lifetime interval $(t_ a, t_ b)$ are computed as the "Equal Precision" case of Nair (1984) and Meeker and Escobar (1998, section 3.8)

\[ [ F_ L, F_ U ] = \left[ \frac{\hat{F}}{\hat{F} + ( 1-\hat{F})w}, \; \; \frac{\hat{F}}{\hat{F} + ( 1-\hat{F})/w} \right] \]

where

\[ w = \exp \left[\frac{e_{a,b,1-\alpha /2}\mbox{se}_{\hat{F}}}{(\hat{F}(1-\hat{F}))}\right] \]

where the factor $x = e_{a,b,1-\alpha /2}$ is the solution of

\[ x\exp (-x^2/2)\log \left[\frac{(1-a)b}{(1-b)a}\right]/\sqrt {8\pi } = \alpha /2 \]

The time interval $(t_ a, t_ b)$ over which the bands are valid depends in a complicated way on the constants a and b defined in Nair (1984), $0 < a < b < 1$. a and b are chosen by default, so that the confidence bands are valid between the lowest and highest times corresponding to failures in the case of multiply censored data, or, to the lowest and highest intervals for which probabilities are computed for arbitrarily censored data. You can optionally specify a and b directly with the NPINTERVALS=SIMULTANEOUS(a,b) option in the PROBPLOT statement.