Structure-preserving Approximate Bayesian Computation (ABC) for stochastic neuronal models
Mar 24, 2022
"Structure-preserving Approximate Bayesian Computation (ABC) for stochastic neuronal models" is the second seminar in the series Pathways to the 2023 IHP thematic project Random Processes in the Brain.
Title: Structure-preserving Approximate Bayesian Computation (ABC) for stochastic neuronal models
Date: 2022-03-29 at 10:00 GMT-3
Speaker: Massimiliano Tamborrino
Affiliation: Department of Statistics, University of Warwick
Abstract: ABC has become one of the major tools for parameter inference in complex mathematical models in the last decade. The method is based on the idea of deriving an approximate posterior density aiming to target the true (unavailable) posterior by running massive simulations from the model for different parameters to replace the intractable likelihood, choosing then those parameters whose simulations are good matches to the observed data. When applying ABC to stochastic models, the derivation of effective summary statistics and proper distances is particularly challenging, since simulations from the model under the same parameter configuration result in different output. Moreover, since exact simulation from complex stochastic models is rarely possible, reliable numerical methods need to be applied. In this talk, we show how to use the underlying structural properties of the model to construct specific ABC summaries that are less sensitive to the intrinsic stochasticity of the model, and the importance of adopting reliable property-preserving numerical (splitting) schemes for the synthetic data generation. Indeed, the commonly used Euler-Maruyama scheme may drastically fail even with very small stepsizes. The proposed approach is illustrated first on the stochastic FitzHugh-Nagumo model, and then on the broad class of partially observed Hamiltonian stochastic differential equations, in particular on the stochastic Jensen-and-Rit neural mass model, both with simulated and with real electroencephalography (EEG) data, for both one neural population [2] and a network of neural populations (ongoing work).
[1] E. Buckwar, A. Samson, M. Tamborrino, I. Tubikanec. A splitting method for SDEs with locally Lipschitz drift: Illustration on the FitzHugh-Nagumo model. ArXiv:2101.01027, https://arxiv.org/abs/2101.01027
[2] E. Buckwar, M. Tamborrino, I. Tubikanec. Spectral density-based and measure-preserving ABC for partially observed diffusion processes. An illustration on Hamiltonian SDEs. Stat. Comput. 30 (3), 627-648, 2020.
Share on Facebook
