DataWhale第23期组队集成学习-Task8

一、bagging的思路

​ 与投票法不同的是，Bagging不仅仅集成模型最后的预测结果，同时采用一定策略来影响基模型训练，保证基模型可以服从一定的假设。在上一章中我们提到，希望各个模型之间具有较大的差异性，而在实际操作中的模型却往往是同质的，因此一个简单的思路是通过不同的采样增加模型的差异性。

二、bagging的原理分析

​ Bagging的核心在于自助采样(bootstrap)这一概念，即有放回的从数据集中进行采样，也就是说，同样的一个样本可能被多次进行采样。

一个自助采样的小例子:

​ 我们希望估计全国所有人口年龄的平均值，那么我们可以在全国所有人口中随机抽取不同的集合（这些集合可能存在交集），计算每个集合的平均值，然后将所有平均值的均值作为估计值。

​ 首先我们随机取出一个样本放入采样集合中，再把这个样本放回初始数据集，重复K次采样，最终我们可以获得一个大小为K的样本集合。同样的方法， 我们可以采样出T个含K个样本的采样集合，然后基于每个采样集合训练出一个基学习器，再将这些基学习器进行结合，这就是Bagging的基本流程。

​ 对回归问题的预测是通过预测取平均值来进行的。对于分类问题的预测是通过对预测取多数票预测来进行的。Bagging方法之所以有效，是因为每个模型都是在略微不同的训练数据集上拟合完成的，这又使得每个基模型之间存在略微的差异，使每个基模型拥有略微不同的训练能力。

​ Bagging同样是一种降低方差的技术，因此它在不剪枝决策树、神经网络等易受样本扰动的学习器上效果更加明显。在实际的使用中，加入列采样的Bagging技术对高维小样本往往有神奇的效果。

三、bagging的案例分析(基于sklearn，介绍随机森林的相关理论以及实例)

Sklearn为我们提供了 BaggingRegressor 与 BaggingClassifier 两种Bagging方法的API，我们在这里通过一个完整的例子演示Bagging在分类问题上的具体应用。这里两种方法的默认基模型是树模型。

3.1 树模型

​ 这里的树模型一般指决策树，它是一种树形结构，树的每个非叶子节点表示对样本在一个特征上的判断，节点下方的分支代表对样本的划分。决策树的建立过程是一个对数据不断划分的过程。每次划分中，首先要选择用于划分的特征，之后要确定划分的方案（类别/阈值）。我们希望通过划分，决策树的分支节点所包含的样本“纯度”尽可能地高。节点划分过程中所用的指标主要是信息增益和GINI系数。

信息增益衡量的是划分前后信息不确定性程度的减小。信息不确定程度一般使用信息熵来度量，其计算方式是：
$$H(Y)=-\sum p_i\log p_i$$
​ 其中i表示样本的标签，p表示该类样本出现的概率。当我们对样本做出划分之后，计算样本的条件熵：
$$H(Y|X)=-\sum _{x\in X}p(X=x)H(Y|X=x)$$
​ 其中x表示用于划分的特征的取值。信息增益定义为信息熵与条件熵的差值：
$$IG=H(Y)-H(Y|X)$$
​ 信息增益IG越大，说明使用该特征划分数据所获得的信息量变化越大，子节点的样本“纯度”越高。

​ 同样的，我们也可以利用Gini指数来衡量数据的不纯度，计算方法如下：
$$Gini=1-\sum p_i^2$$
​ 当我们对样本做出划分后，计算划分后的Gini指数：
$$Gin{i}_x=\sum _{x\in X}p(X=x)[1-\sum p_i^2]$$
​ 一般来说，我们选择使得划分后Gini指数最小的特征（注意这里是直接根据Gini指数进行判断，而 并非其变化量）。

给出了决策树的一个例子：现在需要训练一个模型，根据天气、温度和风力等级来判断是否打网球。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-az42Cj6H-1618606936330)(C:\Users\83769\Desktop\集成学习\中\集成学习（中）教程分任务版本\分任务版本\决策树例子.png)]

1. 通过计算信息增益或者Gini指数确定了先要根据天气情况对样本进行划分；
2. 对于每个分支，再继续考虑除天气之外的其他特征；
3. 直到样本的类别被完全分开，所有特征都已使用，或达到树的最大深度为止。

​ Bagging的一个典型应用是随机森林。顾名思义，“森林”是由许多“树”bagging组成的。在具体实现上， 用于每个决策树训练的样本和构建决策树的特征都是通过随机采样得到的，随机森林的预测结果是多个 决策树输出的组合（投票）。随机森林的示意图如下：

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-XEPkrWhm-1618606936333)(C:\Users\83769\Desktop\集成学习\中\集成学习（中）教程分任务版本\分任务版本\随机森林示意图.png)]

​ Bagging方法有很多种，其主要区别在于随机抽取训练子集的方法不同：

• 若抽取的数据集的随机子集是样例的随机子集，我们就叫做粘贴；
• 若样例抽取是有放回的，我们叫做Bagging；
• 若抽取的数据集的随机子集是特征的随机子集，我们叫做随机子空间；
• 若基估计器构建在对于样本和特征抽取的子集之上时，我们叫做随机补丁。

参考了链接：https://sklearn.apachecn.org/docs/master/12.html

3.2 BaggingClassifier

class sklearn.ensemble.BaggingClassifier(base_estimator=None, n_estimators=10, *, max_samples=1.0, max_features=1.0, bootstrap=True, bootstrap_features=False, oob_score=False, warm_start=False, n_jobs=None, random_state=None, verbose=0)

相关参数：

• base_estimator：一个对象，默认是None；
• n_estimators：数值型，默认值为10。总体的基数估计数；
• max_samples和max_features：控制子集的大小(对于样例和特征)；
• bootstrap和bootstrap_features：控制样例和特征的抽取是有放回还是无放回的。
• 当使用样本子集时，通过设置oob_score=True，可以使用袋外(out-of-bag)样本来评估泛化精度。

简单示例：

from sklearn.svm import SVC
from sklearn.ensemble import BaggingClassifier
from sklearn.datasets import make_classification
X, y = make_classification(n_samples=100, n_features=4,
                        n_informative=2, n_redundant=0,
                        random_state=0, shuffle=False)
clf = BaggingClassifier(base_estimator=SVC(),
                         n_estimators=10, random_state=0).fit(X, y)
clf.predict([[0, 0, 0, 0]])

array([1])

3.3 BaggingRegressior

class sklearn.ensemble.BaggingClassifier(base_estimator=None, n_estimators=10, *, max_samples=1.0, max_features=1.0, bootstrap=True, bootstrap_features=False, oob_score=False, warm_start=False, n_jobs=None, random_state=None, verbose=0)

BaggingRegressior跟BaggingClassifier的参数及其使用方法是类似的，因此就不重复赘述。

简单示例：

from sklearn.svm import SVR
from sklearn.ensemble import BaggingRegressor
from sklearn.datasets import make_regression
X, y = make_regression(n_samples=100, n_features=4,
                        n_informative=2, n_targets=1,
                        random_state=0, shuffle=False)
regr = BaggingRegressor(base_estimator=SVR(),
                         n_estimators=10, random_state=0).fit(X, y)
regr.predict([[0, 0, 0, 0]])

array([-2.84440667])

4. 例子

我们创建一个含有1000个样本20维特征的随机分类数据集：

# test classification dataset
from sklearn.datasets import make_classification
# define dataset
X, y = make_classification(n_samples=1000, n_features=20, n_informative=15, n_redundant=5, random_state=5)
# summarize the dataset
print(X.shape, y.shape)

(1000, 20) (1000,)

我们将使用重复的分层k-fold交叉验证来评估该模型，一共重复3次，每次有10个fold。我们将评估该模型在所有重复交叉验证中性能的平均值和标准差。

# evaluate bagging algorithm for classification
from numpy import mean
from numpy import std
from sklearn.datasets import make_classification
from sklearn.model_selection import cross_val_score
from sklearn.model_selection import RepeatedStratifiedKFold
from sklearn.ensemble import BaggingClassifier
# define dataset
X, y = make_classification(n_samples=1000, n_features=20, n_informative=15, n_redundant=5, random_state=5)
# define the model
model = BaggingClassifier()
# evaluate the model
cv = RepeatedStratifiedKFold(n_splits=10, n_repeats=3, random_state=1)
n_scores = cross_val_score(model, X, y, scoring='accuracy', cv=cv, n_jobs=-1, error_score='raise')
# report performance
print('Accuracy: %.3f (%.3f)' % (mean(n_scores), std(n_scores)))

Accuracy: 0.866(0.031)