|
|
|
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: RELGS8 */
/* TITLE: Getting Started Example 8 for PROC RELIABILITY */
/* PRODUCT: QC */
/* SYSTEM: ALL */
/* KEYS: reliability, lifetime data */
/* PROCS: RELIABILITY */
/* DATA: */
/* */
/* SUPPORT: sasgjj */
/* REF: PROC RELIABILITY, INTRODUCTORY EXAMPLE 8. */
/* MISC: */
/****************************************************************/
data defrost;
t1=t2;
t2+1;
label t1 = 'Months';
input r n;
datalines;
83 22914
35 22914
23 22914
15 22914
22 22914
16 22914
13 22911
12 22875
15 22851
15 22822
24 22785
12 22745
7 2704
11 2690
15 2673
6 2660
8 2632
9 2610
9 2583
5 2519
6 2425
6 2306
6 2188
5 2050
7 862
5 845
5 817
6 718
3 590
;
run;
data d2 ( keep = n rename = ( n = n1 ) );
set defrost;
if _n_ > 1;
run;
data d3;
n1 = 0;
run;
data d2 ;
set d2 d3;
run;
data defrost( keep = t1 t2 r c
rename = ( t1 = Lower t2 = Upper r = Recurrences c = Censored ) );
set defrost;
set d2;
c = n - n1;
run;
proc print data=Defrost;
run;
proc reliability data=defrost;
mcfplot ( interval = Lower Upper
recurrences = Recurrences
censor = Censored ) / plotsymbol = X
vaxis = 0 to .12 by .04
interpolate = join;
run;