sklearn（2）—— 决策树

1 概述

1.1 决策树是如何工作的

决策树（Decision Tree）是一种非参数的有监督学习方法，它能够从一系列有特征和标签的数据中总结出决策规则，并用树状图的结构来呈现这些规则，以解决分类和回归问题。决策树算法容易理解，适用各种数据，在解决各种问题时都有良好表现，尤其是以树模型为核心的各种集成算法，在各个行业和领域都有广泛的应用。

我们来简单了解一下决策树是如何工作的。决策树算法的本质是一种图结构，我们只需要问一系列问题就可以对数据进行分类了。比如说，来看看下面这组数据集，这是一系列已知物种以及所属类别的数据：

我们现在的目标是，将动物们分为哺乳类和非哺乳类。那根据已经收集到的数据，决策树算法为我们算出了下面的这棵决策树：

假如我们现在发现了一种新物种Python，它是冷血动物，体表带鳞片，并且不是胎生，我们就可以通过这棵决策树来判断它的所属类别。
可以看出，在这个决策过程中，我们一直在对记录的特征进行提问。最初的问题所在的地方叫做根节点，在得到结论前的每一个问题都是中间节点，而得到的每一个结论（动物的类别）都叫做叶子节点。

关键概念：节点

• 根节点：没有进边，有出边。包含最初的，针对特征的提问。
• 中间节点：既有进边也有出边，进边只有一条，出边可以有很多条。都是针对特征的提问。
• 叶子节点：有进边，没有出边，每个叶子节点都是一个类别标签。
• 子节点和父节点：在两个相连的节点中，更接近根节点的是父节点，另一个是子节点。

决策树算法的核心是要解决两个问题：

1. 如何从数据表中找出最佳节点和最佳分枝？
2. 如何让决策树停止生长，防止过拟合？

几乎所有决策树有关的模型调整方法，都围绕这两个问题展开。这两个问题背后的原理十分复杂，我们会在讲解模型参数和属性的时候为大家简单解释涉及到的部分。在这门课中，我会尽量避免让大家太过深入到决策树复杂的原理和数学公式中（尽管决策树的原理相比其他高级的算法来说是非常简单了），这门课会专注于实践和应用。如果大家希望理解更深入的细节，建议大家在听这门课之前还是先去阅读和学习一下决策树的原理。

1.2 sklearn中的决策树

模块sklearn.tree
sklearn中决策树的类都在”tree“这个模块之下。这个模块总共包含五个类：

tree.DecisionTreeClassifier	分类树
tree.DecisionTreeRegressor	回归树
tree.export_graphviz	将生成的决策树导出为DOT格式，画图专用
tree.ExtraTreeClassifier	高随机版本的分类树
tree.ExtraTreeRegressor	高随机版本的回归树

我们会主要讲解分类树和回归树，并用图像呈现给大家。

• sklearn的基本建模流程

在那之前，我们先来了解一下sklearn建模的基本流程。

在这个流程下，分类树对应的代码是：

from sklearn import tree            # 导入需要的模块
clf = tree.DecisionTreeClassifier() # 实例化
clf = clf.fit()                     # 用训练集数据训练模型
result = clf.score()                # 导入测试集，从接口中调用需要的信息

回归树中二维数组，分类树中一维数组

2 DecisionTreeClassifier（分类树与红酒数据集）

class sklearn.tree.DecisionTreeClassifier (criterion=’gini’, splitter=’best’, max_depth=None,
min_samples_split=2, min_samples_leaf=1, min_weight_fraction_leaf=0.0, max_features=None,
random_state=None, max_leaf_nodes=None, min_impurity_decrease=0.0, min_impurity_split=None,class_weight=None, presort=False)

2.1 重要参数

2.1.1 criterion

为了要将表格转化为一棵树，决策树需要找出最佳节点和最佳的分枝方法，对分类树来说，衡量这个“最佳”的指标叫做“不纯度”。通常来说，不纯度越低，决策树对训练集的拟合越好。现在使用的决策树算法在分枝方法上的核心大多是围绕在对某个不纯度相关指标的最优化上。

不纯度基于节点来计算，树中的每个节点都会有一个不纯度，并且子节点的不纯度一定是低于父节点的，也就是说，在同一棵决策树上，叶子节点的不纯度一定是最低的。

Criterion这个参数正是用来决定不纯度的计算方法的。sklearn提供了两种选择：
1）输入”entropy“，使用信息熵（Entropy）
2）输入”gini“，使用基尼系数（Gini Impurity）

其中t代表给定的节点，i代表标签的任意分类，$p(i|t)$代表标签分类i在节点t上所占的比例。注意，当使用信息熵时，sklearn实际计算的是基于信息熵的信息增益(Information Gain)，即父节点的信息熵和子节点的信息熵之差。比起基尼系数，信息熵对不纯度更加敏感，对不纯度的惩罚最强。但是在实际使用中，信息熵和基尼系数的效果基本相同。信息熵的计算比基尼系数缓慢一些，因为基尼系数的计算不涉及对数。另外，因为信息熵对不纯度更加敏感，所以信息熵作为指标时，决策树的生长会更加“精细”，因此对于高维数据或者噪音很多的数据，信息熵很容易过拟合，基尼系数在这种情况下效果往往比较好。当然，这不是绝对的。

参数	criterion
如何影响模型?	确定不纯度的计算方法，帮忙找出最佳节点和最佳分枝，不纯度越低，决策树对训练集的拟合越好
可能的输入有哪些？	不填默认基尼系数，填写gini使用基尼系数，填写entropy使用信息增益
怎样选取参数？	通常就使用基尼系数数据维度很大，噪音很大时使用基尼系数 维度低，数据比较清晰的时候，信息熵和基尼系数没区别， 当决策树的拟合程度不够的时候，使用信息熵 两个都试试，不好就换另外一个

到这里，决策树的基本流程其实可以简单概括如下：

直到没有更多的特征可用，或整体的不纯度指标已经最优，决策树就会停止生长。

• 建立一棵树

1. 导入需要的算法库和模块

# 1.导入需要的算法库和模块
from sklearn import tree
from sklearn.datasets import load_wine
from sklearn.model_selection import train_test_split

2. 探索数据

# 2.探索数据
wine = load_wine()
wine.data.shape # (178, 13)

wine.target
'''
输出
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, 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, 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])
'''

如果wine是一张表，应该是这样(结果如下图所示)

# 如果wine是一张表，应该是这样
import pandas as pd
pd.concat([pd.DataFrame(wine.data),pd.DataFrame(wine.target)],axis=1)

wine.feature_names
'''
输出
['alcohol',
 'malic_acid',
 'ash',
 'alcalinity_of_ash',
 'magnesium',
 'total_phenols',
 'flavanoids',
 'nonflavanoid_phenols',
 'proanthocyanins',
 'color_intensity',
 'hue',
 'od280/od315_of_diluted_wines',
 'proline']
'''

wine.target_names
'''
输出
array(['class_0', 'class_1', 'class_2'], dtype='

3. 分训练集和测试集

# 3.分训练集和测试集
# test_size=0.3 表示30%做测试集，70%做训练集
Xtrain,Xtest,Ytrain,Ytest = train_test_split(wine.data,wine.target,test_size=0.3)

Xtrain.shape
'''
输出
(124, 13)
'''

Xtest.shape
'''
输出
(54, 13)
'''

4. 建立模型

# 4. 建立模型
clf = tree.DecisionTreeClassifier(criterion="entropy")
clf = clf.fit(Xtrain,Ytrain)
score = clf.score(Xtest,Ytest) # 返回预测的准确度

score
'''
输出
0.8333333333333334
'''

5. 画出一棵树

# 5.画出一棵树
feature_name = ['酒精','苹果酸','灰','灰的碱性','镁','总酚','类黄酮','非黄烷类酚类','花青素','颜色强度','色调','od280/od315稀释葡萄酒','脯氨酸']

import graphviz
dot_data = tree.export_graphviz(clf,
                                feature_names=feature_name,
                                class_names=["琴酒","雪莉","贝尔摩德"],
                                filled=True,
                                rounded=True
                                )
graph = graphviz.Source(dot_data)
graph
'''
输出
如下图所示
'''

6. 探索决策树

# 6.探索决策树

# 特征重要性
clf.feature_importances_
'''
输出
array([0.        , 0.03472593, 0.02222161, 0.        , 0.        ,
       0.        , 0.44246035, 0.        , 0.        , 0.        ,
       0.        , 0.08561031, 0.41498181])
'''

[*zip(feature_name,clf.feature_importances_)]
'''
输出
[*zip(feature_name,clf.feature_importances_)]

[('酒精', 0.0),
 ('苹果酸', 0.03472592550375031),
 ('灰', 0.022221605749800283),
 ('灰的碱性', 0.0),
 ('镁', 0.0),
 ('总酚', 0.0),
 ('类黄酮', 0.44246034690442393),
 ('非黄烷类酚类', 0.0),
 ('花青素', 0.0),
 ('颜色强度', 0.0),
 ('色调', 0.0),
 ('od280/od315稀释葡萄酒', 0.08561031480649588),
 ('脯氨酸', 0.41498180703552967)]
'''

我们已经在只了解一个参数的情况下，建立了一棵完整的决策树。但是回到步骤4建立模型，score会在某个值附近波动，引起步骤5中画出来的每一棵树都不一样。它为什么会不稳定呢？如果使用其他数据集，它还会不稳定吗？我们之前提到过，无论决策树模型如何进化，在分枝上的本质都还是追求某个不纯度相关的指标的优化，而正如我们提到的，不纯度是基于节点来计算的，也就是说，决策树在建树时，是靠优化节点来追求一棵优化的树，但最优的节点能够保证最优的树吗？集成算法被用来解决这个问题：sklearn表示，既然一棵树不能保证最优，那就建更多的不同的树，然后从中取最好的。怎样从一组数据集中建不同的树？在每次分枝时，不从使用全部特征，而是随机选取一部分特征，从中选取不纯度相关指标最优的作为分枝用的节点。这样，每次生成的树也就不同了。

# random_state代表随机模式的参数
clf = tree.DecisionTreeClassifier(criterion="entropy",random_state=30)
clf = clf.fit(Xtrain,Ytrain)
score = clf.score(Xtest,Ytest) # 返回预测的准确度

score
'''
输出
0.9444444444444444
'''

2.1.2 random_state & splitter

random_state用来设置分枝中的随机模式的参数，默认None，在高维度时随机性会表现更明显，低维度的数据（比如鸢尾花数据集），随机性几乎不会显现。输入任意整数，会一直长出同一棵树，让模型稳定下来。

splitter也是用来控制决策树中的随机选项的，有两种输入值，输入”best"，决策树在分枝时虽然随机，但是还是会优先选择更重要的特征进行分枝（重要性可以通过属性feature_importances_查看），输入“random"，决策树在分枝时会更加随机，树会因为含有更多的不必要信息而更深更大，并因这些不必要信息而降低对训练集的拟合。这也是防止过拟合的一种方式。当你预测到你的模型会过拟合，用这两个参数来帮助你降低树建成之后过拟合的可能性。当然，树一旦建成，我们依然是使用剪枝参数来防止过拟合。

clf = tree.DecisionTreeClassifier(criterion="entropy"
                                  ,random_state=30
                                  ,splitter="random"
                                 )
clf = clf.fit(Xtrain,Ytrain)
score = clf.score(Xtest,Ytest) # 返回预测的准确度
score
'''
输出：0.9814814814814815
'''

dot_data = tree.export_graphviz(clf
                                ,feature_names=feature_name
                                ,class_names=["琴酒","雪莉","贝尔摩德"]
                                ,filled=True
                                ,rounded=True
                                )
graph = graphviz.Source(dot_data)
graph
'''
输出：如下图所示
'''

2.1.3 剪枝参数

在不加限制的情况下，一棵决策树会生长到衡量不纯度的指标最优，或者没有更多的特征可用为止。这样的决策树往往会过拟合，这就是说，它会在训练集上表现很好，在测试集上却表现糟糕。我们收集的样本数据不可能和整体的状况完全一致，因此当一棵决策树对训练数据有了过于优秀的解释性，它找出的规则必然包含了训练样本中的噪声，并使它对未知数据的拟合程度不足。

# 我们的树对训练集的拟合程度如何：

score_train = clf.score(Xtrain,Ytrain)
score_train
'''
输出：1.0
'''

# 我们的树对测试集的拟合程度如何：

score_test = clf.score(Xtest,Ytest)
score_test
'''
输出：0.9814814814814815
'''

为了让决策树有更好的泛化性，我们要对决策树进行剪枝。剪枝策略对决策树的影响巨大，正确的剪枝策略是优化决策树算法的核心。sklearn为我们提供了不同的剪枝策略：

• max_depth
  限制树的最大深度，超过设定深度的树枝全部剪掉这是用得最广泛的剪枝参数，在高维度低样本量时非常有效。决策树多生长一层，对样本量的需求会增加一倍，所以限制树深度能够有效地限制过拟合。在集成算法中也非常实用。实际使用时，建议从=3开始尝试，看看拟合的效果再决定是否增加设定深度。
• min_samples_leaf & min_samples_split
  min_samples_leaf限定，一个节点在分枝后的每个子节点都必须包含至少min_samples_leaf个训练样本，否则分枝就不会发生，或者，分枝会朝着满足每个子节点都包含min_samples_leaf个样本的方向去发生。

一般搭配max_depth使用，在回归树中有神奇的效果，可以让模型变得更加平滑。这个参数的数量设置得太小会引起过拟合，设置得太大就会阻止模型学习数据。一般来说，建议从=5开始使用。如果叶节点中含有的样本量变化很大，建议输入浮点数作为样本量的百分比来使用。同时，这个参数可以保证每个叶子的最小尺寸，可以在回归问题中避免低方差，过拟合的叶子节点出现。对于类别不多的分类问题，=1通常就是最佳选择。

min_samples_split限定，一个节点必须要包含至少min_samples_split个训练样本，这个节点才允许被分枝，否则分枝就不会发生：

'''
max_depth=3
min_samples_leaf=10
min_impurity_split=5
这三个参数逐一验证，当测试集的拟合程度提升时留下该参数，降低时则不添加该参数
'''
clf = tree.DecisionTreeClassifier(criterion="entropy"
                                  ,random_state=30
                                  ,splitter="random"
                                # ,max_depth=3
                                # ,min_samples_leaf=10
                                # ,min_impurity_split=5
                                 )
clf = clf.fit(Xtrain,Ytrain)
dot_data = tree.export_graphviz(clf
                                ,feature_names=feature_name
                                ,class_names=["琴酒","雪莉","贝尔摩德"]
                                ,filled=True
                                ,rounded=True
                                )
graph = graphviz.Source(dot_data)
graph

• max_features & min_impurity_decrease

一般max_depth使用，用作树的”精修“

max_features限制分枝时考虑的特征个数，超过限制个数的特征都会被舍弃。和max_depth异曲同工，max_features是用来限制高维度数据的过拟合的剪枝参数，但其方法比较暴力，是直接限制可以使用的特征数量而强行使决策树停下的参数，在不知道决策树中的各个特征的重要性的情况下，强行设定这个参数可能会导致模型学习不足。如果希望通过降维的方式防止过拟合，建议使用PCA，ICA或者特征选择模块中的降维算法。

min_impurity_decrease限制信息增益的大小，信息增益小于设定数值的分枝不会发生。这是在0.19版本种更新的功能，在0.19版本之前时使用min_impurity_split。

信息增益：父节点和子节点之间 信息熵 的差。子节点的信息熵一定小于父节点的信息熵。当信息增益越大时，代表这一层的分支对决策树的贡献越大。所以可以限制当信息增益小于某个值是，不在发生分支。例<0.001。

• 确认最优的剪枝参数

那具体怎么来确定每个参数填写什么值呢？这时候，我们就要使用确定超参数的曲线来进行判断了，继续使用我们已经训练好的决策树模型clf。超参数的学习曲线，是一条以超参数的取值为横坐标，模型的度量指标为纵坐标的曲线，它是用来衡量不同超参数取值下模型的表现的线。在我们建好的决策树里，我们的模型度量指标就是score

import matplotlib.pyplot as plt
test = []
for i in range(10):
    clf = tree.DecisionTreeClassifier(max_depth=i+1
                                    ,criterion="entropy"
                                    ,random_state=30
                                    ,splitter="random"
                                    )
    clf = clf.fit(Xtrain, Ytrain)
    score = clf.score(Xtest, Ytest)
    test.append(score)
plt.plot(range(1,11),test,color="red",label="max_depth")
plt.legend()
plt.show()
'''
输出：如下图所示
'''

思考：

1. 剪枝参数一定能够提升模型在测试集上的表现吗？

• 调参没有绝对的答案，一切都是看数据本身。

2. 这么多参数，一个个画学习曲线？

无论如何，剪枝参数的默认值会让树无尽地生长，这些树在某些数据集上可能非常巨大，对内存的消耗。所以如果你手中的数据集非常大，你已经预测到无论如何你都是要剪枝的，那提前设定这些参数来控制树的复杂性和大小会比较好。

2.1.4 目标权重参数

• class_weight & min_weight_fraction_leaf

完成样本标签平衡的参数。样本不平衡是指在一组数据集中，标签的一类天生占有很大的比例。比如说，在银行要判断“一个办了信用卡的人是否会违约”，就是是vs否（1%：99%）的比例。这种分类状况下，即便模型什么也不做，全把结果预测成“否”，正确率也能有99%。因此我们要使用class_weight参数对样本标签进行一定的均衡，给少量的标签更多的权重，让模型更偏向少数类，向捕获少数类的方向建模。该参数默认None，此模式表示自动给与数据集中的所有标签相同的权重。

有了权重之后，样本量就不再是单纯地记录数目，而是受输入的权重影响了，因此这时候剪枝，就需要搭配min_weight_fraction_leaf这个基于权重的剪枝参数来使用。另请注意，基于权重的剪枝参数（例如min_weight_fraction_leaf）将比不知道样本权重的标准（比如min_samples_leaf）更少偏向主导类。如果样本是加权的，则使用基于权重的预修剪标准来更容易优化树结构，这确保叶节点至少包含样本权重的总和的一小部分。

2.2 重要属性和接口

属性是在模型训练之后，能够调用查看的模型的各种性质。对决策树来说，最重要的是feature_importances_，能够查看各个特征对模型的重要性。
sklearn中许多算法的接口都是相似的，比如说我们之前已经用到的fit和score，几乎对每个算法都可以使用。除了这两个接口之外，决策树最常用的接口还有apply和predict。apply中输入测试集返回每个测试样本所在的叶子节点的索引，predict输入测试集返回每个测试样本的标签。返回的内容一目了然并且非常容易，大家感兴趣可以自己下去试试看。

在这里不得不提的是。所有接口中要求输入X_train和X_test的部分，输入的特征矩阵必须至少是一个二维矩阵。sklearn不接受任何一维矩阵作为特征矩阵被输入。如果你的数据的确只有一个特征，那必须用 reshape（-1,1）来给矩阵增维；如果你的数据只有一个特征和一个样本，使用 reshape（1,-1）来给你的数据增维.

#apply返回每个测试样本所在的叶子节点的索引
clf.apply(Xtest)
'''
array([ 4, 12,  4,  4, 15, 13, 15, 19, 12, 19, 15, 15, 15,  4,  4, 19, 15,
        7, 15, 20, 15, 12, 19, 19, 12,  4, 15, 15, 15, 19, 12, 15, 13, 19,
       15, 19, 12,  9, 15, 19, 15, 10,  4,  7,  4, 15, 18, 19, 20,  4, 18,
       15, 19,  4], dtype=int64)
'''

#predict返回每个测试样本的分类/回归结果
clf.predict(Xtest)
'''
array([2, 1, 2, 2, 1, 2, 1, 0, 1, 0, 1, 1, 1, 2, 2, 0, 1, 2, 1, 1, 1, 1,
       0, 0, 1, 2, 1, 1, 1, 0, 1, 1, 2, 0, 1, 0, 1, 2, 1, 0, 1, 1, 2, 2,
       2, 1, 1, 0, 1, 2, 1, 1, 0, 2])
'''

至此，我们已经学完了分类树DecisionTreeClassifier和用决策树绘图（export_graphviz）的所有基础。我们讲解了决策树的基本流程，分类树的八个参数，一个属性，四个接口，以及绘图所用的代码。

八个参数：一个决定不纯度计算方法的属性（Criterion），两个随机性相关的参数（random_state，splitter），五个剪枝参数（max_depth, min_sample_split，min_sample_leaf，max_feature，min_impurity_decrease）

一个属性：feature_importances_
四个接口：fit，score，apply，predict

有了这些知识，基本上分类树的使用大家都能够掌握了，接下来再到实例中去磨练就好。

2.3 实例：分类树在合成数集上的表现

我们在红酒数据集上画出了一棵树，并且展示了多个参数会对树形成这样的影响，接下来，我们将在不同结构的数据集上测试一下决策树的效果，让大家更好地理解决策树。

3 Decision Tree Regressor（回归树）

几乎所有参数，属性及接口都和分类树一模一样。需要注意的是，在回归树种，没有标签分布是否均衡的问题，因此没有 class_ weigh这样的参数

3.1 重要参数，属性及接口

criterion

注：与分类树的criterion相比没有不纯度的计算

回归树衡量分枝质量的指标，支持的标准有三种:

1）输入"mse"使用均方误差 mean squared error（MSE），父节点和叶子节点之间的均方误差的差额将被做为特征选择的标准，这种方法通过使用叶子节点的均值来最小化L2损失.

2）输入" friedman mse"使用费尔德曼均方误差，这种指标使用弗里德曼针对潜在分枝中的问题改进后的均方误差.（0.20版本中含有）

3）输入"mae"使用绝对平均误差MAE（ mean absolute error），这种指标使用叶节点的中值来最小化L1损失属性中最重要的依然是 feature_importances_，接口依然是 apply, fit, predict, score最核心.

其中N是样本数量，禔是每一个数据样本，f是模型回归出的数值，y是样本点际的数值标签。所以MSE的本质其实是样本真实数据与回归结果的差异。在回归树中，MSE不只是我们的分枝质量衡量指标，也是我们最常用的衡量回归树回归质量的指标，当我们在使用交叉验证，或者其他方式获取回归树的结果时，我们往往选择均方误差作为我们的评估（在分类树中这个指标是 score代表的预测准确率）。在回归中，我们追求的是，MSE越小越好。

然而，回归树的接口 score返回的是R平方，并不是MSE.R平方被定义如下：

其中u是残差平方和（MSE*N），是总平方和，N是样本数量，禔是每一个数据样本，f是模型回归出的数值，yi是样本点实际的数值标签。y帽是真实数值标签的平均数。R平方可以为正为负（如果模型的残差平方和远远大于模型的总平方和，模型非常糟糕，R平方就会为负），而均方误差永远为正。

R平方越接近1越好，均方误差越小越好（即越接近0越好），则负的均方误差是越大越好（即越接近0越好）

值得一提的是，虽然均方误差永远为正，但是 sklearn当中使用均方误差作为评判标准时，却是计算"负均方误差”（ neg_mean_squared_error）。这是因为 sklearn在计算模型评估指标的时候，会考虑指标本身的性质，均方误差本身是一种误差，所以被 sklearn划分为模型的一种损失（oss，因此在 sklearn当中，都以负数表示。真正的均方误差MSE的数值，其实就是 neg_mean_ _squared_error去掉负号的数字。

交叉验证

简单看看回归树是怎样工作的

'''
回归树的数据集要求是要连续的数据
'''
from sklearn.datasets import load_boston 
from sklearn.model_selection import cross_val_score 
from sklearn.tree import DecisionTreeRegressor 

boston = load_boston()
regressor = DecisionTreeRegressor(random_state=0) # 实例化
'''
交又验证 cross_ val_scorel的用法
参数：
第一个参数：适用于任何算法实例化后的模型，不仅仅适用于决策树
第二个参数：不需要划分测试集和训练集的完整矩阵数据集
第三个参数：完整的数据集标签
第四个参数：交叉验证的次数 （cv=10，代表交叉验证请做10次，一份做测试集，九份做训练集，循环往复）
第五个参数：使用父的均方误差来衡量这个模型（评估模型）
'''
cross_val_score(regressor,boston.data,boston.target,cv=10
				#,scoring="neg_mean_squared_error"
				)


'''
输出：
array([ 0.52939335,  0.60461936, -1.60907519,  0.4356399 ,  0.77280671,
        0.40597035,  0.23656049,  0.38709149, -2.06488186, -0.95162992])
'''

交叉验证是用来观察模型的稳定性的一种方法，我们将数据划分为n份，依次使用其中一份作为测试集，其他n-1份作为训练集，多次计算模型的精确性来评估模型的平均准确程度。训练集和测试集的划分会干扰模型的结果，因此用交叉验证n次的结果求出的平均值，是对模型效果的一个更好的度量

3.2实例：一维回归的图像绘制

接下来我们到二维平面上来观察决策树是怎样拟合一条曲线的。我们用回归树来拟合正弦曲线，并添加一些噪声来观察回归树的表现。

1. 导入需要的库

import numpy as np 
from sklearn.tree import DecisionTreeRegressor 
import matplotlib.pyplot as plt

2. 刨建一条含有噪声的正弦曲线

在这一步，我们的基本思路是，先创建一组随机的，分布在0-5上的横坐标轴的取值（x），然后将这一组值放到sin数中去生成纵坐标的值（y），接着再到y上去添加噪声。全程我们会使用 numpy库来为我们生成这个正弦曲线。

完整代码：

import numpy as np 
from sklearn.tree import DecisionTreeRegressor 
import matplotlib.pyplot as plt

# 生成随机数种子
rng= np.random.RandomState(1) 
'''
rand(80,1) 生成80行1列随机数组（二维数组）分类树中不允许是一维的数组  
5* 数据在0-5之间 
np.sort从小到大排序
'''
x= np.sort(5* rng. rand(80,1),axis=0) 
#生成正弦函数
y =np.sin(x).ravel()  # 标签只能是一维的使用 ravel()降维
# 在正弦曲线上增加噪声
y[::5]+=3*(0.5-rng.rand(16)) # 添加噪声 y[::5] 步长 每5个取出一个，让着数随机加上一个数

#np, random,rand（数组结构），生成随机数组的函数

#了解降维函数rave1O的用法
np. random.random((2,1))
np. random.random((2,1)).ravel()
np. random.random((2,1)).ravel().shape

'''
输出：(2,)
'''

代码解析：

rng= np.random.RandomState(1) # 生成随机数种子
rng.rand(10)
'''
array([4.17022005e-01, 7.20324493e-01, 1.14374817e-04, 3.02332573e-01,
       1.46755891e-01, 9.23385948e-02, 1.86260211e-01, 3.45560727e-01,
       3.96767474e-01, 5.38816734e-01])
'''
rng.rand(2,3)
'''
array([[0.41919451, 0.6852195 , 0.20445225],
       [0.87811744, 0.02738759, 0.67046751]])
'''
5*rng.rand(80,1)
'''
array([[4.00372284],
       [4.84130788],
       [1.56712089],
       [3.46161308],
       [4.38194576],
       [4.47303332],
       [0.42522106],
       [0.19527392],
       [0.8491521 ],
       [4.39071252],
       [0.49173417],
       [2.10553813],
       [4.78944765],
       [2.66582642],
       [3.45938557],
       [1.57757816],
       [3.43250464],
       [4.17312836],
       [0.09144139],
       [3.75072157],
       [4.94430544],
       [3.74082827],
       [1.40221996],
       [3.94639664],
       [0.51613003],
       ...
'''
np.sort(5* rng. rand(80,1),axis=0) 
'''
array([[0.0627799 ],
       ...
       [0.97714741],
       [1.01646617],
       ...
       [1.93930322],
       [2.43172555],
       ...
       [2.92879636],
       [3.03164731],
       ...
       [3.5649449 ],
       [4.00316336],
       ...
'''

rng= np.random.RandomState(1) # 生成随机数种子
x= np.sort(5* rng. rand(80,1),axis=0) 
np.sin(x)
'''
array([[ 5.71874056e-04],
       [ 9.13140084e-02],
       [ 9.66835240e-02],
       [ 1.36510390e-01],
       [ 1.94035253e-01],
       [ 2.47178482e-01],
       [ 2.63658285e-01],
       [ 4.12522163e-01],
       [ 4.45464463e-01],
       [ 4.72155294e-01],
       [ 4.89635918e-01],
       [ 4.93517994e-01],
       [ 6.05300129e-01],
       [ 6.34509327e-01],
       [ 6.41446076e-01],
       [ 6.45696215e-01],
       [ 6.69579781e-01],
       [ 6.69681218e-01],
       ...
'''
np.sin(x).ravel() # 降维
'''
array([ 5.71874056e-04,  9.13140084e-02,  9.66835240e-02,  1.36510390e-01,
        1.94035253e-01,  2.47178482e-01,  2.63658285e-01,  4.12522163e-01,
        ...
       -8.58088588e-01, -9.45898629e-01, -9.48667092e-01, -9.48706725e-01,
       -9.71490934e-01, -9.80974083e-01, -9.85684166e-01, -9.97200340e-01,
       -9.99943982e-01, -9.97032450e-01, -9.91701462e-01, -9.73227696e-01])
'''
np.sin(x).ravel().shape # 降维
'''
(80,)
'''
y=np.sin(x).ravel()
plt.figure()
plt.scatter(x,y,s=20,edgecolors="black",c="darkorange",label="data")

y
'''
array([ 5.71874056e-04,  9.13140084e-02,  9.66835240e-02,  1.36510390e-01,
        ...
       -8.58088588e-01, -9.45898629e-01, -9.48667092e-01, -9.48706725e-01,
       -9.71490934e-01, -9.80974083e-01, -9.85684166e-01, -9.97200340e-01,
       -9.99943982e-01, -9.97032450e-01, -9.91701462e-01, -9.73227696e-01])
'''
y[::5]+=3*(0.5-rng.rand(16))
y[::5]
'''
array([-1.20956669,  0.02614002,  1.98102494,  0.29426147,  1.32246237,
        0.90465031, -0.15784929,  1.30602073, -0.36520521,  0.075758  ,
        1.64636654, -1.57513771, -1.01666557, -2.21266332,  0.03431138,
        0.09139241])
'''
0.5-rng.rand(16) # 取值在-0.5~0.5之间
'''
array([-0.38330609, -0.12367221, -0.25094243,  0.15110166,  0.23007211,
       -0.39588622,  0.07190881, -0.46484005, -0.1634415 , -0.12169572,
        0.38525403, -0.44948926,  0.05008787, -0.07838961,  0.0918632 ,
        0.26297302])
'''
3*(0.5-rng.rand(16))
'''
array([-1.29778639, -0.59045448,  1.30199948, -0.76638916, -0.76162857,
       -1.26907361, -0.63457428,  1.12718711,  1.4403596 ,  1.42136704,
        1.41508054,  0.7613668 , -1.08008385, -0.11649319, -0.15846594,
       -1.02609268])
'''
y
'''
array([-1.20956669,  0.09131401,  0.09668352,  0.13651039,  0.19403525,
        0.02614002,  0.26365828,  0.41252216,  0.44546446,  0.47215529,
        ...
        0.03431138, -0.94870673, -0.97149093, -0.98097408, -0.98568417,
        0.09139241, -0.99994398, -0.99703245, -0.99170146, -0.9732277 ])
'''
plt.figure()
plt.scatter(x,y,s=20,edgecolors="black",c="darkorange",label="data")

3.实例化&训练模型

regr_1= DecisionTreeRegressor(max_depth=2)
regr_2= DecisionTreeRegressor(max_depth=5)
regr_1.fit(x,y)
regr_2.fit(x,y)
'''
DecisionTreeRegressor(criterion='mse', max_depth=5, max_features=None,
           max_leaf_nodes=None, min_impurity_decrease=0.0,
           min_impurity_split=None, min_samples_leaf=1,
           min_samples_split=2, min_weight_fraction_leaf=0.0,
           presort=False, random_state=None, splitter='best')

# np.arange(开始点,结束点,步长)
'''

4.测试集导入模型，预测结果

# np.arange(开始点,结束点,步长)
x_test= np.arange(0.0,5.0,0.01)[:,np.newaxis] # 升维
y_1= regr_1.predict(x_test)
y_2=regr_2.predict(x_test)

#np, arrange（开始点，结束点，步长）生成有序数组的函数
#了解增维切片np. newall的用法
l =np. array([1,2,3,4])
l.shape
'''
(4,)
'''
l[:,np.newaxis]
'''
array([[1],
       [2],
       [3],
       [4]])
'''
l[:,np.newaxis].shape
'''
(4, 1)
'''
l[np.newaxis,:].shape
'''
(1, 4)
'''

5.绘制图像

plt.figure()
# 散点图 edgecolors：边框的颜色 c:点的颜色 label:标签
plt.scatter(x,y,s=20,edgecolors="black",c="darkorange",label="data")

# 折线图 linewidth：线宽
plt.plot(x_test,y_1,color="cornflowerblue",label="max_depth=2",linewidth=2)

plt.plot(x_test,y_2,color="yellowgreen",label="max_depth=5",linewidth=2)

plt.xlabel("data")
plt.ylabel("target")
plt.title("Decision Tree Regresssion")

plt.figure()
plt.scatter(x,y,s=20,edgecolors="black",c="darkorange",label="data")
plt.plot(x_test,y_1,color="cornflowerblue",label="max_depth=2",linewidth=2)
plt.plot(x_test,y_2,color="yellowgreen",label="max_depth=5",linewidth=2)
plt.xlabel("data")
plt.ylabel("target")
plt.title("Decision Tree Regresssion")
plt.legend()
plt.show()

可见，回归树学习了近似正弦曲线的局部线性回归。我们可以看到，如果树的最大深度（由 max_depth参数控制）设置得太高，则决策树学习得太精细，它从训练数据中学了很多细节，包括噪声得呈现，从而使模型偏离真实的正弦曲线，形成过拟合。

4实例：泰坦尼克号幸存者的预测

泰坦尼克号的沉没是世界上最严重的海难事故之一，今天我们通过分类树模型来预测一下哪些人可能成为幸存者。

数据集来着httρs://www.kaggle.com/c/titanic，数据集包含两个csv格式文件，data为我们接下来要使用的数据，test为 haggle提供的测试集接下来我们就来执行我们的代码

# 1.导入所需的库
import pandas as pd
from sklearn.tree import DecisionTreeClassifier
import matplotlib.pyplot as plt
from sklearn.model_selection import GridSearchCV
from sklearn.model_selection import train_test_split
from sklearn.model_selection import cross_val_score 

# 2.导入数据库，探索数据
data = pd.read_csv(r"C:\Users\ASUS\sklearn\mytrain.csv")

data.head() # 默认显示5行 也可以添加自定义显示行数 data.head(10)
# 输出结果如下图

data.info()
'''
输出：

RangeIndex: 891 entries, 0 to 890
Data columns (total 13 columns):
 #   Column       Non-Null Count  Dtype  
---  ------       --------------  -----  
 0   PassengerId  891 non-null    int64  
 1   Survived     891 non-null    int64  
 2   Pclass       891 non-null    int64  
 3   Name         891 non-null    object 
 4   Sex          891 non-null    object 
 5   Age          714 non-null    float64
 6   SibSp        891 non-null    int64  
 7   Parch        891 non-null    int64  
 8   Ticket       891 non-null    object 
 9   Fare         891 non-null    float64
 10  Cabin        204 non-null    object 
 11  Embarked     889 non-null    object 
 12  Unnamed: 12  0 non-null      float64
dtypes: float64(3), int64(5), object(5)
memory usage: 90.6+ KB
'''

# 筛选特征

# 删除与存活不相关的无用信息：仓号，名字，票号，缺失值太多的数据
# inplace=True 表示删除之后覆盖原表，axis=1 表示对列进行操作，若=0 则对行进行操作
data.drop(['Cabin','Name','Ticket','Unnamed: 12'],inplace=True,axis=1)

# data = data.drop(['Cabin','Name','Ticket'],inplace=False)

data.head()
# 输出结果如下图

# 处理缺失值
data["Age"] = data["Age"].fillna(data["Age"].mean())
data.info()
'''
输出：

RangeIndex: 891 entries, 0 to 890
Data columns (total 9 columns):
 #   Column       Non-Null Count  Dtype  
---  ------       --------------  -----  
 0   PassengerId  891 non-null    int64  
 1   Survived     891 non-null    int64  
 2   Pclass       891 non-null    int64  
 3   Sex          891 non-null    object 
 4   Age          891 non-null    float64
 5   SibSp        891 non-null    int64  
 6   Parch        891 non-null    int64  
 7   Fare         891 non-null    float64
 8   Embarked     889 non-null    object 
dtypes: float64(2), int64(5), object(2)
memory usage: 62.8+ KB

'''

# dropna() 表示删除存在缺失值的行，默认asix=0 
data = data.dropna()
data.info()
'''
输出：

Int64Index: 889 entries, 0 to 890
Data columns (total 9 columns):
 #   Column       Non-Null Count  Dtype  
---  ------       --------------  -----  
 0   PassengerId  889 non-null    int64  
 1   Survived     889 non-null    int64  
 2   Pclass       889 non-null    int64  
 3   Sex          889 non-null    object 
 4   Age          889 non-null    float64
 5   SibSp        889 non-null    int64  
 6   Parch        889 non-null    int64  
 7   Fare         889 non-null    float64
 8   Embarked     889 non-null    object 
dtypes: float64(2), int64(5), object(2)
memory usage: 69.5+ KB

'''

data["Embarked"].unique().tolist()
'''
输出：['S', 'C', 'Q']
'''
# 将文字转换为数字
labels = data["Embarked"].unique().tolist()
data["Embarked"] = data["Embarked"].apply(lambda x:labels.index(x))
"S"-->labels.index("S")
data
# 输出如下图所示

(data["Sex"]=="male").astype("int")
'''
输出：
0      1
1      0
2      0
3      0
4      1
      ..
886    1
887    0
888    0
889    1
890    1
Name: Sex, Length: 889, dtype: int32
'''
data.loc[:,"Sex"] = (data["Sex"]=="male").astype("int")
data.head()
# 输出如下图所示

data.columns !="Survived"
'''
输出：array([ True, False,  True,  True,  True,  True,  True,  True,  True])
'''
x = data.iloc[:,data.columns !="Survived"]
x
# 输出如下图所示

y = data.iloc[:,data.columns =="Survived"]
Xtrain,Xtest,Ytrain,Ytest = train_test_split(x,y,test_size=0.3) # 划分数据集
# 将索引排序
for i in [Xtrain,Xtest,Ytrain,Ytest]:
    i.index = range(i.shape[0])
# 3.建立模型
clf = DecisionTreeClassifier(random_state=25)
clf = clf.fit(Xtrain,Ytrain)
score = clf.score(Xtest,Ytest) # 返回预测的准确度

score
'''
输出：0.7191011235955056
'''

clf = DecisionTreeClassifier(random_state=25)
# 交叉验证
score = cross_val_score(clf,x,y,cv=10).mean()

tr = []
te = []
for i in range(10):
    clf = DecisionTreeClassifier(random_state=25
                                ,max_depth=i+1
                                )
    clf = clf.fit(Xtrain,Ytrain)
    score_tr = clf.score(Xtrain,Ytrain)
    score_te = cross_val_score(clf,x,y,cv=10).mean()
    tr.append(score_tr)
    te.append(score_te)
    
print(max(te))
plt.plot(range(1,11),tr,color ="red",label ="train")
plt.plot(range(1,11),te,color ="blue",label ="test")
plt.xticks(range(1,11))
plt.legend()
plt.show()
# 输出结果如下图所示

tr = []
te = []
for i in range(10):
    clf = DecisionTreeClassifier(random_state=25
                                ,max_depth=i+1
                                ,criterion="entropy"
                                )
    clf = clf.fit(Xtrain,Ytrain)
    score_tr = clf.score(Xtrain,Ytrain)
    score_te = cross_val_score(clf,x,y,cv=10).mean()
    tr.append(score_tr)
    te.append(score_te)
    
print(max(te))
plt.plot(range(1,11),tr,color ="red",label ="train")
plt.plot(range(1,11),te,color ="blue",label ="test")
plt.xticks(range(1,11))
plt.legend()
plt.show()
# 输出结果如下图所示

# 4.网格搜索：能够帮助我们调整多个参数的技术，枚举技术
# 缺点：不能舍弃给定的参数，只能在给定参数的基础上找到最优
import numpy as np
gini_threholds = np.linspace(0,0.5,20)

#entropy_threholds = np.linspace(0,1,50)


# 一串参数和这些参数对应的，我们希望网格搜索来搜索的参数的取值范围
parmeters = {"criterion":("gini","entropy")
            ,"splitter":("best","random")
            ,"max_depth":[*range(1,10)]
            ,"min_samples_leaf":[*range(1,50,5)]
            ,"min_impurity_decrease":[*np.linspace(0,0.5,20)]
            }

clf = DecisionTreeClassifier(random_state=25)
GS = GridSearchCV(clf,parmeters,cv=10)
GS = GS.fit(Xtrain,Ytrain)

GS.best_params_  # 从我们输入的参数和参数取值的列表中返回最佳组合
'''
输出：
{'criterion': 'gini',
 'max_depth': 6,
 'min_impurity_decrease': 0.0,
 'min_samples_leaf': 1,
 'splitter': 'best'}
'''
GS.best_score_ #网格搜索后的模型评价标准
'''
输出：0.8215434083601286

'''

完整代码：

# 1.导入所需的库
import numpy as np
import pandas as pd
from sklearn.tree import DecisionTreeClassifier
import matplotlib.pyplot as plt
from sklearn.model_selection import GridSearchCV
from sklearn.model_selection import train_test_split
from sklearn.model_selection import cross_val_score 

# 2.导入数据库，探索数据
data = pd.read_csv(r"C:\Users\ASUS\sklearn\mytrain.csv")

# 筛选特征

# 删除与存活不相关的无用信息：仓号，名字，票号，缺失值太多的数据
# inplace=True 表示删除之后覆盖原表，axis=1 表示对列进行操作，若=0 则对行进行操作
data.drop(['Cabin','Name','Ticket','Unnamed: 12'],inplace=True,axis=1)

# data = data.drop(['Cabin','Name','Ticket'],inplace=False)

# 处理缺失值
data["Age"] = data["Age"].fillna(data["Age"].mean())
# dropna() 表示删除存在缺失值的行，默认asix=0 
data = data.dropna()

# 将文字转换为数字
labels = data["Embarked"].unique().tolist()
data["Embarked"] = data["Embarked"].apply(lambda x:labels.index(x))
"S"-->labels.index("S")
data.loc[:,"Sex"] = (data["Sex"]=="male").astype("int")

x = data.iloc[:,data.columns !="Survived"]
y = data.iloc[:,data.columns =="Survived"]
Xtrain,Xtest,Ytrain,Ytest = train_test_split(x,y,test_size=0.3) # 划分数据集
# 将索引排序
for i in [Xtrain,Xtest,Ytrain,Ytest]:
    i.index = range(i.shape[0])

# 3.建立模型
clf = DecisionTreeClassifier(random_state=25)
clf = clf.fit(Xtrain,Ytrain)
score = clf.score(Xtest,Ytest) # 返回预测的准确度
score

tr = []
te = []
for i in range(10):
    clf = DecisionTreeClassifier(random_state=25
                                ,max_depth=i+1
                                ,criterion="entropy"
                                )
    clf = clf.fit(Xtrain,Ytrain)
    score_tr = clf.score(Xtrain,Ytrain)
    score_te = cross_val_score(clf,x,y,cv=10).mean()
    tr.append(score_tr)
    te.append(score_te)
    
print(max(te))
plt.plot(range(1,11),tr,color ="red",label ="train")
plt.plot(range(1,11),te,color ="blue",label ="test")
plt.xticks(range(1,11))
plt.legend()
plt.show()


# 4.网格搜索：能够帮助我们调整多个参数的技术，枚举技术
# 缺点：不能舍弃给定的参数，只能在给定参数的基础上找到最优

gini_threholds = np.linspace(0,0.5,20)

#entropy_threholds = np.linspace(0,1,50)


# 一串参数和这些参数对应的，我们希望网格搜索来搜索的参数的取值范围
parmeters = {"criterion":("gini","entropy")
            ,"splitter":("best","random")
            ,"max_depth":[*range(1,10)]
            ,"min_samples_leaf":[*range(1,50,5)]
            ,"min_impurity_decrease":[*np.linspace(0,0.5,20)]
            }

clf = DecisionTreeClassifier(random_state=25)
GS = GridSearchCV(clf,parmeters,cv=10)
GS = GS.fit(Xtrain,Ytrain)

GS.best_params_  # 从我们输入的参数和参数取值的列表中返回最佳组合
GS.best_score_ #网格搜索后的模型评价标准