Как заполнить AUC участка ROC в R?

Question

У меня есть этот ROC сюжет:

library(Epi) 
library(pROC) 
data(aSAH) 
data <- aSAH 
plot.roc(data$outcome, data$s100b) 

Я хочу, чтобы цвет Area under the curve: 0.7314.

Я попытался вариации на ...

x <- seq(1, 0, by = - .001) 
polygon(x, roc(data$outcome, data$s100b), 
     col = rgb(.35,0.31,0.61, alpha = 0.4), 
     border = rgb(.35,0.31,0.61, 0.4), 
     lwd=2) 

... и получаю сообщение об ошибке: Error in xy.coords(x, y) : 'x' and 'y' lengths differ.

Как можно отрегулировать длину?

Answer 1

plot.roc функция от pROC имеет встроенные параметры для включения (auc.polygon) и настройки (auc.polygon.col, auc.polygon.border и т.д., см ?plot.roc) отображение ППК. В вашем случае вы можете использовать:

plot.roc(data$outcome, data$s100b, 
    auc.polygon = TRUE, 
    auc.polygon.col=rgb(.35,0.31,0.61, alpha = 0.4), 
    auc.polygon.border=rgb(.35,0.31,0.61, 0.4)) 

Answer 2

Попробуйте

library(Epi) 
library(pROC) 
data(aSAH) 
data <- aSAH 
res <- plot.roc(data$outcome, data$s100b) 
polygon(with(res, cbind(specificities, sensitivities)), 
     col = rgb(.35,0.31,0.61, alpha = 0.4), 
     border = rgb(.35,0.31,0.61, 0.4), 
     lwd=2)