Abstract:Observational astrophysics consists of making inferences about the Universe by comparing data and models. The credible intervals placed on model parameters are often as important as the maximum a posteriori probability values, as the intervals indicate concordance or discordance between models and with measurements from other data. Intermediate statistics (e.g. the power spectrum) are usually measured and inferences made by fitting models to these rather than the raw data, assuming that the likelihood for these statistics has multivariate Gaussian form. The covariance matrix used to calculate the likelihood is often estimated from simulations, such that it is itself a random variable. This is a standard problem in Bayesian statistics, which requires a prior to be placed on the true model parameters and covariance matrix, influencing the joint posterior distribution. As an alternative to the commonly-used Independence-Jeffreys prior, we introduce a prior that leads to a posterior that has approximately frequentist matching coverage. This is achieved by matching the covariance of the posterior to that of the distribution of true values of the parameters around the maximum likelihood values in repeated trials, under certain assumptions. Using this prior, credible intervals derived from a Bayesian analysis can be interpreted approximately as confidence intervals, containing the truth a certain proportion of the time for repeated trials. Linking frequentist and Bayesian approaches that have previously appeared in the astronomical literature, this offers a consistent and conservative approach for credible intervals quoted on model parameters for problems where the covariance matrix is itself an estimate.

Robust Bayesian regression in astronomy

A Novel Bayesian Robust Model and Its Application for Fault Detection and Automatic Supervision of Nonlinear Process

Robust Bayesian Filtering and Smoothing Using Student's t Distribution

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

Regression for Astronomical Data with Realistic Distributions, Errors and Non-linearity

Bayesian Inference and Application of Robust Growth Curve Models Using Student'S T Distribution

Robust Inference and Model Criticism Using Bagged Posteriors

Semi-supervised Variational Bayesian Student’s T Mixture Regression and Robust Inferential Sensor Application

Robust Radio Interferometric Calibration Using the t-Distribution

Robust LS-SVR Based on Variational Bayesian and Its Applications

Robust Bayesian Inference for Berkson and Classical Measurement Error Models

Bayesian inference of T Tauri star properties using multi-wavelength survey photometry

Practical Guidance for Bayesian Inference in Astronomy

Robust Bayesian linear regression with application to an analysis of the CODATA values for the Planck constant

Bayesian Analysis on Limiting the Student-$t$ Linear Regression Model

Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy

Bayesian Restricted Likelihood Methods: Conditioning on Insufficient Statistics in Bayesian Regression

Growing Pains: Understanding the Impact of Likelihood Uncertainty on Hierarchical Bayesian Inference for Gravitational-Wave Astronomy

Bayesian astrostatistics: a backward look to the future

Robustness and accuracy of methods for high dimensional data analysis based on Student's t statistic

Bayesian Robustness to Outliers in Linear Regression and Ratio Estimation