Bayesian Analysis for the Ranking of Transits

Similar documents
Approximate Bayesian Computation. Alireza Shafaei - April 2016

Hierarchical Bayesian Modeling with Ensemble MCMC. Eric B. Ford (Penn State) Bayesian Computing for Astronomical Data Analysis June 12, 2014

An Introduction to Markov Chain Monte Carlo

Nested Sampling: Introduction and Implementation

Discussion on Bayesian Model Selection and Parameter Estimation in Extragalactic Astronomy by Martin Weinberg

Particle Filters for Visual Tracking

Overview. Monte Carlo Methods. Statistics & Bayesian Inference Lecture 3. Situation At End Of Last Week

MCMC Methods for data modeling

Chapter 1. Introduction

Probabilistic Graphical Models

1

Issues in MCMC use for Bayesian model fitting. Practical Considerations for WinBUGS Users

Quantitative Biology II!

Tutorial using BEAST v2.4.1 Troubleshooting David A. Rasmussen

10.4 Linear interpolation method Newton s method

Markov Chain Monte Carlo (part 1)

Computer vision: models, learning and inference. Chapter 10 Graphical Models

From Bayesian Analysis of Item Response Theory Models Using SAS. Full book available for purchase here.

Metropolis Light Transport

CSCI 599 Class Presenta/on. Zach Levine. Markov Chain Monte Carlo (MCMC) HMM Parameter Es/mates

Bayesian Statistics Group 8th March Slice samplers. (A very brief introduction) The basic idea

STAT 203 SOFTWARE TUTORIAL

WALLABY kinematic parametre extraction:

GT "Calcul Ensembliste"

Using the package glmbfp: a binary regression example.

Modified Metropolis-Hastings algorithm with delayed rejection

Approximate (Monte Carlo) Inference in Bayes Nets. Monte Carlo (continued)

Level-set MCMC Curve Sampling and Geometric Conditional Simulation

UCT Algorithm Circle: Probabilistic Algorithms

Optimization Methods III. The MCMC. Exercises.

Linear Modeling with Bayesian Statistics

Laboratorio di Problemi Inversi Esercitazione 4: metodi Bayesiani e importance sampling

Bayesian Estimation for Skew Normal Distributions Using Data Augmentation

CSC412: Stochastic Variational Inference. David Duvenaud

1 Methods for Posterior Simulation

GAMES Webinar: Rendering Tutorial 2. Monte Carlo Methods. Shuang Zhao

SysML model of Exoplanet Archive Functionality and Activities Solange Ramirez a, and the NASA Exoplanet Archive Team a

Bayesian estimation of optical properties of the human head via 3D structural MRI p.1

WRITING, STRUCTURING AND DEVELOPING MONTE CARLO RADIATIVE TRANSFER CODES ANTONIA BEVAN, UCL ST ANDREWS MONTE CARLO SUMMER SCHOOL 2017

AMCMC: An R interface for adaptive MCMC

A Short History of Markov Chain Monte Carlo

Clustering Relational Data using the Infinite Relational Model

Image analysis. Computer Vision and Classification Image Segmentation. 7 Image analysis

GNU MCSim Frederic Yves Bois Chair of mathematical modeling and system biology for predictive toxicology

MSA101/MVE Lecture 5

Parametric Approaches for Refractivity-From-Clutter Inversion

Short-Cut MCMC: An Alternative to Adaptation

PROBABILISTIC IMAGE SEGMENTATION FOR LOW-LEVEL MAP BUILDING IN ROBOT NAVIGATION. Timo Kostiainen, Jouko Lampinen

Sampling Using GPU Accelerated Sparse Hierarchical Models

RJaCGH, a package for analysis of

Bayesian Modelling with JAGS and R

Monte Carlo Methods and Statistical Computing: My Personal E

Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri

Samuel Coolidge, Dan Simon, Dennis Shasha, Technical Report NYU/CIMS/TR

( ylogenetics/bayesian_workshop/bayesian%20mini conference.htm#_toc )

MCMC Diagnostics. Yingbo Li MATH Clemson University. Yingbo Li (Clemson) MCMC Diagnostics MATH / 24

Learning the Structure of Deep, Sparse Graphical Models

Multi-Modal Metropolis Nested Sampling For Inspiralling Binaries

Reddit Recommendation System Daniel Poon, Yu Wu, David (Qifan) Zhang CS229, Stanford University December 11 th, 2011

Autoencoder. Representation learning (related to dictionary learning) Both the input and the output are x

Tracking Algorithms. Lecture16: Visual Tracking I. Probabilistic Tracking. Joint Probability and Graphical Model. Deterministic methods

Markov chain Monte Carlo methods

Graphical Models, Bayesian Method, Sampling, and Variational Inference

BeviMed Guide. Daniel Greene

A Heuristic Robust Approach for Real Estate Valuation in Areas with Few Transactions

Monte Carlo for Spatial Models

Warm-up as you walk in

Bayesian Multiple QTL Mapping

PRINCIPLES OF PHYLOGENETICS Spring 2008 Updated by Nick Matzke. Lab 11: MrBayes Lab

Climate model-based probabilistic wind risk assessment under future climate

Package beast. March 16, 2018

Bayesian Robust Inference of Differential Gene Expression The bridge package

Introduction to WinBUGS

Analysis of Incomplete Multivariate Data

Practical Course WS12/13 Introduction to Monte Carlo Localization

Infectious Disease Models. Angie Dobbs. A Thesis. Presented to. The University of Guelph. In partial fulfilment of requirements.

Markov chain Monte Carlo sampling

Automated Video Analysis of Crowd Behavior

Variational Methods for Graphical Models

rendering equation camera all

Debugging MCMC Code. Charles J. Geyer. April 15, 2017

Recap: The E-M algorithm. Biostatistics 615/815 Lecture 22: Gibbs Sampling. Recap - Local minimization methods

INTRO TO THE METROPOLIS ALGORITHM

Adaptive spatial resampling as a Markov chain Monte Carlo method for uncertainty quantification in seismic reservoir characterization

Bayesian Computation: Overview and Methods for Low-Dimensional Models

Package INLABMA. February 4, 2017

Human Upper Body Pose Estimation in Static Images

CS281 Section 9: Graph Models and Practical MCMC

Robust False Positive Detection for Real-Time Multi-Target Tracking

A Sample of Monte Carlo Methods in Robotics and Vision. Credits. Outline. Structure from Motion. without Correspondences

Hybrid Inference for Sensor Network Localization using a Mobile Robot

Monte Carlo Methods. CS1538: Introduction to Simulations

Package dirmcmc. R topics documented: February 27, 2017

CHAPTER 1 INTRODUCTION

Mesh segmentation. Florent Lafarge Inria Sophia Antipolis - Mediterranee

Machine Learning (BSMC-GA 4439) Wenke Liu

Simulated Annealing. Slides based on lecture by Van Larhoven

Computer Vision 2 Lecture 8

Introduction to Machine Learning Prof. Mr. Anirban Santara Department of Computer Science and Engineering Indian Institute of Technology, Kharagpur

A guided walk Metropolis algorithm

Transcription:

Bayesian Analysis for the Ranking of Transits Pascal Bordé, Marc Ollivier, Alain Léger

Layout I. What is BART and what are its objectives? II. Description of BART III.Current results IV.Conclusions and outlooks 2

I. What is BART and what are its objectives? I.1) Issues 3

I.1) Issues of the transits detected by the detection teams turn out to be false positives; 4

I.1) Issues of the transits detected by the detection teams turn out to be false positives; The follow-up observations requires lots of time, money and man-power. CoRoT-7b : 70h (100 nigths) on HARPS; 5

I.1) Issues of the transits detected by the detection teams turn out to be false positives; The follow-up observations requires lots of time, money and man-power. CoRoT-7b : 70h (100 nigths) on HARPS; Most targets with magr can t be followed (especially in RV). 6

I.1) Issues of the transits detected by the detection teams turn out to be false positives; The follow-up observations requires lots of time, money and man-power. CoRoT-7b : 70h (100 nigths) on HARPS; Most targets with magr can t be followed (especially in RV). Probabilistic validation (BLENDER, PASTIS) is a long (computer time) and complex task (CoRoT-22b). 7

I.2) Objectives of BART 8

I.2) Objectives of BART Link between detection and follow-up teams 9

I.2) Objectives of BART Link between detection and follow-up teams Point at the most promising transits for follow-up observations and probabilistic validation; 10

I.2) Objectives of BART Link between detection and follow-up teams Point at the most promising transits for follow-up observations and probabilistic validation; Make a report on what we can infer with the light curve alone on each transit. 11

I.3) What is BART Software developed in Python; 12

I.3) What is BART Software developed in Python; Automated and quick enough to be applied on all the transit of a run within a week; 13

I.3) What is BART Software developed in Python; Automated and quick enough to be applied on all the transit of a run within a week; Use bayesian model comparison to infer the probabilities for a transit to be caused by PS, EB, CPS, CEB; 14

I.3) What is BART Software developed in Python; Automated and quick enough to be applied on all the transit of a run within a week; Use bayesian model comparison to infer the probabilities for a transit to be caused by PS, EB, CPS, CEB; Produce a ranking of transit based on the probability of each transit to a be PS; 15

I.3) What is BART Software developed in Python; Automated and quick enough to be applied on all the transits of a run within a week; Use bayesian model comparison to infer the probabilities for a transit to be caused by PS, EB, CPS, CEB; Produce a ranking of transit based on the probability of each transit to a be PS; Use bayesian inference of parameters to obtain the probability density function of all the parameters (density, radius ratio,... ) of all the scenari (PS, EB,...) 16

II) Description of BART 17

II.1) Bayesian model comparison Posterior probability Likelihood Prior probability M : model (PS, EB, CEB ) d : data I : context of the measurement Marginalized likelihood (Normalization factor) 18

II.1) Bayesian model comparison Posterior probability Likelihood Prior probability M : model (PS, EB, CEB ) d : data I : context of the measurement Marginalized likelihood (Normalization factor) Eclipsing Binary single and double periode Models, M, Planetary System, Contaminating Eclipsing Binary single and double periode Contaminating Planetary System 19

II.1) Bayesian model comparison Posterior probability Likelihood Prior probability M : model (PS, EB, CEB ) d : data I : context of the measurement Marginalized likelihood (Normalization factor) Time Data, d To take advantage of the work done by the detection teams, I use as data measured on the Normalized flux transit Phase angle 20

II.2) Bayesian inference parameters Posterior probability density Likelihood Prior probability density : Parameters PS impact parameter radius ratio Marginalized likelihood (Normalization factor) CPS impact parameter radius ratio EB (P, 2P) impact parameter R prim R sec som radius ratio system density CEB (P, 2P) impact parameter R prim R sec radius ratio som system density stellar density stellar density orbital periode orbital periode orbital periode orbital periode effective effective contamination factor temperature ratio temperature ratio contamination factor 21

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 22

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 23

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 24

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 25

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 26

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 27

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 28

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 29

II.3) Exploration of parameter space and integration Fasten execution time : MCMC exploration Monte carlo integration Automatisation : Personal adaptation of P.C. Gregory (2005) Adaptative Metropolis-Hasting algorithm (acceptance rate and trace correlation feedback) more efficient with correlated parameters 30

III) Current results 31

III.1) Current results : Ranking of 22 candidates from LRc02 32

III.2) Current results : Parameters inference CoRoT-6b -- PS 33

III.3) Current results : Correlation diagrams CoRoT-6b -- PS 34

IV) Conclusions and outlook The ranking is already quite good (except for high impact parameter planets) ; I think that it could be even better with a more detailed description of contamination (contaminant list and colors); Development of more detailed prior probabilities on models and parameters 35

Thank you for attention! 36