使用R中的pROC包从特定阈值获取SP和SN的CI

时间:2012-08-28 15:01:58

标签: r roc

使用pROC包,是否有人知道如何通过输入特定阈值而不是特定SP或SN值来获取SP和SN的置信区间? ci()似乎仅使用其对应的SN或SP作为参数显示SP和SN值;我可以在SN为0.8时要求SP的CI,但对于阈值为9.0的SP的CI则不能。

我知道我可以找到阈值的SP / SN值,然后将它们用作ci()的参数,但这看起来很笨拙。当然,这可能最终成为包开发人员的问题。

提前感谢您提供的任何帮助。

1 个答案:

答案 0 :(得分:4)

当首次创建roc对象时,可以计算指定阈值的SN和SP的置信区间,以及ci.thresholds(rocobj, thresholds=0.9)

使用pROC的帮助文件中的示例数据:

library(pROC)
data(aSAH)

rocobj <- roc(aSAH$outcome, aSAH$s100b, ci=TRUE, of="thresholds", thresholds=0.9)

95% CI (2000 stratified bootstrap replicates):
 thresholds sp.low sp.median sp.high se.low se.median se.high
        0.9      1         1       1      0   0.04878   0.122