2016-01-13 2801 views
0

我要绘制ROC曲线以下Cox比例风险模型。如何绘制一个coxmodel(存活出现分析)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) 

有什么建议吗?你知道的任何其他包

+1

[mcve]的任何数据? – Tensibai

+1

“我试过了survivalROC软件包。”发生什么了? – doctorlove

+0

预测模型上具有时变变量需要为每个协在要求预测的持续时间的场景。你有没有建立这样一套完整的方案? –

回答

0

使用ROCR包。它有许多很好的模型诊断功能,包括AUC。

ROCR Website