Resources

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;