Bayesian modelling and analysis of spatio-temporal neuronal networks

This paper illustrates a novel hierarchical dynamic Bayesian network modelling the spiking patterns of neuronal ensembles over time. We introduce, at separate model stages, the parameters characterizing the discrete-time spiking process, the unknown structure of the functional connections among the analysed neurons and its dependence on their spatial arrangement. Estimates for all model parameters and predictions for future spiking states are computed under the Bayesian paradigm via the standard Gibbs sampler using a shrinkage prior. The adequacy of the model is investigated by plotting the residuals and by applying the time-rescaling theorem. We analyse a simulated dataset and a set of experimental multiple spike trains obtained from a culture of neurons in vitro. For the latter data, we nd that one neuron plays a pivotal role for the initiation of each cycle of network activity and that the estimated network structure signicantly depends on the spatial arrangement of the neurons. © 2006 International Society for Bayesian Analysis.

[1]  S E Fienberg,et al.  Stochastic models for single neuron firing trains: a survey. , 1974, Biometrics.

[2]  Michael A. West,et al.  Deconvolution of Mixtures in Analysis of Neural Synaptic Transmission , 1994 .

[3]  Nir Friedman,et al.  Inferring Cellular Networks Using Probabilistic Graphical Models , 2004, Science.

[4]  M. Wilson,et al.  Analyzing Functional Connectivity Using a Network Likelihood Model of Ensemble Neural Spiking Activity , 2005, Neural Computation.

[5]  Wang,et al.  Nonuniversal critical dynamics in Monte Carlo simulations. , 1987, Physical review letters.

[6]  Yuzuru Takamura,et al.  Investigating neuronal activity with planar microelectrode arrays: achievements and new perspectives. , 2005, Journal of bioscience and bioengineering.

[7]  Wulfram Gerstner,et al.  Spiking Neuron Models , 2002 .

[8]  David J. Spiegelhalter,et al.  Local computations with probabilities on graphical structures and their application to expert systems , 1990 .

[9]  P. Gehler,et al.  An introduction to graphical models , 2001 .

[10]  P. Green,et al.  Bayesian Variable Selection and the Swendsen-Wang Algorithm , 2004 .

[11]  David R. Brillinger,et al.  Nerve Cell Spike Train Data Analysis: A Progression of Technique , 1992 .

[12]  R. Kass,et al.  Bayesian curve-fitting with free-knot splines , 2001 .

[13]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

[14]  Robert E. Kass,et al.  A Spike-Train Probability Model , 2001, Neural Computation.

[15]  David R. Brillinger,et al.  Assessing Connections in Networks of Biological Neurons , 2007 .

[16]  Simon J. Godsill,et al.  On the relationship between MCMC model uncertainty methods , 1997 .

[17]  R. Kass,et al.  Multiple neural spike train data analysis: state-of-the-art and future challenges , 2004, Nature Neuroscience.

[18]  Mark W. Woolrich,et al.  Fully Bayesian spatio-temporal modeling of FMRI data , 2004, IEEE Transactions on Medical Imaging.

[19]  Kathryn B. Laskey,et al.  Neural Coding: Higher-Order Temporal Patterns in the Neurostatistics of Cell Assemblies , 2000, Neural Computation.

[20]  Eugene M. Izhikevich,et al.  Resonate-and-fire neurons , 2001, Neural Networks.

[21]  D. Pregibon,et al.  Graphical Methods for Assessing Logistic Regression Models , 1984 .

[22]  L. Tierney.  Markov Chains for Exploring Posterior Distributions , 1994 .

[23]  R. Kass,et al.  Statistical smoothing of neuronal data. , 2003, Network.

[24]  Marco Gori,et al.  Adaptive Processing of Sequences and Data Structures , 1998, Lecture Notes in Computer Science.

[25]  David R. Brillinger,et al.  Empirical examination of the threshold model of neuron firing , 1979, Biological Cybernetics.

[26]  Kevin Murphy,et al.  Modelling Gene Expression Data using Dynamic Bayesian Networks , 2006 .

[27]  Satish Iyengar,et al.  The Analysis of Multiple Neural Spike Trains , 2003, Advances on Methodological and Applied Aspects of Probability and Statistics.

[28]  D. Brillinger.  Maximum likelihood analysis of spike trains of interacting nerve cells , 2004, Biological Cybernetics.

[29]  Emery N. Brown,et al.  The Time-Rescaling Theorem and Its Application to Neural Spike Train Data Analysis , 2002, Neural Computation.

[30]  A. Gelfand,et al.  Bayesian Model Choice: Asymptotics and Exact Calculations , 1994 .

[31]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[32]  E. George,et al.  Journal of the American Statistical Association is currently published by American Statistical Association. , 2007 .

[33]  M. Escobar,et al.  Bayesian Density Estimation and Inference Using Mixtures , 1995 .

[34]  P. Dellaportas,et al.  Bayesian variable selection using the Gibbs sampler , 2000 .

[35]  Robert E Kass,et al.  Statistical issues in the analysis of neuronal data. , 2005, Journal of neurophysiology.

[36]  Wim L. C. Rutten,et al.  Long-term characterization of firing dynamics of spontaneous bursts in cultured neural networks , 2004, IEEE Transactions on Biomedical Engineering.

[37]  Kevin P. Murphy,et al.  Learning the Structure of Dynamic Probabilistic Networks , 1998, UAI.

[38]  C. Geyer.  Markov Chain Monte Carlo Maximum Likelihood , 1991 .

[39]  G. C. Tiao,et al.  A bayesian approach to some outlier problems. , 1968, Biometrika.

[40]  Michael I. Jordan.  Graphical Models , 2003 .

[41]  Dean Phillips Foster,et al.  Calibration and Empirical Bayes Variable Selection , 1997 .

[42]  E. George,et al.  APPROACHES FOR BAYESIAN VARIABLE SELECTION , 1997 .

[43]  T. Ferguson.  BAYESIAN DENSITY ESTIMATION BY MIXTURES OF NORMAL DISTRIBUTIONS , 1983 .

[44]  H. L. Bryant,et al.  Identification of synaptic interactions , 1976, Biological Cybernetics.

[45]  Rajesh P. N. Rao.  Hierarchical Bayesian Inference in Networks of Spiking Neurons , 2004, NIPS.

[46]  E. S. Chornoboy,et al.  Maximum likelihood identification of neural point process systems , 1988, Biological Cybernetics.

[47]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[48]  Michael A. West,et al.  Hierarchical Mixture Models in Neurological Transmission Analysis , 1997 .

[49]  David Draper,et al.  Assessment and Propagation of Model Uncertainty , 2011 .

[50]  D. Brillinger.  Some Statistical Methods for Random Process Data from Seismology and Neurophysiology , 1988 .

[51]  Hani Doss,et al.  On estimating the dependence between two point processes , 1989 .

[52]  D. Pregibon.  Logistic Regression Diagnostics , 1981 .

[53]  S. Chib,et al.  Bayesian residual analysis for binary response regression models , 1995 .

[54]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[55]  M. Corner,et al.  Dynamics and plasticity in developing neuronal networks in vitro. , 2005, Progress in brain research.

[56]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[57]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[58]  Peter Spirtes,et al.  Directed Cyclic Graphical Representations of Feedback Models , 1995, UAI.

[59]  F. Papangelou.  Integrability of expected increments of point processes and a related random change of scale , 1972 .

[60]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.