2016-02-26 56 views
2

我是新來JAGS(只是另一個吉布斯採樣),我想知道是否有可能延長它從隨機排列的空間來樣?隨機置換爲MCMC

我問的原因是,我遇到了關於擴展JAGS的本教程,其中一個所需函數「logDensity」具有所需的輸入作爲double,但要使用排列我需要輸入double。

http://www.cidlab.com/prints/wabersich2013extending.pdf

我也沒有隻JAGS工作,如果任何人有另一項建議,將工作我將不勝感激任何輸入。

+0

我可以編寫採樣的自己,我意識到這是一個簡單的任務。然而,我想用dclone(與尖齒使用R R庫),以使並行處理。 – steve3051980

回答

0

如果置換超過18個或更少元素,則可以通過計算其lexicographic rank來將置換表示爲雙精度。 (這是可能通過計算排名爲64位整數,然後重新詮釋位作爲雙重擠壓出一對夫婦更多的元素。)

+0

謝謝你,我也只是看萊默代碼也。我認爲這會工作。謝謝 – steve3051980