我有这个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
。
如何确定长度?
在plot.roc
从功能pROC
已经建立参数,以使(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))
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句