码迷,mamicode.com
首页 > 其他好文 > 详细

ROC与AUC曲线绘制

时间:2019-07-12 15:44:17      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:windows   png   lis   rgs   高度   inf   ted   面积   int   

由于ROC曲线面积比较难求得,所以判断模型好坏一般使用AUC曲线

技术图片

 

关于AUC曲线的绘制,西瓜书上写得比较学术,不太能理解,假设有这么一个样本集:

技术图片

假设预测样本为20个,预测为正类的概率已经进行了排序,得分递减,画图步骤为:

(1) 在所排序的样本最左边,画一条线即  无 | 1 2 3 4 5 …,线左边的认为是正类,右边认为是负类,可以算出,TP(实际为正,预测为正)=0,FN(实际为正,预测为负)=10,TN(实际为负,预测为负)=10,FP(实际为负,预测为正)=0,故而TPR=0,FPR=0

(2) 在第一个样本后面,画一条线即 1 | 2 3 4 …,线左边的认为是正类,右边认为是负类,可以算出,TP(实际为正,预测为正)=1,FN(实际为正,预测为负)=9,TN(实际为负,预测为负)=10,FP(实际为负,预测为正)=0,故而TPR=0.1,FPR=0

(3) 在第二个样本后面,画一条线即 1 2| 3 4 …,线左边的认为是正类,右边认为是负类,可以算出,TP(实际为正,预测为正)=2,FN(实际为正,预测为负)=8,TN(实际为负,预测为负)=10,FP(实际为负,预测为正)=0,故而TPR=0.2,FPR=0

(4) 在第三个样本后面,画一条线即 1 2 3| 4 5…,线左边的认为是正类,右边认为是负类,可以算出,TP(实际为正,预测为正)=2,FN(实际为正,预测为负)=8,TN(实际为负,预测为负)=10,FP(实际为负,预测为正)=1,故而TPR=0.2,FPR=0.1

... ...

依次例推:起点一定为(0,0),终点一定为(1,1),由此即可画出AUC曲线

但这样做比较繁琐,可以总结规律为:

按照预测为正类的概率值由大到小排序,首先起点一定是(0,0),然后使用上述划线法进行计算,每次碰到一个预测为正的样本,则在Y方向加一个Y步长,否则在X方向加一个X步长,其中X和Y的步长,计算方法为X: 负类样本个数分之1,Y:正类样本个数分之1,如果是降序排列,则是减一个对应方向的步长即可

代码部分:

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.font_manager import FontProperties

def plotROC(predStrengths,classLabels):
    font = FontProperties(fname = r"c:\windows\fonts\simsun.ttc",size=14)
    cur=(1.0,1.0)             #绘制光标位置
    ySum=0.0                  #用于计算AUC
    numPosClas=np.sum(np.array(classLabels)==1.0)  #正类样本的个数
    yStep = 1/float(numPosClas)                    #Y轴步长
    xStep = 1/float(len(classLabels)-numPosClas)   #X轴步长

    sortedIndicies = predStrengths.argsort()
    fig =plt.figure()
    fig().clf()
    ax=plt.figure()
    for inddex in sortedIndicies.tolist()[0]:
        if classLabels[index]==1.0
            delX=0; delY = yStep   #碰到一个正类,Y的步长减 1
        else:
            delX = xStep; delY=0   #碰到一个负类,X的步长减1
            ySum += cur[1]         #每次cur[1]都可能变化了       #高度累加
        ax.plot([cur[0],cur[0]-delX],[cur[1],cur[1]-delY],c = b)   #绘制ROC
        cur = (cur[0]-delX,cur[1]-delY)                         #更新绘制光标的位置
    ax.plot([0,1],[0,1],b--)
    plt.title(ROC曲线,FontProperties =font)
    plt.xlabel(假正例率,FontProperties = font)
    plt.ylabel((真正例率,FontProperties = font))
    ax.axis([0,1,0,1])
    print(AUC面积:,ySum*xStep)                         #计算AUC
    plt.show()

 

ROC与AUC曲线绘制

标签:windows   png   lis   rgs   高度   inf   ted   面积   int   

原文地址:https://www.cnblogs.com/roscangjie/p/11173900.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!