设为首页收藏本站获取积分奖励

R语言中文网

 找回密码
 立即注册
查看: 2991|回复: 0

[实例投稿提交] 各位大佬 问个问题呗?在运行R时报了一个这样的错,这是啥原因呢

[复制链接]
发表于 2021-4-4 10:08:53 | 显示全部楼层 |阅读模式
各位大佬 问个问题呗?在运行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)

微信图片_20210404100717.png
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|Archiver|R语言中文网    

GMT+8, 2024-5-1 10:03 , Processed in 0.029103 second(s), 25 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表