原文链接：http://tecdat.cn/?p=15508

绘制ROC曲线通过Logistic回归进行分类

加载样本数据。

load fisheriris

通过使用与versicolor和virginica物种相对应的度量来定义二元分类问题。

pred = meas(51:end,1:2);

定义二进制响应变量。

resp = (1:100)'>50;  % Versicolor = 0, virginica = 1

拟合逻辑回归模型。

mdl = fitglm(pred,resp,'Distribution','binomial','Link','logit');

计算ROC曲线。使用逻辑回归模型中的概率估计值作为得分。

perfcurve 将阈值存储在数组中。

显示曲线下的面积。

1.  
   AUC
2.  
   AUC = 0.7918

曲线下的面积为0.7918。最大AUC为1，对应于理想分类器。较大的AUC值表示更好的分类器性能。

1.  
   plot(X,Y)
2.  
   xlabel('False positive rate')
3.  
   ylabel('True positive rate')
4.  
   title('ROC for Classification by Logistic Regression')

使用ROC曲线比较分类方法

加载样本数据

load ionosphere

X 是351x34预测变量的矩阵。 Y 是类别标签的字符数组：  'b' 不良雷达回波和  'g' 良好雷达回波。

重新格式化因变量以适合逻辑回归。

拟合一个逻辑回归模型来估计雷达返回的后验概率是一个不好的概率。

1.  
   mdl = fitglm(pred,resp,'Distribution','binomial','Link','logit');
2.  
   score_log = mdl.Fitted.Probability; % Probability estimates

使用得分的概率计算标准ROC曲线。

在相同的样本数据上训练SVM分类器标准化数据。

mdlSVM = fitcsvm(pred,resp,'Standardize',true);

计算后验概率。

第二列  score_svm 包含不良雷达收益的后验概率。

使用SVM模型的分数计算标准ROC曲线。

在同一样本数据上拟合朴素贝叶斯分类器。

计算后验概率（分数）

[~,score_nb] = resubPredict(mdlNB);

使用朴素贝叶斯分类的分数计算标准ROC曲线。

将ROC曲线绘制在同一张图上。

尽管对于较高的阈值，SVM可以产生更好的ROC值，但逻辑回归通常更擅长区分不良雷达收益与良好雷达。朴素贝叶斯的ROC曲线通常低于其他两个ROC曲线，这表明样本内性能比其他两个分类器方法差。

比较所有三个分类器的曲线下面积。

1.  
   AUClog
2.  
   AUClog = 0.9659
3.  
   AUCsvm
4.  
   AUCsvm = 0.9489
5.  
   AUCnb
6.  
   AUCnb = 0.9393

Logistic回归的AUC度量最高，而朴素的贝叶斯则最低。该结果表明，逻辑回归对此样本数据具有更好的样本内平均性能。

确定自定义内核功能的参数值

本示例说明如何使用ROC曲线为分类器中的自定义内核函数确定更好的参数值。

在单位圆内生成随机的一组点。

定义预测变量。将第一象限和第三象限中的点标记为属于正类别，而将第二象限和第二象限中的点标记为负类。

1.  
   pred = [X1; X2];
2.  
   resp = ones(4*n,1);
3.  
   resp(2*n + 1:end) = -1; % Labels

创建函数mysigmoid.m ，该函数 接受要素空间中的两个矩阵作为输入，并使用S形内核将其转换为Gram矩阵。

使用Sigmoid内核函数训练SVM分类器。使用标准化数据。

设置  gamma = 0.5 ，使用调整后的S形核训练SVM分类器。

1.  
    
2.  
    
3.  
   SVMModel2 = fitPosterior(SVMModel2);
4.  
   [~,scores2] = resubPredict(SVMModel2);

计算两个模型的ROC曲线和曲线下面积（AUC）。

绘制ROC曲线。

1.  
   plot(x1,y1)
2.  
   hold on
3.  
   plot(x2,y2)
4.  
   hold off
5.  
    
6.  
   title('ROC for classification by SVM');

将gamma参数设置为0.5的内核函数可提供更好的样本内结果。

比较AUC度量。

1.  
   auc1
2.  
   auc2
3.  
   auc1 =
4.  
    
5.  
   0.9518
6.  
    
7.  
    
8.  
   auc2 =
9.  
    
10.  
    0.9985

伽玛设置为0.5时曲线下的面积大于伽玛设置为1时曲线下的面积。这也证实了伽玛参数值为0.5会产生更好的结果。为了直观比较这两个伽玛参数值的分类性能。

绘制分类树的ROC曲线

加载样本数据。

load fisheriris

列向量  species由三种不同物种的鸢尾花组成。双矩阵  meas 包含对花朵的四种测量类型：萼片长度，萼片宽度，花瓣长度和花瓣宽度。所有度量单位均为厘米。

使用萼片的长度和宽度作为预测变量训练分类树。

根据树预测物种的分类标签和分数 。

[~,score] = resubPredict(Model);

分数是观察值（数据矩阵中的一行）所属类别的后验概率。列  score 对应于所指定的类  'ClassNames'。

由于这是一个多类问题，因此不能仅将其  score(:,2) 作为输入。这样做将无法提供  perfcurve 有关两个阴性类别（setosa和virginica）分数的足够信息。此问题与二元分类问题不同，在二元分类问题中，知道一个类别的分数就足以确定另一个类别的分数。因此，必须提供  perfcurve 将两个否定类的得分纳入考虑范围的函数。一种函数是score（：，2）-max（score（：，1），score（：，3））。

X，默认为假阳性率，  Y，默认为真阳性率（召回率或敏感性）。正类标签为  versicolor。由于未定义否定类别，因此  perfcurve 假设不属于肯定类别的观测值属于一个类别。该函数将其接受为否定类。

1.  
    
2.  
   suby = 12×2
3.  
    
4.  
   0 0
5.  
   0.1800 0.1800
6.  
   0.4800 0.4800
7.  
   0.5800 0.5800
8.  
   0.6200 0.6200
9.  
   0.8000 0.8000
10.  
    0.8800 0.8800
11.  
    0.9200 0.9200
12.  
    0.9600 0.9600
13.  
    0.9800 0.9800
14.  
    ⋮
15.  
     
16.  
     
17.  
    subnames = 1x2 cell
18.  
    {'setosa'} {'virginica'}

在ROC曲线上绘制ROC曲线和最佳工作点。

找到与最佳工作点相对应的阈值。

1.  
   T((X==OPTROCPT(1))&(Y==OPTROCPT(2)))
2.  
   ans = 0.2857

指定  virginica 为否定类，并计算和绘制ROC曲线  versicolor。

同样，必须提供  perfcurve 将否定类分数纳入考量的函数。要使用的函数的一个示例是score（：，2）-score（：，3）。

计算ROC曲线的逐点置信区间

加载样本数据。

load fisheriris

仅将前两个变量用作预测变量，来定义二元问题。

pred = meas(51:end,1:2);

定义二进制因变量。

resp = (1:100)'>50;  % Versicolor = 0, virginica = 1

拟合逻辑回归模型。

通过垂直平均（VA）和使用bootstrap进行采样，计算真实正率（TPR）上的逐点置信区间。

'NBoot',1000 将引导样本的数量设置为1000。  'XVals','All' 提示  perfcurve 返回  X，  Y和  T 所有分数的Y 值，并X 使用垂直平均将所有值的值（真阳性率）  平均  （假阳性率）。 默认情况下将使用阈值平均来计算置信范围。

绘制逐点置信区间。

errorbar(X,Y(:,1),Y(:,1)-Y(:,2),Y(:,3)-Y(:,1));

不一定总是可以控制误报率（FPR，X 此示例中的  值）。因此，可能希望通过阈值平均来计算真实正利率（TPR）的逐点置信区间。

绘制置信区间。

1.  
   figure()
2.  
   errorbar(X1(:,1),Y1(:,1),Y1(:,1)-Y1(:,2),Y1(:,3)-Y1(:,1));

指定阈值计算ROC曲线。然后绘制曲线。

1.  
    
2.  
   figure()
3.  
   errorbar(X1(:,1),Y1(:,1),Y1(:,1)-Y1(:,2),Y1(:,3)-Y1(:,1));

参考文献

最受欢迎的见解

1.R语言多元Logistic逻辑回归 应用案例

2.面板平滑转移回归(PSTR)分析案例实现

3.matlab中的偏最小二乘回归（PLSR）和主成分回归（PCR）

4.R语言泊松Poisson回归模型分析案例

5.R语言回归中的Hosmer-Lemeshow拟合优度检验

6.r语言中对LASSO回归，Ridge岭回归和Elastic Net模型实现

7.在R语言中实现Logistic逻辑回归

8.python用线性回归预测股票价格

9.R语言如何在生存分析与Cox回归中计算IDI，NRI指标