Skip to content

Commit 1a31c58

Browse files
committed
update 1
1 parent 4844e48 commit 1a31c58

24 files changed

+254
-0
lines changed

docs/about.md

Lines changed: 35 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,35 @@
1+
# About Me
2+
3+
I'm an exoplanet researcher focused on detecting and characterising long-period giant planets by combining precision radial-velocities with astrometry. I develop Bayesian methods and open-source tools to make model selection and evidence estimation more reliable in challenging multi-modal posteriors.
4+
5+
## Highlights
6+
7+
I built [`reddemcee`](https://reddemcee.readthedocs.io/en/latest), an adaptive parallel-tempering MCMC sampler in Python with original implementations of temperature ladders, and next-gen Bayesian evidence estimation (including a novel hybrid method which combines the best of thermodynamical integration and the stepping-stones algorithm); and [`EMPEROR`](https://astroemperor.readthedocs.io/en/latest), a Python library for joint RV+astrometry orbit fitting (multi-instrument offsets and jitters, GP noise models, stellar activity indices, etc). On the observational side, I've led and analysed programs on CORALIE and worked with HARPS data, extending the CHEPS sample to uncover Jupiter analogues (using combined RV and Hipparcos/Gaia data).
8+
9+
## Methods
10+
Computationally, I specialise in APT MCMC, evidence-driven model comparison, and stellar activity noise modelling.
11+
12+
Data-wise, my expertise lays in high-precision radial velocities and astrometry data.
13+
14+
## Now seeking
15+
Currently looking for postdoctoral positions where I can combine algorithm development with bulk data analysis (in preparation Gaia DR4/Roman/PLATO), to tackle down on important exoplanet population questions, like how common true Jupiter analogues are, to answer exoplanet formation questions, like how outer giants help shaping the inner architecture of stellar systems.
16+
17+
18+
## Relevant Links
19+
20+
=== "Pip installable libraries"
21+
* [GitHub](https://github.com/ReddTea)
22+
* [`reddemcee`](https://github.com/ReddTea/reddemcee)
23+
* [`EMPEROR`](https://github.com/ReddTea/astroemperor)
24+
* [Astronomy Utilities](https://github.com/ReddTea/reddutils)
25+
* [kepler solver in C](https://github.com/ReddTea/fast_kepler)
26+
27+
=== "Documentation"
28+
* [`reddemcee`](https://reddemcee.readthedocs.io/en/latest)
29+
* [`EMPEROR`](https://astroemperor.readthedocs.io/en/latest)
30+
31+
=== "Publications"
32+
* [`ORCID`](https://orcid.org/my-orcid?orcid=0000-0002-8770-4398)
33+
* [`ADS`](https://ui.adsabs.harvard.edu/user/libraries/Uj9qAAVWSzS_Z9ulrRovZQ)
34+
35+

docs/emperor.md

Lines changed: 62 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,62 @@
1+
[`EMPEROR`](https://astroemperor.readthedocs.io/en/latest) (Exoplanet Mcmc Parallel tEmpering for Rv Orbit Retrieval), is an open-source Python-based framework designed for the efficient detection and characterisation of exoplanets by using radial velocity (RV) methods and astrometry.
2+
3+
Its combination of performance, flexibility, and ease of use makes it a robust tool for any exoplanet detection endeavour. `EMPEROR` integrates Dynamic Nested Sampling (DNS) and Adaptive Parallel Tempering (APT) Markov Chain Monte Carlo (MCMC), supporting multiple noise models such as Gaussian Processes (GPs) and Moving Averages (MA). The framework facilitates systematic model comparison using statistical metrics, including Bayesian evidence and Bayes Information Criterion (BIC), while providing automated, publish-ready visualisations.
4+
5+
The code is easy to install and easy to use, providing full publication-ready plots and copy-pastable latex tables:
6+
7+
```python
8+
9+
import astroemperor as emp
10+
import numpy as np
11+
np.random.seed(1234)
12+
13+
sim = emp.Simulation()
14+
sim.load_data('51Peg') # folder read from /datafiles/
15+
16+
sim.engine_config['setup'] = [10, 500, 3000, 1] # ntemps, nwalkers, nsweeps, nsteps
17+
18+
sim.add_condition(['Period 1', 'limits', [3, 5]]) # short chain
19+
```
20+
21+
Some Automatic Plots:
22+
=== "Best Fit"
23+
Full model, and phase-folded per planet
24+
![Best Fit](img/emperor/model/k1.png)
25+
=== "Arviz: Trace and More"
26+
For every model block:
27+
=== "Corner"
28+
![Corner](img/emperor/trace/corner_kep.png)
29+
=== "Trace"
30+
![Trace](img/emperor/trace/trace_kep.png)
31+
=== "HDI Intervals"
32+
![hdi-intervals](img/emperor/trace/hdi_off.png)
33+
=== "Normalised Posteriors"
34+
![norm-post](img/emperor/trace/normpost_amp.png)
35+
36+
37+
=== "Posteriors"
38+
For every model block, and temperature:
39+
=== "Scatter"
40+
![Posteriors1](img/emperor/posteriors/scatter_kep_t0.png)
41+
=== "Higher temperatures!"
42+
![Posteriors2](img/emperor/posteriors/scatter_kep_t8.png)
43+
=== "Different styles 1"
44+
![Posteriors3](img/emperor/posteriors/gau_offset_t0.png)
45+
=== "Different styles 1"
46+
![Posteriors4](img/emperor/posteriors/hex_offset_t0.png)
47+
48+
=== "Estimates"
49+
For every parameter:
50+
=== "GME"
51+
![GME](img/emperor/gm/kep_amp.png)
52+
=== "Histograms"
53+
![histo](img/emperor/histograms/OffsetBlock.png)
54+
55+
56+
=== "Temperature Ladder"
57+
=== "Swap Acceptance"
58+
![temp_rates](img/emperor/betas/rates.png)
59+
=== "Ladder"
60+
![temp_ladder](img/emperor/betas/beta_ladder.png)
61+
=== "Density"
62+
![temp_dens](img/emperor/betas/density.png)

0 commit comments

Comments
 (0)