Gibbs采样
(学习这部分内容大约需要50分钟)
摘要
Gibbs采样是一种马尔科夫连蒙特卡洛(Markov Chain Monte Carlo, MCMC)算法, 其中每个随机变量从给定剩余变量的条件分布迭代地重新采样. 它是在概率模型中执行后验推理的简单且常用的高效方法.
预备知识
学习Gibbs采样需要以下预备知识
- 条件分布: Gibbs采样根据条件分布定义.
- 马尔科夫蒙特卡洛(Markov chain Monte Carlo, MCMC): Gibbs采样是一种MCMC算法.
- 马尔科夫随机场(Markov random fields): 图模型的结构表明为什么Gibbs更新计算会高效.
学习目标
- 知道Gibbs采样更新规则
- 明白为什么 Gibbs 采样的稳态分布是模型分布.
- 明白为什么 Gibbs 采样在变量紧密耦合时会变得低效.
核心资源
(阅读/观看其中一个资源)
免费
Coursera: Probabilistic Graphical Models (2013)
简介: 概率图模型的一门在线课程
位置: Lecture "Gibbs sampling" ](https://class.coursera.org/pgm/lecture/74)
网站
作者: Daphne Koller
其他依赖知识:备注:
- 点击"Preview"观看视频
Bayesian Reasoning and Machine Learning
简介: 一门研究生机器学习课程- Section 27.3, Gibbs sampling, up to 27.3.2, "Structured Gibbs sampling" (pages 547-548)
- Section 27.3.3, "Remarks" (pages 549-550)
作者: David Barber
其他依赖知识:
付费
- Pattern Recognition and Machine Learning(PRML)
简介: 一本研究生机器学习课程, 专注于贝叶斯方法
位置: Section 11.3, pages 542-546
网站
作者: Christopher M. Bishop
其他依赖知识: - Probabilistic Graphical Models: Principles and Techniques
简介: 一本非常全面的概率AI研究生教材
位置: Sections 12.3.1 (pages 505-507) and 12.3.3 (pages 512-515)
网站
作者: Daphne Koller,Nir Friedman
其他依赖知识:
增补资源
(以下为可选内容, 你可能会发现它们很有用)
免费
- Information Theory, Inference, and Learning Algorithms
简介: 一本机器学习和信息论的研究生教材
位置: Section 29.5, "Gibbs sampling," pages 370-371
网站
作者: David MacKay
其他依赖知识:- Metropolis-Hastings算法
- Machine learning summer school: Markov chain Monte Carlo (2009)
简介: 一个 MCMC 方法视频教程
位置: 69:40 to 77:34
网站
作者: Iain Murray - Computational Cognition Cheat Sheets (2013)
简介: 认知科学家写的一些笔记
位置: Bayesian Inference: Gibbs Sampling
网站
付费
- Machine Learning: a Probabilistic Perspective(MLAPP)
简介: 一本非常全面的研究生机器学习教材
位置: Section 24.2-24.2.2, pages 838-839
网站
作者: Kevin P. Murphy
相关知识
Gibbs采样可以被看作 Metropolis-Hastings 算法的特例
简单 Gibbs 采样常常 mix 得很慢. 以下是一些改进版本:
- block Gibbs sampling: 一次采样多个随机变量.
- collapsed Gibbs sampling: 一些随机变量被解析地积分掉了
Slice sampling 是一类特殊的 Gibbs 采样, 适合从没有解析形式的单变量分布中采样.
我们可以使用spectral graph理论分析 mixing 率.
Gibbs采样的更多相关文章
- MC, MCMC, Gibbs采样 原理&实现(in R)
本文用讲一下指定分布的随机抽样方法:MC(Monte Carlo), MC(Markov Chain), MCMC(Markov Chain Monte Carlo)的基本原理,并用R语言实现了几个例 ...
- 如何做Gibbs采样(how to do gibbs-sampling)
原文地址:<如何做Gibbs采样(how to do gibbs-sampling)> 随机模拟 随机模拟(或者统计模拟)方法最早有数学家乌拉姆提出,又称做蒙特卡洛方法.蒙特卡洛是一个著名 ...
- MCMC(四)Gibbs采样
MCMC(一)蒙特卡罗方法 MCMC(二)马尔科夫链 MCMC(三)MCMC采样和M-H采样 MCMC(四)Gibbs采样 在MCMC(三)MCMC采样和M-H采样中,我们讲到了M-H采样已经可以很好 ...
- 文本主题模型之LDA(二) LDA求解之Gibbs采样算法
文本主题模型之LDA(一) LDA基础 文本主题模型之LDA(二) LDA求解之Gibbs采样算法 文本主题模型之LDA(三) LDA求解之变分推断EM算法(TODO) 本文是LDA主题模型的第二篇, ...
- 关于LDA的gibbs采样,为什么可以获得正确的样本?
算法里面是随机初始了一个分布,然后进行采样,然后根据每次采样的结果去更新分布,之后接着采样直到收敛. 1.首先明确一下MCMC方法. 当我们面对一个未知或者复杂的分布时,我们经常使用MCMC方法来进行 ...
- Gibbs 采样的应用
Gibbs 采样的最大作用在于使得对高维连续概率分布的抽样由复杂变得简单. 可能的应用: 计算高维连续概率分布函数的数学期望, Gibbs 采样得到 n 个值,再取均值: 比如用于 RBM:
- Gibbs 采样定理的若干证明
坐标平面上的三点,A(x1,y1),B(x1,y2),C(x2,y1),假设有概率分布 p(x,y)(P(X=x,Y=y) 联合概率),则根据联合概率与条件概率的关系,则有如下两个等式: {p(x1, ...
- MCMC:Gibbs 采样(matlab 实现)
MCMC: The Gibbs Sampler 多元高斯分布的边缘概率和条件概率 Marginal and conditional distributions of multivariate norm ...
- 随机采样和随机模拟:吉布斯采样Gibbs Sampling实现文档分类
http://blog.csdn.net/pipisorry/article/details/51525308 吉布斯采样的实现问题 本文主要说明如何通过吉布斯采样进行文档分类(聚类),当然更复杂的实 ...
随机推荐
- ArcEngine下SceneControl叠加影像数据(构建三维地形)
载入完TIN数据后.须要在三维物体上描绘细节.即纹理:建立DEM表面点与纹理空间点的关系,即纹理映射.叠加影像数据就是把影像看作纹理.将其贴在地形表面,让其具有地形起伏的三维效果. 这里与GlobeC ...
- C#内置泛型委托:Action委托
1.什么是Action泛型委托 Action<T>是.NET Framework内置的泛型委托,可以使用Action<T>委托以参数形式传递方法,而不用显示声明自定义的委托.封 ...
- Spring Cloud搭建手册(2)——Spring Cloud Config
※在Dalston.SR2版本以后,均不能正常加密,如果必须使用此功能,需要降级到SR1或Camden SR7. 1.首先需要创建一个config-server工程,作为配置中心的服务器,用来与git ...
- 微信小程序——修改data里面数组某一个值
比如我现在有个data数据如下: data: { playIndex: null, courseList: [{ videoId: '0', isPlaying: false, }, { videoI ...
- Java设计模式(16)中介模式(Mediator模式)
Mediator定义:用一个中介对象来封装一系列关于对象交互行为. 为何使用Mediator模式/中介模式 各个对象之间的交互操作非常多,每个对象的行为操作都依赖彼此对方,修改一个对象的行为,同时会涉 ...
- 如何利用pyCharm编写和运行python文件
在安装python环境后,通常可以利用IDE pyCharm来编译我们的python文件.创建一个python文件夹,用pyCharm打开文件夹,在文件夹中新建一个python文件demo.py 也许 ...
- (笔记)Mysql实例:建库建表并插入数据2
drop database if exists school; // 如果存在SCHOOL则删除create database school; // 建立库SCHOOLuse school; / ...
- (笔记)Linux下检测网卡与网线连接状态
http://blog.chinaunix.net/space.php?uid=20357359&do=blog&cuid=1798479 Linux下检测网卡与网线连接状态,使用io ...
- 层层递进——宽度优先搜索(BFS)
问题引入 我们接着上次“解救小哈”的问题继续探索,不过这次是用宽度优先搜索(BFS). 注:问题来源可以点击这里 http://www.cnblogs.com/OctoptusLian/p/74296 ...
- CentOS 7 之前好好的,突然一天启动时黑屏,没有登陆界面了(配置 network-scripts 连网)
原因: 百度大神说是Gnome(一套纯粹自由的计算机软件,运行在操作系统上,提供图形桌面环境)不行了. 解决方法: 1. 重启系统,ctrl + alt + F2 进入命令行界面. 2. sudo s ...