A Numerical Study on Optimal Monte Carlo Algorithm for Multidimensional Integrals

Author(s):  
Venelin Todorov ◽  
Stoyan Apostolov ◽  
Ivan Dimov ◽  
Yuri Dimitrov ◽  
Stoyan Poryazov ◽  
...  
2020 ◽  
Vol 26 (3) ◽  
pp. 223-244
Author(s):  
W. John Thrasher ◽  
Michael Mascagni

AbstractIt has been shown that when using a Monte Carlo algorithm to estimate the electrostatic free energy of a biomolecule in a solution, individual random walks can become entrapped in the geometry. We examine a proposed solution, using a sharp restart during the Walk-on-Subdomains step, in more detail. We show that the point at which this solution introduces significant bias is related to properties intrinsic to the molecule being examined. We also examine two potential methods of generating a sharp restart point and show that they both cause no significant bias in the examined molecules and increase the stability of the run times of the individual walks.


Biometrika ◽  
2020 ◽  
Vol 107 (3) ◽  
pp. 745-752 ◽  
Author(s):  
Sirio Legramanti ◽  
Daniele Durante ◽  
David B Dunson

Summary The dimension of the parameter space is typically unknown in a variety of models that rely on factorizations. For example, in factor analysis the number of latent factors is not known and has to be inferred from the data. Although classical shrinkage priors are useful in such contexts, increasing shrinkage priors can provide a more effective approach that progressively penalizes expansions with growing complexity. In this article we propose a novel increasing shrinkage prior, called the cumulative shrinkage process, for the parameters that control the dimension in overcomplete formulations. Our construction has broad applicability and is based on an interpretable sequence of spike-and-slab distributions which assign increasing mass to the spike as the model complexity grows. Using factor analysis as an illustrative example, we show that this formulation has theoretical and practical advantages relative to current competitors, including an improved ability to recover the model dimension. An adaptive Markov chain Monte Carlo algorithm is proposed, and the performance gains are outlined in simulations and in an application to personality data.


Sign in / Sign up

Export Citation Format

Share Document