timeROC出图时显示不全

 右侧被占用

原因:Rstudio右下方的plots窗口太小

解决方法:向上向左增大窗口

 Time_ROC代码

####timeROC####
setwd("C:\\Users\\86189\\Desktop\\Time_ROC")
#R包
#install.packages("timeROC")
#install.packages("survival")
library(timeROC)
library(survival)
library(tidyverse)

#2.2 数据的整理与载入
exp_sur <- read.table("exp_sur.txt", header=T,sep="\t", check.names=F, row.names=1)
exp_sur$OS.time <- exp_sur$OS.time/365 #按年计算
exp_sur_01A <- exp_sur[substr(rownames(exp_sur),14,16) == "01A",]#将行名第14到16位为“01A”的取出来
write.table(exp_sur_01A, file = "exp_sur_01A.txt",sep = "\t",row.names = T,col.names = NA,quote = F)
#2.3 构建ROC曲线函数
ROC3 <- timeROC(T=exp_sur_01A$OS.time,   #结局时间
                delta=exp_sur_01A$OS,   #结局指标,delta:事件状态. 删失数据编码为0
                marker=exp_sur_01A$SPP1,   #预测变量,marker :计算ROC的biomaker,默认是marker值越大,事件越可能发生;反之的话,前面加-号
                cause=1,   #阳性结局指标数值,cause:事件结局。没有竞争风险(Without competing risks)中,必须是非删失数据的编码方式,一般为1
                weighting="marginal",   #weighting:计算方法,默认是"marginal",KM模型;weighting="cox" 和weighting="aalen"分别为COX模型和additive Aalen 模型
                times=c(1, 3, 5),   #times:想计算的ROC曲线的时间节点,选取1年,3年和5年的生存率
                iid=TRUE)#: 默认值iid = FALSE。iid = TRUE 才会保存置信区间,但是样本量大了后,耗时耗资源
ROC3   #查看模型变量信息
ROC3$AUC #查看AUC值
confint(ROC3)$CI_AUC #使用confint计函数计算95%的置信区间

#2.4 绘制ROC曲线
plot(ROC3,#time是时间点,col是线条颜色
     time=1, col="red", lwd=2,title = "") #title = "",图片的title为空 
plot(ROC3,
     time=3, col="green", add=TRUE, lwd=2)   
plot(ROC3,#add指是否添加在上一张图中
     time=5, col="blue", add=TRUE, lwd=2)
legend("bottomright",c("Y-1 AUC=0.69","Y-3 AUC=0.63","Y-5 AUC=0.65"),
       col=c("red","green","blue"),
       lty=1,lwd=2)   #添加标签信息.
#lwd:line width,线条宽度,默认为1.
#lty:;line type,线条类型,1为实线
dev.off()
 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值