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.

Probabilities in multimatrix variate distributions: an application in SARS-CoV-2

Multimatrix variate distributions

Multimatricvariate and multimatrix variate distributions based on elliptically contoured laws under real normed division algebras

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

Matrix‐variate risk measures under Wishart and gamma distributions

A Vine-Copula Based Adaptive Mcmc Sampler for Efficient Inference of Dynamical Systems

Multivariate Distributions in Non-Stationary Complex Systems I: Random Matrix Model and Formulae for Data Analysis

Conjugacy properties of multivariate unified skew-elliptical distributions

Bayesian Nonparametric Modelling of Conditional Multidimensional Dependence Structures

Matrix-Variate Dirichlet Process Priors with Applications

Monte Carlo simulation of COVID-19 pandemic using statistical physics-inspired probabilities

Multiscale statistical physics of the Human-SARS-CoV-2 interactome

Fortuitous Correlations in Molecular Dynamics Simulations: Their Harmful Influence on the Probability Distributions of the Main Principal Components

Multivariate nonnegative trigonometric sums distributions for high-dimensional multivariate circular data

Multivariate Generalized Marshall–Olkin Distributions and Copulas

Matrix variate p-value in MANOVA

Total probability and Bayes formulae for joint multidimensional-matrix Gaussian distributions

Multivariate Subexponential Distributions and Their Applications

Practical Notes on Multivariate Modeling Based on Elliptical Copulas

Compatibility and attainability of matrices of correlation-based measures of concordance

Kernel Density Matrices for Probabilistic Deep Learning