*********************************************;
*** Competing failure mode data example 5 ***;
*** From Nelson handout (page 25)         ***;
*********************************************;

OPTIONS NODATE NONUMBER LS=72 PS=80;

TITLE F=SWISSB H=6 CM 'RELIABILITY ANALYSIS -- EXAMPLE 5';

DATA example5;
     INPUT strength mode $ n @@;
     censor = 0;
LABEL strength = 'BREAKING STRENGTH AT FAILURE';
CARDS;
 550 B 1   750 W 1   950 B 1   950 W 1  1150 W 3  1150 B 2
1250 B 2  1350 W 1  1450 B 2  1450 W 1  1550 B 1  1550 W 2
1850 W 1  2050 B 1
 ;
PROC RELIABILITY DATA=example5;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'DATA SET WITH BOTH FAILURE MODES';

PROC RELIABILITY DATA=example5;
   DISTRIBUTION WEIBULL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING A WEIBULL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'DATA SET WITH BOTH FAILURE MODES';

DATA wire; SET example5;
     censor = (mode='B');
PROC RELIABILITY DATA=wire;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'WITH BOND FAILURE DATA CENSORED';
PROC RELIABILITY DATA=wire;
   DISTRIBUTION WEIBULL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING A WEIBULL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'WITH BOND FAILURE DATA CENSORED';

DATA bond; SET example5;
     censor = (mode='W');
PROC RELIABILITY DATA=bond;
   DISTRIBUTION EXPONENTIAL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING AN EXPONENTIAL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'WITH WIRE FAILURE DATA CENSORED';

PROC RELIABILITY DATA=bond;
   DISTRIBUTION WEIBULL;
   PROBPLOT strength*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,25)PCT 'CONNECTION';
   NOTE F=SWISSB H=.43 CM MOVE=(58,23)PCT 'BREAKING';
   NOTE F=SWISSB H=.43 CM MOVE=(58,21)PCT 'STRENGTH DATA';
TITLE2 F=SWISSB H=.4 CM 'FITTING A WEIBULL DISTRIBUTION';
TITLE3 F=SWISSB H=.4 CM 'WITH WIRE FAILURE DATA CENSORED';
RUN;