install.packages("pROC")
macs.labels = read.table("../../tests/mhc_test_data/macs_sequences.fasta.roc", sep=T)
macs.labels = read.table("../../tests/mhc_test_data/macs_sequences.fasta.roc", sep=",")
macs.labels
macs.roc = roc.curve(macs.labels)
plot(macs.roc)
macs.roc
macs.roc$TPR
macs.labels[0,]
macs.labels[1,]
labels = as.vector(macs.labels[0,])
labels
labels = as.vector(macs.labels[1,])
labels
labels = as.numeric(macs.labels[1,])
labels
macs.roc = roc.curve(labels)
mcas.roc$TPR
macs.roc$TPR
plot(macs.roc$FPR, macs.roc$TPR)
macs.labels = read.table("../../tests/mhc_test_data/macs_sequences.fasta.roc", sep=",")

compare.roc("../../tests/mhc_test_data/macs_sequences.fasta.roc", "../../tests/mhc_test_data/test_sequences.fasta.roc")
source("roc_analyisis.r")
source("roc_analyisis.r")
source("roc_analyisis.r")
compare.roc("../../tests/mhc_test_data/macs_sequences.fasta.roc", "../../tests/mhc_test_data/test_sequences.fasta.roc")
source("roc_analyisis.r")
compare.roc("../../tests/mhc_test_data/macs_sequences.fasta.roc", "../../tests/mhc_test_data/test_sequences.fasta.roc")
source("roc_analyisis.r")
compare.roc("../../tests/mhc_test_data/macs_sequences.fasta.roc", "../../tests/mhc_test_data/test_sequences.fasta.roc")
labels
len(labels)
length(labels)
 library(pROC)
roc(labels, 1:215)
1:215
roc(c(0,1,0,0), c(1, 2, 3, 4))
rocobj <- roc(c(0,1,0,0), c(1, 2, 3, 4))
rocobj <- roc(c(0,1,0,0), c(1, 2, 3, 4))
help(roc)
roc(c(1,1,1,0,0), c(1, 2, 3, 4, 5))
roc(c(1,1,1,0,0), c(1, 2, 3, 4, 5), quiet=FALSE)
roc(response=c(1,1,1,0,0), predictor=c(1, 2, 3, 4, 5))
roc(response=c(1,1,1,0,0), predictor=c(1, 2, 3, 4, 5), auc=FALSE)
library(pROC)
source("roc_analyisis.r")
compare.roc("../../tests/mhc_test_data/macs_sequences.fasta.roc", "../../tests/mhc_test_data/test_sequences.fasta.roc")
roc(response=c(1,1,1,0,0), predictor=c(1, 2, 3, 4, 5), auc=FALSE)
 df = data.frame
df$response = c(1, 1, 1, 0, 0)
df$response = l(1, 1, 1, 0, 0)
c(1, 1,1)
a = c(1, 1, 1)
a
roc(c(1,1,1,0,0) ~c(1, 2, 3, 4, 5), auc=FALSE)
roc(c(1,1,1,0,0)~c(1, 2, 3, 4, 5), auc=FALSE)
roc(c(1,1,1,0,0)~c(1, 2, 3, 4, 5))
c
response = c(1, 0, 1, 0)
predictor = c(1, 2, 3, 4)
roc(response~predictor)
 response
auc(response~predictor)
quit
library(pROC)
roc(c(1,0,1, 0)~c(1,2,3,4))
install.packages("pROC")
install.packages("plyr")
