/* Programs for Lab note 3 */
/* Prepared by Sheng Luo */

DATA recid;
	SET "C:\courses\Survival\2006\labs\lab3\recid";
RUN;

proc print;
run;

proc phreg data=recid;
   model week*arrest(0)=fin age race wexp mar paro prio;
run;

proc phreg data=recid;
   model week*arrest(0)=fin age race wexp mar paro prio employed
         / ties=efron;
   array emp{*} emp1-emp52;
   employed=emp{week};
run;

proc phreg data=recid;
   where week>1;
   model week*arrest(0)=fin age race wexp mar paro prio employed
         / ties=efron;
   array emp{*} emp1-emp52;
   employed=emp{week-1};
run;


proc phreg data=recid;
   model week*arrest(0)=fin age race wexp mar paro prio fintime
         / ties=efron;
   fintime=fin*week;
run;

proc phreg data=recid;
   where week>1;
   model week*arrest(0)=fin age wexp mar paro prio employed
         / ties=efron;
   array emp{*} emp1-emp52;
   employed=emp{week-1};
   strata race;
run;


proc phreg data=recid;
   where week>1;
   model week*arrest(0)=fin age race wexp mar paro prio / ties=efron;
   baseline out=a survival=s logsurv=ls loglogs=lls;
run;

proc print data=a;
run;

proc gplot data=a;
   symbol1 value=none interpol=join;
   plot lls*week;
run;