**************************************;
*** Singly censored data example 4 ***;
*** From Nelson handout (page 25)  ***;
**************************************;

OPTIONS NODATE NONUMBER LS=72 PS=80;

TITLE F=SWISSB H=6 CM 'RELIABILITY ANALYSIS -- EXAMPLE 4';
DATA example4;
     INPUT miles @@;
     k=_N_;
     IF k LE 37 THEN DO;  n=1; censor=0; END;
                ELSE DO; n=59; censor=1; END;
LABEL MILES = 'NUMBER OF MILES TO FAILURE';
CARDS;
 22.5   37.5   46.0   48.5   51.5   53.0   54.5
 57.5   66.5   68.0   69.5   76.5   77.0   78.5
 80.0   81.5   82.0   83.0   84.0   91.5   93.5
102.5  107.0  108.5  112.5  113.5  116.0  117.0
118.5  119.0  120.0  122.5  123.0  127.5  131.0
132.5  134.0  135.0
;
PROC RELIABILITY DATA=example4;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT miles*censor(1) / WAXIS=2 WFIT=2 FONT=SWISSB;
   FREQ n;
   SYMBOL1 H=1.5 V=CIRCLE W=2;
   NOTE F=SWISSB H=.43 CM MOVE=(58,26)PCT 'LOCOMOTIVE';
   NOTE F=SWISSB H=.43 CM MOVE=(58,24)PCT 'CONTROLS:';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'NUMBER OF MILES';
   NOTE F=SWISSB H=.43 CM MOVE=(58,20)PCT 'UNTIL FAILURE';
TITLE2 F=SWISSB H=.5 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';

PROC RELIABILITY DATA=example4;
   DISTRIBUTION WEIBULL;
   PROBPLOT miles*censor(1) / WAXIS=2 WFIT=2 FONT=SWISSB;
   FREQ n;
   SYMBOL1 H=1.5 V=CIRCLE W=2;
   NOTE F=SWISSB H=.43 CM MOVE=(58,26)PCT 'LOCOMOTIVE';
   NOTE F=SWISSB H=.43 CM MOVE=(58,24)PCT 'CONTROLS:';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'NUMBER OF MILES';
   NOTE F=SWISSB H=.43 CM MOVE=(58,20)PCT 'UNTIL FAILURE';
TITLE2 F=SWISSB H=.5 CM 'FITTING A WEIBULL DISTRIBUTION';

RUN;