data surv;
   set restat.t3data;
   if surv=t then exit=0;
   if t=117 then exit=0;
   
proc sort; by c t; run;
/*
proc print; var country c t xcumimf cumimf dprog enter exit surv; run;
*/
  /*  This set includes all variables  */

proc phreg data=surv;
 model surv*exit(0)=xcumimf rrat0  rir rxr0 tot0 
   f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02    / ties=exact;
 strata c;
 output out=resid resdev=dev;
 run;


 data resid;
   set resid;
   if dev=. then delete;
   ldev=lag(dev);
   lsurv=lag(surv);
   if surv=lsurv+1;


 proc print data=resid; var c surv dev; run;

 proc autoreg data=resid;
   model dev=ldev /noint;
   by c;
   run;


proc phreg data=surv;
 model surv*exit(0)=xcumimf rrat0  rir rxr0 tot0 
   f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02   / ties=exact;

 run;

 proc phreg data=surv;
 model surv*exit(0)=xcumimf rrat0  rir rxr0 tot0 
      / ties=exact;
 strata c;
 run;

 
proc phreg data=surv;
 model surv*exit(0)=xcumimf rrat0  rir rxr0 tot0 
    / ties=exact;

 run;

 /*      This is the set excluding rxr0   */

 proc phreg data=surv;
  model surv*exit(0)=xcumimf rrat0  rir tot0 
    f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02   / ties=exact;
  strata c;
  run;

  proc phreg data=surv;
  model surv*exit(0)=xcumimf rrat0  rir tot0 
    f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02   / ties=exact;
    run;

  proc phreg data=surv;
  model surv*exit(0)=xcumimf rrat0  rir tot0 
     / ties=exact;
  strata c;
  run;

   proc phreg data=surv;
  model surv*exit(0)=xcumimf rrat0  rir tot0 
     / ties=exact;
  
  run;

 
/*  This set just includes the cumulative IMF variable  */

  proc phreg data=surv;
 model surv*exit(0)=xcumimf 
    f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02   / ties=exact;
 strata c;
 run;

   proc phreg data=surv;
 model surv*exit(0)=xcumimf 
    f79 f80 f81 f82 f83 f84 f85 f86 f87 f88 f89 f90 f91 f92
   f93 f94 f95 f96 f97 f98 f99 f00 f01 f02   / ties=exact;
  run;

proc phreg data=surv;
 model surv*exit(0)=xcumimf 
    / ties=exact;
 strata c;
 run;
 

proc phreg data=surv outest=axx;
   model surv*exit(0)=xcumimf       / ties=exact;

run;