决策树(Decison Tree)是一种基于树结构决策的机器学习算法。
决策树采用分而治之(Divide and Conquer)策略,以一系列的子决策决定分类结果。
一般的,一颗决策树包含一个根结点、若干个子结点和若干个叶结点。根结点包含样本全集;子结点对应属性划分,包含划分样本;叶结点对应决策结果,包含决策样本。从根结点到每个叶结点的路径对应一个判定测试序列。
学习目的:决策树学习的目的是为了产生一颗泛化能力强,即处理未见示例能力强的决策树。 显然,决策树的生成是一个递归过程。其核心是最优划分属性的选择,在决策树基本算法中,有以下三种情形导致递归返回:
当前结点包含的样本全属于同一类别,无需划分。
所有样本在所有属性值相同,或属性集为空,无法划分,该结点类别设定为所含样本最多的类别。
当前结点包含的样本集合为空,不能划分。该节点类别设置为其父节点所含样本最多的类别。
一般而言,随着划分过程不断进行,我们希望决策树的分支结点所包含的样本尽可能属于同一类别,即结点的“纯度”(purity)越来越高
“信息熵”(information entropy):度量样本集合纯度最常用的一种最常用的指标(可以狭隘的理解为一种混乱程度的度量) 假定当前样本集合 D D D中第 k k k类样本所占的比例 p k ( k = 1 , 2 , … , ∣ Y ∣ ) p_{k}(k=1,2,…,∣Y∣) pk(k=1,2,…,∣Y∣) 则D的信息熵定义为:
E n t ( D ) = − ∑ k = 1 ∣ Y ∣ p k l o g 2 p k Ent(D)=−∑^{∣Y∣}_{k=1}p_{k}log_{2}{p_k} Ent(D)=−k=1∑∣Y∣pklog2pk
E n t ( D ) Ent(D) Ent(D)值越小,则 D D D的纯度越高 计算信息熵约定:若 p = 0 p=0 p=0,则 p l o g 2 p = 0 p log 2 p = 0 plog_2p=0 p\log_2p=0 plog2p=0plog2p=0 E n t ( D ) Ent(D) Ent(D)最小值为 0 0 0,最大值为 l o g 2 ∣ Y ∣ log2∣\mathcal{Y}∣ log2∣Y∣
假定离散属性 a a a有 V V V个可能取值 { a 1 , a 2 , … , a V } \{a1,a2,…,aV\} {a1,a2,…,aV} 若使用 a a a来对样本集 D D D划分,则会产生 V V V个分支结点,其中第 v v v个分支结点包含了 D D D中所有在属性 a a a上取值为 a v a^{\mathcal{v}} av的样本,记为 D v D^{\mathcal{v}} Dv。再考虑到不同的分支结点所包含的样本数不同,给分支结点赋予权重 D v D \dfrac{D^{\mathcal{v}}}{D} DDv,即样本数越多的分支结点的影响越大,于是可计算出用属性 a a a对样本集 D D D进行划分所获得的“信息增益”(information gain): G a i n ( D , a ) = E n t ( D ) − ∑ v = 1 V ∣ D v ∣ ∣ D ∣ E n t ( D v ) Gain(D,a)=Ent(D)−∑_{v=1}^{V}\frac{|D_v∣}{∣D∣}Ent(D^v) Gain(D,a)=Ent(D)−v=1∑V∣D∣∣Dv∣Ent(Dv)
一般而言,信息增益越大,则意味着使用属性a来进行划分所获得的的“纯度提升”越大
因此,我们可以根据信息增益来选择结点,一般选择信息增益最大的属性作为划分结点,即: a ∗ = a r g m a x a ϵ A G a i n ( D , a ) a∗=argmax_{aϵA} Gain(D,a) a∗=argmaxaϵAGain(D,a)
著名的ID3决策树学习算法就是以信息增益为准则来选择划分属性。
实际上,由于信息增益准则对可取值数目较多的属性有所偏好,为减少这种偏好可能带来的不利影响,著名的C4.5决策树算法不直接使用信息增益,而是使用“增益率”来选择最优划分属性。增益率定义为: G a i n r a t i o ( D , a ) = G a i n ( D , a ) I V ( a ) Gain_{ratio}(D,a)=\frac{Gain(D,a)}{IV(a)} Gainratio(D,a)=IV(a)Gain(D,a) 其中: I V ( a ) = − ∑ v = 1 V ∣ D v ∣ ∣ D ∣ l o g 2 ∣ D v ∣ ∣ D ∣ IV(a)=−∑_{v=1}^{V} \frac{∣Dv∣}{∣D∣}log_2\frac{∣Dv∣}{∣D∣} IV(a)=−v=1∑V∣D∣∣Dv∣log2∣D∣∣Dv∣ 称为属性a的“固有值”( i n t r i n s i c v a l u e intrinsic value intrinsicvalue).
须注意的是,由于增益率准则对可取值数目较少的属性有所偏好,因此C4.5算法采用先从候选划分属性中找出信息增益高于平均水平的属性,再从中选择增益率最高的方法来选择最优划分属性
CART(Classification and Regression)决策树使用基尼指数(Gini Index)来选择划分属性
数据集D的纯度可用基尼值来度量: G i n i ( D ) = ∑ k = 1 ∣ Y ∣ ∑ k ′ ≠ k p k p k ′ Gini(D)=∑_{k=1}^{\mathcal{|Y|}}∑_{k′ \neq k}p_{k}p_{k′} Gini(D)=k=1∑∣Y∣k′̸=k∑pkpk′
直观来说, G i n i ( D ) Gini(D) Gini(D)反映了从数据集 D D D中随机抽取两个样本,其类别标记不一致的概率。因此, G i n i ( D ) Gini(D) Gini(D)越小,则数据集 D D D的纯度越高。
属性 a a a的基尼指数定义为: G i n i i n d e x ( D , a ) = ∑ v = 1 V ∣ D v ∣ ∣ D ∣ G i n i ( D v ) Gini_{index}(D,a)=∑_{v=1}^{V}\frac{∣Dv∣}{∣D∣}Gini(D^v) Giniindex(D,a)=v=1∑V∣D∣∣Dv∣Gini(Dv) 于是,我们在划分属性集合A中,选择那个使得划分后基尼指数最小的属性作为最优划分属性,即: a ∗ = a r g m i n a ϵ A G i n i _ i n d e x ( D , a ) a∗=argmin_{aϵA}Gini\_index(D,a) a∗=argminaϵAGini_index(D,a)
信息熵代表了混乱程度。信息熵越小,信息增益越大,纯度越大 基尼值表示了类别不一致的概率,基尼值越小,纯度越大
“剪枝”(Pruning)是决策树学习算法用于防止模型过拟合(Overfitting)的主要手段。
决策树剪枝的基本策略有"预剪枝"(prepruning)和"后剪枝"(post-pruning)。
预剪枝是指在决策树生成过程中,对每个结点在划分前先进行估计,若当前结点不能带来泛化能力的提升,则停止划分并将该结点标记为叶结点。
仅有一层划分的决策树,称为“决策树桩(decision stump)”
预剪枝基于贪心本质,预划分当前结点,减少了决策树的分支
优点:
显著降低了过拟合的风险显著减少了决策树的训练时间开销和测试时间开销缺点:
数据集可能存在当前划分验证集精度低,但后续划分显著提高的情形,无法得到最优决策树;给预剪枝决策树带来了欠拟合的风险;后剪枝就是先从训练集中生成一颗完整的决策树,然后自下向上对非叶结点进行考察,若将该结点对应的子树替换为叶结点能带来决策树泛化能力的提升,则将该子树替换为叶结点。
优点:
保留了更多分支,泛化性能增强降低了欠拟合的风险缺点:
先从训练集生成一颗完整的决策树,训练时间开销和测试时间开销大到目前为止我们仅讨论了基于离散属性来生成决策树。但是在现实学习任务中常会遇到连续属性, 由于连续属性的可取数目不再有限,因此不能根据连续属性的可取值来对结点进行划分。此时,连续属性离散化技术派上用场。最简单的策略是采用二分法(bi-partition)对连续属性进行处理,这正是C4.5算法中采用的机制 。
给定样本集D和属性a,假定在D上出现了n个不同的取值,将这些值从小到大进行排序,记为 { a 1 , a 2 , … , a n } \{a1,a2,…,an \} {a1,a2,…,an},基于划分点 t t t可将 D D D分为子集 D t − D_{t}^{-} Dt−和 D t + D_{t}^{+} Dt+,其中, D t − D_{t}^{-} Dt−包含那些在属性 a a a取值上不大于 t t t的样本,而 D t + D_{t}^{+} Dt+则包含那些在属性a上大于t的样本。
显然,对相邻的属性取值 a i a^i ai和 a i + 1 a^{i+1} ai+1来说, t t t在区间 [ a i , a i + 1 ] [ai,ai+1] [ai,ai+1]中取任意值所产生的划分结果相同。因此,对于连续属性 a a a,我们可考察包含 n − 1 n-1 n−1个元素的候选划分点集合: T a = { a i + a i + 1 2 ∣ 1 ≤ i ≤ n − 1 } Ta=\{\frac{a^i+a^{i+1}}{2}∣1≤i≤n−1\} Ta={2ai+ai+1∣1≤i≤n−1}
计算纯度的方式跟之前一致,但是将中位点值替换为划分属性值
现实任务中常会遇到不完整样本,即样本的某些属性值缺失
考虑利用有缺失属性值的不完整样本来进行学习,需解决以下两个问题:
如何在属性值缺失的情况下进行划分属性选择给定划分属性,若样本在该属性的值缺失,如何对该样本进行划分问题一 对于训练集 D D D和属性 a a a,令 D ~ \tilde{D} D~表示 D D D中在属性 a a a上没有缺失值的样本子集,使用 D ~ \tilde{D} D~来进行划分属性选择,假定属性 a a a有 V V V个可取值 { a 1 , a 2 , ⋯   , a V } \{a^1, a^2, \cdots, a^V\} {a1,a2,⋯,aV},则:
D ~ v \tilde{D}^{v} D~v表示 D ~ \tilde{D} D~在属性a上取值为a^{v}的样本子集(共有V个值) D ~ v \tilde{D}^{v} D~v表示 D ~ \tilde{D} D~中属于第 K K K类的样本子集(共 Y \mathcal{Y} Y类) 显然有 D ~ ⋃ k = 1 Y D ~ k = ⋃ k = 1 V D ~ v \tilde{D} \bigcup_{k=1}^{\mathcal{Y}} \tilde{D}_k = \bigcup_{k=1}^{\mathcal{V}} \tilde{D}^{v} D~k=1⋃YD~k=k=1⋃VD~v 假定我们为每个样本 x x x赋予一个权重 w x w_x wx ,并定义: ρ = ∑ x ϵ D ~ w x ∑ x ϵ D w x \rho=\frac{\sum_{x{\epsilon}{\tilde{D}}}w_x}{\sum_{x{\epsilon}{D}}w_x} ρ=∑xϵDwx∑xϵD~wxp ~ k = ∑ x ϵ D ~ k w x ∑ x ϵ D ~ w x ( 1 ≤ k ≤ ∣ Y ∣ ) \tilde{\mathcal{p}}_k=\frac{\sum_{x{\epsilon}{\tilde{D}}_k}w_x}{\sum_{x{\epsilon}{\tilde{D}}}w_x} (1{\leq}k{\leq}{|\mathcal{Y}|}) p~k=∑xϵD~wx∑xϵD~kwx(1≤k≤∣Y∣)
r ~ v = ∑ x ϵ D ~ v w x ∑ x ϵ D ~ w x ( 1 ≤ v ≤ V ) \tilde{\mathcal{r}}_{\mathcal{v}}=\frac{\sum_{x{\epsilon}{\tilde{D}}^{\mathcal{v}}}w_x}{\sum_{x{\epsilon}{\tilde{D}}}w_x}(1{\leq}{v}{\leq}{V}) r~v=∑xϵD~wx∑xϵD~vwx(1≤v≤V)
直观的看, 对属性 a , ρ a, \rho a,ρ表示无缺失值样本所占的比例; p k ~ \tilde{p_k} pk~表示无缺失值样本中第k类所占的比例; r v ~ \tilde{r_v} rv~表示无缺失值样本中在属性 a a a上取值 a v a^{v} av的样本所占的比例; 为权重值,默认为1,在属性缺失的样本会同时进入所有分支,并将权重调整为各分支占的比重。
基于上述定义,我们可将信息增益的计算式推广为:
G a i n ( D , a ) = ρ × G a i n ( D ~ , a ) Gain(D,a)=ρ×Gain(\tilde{D},a) Gain(D,a)=ρ×Gain(D~,a) 其中: E n t ( D ~ ) = − ∑ k = 1 ∣ Y ∣ p ~ k l o g 2 p ~ k Ent(\tilde{D})=−∑_{k=1}^{∣Y∣}\tilde{p}_{k}log_2\tilde{p}_k Ent(D~)=−k=1∑∣Y∣p~klog2p~k 问题二
若样本 x x xxx在划分属性 a a a上的取值已知,则将 x x x划入与其取值对应的子结点,且样本权值在子结点中保持为 w x w_x wx。若样本 x x x在划分属性 a a a上的取值未知,则将 x x x同时划入所有结点,且样本权值在与属性值 a v a^v av 对应的子结点中调整为 r ~ v ⋅ w x \tilde{r}_v⋅w_x r~v⋅wx ;直观的看,这是让同一个样本以不同的概率划入到不同的子结点中去。
将样本集合对应多维空间,每个属性对应一个维度,分类就是在不同类空间寻找边界。单变量决策树的分类边界是由若干个与坐标轴平行的分段组成。
若能使用斜的划分边界,则决策树模型将大为简化, “多变量决策树“(multivariate decision tree)就是能实现这样的“斜划分”,甚至更复杂划分的决策树 在多变量决策树的学习过程中,不是为每个非叶结点寻找一个最优划分属性,而是试图建立一个合适的线性分类器。