****************************************;
*** Multiply censored data example 3 ***;
*** From Nelson handout (page 21)    ***;
****************************************;

OPTIONS NODATE NONUMBER LS=72 PS=40;

DATA example3;
     INPUT hours censor n @@;
LABEL hours = 'NUMBER OF HOURS TO FAILURE';
CARDS;
  450 0 1    460 1 1   1150 0 2   1560 1 1   1600 0 1   1660 1 1
 1850 1 5   2030 1 3   2070 0 2   2080 0 1   2200 1 1   3000 1 4
 3100 0 1   3200 1 1   3450 0 1   3750 1 2   4150 1 4   4300 1 4
 4600 0 1   4850 1 4   5000 1 3   6100 1 3   6100 0 1   6300 1 1
 6450 1 2   6700 1 1   7450 1 1   7800 1 2   8100 1 2   8200 1 1
 8500 1 3   8750 1 2   8750 0 1   9400 1 1   9900 1 1  10100 1 3
11500 1 1
 ;
PROC LIFETEST DATA= example3 PLOTS=(LS,LLS,S) OUTSURV=survive GRAPHICS;
   TITLE F=SWISSB H=.6 CM 'RELIABILITY ANALYSIS: EXAMPLE 3';
   TIME hours*censor(1);
   FREQ n;
   SYMBOL1 H=1 V=DOT W=2;
   NOTE F=SWISSB H=.45 CM MOVE=(18,67)PCT 'FAN FAILURES';
   NOTE F=SWISSB H=.45 CM MOVE=(18,65)PCT 'LIFETIME DATA';
PROC PRINT DATA=survive;
RUN;