data melanoma_1684;
 infile 'H:\Bayes Workshop 2008051415\link\e1684ws.dat';
 input  case study age trt sex perform nodes breslow failtime rfscens survtime scens;
 run;
ods listing close;
ods html file="junk2.htm";
/* This statement sets the device to ActiveX */
/* to produce an ActiveX graph that will     */
/* appear in your web page. It also sets up  */
/* customizations for the graph.             */
goptions reset   = all
         device  = activex
         gunit   = pct
         border
         xpixels = 640
         ypixels = 480
         htitle  = 4
         htext   = 4
         cback   = CXFFF7CE
         colors  = (CX00FF94 CX00C4FF CX001FFF CX8400FF CXFF00D4);
ods graphics on;
proc lifereg;
 model failtime*rfscens(0) = trt age sex /dist=exponential;
 bayes seed = 5432 nbi=3000 nmc = 20000 
	   coeffprior = uniform
       diagnostics = all plots=all;
       
run;

proc lifereg;
 model failtime*rfscens(0) = trt age sex /dist=weibull;
 bayes seed = 5432 nbi=3000 nmc = 20000 
       WeibullShapePrior=gamma 
	   coeffprior = uniform
       diagnostics = all plots=all;
       
run;

ods html close;
ods listing;

proc lifereg;
 model failtime*rfscens(0) = trt age sex /dist=weibull;
 bayes seed = 5432 nbi=3000 nmc = 20000 
       WeibullShapePrior=gamma(shape=0.01,iscale=.01)
	   coeffprior = normal
       diagnostics = all plots=all;
run;

proc lifereg;
 model failtime*rfscens(0) = trt age sex /dist=Lognormal;
 bayes seed = 5432 nbi=3000 nmc = 20000
       coeffprior = uniform
       diagnostics = all plots=all;
       
run;