AU Mic b Samples

This folder contains data needed to recreate figures from the AU Mic b discovery paper.

Files and descriptions:

Figure 1, top panel:

F1_a.csv: two columns of TESS data, covering the first transit (green):
time, flux

F1_b.csv: two columns of TESS data, covering the second transit (red):
time, flux

F1_c.csv: transit model for TESS data (orange data)
time, model median, model 5th percentile, model 95th percentile

F1_d.csv: two columns of Spitzer data (purple dots):
time, flux

F1_e.csv: transit model for Spitzer data (orange data)
time, model median, model 5th percentile, model 95th percentile


Figure 1, lower panel:

F1_f.csv two columns of TESS data, covering the candidate planet transit (green):
time, flux

F1_g.csv transit model for TESS data (orange data)
time, model median, model 5th percentile, model 95th percentile



Extended Data Figure 2, top two panels:

ED2_a.csv: two columns of TESS data (black dots)
time, flux

ED2_b.csv: transit model (orange data)
time, model median, model 5th percentile, model 95th percentile

ED2_c.csv: GP model (green data)
time, model median, model 5th percentile, model 95th percentile

ED2_d.csv: combined model (red data)
time, model median, model 5th percentile, model 95th percentile


Extended Data Figure 2, third panel:

ED2_e.csv: two columns of Spitzer data (black dots)
time, flux

ED2_f.csv: transit model (orange data)
time, model median, model 5th percentile, model 95th percentile

ED2_g.csv: GP model (green data)
time, model median, model 5th percentile, model 95th percentile

ED2_h.csv: combined model (red data)
time, model median, model 5th percentile, model 95th percentile


Extended Data Figure 2, lower panel:

ED2_i.csv: two columns of TESS data (black dots)
time, flux

ED2_j.csv: transit model (orange data)
time, model median, model 5th percentile, model 95th percentile

ED2_k.csv: GP model (green data)
time, model median, model 5th percentile, model 95th percentile

ED2_l.csv: combined model (red data)
time, model median, model 5th percentile, model 95th percentile


Extended Data Figure 3:
Samples from the MCMC model of AU Mic b. Samples are stored in a pymc3 trace file called aumicb_pymc3.tgz. The tile will need to be untarred first you can use tar -xzvf aumicb_pymc3.tgz

This is a custom data format for PyMC3 traces. Each chain goes inside a directory, and each directory contains a metadata json file, and a numpy compressed file.
File can be read using sample code supplied. The code is in a jupyter notebook called AU_Mic_read_samples.ipynb. The full file will need to be run because the samples rely on the model being set up correctly.
Several python packages are needed to run the notebook:
numpy, matplotlib, lightkurve, exoplanet, pymc3, theano, scipy, corner, pandas, and astropy