In am using the Aalen Johansen estimator to plot the cumulative incidence of my specific outcome over time in the competing risk setting. The graph is shown below.
[ATTACH=CONFIG]temp_22802_1622114226758_319[/ATTACH]
For this graph, the journal requires that we include the number of individuals at risk at each of the 5-year age intervals displayed on the x-axis.
In a normal Kaplan-Meier curve, the STATA syntax could easier provide me with such information using a syntax such as:
Code:
sts graph, by(group) risktable
Does anybody have experience with this and can guide me?
A brief description of the STATA syntax that I have used to create the Aalen Johansen graph is displayed in the box below.
Code:
stset age, failure(event=1) id (id) stcompet CumInc_group1=ci if(group==1), compet1(2) stcompet CumInc_group2=ci if(group==2), compet1(2) stcompet CumInc_group3=ci if(group==3), compet1(2) gen CI_group1=CumInc_group1 if event==1 gen CI_group2=CumInc_group2 if event==1 gen CI_group3=CumInc_group3 if event==1 sort age twoway (scatter CI_group1 age, connect(J)) /// (scatter CI_group2 age, connect(J)) /// (scatter CI_group3 age, connect(J))
Any help would be much appreciated.
Thank you,
Line
0 Response to Aalen Johansen estimates of the cumulative incidence - how to include numbers at risk in figure
Post a Comment