各位大佬 问个问题呗?在运行R时报了一个这样的错,这是啥原因呢?
library(survivalROC)
setwd("F:\\TCGA_0.5\\26.model\\train")
rt=read.table("testRisk.txt", header=T, sep="\t", check.names=F, row.names=1)
rocCol=c("red", "green", "blue")
aucText=c()
pdf(file="ROC.multiTime_test.pdf",width=6,height=6)
predictTime=1
par(oma=c(0.5,1,0,1),font.lab=1.5,font.axis=1.5)
roc=survivalROC(Stime=rt$futime, status=rt$fustat, marker=rt$riskScore, predict.time=predictTime, method="KM")
plot(roc$FP, roc$TP, type="l", xlim=c(0,1), ylim=c(0,1),col=rocCol[1],
xlab="False positive rate", ylab="True positive rate",
lwd = 2, cex.main=1.3, cex.lab=1.2, cex.axis=1.2, font=1.2)
aucText=c(aucText,paste0("one year"," (AUC=",sprintf("%.3f",roc$AUC),")"))
abline(0,1)
|