提取用于计算R中ROC的阈值或截止水平（“Epi”）

Question

有些解决方案可以在packages pROC and ROCR中提取这些值，但我不知道如何对Epi包执行相同的操作。

以下是我正在尝试解决的例子：

library(pROC)
library(Epi)
data(aSAH)
ROC(form = outcome ~ s100b, data=aSAH, plot = "ROC", MX = T)

我尝试了以下内容：

    rc <- ROC(form = outcome ~ s100b, data=aSAH, plot="sp" )
    rc$lr$weights
    str(rc$lr$weights)
 Named num [1:113] 0.185 0.19 0.171 0.143 0.185 ...
 - attr(*, "names")= chr [1:113] "29" "30" "31" "32" ..

这可能是正确的，因为documentation for the package将weights链接到模型矩阵，但我不确定，特别是因为他weights似乎从29开始。< / p>