python计算ROC曲线和面积AUC

ROC曲线是根据一系列不同的二分类方式(分界值或决定阈),以真正率(也就是灵敏度)(True Positive Rate,TPR)为纵坐标,假正率(1-特效性)(False Positive Rate,FPR)为横坐标绘制的曲线。通过将连续变量设定出多个不同的临界值,从而计算出一系列敏感性和特异性,从而可以绘制ROC曲线。

纵坐标:真正率(True Positive Rate , TPR)或灵敏度(sensitivity)

TPR = TP /(TP + FN) (正样本预测结果数 / 正样本实际数)

横坐标:假正率(False Positive Rate , FPR)

FPR = FP /(FP + TN) (被预测为正的负样本结果数 /负样本实际数)

利用sklearn.metrics.roc_curve可以计算ROC曲线

from sklearn.metrics import roc_curve, aucy_true = [0, 1, 1]
y_score = [0.1, 0.8, 0.7]fpr, tpr, thresholds = roc_curve(y_true, y_score)
print(fpr, tpr, thresholds)
"""
[0. 0. 0. 1.] 
[0.  0.5 1.  1. ] 
[1.8 0.8 0.7 0.1]
"""

其中y_true是真实标签,y_score是预测概率,fpr是假正率,tpr是真正率,thresholds是阈值(thresholds[0]=max(y_score)+1=1.8)。thresholds是y_score去重复元素后加上1.8所以一共4个元素。

在这个例子中,正样本实际数(TP + FN)=2,负样本实际数(FP + TN)=1

  • 当index=0,thresholds[0]=1.8,大于等于1.8的预测为正样本,小于1.8的预测为负样本,此时预测标签为[0,0,0],此时TP=0,FP=0,fpr=0,tpr=0
  • 当index=1,thresholds[1]=0.8,大于等于0.8的预测为正样本,小于0.8的预测为负样本,此时预测标签为[0,1,0],此时TP=1,FP=0,fpr=0,tpr=0.5
  • 当index=2,thresholds[1]=0.7,大于等于0.7的预测为正样本,小于0.7的预测为负样本,此时预测标签为[0,1,1],此时TP=2,FP=0,fpr=0,tpr=1
  • 当index=3,thresholds[1]=0.1,大于等于0.1的预测为正样本,小于0.1的预测为负样本,此时预测标签为[1,1,1],此时TP=2,FP=1,fpr=1,tpr=1

计算出ROC曲线后,可以利用sklearn.metrics.auc计算AUC:

from sklearn.metrics import roc_curve, aucy_true = [0, 1, 1]
y_score = [0.1, 0.8, 0.7]fpr, tpr, thresholds = roc_curve(y_true, y_score)
print(auc(fpr, tpr))

 

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/490792.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

python from __future__ import division

python from __future__ import division 之前一直很困惑,为什么这个模块叫future呢,难道有什么特殊功能能够让人们想到未来吗,最近才恍然大悟。 python的更新和前进是由社区进行推动的,而且是免费开源的,不受大型…

【VS开发】CTimeSpan类

CTimeSpan类。 日期和时间类简介 CTime类的对象表示的时间是基于格林威治标准时间(GMT)的。CTimeSpan类的对象表示的是时间间隔。 CTime类和CTimeSpan类一般不会被继承使用。两者对象的大小都是8个字节。 CTime表示的日期上…

搅动世界的两大因素

原创:张晓峰提要:移动互联、云计算、大数据、人工智能等技术因素逐步成为新基础设施,而泛连接、泛共享、泛融合与泛协同为代表的非技术因素正在重构这个世界。二者叠加融汇、相因相生。每个人都渐进或主动或被动地“被”函数化、数字化、孪生…

python计算PR曲线sklearn.metrics.precision_recall_curve

PR曲线实则是以precision(精准率)和recall(召回率)这两个为变量而做出的曲线,其中recall为横坐标,precision为纵坐标。设定一系列阈值,计算每个阈值对应的recall和precision,即可计算…

amigo幸运字符什么意思_转载 | 史上最全 python 字符串操作指南

点击蓝字关注,创智助你长姿势【本文已由 清风Python 授权转载(原创)作者:王翔,转载请联系出处】字符串的定义完了,估计很多人看到这个标题就要关网页了,稍等不妨再往下看看?python 定义字符、字符串没有 j…

在物理学的语言里,“生命”是什么?

转自:Darthusian“想象一种语言就像想象一种形式的生命。”--- 路德维希.维特根斯坦,《哲学研究》当今世界人们使用大约6,800种不同的语言。不是每个词都能在不同的语言之间完美地翻译,意义有时会落入语义的裂缝。例如,日语词wabi…

HTML5中的时间类型,另外EL表达式的时间值来读取时间,并且还可以更改时间

HTML5规范里只规定date新型input输入类型,并没有规定日历弹出框的实现和样式。所以,各浏览器根据自己的设计实现日历。目前只有谷歌浏览器完全实现日历功能。相信这种局面很快就会结束,所有的浏览器最终都将会提供原生的日历组件。 如果你使用…

python 中 yield 的使用

python中yield在函数中的使用 正是因为函数含有这个yield,所以,该函数不再是普通的函数,而是生成器函数(generator function)。下面通过小例子来说明一下这个内置函数的特性: 1、直接运行函数并不会执行。…

支付接口申请_2020新版去水印小程序搭建教程(附源码+对接支付)

可以设置小程序跳转,引流也不错支持小程序后台流量主激励视频广告用户下载一次观看一次视频广告,收入非常可观支持小程序后台流量主banner广告支持全网短视频解析,苹果安卓通用支持后台无限生成卡密支持自定义文字广告支持图片广告支持小程序…

python 绘制功率谱密度图 plt.psd

功率谱密度图以横轴为频率,纵轴为功率密度,表示信号功率密度随着频率的变化情况 python绘制功率谱密度: matplotlib.pyplot.psd(x, NFFT256, Fs2, Fc0, detrendmlab.detrend_none,windowmlab.window_hanning, noverlap0, pad_toNone,sidesd…

概述自动机器学习(AutoML)

来源:混沌巡洋舰华为年薪百万招聘的机器学习博士,其中一个的研究方向就是AutoML,今天借着最新的综述,来对AutoML进行简要的概述。机器学习是让算法自动的从数据中找出一组规则,从而提取数据中对分类/聚类/决策有帮助的特征&#x…

python assert 与 slice

assert: assert 断言的作用:assert是断言的意思,python assert 断言语句格式及用法很简单。assert断言是声明其布尔值必须为真的判定,如果发生异常就说明表达式为假,用来测试表达式,如果表达式的返回值为假…

linux-ntpdate同步更新时间

Linux服务器运行久时,系统时间就会存在一定的误差,一般情况下可以使用date命令进行时间设置,但在做数据库集群分片等操作时对多台机器的时间差是有要求的,此时就需要使用ntpdate进行时间同步 安装: yum install ntpdat…

卸载52好压,极速输入法,手机模拟大师这些流氓软件

卸载52好压,极速输入法的方法: 通过控制面板里卸载了之后发现这两个流氓软件居然还在?于是换了一种卸载方法: 先在官网下载并重新安装这两个软件,如果电脑已经安装了旧版本会提示安装的新版本会覆盖旧的版本&#xf…

h5 bootstrap 小程序模板_软件测试人员必知H5/小程序测试点

最近接触了较多关于H5页面的测试,H5页面的测试除了业务逻辑功能测试外,其他部分的测试方法基本是可以通用的,在此对H5页面和小程序的一些通用测试方法进行总结分享给大家01H5相关测试H5优势:1.H5可以跨平台,开发成本相对较低&…

【深度学习】这才是深度学习的本源

内容来源:2019年7月26日,在长城会主办的“GMIC之院士AI论坛”上,索尔克生物研究所计算神经生物学实验室主任特伦斯进行了以“人工智能,将成科技的最大风口”为主题的精彩分享。笔记侠作为合作方,未经主办方和讲者审核。…

KALI Linux中GURB安装失败如何处理

KALI Linux中GURB安装失败如何处理 # GRUB安装失败 # 无法将grub-pc软件包安装到/target/中,如果没有GRUB启动引导器,所安装的系统将无法启动。 解决办法: 点击返回 点击配置软件包管理 将之前设置的网络镜像“是”改为“否”&#xff0c…

IEEE选择和下载投稿期刊模板(LaTeX和Word)网站

https://template-selector.ieee.org/ (复制网址在新标签页打开,不要直接点)

d3.js 旋转图形_苏教版三年级数学上册6.1平移和旋转微课视频 | 练习

本公众号学习资源丰富!欢迎分享到班群、朋友圈、帮助更多孩子。微课视频第一课时微课视频第二课时同步练习1.连一连。升旗时国旗的运动 钟摆的运动风扇叶片的运动 平移 电梯的运动光盘在电脑里的运动 旋转 火车的…

任正非签发最新电邮:过去我们是为了赚点小钱,现在是要战胜美国

来源:蓝血研究(lanxueyanjiu)2019年7月31日,华为举行“千疮百孔的烂伊尔2飞机”战旗交接仪式,任正非在仪式上做了题为《钢铁是怎么炼成的》的讲话。他表示:人工智能才是未来的大产业,要让“人才…