matlab 的plotroc函数怎么调用

ROC曲线是通用的分类器评价工具，matlab函数中自带了绘制该曲线的函数plotroc。转载

plotroc函数的原型为：plotroc(targets, outputs)

其中参数targets是一个矩阵，代表测试集，每一列表示一个测试样本的标签

如果有两类样本，比如第1,2,5个样本属于第1类，第3,4,6个样本属于第2类....则targets应为：

1 1 0 0 1 0 ...

0 0 1 1 0 1 ...

如果只有一类样本，包含了负样本，则只要一行，用1表示正样本，0表示负样本即可,比如targets为：

1 0 1 1 0 0 0 0 1 ...

参数outputs也是一个矩阵，代表分类结果，同样每一列表示一个测试样本的分类结果

同样如果有两类样本，则应有两个分类器，每一列记录了每个测试样本在两个分类器上的得分，此时outputs为：

0.8 0.85 0.2 0.75 0.21 ...

0.8 0.01 0.9 0.23 0.67 ...

如果只有一类，则outputs只有一行，如：

0.8 0.6 0.8 0.7 0.05 0.3 0.03 ...

注意，得分必须在[0, 1]的区间内，可以自己规约一下。

我们将相应的测试标签targets和对应的分类得分outputs输入plotroc中就可以绘制出相应的ROC曲线了。

有人问起，我也就在网上搜了一下，发现还有很多人不会用，写下来以供参考，欢迎指正。

补记：似乎使用matlab中的plot()与roc()组合也能完成ROC曲线绘制。

ROC曲线（Receiver Operating Characteristic Curve）是利用Classification模型真正率（True Positive Rate）和假正率（False Positive Rate）作为坐标轴，图形化表示分类方法的准确率的高低。

ROC图的一些概念定义:：

真正（True Positive , TP）被模型预测为正的正样本

假负（False Negative , FN）被模型预测为负的正样本

假正（False Positive , FP）被模型预测为正的负样本

真负（True Negative , TN）被模型预测为负的负样本

真正率（TPR）

TPR = TP /（TP + FN）

正样本预测结果数 / 正样本实际数 。在ROC曲线中，TPR作为Y轴

假正率（ FPR）

FPR = FP /（FP + TN）

被预测为正的负样本结果数 /负样本实际数 。在ROC曲线中，FPR作为X轴

关于ROC曲线的相关概念的解释，请参阅cowboy_wz的博客ROC曲线：http://blog.csdn.net/chl033/article/details/4507544

我在此主要做的事情是画出ROC曲线工程实现方面的一些解释。我们设计一个函数，此函数需要有一个模型预测值predict和数据标签值ground_truth作为输入参数。分为几步进行实现。

（1）统计数据标签值ground_truth(及y)中分类为0和分类为1的数据数目：

pos_num=sum(ground_truth==1)neg_num=sum(ground_truth)

(2)对利用模型求出的预测值predict由低到高进行排序对应数据原来所在位置进行索引记录，用于重新排序ground_truth.利用函数sort实现，sort详情请查看help文档：

[pre,Index]=sort(predict)ground_truth=ground_truth(Index)

(3)对ground_truth和predict遍历i=1:n，n是测试集数目。其目的是随着predict中概率的增加，随着增加判断正负样本的阈值；也就是说取遍历到的predict值为阈值，大于阈值的假设预测为正样本（阈值右边），小于阈值的假设预测为负样本（阈值左边）。

所以同时我们可得到真正TP 和假正FP值：TP=sum(ground_truth(i:n)==1)FP=sum(ground_truth(i:n)==0)

这时我们就可以求取TPR 和FPR了：TPR=TP/pos_num FPR=FP/neg_num 把求取到的值保存起来（x(i),y(i)）,因为这就是我们要在图上画的点。

（4）返回曲线与坐标轴间的面积auc。我们的目的是测量数据的准确率，这个面积就是一个量度，auc越大，准确率越高。

auc=auc+(y(i)+y(i-1))*(x(i-1)-x(i))/2；

实现ROC曲线的代码如下：

% predict - 分类器对测试集的分类结果

% ground_truth - 测试集的正确标签,这里只考虑二分类，即0和1

% auc- 返回ROC曲线的曲线下的面积

function auc = plot_roc( predict, ground_truth )

%初始点为（1.0, 1.0）

%计算出ground_truth中正样本的数目pos_num和负样本的数目neg_num

pos_num = sum(ground_truth==1)

neg_num = sum(ground_truth==0)

m=size(ground_truth,1)

[pre,Index]=sort(predict)

ground_truth=ground_truth(Index)

x=zeros(m+1,1)

y=zeros(m+1,1)

auc=0

x(1)=1y(1)=1

for i=2:m

TP=sum(ground_truth(i:m)==1)FP=sum(ground_truth(i:m)==0)

x(i)=FP/neg_num

y(i)=TP/pos_num

auc=auc+(y(i)+y(i-1))*(x(i-1)-x(i))/2

end

x(m+1)=0y(m+1)=0

auc=auc+y(m)*x(m)/2

plot(x,y)

end

---