本文转载自CSDN aliceyangxi1987的博客
通过本篇文章可以对机器学习(machine learning, ML)的常用算法有个常识性的认识,没有代码,没有复杂的理论推导,就是图解一下,知道这些算法是什么,它们怎么做的,怎么应用的,例子主要是分类问题。
每个算法都进行了简洁的讲解,便于科普,以后有时间再对单个算法做深入地解析。
今天的算法如下:
决策树
随机森林算法
逻辑回归
SVM
朴素贝叶斯
K最近邻算法
K均值算法
Adaboost 算法
神经网络
马尔可夫
根据一些 feature 进行分类,每个节点提一个问题,通过判断,将数据分为两类,再继续提问。这些问题是根据已有数据学习出来的,再投入新数据的时候,就可以根据这棵树上的问题,将数据划分到合适的叶子上。
在源数据中随机选取数据,组成几个子集:
S 矩阵是源数据,有 1-N 条数据,A B C 是feature,最后一列C是类别。
由 S 随机生成 M 个子矩阵:
这 M 个子集得到 M 个决策树
将新数据投入到这 M 个树中,得到 M 个分类结果,计数看预测成哪一类的数目最多,就将此类别作为最后的预测结果。
此算法,在扩增子、宏基因组数据分析中很常用,用于找出分组贡献率较大的OTUs或功能分类,同时方便进一步建立分类器,用于指导诊断。
当预测目标是概率这样的,值域需要满足大于等于0,小于等于1的,这个时候单纯的线性模型是做不到的,因为在定义域不在某个范围之内时,值域也超出了规定区间。
所以此时需要这样的形状的模型会比较好。
那么怎么得到这样的模型呢?
这个模型需要满足两个条件 大于等于0,小于等于1 。
大于等于0 的模型可以选择 绝对值,平方值,这里用 指数函数,一定大于0 。
小于等于1 用除法,分子是自己,分母是自身加上1,那一定是小于1的了。
再做一下变形,就得到了 logistic regression 模型
通过源数据计算可以得到相应的系数了
最后得到 logistic 的图形
support vector machine
要将两类分开,想要得到一个超平面,最优的超平面是到两类的 margin 达到最大,margin就是超平面与离它最近一点的距离,如下图,Z2>Z1,所以绿色的超平面比较好。
将这个超平面表示成一个线性方程,在线上方的一类,都大于等于1,另一类小于等于-1。
点到面的距离根据图中的公式计算。
所以得到 total margin 的表达式如下,目标是最大化这个 margin,就需要最小化分母,于是变成了一个优化问题。
举个栗子,三个点,找到最优的超平面,定义了 weight vector=(2,3)-(1,1)。
得到 weight vector 为(a,2a),将两个点代入方程,代入(2,3)另其值=1,代入(1,1)另其值=-1,求解出 a 和 截矩 w0 的值,进而得到超平面的表达式。
a 求出来后,代入(a,2a)得到的就是 support vector。
a 和 w0 代入超平面的方程就是 support vector machine。
举个在 NLP 的应用:
给一段文字,返回情感分类,这段文字的态度是positive,还是negative。
为了解决这个问题,可以只看其中的一些单词。
这段文字,将仅由一些单词和它们的计数代表。
原始问题是:给你一句话,它属于哪一类
通过 bayes rules 变成一个比较简单容易求得的问题
问题变成,这一类中这句话出现的概率是多少,当然,别忘了公式里的另外两个概率。
栗子:单词 love 在 positive 的情况下出现的概率是 0.1,在 negative 的情况下出现的概率是 0.001。
k nearest neighbours
给一个新的数据时,离它最近的 k 个点中,哪个类别多,这个数据就属于哪一类。
栗子:要区分 猫 和 狗,通过 claws 和 sound 两个feature来判断的话,圆形和三角形是已知分类的了,那么这个 star 代表的是哪一类呢?
k=3时,这三条线链接的点就是最近的三个点,那么圆形多一些,所以这个star就是属于猫。
想要将一组数据,分为三类,粉色数值大,黄色数值小 。
最开心先初始化,这里面选了最简单的 3,2,1 作为各类的初始值 。
剩下的数据里,每个都与三个初始值计算距离,然后归类到离它最近的初始值所在类别。
分好类后,计算每一类的平均值,作为新一轮的中心点。
几轮之后,分组不再变化了,就可以停止了。
adaboost 是 bosting 的方法之一。
bosting就是把若干个分类效果并不好的分类器综合起来考虑,会得到一个效果比较好的分类器。
下图,左右两个决策树,单个看是效果不怎么好的,但是把同样的数据投入进去,把两个结果加起来考虑,就会增加可信度。
adaboost 的栗子,手写识别中,在画板上可以抓取到很多 features,例如 始点的方向,始点和终点的距离等等。
training 的时候,会得到每个 feature 的 weight,例如 2 和 3 的开头部分很像,这个 feature 对分类起到的作用很小,它的权重也就会较小。
而这个 alpha 角 就具有很强的识别性,这个 feature 的权重就会较大,最后的预测结果是综合考虑这些 feature 的结果。
Adaboost算法在宏基因组领域也很常用,昨天推送的Nature文章《Salt-responsive gut commensal modulates TH17 axis and disease》中找到高盐调控高血压的菌,就是采用此算法分析得到贡献度最大的菌。
Neural Networks 适合一个input可能落入至少两个类别里。
NN 由若干层神经元,和它们之间的联系组成 。
第一层是 input 层,最后一层是 output 层。在 hidden 层 和 output 层都有自己的 classifier。
input 输入到网络中,被激活,计算的分数被传递到下一层,激活后面的神经层,最后output 层的节点上的分数代表属于各类的分数,下图例子得到分类结果为 class 1。
同样的 input 被传输到不同的节点上,之所以会得到不同的结果是因为各自节点有不同的weights 和 bias。
这也就是 forward propagation。
Markov Chains 由 state 和 transitions 组成
栗子,根据这一句话 ‘the quick brown fox jumps over the lazy dog’,要得到 markov chain。
步骤:先给每一个单词设定成一个状态,然后计算状态间转换的概率
这是一句话计算出来的概率,当你用大量文本去做统计的时候,会得到更大的状态转移矩阵,例如 the 后面可以连接的单词,及相应的概率。
生活中,键盘输入法的备选结果也是一样的原理,模型会更高级
热文:1高分文章 2不可或缺的人 3图表规范
一文读懂:1微生物组 2寄生虫益处 3进化树
必备技能:1提问 2搜索 3Endnote
文献阅读 1热心肠 2SemanticScholar 3geenmedical
扩增子分析:1图表解读 2分析流程 3统计绘图 4功能预测
科研经验:1云笔记 2云协作 3公众号
系列教程:1Biostar 2微生物组 3宏基因组
生物科普 1肠道细菌 2人体上的生命 3生命大跃进 4细胞的暗战 5人体奥秘
为鼓励读者交流、快速解决科研困难,我们建立了“宏基因组”专业讨论群,目前己有国内外100+ PI,1000+ 一线科研人员加入。参与讨论,获得专业解答,欢迎分享此文至朋友圈,并扫码加主编好友带你入群,务必备注“姓名-单位-研究方向-职称/年级”。技术问题寻求帮助,首先阅读《如何优雅的提问》学习解决问题思路,仍末解决群内讨论,问题不私聊,帮助同行。
学习16S扩增子、宏基因组科研思路和分析实战,关注“宏基因组”