Chapter Contents |
Previous |
Next |
The RELIABILITY Procedure |
Inspection (Months) | Number | ||
Start | End | Cracked | Cumulative |
0 | 6.12 | 5 | 5 |
6.12 | 19.92 | 16 | 21 |
19.92 | 29.64 | 12 | 33 |
29.64 | 35.40 | 18 | 51 |
35.40 | 39.72 | 18 | 69 |
39.72 | 45.24 | 2 | 71 |
45.24 | 52.32 | 6 | 77 |
52.32 | 63.48 | 17 | 94 |
63.48 | Survived | 73 | 167 |
Table 30.1 shows the time in months of each inspection period and the number of cracked parts found in each period. These data are said to be interval censored since only the time interval in which failures occurred is known, not the exact failure times. Seventy-three parts had not cracked at the last inspection, which took place at 63.48 months. These 73 lifetimes are right censored, since the lifetimes are known only to be greater than 63.48 months.
The interval data in this example is read from a SAS data set with a special structure. All units must have a common inspection schedule. This type of interval data is called readout data. The SAS data set named CRACKS, shown in Figure 30.8, provides the data in Table 30.1 with this structure. The variable TIME is the inspection time, that is, the upper endpoint of each interval. The variable UNITS is the number of unfailed units at the beginning of each interval, and the variable FAIL is the number of units with cracks at the inspection time.
|
The following statements use the RELIABILITY procedure to produce the probability plot in Figure 30.9 for the data in the data set CRACKS.
symbol v=plus c=blue; proc reliability data=cracks; freq fail; nenter units; distribution weibull; probplot time / readout pconfplt noconf cframe = ligr; inset / ctext = black cfill = ywh ; run;The FREQ statement specifies that the variable FAIL provides the number of failures in each interval. The NENTER statement specifies that the variable UNITS provides the number of unfailed units at the beginning of each interval. The DISTRIBUTION statement specifies that the Weibull distribution is used for parameter estimation and probability plotting. The PROBPLOT statement requests a probability plot of the data.
The PROBPLOT statement option READOUT indicates that the data in the CRACKS data set are readout (or interval) data. The option PCONFPLT specifies that confidence intervals for the cumulative probability of failure are to be plotted. The confidence intervals for the cumulative probability are based on the binomial distribution for time intervals until right censoring occurs. For time intervals after right censoring occurs, the binomial distribution is not valid, and a normal approximation is used to compute confidence intervals.
The option NOCONF suppresses the display of confidence intervals for distribution percentiles in the probability plot.
A listing of the tabular output produced by the preceding SAS statements is shown in Figure 30.10. By default, the specified Weibull distribution is fitted by maximum likelihood. The line plotted on the probability plot and the tabular output summarize this fit. For interval data, the estimated cumulative probabilities and associated confidence intervals are tabulated. In addition, general fit information, parameter estimates, percentile estimates, standard errors, and confidence intervals are tabulated.
In this example, the number of unfailed units at the beginning of an interval minus the number failing in the interval is equal to the number of unfailed units entering the next interval. This is not always the case since some unfailed units might be removed from the test at the end of an interval; that is, they might be right censored. The special structure of the input SAS data set required for interval data enables the RELIABILITY procedure to analyze this more general case.
Chapter Contents |
Previous |
Next |
Top |
Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.