library(pROC)
data(aSAH)
fit1 <- roc(aSAH$outcome, aSAH$s100b, levels = c("Good", "Poor"))
fit2 <- roc(aSAH$outcome, aSAH$ndka, levels = c("Good", "Poor"))
plot.roc(fit1, print.auc=T, print.auc.x=0.8, print.auc.y=0.8)
plot.roc(fit2, print.auc=T, print.auc.x=0.5, print.auc.y=0.5, add=T, col=2)
bGlicmFyeShwUk9DKQpkYXRhKGFTQUgpCmZpdDEgPC0gcm9jKGFTQUgkb3V0Y29tZSwgYVNBSCRzMTAwYiwgbGV2ZWxzID0gYygiR29vZCIsICJQb29yIikpCmZpdDIgPC0gcm9jKGFTQUgkb3V0Y29tZSwgYVNBSCRuZGthLCBsZXZlbHMgPSBjKCJHb29kIiwgIlBvb3IiKSkKcGxvdC5yb2MoZml0MSwgcHJpbnQuYXVjPVQsIHByaW50LmF1Yy54PTAuOCwgcHJpbnQuYXVjLnk9MC44KQpwbG90LnJvYyhmaXQyLCBwcmludC5hdWM9VCwgcHJpbnQuYXVjLng9MC41LCBwcmludC5hdWMueT0wLjUsIGFkZD1ULCBjb2w9MikK