請用此 Handle URI 來引用此文件:
http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/15998
標題: | 探討結合動態漸近蒙地卡羅演算法與平行調整演算法在多峰分布抽樣的表現 Combining stochastic approximation monte carlo and parallel tempering algorithms in sampling multimodal distributions |
作者: | Kuei-Ling Huang 黃貴鈴 |
指導教授: | 張淑惠 |
關鍵字: | Metropolis-Hastings演算法,平行調整演算法,動態漸近蒙地卡羅演算法,多峰分布, Metropolis-Hastings algorithms,parallel tempering,stochastic approximation Monte Carlo algorithm,multimodal distribution, |
出版年 : | 2012 |
學位: | 碩士 |
摘要: | Metropolis-Hastings演算法是基於建構馬可夫鏈自多變量分布產生一序列隨機樣本的方法,當分布為崎嶇不平或是多峰分布的維度很高時,Metropolis-Hastings演算法很容易會陷在某個區域性的單峰分布,在文獻中,有數種改善Metropolis-Hasting演算法的演算法被提出。舉例來說,平行調整演算法是利用輔助變數改善Metropolis-Hastings演算法的模擬方法;而動態漸近蒙地卡羅演算法利用過去樣本資訊改造Metropolis-Hastings演算法。在此研究中,藉由結合此兩種方法,同時利用輔助變數與過去樣本資訊提出新的演算法。模擬研究為比較新的演算法與上述兩種演算法的表現。模擬的結果顯示動態漸近蒙地卡羅演算法在多峰分布的峰有覆蓋時表現不佳,多峰分布的峰無覆蓋時其表現佳;平行調整演算法在多峰分布的峰有無覆蓋皆表現佳,而結合兩種方法的表現依賴交換發生率及候選函數。 Metropolis-Hastings algorithm is established based on a Markov chain method to generate a series of random samples from multivariate distributions. When the distributions are rugged or the number of dimensions in multimodal distributions is high, Metropolis-Hastings algorithm is likely to be trapped locally by a certain unimodal distributions. There are several algorithms proposed to improve Metropolis-Hastings algorithm in literature. For example, parallel tempering is a simulation method which uses auxiliary variables to modify Metropolis-Hastings algorithm. Alternatively, the stochastic approximation Monte Carlo algorithm exploits the past sample information to adapt Metropolis-Hastings algorithm. In this study, a new algorithm is proposed by combing these two methods for using both information from auxiliary variables and past samples. A simulation study is conducted to investigate and compare the performance of the new algorithm and the abovementioned algorithms. The simulation results show that the performance of stochastic approximation Monte Carlo algorithm for the multimodal distribution modes coverage is poor but its performance is better for modes without covering. Parallel tempering performs well for both situations, while performance of the combined method is dependent on the exchange of incidence and proposal function. |
URI: | http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/15998 |
全文授權: | 未授權 |
顯示於系所單位: | 流行病學與預防醫學研究所 |
文件中的檔案:
檔案 | 大小 | 格式 | |
---|---|---|---|
ntu-101-1.pdf 目前未授權公開取用 | 1.95 MB | Adobe PDF |
系統中的文件,除了特別指名其著作權條款之外,均受到著作權保護,並且保留所有的權利。