ravest

Radial Velocity fitting with Bayesian model comparison

Features:

  • Model exoplanets and host stars, to simulate RV data for given orbital and instrumental parameters

  • Fit RV data with MCMC to explore posterior distributions for parameters - including Gaussian Processes for stellar activity

  • Bayesian Model Comparison using the Learned Harmonic Mean Estimator from harmonic

  • Visualise/animate the star’s orbit (coming soon!)

Check out the tutorial notebooks at the online documentation to see examples of how to do all of these!

Installation

It should be as simple as

$ pip install ravest

JAX is a requirement for ravest (and harmonic), so you may want to consult the JAX installation docs if you want GPU or TPU support (tl;dr: install JAX first according to those instructions, then install ravest on top).

Usage

For an introduction to modelling planetary and stellar data, see the example modelling notebook for ravest.model.

For an example of how to fit a model to RV data, see the example fitting notebook where we fit some ELODIE data for 51 Peg b.

For an example of how to use a Gaussian Process to mitigate stellar variability, see the example GP notebook where we use a quasiperiodic kernel on HARPS data for K2-229.

For an example on using the Learned Harmonic Mean Estimator from harmonic to compare two competing RV models by estimating the Bayesian evidence \(\mathcal{Z}\) and Bayes Factors, see the example harmonic notebook where we compare a one-planet and two-planet fit for TOI-544.

If you have any questions, check the Frequently Asked Questions, raise an issue on Github, or email me and I’ll be happy to help.

Contributing

Interested in contributing? Check out the contributing guidelines. Please note that this project is released with a Code of Conduct. By contributing to this project, you agree to abide by its terms.

Acknowledgements

Ravest makes use of the following open-source packages:

  • NumPy for numerical computing

  • SciPy for scientific computing algorithms

  • Matplotlib for plotting and visualisation

  • Astropy for astronomical calculations and utilities

  • pandas for data manipulation

  • tqdm for progress bars

  • emcee for MCMC sampling

  • corner for visualising posterior distributions

  • tinygp for Gaussian Process modelling, which requires JAX

  • harmonic for Bayesian evidence estimation via the Learned Harmonic Mean Estimator

License

ravest was created by Ross Dobson. It is licensed under the terms of the MIT license.