Data visualization with SAS programming

How do I add censoring marks to cumulative incidence function plots?

Reply
Occasional Contributor
Posts: 5

How do I add censoring marks to cumulative incidence function plots?

Hi,

 

I want to show the difference in cumulative incidence of tumor recurrence between two groups (strata "variable" in the code below), treating death as a competing event. I am having trouble displaying censoring marks on cumulative incidence function plots when using competing risks analysis. I am also having trouble displaying the number of patients at risk at each time point (but I realize that the OUTCIF function can compute a table with number at risk at each time point which can then be manually entered underneath the CIF plot). I have no trouble showing censoring marks using the cumulative incidence capabilities of the Kaplan-Meier analysis; however, Kaplan-Meier overestimates cumulative incidence because it censors death instead of treating it as a competing event. I would therefore prefer not to use Kaplan-Meier. Below is the code for CIF and for the Kaplan-Meier method, respectively. Any help to display censoring marks for CIF in SAS 9.4 would be greatly appreciated.

 

proc lifetest data=filename 

plots=cif(test);
time Months*censor(0) / eventcode=1;
strata variable;
run;

 

proc lifetest data=filename
plots=(survival(outside f ATRISK PVALUE));

strata variable;
time Months*censor(0);
run;

 

Thanks!

SAS Super FREQ
Posts: 496

Re: How do I add censoring marks to cumulative incidence function plots?

Posted in reply to sasuser110

If no one responds, you should probably ask technical support.  We have many options for customizing the survival plot.  We have an entire chapter devoted to it https://support.sas.com/documentation/onlinedoc/stat/142/kaplan.pdf  We don't provide that type of customization for any other graph.  You have unlimited ways to customize if you can make a data set with the right form.  Sometimes this involves running the procedure multiple times or running multiple procedures and then combining ODS output data set or other output data sets.  Then you can customize the template.  The default template for the CIF plot does not contain the axis table needed to display the number of subjects at risk. I don't know enough about the CIF plot to help you with this.

Ask a Question
Discussion stats
  • 1 reply
  • 155 views
  • 0 likes
  • 2 in conversation