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

OPTIONS NODATE NONUMBER LS=72 PS=80;

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 LIFETEST DATA= example5 PLOTS=(LS,LLS,S) OUTSURV=survive GRAPHICS;
   TITLE  F=SWISSB H=.6 CM 'RELIABILITY ANALYSIS: EXAMPLE 5';
   TITLE2 F=SWISSB H=.4 CM 'DATA SET WITH BOTH FAILURE MODES';
   TIME strength*censor(1);
   FREQ n;
   SYMBOL1 H=1 V=DOT W=2;
   NOTE F=SWISSB H=.45 CM MOVE=(18,65)PCT 'CONNECTION';
   NOTE F=SWISSB H=.45 CM MOVE=(18,63)PCT 'BREAKING';
   NOTE F=SWISSB H=.45 CM MOVE=(18,61)PCT 'STRENGTH DATA';
PROC PRINT DATA=survive;

DATA wire; SET example5;
     censor = (mode='B');
PROC LIFETEST DATA= wire PLOTS=(LS,LLS,S) OUTSURV=survive GRAPHICS;
   TITLE  F=SWISSB H=.6 CM 'RELIABILITY ANALYSIS: EXAMPLE 5';
   TITLE2 F=SWISSB H=.4 CM 'WITH BOND FAILURE DATA CENSORED';
   TIME strength*censor(1);
   FREQ n;
   SYMBOL1 H=1 V=DOT W=2;
   NOTE F=SWISSB H=.45 CM MOVE=(18,65)PCT 'CONNECTION';
   NOTE F=SWISSB H=.45 CM MOVE=(18,63)PCT 'BREAKING';
   NOTE F=SWISSB H=.45 CM MOVE=(18,61)PCT 'STRENGTH DATA';
PROC PRINT DATA=survive;

DATA bond; SET example5;
     censor = (mode='W');
PROC LIFETEST DATA= bond PLOTS=(LS,LLS,S) OUTSURV=survive GRAPHICS;
   TITLE  F=SWISSB H=.6 CM 'RELIABILITY ANALYSIS: EXAMPLE 5';
   TITLE2 F=SWISSB H=.4 CM 'WITH WIRE FAILURE DATA CENSORED';
   TIME strength*censor(1);
   FREQ n;
   SYMBOL1 H=1 V=DOT W=2;
   NOTE F=SWISSB H=.45 CM MOVE=(18,65)PCT 'CONNECTION';
   NOTE F=SWISSB H=.45 CM MOVE=(18,63)PCT 'BREAKING';
   NOTE F=SWISSB H=.45 CM MOVE=(18,61)PCT 'STRENGTH DATA';
PROC PRINT DATA=survive;
RUN;