《贝叶斯统计》ppt课件_第1页
《贝叶斯统计》ppt课件_第2页
《贝叶斯统计》ppt课件_第3页
《贝叶斯统计》ppt课件_第4页
《贝叶斯统计》ppt课件_第5页
已阅读5页,还剩40页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、MCMC方法方法一、一、贝贝叶斯叶斯统计统计的框架分析的框架分析困难困难: 后验分布是复杂的、高维的分布后验分布是复杂的、高维的分布处理方法处理方法:马尔可夫链蒙特卡罗马尔可夫链蒙特卡罗MCMC方法方法后验分布后验分布 先验信息先验信息 似然函数似然函数 目前,MCMC曾经成为一种处置复杂统计问题的特别流行的工具,尤其在经常需求复杂的高维积分运算的贝叶斯分析领域更是如此。在那里,高维积分运算主要是用来求取普通方法无法得到的后验分布密度。假设合理的定义和实施,MCMC总能得到一条或几条收敛的马尔可夫链,该马尔可夫链的极限分布就是所需的后验分布一预备知识二根本思想三常用MCMC算法Gibbs抽样吉

2、布斯采样算法每次迭带的时候的一些元素曾经被跟新了,假设在更新其他的元素时不运用这些更新后的元素会呵斥一定程度新其他的元素时不运用这些更新后的元素会呵斥一定程度的浪费。现实上,的浪费。现实上, Gibbs抽样抽样 可经过在每一步都利用近似可经过在每一步都利用近似得到的其他元素的值来获得更好的效果。这种方法改良了得到的其他元素的值来获得更好的效果。这种方法改良了练的混合,换句话说,链能更加迅速,更加详尽的搜索目练的混合,换句话说,链能更加迅速,更加详尽的搜索目的分布的支撑空间。的分布的支撑空间。)(0)(02逐个生成。3添加m,前往第2步。Metropolis-Hastings抽样00000001

3、N,401,( )()pmin(1,)( )aa 假设数据是 (1 )的1000个随机数;和的初值是 和 ,用随机移动的正态分布作为建议分布做法就应该是,=建议分布为N(,I),再由它生成一个随机向量作为,然后看接受概率设先验为均匀分布,设p(x,x )= (x,x),则诊断方法1察看样本途径 产生多条马尔可夫链,察看样本途径对多个初始值产生多个马尔可夫链 样本途径是一个描画迭代数对应 的实现图。样本途径有时也称为历史图。假设链的混合不是很好,那么在很多次迭代中它会取 一样或者相近的数值。一个好的链可以快速地远离初始值,无论以何值开场。 )(m历史迭代图不收敛收敛(2)察看自相关性图自相关性图用于描画 序列在不同迭代延迟下的相关性,延迟i的自相关性是指相距i步的两迭代之间的相关性。具有较差的性质的链随着迭代延迟的添加会表现出较慢的自相关衰弱。)(m)(m)()()()()(xLxfxh)()(supxhxhMDMDdxhxEE)()(E5 . 0)(5 . 0dxhuMeMe5 . 0uMeB),()(21nxxxx1)(ULP1)()(dxhxPULUL1)()(dxhxIPIII21,)()(21xhxhdxhxPxdxhxPx10)()()()()()(111

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论