I used the following SAS code to generate cumulative incidence functions (CIF). However, I'd like to know how to test the heterogeneity on the cumulative hazards of these multiple trts for recurrent events. Now my plan is to use some specific tests (i.e.,Mann-Kendall test ) to compare the generated CIFs of these trts. Or I can directly use the p-value of the HR of the phreg model to compare these groups? Any advice? Thanks.
title2 'PWP Total Time Model with Common Effects'; proc phreg data=Bladder2; model (tstart,tstop) * status(0) = Trt Number Size; strata Visit; run;
Original code can be access from
https://documentation.sas.com/doc/en/statug/15.2/statug_phreg_examples10.htm
1) You can use PROC LIFETEST to do Gray Test.
proc lifetest data=sashelp.bmt plots=cif(test ) notable ;
time t*status(0)/eventcode=1;
strata group/order=internal;
run;
2) Or did you try to do "Assessing the proportional hazards assumption"? Check this at bottom of paper.
https://stats.oarc.ucla.edu/sas/seminars/sas-survival/
Catch the best of SAS Innovate 2025 — anytime, anywhere. Stream powerful keynotes, real-world demos, and game-changing insights from the world’s leading data and AI minds.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
Ready to level-up your skills? Choose your own adventure.