Code for analysis

*** code for survival analysis
 
stset exitdate16, fail(tbCP16) origin(skint_date) enter(skint_date) id(ident) scale(365.25)
** and similarly for the other outcomes
 
stsplit timeband4, after(skint_date) at (2,10,20)
strate timeband4 large17 , per(1000)
 
xi: streg i.large17 if timeband4==0, dist(exp)nolog
xi: streg i.large17 i.age3 if timeband4==0, dist(exp)nolog
xi: streg i.large17 i.age3 i.sex i.area if timeband4==0, dist(exp)nolog
 
** and similarly for the other timebands