1. 随机森林使用背景
1.1 随机森林定义
随机森林是一种比较新的机器学习模型。经典的机器学习模型是神经网络,有半个多世纪的历史了。神经网络预测精确,但是计算量很大。上世纪八十年代Breiman等人发明分类树的算法(Breiman et al. 1984),通过反复二分数据进行分类或回归,计算量大大降低。2001年Breiman把分类树组合成随机森林(Breiman 2001a),即在变量(列)的使用和数据(行)的使用上进行随机化,生成很多分类树,再汇总分类树的结果。随机森林在运算量没有显著提高的前提下提高了预测精度。随机森林对多元共线性不敏感,结果对缺失数据和非平衡的数据比较稳健,可以很好地预测多达几千个解释变量的作用(Breiman 2001b),被誉为当前最好的算法之一(Iverson et al. 2008)。
随机森林顾名思义,是用随机的方式建立一个森林,森林里面有很多的决策树组成,随机森林的每一棵决策树之间是没有关联的。在得到森林之后,当有一个新的输入样本进入的时候,就让森林中的每一棵决策树分别进行一下判断,看看这个样本应该属于哪一类(对于分类算法),然后看看哪一类被选择最多,就预测这个样本为那一类。
1.2 随机森林优点
随机森林是一个最近比较火的算法,它有很多的优点:
a. 在数据集上表现良好,两个随机性的引入,使得随机森林不容易陷入过拟合
b. 在当前的很多数据集上,相对其他算法有着很大的优势,两个随机性的引入,使得随机森林具有很好的抗噪声能力
c. 它能够处理很高维度(feature很多)的数据,并且不用做特征选择,对数据集的适应能力强:既能处理离散型数据,也能处理连续型数据,数据集无需规范化
d. 可生成一个Proximities=(pij)矩阵,用于度量样本之间的相似性: pij=aij/N, aij表示样本i和j出现在随机森林中同一个叶子结点的次数,N随机森林中树的颗数
e. 在创建随机森林的时候,对generlization error使用的是无偏估计
f. 训练速度快,可以得到变量重要性排序(两种:基于OOB误分率的增加量和基于分裂时的GINI下降量
g. 在训练过程中,能够检测到feature间的互相影响
h. 容易做成并行化方法
i. 实现比较简单
1.3 随机森林应用范围
随机森林主要应用于回归和分类。本文主要探讨基于随机森林的分类问题。随机森林和使用决策树作为基本分类器的(bagging)有些类似。以决策树为基本模型的bagging在每次bootstrap放回抽样之后,产生一棵决策树,抽多少样本就生成多少棵树,在生成这些树的时候没有进行更多的干预。而随机森林也是进行bootstrap抽样,但它与bagging的区别是:在生成每棵树的时候,每个节点变量都仅仅在随机选出的少数变量中产生。因此,不但样本是随机的,连每个节点变量(Features)的产生都是随机的。
许多研究表明, 组合分类器比单一分类器的分类效果好,随机森林(random forest)是一种利用多个分类树对数据进行判别与分类的方法,它在对数据进行分类的同时,还可以给出各个变量(基因)的重要性评分,评估各个变量在分类中所起的作用。
2. 随机森林方法理论介绍
2.1 随机森林基本原理
随机森林由LeoBreiman(2001)提出,它通过自助法(bootstrap)重采样技术,从原始训练样本集N中有放回地重复随机抽取k个样本生成新的训练样本集合,然后根据自助样本集生成k个分类树组成随机森林,新数据的分类结果按分类树投票多少形成的分数而定。其实质是对决策树算法的一种改进,将多个决策树合并在一起,每棵树的建立依赖于一个独立抽取的样品,森林中的每棵树具有相同的分布,分类误差取决于每一棵树的分类能力和它们之间的相关性。特征选择采用随机的方法去分裂每一个节点,然后比较不同情况下产生的误差。能够检测到的内在估计误差、分类能力和相关性决定选择特征的数目。单棵树的分类能力可能很小,但在随机产生大量的决策树后,一个测试样品可以通过每一棵树的分类结果经统计后选择最可能的分类。
2.2 随机森林算法
2.2.1 决策树
决策树(decision tree)是一个树结构(可以是二叉树或非二叉树)。其每个非叶节点表示一个特征属性上的测试,每个分支代表这个特征属性在某个值域上的输出,而每个叶节点存放一个类别。使用决策树进行决策的过程就是从根节点开始,测试待分类项中相应的特征属性,并按照其值选择输出分支,直到到达叶子节点,将叶子节点存放的类别作为决策结果。
随机森林是用随机的方式建立一个森林,森林里面有很多的决策树组成,随机森林的每一棵决策树之间是没有关联的。在得到森林之后,当有一个新的输入样本进入的时候,就让森林中的每一棵决策树分别进行一下判断,看看这个样本应该属于哪一类,然后看看哪一类被选择最多,就预测这个样本为那一类。
在建立每一棵决策树的过程中,有两点需要注意采样与完全分裂。首先是两个随机采样的过程,random forest对输入的数据要进行行、列的采样。对于行采样,采用有放回的方式,也就是在采样得到的样本集合中,可能有重复的样本。假设输入样本为N个,那么采样的样本也为N个。这样使得在训练的时候,每一棵树的输入样本都不是全部的样本,使得相对不容易出现over-fitting。然后进行列采样,从M个feature中,选择m个(m << M)。之后就是对采样之后的数据使用完全分裂的方式建立出决策树,这样决策树的某一个叶子节点要么是无法继续分裂的,要么里面的所有样本的都是指向的同一个分类。一般很多的决策树算法都一个重要的步骤——剪枝,但是这里不这样干,由于之前的两个随机采样的过程保证了随机性,所以就算不剪枝,也不会出现over-fitting。
决策树中分裂属性的两个选择度量:
1)信息增益
随机森林模型任意样本分类的期望信息:
a) I(s1,s2,……,sm)=
∑Pi log2(pi)(i=1..m)
其中,数据集为S,m为S的分类数目,Pi≈|Si/|S|,Ci为某分类标号,Pi为任意样本属于Ci的概率,si为分类Ci上的样本数
b) I(s1,s2,……,sm)越小,s1,s2,……,sm就越有序(越纯),分类效果就越好。
c) 由属性A划分为子集的熵:
A为属性,具有V个不同的取值, S被A 划分为V 个子集s1,s2,……,sv,sij是子集sj中类Ci的样本数。E(A)= ∑(s1j+ ……+smj)/s * I(s1j,……,smj)
d) 信息增益:Gain(A)= I(s1,s2,……,sm)
E(A)
e) 分裂属性选择规则:选择具有最大信息增益的属性为分裂属性
2)基尼指数
a) 集合T包含N个类别的记录,那么其Gini指标就是pj 类别j出现的频率
b) 如果集合T分成m部分 N1 , N2 ,…, Nm 。那么这个分割的Gini就是
c)分裂属性选择规则:选择具有最小Ginisplit的属性为分裂属性(对于每个属性都要遍历所有可能的分割方法)。
2.2.3 随机森林模型的注意点
设有N个样本,每个样本有M个features,决策树们其实都是随机地接受n个样本(对行随机取样)的m个feature(对列进行随机取样),每颗决策树的m个feature相同。每颗决策树其实都是对特定的数据进行学习归纳出分类方法,而随机取样可以保证有重复样本被不同决策树分类,这样就可以对不同决策树的分类能力做个评价。
2.2.4随机森林实现过程
随机森林中的每一棵分类树为二叉树,其生成遵循自顶向下的递归分裂原则,即从根节点开始依次对训练集进行划分;在二叉树中,根节点包含全部训练数据, 按照节点
纯度最小原则,分裂为左节点和右节点,它们分别包含训练数据的一个子集,按照同样的规则节点继续分裂,直到满足分支停止规则而停止生长。若节点n上的分类数据全部来自于同一类别,则此节点的纯度I(n)=0,
纯度度量方法是Gini准则,即假设P(Xj)是节点n上属于Xj 类样本个数占训练。
具体实现过程如下:
(1)原始训练集为N,应用bootstrap法有放回地随机抽取k个新的自助样本集,并由此构建k棵分类树,每次未被抽到的样本组成了k个袋外数据;
(2)设有mall个变量,则在每一棵树的每个节点处随机抽取mtry个变量(mtry n mall),然后在mtry中选择一个最具有分类能力的变量,变量分类的阈值通过检查每一个分类点确定;
(3)每棵树最大限度地生长, 不做任何修剪;
(4)将生成的多棵分类树组成随机森林,用随机森林分类器对新的数据进行判别与分类,分类结果按树分类器的投票多少而定。
3. 随机森林应用
由于R中早就出现randomForest包了,本文主要讨论R中随机森林的应用。两个主要函数比较重要:randomForest用来构建随机森林模型,predict()使用训练后的随机森林对新数据进行预测。
3.1目标
通过随机森林的算法,根据一些特征,例如花瓣的长,宽,花萼的长宽。来预测植株的种类。
3.2 准备的数据集
iris数据集,是R语言自带的数据集。
- Sepal.Length Sepal.Width Petal.Length Petal.Width Species
- 1 5.1 3.5 1.4 0.2 setosa
- 2 4.9 3.0 1.4 0.2 setosa
- 3 4.7 3.2 1.3 0.2 setosa
- 4 4.6 3.1 1.5 0.2 setosa
- 5 5.0 3.6 1.4 0.2 setosa
- 6 5.4 3.9 1.7 0.4 setosa
- 7 4.6 3.4 1.4 0.3 setosa
- 8 5.0 3.4 1.5 0.2 setosa
- 9 4.4 2.9 1.4 0.2 setosa
- 10 4.9 3.1 1.5 0.1 setosa
- 11 5.4 3.7 1.5 0.2 setosa
- 12 4.8 3.4 1.6 0.2 setosa
- 13 4.8 3.0 1.4 0.1 setosa
- 14 4.3 3.0 1.1 0.1 setosa
- 15 5.8 4.0 1.2 0.2 setosa
- 16 5.7 4.4 1.5 0.4 setosa
- 17 5.4 3.9 1.3 0.4 setosa
- 18 5.1 3.5 1.4 0.3 setosa
- 19 5.7 3.8 1.7 0.3 setosa
- 20 5.1 3.8 1.5 0.3 setosa
- 21 5.4 3.4 1.7 0.2 setosa
- 22 5.1 3.7 1.5 0.4 setosa
- 23 4.6 3.6 1.0 0.2 setosa
- 24 5.1 3.3 1.7 0.5 setosa
- 25 4.8 3.4 1.9 0.2 setosa
- 26 5.0 3.0 1.6 0.2 setosa
- 27 5.0 3.4 1.6 0.4 setosa
- 28 5.2 3.5 1.5 0.2 setosa
- 29 5.2 3.4 1.4 0.2 setosa
- 30 4.7 3.2 1.6 0.2 setosa
- 31 4.8 3.1 1.6 0.2 setosa
- 32 5.4 3.4 1.5 0.4 setosa
- 33 5.2 4.1 1.5 0.1 setosa
- 34 5.5 4.2 1.4 0.2 setosa
- 35 4.9 3.1 1.5 0.2 setosa
- 36 5.0 3.2 1.2 0.2 setosa
- 37 5.5 3.5 1.3 0.2 setosa
- 38 4.9 3.6 1.4 0.1 setosa
- 39 4.4 3.0 1.3 0.2 setosa
- 40 5.1 3.4 1.5 0.2 setosa
- 41 5.0 3.5 1.3 0.3 setosa
- 42 4.5 2.3 1.3 0.3 setosa
- 43 4.4 3.2 1.3 0.2 setosa
- 44 5.0 3.5 1.6 0.6 setosa
- 45 5.1 3.8 1.9 0.4 setosa
- 46 4.8 3.0 1.4 0.3 setosa
- 47 5.1 3.8 1.6 0.2 setosa
- 48 4.6 3.2 1.4 0.2 setosa
- 49 5.3 3.7 1.5 0.2 setosa
- 50 5.0 3.3 1.4 0.2 setosa
- 51 7.0 3.2 4.7 1.4 versicolor
- 52 6.4 3.2 4.5 1.5 versicolor
- 53 6.9 3.1 4.9 1.5 versicolor
- 54 5.5 2.3 4.0 1.3 versicolor
- 55 6.5 2.8 4.6 1.5 versicolor
- 56 5.7 2.8 4.5 1.3 versicolor
- 57 6.3 3.3 4.7 1.6 versicolor
- 58 4.9 2.4 3.3 1.0 versicolor
- 59 6.6 2.9 4.6 1.3 versicolor
- 60 5.2 2.7 3.9 1.4 versicolor
- 61 5.0 2.0 3.5 1.0 versicolor
- 62 5.9 3.0 4.2 1.5 versicolor
- 63 6.0 2.2 4.0 1.0 versicolor
- 64 6.1 2.9 4.7 1.4 versicolor
- 65 5.6 2.9 3.6 1.3 versicolor
- 66 6.7 3.1 4.4 1.4 versicolor
- 67 5.6 3.0 4.5 1.5 versicolor
- 68 5.8 2.7 4.1 1.0 versicolor
- 69 6.2 2.2 4.5 1.5 versicolor
- 70 5.6 2.5 3.9 1.1 versicolor
- 71 5.9 3.2 4.8 1.8 versicolor
- 72 6.1 2.8 4.0 1.3 versicolor
- 73 6.3 2.5 4.9 1.5 versicolor
- 74 6.1 2.8 4.7 1.2 versicolor
- 75 6.4 2.9 4.3 1.3 versicolor
- 76 6.6 3.0 4.4 1.4 versicolor
- 77 6.8 2.8 4.8 1.4 versicolor
- 78 6.7 3.0 5.0 1.7 versicolor
- 79 6.0 2.9 4.5 1.5 versicolor
- 80 5.7 2.6 3.5 1.0 versicolor
- 81 5.5 2.4 3.8 1.1 versicolor
- 82 5.5 2.4 3.7 1.0 versicolor
- 83 5.8 2.7 3.9 1.2 versicolor
- 84 6.0 2.7 5.1 1.6 versicolor
- 85 5.4 3.0 4.5 1.5 versicolor
- 86 6.0 3.4 4.5 1.6 versicolor
- 87 6.7 3.1 4.7 1.5 versicolor
- 88 6.3 2.3 4.4 1.3 versicolor
- 89 5.6 3.0 4.1 1.3 versicolor
- 90 5.5 2.5 4.0 1.3 versicolor
- 91 5.5 2.6 4.4 1.2 versicolor
- 92 6.1 3.0 4.6 1.4 versicolor
- 93 5.8 2.6 4.0 1.2 versicolor
- 94 5.0 2.3 3.3 1.0 versicolor
- 95 5.6 2.7 4.2 1.3 versicolor
- 96 5.7 3.0 4.2 1.2 versicolor
- 97 5.7 2.9 4.2 1.3 versicolor
- 98 6.2 2.9 4.3 1.3 versicolor
- 99 5.1 2.5 3.0 1.1 versicolor
- 100 5.7 2.8 4.1 1.3 versicolor
- 101 6.3 3.3 6.0 2.5 virginica
- 102 5.8 2.7 5.1 1.9 virginica
- 103 7.1 3.0 5.9 2.1 virginica
- 104 6.3 2.9 5.6 1.8 virginica
- 105 6.5 3.0 5.8 2.2 virginica
- 106 7.6 3.0 6.6 2.1 virginica
- 107 4.9 2.5 4.5 1.7 virginica
- 108 7.3 2.9 6.3 1.8 virginica
- 109 6.7 2.5 5.8 1.8 virginica
- 110 7.2 3.6 6.1 2.5 virginica
- 111 6.5 3.2 5.1 2.0 virginica
- 112 6.4 2.7 5.3 1.9 virginica
- 113 6.8 3.0 5.5 2.1 virginica
- 114 5.7 2.5 5.0 2.0 virginica
- 115 5.8 2.8 5.1 2.4 virginica
- 116 6.4 3.2 5.3 2.3 virginica
- 117 6.5 3.0 5.5 1.8 virginica
- 118 7.7 3.8 6.7 2.2 virginica
- 119 7.7 2.6 6.9 2.3 virginica
- 120 6.0 2.2 5.0 1.5 virginica
- 121 6.9 3.2 5.7 2.3 virginica
- 122 5.6 2.8 4.9 2.0 virginica
- 123 7.7 2.8 6.7 2.0 virginica
- 124 6.3 2.7 4.9 1.8 virginica
- 125 6.7 3.3 5.7 2.1 virginica
- 126 7.2 3.2 6.0 1.8 virginica
- 127 6.2 2.8 4.8 1.8 virginica
- 128 6.1 3.0 4.9 1.8 virginica
- 129 6.4 2.8 5.6 2.1 virginica
- 130 7.2 3.0 5.8 1.6 virginica
- 131 7.4 2.8 6.1 1.9 virginica
- 132 7.9 3.8 6.4 2.0 virginica
- 133 6.4 2.8 5.6 2.2 virginica
- 134 6.3 2.8 5.1 1.5 virginica
- 135 6.1 2.6 5.6 1.4 virginica
- 136 7.7 3.0 6.1 2.3 virginica
- 137 6.3 3.4 5.6 2.4 virginica
- 138 6.4 3.1 5.5 1.8 virginica
- 139 6.0 3.0 4.8 1.8 virginica
- 140 6.9 3.1 5.4 2.1 virginica
- 141 6.7 3.1 5.6 2.4 virginica
- 142 6.9 3.1 5.1 2.3 virginica
- 143 5.8 2.7 5.1 1.9 virginica
- 144 6.8 3.2 5.9 2.3 virginica
- 145 6.7 3.3 5.7 2.5 virginica
- 146 6.7 3.0 5.2 2.3 virginica
- 147 6.3 2.5 5.0 1.9 virginica
- 148 6.5 3.0 5.2 2.0 virginica
- 149 6.2 3.4 5.4 2.3 virginica
- 150 5.9 3.0 5.1 1.8 virginica
R 源代码:
- library( ”randomForest” )
- data(iris)
- set.seed(100)
- ind=sample(2,nrow(iris),replace=TRUE,prob=c(0.8,0.2))
- iris.rf=randomForest(Species~.,iris[ind==1,],ntree=50,nPerm=10,mtry=3,proximity=TRUE,importance=TRUE)
- print(iris.rf)
- iris.pred=predict( iris.rf,iris[ind==2,] )
- table(observed=iris[ind==2,”Species”],predicted=iris.pred )
3.4 一些重要参数说明
randomForest()对训练集的数据进行处理,生成决策树
iris.rf=randomForest(Species~.,iris[ind==1,],ntree=50,nPerm=10,mtry=3,proximity=TRUE,importance=TRUE)
Species~.:代表需要预测的列,species是列的名称。
iris[ind==1,]:生成决策树的训练集
ntree:生成决策树的数目
nperm:计算importance时的重复次数
mtry:选择的分裂属性的个数
proximity=TRUE:表示生成临近矩阵
importance=TRUE:输出分裂属性的重要性
predict()
iris.pred=predict( iris.rf,iris[ind==2,] )
iris.rf:表示生成的随机森林模型
iris[ind==2,] :进行预测的测试集
3.5预测结果
- predicted
- served setosa versicolor virginica
- setosa 35 0 0
- versicolor 0 37 1
- virginica 0 3 33
- virginica 0 3 33