机器学习 —— 概率图模型(推理:MAP)
MAP 是最大后验概率的缩写。后验概率指的是当有一定观测结果的情况下,对其他随机变量进行推理。假设随机变量的集合为X ,观察到的变量为 e, W = X-e , AP = P(W|e). 后验概率和联合概率是不同的两个概念。事实上,后验概率更接近推理本身的“意义”,并且被越来越多的用于诊断系统中。在医疗诊断系统中,存在包括病症,症状等许多随机变量,使用VE或者消息传递之类的推理手段确实可以获得每个随机变量的概率以及某些随机变量的联合概率(一个Scope的概率)。但实际上,如果面对某些很少出现的症状以及对应病症,联合概率密度函数并不合适,比如病人红细胞减少的概率非常小,但红细胞减少很大几率上对应“白血病”。如果求联合分布则会得到一个较小的概率值(因为得白血病的人本来就不多,加上有其他症状干扰),但是如果使用后验概率,则能排除此干扰。
1. 积化和的最大化
最大后验概率是一种推理手段。w = argmax P(W|e)。 W是尚未观测到的随机变量,使得此概率最大的意义是在获得某观测后,推断最可能发生了什么。这个公式把MAP变成了一个优化问题。P(X)实际上是一系列 P(scope)的乘积。在取对数的情况下,积就变成了和,对数的底是可以随意选择的。demo example 如下图所示.
如果 PHI_k (Dk) 是链状分解的情况下,可以采用变量分离最大化的方式来求取其最大值。如图所示:
由变量边际类比可知,还可由消息传递的方式来计算最终结果。实际上这里的最大值代替了边际。利用消息传递的方式计算最大后验概率如图所示:
最大化执行完毕后,得到的是各个变量的“势”,以及使得“势”最大的变量组合取值。简而言之,就是一组推断的结果。
2. NP完备的MAP问题
2.1 对应问题
对应问题是在工程中经常碰到的问题。例如多目视觉中的配准,同一个物体被不同像素观测到。那么我们关心的两个不同的像素值各是多少。这个像素值本质是随机变量,物体是观测量。该问题可以直观的表示为下图:
2.2 联系势能与最小割
联系势能往往用于一种特殊的MRF中,例如 pair wise MRF(成对的马尔科夫网络)。这种网络有一个特点就是随机变量之间两两相连,随机变量(用于图像分割问题则是label,三维重构中则是深度)可以分为两个部分,第一是图的顶点,也就是图片的像素,第二是图的边,也就是两个相连像素的邻接关系(xixj)。pair wise 马尔科夫场就是这样定义变量scope的,原因不明。
针对相邻的两个变量,可以定义Phi_(xi,xj)为联系势能,其表述如下图所示。
其中,a+d>b+c表示相邻的两个像素应该相近。此外再利用罚函数防止边缘误判。最后通过最小割算法对连接势能进行分割则可达成图像分割的目的。
2.3 基数因子,离散因子
因子法比较简单粗暴,直接对某种随机变量的特定组合,一般是scope的组合,定义一个因子,希望可能发生的就把因子定义的尽可能大。其中基数因子利用直接统计随机变量取值的方式,如相邻4像素均为蓝色,则label可能是天空的概率就极高。
离散因子法是针对某些特定组合定义较高的势能,而其他组合则给较低的势能,这种势能定义方法往往用于输入判断。如果输入某个拼音组合和字典中的某个词对应,则提示改词。
2.4 凸因子
凸因子是在一串随机变量中给出一个凸紧集的方式,本质上是强化目标的连续性。例如图像分割问题中,相邻像素的label极大可能是相同的。以此为约束可提高运算速度。
机器学习 —— 概率图模型(推理:MAP)的更多相关文章
- 机器学习 —— 概率图模型(Homework: MCMC)
除了精确推理之外,我们还有非精确推理的手段来对概率图单个变量的分布进行求解.在很多情况下,概率图无法简化成团树,或者简化成团树后单个团中随机变量数目较多,会导致团树标定的效率低下.以图像分割为例,如果 ...
- 机器学习 —— 概率图模型(Homework: Exact Inference)
在前三周的作业中,我构造了概率图模型并调用第三方的求解器对器进行了求解,最终获得了每个随机变量的分布(有向图),最大后验分布(双向图).本周作业的主要内容就是自行编写概率图模型的求解器.实际上,从根本 ...
- 机器学习 —— 概率图模型(Homework: Representation)
前两周的作业主要是关于Factor以及有向图的构造,但是概率图模型中还有一种更强大的武器——双向图(无向图.Markov Network).与有向图不同,双向图可以描述两个var之间相互作用以及联系. ...
- 机器学习 —— 概率图模型(Homework: StructuredCPD)
Week2的作业主要是关于概率图模型的构造,主要任务可以分为两个部分:1.构造CPD;2.构造Graph.对于有向图而言,在获得单个节点的CPD之后就可依据图对Combine CPD进行构造.在获得C ...
- 机器学习 —— 概率图模型(Homework: CRF Learning)
概率图模型的作业越往后变得越来越有趣了.当然,难度也是指数级别的上涨啊,以至于我用了两个周末才完成秋名山神秘车牌的寻找,啊不,CRF模型的训练. 条件随机场是一种强大的PGM,其可以对各种特征进行建模 ...
- 机器学习 —— 概率图模型(Homework: Factors)
Talk is cheap, I show you the code 第一章的作业主要是关于PGM的因子操作.实际上,因子是整个概率图的核心.对于有向图而言,因子对应的是CPD(条件分布):对无向图而 ...
- 机器学习 —— 概率图模型(学习:CRF与MRF)
在概率图模型中,有一类很重要的模型称为条件随机场.这种模型广泛的应用于标签—样本(特征)对应问题.与MRF不同,CRF计算的是“条件概率”.故其表达式与MRF在分母上是不一样的. 如图所示,CRF只对 ...
- 机器学习 —— 概率图模型(CPD)
CPD是conditional probability distribution的缩写,翻译成中文叫做 条件概率分布.在概率图中,条件概率分布是一个非常重要的概念.因为概率图研究的是随机变量之间的练习 ...
- 机器学习 —— 概率图模型(Homework: Structure Learning)
概率图的学习真的要接近尾声了啊,了解的越多越发感受到它的强大.这周的作业本质上是data mining.从数据中学习PGM的结构和参数,完全使用数据驱动 —— No structure, No par ...
随机推荐
- edge.js架起node.js和.net互操作桥梁
今天要介绍的是edge.js这个github上刚兴起的开源项目,它可以让node.js和.net之间在in-process下互操作..net版本在4.5及以上,因为.net4.5带来的Task,asy ...
- 升级Mac osx 10.10 后安装php 的mcrypt extension
昨天折腾到凌晨四点多,今天又到这个点,终于成功了,心情比较复杂,先记录一下,改日再整理,此刻头昏脑胀,只想蒙头大睡. 发的问题,还没人回复,自己找到答案了,英语比较蹩脚,凑合看.. http://st ...
- VB数据库经典实例总结(一)
先让大家看一张图.随后讲解..... 敲完五个例子之后的心情是非常好的.并没有想象中的那么难,深究它的话大致思路就是.: 建立数据库 --->利用VB导出数据 --->供人们使用.. 因为 ...
- IIS 分析器错误消息: 未能加载类型“_Default”
还会出现不能加载程序集,如: using Model;using BLL; 均报错 但是在vs2012 iis express中调试正常, 按照网上方法: 我将vs web项目下的目录全部cop ...
- ubuntu下修改ip重启系统ip不变
今天同学问我ubuntu下ip如何写死,我想起这周在公司我们队长也问过我,我就在这把我实验的方法说一下. 打开终端: sudo vim /etc/network/interfaces 然后按如下修改: ...
- 【UVA】【11427】玩纸牌
数学期望 也是刘汝佳老师白书上的例题……感觉思路很神奇啊 //UVA 11427 #include<cmath> #include<cstdio> #include<cs ...
- leetcode Largest Rectangle in Histogram 解法二
上一篇文章讲了该题的一个解法.后来又发现一个更好的解法. 首先依旧考虑一个升序的数列,例如1,2,3,4,5.那么它的最大矩形显然是有5种可能,即 1*5,2*4,3*3,4*2,1*5.所以最大的矩 ...
- django在nginx uwsgi和tornado异步方案在项目中的体验
原创作品,允许转载,转载时请务必以超链接形式标明文章 原始出处 .作者信息和本声明.否则将追究法律责任.http://rfyiamcool.blog.51cto.com/1030776/1397495 ...
- [转载]Spring Bean Configuration Inheritance
转自: http://www.mkyong.com/spring/spring-bean-configuration-inheritance/ In Spring, the inheritance i ...
- POJ 2771 Guardian of Decency(求最大点独立集)
该题反过来想:将所有可能发生恋爱关系的男女配对,那么可以带出去的人数应该等于这个二分图的最大独立集 先要做一下预处理,把不符合要求的双方先求出来, company[i][j]表示i.j四个标准都不符合 ...