2014/05/05 by Kushal Kr. Dey, Kushal Kr Dey, Dey, Kushal Kr +2
Mathematics · Physics and Astronomy · #60H10 #60J05 #60J10 #60J60 #60J70 #Computation (stat.CO) #FOS: Computer and information sciences #Markov Chains and Monte Carlo Methods #Stochastic processes and statistical mechanics #Theoretical and Computational Physics #msc:60H10 #msc:60J05 #msc:60J10 #msc:60J60 #msc:60J70 #stat.CO
paper · pdf · doi:10.48550/arxiv.1405.0913
38 pages, this version will appear in Brazilian Journal of Probability and Statistics
openalex publication_date 2014/05/05 · arxiv created 2017/07/25 · arxiv updated 2017/07/26 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
Very recently, Transformation based Markov Chain Monte Carlo (TMCMC) was proposed by Dutta and Bhattcharya (2013) as a much efficient alternative to the Metropolis-Hastings algorithm, Random Walk Metropolis (RWM) algorithm, especially in high dimensions. The main advantage of this algorithm is that it simultaneously updates all components of a high dimensional parameter by some appropriate deterministic transformation of a single random variable, thereby reducing time complexity and enhancing the acceptance rate. The optimal scaling of the additive TMCMC approach has already been studied for the Gaussian proposal density by Dey and Bhattacharya(2013). In this paper, we discuss diffusion-based optimal scaling behavior for non-Gaussian proposal densities - in particular, uniform, Student's t and Cauchy proposals. We also consider diffusion based optimal scaling for non-Gaussian proposals when the target density is discontinuous. In the case of the Random Walk metropolis (RWM) algorithm these non-regular situations have been studied by Neal and Roberts (2011) in terms of expected squared jumping distance (ESJD), but the diffusion based approach has not been considered. Although we could not formally prove our diffusion result for the Cauchy proposal, simulation based results led us to a conjecture that the diffusion result still holds for the Cauchy case. We compare our diffusion based TMCMC approach with that of ESJD based RWM approach for the very challenging Cauchy proposal case, showing that our former approach clearly outperforms the latter.