使用pROC
包,是否有人知道如何通过输入特定阈值而不是特定SP或SN值来获取SP和SN的置信区间? ci()
似乎仅使用其对应的SN或SP作为参数显示SP和SN值;我可以在SN为0.8时要求SP的CI,但对于阈值为9.0的SP的CI则不能。
我知道我可以找到阈值的SP / SN值,然后将它们用作ci()
的参数,但这看起来很笨拙。当然,这可能最终成为包开发人员的问题。
提前感谢您提供的任何帮助。
答案 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