A pragmatic Bayesian perspective on correlation analysis: The exoplanetary gravity - stellar activity case

Abstract

We apply the Bayesian framework to assess the presence of a correlation between two quantities. To do so, we estimate the probability distribution of the parameter of interest, $\rho$, characterizing the strength of the correlation. We provide an implementation of these ideas and concepts using python programming language and the pyMC module in a very short ($\sim$130 lines of code, heavily commented) and user-friendly program. We used this tool to assess the presence and properties of the correlation between planetary surface gravity and stellar activity level as measured by the log($R_{\mathrm{HK}}^{\prime }$) indicator. The results of the Bayesian analysis are qualitatively similar to those obtained via p-value analysis, and support the presence of a correlation in the data. The results are more robust in their derivation and more informative, revealing interesting features such as asymmetric posterior distributions or markedly different credible intervals, and allowing for a deeper exploration. We encourage the reader interested in this kind of problem to apply our code to his/her own scientific problems. The full understanding of what the Bayesian framework is can only be gained through the insight that comes by handling priors, assessing the convergence of Monte Carlo runs, and a multitude of other practical problems. We hope to contribute so that Bayesian analysis becomes a tool in the toolkit of researchers, and they understand by experience its advantages and limitations.

10 pages, 1 figure, accepted for publication in the peer-reviewed proceedings of the conference "Habitability in the Universe: From the Early Earth to Exoplanets"

References in corpus (3)

• Objective priors for the bivariate normal model
• Revisiting the correlation between stellar activity and planetary surface gravity
• On the correlation between stellar chromospheric flux and the surface gravity of close-in planets

Cited by in corpus (19)

• Hierarchical Bayesian calibration of tidal orbit decay rates among hot Jupiters
• The HARPS search for southern extra-solar planets. XLIII. A compact system of four super-Earth planets orbiting HD 215152
• A pair of Sub-Neptunes transiting the bright K-dwarf TOI-1064 characterised with CHEOPS
• Asteroseismic modelling of solar-type stars: internal systematics from input physics and surface correction methods
• The discovery of WASP-151b, WASP-153b, WASP-156b: Insights on giant planet migration and the upper boundary of the Neptunian desert
• The correlation between photometric variability and radial velocity jitter, based on TESS and HARPS observations
• Understanding stellar activity-induced radial velocity jitter using simultaneous K2 photometry and HARPS RV measurements
• Activity induced variation in spin-orbit angles as derived from Rossiter-McLaughlin measurements
• The HD 137496 system: A dense, hot super-Mercury and a cold Jupiter
• K2-291 b: A rocky super-Earth in a 2.2 day orbit
• HD 38858: a solar type star with a $\sim$10.8 yr activity cycle. Searching for Balmer and metallic line variations
• The SOPHIE search for northern extrasolar planets. XVIII: Six new cold Jupiters, including one of the most eccentric exoplanet orbits
• Ross 128 -- GL 447A possible activity cycle for a slow-rotating fully-convective star
• CHEOPS and TESS view of the ultra-short period super-Earth TOI-561 b
• Correlation between activity indicators: H$\alpha$ and Ca II lines in M-dwarf stars
• Evidence of a substellar companion around a very young T Tauri star
• The inflated, eccentric warm Jupiter TOI-4914 b orbiting a metal-poor star, and the hot Jupiters TOI-2714 b and TOI-2981 b
• Tracking Advanced Planetary Systems (TAPAS) with HARPS-N. V.: A Massive Jupiter orbiting the very low metallicity giant star BD+03 2562 and a possible planet around HD~103485
• New Photometry and Spectroscopy of DW Cnc