pROC包绘制诊断试验ROC曲线

gja <- read.csv("C:\\Users\\Administrator\\Desktop\\gja2.csv",header = TRUE)
View(gja)
str(gja)
install.packages("pROC")
install.packages("ggplot2")

library(pROC)
library(ggplot2)
#建立曲线

rocobj1<-roc(gja$dia,gja$dia_tbs)
rocobj2<-roc(gja$dia,gja$dia_hpv)
rocobj3<-roc(gja$dia,gja$dia_ydj)
rocobj4<-roc(gja$dia,gja$lh)
rocobj5<-roc(gja$dia,gja$loh)

#计算full AUC
auc(rocobj1)
auc(rocobj2)
auc(rocobj3)
auc(rocobj4)
auc(rocobj5)

#绘制曲线
plot(rocobj1,print.auc=TRUE,legacy.axes=T)
plot(rocobj2,print.auc=TRUE,legacy.axes=T)
plot(rocobj3,print.auc=TRUE,legacy.axes=T)
plot(rocobj4,print.auc=TRUE,legacy.axes=T)
plot(rocobj5,print.auc=TRUE,legacy.axes=T)

#将多条ROC曲线放在一张图上
plot(rocobj1,col="blue",legacy.axes=T)#绘制第一条曲线
lines(rocobj2,col="red")#添加第二条
text(0.3,0.7,"AUC:0.679",col="red")#添加文本（位置，文本，颜色）
lines(rocobj3,col="brown