Gibbs采样

(学习这部分内容大约需要50分钟)

摘要

Gibbs采样是一种马尔科夫连蒙特卡洛(Markov Chain Monte Carlo, MCMC)算法, 其中每个随机变量从给定剩余变量的条件分布迭代地重新采样. 它是在概率模型中执行后验推理的简单且常用的高效方法.

预备知识

学习Gibbs采样需要以下预备知识

学习目标

  • 知道Gibbs采样更新规则
  • 明白为什么 Gibbs 采样的稳态分布是模型分布.
  • 明白为什么 Gibbs 采样在变量紧密耦合时会变得低效.

核心资源

(阅读/观看其中一个资源)

免费

付费

  • 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 率.


返回贝叶斯机器学习路线图

转载于:https://www.cnblogs.com/bayesianML/p/6397780.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Gibbs采样算法是一种马尔可夫链蒙特卡洛(MCMC)方法,用于从多维概率分布采样。其代码如下: 1. 初始化样本:给定多维概率分布P(X₁, X₂, ..., Xₙ),随机初始化样本X₀ = (x₀₁, x₀₂, ..., x₀ₙ)。 2. 设置迭代次数和burn-in期:设定迭代次数T和burn-in期B。 3. 开始迭代:对于t = 1到T+B: 3.1 对于每个维度i = 1到n: a. 计算条件概率分布P(Xᵢ|X₁(t-1), X₂(t-1), ..., Xₙ(t-1)),即给定其他维度的采样样本情况下,当前维度的概率分布。 b. 从条件概率分布抽样得到新的样本值xᵢ(t)。 3.2 将得到的样本X(t)加入采样结果集中,如果t > B。 4. 返回采样结果集:返回采样得到的样本集合X = {X(B+1), X(B+2), ..., X(T+B)}。 Gibbs采样算法的核心思想是通过给定其他维度的样本情况下,逐个采样每个维度的值。每个维度的采样值都是依赖于其他维度的值,通过迭代逼近真实的多维概率分布。迭代过程中的burn-in期是为了达到一个稳定的采样状态,保证得到的样本更接近于真实概率分布。 需要注意的是,具体实现时还需要考虑到条件概率分布的计算和抽样方法。对于简单的情况,条件概率分布可以直接通过概率计算得到。对于复杂的情况,可以借助数值计算方法如Metropolis-Hastings算法等来估计条件概率分布。此外,根据具体问题,还可以对采样结果进行适当的调整和优化,以提高采样效率和精度。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值