Re: [R] Cox model+ROCR

From: Terry Therneau <>
Date: Tue, 04 Mar 2008 08:10:05 -0600 (CST)

> I am trying to build a cox model and then perform ROC analysis in order to
> retrieve some genes that are correlated with breast cancer. When I calculate

 Extension of ROC values to the censored data case is handled by the rcorr.cens  function found in the Hmisc library. See the references and examples there.

    In short, one interpretation of the ROC statistic is in terms of concordance. Consider all pairs of "comparable" subjects. For a yes/no outcome, this is all pairs where one of the subjects is a yes and the other a no. The c-statistic = % of all pairs where the model's score correctly ordered them (higher score for the yes). For a yes/no outcome this = the ROC. For your endpoint it is just somewhat harder to count up "comparables": they are all pairs for which I can determine which relapsed first. Relapse at 20 vs censor at 15, for instance, is not comparable.     

    Terry Therneau mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Tue 04 Mar 2008 - 14:12:19 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 04 Mar 2008 - 14:30:20 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive