roc<-read.csv("경로/roc1.csv")
head(roc)

#진단검사
install.packages("reportROC")
library(reportROC)

reportROC(roc$smoking,roc$test)

#결과 확인
roc$Nicotinetest <- factor(roc$test>0.5,levels=c(TRUE,FALSE),labels=c("smoking","non-smoking"))
roc$smoking<-factor(roc$smoking>=1,levels=c(TRUE,FALSE),labels=c("smoking","non-smoking"))
table(roc$smoking,roc$Nicotinetest)

install.packages("Epi")
library(Epi)

result=ROC(form= smoking~test,data=roc,plot="ROC")