*********************************************; *** Multiply censored data -- Homework #7 ***; *********************************************; DM 'LOG;CLEAR;OUT;CLEAR;'; options nodate nonumber ls=76 ps=54; data hw7; input Y censor @@; logY = log(Y); lines; 1.151 0 1.171 0 1.248 0 1.331 0 1.381 0 1.499 1 1.508 0 1.534 0 1.577 0 1.584 0 1.667 1 1.695 1 1.710 1 1.955 0 1.965 1 2.012 0 2.051 0 2.076 0 2.109 1 2.116 0 2.119 0 2.135 1 2.197 1 2.199 0 2.227 1 2.250 0 2.254 1 2.261 0 2.349 0 2.369 1 2.547 1 2.548 1 2.738 0 2.794 1 2.883 0 2.883 0 2.910 1 3.015 1 3.017 1 3.793 0 ; /* PROC LIFETEST DATA= hw7 PLOTS=(LS,LLS,S) OUTSURV=survive; TITLE 'RELIABILITY STUDY OF Switch Lifetimes'; TIME Y*censor(0); SYMBOL1 H=1 V=DOT W=1.5; NOTE F=SWISS H=.45 CM MOVE=(18,63)PCT 'Switch Lifetimes'; */ *****************************************************************; /* PROC RELIABILITY DATA=hw7; DISTRIBUTION EXPONENTIAL; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING AN EXPONENTIAL DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION EXTREME; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING AN EXTREME VALUE DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION LOGISTIC; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING A LOGISTIC DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION LOGLOGISTIC; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING A LOGLOGISTIC DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION LOGNORMAL; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING A LOGNORMAL DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION NORMAL; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING A NORMAL DISTRIBUTION'; */ /* PROC RELIABILITY DATA=hw7; DISTRIBUTION WEIBULL; PROBPLOT Y*censor(0) / WAXIS=1.5 WFIT=1.5 FONT=SWISSB; SYMBOL1 H=1.5 V=CIRCLE W=2; TITLE F=SWISSB H=.5 CM 'FITTING A WEIBULL DISTRIBUTION'; */ RUN;