Geodesic Lagrangian Monte Carlo over the space of positive definite matrices: with application to Bayesian spectral density estimation

Andrew Holbrook, Shiwei Lan, Alexander Vandenberg-Rodes, Babak Shahbaba

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

Abstract

We present geodesic Lagrangian Monte Carlo, an extension of Hamiltonian Monte Carlo for sampling from posterior distributions defined on general Riemannian manifolds. We apply this new algorithm to Bayesian inference on symmetric or Hermitian positive definite (PD) matrices. To do so, we exploit the Riemannian structure induced by Cartan's canonical metric. The geodesics that correspond to this metric are available in closed-form and–within the context of Lagrangian Monte Carlo–provide a principled way to travel around the space of PD matrices. Our method improves Bayesian inference on such matrices by allowing for a broad range of priors, so we are not limited to conjugate priors only. In the context of spectral density estimation, we use the (non-conjugate) complex reference prior as an example modelling option made available by the algorithm. Results based on simulated and real-world multivariate time series are presented in this context, and future directions are outlined.

Original languageEnglish (US)
Pages (from-to)982-1002
Number of pages21
JournalJournal of Statistical Computation and Simulation
Volume88
Issue number5
DOIs
StatePublished - Mar 24 2018

Keywords

  • HMC
  • Riemannian geometry
  • spectral analysis

ASJC Scopus subject areas

  • Statistics and Probability
  • Modeling and Simulation
  • Statistics, Probability and Uncertainty
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Geodesic Lagrangian Monte Carlo over the space of positive definite matrices: with application to Bayesian spectral density estimation'. Together they form a unique fingerprint.

Cite this