R логистическая область регрессии под кривой

Question

Выполняю логистическую регрессию с использованием этого page. Мой код выглядит следующим образом.

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") 
mylogit <- glm(admit ~ gre, data = mydata, family = "binomial") 
summary(mylogit) 
prob=predict(mylogit,type=c("response")) 
mydata$prob=prob 

После запуска этого кода mydata dataframe имеет две колонки - «признать» и «prob». Должны ли эти два столбца быть достаточными для получения кривой ROC?

Как получить кривую ROC.

Во-вторых, при поиске в mydata, кажется, что модель прогнозирует вероятность admit=1.

Верно ли это?

Как определить, какое конкретное событие предсказывает модель?

Благодаря

UPDATE: кажется, что ниже трех команд очень полезны. Они обеспечивают отсечку, которая будет иметь максимальную точность, а затем поможет получить кривую ROC.

coords(g, "best") 

mydata$prediction=ifelse(prob>=0.3126844,1,0) 

confusionMatrix(mydata$prediction,mydata$admit 

Answer 1

Кривая ROC сравнивает ранг прогнозирования и ответа. Таким образом, можно оценить кривую ROC с пакетом pROC следующим образом:

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") 
mylogit <- glm(admit ~ gre, data = mydata, family = "binomial") 
summary(mylogit) 
prob=predict(mylogit,type=c("response")) 
mydata$prob=prob 
library(pROC) 
g <- roc(admit ~ prob, data = mydata) 
plot(g)  

Answer 2

Другой способ построения кривой ROC ...

library(Deducer) 
modelfit <- glm(formula=admit ~ gre + gpa, family=binomial(), data=mydata, na.action=na.omit) 
rocplot(modelfit) 

Answer 3

#Another way to plot ROC 

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") 
mylogit <- glm(admit ~ gre, data = mydata, family = "binomial")  
summary(mylogit)  
prob=predict(mylogit,type=c("response"))  
library("ROCR")  
pred <- prediction(prob, mydata$admit)  
perf <- performance(pred, measure = "tpr", x.measure = "fpr")  
plot(perf, col=rainbow(7), main="ROC curve Admissions", xlab="Specificity", 
    ylab="Sensitivity")  
abline(0, 1) #add a 45 degree line