 foreach age in 50 65 80 {\lbr}
  2.   rcsgen, scalar(`age') knots($knots) rmatrix(Rage) gen(c) center(45)
{\smallskip}
  3.   global cancer_mostdep_`age'_male sex_cancer 0 sex_cvd 0 sex_other 0 ///
>         mostdep_cancer 1  mostdep_cvd 0  mostdep_other 0   ///
>         rcsage1_cancer `=c1' rcsage2_cancer `=c2' rcsage3_cancer `=c3' ///
>         rcsage1_other 0 rcsage2_other 0 rcsage3_other 0 ///
>         rcsage1_cvd 0 rcsage2_cvd 0 rcsage3_cvd 0 _cancer 1 _cvd 0 _other 0
{\smallskip}
  4.   global cancer_leastdep_`age'_male sex_cancer 0 sex_cvd 0 sex_other 0 ///
>         mostdep_cancer 0 mostdep_cvd 0 mostdep_other 0 ///
>         rcsage1_cancer `=c1' rcsage2_cancer `=c2' rcsage3_cancer `=c3' ///
>         rcsage1_other 0 rcsage2_other 0 rcsage3_other 0 ///
>         rcsage1_cvd 0 rcsage2_cvd 0 rcsage3_cvd 0 _cancer 1 _cvd 0 _other 0
  5. {\rbr}