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

OPTIONS NODATE NONUMBER LS=72 PS=40;

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 LIFETEST DATA= example4 PLOTS=(LS,LLS,S) OUTSURV=survive GRAPHICS;
   TITLE F=SWISSB H=.6 CM 'RELIABILITY ANALYSIS: EXAMPLE 4';
   TIME miles*censor(1);
   FREQ n;
   SYMBOL1 H=1 V=DOT W=2;
   NOTE F=SWISSB H=.45 CM MOVE=(18,67)PCT 'LOCOMOTIVE CONTROLS:';
   NOTE F=SWISSB H=.45 CM MOVE=(18,65)PCT 'NUMBER OF MILES';
   NOTE F=SWISSB H=.45 CM MOVE=(18,63)PCT 'UNTIL FAILURE';
PROC PRINT DATA=survive;
RUN;