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)
有什麼建議嗎?你知道的任何其他包
[mcve]的任何數據? – Tensibai
「我試過了survivalROC軟件包。」發生什麼了? – doctorlove
預測模型上具有時變變量需要爲每個協在要求預測的持續時間的場景。你有沒有建立這樣一套完整的方案? –