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