Home News About Us Contact Contributors Disclaimer Privacy Policy Help FAQ

Quick Search
My eDoc
Session History
Support Wiki
Direct access to
document ID:

          Institute: MPI für Astronomie     Collection: Publikationen_mpia     Display Documents

ID: 742491.0, MPI für Astronomie / Publikationen_mpia
The joker: A custom Monte Carlo sampler forbinary-star and exoplanet radial velocity data
Authors:Price-Whelan, A. M.; Hogg, D. W.; Foreman-Mackey, D.; Rix, H.-W.
Date of Publication (YYYY-MM-DD):2017
Title of Journal:The Astrophysical Journal
Issue / Number:1
Start Page:id. 20 (14 pp)
Audience:Not Specified
Abstract / Description:Given sparse or low-quality radial velocity measurements of a star, there are often many qualitatively different stellar or exoplanet companion orbit models that are consistent with the data. The consequent multimodality of the likelihood function leads to extremely challenging search, optimization, and Markov chain Monte Carlo (MCMC) posterior sampling over the orbital parameters. Here we create a custom Monte Carlo sampler for sparse or noisy radial velocity measurements of two-body systems that can produce posterior samples for orbital parameters even when the likelihood function is poorly behaved. The six standard orbital parameters for a binary system can be split into four nonlinear parameters (period, eccentricity, argument of pericenter, phase) and two linear parameters (velocity amplitude, barycenter velocity). We capitalize on this by building a sampling method in which we densely sample the prior probability density function (pdf) in the nonlinear parameters and perform rejection sampling using a likelihood function marginalized over the linear parameters. With sparse or uninformative data, the sampling obtained by this rejection sampling is generally multimodal and dense. With informative data, the sampling becomes effectively unimodal but too sparse: in these cases we follow the rejection sampling with standard MCMC. The method produces correct samplings in orbital parameters for data that include as few as three epochs. The Joker can therefore be used to produce proper samplings of multimodal pdfs, which are still informative and can be used in hierarchical (population) modeling. We give some examples that show how the posterior pdf depends sensitively on the number and time coverage of the observations and their uncertainties.
Free Keywords:binaries: spectroscopic; methods: data analysis; methods: statistical; planets and satellites: fundamental parameters; surveys; techniques: radial velocities
External Publication Status:published
Document Type:Article
Communicated by:N. N.
Affiliations:MPI für Astronomie
Identifiers:ISSN:0004-637X %R 10.3847/1538-4357/aa5e50
The scope and number of records on eDoc is subject to the collection policies defined by each institute - see "info" button in the collection browse view.