stset exit2, failure(cod=1,2,3) id(id) scale(365.25) origin(dx) ///
>   exit(time min(dx + 365.25*10.01,mdy(12,31,2013))) 
{\smallskip}
//must ensure that those that die after follow-up time 
//   are administratively censored 
gen cod2 = cond(_d==0,0,cod) 
{\smallskip}
//center non-linear age (rcsage) at 45 years old
rcsgen age, gen(rcsage) df(3) orthog center(45)
//store knot positions in global macro
global knots `r(knots)'
//save matrix for orthogonalization 
matrix Rage = r(R) 
