2016-01-13 59 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