算法工程师是伴随着人工智能火起来的一个领域。听着名字似乎门槛很高。但是,得益于Python生态下的包共享机制,机器模型构建的过程其实已经变得非常简单了,很多听起来牛逼的算法,其实根本不需要自己实现,甚至都不需要知道这些算法的具体原理。
你只需要两步就能构建起自己的机器学习模型:
- 明确你需要解决的问题是什么类型,以及知道解决该类型问题所对应的算法。
- 从skicit-learn中调用相应的算法构建模型即可。是的!在机器学习领域,如果你只是抱着体验机器学习的心态,实现起来就是这么简单。
第一步很好解决
常见的问题类型只有三种:分类、回归、聚类。而明确具体问题对应的类型也很简单。比如,如果你需要通过输入数据得到一个类别变量,那就是分类问题。分成两类就是二分类问题,分成两类以上就是多分类问题。常见的有:判别一个邮件是否是垃圾邮件、根据图片分辩图片里的是猫还是狗等等。
如果你需要通过输入数据得到一个具体的连续数值,那就是回归问题。比如:预测某个区域的房价等。
常用的分类和回归算法算法有:SVM (支持向量机) 、xgboost、, KNN、LR算法、SGD (随机梯度下降算法)、Bayes (贝叶斯估计)以及随机森林等。这些算法大多都既可以解分类问题,又可以解回归问题。
如果你的数据集并没有对应的属性标签,你要做的,是发掘这组样本在空间的分布, 比如分析哪些样本靠的更近,哪些样本之间离得很远, 这就是属于聚类问题。常用的聚类算法有k-means算法。
在本文中,我们主要解决第二步:通过skicit-learn构建模型。告诉你你一套让你简单到想笑的通用模型构建模板。只要scikit-learn实现的算法,都可以通过这种方式快速调用。牢记这个万能模板,你就能轻松构建起自己的机器学习模型。
预备工作
在介绍万能模板之前,为了能够更深刻地理解这三个模板,我们加载一个Iris(鸢尾花)数据集来作为应用万能模板的小例子,Iris数据集在前边的文章中已经提到过多次了,这里不再赘述。它是一个典型的多分类问题。加载步骤如下:
1、加载数据集
因为原始的数据集中包含很多空值,而且类别特征用英文名表示各个花的名字,也需要我们转换成数字。
在scikit-learn下的datasets子包里,也自带了一个Iris数据集,这个数据集和原始数据集的区别就是scikit-learn已经帮我们提前处理好了空值等问题,可以直接输入模型用来训练。所以为了方便起见,我们直接使用scikit-learn的数据集。加载方法如下:
from sklearn.datasets import load_iris data = load_iris() x = data.data y = data.target
x值如下,可以看到scikit-learn把数据集经过去除空值处理放在了array里,所以x是一个(150,4)的数组,保存了150个数据的4个特征:
array([[5.1, 3.5, 1.4, 0.2], [4.9, 3. , 1.4, 0.2], [4.7, 3.2, 1.3, 0.2], [4.6, 3.1, 1.5, 0.2], [5. , 3.6, 1.4, 0.2], [5.4, 3.9, 1.7, 0.4], [4.6, 3.4, 1.4, 0.3], [5. , 3.4, 1.5, 0.2], [4.4, 2.9, 1.4, 0.2], [4.9, 3.1, 1.5, 0.1], [5.4, 3.7, 1.5, 0.2], [4.8, 3.4, 1.6, 0.2], [4.8, 3. , 1.4, 0.1], [4.3, 3. , 1.1, 0.1], …………
y值如下,共有150行,其中0、1、2分别代表三类花:
array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2])
2、数据集拆分
数据集拆分是为了验证模型在训练集和测试集是否过拟合,使用train_test_split
的目的是保证从数据集中均匀拆分出测试集。这里,简单把10%的数据集拿出来用作测试集。
from sklearn.model_selection import train_test_split train_x,test_x,train_y,test_y = train_test_split(x,y,test_size=0.1,random_state=0)
万能模板
助你快速构建一个基本的算法模型
不同的算法只是改变了名字,以及模型的参数不同而已。
有了这个万能模板,接下来就是简单的复制粘贴改名字了:
而且在scikit-learn中,每个包的位置都是有规律的,比如:随机森林就是在集成学习文件夹下。
模板应用案例:
1、构建SVM分类模型
通过查阅资料,我们知道svm算法在scikit-learn.svm.SVC
下,所以:
- 算法位置填入:
svm
- 算法名填入:
SVC()
- 模型名自己起,这里我们就叫
svm_model
套用模板得到程序如下:
# svm分类器 from sklearn.svm import SVC from sklearn.metrics import accuracy_score svm_model = SVC() svm_model.fit(train_x,train_y) pred1 = svm_model.predict(train_x) accuracy1 = accuracy_score(train_y,pred1) print('在训练集上的精确度: %.4f'%accuracy1) pred2 = svm_model.predict(test_x) accuracy2 = accuracy_score(test_y,pred2) print('在测试集上的精确度: %.4f'%accuracy2)
输出:
在训练集上的精确度: 0.9810 在测试集上的精确度: 0.9778
2、构建LR分类模型
同理,找到LR算法在sklearn.linear_model.LogisticRegression
下,所以:
- 算法位置填入:
linear_model
- 算法名填入:
LogisticRegression
- 模型名叫做:lr_model。
程序如下:
套用模板得到程序如下:
# LogisticRegression分类器 from sklearn.linear_model import LogisticRegression from sklearn.metrics import accuracy_score #评分函数用精确度评估 lr_model = LogisticRegression() lr_model.fit(train_x,train_y) pred1 = lr_model.predict(train_x) accuracy1 = accuracy_score(train_y,pred1) print('在训练集上的精确度: %.4f'%accuracy1) pred2 = lr_model.predict(test_x) accuracy2 = accuracy_score(test_y,pred2) print('在测试集上的精确度: %.4f'%accuracy2)
输出:
在训练集上的精确度: 0.9429
在测试集上的精确度: 0.8889
3、构建随机森林分类模型
随机森林算法在sklearn.ensemble.RandomForestClassifier
下,好了,现在你应该可以自己写了,这个作为本文的一个小测试,欢迎在评论区写下你的答案。
万能模板改进版
加入交叉验证,让算法模型评估更加科学
在1.0版的模板中,当你多次运行同一个程序就会发现:每次运行得到的精确度并不相同,而是在一定范围内浮动,这是因为数据输入模型之前会进行选择,每次训练时数据输入模型的顺序都不一样。所以即使是同一个程序,模型最后的表现也会有好有坏。
更糟糕的是,有些情况下,在训练集上,通过调整参数设置使模型的性能达到了最佳状态,但在测试集上却可能出现过拟合的情况。这个时候,我们在训练集上得到的评分不能有效反映出模型的泛化性能。
为了解决上述两个问题,还应该在训练集上划分出验证集(validation set)并结合交叉验证来解决。首先,在训练集中划分出不参与训练的验证集,只是在模型训练完成以后对模型进行评估,接着再在测试集上进行最后的评估。
但这样大大减少了可用于模型学习的样本数量,所以还需要采用交叉验证的方式多训练几次。比如说最常用的k-折交叉验证如下图所示,它主要是将训练集划分为 k 个较小的集合。然后将k-1份训练子集作为训练集训练模型,将剩余的 1 份训练集子集作为验证集用于模型验证。这样需要训练k次,最后在训练集上的评估得分取所有训练结果评估得分的平均值。
这样一方面可以让训练集的所有数据都参与训练,另一方面也通过多次计算得到了一个比较有代表性的得分。唯一的缺点就是计算代价很高,增加了k倍的计算量。
原理就是这样,但理想很丰满,现实很骨干。在自己实现的时候却有一个很大的难题摆在面前:怎么能够把训练集均匀地划分为K份?
这个问题不用思考太多,既然别忘了,我们现在是站在巨人的肩膀上,scikit-learn已经将优秀的数学家所想到的均匀拆分方法和程序员的智慧融合在了cross_val_score() 这个函数里了,只需要调用该函数即可,不需要自己想什么拆分算法,也不用写for循环进行循环训练。
万能模板改进如下:
把模型、数据、划分验证集的个数一股脑输入函数,函数会自动执行上边所说的过程。
在求精确度的时候,我们可以简单地输出平均精确度:
# 输出精确度的平均值 print("训练集上的精确度: %0.2f " % scores1.mean())
改进模板应用案例:
1、构建SVM分类模型
程序如下:
### svm分类器 from sklearn.model_selection import cross_val_score from sklearn.svm import SVC svm_model = SVC() svm_model.fit(train_x,train_y) scores1 = cross_val_score(svm_model,train_x,train_y,cv=5, scoring='accuracy') # 输出精确度的平均值和置信度区间 print("训练集上的精确度: %0.2f (+/- %0.2f)" % (scores1.mean(), scores1.std() * 2)) scores2 = cross_val_score(svm_model,test_x,test_y,cv=5, scoring='accuracy') # 输出精确度的平均值和置信度区间 print("测试集上的平均精确度: %0.2f (+/- %0.2f)" % (scores2.mean(), scores2.std() * 2)) print(scores1) print(scores2)
输出:
输出:
训练集上的精确度: 0.97 (+/- 0.08)
测试集上的平均精确度: 0.91 (+/- 0.10)
[1. 1. 1. 0.9047619 0.94736842]
[1. 0.88888889 0.88888889 0.875 0.875 ]
2、构建LR分类模型
代码如下
# LogisticRegression分类器 from sklearn.model_selection import cross_val_score from sklearn.linear_model import LogisticRegression lr_model = LogisticRegression() lr_model.fit(train_x,train_y) scores1 = cross_val_score(lr_model,train_x,train_y,cv=5, scoring='accuracy') # 输出精确度的平均值和置信度区间 print("训练集上的精确度: %0.2f (+/- %0.2f)" % (scores1.mean(), scores1.std() * 2)) scores2 = cross_val_score(lr_model,test_x,test_y,cv=5, scoring='accuracy') # 输出精确度的平均值和置信度区间 print("测试集上的平均精确度: %0.2f (+/- %0.2f)" % (scores2.mean(), scores2.std() * 2)) print(scores1) print(scores2)
输出:
训练集上的精确度: 0.94 (+/- 0.07) 测试集上的平均精确度: 0.84 (+/- 0.14) [0.90909091 1. 0.95238095 0.9047619 0.94736842] [0.90909091 0.88888889 0.88888889 0.75 0.75 ]
随机森林依旧留作小测试。
注: 如果想要一次性评估多个指标,也可以使用可以一次性输入多个评估指标的 cross_validate()
函数。