Programming the statistical procedures from SAS

Competing risk graph question

Reply
New Contributor
Posts: 2

Competing risk graph question

Hi 

 

I am trying to perform competing risk analysis for outcomes using a patient cohort data on SAS studio release 9.04. I checked codes on line and found that we could use proc phreg or %macro cuminc. I am planning on using proc phreg using these codes on sas support site

ods graphics on;
proc phreg data=Bmt plots(overlay=stratum)=cif;
   class Disease (order=internal ref=first);
   model T*Status(0)=Disease / eventcode=1;
   Hazardratio 'Pairwise' Disease / diff=pairwise;
   baseline covariates=Risk out=out1 cif=_all_ / seed=191;
run;

 My question is how do I code so that the graph shows all outcome status. This code only plots for the outcome of interest (eventcode =1). Thanks

Super Contributor
Posts: 287

Re: Competing risk graph question

[ Edited ]

You can not calculate all the survival curves in only one run of phreg. You can of course run phreg for each outcome. The Fine and Gray method which proc phreg use only take account for one outcome at a time.

I am not sure that if you run phreg for each outcome that the total sum will add up to one (exactly). 
If you like, I can give you a macro which calculate all curves in one time using the Aalen Johansen estimator.

New Contributor
Posts: 2

Re: Competing risk graph question

What you say is accurate - I tried to do individual phregs and it did not add up to one. Others suggested individual KM for each outcome, which will work but is obviously univariate. I would appreciate if you can send me your macro. Thanks.
Super Contributor
Posts: 287

Re: Competing risk graph question

Here you have the program that calculate Aalen Johansen estimates. It was not written as a macro, so bit editing will perhaps be neccessary. I hope you find it useful.

Attachment
Ask a Question
Discussion stats
  • 3 replies
  • 245 views
  • 1 like
  • 2 in conversation