广州企业建设网站,云南手机网站建设,网站导航栏下面的文章,公司想建个网站文章目录一、多分类评估指标的 macro 和 weighted 过程1. 多分类 F1-Score 评估指标2. 多分类 ROC-AUC 评估指标二、借助机器学习流构建全域参数搜索空间三、优化评估指标选取1. 高级评估指标的选用方法2. 同时输入多组评估指标四、优化后建模流程在正式讨论关于网格搜索的进阶…
文章目录一、多分类评估指标的 macro 和 weighted 过程1. 多分类 F1-Score 评估指标2. 多分类 ROC-AUC 评估指标二、借助机器学习流构建全域参数搜索空间三、优化评估指标选取1. 高级评估指标的选用方法2. 同时输入多组评估指标四、优化后建模流程在正式讨论关于网格搜索的进阶使用方法之前我们需要先补充一些关于多分类问题的评估指标计算过程。在此前的内容中我们曾经介绍过分类模型在解决多分类问题时的不同策略同时也介绍过二分类问题的更高级评估指标如 F1-Score 和 ROC-AUC 等。接下来我们将详细讨论关于多分类预测结果在 F1-Socre 和 ROC-AUC 中的评估过程以及 在sklearn 中如何调用函数进行计算。
# 科学计算模块
import numpy as np
import pandas as pd
# 绘图模块
import matplotlib as mpl
import matplotlib.pyplot as plt
# 自定义模块
from ML_basic_function import *
# Scikit-Learn相关模块
# 评估器类
from sklearn.preprocessing import StandardScaler
from sklearn.preprocessing import PolynomialFeatures
from sklearn.linear_model import LogisticRegression
from sklearn.pipeline import make_pipeline
from sklearn.model_selection import GridSearchCV
# 实用函数
from sklearn.model_selection import train_test_split
from sklearn.metrics import accuracy_score
# 数据准备
from sklearn.datasets import load_iris一、多分类评估指标的 macro 和 weighted 过程
1. 多分类 F1-Score 评估指标
首先导入和 F1-Score 相关的评估指标计算函数。
from sklearn.metrics import precision_score,recall_score,f1_score然后简单查看相关说明文档发现这几组和混淆矩阵相关的评估指标基本是共用了一套参数命名并且大多数参数其实都是作用于多分类问题对于二分类问题我们可以简单调用相关函数直接计算
y_true np.array([1, 0, 0, 1, 0, 1])
y_pred np.array([1, 1, 0, 1, 0, 1])precision_score(y_true, y_pred), recall_score(y_true, y_pred), f1_score(y_true, y_pred)
#(0.75, 1.0, 0.8571428571428571)precision_score?具体参数含义解释如下
NameDescriptiony_true数据集真实标签y_pred标签预测结果labels允许以列表形式输入其他形态的标签一般不进行修改pos_labelpositive类别标签average多分类时指标计算方法sample_weight不同类别的样本权重zero_division当分母为0时返回结果
其中需要重点介绍多分类问题时 average 参数不同取值时的计算方法。此处以 recall 为例进行计算重点介绍当 average 取值为 ‘macro’、‘micro’ 和 ‘weighted’ 的情况其他指标也类似例如有简单多分类问题如下 我们令 1 类标签为 0、2 类标签为 1、3 类标签为 2则上述数据集真实标签为
y_true np.array([0, 1, 2, 2, 0, 1, 1, 2, 0, 2])并且最终分类预测结果为
y_pred np.array([0, 1, 0, 2, 2, 1, 2, 2, 0, 2])据此我们可以构造多分类混淆矩阵如下 据此我们可以计算三个类别的 TP 和 FN
tp1 2
tp2 2
tp3 3fn1 1
fn2 1
fn3 1接下来有两种计算 recall 的方法其一是先计算每个类别的 recall然后求均值
re1 2/3
re2 2/3
re3 3/4np.mean([re1, re2, re3])
#0.6944444444444443这也就是 average 参数取值为 macro 时的计算结果
recall_score(y_true, y_pred, averagemacro)
#0.6944444444444443当然如果上述手动实现过程不求均值而是根据每个类别的数量进行加权求和则就是参数 average 参数取值为 weighted 时的结果
re1 * 3/10 re2 * 3/10 re3 * 4/10
#0.7recall_score(y_true, y_pred, averageweighted)
#0.7当然还有另外一种计算方法那就是先计算整体的 TP 和 FN然后根据整体 TP 和 FN 计算 recall
tp tp1 tp2 tp3
fn fn1 fn2 fn3tp / (tpfn)
#0.7该过程也就是 average 参数取值 micro 时的计算结果
recall_score(y_true, y_pred, averagemicro)
#0.7对于上述三个不同参数的选取首先如果是样本不平衡问题如果是要侧重训练模型判别小类样本的能力的情况下、则应排除 weighted 参数以避免赋予大类样本更高的权重。除此以外在大多数情况下这三个不同的参数其实并不会对最后评估器的选取结果造成太大影响只是在很多要求严谨的场合下需要说明多分类的评估结果的计算过程此时需要简单标注下是按照何种方法进行的计算。不过如果是混淆矩阵中相关指标和 ROC-AUC 指标放在一起讨论由于新版 sklearn 中 ROC-AUC 本身不支持在多分类时按照 micro 计算、只支持 macro 计算因此建议混淆矩阵的多分类计算过程也选择 macro 过程以保持一致。不过值得注意的是还有一种观点尽管 micro 和 macro 方法在混淆矩阵相关指标的计算过程中差别不大在 ROC-AUC 中macro 指标并不利于非平衡样本的计算混淆矩阵中可以通过 positive 的类别选择来解决这一问题需要配合 OVR 分类方法才能够有所改善。
2. 多分类 ROC-AUC 评估指标
接下来继续讨论关于多分类的 ROC-AUC 评估指标的相关问题
from sklearn.metrics import roc_auc_score能够发现roc_auc_score 评估指标函数中大多数参数都和此前介绍的混淆矩阵中评估指标类似。接下来我们简单尝试使用 ROC-AUC 函数进行评估指标计算根据 ROC-AUC 的计算流程可知此处我们需要在 y_pred 参数位中输入模型概率预测结果
y_true np.array([1, 0, 0, 1, 0, 1])
y_pred np.array([0.9, 0.7, 0.2, 0.7, 0.4, 0.8])roc_auc_score(y_true, y_pred)
#0.9444444444444444当然如果我们在 y_pred 参数中输入分类结果该函数也能计算出最终结果
y_true np.array([1, 0, 0, 1, 0, 1])
y_pred np.array([1, 1, 0, 1, 0, 1])roc_auc_score(y_true, y_pred)
#0.8333333333333334不过此时模型会默认预测标签为 0 的概率结果为 0.4、预测标签为 1 的概率预测结果为 0.6即上述结果等价于
y_true np.array([1, 0, 0, 1, 0, 1])
y_pred np.array([0.6, 0.6, 0.4, 0.6, 0.4, 0.6])roc_auc_score(y_true, y_pred)
#0.8333333333333334即计算过程会默认模型概率预测结果更差。接下来详细解释 ROC-AUC 中其他参数
roc_auc_score?NameDescriptionmax_fprfpr最大值fpr是roc曲线的横坐标multi_class分类器在进行多分类时进行的多分类问题处理策略
此处需要注意的是关于 multi_class 参数的选择。一般来说 sklearn 中的 multi_class 参数都是二分类器中用于解决多元分类问题时的参数如逻辑回归。由于 ROC-AUC 需要分类结果中的概率来完成最终计算因此需要知道概率结果对应分类标签——即到底是以 OVO 还是 OVR 模式在进行多分类因此如果是进行多分类 ROC-AUC 计算时需要对其进行明确说明。不过对于多分类逻辑回归来说无论是 OVR 还是 MVM 策略最终分类结果其实都可以看成是 OVR 分类结果因此如果是多分类逻辑回归计算 ROC-AUC 需要设置 multi_class 参数为 OVR 。同时由于根据 ROC-AUC 的函数参数说明可知在 multi_class 参数取为 OVR 时average 参数取值为 macro 时能够保持一个较高的偏态样本敏感性因此对于 ROC-AUC 来说大多数时候 average 参数建议取值为 macro。总结一下对于 ROC-AUC 进行多分类问题评估时建议选择的参数组合是 OVR/OVOmacro而 OVR/OVO 的参数选择需要根据具体的多分类模型来定如果是围绕逻辑回归多分类评估器来进行结果评估则建议 ROC-AUC 和逻辑回归评估器的 multi_class 参数都选择 OVR。在新版的 sklearn 中 ROC-AUC 函数的multi_class参数已不支持 micro 参数面对多分类问题该参数只能够在 macro 和 weighted 中进行选择。接下来我们简单测算 average 参数中 macro 和 weighted 的计算过程。还是围绕上述数据集进行计算 据此我们可以计算每个类别单独的 ROC-AUC 值
y_true_1 np.array([1, 0, 0, 0, 1, 0, 0, 0, 1, 0])
y_pred_1 np.array([0.8, 0.2, 0.5, 0.2, 0.3, 0.1, 0.3, 0.3, 0.9, 0.3])r1 roc_auc_score(y_true_1, y_pred_1)
r1
#0.8809523809523809y_true_2 np.array([0, 1, 0, 0, 0, 1, 1, 0, 0, 0])
y_pred_2 np.array([0.2, 0.6, 0.3, 0, 0.2, 0.8, 0.2, 0.3, 0, 0.1])r2 roc_auc_score(y_true_2, y_pred_2)
r2
#0.8571428571428571y_true_3 np.array([0, 0, 1, 1, 0, 0, 0, 1, 0, 1])
y_pred_3 np.array([0, 0.2, 0.2, 0.8, 0.5, 0.1, 0.5, 0.4, 0.1, 0.6])r3 roc_auc_score(y_true_3, y_pred_3)
r3
#0.8125此时 r1、r2、r3 的均值如下
np.mean([r1, r2, r3])
#0.8501984126984127该结果应当和 macromulti_class 参数计算结果相同
y_pred np.concatenate([y_pred_1.reshape(-1, 1), y_pred_2.reshape(-1, 1), y_pred_3.reshape(-1, 1)], 1)
y_pred
#array([[0.8, 0.2, 0. ],
# [0.2, 0.6, 0.2],
# [0.5, 0.3, 0.2],
# [0.2, 0. , 0.8],
# [0.3, 0.2, 0.5],
# [0.1, 0.8, 0.1],
# [0.3, 0.2, 0.5],
# [0.3, 0.3, 0.4],
# [0.9, 0. , 0.1],
# [0.3, 0.1, 0.6]])y_true np.array([0, 1, 2, 2, 0, 1, 1, 2, 0, 2])roc_auc_score(y_true, y_pred, averagemacro, multi_classovr)
#0.8501984126984127当然如果 ROC-AUC 函数的参数是 OVRweighted则计算结果过程验证如下
r1 * 3/10 r2 * 3/10 r3 * 4/10
#0.8464285714285713roc_auc_score(y_true, y_pred, averageweighted, multi_classovr)
#0.8464285714285713至此我们就能够较为清楚的了解关于 F1-Score 和 ROC-AUC 评估指标在调用 sklearn 中相关函数解决多分类问题评估的具体方法。
二、借助机器学习流构建全域参数搜索空间
在 Lesson 6.5 中我们已经完整总结了机器学习调参的基本理论同时介绍了 sklearn 中网格搜索GridSearchCV评估器的参数及基本使用方法。我们将进一步介绍网格搜索的进阶使用方法并同时补充多分类问题评估指标在 sklearn 中实现的相关方法然后围绕 Lesson 6.4 中提出的问题给出一个基于网格搜索的解决方案。首先是关于评估器全参数的设置方法。在此前的实验中我们只是保守的选取了部分我们觉得会对模型产生比较大影响的超参数来构建参数空间但在实际场景中调参应该是纳入所有对模型结果有影响的参数进行搜索、并且是全流程中的参数来进行搜索。也就是说我们设置参数的空间的思路不应该更加“激进”一些首先是对逻辑回归评估器来说应该是排除无用的参数外纳入所有参数进行调参并且就逻辑回归模型来说往往我们需要在模型训练前进行特征衍生以增强模型表现。因此我们应该先构建一个包含多项式特征衍生的机器学习流、然后围绕这个机器学习流进行参数搜索这才是一个更加完整的调参过程。首先仿造 Lesson 6.4 中展示过程创造数据集如下
np.random.seed(24)
X np.random.normal(0, 1, size(1000, 2))
y np.array(X[:,0]X[:, 1]**2 1.5, int)np.random.seed(24)
for i in range(200):y[np.random.randint(1000)] 1y[np.random.randint(1000)] 0plt.scatter(X[:, 0], X[:, 1], cy)X_train, X_test, y_train, y_test train_test_split(X, y, train_size0.7, random_state 42)然后开始构造机器学习流。
# 构造机器学习流
pipe make_pipeline(PolynomialFeatures(), StandardScaler(), LogisticRegression(max_iterint(1e6)))# 查看参数
pipe.get_params()# 构造参数空间
param_grid [{polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [l1], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__solver: [saga]}, {polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [l2], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__solver: [lbfgs, newton-cg, sag, saga]},{polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [elasticnet], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__l1_ratio: np.arange(0.1, 1, 0.1).tolist(), logisticregression__solver: [saga]}
]三、优化评估指标选取
1. 高级评估指标的选用方法
根据此前介绍如果需要更好的验证模型本身泛化能力建议使用 F1-Score 或者 ROC-AUC当然调整网格搜索过程的模型评估指标过程其实并不难理解核心就是修改 scoring 参数取值。但由于涉及到在参数中调用评估函数因此需要补充一些关于常用分类评估指标在 sklearn 中的函数使用方法以及不同评估指标函数在不同参数取值时在网格搜索评估器中的调用方法。
GridSearchCV?从评估器的说明文档中能够看出scoring 参数最基础的情况下可以选择输入 str字符串或者 callable可调用对象也就是可以输入指代某个评估过程的字符串一个字符串代表不同参数取值下的某评估函数或者直接输入某评估指标函数或者通过 make_score 函数创建的函数来进行模型结果的评估。当然也可以在该参数位上直接输入一个字典或者 list其中如果是字典的话字典的 value 需要是 str字符串或者 callable可调用对象。由于 sklearn 中的评估指标函数一般都是有多个不同参数而不同参数代表不同的计算过程因此这些评估指标函数作为参数输入网格搜索评估器中的时候必须通过“某种方式”确定这些参数取值.因此就有了如下方法即通过字符串对应表来查看不同字符串所代表的不同参数取值下的评估指标函数如下所示 不难看出在网格搜索中输出评估指标参数和调用评估指标函数进行数据处理还是有很大的区别。例如metrics.roc_auc_score 函数能够同时处理多分类问题和二分类问题但如果作为参数输入到网格搜索中roc_auc 参数只能指代 metrics.roc_auc_score 函数的二分类功能.如果需要进行多分类则需要在 scoring 参数中输入 roc_auc_ovr、roc_auc_ovo 或者 roc_auc_ovr_weighted、roc_auc_ovo_weighted。我们先简单尝试在 scoring 中输入字符串的基本操作然后在深入解释 ROC-AUC 评估指标的使用方法。同时该参数列表也可以通过如下方式获得
import sklearn
sorted(sklearn.metrics.SCORERS.keys())from sklearn.metrics import roc_auc_scoreroc_auc_score?例如字符串 roc_auc_ovr 就代表 roc_auc_score 函数中 multi_class 参数取值为 ovr 时的计算流程。也就是说当网格搜索的 scoring 参数取值为字符串 roc_auc_ovr 时就代表调用了 multi_classovr、而其他参数选用默认参数的 roc_auc_score 函数作为模型评估函数。
GridSearchCV(estimatorpipe,param_gridparam_grid,scoringroc_auc_ovr)当然scoring 参数还支持直接输入可调用对象即支持输入经过 make_scorer 函数转化之后的评估指标函数
from sklearn.metrics import make_scoreracc make_scorer(roc_auc_score)GridSearchCV(estimatorpipe,param_gridparam_grid,scoringacc)但此时我们无法修改评估指标函数的默认参数。值得注意的是此处 make_scorer 函数实际上会将一个简单的评估指标函数转化为评估器结果评估函数。对于评估指标函数来说只需要输入标签的预测值和真实值即可进行计算例如
accuracy_score([1, 1, 0], [1, 1, 1])
#0.6666666666666666而评估器结果评估函数则需要同时输入评估器、特征矩阵以及对应的真实标签其执行过程是先将特征矩阵输入评估器、然后将输出结果和真实标签进行对比
acc make_scorer(accuracy_score)acc(search.best_estimator_, X_train, y_train)
#0.9732142857142857search.score(X_train, y_train)
#0.9732142857142857而在网格搜索或者交叉验证评估器中只支持输入经过 make_scorer 转化后的评估指标函数。
2. 同时输入多组评估指标
当然有的时候我们可能需要同时看不同参数下多项评估指标的结果此时我们就可以在 scoring 中输入列表、元组或者字典当然字典对象会较为常用。例如如果我们需要同时选用 roc-auc 和 accuracy 作为模型评估指标则需要创建如下字典
scoring {AUC: roc_auc, Accuracy: make_scorer(accuracy_score)}然后将其作为参数传入网格搜索评估器内
GridSearchCV?search GridSearchCV(estimatorclf,param_gridparam_grid_simple, scoringscoring, refitAccuracy)当然roc-auc 指标也可以用 make_score 来传输accuracy 也可以用字符串形式来传输即来上述多评估指标的字典等价于
{AUC: make_scorer(roc_auc_score), Accuracy: accuracy}
#{AUC: make_scorer(roc_auc_score), Accuracy: accuracy}不过需要注意的是尽管此时网格搜索评估器将同时计算一组参数下的多个评估指标结果并输出但我们只能选取其中一个评估指标作为挑选超参数的依据而其他指标尽管仍然会计算但结果只作参考。而 refit 参数中输入的评估指标就是最终选择参数的评估指标尽管网格搜索支持依据不同的评估指标进行参数搜索但最终选择何种参数可以参考如下依据有明确模型评估指标的在很多竞赛或者项目算法验收环节可能都会存在较为明确的模型评估指标例如模型排名根据 f1-score 计算结果得出等。在该情况下应当尽量选择要求的评估指标。没有明确模型评估指标的但是如果没有明确的评估指标要求则选择评估指标最核心的依据就是尽可能提升/确保模型的泛化能力。此时根据 Lesson 5 中对各评估指标的讨论结果如果数据集的各类别并没有明确的差异在算力允许的情况下应当优先考虑 roc-auc。而如果希望重点提升模型对类别 1或者某类别的识别能力则可以优先考虑 f1-score 作为模型评估指标。
四、优化后建模流程
接下来依据上述优化后的过程来执行网格搜索。完整流程如下1 构造机器学习流
#构造机器学习流
pipe make_pipeline(PolynomialFeatures(), StandardScaler(), LogisticRegression(max_iterint(1e6)))2 构造包含多个评估器的机器学习流的参数空间
# 构造参数空间
param_grid [{polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [l1], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__solver: [saga]}, {polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [l2], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__solver: [lbfgs, newton-cg, sag, saga]},{polynomialfeatures__degree: np.arange(2, 10).tolist(), logisticregression__penalty: [elasticnet], logisticregression__C: np.arange(0.1, 2, 0.1).tolist(), logisticregression__l1_ratio: np.arange(0.1, 1, 0.1).tolist(), logisticregression__solver: [saga]}
]3 实例化网格搜索评估器考虑到实际参数空间较大网格搜索需要耗费较长时间此处使用单一指标 roc 作为参数选取指标进行搜索
search GridSearchCV(estimatorpipe,param_gridparam_grid, scoringroc_auc, n_jobs5)4 执行训练
search.fit(X_train, y_train)5 查看结果
search.best_score_
#0.7879905483853072search.best_params_
#{logisticregression__C: 0.2,
# logisticregression__penalty: l1,
# logisticregression__solver: saga,
# polynomialfeatures__degree: 3}需要注意的是上述 best_score_ 属性查看的结果是在 roc-auc 评估指标下默认五折交叉验证时验证集上的 roc-auc 的平均值。但如果我们对训练好的评估器使用 .socre 方法查看的仍然是 pipe 评估器默认的结果评估方式也就是准确率计算结果
search.best_estimator_.score(X_train,y_train)
#0.7857142857142857search.best_estimator_.score(X_test,y_test)
#0.7866666666666666验证准确率计算结果
accuracy_score(search.best_estimator_.predict(X_train), y_train)
#0.7857142857142857accuracy_score(search.best_estimator_.predict(X_test), y_test)
#0.78666666666666666 结果分析最终模型结果准确率在 78% 上下。当然如果只看模型准确率结果我们发现该结果相比 Lesson 6.4 中结果较差Lesson 6.4 中测试集最高得分达到 0.8。但是该模型相比 Lesson 6.4 中模型来看该模型基本没有过拟合隐患测试集分数甚至高于训练集因此该模型在未来的使用过程中更有可能能够确保一个稳定的预测输出结果泛化能力更强。这也是交叉验证和 roc-auc 共同作用的结果。当然如果有明确要求根据准确率判断模型效果则上述过程应该选择准确率同时如果算力允许也可以近一步扩大搜索空间Lesson 6.4 中 0.8 的准确率就是在 15 阶多项式特征衍生基础上进行的运算。至此我们就完成了在实验数据上的建模调优。
