Abstract:Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics. A central problem in computational neural modeling is estimating model parameters that can account for observed activity patterns. While several techniques exist to perform parameter inference in special classes of abstract neural models, there are comparatively few approaches for large scale biophysically detailed neural models. In this work, we describe challenges and solutions in applying a deep learning based statistical framework to estimate parameters in a biophysically detailed large scale neural model, and emphasize the particular difficulties in estimating parameters for time series data. Our example uses a multi-scale model designed to connect human MEG/EEG recordings to the underlying cell and circuit level generators. Our approach allows for crucial insight into how cell-level properties interact to produce measured neural activity, and provides guidelines for diagnosing the quality of the estimate and uniqueness of predictions for different MEG/EEG biomarkers.

Neural Posterior Estimation of hierarchical models in neuroscience

HNPE: Leveraging global parameters for neural posterior estimation

Leveraging Global Parameters for Flow-based Neural Posterior Estimation

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Hierarchical Neural Simulation-Based Inference Over Event Ensembles

Joint Learning of Full-structure Noise in Hierarchical Bayesian Regression Models

Neural Sampling in Hierarchical Exponential-family Energy-based Models

Posterior Estimation Using Deep Learning: A Simulation Study of Compartmental Modeling in Dynamic PET

An amortized approach to non-linear mixed-effects modeling based on neural posterior estimation

Hierarchical models for neural population dynamics in the presence of non-stationarity

A Self-Organizing State-Space-Model Approach for Parameter Estimation in Hodgkin-Huxley-Type Models of Single Neurons

Network inference in a stochastic multi-population neural mass model via approximate Bayesian computation

Posterior estimation using deep learning: a simulation study of compartmental modeling in dynamic positron emission tomography

Deep Learning-Based Parameter Estimation for Neurophysiological Models of Neuroimaging Data

Bayesian Algorithms for Joint Estimation of Brain Activity and Noise in Electromagnetic Imaging

Bayesian Joint Modeling of Multiple Brain Functional Networks

Hierarchical Bayesian Inference and Learning in Spiking Neural Networks.

Hierarchical Bayesian modeling of multi-region brain cell count data

Bayesian Inference for Brain Activity from Functional Magnetic Resonance Imaging Collected at Two Spatial Resolutions

Preconditioned Neural Posterior Estimation for Likelihood-free Inference

Implementation of Bayesian Inference in Distributed Neural Networks