Bayesian Lesion Estimation with a Structured Spike-and-Slab Prior

by   Anna Menacher, et al.

Neural demyelination and brain damage accumulated in white matter appear as hyperintense areas on MRI scans in the form of lesions. Modeling binary images at the population level, where each voxel represents the existence of a lesion, plays an important role in understanding aging and inflammatory diseases. We propose a scalable hierarchical Bayesian spatial model, called BLESS, capable of handling binary responses by placing continuous spike-and-slab mixture priors on spatially-varying parameters and enforcing spatial dependency on the parameter dictating the amount of sparsity within the probability of inclusion. The use of mean-field variational inference with dynamic posterior exploration, which is an annealing-like strategy that improves optimization, allows our method to scale to large sample sizes. Our method also accounts for underestimation of posterior variance due to variational inference by providing an approximate posterior sampling approach based on Bayesian bootstrap ideas and spike-and-slab priors with random shrinkage targets. Besides accurate uncertainty quantification, this approach is capable of producing novel cluster size based imaging statistics, such as credible intervals of cluster size, and measures of reliability of cluster occurrence. Lastly, we validate our results via simulation studies and an application to the UK Biobank, a large-scale lesion mapping study with a sample size of 40,000 subjects.


page 3

page 27

page 30


Bayesian Bootstrap Spike-and-Slab LASSO

The impracticality of posterior sampling has prevented the widespread ad...

Sparse Dynamic Factor Models with Loading Selection by Variational Inference

In this paper we develop a novel approach for estimating large and spars...

Flexible mean field variational inference using mixtures of non-overlapping exponential families

Sparse models are desirable for many applications across diverse domains...

A Horseshoe Pit mixture model for Bayesian screening with an application to light sheet fluorescence microscopy in brain imaging

Finding parsimonious models through variable selection is a fundamental ...

ADAVI: Automatic Dual Amortized Variational Inference Applied To Pyramidal Bayesian Models

Frequently, population studies feature pyramidally-organized data repres...

Fast Exact Bayesian Inference for Sparse Signals in the Normal Sequence Model

We consider exact algorithms for Bayesian inference with model selection...

A Variational Inference Algorithm for BKMR in the Cross-Sectional Setting

The identification of pollutant effects is an important task in environm...