马尔可夫链蒙特卡罗方法及其R实现

来源 :南京大学 | 被引量 : 0次 | 上传用户:choasy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
普通Monte Carlo方法通过随机模拟为一些复杂积分的计算提供了便利,但是当碰到不容易直接采样的时候,将马尔可夫链的概念引入到采样过程中,利用非周期且不可约的马氏链可收敛到平稳不变分布的特性,在马氏链运行足够长之后,认为产生的马氏链序列来源于目标分布的随机样本,这样便产生了Markov Chain Monte Carlo方法,即所谓的MCMC。本文重点梳理和研究了MCMC的理论基础和构造原理,使得该方法能被更好的认识和接受。本文第一部分从简单Monte Carlo模拟出发,介绍了其原理,并分析了缩小Monte Carlo积分误差的方法。第二部分介绍MCMC方法的理论基础、构建方法及收敛性分析。其中第一块内容结合马氏链的数学知识和推理,分析了MCMC方法的理论基础,第二块主要梳理和分析了MCMC方法的构造原理,重点介绍了Metropolis-Hastings算法和Gibbs抽样。在分析MCMC算法的同时,使用开源的统计软件R语言编写了算法的实现代码,得到了较好的分析结果,包括图形诊断和参数估计结果。第三块探讨了MCMC方法的收敛性和误差问题。第三部分介绍了MCMC在贝叶斯统计中的应用,并推导了Logistic回归模型中参数关于样本的后验分布f(β|y)对应Metropolis-Hastings算法的形式,并用R语言对一个Logistic案例和一个更一般化的贝叶斯模型案例进行计算,得到较好的诊断图形和参数估计。本文最后对MCMC的并行化算法进行了展望和尝试,以期对MCMC算法在大样本或大数据形式下有好的应用发展。
其他文献
自2009年乘用车销量跃居世界第一以来,中国已连续九年蝉联全球最大汽车市场。伴随乘用车市场进入成熟期阶段,SUV及豪华车市场作为发展较晚的车型细分市场,开始进入快速发展的
有责任透彻地理解社会主义理论问题──学习马克思恩格斯社会主义观札记之一严枝马克思和恩格斯的社会主义观,是以他们所创立的科学社会主义理论为基础的,同时也是这一理论的基
案件质量是检察机关执法办案的生命线,而案件质量评查是评价案件质量优劣的标尺,是权力相互监督制约的有效方式,是保障司法公正的有力举措,是契合检察改革发展的现实需要。评
随着我国经济的快速发展,特别是”金工业的快速发展,大量的工业废渣对环境产生了重大的影响。工业废渣在道路工程中的应用不但可以为道路工程提供大量原材料,而且还对废渣的处理
当代科技发展迅速,物联网服务不断增多。要想实现物联网的大规模应用,将信息进行聚合已变得越来越重要。服务执行环境将从大而全的个体服务模式转向泛在多终端节点设备多而广