# dataset is assigned to jv # install.packages("pROC") library(pROC) rocjv <- roc(jv$intubated, jv$wbc1) coords(rocjv, "best") coords(rocjv, x = "best", input = "threshold", best.method = "youden")