****************************************;
*** Multiply censored data example 1 ***;
*** From Nelson handout (page 13)    ***;
****************************************;

OPTIONS NODATE NONUMBER LS=72 PS=80;

DATA example1;
     INPUT cycles censor n @@;
LABEL CYCLES = 'NUMBER OF CYCLES TO FAILURE';
TITLE F=SWISSB H=6 CM 'RELIABILITY ANALYSIS -- EXAMPLE 1';

CARDS;
  45 1 1     47 0 1     73 0 1    136 1 5    145 0 1
 190 1 2    281 1 1    311 0 1    417 1 1    485 1 2
 490 0 1    569 1 1    571 1 1    571 0 1    575 0 1
 608 1 12   608 0 2    630 0 1    670 0 2    731 1 1
 838 0 1    964 0 2   1164 1 7   1198 1 1   1198 0 1
1300 1 3

PROC RELIABILITY DATA=example1;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT cycles*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 'LIFETIMES OF';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'APPLIANCE';
   NOTE F=SWISSB H=.43 CM MOVE=(58,20)PCT 'COMPONENTS';
TITLE2 F=SWISSB H=.5 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';

PROC RELIABILITY DATA=example1;
   DISTRIBUTION WEIBULL;
   PROBPLOT cycles*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 'LIFETIMES OF';
   NOTE F=SWISSB H=.43 CM MOVE=(58,22)PCT 'APPLIANCE';
   NOTE F=SWISSB H=.43 CM MOVE=(58,20)PCT 'COMPONENTS';
TITLE2 F=SWISSB H=.5 CM 'FITTING A WEIBULL DISTRIBUTION';

RUN;