library(survival)
d=sampleData(100,outcome="survival")
nd=sampleData(100,outcome="survival")
f1=coxph(Surv(time,event)~X1+X6+X8,data=d,x=TRUE,y=TRUE)
f2=coxph(Surv(time,event)~X2+X5+X9,data=d,x=TRUE,y=TRUE)
xx=Score(list(f1,f2), formula=Surv(time,event)~1,
data=nd, metrics="auc", nullModel=FALSE, times=seq(3:10))
aucgraph <- qplotAUC(xx)
qplotAUC(xx,confint=TRUE)+ggtitle("AUC")+theme_classic()
qplotAUC(xx,type="contrasts")
a=qplotAUC(xx,type="contrasts",confint=TRUE)
a+theme_bw()
Run the code above in your browser using DataLab