I want to compare cumulative incidence for two arms of a cohort at time _t and generate confidence interval for the difference.

I have used the following commands to obtain the cumulative incidence for each cohort:

stset timeelapsed, failure(failure) id(hcno)
sts list, failure by(Cohort)

How can I calculate the difference in cumulative incidence (that is the failure function) for Cohort==1 and Cohort==0 and generate a confidence interval for that difference using 500 resamples? Any advice is greatly appreciated