Getting Started Example 4 for PROC RELIABILITY


/****************************************************************/
/*          S A S   S A M P L E   L I B R A R Y                 */
/*                                                              */
/*    NAME: RELGS4                                              */
/*   TITLE: Getting Started Example 4 for PROC RELIABILITY      */
/* PRODUCT: QC                                                  */
/*  SYSTEM: ALL                                                 */
/*    KEYS: reliability, lifetime data                          */
/*   PROCS: RELIABILITY                                         */
/*    DATA:                                                     */
/*                                                              */
/* SUPPORT: sasgjj                                              */
/*     REF: PROC RELIABILITY, INTRODUCTORY EXAMPLE 4.           */
/*    MISC:                                                     */
/****************************************************************/


data cracks;
   input Time units fail;
   datalines;
6.12  167 5
19.92 162 16
29.64 146 12
35.4  134 18
39.72 116 18
45.24 98  2
52.32 96  6
63.48 90 17
;

proc print data=cracks;
run;

proc reliability data=cracks;
   freq fail;
   nenter units;
   distribution Weibull;
   probplot time / readout
                   ppout
                   pconfplt
                   noconf;
run;

data turbine;
   label t1 = 'Time of Cracking (Hours x 100 )';
   input t1 t2 f;
   datalines;
.   4  0
4   .  39
.  10  4
10  .  49
.  14  2
14  .  31
.  18  7
18  .  66
.  22  5
22  .  25
.  26  9
26  .  30
.  30  9
30  .  33
.  34  6
34  .  7
.  38  22
38  .  12
.  42  21
42  .  19
.  46  21
46  .  15
;

proc reliability data = turbine;
   distribution lognormal;
   freq f;
   pplot ( t1 t2 ) / maxitem = 5000
                     ppout;
run;