Efficient cosmological parameter estimation with Hamiltonian Monte Carlo technique
- Department of Physics, Jadwin Hall, Princeton University, P.O. Box 708, Princeton, New Jersey 08542 (United States)
Traditional Markov Chain Monte Carlo methods suffer from low acceptance rate, slow mixing, and low efficiency in high dimensions. Hamiltonian Monte Carlo resolves this issue by avoiding the random walk. Hamiltonian Monte Carlo (HMC) is a Markov Chain Monte Carlo (MCMC) technique built upon the basic principle of Hamiltonian mechanics. Hamiltonian dynamics allows the chain to move along trajectories of constant energy, taking large jumps in the parameter space with relatively inexpensive computations. This new technique improves the acceptance rate by a factor of 4 while reducing the correlations and boosts up the efficiency by almost a factor of D in a D-dimensional parameter space. Therefore shorter chains will be needed for a reliable parameter estimation comparing to a traditional MCMC chain yielding the same performance. Besides that, the HMC is well suited for sampling from non-Gaussian and curved distributions which are very hard to sample from using the traditional MCMC methods. The method is very simple to code and can be easily plugged into standard parameter estimation codes such as CosmoMC. In this paper we demonstrate how the HMC can be efficiently used in cosmological parameter estimation. Also we discuss possible ways of getting good estimates of the derivatives of (the log of) posterior which is needed for HMC.
- OSTI ID:
- 21020367
- Journal Information:
- Physical Review. D, Particles Fields, Vol. 75, Issue 8; Other Information: DOI: 10.1103/PhysRevD.75.083525; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA); ISSN 0556-2821
- Country of Publication:
- United States
- Language:
- English
Similar Records
MARKOV CHAIN MONTE CARLO POSTERIOR SAMPLING WITH THE HAMILTONIAN METHOD
Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo