I have estimated the overall crude incidence rate of diabetes in people with hypertension (the study population) between 2010-2019.

I used the following command to calculate the overall incidence rate: -stset enddate, failure(diabetes==1) enter(startdate) id(id) scale(365.25)-
enddate = earliest of diabetes diagnosis, death, exit from database or 31st Dec 2019.
startdate = date of first hypertension diagnosis

I now would like to stratify the incidence rate by calendar year and produce a graph with calendar year on the x-axis and incidence rate per 1000 on the y-axis.

I attempted to use the command: -stptime, by (startyear) per(100) dd(2)- where (startyear) is the year hypertension was diagnosed for each id. However, I do not think this is correct.

I have done some reading and come across -stsplit- and -strate-. How would I go about stratifying incidence rate and graph it by calendar year? I've put a dummy example of my data below.

Thank you.

Code:
* Example generated by -dataex-. For more info, type help dataex
clear
input str2 id byte diabetes int(startdate enddate) float survival_time
"1"  0 19304 21816  6.877481
"2"  0 20021 20973  2.606434
"3"  0 18497 21796  9.032169
"4"  0 18883 21762  7.882272
"5"  1 20273 21209  2.562628
"6"  0 18868 21803  8.035592
"7"  1 18512 20755  6.140999
"8"  0 19908 20408 1.3689255
"9"  0 19547 20828  3.507187
"10" 0 20409 20983 1.5715264
"11" 1 19643 21817  5.952087
"12" 0 21019 21638 1.6947297
"13" 0 20614 21091 1.3059548
"14" 1 19207 21817  7.145791
"15" 0 20619 21817  3.279945
"16" 1 18517 20078  4.273785
"17" 0 19036 19591  1.519507
"18" 0 19598 20752   3.15948
"19" 1 20988 21628 1.7522244
"20" 0 18713 21816  8.495551
end
format %td startdate
format %td enddate