MCMC ，M-H採樣 示例 （pyhton版本）

初看MCMC採樣的時候，一腦子的霧水。不盡會問，馬爾科夫鏈穩定後怎麼採樣啊，都穩定了每次採樣出來的東西不是一樣的嗎？這時候特別希望有個例子能看看，怎麼採樣的，採出來的是什麼東西。很不幸，網上的例子全是連續函數的版本的，沒有離散版本的。那麼今天，我就把我的代碼分享給大家。 ''' Created on 2018年5月16日 p:輸入的概率分佈，離散情況採用元素爲概率值的數組表示 N:認爲迭代N次馬爾

>>阅读原文<<