Hi Stata forum,
I try to post once more regarding my problem - now hopefully its easier to read my codes, as I corrected my post.

I would like to estimate the 5-year cumulative incidence proportions with a competing risk (CIF)(and 95%CI) among my exposed(exp==1) and non-exposed(exp==0) according to a list of 22 covariates. My event of interest is event==1, whereas my competing risk is event==2.
So far I have used the following syntax to make up a cumulative incidence graph(preceeded by setting up my data for survival analysis):

Code:
Stcompet cif=ci hilim=hi lowli=lo, compet(2) by(exp)
gen cif_exp1=cif if event==1 & exp==1
gen cif_exp0=cif if event==1 & exp==0
graph twoway line cif_exp* _t, connect(JJ)
Lets say I would like to calculate the 5-year CIF for males only(sex==1) among the exposed. I have done this the following way
Code:
gen c2=cif if event==1 & exp==1 & sex==1
sort _t
replace c2=c2[_n-1] if mi(c2)
list _t c2 if _t<=5 & _t[_n+1]>5
The same is done for the lower and higher 95%confidence limits.
First, I am unsure if I am allowed to stratify on sex, since my variable cif is calculated for the entire cohort by exposure group.


I have also used the stcomlist command to do the same:

stcomlist if sex==1, compet1(2) at (1 5) by(eks)

Which procedure is the better or does there exist another less time-consuming way to do it?

All help will be much appreciated.

Best regards Kirsten