R ROC, AUC

source("http://bioconductor.org/getBioC.R")
getBioC("ROCR")
getBioC("GEOquery")
install.packages("e1071")

http://rocr.bioinf.mpi-sb.mpg.de/ 
data(ROCR.simple) pred <- ROCR::prediction(ROCR.simple$predictions, ROCR.simple$labels)
perf <- performance(pred,"tpr","fpr")
plot(perf,colorize=TRUE)
plot(performance(prediction(c(0,1,-1), labels=c("a","b","a")), "tpr", "fpr"))
performance(prediction(c(0,1,1), labels=c("a","b","a")), "auc")@y.values