0

我想为以下 Cox 比例风险模型绘制 ROC 曲线。

Cox 比例风险模型 - 使用时变协变量

coxph_tvc <- coxph(Surv(time0, time1, event_tvc) ~ channelname)

我尝试了survivalROC 包。

    library(survivalROC)    
    predsurv<- predict(coxph_tvc, type = "lp")    
    nobs <- NROW(data_tvc)
    cutoff <- 3650      
    rocfit <- survivalROC.C( Stime = time1,
                                status = event_tvc,
                                marker = channelname,
                                predict.time = cutoff,
                                span = 1500)
   plot(rocfit$FP, rocfit$TP, type = "l",
         xlim = c(0,1), ylim = c(0,1),
         xlab = paste( "FP \n AUC =",round(rocfit$AUC,3)),
         ylab = "TP",main = "Year = 1" )
    abline(0,1)

有什么建议吗？你知道其他的包吗

4

1 回答 1

0

使用 ROCR 包。它有很多很好的模型诊断绘图功能，包括 AUC。

ROCR 网站

于 2016-06-20T05:18:35.353 回答