其他
一行代码输出ROC曲线的截断值、AUC、准确度、SE、SP、PLR、NLR、PPV、NPV、TPA、KAPPA等详细统计信息
今天来学习一个新R包——reportROC,这个R包可以一行代码输出ROC曲线的各项统计数值及ROC曲线,比如截断值、AUC及置信区间、灵敏度及置信区间、特异度及置信区间、阳性似然比、阴性似然比、阴性预测值等等。
安装R包可以直接从CRAN上安装。
install.packages("reportROC")
library(reportROC)
使用官方的示例进行演示。
data(aSAH)
#for continuous variables
reportROC(gold=aSAH$outcome,
predictor=aSAH$s100b,
important="se",
plot=TRUE)
直接输出ROC曲线的各种统计信息。
可以将上述的信息存储为数据框,然后导出来。
cc <- reportROC(gold=aSAH$outcome,
predictor=aSAH$s100b,
important="se",
plot=TRUE)
write.csv(cc, "cc.csv")
这个R包函数不多,有兴趣的可以具体看下R包帮助文档。
关注下方公众号,分享更多更好玩的R语言知识。
如果需要进群,请在菜单栏处添加作者微信,添加入群。
觉得有帮助的请点赞、分享、在看走起!
点个在看,SCI马上发表。