On a framework of data assimilation for hyperparameter estimation of spiking neuronal networks

When handling real-world data modeled by a complex network dynamical system, the number of the parameters is often much more than the size of the data. Therefore, in many cases, it is impossible to estimate these parameters and the exact value of each parameter is frequently less interesting than the distribution of the parameters. In this paper, we aim to estimate the distribution of the parameters in the mesoscopic neuronal network model from the macroscopic experimental data, for example, the BOLD (blood oxygen level dependent) signal. Herein, we assume that the parameters of the neurons and synapses are inhomogeneous but independently and identically distributed from certain distributions with unknown hyperparameters. Thus, we estimate these hyperparameters of the distributions of the parameters, instead of estimating the parameters themselves. We formulate this problem under the framework of data assimilation and hierarchical Bayesian method and present an efficient method named Hierarchical Data Assimilation (HDA) to conduct the statistical inference on the neuronal network model with the BOLD signal data simulated by the hemodynamic model. We consider the Leaky Integral-Fire (LIF) neuronal networks with four synapses and show that the proposed algorithm can estimate the BOLD signals and the hyperparameters with high preciseness. In addition, we discuss the influence on the performance of the algorithm configuration and the LIF network model setup.