options pagesize=47 linesize=64;
%include 'l11m.sas';
filename grafout 'l11.ps';
goptions device=pslmono gsfname=grafout gsfmode=append;
data emerson ; infile 'emerson.dat' ;
   input age ttt crind crtime transs transp time status fab 
      gender wbc kscore chem; control=0; 
   if _n_<2 then control=1; zero=0; run;
data fake; drop n; do n=1 to 1000;
   time=ranexp(0); cent=ranexp(0); time1=time;
   status=0; if time<cent then status=1; else time1=cent; 
   time2=time1**2; time3=time2*exp(1.0);
   output; end; run;
title1 'Simulated data with intercept 0, scale 1';
proc lifereg data=fake; model time1*status(0)=; run;
title1 'Simulated data with intercept 0, scale 2';
proc lifereg data=fake; model time2*status(0)=; run;
title1 'Simulated data with intercept 1, scale 2';
proc lifereg data=fake; model time3*status(0)=; run;
title1 'Regression of remission time on treatment, kscore';
proc lifereg data=emerson covout outest=out1;
   model crtime*crind(0)=ttt kscore/covb; 
   output out=pred p=acf control=control std=stdacf;
   run;
/* The following macro is no longer necessary */
/* %acse(out1,0 0 70 0,crtime);
%acse(out1,0 0 10 0,crtime); */
data pred; set pred; lp=log(acf); sdlp=stdacf/acf; run;
title1 'Estimates and standard errors of accel. factors';
proc print data=pred; var ttt kscore lp sdlp acf stdacf; run;
proc phreg data=emerson ;
  title1 'Regression of remission time on treatment, kscore';
   model crtime*crind(0)=ttt kscore/rl;
run;
%lifediag(emerson,crtime,crind,zero,ttt)
%lifediag(emerson,time,status,zero,ttt)
%lifediag(emerson,crtime,crind,kscore,ttt)
%lifediag(emerson,time,status,kscore,ttt)