如何在R中填寫ROC圖的AUC？

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

從pROC的plot.roc功能已經建立參數，以使（auc.polygon）和曲調（auc.polygon.col，auc.polygon.border等，見?plot.roc）的AUC的顯示。在你的情況下，你可以使用：

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)