Dan Foreman-Mackey SEARCH CHARACTERIZATION

3 downloads 157 Views 2MB Size Report
100. 1000 ms s COMPUTATION TIME. NUMBER OF DATAPOINTS. YOU DON'T NEED TO DETREND YOUR DATA. Dan Foreman-Mackey. Center for ...
YOU DON’T NEED TO DETREND YOUR DATA Dan Foreman-Mackey

Center for Cosmology & Particle Physics, NYU dfm.io — github.com/dfm — @exoplaneteer

Gaussian process

likelihood function

p(data | model) = N [f ; µ(✓), ⌃(↵)] 3-parameter kernel function

explicit sparsity SuiteSparse1

covariance includes: observational uncertainties systematics (trends, discontinuities, etc.)

computationally

s

tractable George2

state-of-the art sparse linear algebra

an open source implementation C (+Python bindings)

COMPUTATION TIME

a generalization of standard χ² likelihood

model includes: physics (Kepler’s equation, N-body, etc.) geometry (limb darkening, etc.)



observations includes: “raw” photometry auxiliary measurements (optional)

3

O(N )

O(N

ms

)

NUMBER OF DATAPOINTS

100

t

SEARCH

⇠2

1000

CHARACTERIZATION Δ log(likelihood)

search is the process of mining light curves for new candidates. Writing this question as a hypothesis test lets us incorporate a better likelihood function and search without median filtering the data. We find: Kepler is sensitive to Earth-analog transits! anywhere that one might use a χ² likelihood function, this Gaussian process is a drop-in replacement. For characterizing planet candidates, all we need to do is trivially update the computation of the probability function in the MCMC loop. It’s even easy to marginalize the kernel parameters!

phase

3

log(period)

check out code.dfm.io for open source code and demos

2 1 0 1 2 3 time since transit [days]

3

2 1 0 1 2 3 time since transit [days]

1 cise.ufl.edu/research/sparse/SuiteSparse 2 github.com/dfm/george