右侧被占用

原因：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()