Slice Smapling算法

今天看到一種新的Sampling技術，頗爲興奮。配上adaptive rejection sampling, Gibbs抽樣就可以「無敵」了。   大致說說Slice sampling的算法： 我們想從一個密度分佈爲f(x)的函數中抽樣 初始化：任意選擇x0，步長w for i= 1:n  固定xi,在[0,f(xi)]區間上均勻抽樣，得到yi  固定yi，在函數值大於yi的x區間（{x: f(x

>>阅读原文<<