Journal
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
Volume 113, Issue 524, Pages 1431-1442Publisher
TAYLOR & FRANCIS INC
DOI: 10.1080/01621459.2017.1409121
Keywords
Bayesian methods; Computationally intensive methods; Finite precision; Gibbs sampling
Categories
Funding
- NIH [GM109452]
Ask authors/readers for more resources
Three-dimensional confocal scanning laser microscope images offer dramatic visualizations of living biofilms before and after interventions. Here, we use confocal microscopy to study the effect of a treatment over time that causes a biofilm to swell and contract due to osmotic pressure changes. From these data (the video is provided in the supplementary materials), our goal is to reconstruct biofilm surfaces, to estimate the effect of the treatment on the biofilm's volume, and to quantify the related uncertainties. We formulate the associated massive linear Bayesian inverse problem and then solve it using iterative samplers from large multivariate Gaussians that exploit well-established polynomial acceleration techniques from numerical linear algebra. Because of a general equivalence with linear solvers, these polynomial accelerated iterative samplers have known convergence rates, stopping criteria, and perform well in finite precision. An explicit algorithm is provided, for the first time, for an iterative sampler that is accelerated by the synergistic implementation of preconditioned conjugate gradient and Chebyshev polynomials. Supplementary materials for this article are available online.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available