Как получить кривую ROC в R / plotROC? - PullRequest
Новая
       23

Как получить кривую ROC в R / plotROC?

0 голосов
/ 17 декабря 2018

Я только начал изучать кривых.Я пытаюсь создать ROC plot, но кажется, что кривая изгибается в неправильном направлении, чем обычно - см. Прилагаемый

Можете ли вы помочь повернуть кривую так, чтобы она изгибалась в «обычном» направлении?

enter image description here

Мои данные 

p <- structure(list(t = c(29354L, 7445L, 22309L, 29699L, 29711L, 14765L, 22257L, 
29715L, 29772L, 13320L, 20905L, 12950L, 3400L, 14800L,7400L, 21890L, 19400L, 14800L, 14700L, 22200L, 1688L, 4500L, 8438L, 13500L, 14800L, 
12580L, 12950L, 13320L, 11840L, 13320L, 14800L, 13690L, 11250L, 12210L, 13320L, 13320L, 14800L, 12580L,20720L, 11840L, 14800L, 7030L, 14800L, 
14800L, 8325L, 11100L,10730L, 13690L, 12210L, 14800L), a = c(0L, 1L, 1L, 0L, 0L, 
1L,0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,1L, 0L, 0L, 0L, 
0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L,0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 
0L, 0L, 0L)), .Names = c("t","a"), class = "data.frame", row.names = c(NA, 50L))

И я использовал следующие 

library(plotROC)
basicplot <- ggplot(p, aes(d = p$a, m = p$t)) + geom_roc() + theme_bw()
basicplot + style_roc()

Бонусвопросы

Я использовал Youden's, чтобы определить наилучшее значение отсечения, которое было 13410. Есть идеи, как привлечь особое внимание к этой точке - например, выделив красную точку?

1 Ответ

0 голосов
/ 17 декабря 2018

Вы можете сделать две вещи, поменять местами обозначения 1 и 0 или использовать geom_roc(increasing = FALSE).См. PlotROC виньетка .Это предполагает, что ваша модель действительно работает таким образом, и на самом деле не имеет низкой чувствительности и высокого уровня ложных срабатываний, а вы просто инвертируете ее. 

basicplot <- ggplot(p, aes(d = a, m = t)) + geom_roc(increasing = FALSE) + theme_bw()
basicplot + style_roc()

...