The Efficiency of Geometric Samplers for Exoplanet Transit Timing Variation Models

07/30/2018
by   Noah W. Tuchow, et al.
0

Transit timing variations (TTVs) are a valuable tool to determine the masses and orbits of transiting planets in multi-planet systems. TTVs can be readily modeled given knowledge of the interacting planets' orbital configurations and planet-star mass ratios, but such models are highly nonlinear and difficult to invert. Markov chain Monte Carlo (MCMC) methods are often used to explore the posterior distribution for model parameters, but, due to the high correlations between parameters, nonlinearity, and potential multi-modality in the posterior, many samplers perform very inefficiently. Therefore, we assess the performance of several MCMC samplers that use varying degrees of geometric information about the target distribution. We generate synthetic datasets from multiple models, including the TTVFaster model and a simple sinusoidal model, and test the efficiencies of various MCMC samplers. We find that sampling efficiency can be greatly improved for all models by sampling from a parameter space transformed using an estimate of the covariance and means of the target distribution. No one sampler performs the best for all datasets, but several samplers, such as Differential Evolution Monte Carlo and Geometric adaptive Monte Carlo, have consistently efficient performance. For datasets with near Gaussian posteriors, Hamiltonian Monte Carlo samplers with 2 or 3 leapfrog steps obtained the highest efficiencies. Based on differences in effective sample sizes per time, we show that the right choice of sampler can improve sampling efficiencies by several orders of magnitude.

READ FULL TEXT

page 11

page 12

page 13

research
08/05/2021

The No-U-Turn Sampler as a Proposal Distribution in a Sequential Monte Carlo Sampler with a Near-Optimal L-Kernel

Markov Chain Monte Carlo (MCMC) is a powerful method for drawing samples...
research
06/05/2017

Stochastic Gradient Monomial Gamma Sampler

Recent advances in stochastic gradient techniques have made it possible ...
research
06/19/2018

Bayesian Sequential Inference in Dynamic Survival Models

Dynamic hazard models are applied to analyze time-varying effects of cov...
research
10/08/2014

Bayesian tracking and parameter learning for non-linear multiple target tracking models

We propose a new Bayesian tracking and parameter learning algorithm for ...
research
08/29/2016

On the Computational Complexity of Geometric Langevin Monte Carlo

Manifold Markov chain Monte Carlo algorithms have been introduced to sam...
research
10/22/2021

Focusing on Difficult Directions for Learning HMC Trajectory Lengths

Hamiltonian Monte Carlo (HMC) is a premier Markov Chain Monte Carlo (MCM...
research
04/04/2018

Posterior Inference for Sparse Hierarchical Non-stationary Models

Gaussian processes are valuable tools for non-parametric modelling, wher...

Please sign up or login with your details

Forgot password? Click here to reset