决策树(decision tree)是一个树结构(可以是二叉树或非二叉树)。其每个非叶节点表示一个特征属性上的测试,每个分支代表这个特征属性在某个值域上的输出,而每个叶节点存放一个类别。使用决策树进行决策的过程就是从根节点开始,测试待分类项中相应的特征属性,并按照其值选择输出分支,直到到达叶子节点,将叶子节点存放的类别作为决策结果。
我们构建一棵决策树的基本想法就是,我们希望决策树每个叶子节点包含的样本尽可能属于同一个类别,即结点的“纯度”越来越高
根据构建树的原则来看,即使得每个结点的纯度尽可能小,那么我们需要一些指标评价“纯度”这个概念。信息熵和基尼指数是两个常用的指标。
信息熵(information entropy)是度量样本集合纯度的常用指标;
假定当前样本集合D中第k类样本所占的比例为 pk(k=1,2,…,|Y|) ,
则D的信息熵为:
Ent(D)的值越小,D的纯度越高(约定:若p=0则plog2p=0)
一般而言,信息增益越大,则意味着用属性a来进行划分所获得的纯度提升越大:
ID3就是以信息增益为准则来选择划分属性的
实际上,信息增益对可取值数目较多的属性有所偏好(如编号,在西瓜集中若以编号为划分属性,则其信息增益最大),为减少由于偏好而带来的不利影响,C4.5算法使用增益率(gain ratio)来选择最优划分属性
其中:
称为属性a的固有值(intrinsic value),属性a的可能数目越多,则IV(a)的值通常越大
然而,增益率准则对可取值数目较少的属性有所偏好,C4.5采用的是先从候选划分属性中找出信息增益高于平均水平的属性,再从中选择增益率最高的
CART(Classification and Regression Tree)使用基尼指数(Gini index)来选择划分属性,数据集的纯度可用基尼值来度量
属性a的基尼指数定义为:
在属性集合A中寻找:
CART决策树使用基尼指数作为属性划分的标准
剪枝(pruning)是决策树学习算法对付过拟合的主要手段,基本策略有预剪枝(prepruning)和后剪枝(post-pruning)
在C4.5决策树算法当中,使用二分法对连续的数值进行处理:我们可以考察包含n-1个元素的候选划分点集合
我们将每个区间的中位点作为候选划分点,然后我们使用想离散值属性一样来考察这些划分点,选取最优的划分点进行样本集合的划分,例如:
对上图表格当中的例子而言,设置密度为:
根据Gain的计算公式可以得到属性”密度“的信息增益位0.262,对应于划分点0.381。同时按照之前的离散值的计算方法,计算离散属性的信息增益的值:
Gain(D ,色泽) = 0.109; Gain(D ,根蒂) = 0.143;
Gain(D ,敲声) = 0.141; Gain(D ,纹理) = 0.381;
Gain(D ,脐部) = 0.289; Gain(D , 触感) = 0.006;
Gain(D ,密度) = 0.262; Gain(D ,含糖率) = 0.349.
可以发现纹理的信息增益是最大的,所以我们选择”纹理“作为根节点作为划分属性,然后每个结点划分过程递归进行,最终生成如图所示的决策树:
一些数据由于敏感等原因,部分数据可能会出现缺失的情况,例如下面的情况:
在决策树的C4.5算法当中,我们使用了没有缺失值的样本子集进行树的构建。以上述表格为例子举例,没有缺失值的样例子集包含编号为{2,3,4,6,7,8,9,10,11,12,14,15,16,17}的14个样例(总共有17个样例)。那么相应的信息熵为:
其分别在”色泽“属性上取值为”青绿“,”乌黑“以及”浅白“的样本子集,那么有:
因此在样本子集上,其信息增益为:
那么在样本集上的”色泽“的信息增益为,要乘以其没有缺失的样例数量除以全部的样例数量:
在上述文章提及的变量为,其中每个样本的权重 wk为1:
1. 很可能对噪声数据产生过拟合。决策树越深,由噪声产生过拟合的可能性就越大。一种解决方案是对决策树进行剪枝。
# -*- coding: utf-8 -*-
"""
Created on Wed Jul 31 16:51:08 2019
@author: 86182
"""
from sklearn.datasets import load_iris
from sklearn import tree
import pydotplus
#用于划分训练集与测试集
from sklearn.model_selection import train_test_split
from sklearn.metrics import classification_report
#加载数据
iris = load_iris()
#划分训练集与测试集
(training_inputs, testing_inputs, training_classes, testing_classes)
=train_test_split(iris.data, iris.target,test_size=0.4, random_state=1)
# 构建模型
clf = tree.DecisionTreeClassifier()
clf = clf.fit(training_inputs, training_classes)
#测试值预测
y_predict = clf.predict(testing_inputs)
#预测值和测试值打分
score = classification_report(testing_classes, y_predict)
print(score)
# 保存模型
with open("iris.dot", 'w') as f:
f = tree.export_graphviz(clf, out_file=f)
# 画图,保存到pdf文件
# 设置图像参数
dot_data = tree.export_graphviz(clf, out_file=None,
feature_names=iris.feature_names,
class_names=iris.target_names,
filled=True, rounded=True,
special_characters=True)
graph = pydotplus.graph_from_dot_data(dot_data)
# 保存图像到pdf文件
graph.write_pdf("irsi.pdf")
阅读量:81
点赞量:0
收藏量:0