Probability Plot for Engine Fan Data
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: RELPPLT1 */
/* TITLE: Probability Plot for Engine Fan Data */
/* PRODUCT: QC */
/* SYSTEM: ALL */
/* KEYS: */
/* PROCS: RELIABILITY */
/* DATA: Nelson, W. Applied Life Data Analysis, p. 336 */
/* REF: SAS/QC Software: The RELIABILITY Procedure */
/* */
/* */
/****************************************************************/
options nodate nostimer nonumber source2 ls=76 ps=80;
data fan;
input lifetime censor @@;
lifetime = lifetime/1000;
cards;
2070 0 2080 0 2200 1 3000 1
3000 1 3000 1 3000 1 3100 0
3200 1 3450 0 3750 1 3750 1
4150 1 4150 1 4150 1 4150 1
4300 1 4300 1 4300 1 4300 1
450 0 460 1 1150 0 1150 0
1560 1 1600 0 1660 1 1850 1
1850 1 1850 1 1850 1 1850 1
2030 1 2030 1 2030 1 2070 0
7450 1 7800 1 7800 1 8100 1
8100 1 8200 1 8500 1 8500 1
8500 1 8750 1 8750 0 8750 1
9400 1 9900 1 10100 1 10100 1
10100 1 11500 1 4600 0 4850 1
4850 1 4850 1 4850 1 5000 1
5000 1 5000 1 6100 1 6100 0
6100 1 6100 1 6300 1 6450 1
6450 1 6700 1
;
title 'Probability Plot for Fan Data';
symbol c=blue ;
proc reliability;
distribution weibull;
pplot lifetime*censor( 1 ) / cfit = yellow
cframe = megr
ccensor = ligr
vref (intersect)= 15
vreflabel = '15th Percentile'
cvref = orange;
inset n nfail confidence / cfill = ligr
ctext = yellow
cframe = green ;
run;