Abstract:Abstract In diffusion MRI analysis, advances in biophysical multi-compartment modeling have gained popularity over the conventional Diffusion Tensor Imaging (DTI), because they possess greater specificity in relating the dMRI signal to underlying cellular microstructure. Biophysical multi-compartment models require parameter estimation, typically performed using either Maximum Likelihood Estimation (MLE) or using Monte Carlo Markov Chain (MCMC) sampling. Whereas MLE provides only a point estimate of the fitted model parameters, MCMC recovers the entire posterior distribution of the model parameters given the data, providing additional information such as parameter uncertainty and correlations. MCMC sampling is currently not routinely applied in dMRI microstructure modeling because it requires adjustments and tuning specific to each model, particularly in the choice of proposal distributions, burn-in length, thinning and the number of samples to store. In addition, sampling often takes at least an order of magnitude more time than non-linear optimization. Here we investigate the performance of MCMC algorithm variations over multiple popular diffusion microstructure models to see whether a single well performing variation could be applied efficiently and robustly to many models. Using an efficient GPU-based implementation, we show that run times can be removed as a prohibitive constraint for sampling of diffusion multi-compartment models. Using this implementation, we investigated the effectiveness of different adaptive MCMC algorithms, burn-in, initialization and thinning. Finally we apply the theory of Effective Sample Size to diffusion multi-compartment models as a way of determining a relatively general target for the number of samples needed to characterize parameter distributions for different models and datasets. We conclude that robust and fast sampling is achieved in most diffusion microstructure models with the Adaptive Metropolis-Within-Gibbs (AMWG) algorithm initialized with an MLE point estimate, in which case 100 to 200 samples are sufficient as a burn-in and thinning is mostly unnecessary. As a relatively general target for the number of samples, we recommend a multivariate Effective Sample Size of 2200.

Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images

Realistic voxel sizes and reduced signal variation in Monte-Carlo simulation for diffusion MR data synthesis

Portable simulation framework for diffusion MRI

A generative model of realistic brain cells with application to numerical simulation of diffusion-weighted MR signal

CMRsim–A python package for cardiovascular MR simulations incorporating complex motion and flow

Robust and fast Monte Carlo Markov Chain sampling of diffusion MRI microstructure models

Particle‐based MR modeling with diffusion, microstructure, and enzymatic reactions

Simulation of Water Diffusion Through a Simple Virtual Cardiac Cell Model: Optimization of Monte Carlo Parameters and Observation by Simulated MRI

Monte Carlo Simulation of Water Diffusion Through Cardiac Tissue Models.

Submillimeter diffusion MRI using an in-plane segmented 3D multi-slab acquisition and denoiser-regularized reconstruction

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models

High-Frequency Space Diffusion Models for Accelerated MRI

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems

Practical computation of the diffusion MRI signal of realistic neurons based on Laplace eigenfunctions

Microstructure estimation from diffusion-MRI: Compartmentalized models in permeable cellular tissue

High-Frequency Space Diffusion Model for Accelerated MRI

A diffusion MRI model for random walks confined on cylindrical surfaces: Towards non-invasive quantification of myelin sheath radius

PhoenixMR: A GPU-based MRI simulation framework with runtime-dynamic code execution

Histology-informed microstructural diffusion simulations for MRI cancer characterisation - the Histo-microSim framework

Characterization of Anomalous Diffusion in Porous Biological Tissues Using Fractional Order Derivatives and Entropy

High-resolution whole-brain diffusion MRI at 3T using simultaneous multi-slab (SMSlab) acquisition