r - 如何在逻辑回归中以数字的形式获得ROC的最佳分界点

Question

我想在逻辑回归中获得ROC的最佳分界点，而不是两个交叉曲线。
使用下面的代码，我可以得到显示最佳点的图，但是在某些情况下，我只需要将该点作为可用于其他计算的数字即可。
以下是代码行:

library(Epi)
ROC( form = IsVIP ~ var1+var2+var3+var4+var5, plot="sp", data=vip_data ) 

谢谢

Answer 1

根据文档，最佳截止点定义为灵敏度+特异性最大的点(请参见MX中的?ROC参数)。您可以按以下方式获得相应的值(请参阅?ROC中的示例):

x <- rnorm(100)
z <- rnorm(100)
w <- rnorm(100)
tigol <- function(x) 1 - (1 + exp(x))^(-1)
y <- rbinom(100, 1, tigol(0.3 + 3*x + 5*z + 7*w))
rc <- ROC(form = y ~ x + z, plot="sp") 
## optimal combination
opt <- which.max(rowSums(rc$res[, c("sens", "spec")]))
## optimal cut-off point 
rc$res$lr.eta[opt]

这是您运行时将显示的要点

ROC(form = y ~ x + z, plot = "ROC", MX = TRUE)