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

OPTIONS NODATE NONUMBER LS=72 PS=80;

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
 ;
TITLE F=SWISSB H=6 CM 'RELIABILITY ANALYSIS -- EXAMPLE 3';

PROC RELIABILITY DATA=example3;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT hours*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,24)PCT 'FAN FAILURES';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'LIFETIME DATA';
TITLE2 F=SWISSB H=.5 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';

PROC RELIABILITY DATA=example3;
   DISTRIBUTION WEIBULL;
   PROBPLOT hours*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,24)PCT 'FAN FAILURES';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'LIFETIME DATA';
TITLE2 F=SWISSB H=.5 CM 'FITTING A WEIBULL DISTRIBUTION';

RUN;