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):
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
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
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)
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
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
