MCMC ,M-H采样 示例 (pyhton版本)
初看MCMC采样的时候,一脑子的雾水。不尽会问,马尔科夫链稳定后怎么采样啊,都稳定了每次采样出来的东西不是一样的吗?这时候特别希望有个例子能看看,怎么采样的,采出来的是什么东西。很不幸,网上的例子全是概率密度函数的版本的,没有多项分布版本的。那么今天,我就把我的代码分享给大家。
'''
Created on 2018年5月16日
p:输入的概率分布,离散情况采用元素为概率值的数组表示
N:认...
原创
2019-04-07 19:12:49 ·
5263 阅读 ·
1 评论