topic Re: Concordance index for Competing risks (PROC PHREG) in Statistical Procedures

Concordance index for Competing risks (PROC PHREG)

ks17 — Wed, 21 Feb 2018 19:51:07 GMT

I am fitting a Fine and Gray regression model to competing risks data. I am trying to compute the time-dependent AUC for this model, but the "plots=auc rocoptions(method=ipcw(cl seed=1234))" option does not work. I get the note "Neither concordance nor ROC analysis is performed when the EVENTCODE option in the MODEL statement is specified". However, I need that option to specify fitting a Fine and Gray regression. The data is coded as having Status 0 if censored, 1 if experience the event of interest, 2 if they experience the competing event. Here is my SAS code:

proc phreg data=train plots=auc rocoptions(method=ipcw(cl seed=1234));
class X1 X2;
model Time*Status(0)=X1 X2/eventcode=1;
run;

Any suggestions for how I can obtain the concordance index from a Fine and Gray regression model? Thanks in advance!

Re: Concordance index for Competing risks (PROC PHREG)

k-san — Tue, 08 May 2018 20:18:56 GMT

Did you ever find a solution to this? We are having the same issue.

Re: Concordance index for Competing risks (PROC PHREG)

ks17 — Tue, 08 May 2018 20:45:37 GMT

I emailed one of the authors of "Analyzing Survival Data with Competing Risks using SAS Software" in February 2018. They were quick to respond and said that the AUC and time dependent ROC method (using the IPCW) could be carried over to the Fine-Gray model, but that it was not yet available in PHREG and that it would be a future project.

For now, I am using the FGR function in the riskRegression package in R (https://cran.r-project.org/web/packages/riskRegression/riskRegression.pdf) for the Fine-Gray model, and the timeROC package in R to get the time-dependent AUC.

Re: Concordance index for Competing risks (PROC PHREG)

Bowen2006 — Fri, 21 Dec 2018 20:51:19 GMT

I used cause specific model, and created AUC at different time point, it works.

Try this:

proc phreg data=train plots=auc rocoptions(method=ipcw(cl seed=1234));
class X1 X2;
model Time*Status(0, 2)=X1 X2;
run;

Re: Concordance index for Competing risks (PROC PHREG)

FP_Italy — Wed, 22 Sep 2021 10:34:15 GMT

It is not a Fine and Gray model, so it works

Re: Concordance index for Competing risks (PROC PHREG)

Season — Sun, 08 Dec 2024 06:59:20 GMT

There is a paper on the concordance statistic for competing risks models: Concordance for prognostic models with competing risks | Biostatistics | Oxford Academic. The authors themselves provided R code in the supplemental material of the paper: Supplementary material to “Concordance for prognostic models with competing risks”. Another person has also developed an R package for this calculation based upon the paper mentioned above: R: Concordance index in the Presence of Competing Risks. An issue to note is that independence of observations was assumed in the paper. Therefore, this method is not applicable in the analysis of clustered data with competing risks.