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 phreg;
 model failtime*rfscens(0) = trt age sex;
 bayes seed = 532 nbi=3000 nmc = 20000 
	   coeffprior = uniform
       diagnostics = all plots=all;
       
run;

proc phreg;
 model failtime*rfscens(0) = trt age sex;
 bayes seed = 5432 nbi=3000 nmc = 20000 
       coeffprior = normal
       diagnostics = all plots=all;
       
run;

ods html close;
ods listing;

proc phreg;
 model failtime*rfscens(0) = trt age sex;
 bayes seed = 5432 nbi=3000 nmc = 20000
       coeffprior = normal
       piecewise = hazard(N=15 prior=gamma)
       diagnostics = all plots=all;
run;

proc phreg;
 model failtime*rfscens(0) = trt age sex;
 bayes seed = 5432 nbi=3000 nmc = 20000
       coeffprior = normal
       piecewise = hazard(N=15 prior=improper)
       diagnostics = all plots=all;
run;


       
proc phreg;
 model failtime*rfscens(0) = trt age sex;
 bayes seed = 5432 nbi=3000 nmc = 20000
       coeffprior = normal
       piecewise = loghazard(N=15 prior=normal)
       diagnostics = all plots=all;
run;