2021/07/29 by Fedor Goncharov, Goncharov, Fedor, Éric Barat +3
Computer Science · Mathematics · Medicine · #Applications (stat.AP) #FOS: Computer and information sciences #FOS: Physical sciences #Machine Learning (stat.ML) #Machine Learning and Algorithms #Mathematical Physics (math-ph) #Medical Imaging Techniques and Applications #Methodology (stat.ME) #Statistical Methods and Inference
paper · pdf · doi:10.48550/arxiv.2108.00866
openalex publication_date 2021/07/29 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We continue studies of the uncertainty quantification problem in emission tomographies such as PET or SPECT when additional multimodal data (e.g., anatomical MRI images) are available. To solve the aforementioned problem we adapt the recently proposed nonparametric posterior learning technique to the context of Poisson-type data in emission tomography. Using this approach we derive sampling algorithms which are trivially parallelizable, scalable and very easy to implement. In addition, we prove conditional consistency and tightness for the distribution of produced samples in the small noise limit (i.e., when the acquisition time tends to infinity) and derive new geometrical and necessary condition on how MRI images must be used. This condition arises naturally in the context of identifiability problem for misspecified generalized Poisson models. We also contrast our approach with Bayesian Markov Chain Monte Carlo sampling based on one data augmentation scheme which is very popular in the context of Expectation-Maximization algorithms for PET or SPECT. We show theoretically and also numerically that such data augmentation significantly increases mixing times for the Markov chain. In view of this, our algorithms seem to give a reasonable trade-off between design complexity, scalability, numerical load and assessment for the uncertainty.