****************************************;
*** Multiply censored data example 2 ***;
****************************************;

OPTIONS NODATE NONUMBER LS=72 PS=80;

DATA example2;
     INPUT failtime censor @@;
LABEL failtime = 'TIME TO FAILURE IN HOURS';
TITLE F=SWISSB H=6 CM 'RELIABILITY ANALYSIS -- EXAMPLE 2';
CARDS;
 42.1 1   77.8 0   83.3 1   88.7 0  101.1 0
105.9 0  117.0 0  126.9 0  138.7 0  148.9 0
151.3 1  157.3 0  163.8 0  177.2 1  194.3 1
195.6 1  207.0 0  215.3 1  217.4 0  258.8 1
;
PROC RELIABILITY DATA=example2;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT failtime*censor(1) / WAXIS=2 WFIT=2 FONT=SWISSB;
   SYMBOL1 H=1.5 V=CIRCLE W=2;
   NOTE F=SWISSB H=.4 CM MOVE=(58,22)PCT 'GRINDER LIFETIMES';
TITLE2 F=SWISSB H=.5 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';

PROC RELIABILITY DATA=example2;
   DISTRIBUTION WEIBULL;
   PROBPLOT failtime*censor(1) / WAXIS=2 WFIT=2 FONT=SWISSB;
   SYMBOL1 H=1.5 V=CIRCLE W=2;
   NOTE F=SWISSB H=.4 CM MOVE=(58,22)PCT 'GRINDER LIFETIMES';
TITLE2 F=SWISSB H=.5 CM 'FITTING A WEIBULL DISTRIBUTION';
RUN;