(ROCR tutorial)[https://www.r-bloggers.com/a-small-introduction-to-the-rocr-package/] prediction=list(data[1:16517,”SCORESUM”],data[16518:(165172),”SCORESUM”],data[(165172+1):(165173),”SCORESUM”],data[(165173+1):(165174),”SCORESUM”],data[(165174+1):82581,”SCORESUM”]) labels=list(data[1:16517,”myopia”],data[16518:(165172),”myopia”],data[(165172+1):(165173),”myopia”],data[(165173+1):(165174),”myopia”],data[(165174+1):82581,”myopia”]) pred=prediction(prediction, labels) pred
perf<- performance(pred, “tpr”, “fpr”) plot(perf,col=as.list(rainbow(5))) #col=1:10 abline(a=0,b=1)
#area under the curve auc.perf = performance(pred, measure = “auc”) auc.perf@y.values
#
