如何生成阳性预测值 (PPV) 与各种分类截止点的关系图？

我生成了一些分数来帮助预测某些内容是“是”(1) 还是“否”(0)，假设数据包括：

scores = c(10:20)

response = c(0,0,1,0,1,0,1,1,0,1,1)

mydata = data.frame(scores, response)

我可以进行 ROC 分析，得出 AUC 为 0.77：

roc(response = mydata$response, predictor = mydata$scores) 

现在，我如何准确地看到选择不同的截止值时会发生什么？我想在 x 轴上设置截止值（比方说 13,14,15,16,17），在 y 轴上设置 PPV。有什么好的方法可以做到这一点？我需要什么功能/包？

我将根据 pROC 包*给出答案。使用 ROCR 包也可以获得类似的结果。

您想使用coords https://www.rdocumentation.org/packages/pROC/versions/1.9.1/topics/coords函数，它可以在某些给定阈值下计算一些常见的统计数据。例如，为了获得所有阈值的 PPV，您可以执行以下操作：

library(pROC)
r <- roc(response = response, predictor = scores)
coordinates <- coords(r, x = "all", input = "threshold", ret = c("threshold", "ppv"))

然后您可以绘制这些值：

plot(t(coordinates))

Replace "all"按兴趣阈值：

 coordinates <- coords(r, x = c(13, 14, 15, 16, 17), input = "threshold", ret = c("threshold", "ppv"))

_{* Disclaimer: I am the author of the pROC package.}