1
|
Weng G, Clark K, Akbarian A, Noudoost B, Nategh N. Time-varying generalized linear models: characterizing and decoding neuronal dynamics in higher visual areas. Front Comput Neurosci 2024; 18:1273053. [PMID: 38348287 PMCID: PMC10859875 DOI: 10.3389/fncom.2024.1273053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2023] [Accepted: 01/09/2024] [Indexed: 02/15/2024] Open
Abstract
To create a behaviorally relevant representation of the visual world, neurons in higher visual areas exhibit dynamic response changes to account for the time-varying interactions between external (e.g., visual input) and internal (e.g., reward value) factors. The resulting high-dimensional representational space poses challenges for precisely quantifying individual factors' contributions to the representation and readout of sensory information during a behavior. The widely used point process generalized linear model (GLM) approach provides a powerful framework for a quantitative description of neuronal processing as a function of various sensory and non-sensory inputs (encoding) as well as linking particular response components to particular behaviors (decoding), at the level of single trials and individual neurons. However, most existing variations of GLMs assume the neural systems to be time-invariant, making them inadequate for modeling nonstationary characteristics of neuronal sensitivity in higher visual areas. In this review, we summarize some of the existing GLM variations, with a focus on time-varying extensions. We highlight their applications to understanding neural representations in higher visual areas and decoding transient neuronal sensitivity as well as linking physiology to behavior through manipulation of model components. This time-varying class of statistical models provide valuable insights into the neural basis of various visual behaviors in higher visual areas and hold significant potential for uncovering the fundamental computational principles that govern neuronal processing underlying various behaviors in different regions of the brain.
Collapse
Affiliation(s)
- Geyu Weng
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, United States
- Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, UT, United States
| | - Kelsey Clark
- Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, UT, United States
| | - Amir Akbarian
- Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, UT, United States
| | - Behrad Noudoost
- Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, UT, United States
| | - Neda Nategh
- Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, UT, United States
- Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, UT, United States
| |
Collapse
|
2
|
Kass RE, Bong H, Olarinre M, Xin Q, Urban KN. Identification of interacting neural populations: methods and statistical considerations. J Neurophysiol 2023; 130:475-496. [PMID: 37465897 PMCID: PMC10642974 DOI: 10.1152/jn.00131.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 07/17/2023] [Accepted: 07/17/2023] [Indexed: 07/20/2023] Open
Abstract
As improved recording technologies have created new opportunities for neurophysiological investigation, emphasis has shifted from individual neurons to multiple populations that form circuits, and it has become important to provide evidence of cross-population coordinated activity. We review various methods for doing so, placing them in six major categories while avoiding technical descriptions and instead focusing on high-level motivations and concerns. Our aim is to indicate what the methods can achieve and the circumstances under which they are likely to succeed. Toward this end, we include a discussion of four cross-cutting issues: the definition of neural populations, trial-to-trial variability and Poisson-like noise, time-varying dynamics, and causality.
Collapse
Affiliation(s)
- Robert E Kass
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
- Department of Statistics & Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
| | - Heejong Bong
- Department of Statistics & Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
| | - Motolani Olarinre
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
- Department of Statistics & Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
| | - Qi Xin
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
- Department of Statistics & Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
| | - Konrad N Urban
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
- Department of Statistics & Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States
| |
Collapse
|
3
|
In vivo extracellular recordings of thalamic and cortical visual responses reveal V1 connectivity rules. Proc Natl Acad Sci U S A 2022; 119:e2207032119. [PMID: 36191204 PMCID: PMC9564935 DOI: 10.1073/pnas.2207032119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
The brain's connectome provides the scaffold for canonical neural computations. However, a comparison of connectivity studies in the mouse primary visual cortex (V1) reveals that the average number and strength of connections between specific neuron types can vary. Can variability in V1 connectivity measurements coexist with canonical neural computations? We developed a theory-driven approach to deduce V1 network connectivity from visual responses in mouse V1 and visual thalamus (dLGN). Our method revealed that the same recorded visual responses were captured by multiple connectivity configurations. Remarkably, the magnitude and selectivity of connectivity weights followed a specific order across most of the inferred connectivity configurations. We argue that this order stems from the specific shapes of the recorded contrast response functions and contrast invariance of orientation tuning. Remarkably, despite variability across connectivity studies, connectivity weights computed from individual published connectivity reports followed the order we identified with our method, suggesting that the relations between the weights, rather than their magnitudes, represent a connectivity motif supporting canonical V1 computations.
Collapse
|
4
|
Hu Q, Zheng Z, Sui X, Li L, Chai X, Chen Y. Spatial Attention Modulates Spike Count Correlations and Granger Causality in the Primary Visual Cortex. Front Cell Neurosci 2022; 16:838049. [PMID: 35783091 PMCID: PMC9246483 DOI: 10.3389/fncel.2022.838049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 05/23/2022] [Indexed: 11/16/2022] Open
Abstract
The influence of spatial attention on neural interactions has been revealed even in early visual information processing stages. It resolves the process of competing for sensory information about objects perceived as targets and distractors. However, the attentional modulation of the interaction between pairs of neurons with non-overlapping receptive fields (RFs) is not well known. Here, we investigated the activity of anatomically distant neurons in two behaving monkeys’ primary visual cortex (V1), when they performed a spatial attention task detecting color change. We compared attentional modulation from the perspective of spike count correlations and Granger causality among simple and complex cells. An attention-related increase in spike count correlations and a decrease in Granger causality were found. The results showed that spatial attention significantly influenced only the interactions between rather than within simple and complex cells. Furthermore, we found that the attentional modulation of neuronal interactions changed with neuronal pairs’ preferred directions differences. Thus, we found that spatial attention increased the functional communications and competing connectivities when attending to the neurons’ RFs, which impacts the interactions only between simple and complex cells. Our findings enrich the model of simple and complex cells and further understand the way that attention influences the neurons’ activities.
Collapse
Affiliation(s)
- Qiyi Hu
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Zhiyan Zheng
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Xiaohong Sui
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Liming Li
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Xinyu Chai
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Yao Chen
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
5
|
She X, Berger TW, Song D. A Double-Layer Multi-Resolution Classification Model for Decoding Spatiotemporal Patterns of Spikes With Small Sample Size. Neural Comput 2021; 34:219-254. [PMID: 34758485 DOI: 10.1162/neco_a_01459] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 08/19/2021] [Indexed: 11/04/2022]
Abstract
We build a double-layer, multiple temporal-resolution classification model for decoding single-trial spatiotemporal patterns of spikes. The model takes spiking activities as input signals and binary behavioral or cognitive variables as output signals and represents the input-output mapping with a double-layer ensemble classifier. In the first layer, to solve the underdetermined problem caused by the small sample size and the very high dimensionality of input signals, B-spline functional expansion and L1-regularized logistic classifiers are used to reduce dimensionality and yield sparse model estimations. A wide range of temporal resolutions of neural features is included by using a large number of classifiers with different numbers of B-spline knots. Each classifier serves as a base learner to classify spatiotemporal patterns into the probability of the output label with a single temporal resolution. A bootstrap aggregating strategy is used to reduce the estimation variances of these classifiers. In the second layer, another L1-regularized logistic classifier takes outputs of first-layer classifiers as inputs to generate the final output predictions. This classifier serves as a meta-learner that fuses multiple temporal resolutions to classify spatiotemporal patterns of spikes into binary output labels. We test this decoding model with both synthetic and experimental data recorded from rats and human subjects performing memory-dependent behavioral tasks. Results show that this method can effectively avoid overfitting and yield accurate prediction of output labels with small sample size. The double-layer, multi-resolution classifier consistently outperforms the best single-layer, single-resolution classifier by extracting and utilizing multi-resolution spatiotemporal features of spike patterns in the classification.
Collapse
Affiliation(s)
- Xiwei She
- Neuroscience Graduate Program, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Theodore W Berger
- Department of Biomedical Engineering and Neuroscience Graduate Program, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Dong Song
- Department of Biomedical Engineering and Neuroscience Graduate Program, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| |
Collapse
|
6
|
Keeley SL, Zoltowski DM, Aoi MC, Pillow JW. Modeling statistical dependencies in multi-region spike train data. Curr Opin Neurobiol 2020; 65:194-202. [PMID: 33334641 PMCID: PMC7769979 DOI: 10.1016/j.conb.2020.11.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Revised: 11/10/2020] [Accepted: 11/10/2020] [Indexed: 11/17/2022]
Abstract
Neural computations underlying cognition and behavior rely on the coordination of neural activity across multiple brain areas. Understanding how brain areas interact to process information or generate behavior is thus a central question in neuroscience. Here we provide an overview of statistical approaches for characterizing statistical dependencies in multi-region spike train recordings. We focus on two classes of models in particular: regression-based models and shared latent variable models. Regression-based models describe interactions in terms of a directed transformation of information from one region to another. Shared latent variable models, on the other hand, seek to describe interactions in terms of sources that capture common fluctuations in spiking activity across regions. We discuss the advantages and limitations of each of these approaches and future directions for the field. We intend this review to be an introduction to the statistical methods in multi-region models for computational neuroscientists and experimentalists alike.
Collapse
Affiliation(s)
- Stephen L Keeley
- Princeton Neuroscience Institute, Princeton University, Princeton, NJ, USA
| | - David M Zoltowski
- Princeton Neuroscience Institute, Princeton University, Princeton, NJ, USA
| | - Mikio C Aoi
- Princeton Neuroscience Institute, Princeton University, Princeton, NJ, USA
| | - Jonathan W Pillow
- Princeton Neuroscience Institute, Princeton University, Princeton, NJ, USA.
| |
Collapse
|
7
|
Tauste Campo A. Inferring neural information flow from spiking data. Comput Struct Biotechnol J 2020; 18:2699-2708. [PMID: 33101608 PMCID: PMC7548302 DOI: 10.1016/j.csbj.2020.09.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 09/05/2020] [Accepted: 09/07/2020] [Indexed: 01/02/2023] Open
Abstract
The brain can be regarded as an information processing system in which neurons store and propagate information about external stimuli and internal processes. Therefore, estimating interactions between neural activity at the cellular scale has significant implications in understanding how neuronal circuits encode and communicate information across brain areas to generate behavior. While the number of simultaneously recorded neurons is growing exponentially, current methods relying only on pairwise statistical dependencies still suffer from a number of conceptual and technical challenges that preclude experimental breakthroughs describing neural information flows. In this review, we examine the evolution of the field over the years, starting from descriptive statistics to model-based and model-free approaches. Then, we discuss in detail the Granger Causality framework, which includes many popular state-of-the-art methods and we highlight some of its limitations from a conceptual and practical estimation perspective. Finally, we discuss directions for future research, including the development of theoretical information flow models and the use of dimensionality reduction techniques to extract relevant interactions from large-scale recording datasets.
Collapse
Affiliation(s)
- Adrià Tauste Campo
- Centre for Brain and Cognition, Universitat Pompeu Fabra, Ramon Trias Fargas 25, 08018 Barcelona, Spain
| |
Collapse
|
8
|
Qian C, Sun X, Wang Y, Zheng X, Wang Y, Pan G. Binless Kernel Machine: Modeling Spike Train Transformation for Cognitive Neural Prostheses. Neural Comput 2020; 32:1863-1900. [PMID: 32795229 DOI: 10.1162/neco_a_01306] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
Modeling spike train transformation among brain regions helps in designing a cognitive neural prosthesis that restores lost cognitive functions. Various methods analyze the nonlinear dynamic spike train transformation between two cortical areas with low computational eficiency. The application of a real-time neural prosthesis requires computational eficiency, performance stability, and better interpretation of the neural firing patterns that modulate target spike generation. We propose the binless kernel machine in the point-process framework to describe nonlinear dynamic spike train transformations. Our approach embeds the binless kernel to eficiently capture the feedforward dynamics of spike trains and maps the input spike timings into reproducing kernel Hilbert space (RKHS). An inhomogeneous Bernoulli process is designed to combine with a kernel logistic regression that operates on the binless kernel to generate an output spike train as a point process. Weights of the proposed model are estimated by maximizing the log likelihood of output spike trains in RKHS, which allows a global-optimal solution. To reduce computational complexity, we design a streaming-based clustering algorithm to extract typical and important spike train features. The cluster centers and their weights enable the visualization of the important input spike train patterns that motivate or inhibit output neuron firing. We test the proposed model on both synthetic data and real spike train data recorded from the dorsal premotor cortex and the primary motor cortex of a monkey performing a center-out task. Performances are evaluated by discrete-time rescaling Kolmogorov-Smirnov tests. Our model outperforms the existing methods with higher stability regardless of weight initialization and demonstrates higher eficiency in analyzing neural patterns from spike timing with less historical input (50%). Meanwhile, the typical spike train patterns selected according to weights are validated to encode output spike from the spike train of single-input neuron and the interaction of two input neurons.
Collapse
Affiliation(s)
- Cunle Qian
- College of Computer Science, Zhejiang University, Hangzhou 310027, P.R.C., and Department of Electronic and Computer Engineering, Hong Kong University of Science and Technology, Kowloon, Hong Kong SAR 99077, P.R.C.
| | - Xuyun Sun
- College of Computer Science, Zhejiang University, Hangzhou 310027, P.R.C.
| | - Yueming Wang
- Qiushi Academy for Advanced Studies, Zhejiang University, Hangzhou 310027, P.R.C.
| | - Xiaoxiang Zheng
- Qiushi Academy for Advanced Studies, Zhejiang University, Hangzhou 310027, P.R.C.
| | - Yiwen Wang
- Department of Electronic and Computer Engineering and Department of Chemical and Biological Engineering, Hong Kong University of Science and Technology, Kowloon, Hong Kong SAR 99077, P.R.C.
| | - Gang Pan
- College of Computer Science and State Key Laboratory of CAD&CG, Zhejiang University, Hangzhou 310027, P.R.C.
| |
Collapse
|
9
|
Inference of synaptic connectivity and external variability in neural microcircuits. J Comput Neurosci 2020; 48:123-147. [PMID: 32080777 DOI: 10.1007/s10827-020-00739-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2019] [Revised: 11/15/2019] [Accepted: 01/03/2020] [Indexed: 10/25/2022]
Abstract
A major goal in neuroscience is to estimate neural connectivity from large scale extracellular recordings of neural activity in vivo. This is challenging in part because any such activity is modulated by the unmeasured external synaptic input to the network, known as the common input problem. Many different measures of functional connectivity have been proposed in the literature, but their direct relationship to synaptic connectivity is often assumed or ignored. For in vivo data, measurements of this relationship would require a knowledge of ground truth connectivity, which is nearly always unavailable. Instead, many studies use in silico simulations as benchmarks for investigation, but such approaches necessarily rely upon a variety of simplifying assumptions about the simulated network and can depend on numerous simulation parameters. We combine neuronal network simulations, mathematical analysis, and calcium imaging data to address the question of when and how functional connectivity, synaptic connectivity, and latent external input variability can be untangled. We show numerically and analytically that, even though the precision matrix of recorded spiking activity does not uniquely determine synaptic connectivity, it is in practice often closely related to synaptic connectivity. This relation becomes more pronounced when the spatial structure of neuronal variability is jointly considered.
Collapse
|
10
|
Qian C, Sun X, Yang Z, Pan G, Wang Y. A K-Medoids based Point-Process Modeling on Neural Spike Transformation using Binless Kernel. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020; 2019:4387-4390. [PMID: 31946839 DOI: 10.1109/embc.2019.8856479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
A neural prosthesis is designed to compensate for cognitive functional losses by modeling the information transmission among cortical areas. Existing methods generally build a generalized linear model to approximate the nonlinear transformation among two areas, and use the temporal information of the neural spike with low efficiency. It is essential to efficiently model the nonlinearity embedded in spike generation and transmission for the real-time. This paper proposes a nonlinear point-process model to describe spike-in and spike-out transformation using the theory of reproducing kernel Hilbert space (RKHS) and the binless kernel on spike trains. The binless kernel efficiently maps exact spike timing information to the RKHS to describe nonlinear transformations with global minimum regardless of the weight initialization. A streaming K-medoids algorithm is introduced to select typical and important features in this nonlinear binless kernel for further modeling. We test our model on the nonlinearly generated synthetic neural spike trains, and compare with the existing spike transformation methods, such as Volterra model and staged point-process model. The results show that our model has higher goodness-of-fit evaluated by Kolmogorov-Smirnov test and less variance on the prediction results, which indicates the potential better modeling approach for neural prosthesis application.
Collapse
|
11
|
Kobayashi R, Kurita S, Kurth A, Kitano K, Mizuseki K, Diesmann M, Richmond BJ, Shinomoto S. Reconstructing neuronal circuitry from parallel spike trains. Nat Commun 2019; 10:4468. [PMID: 31578320 PMCID: PMC6775109 DOI: 10.1038/s41467-019-12225-2] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 08/27/2019] [Indexed: 11/23/2022] Open
Abstract
State-of-the-art techniques allow researchers to record large numbers of spike trains in parallel for many hours. With enough such data, we should be able to infer the connectivity among neurons. Here we develop a method for reconstructing neuronal circuitry by applying a generalized linear model (GLM) to spike cross-correlations. Our method estimates connections between neurons in units of postsynaptic potentials and the amount of spike recordings needed to verify connections. The performance of inference is optimized by counting the estimation errors using synthetic data. This method is superior to other established methods in correctly estimating connectivity. By applying our method to rat hippocampal data, we show that the types of estimated connections match the results inferred from other physiological cues. Thus our method provides the means to build a circuit diagram from recorded spike trains, thereby providing a basis for elucidating the differences in information processing in different brain regions. Current techniques have enabled the simultaneous collection of spike train data from large numbers of neurons. Here, the authors report a method to infer the underlying neural circuit connectivity diagram based on a generalized linear model applied to spike cross-correlations between neurons.
Collapse
Affiliation(s)
- Ryota Kobayashi
- National Institute of Informatics, Tokyo, 101-8430, Japan.,Department of Informatics, SOKENDAI (The Graduate University for Advanced Studies), Tokyo, 101-8430, Japan
| | - Shuhei Kurita
- Center for Advanced Intelligence Project, RIKEN, Tokyo, 103-0027, Japan
| | - Anno Kurth
- Institute of Neuroscience and Medicine (INM-6) and Institute for Advanced Simulation (IAS-6) and JARA-Institute Brain Structure-Function Relationships (INM-10), Jülich Research Centre, 52425, Jülich, Germany.,Department of Physics, Faculty 1, RWTH Aachen University, Aachen, Germany
| | - Katsunori Kitano
- Department of Information Science and Engineering, Ritsumeikan University, Kusatsu, 525-8577, Japan
| | - Kenji Mizuseki
- Department of Physiology, Osaka City University Graduate School of Medicine, Osaka, 545-8585, Japan
| | - Markus Diesmann
- Institute of Neuroscience and Medicine (INM-6) and Institute for Advanced Simulation (IAS-6) and JARA-Institute Brain Structure-Function Relationships (INM-10), Jülich Research Centre, 52425, Jülich, Germany.,Department of Physics, Faculty 1, RWTH Aachen University, Aachen, Germany.,Department of Psychiatry, Psychotherapy and Psychosomatics, School of Medicine, RWTH Aachen University, Aachen, Germany
| | - Barry J Richmond
- Laboratory of Neuropsychology, NIMH/NIH/DHHS, Bethesda, MD, 20814, USA
| | - Shigeru Shinomoto
- Department of Physics, Kyoto University, Kyoto, 606-8502, Japan. .,Brain Information Communication Research Laboratory Group, ATR Institute International, Kyoto, 619-0288, Japan.
| |
Collapse
|
12
|
Verzelli P, Sacerdote L. A study of dependency features of spike trains through copulas. Biosystems 2019; 184:104014. [PMID: 31401080 DOI: 10.1016/j.biosystems.2019.104014] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Revised: 06/15/2019] [Accepted: 08/05/2019] [Indexed: 11/27/2022]
Abstract
Despite the progresses of statistical and machine learning techniques, simultaneous recordings from many neurons hide important information and the connections characterizing the network remain generally undiscovered. Discerning the presence of direct links between neurons from data is still a not completely solved problem. We propose the use of copulas, to enlarge the number of tools for detecting the network structure, pursuing on a research direction we started in Sacerdote et al. (2012). Here, our aim is to distinguish different types of connections on a very simple network. Our proposal consists in choosing suitable random intervals in pairs of spike trains determining the shapes of their copulas. We show that this approach allows to detect different types of dependencies. We illustrate the features of the proposed method on synthetic data from suitably connected networks of two or three formal neurons directly connected or influenced by the surrounding network. We show how a smart choice of pairs of random times together with the use of empirical copulas allows to discern between direct and indirect interactions.
Collapse
|
13
|
Liu H, Bridges D, Randall C, Solla SA, Wu B, Hansma P, Yan X, Kosik KS, Bouchard K. In vitro validation of in silico identified inhibitory interactions. J Neurosci Methods 2019; 321:39-48. [PMID: 30965073 DOI: 10.1016/j.jneumeth.2019.04.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Accepted: 04/01/2019] [Indexed: 10/27/2022]
Abstract
BACKGROUND Understanding how neuronal signals propagate in local network is an important step in understanding information processing. As a result, spike trains recorded with multi-electrode arrays (MEAs) have been widely used to study the function of neural networks. Studying the dynamics of neuronal networks requires the identification of both excitatory and inhibitory connections. The detection of excitatory relationships can robustly be inferred by characterizing the statistical relationships of neural spike trains. However, the identification of inhibitory relationships is more difficult: distinguishing endogenous low firing rates from active inhibition is not obvious. NEW METHOD In this paper, we propose an in silico interventional procedure that makes predictions about the effect of stimulating or inhibiting single neurons on other neurons, and thereby gives the ability to accurately identify inhibitory effects. COMPARISON To experimentally test these predictions, we have developed a Neural Circuit Probe (NCP) that delivers drugs transiently and reversibly on individually identified neurons to assess their contributions to the neural circuit behavior. RESULTS Using the NCP, putative inhibitory connections identified by the in silico procedure were validated through in vitro interventional experiments. CONCLUSIONS Together, these results demonstrate how detailed microcircuitry can be inferred from statistical models derived from neurophysiology data.
Collapse
Affiliation(s)
- Honglei Liu
- Department of Computer Science, University of California, Santa Barbara, CA, USA
| | - Daniel Bridges
- Department of Physics, University of California, Santa Barbara, CA, USA
| | - Connor Randall
- Department of Physics, University of California, Santa Barbara, CA, USA
| | - Sara A Solla
- Department of Physiology, Northwestern University, Chicago, IL, USA
| | - Bian Wu
- Neuroscience Research Institute, University of California, Santa Barbara, CA, USA; Department of Molecular Cellular and Developmental Biology, University of California, Santa Barbara, CA, USA
| | - Paul Hansma
- Department of Physics, University of California, Santa Barbara, CA, USA
| | - Xifeng Yan
- Department of Computer Science, University of California, Santa Barbara, CA, USA
| | - Kenneth S Kosik
- Neuroscience Research Institute, University of California, Santa Barbara, CA, USA.
| | - Kristofer Bouchard
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA; Helen Wills Neuroscience Institute and Redwood Center for Theoretical Neuroscience, University of California, Berkeley, CA, USA.
| |
Collapse
|
14
|
Geng K, Shin DC, Song D, Hampson RE, Deadwyler SA, Berger TW, Marmarelis VZ. Multi-Input, Multi-Output Neuronal Mode Network Approach to Modeling the Encoding Dynamics and Functional Connectivity of Neural Systems. Neural Comput 2019; 31:1327-1355. [PMID: 31113305 DOI: 10.1162/neco_a_01204] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
This letter proposes a novel method, multi-input, multi-output neuronal mode network (MIMO-NMN), for modeling encoding dynamics and functional connectivity in neural ensembles such as the hippocampus. Compared with conventional approaches such as the Volterra-Wiener model, linear-nonlinear-cascade (LNC) model, and generalized linear model (GLM), the NMN has several advantages in terms of estimation accuracy, model interpretation, and functional connectivity analysis. We point out the limitations of current neural spike modeling methods, especially the estimation biases caused by the imbalanced class problem when the number of zeros is significantly larger than ones in the spike data. We use synthetic data to test the performance of NMN with a comparison of the traditional methods, and the results indicate the NMN approach could reduce the imbalanced class problem and achieve better predictions. Subsequently, we apply the MIMO-NMN method to analyze data from the human hippocampus. The results indicate that the MIMO-NMN method is a promising approach to modeling neural dynamics and analyzing functional connectivity of multi-neuronal data.
Collapse
Affiliation(s)
- Kunling Geng
- Department of Biomedical Engineering and Biomedical Simulations Resource Center, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Dae C Shin
- Department of Biomedical Engineering and Biomedical Simulations Resource Center, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Dong Song
- Department of Biomedical Engineering and Biomedical Simulations Resource Center, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Robert E Hampson
- Department of Physiology and Pharmacology, Wake Forest School of Medicine, Winston-Salem, NC, 27157, U.S.A.
| | - Samuel A Deadwyler
- Department of Physiology and Pharmacology, Wake Forest School of Medicine, Winston-Salem, NC, 27157, U.S.A.
| | - Theodore W Berger
- Department of Biomedical Engineering and Biomedical Simulations Resource Center, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| | - Vasilis Z Marmarelis
- Department of Biomedical Engineering and Biomedical Simulations Resource Center, University of Southern California, Los Angeles, CA, 90089, U.S.A.
| |
Collapse
|
15
|
Zoltowski DM, Pillow JW. Scaling the Poisson GLM to massive neural datasets through polynomial approximations. ADVANCES IN NEURAL INFORMATION PROCESSING SYSTEMS 2018; 31:3517-3527. [PMID: 31244513 PMCID: PMC6594562] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Recent advances in recording technologies have allowed neuroscientists to record simultaneous spiking activity from hundreds to thousands of neurons in multiple brain regions. Such large-scale recordings pose a major challenge to existing statistical methods for neural data analysis. Here we develop highly scalable approximate inference methods for Poisson generalized linear models (GLMs) that require only a single pass over the data. Our approach relies on a recently proposed method for obtaining approximate sufficient statistics for GLMs using polynomial approximations [7], which we adapt to the Poisson GLM setting. We focus on inference using quadratic approximations to nonlinear terms in the Poisson GLM log-likelihood with Gaussian priors, for which we derive closed-form solutions to the approximate maximum likelihood and MAP estimates, posterior distribution, and marginal likelihood. We introduce an adaptive procedure to select the polynomial approximation interval and show that the resulting method allows for efficient and accurate inference and regularization of high-dimensional parameters. We use the quadratic estimator to fit a fully-coupled Poisson GLM to spike train data recorded from 831 neurons across five regions of the mouse brain for a duration of 41 minutes, binned at 1 ms resolution. Across all neurons, this model is fit to over 2 billion spike count bins and identifies fine-timescale statistical dependencies between neurons within and across cortical and subcortical areas.
Collapse
Affiliation(s)
- David M Zoltowski
- Princeton Neuroscience Institute, Princeton University; Princeton, NJ 08544
| | - Jonathan W Pillow
- Princeton Neuroscience Institute & Psychology, Princeton University; Princeton, NJ 08544
| |
Collapse
|
16
|
Capone C, Gigante G, Del Giudice P. Spontaneous activity emerging from an inferred network model captures complex spatio-temporal dynamics of spike data. Sci Rep 2018; 8:17056. [PMID: 30451957 PMCID: PMC6242821 DOI: 10.1038/s41598-018-35433-0] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2018] [Accepted: 10/28/2018] [Indexed: 11/09/2022] Open
Abstract
Inference methods are widely used to recover effective models from observed data. However, few studies attempted to investigate the dynamics of inferred models in neuroscience, and none, to our knowledge, at the network level. We introduce a principled modification of a widely used generalized linear model (GLM), and learn its structural and dynamic parameters from in-vitro spike data. The spontaneous activity of the new model captures prominent features of the non-stationary and non-linear dynamics displayed by the biological network, where the reference GLM largely fails, and also reflects fine-grained spatio-temporal dynamical features. Two ingredients were key for success. The first is a saturating transfer function: beyond its biological plausibility, it limits the neuron's information transfer, improving robustness against endogenous and external noise. The second is a super-Poisson spikes generative mechanism; it accounts for the undersampling of the network, and allows the model neuron to flexibly incorporate the observed activity fluctuations.
Collapse
Affiliation(s)
- Cristiano Capone
- Physics department, "Sapienza" University, Rome, Italy
- INFN, Sezione di Roma, Rome, Italy
| | | | | |
Collapse
|
17
|
Qian C, Sun X, Zhang S, Xing D, Li H, Zheng X, Pan G, Wang Y. Nonlinear Modeling of Neural Interaction for Spike Prediction Using the Staged Point-Process Model. Neural Comput 2018; 30:3189-3226. [PMID: 30314427 DOI: 10.1162/neco_a_01137] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Neurons communicate nonlinearly through spike activities. Generalized linear models (GLMs) describe spike activities with a cascade of a linear combination across inputs, a static nonlinear function, and an inhomogeneous Bernoulli or Poisson process, or Cox process if a self-history term is considered. This structure considers the output nonlinearity in spike generation but excludes the nonlinear interaction among input neurons. Recent studies extend GLMs by modeling the interaction among input neurons with a quadratic function, which considers the interaction between every pair of input spikes. However, quadratic effects may not fully capture the nonlinear nature of input interaction. We therefore propose a staged point-process model to describe the nonlinear interaction among inputs using a few hidden units, which follows the idea of artificial neural networks. The output firing probability conditioned on inputs is formed as a cascade of two linear-nonlinear (a linear combination plus a static nonlinear function) stages and an inhomogeneous Bernoulli process. Parameters of this model are estimated by maximizing the log likelihood on output spike trains. Unlike the iterative reweighted least squares algorithm used in GLMs, where the performance is guaranteed by the concave condition, we propose a modified Levenberg-Marquardt (L-M) algorithm, which directly calculates the Hessian matrix of the log likelihood, for the nonlinear optimization in our model. The proposed model is tested on both synthetic data and real spike train data recorded from the dorsal premotor cortex and primary motor cortex of a monkey performing a center-out task. Performances are evaluated by discrete-time rescaled Kolmogorov-Smirnov tests, where our model statistically outperforms a GLM and its quadratic extension, with a higher goodness-of-fit in the prediction results. In addition, the staged point-process model describes nonlinear interaction among input neurons with fewer parameters than quadratic models, and the modified L-M algorithm also demonstrates fast convergence.
Collapse
Affiliation(s)
- Cunle Qian
- College of Computer Science, Zhejiang University, Hangzhou, 310027, China
| | - Xuyun Sun
- College of Computer Science, Zhejiang University, Hangzhou, 310027, China
| | - Shaomin Zhang
- Qiushi Academy for Advanced Studies, Zhejiang University, Hangzhou, 310027, China
| | - Dong Xing
- College of Computer Science, Zhejiang University, Hangzhou, 310027, China
| | - Hongbao Li
- Qiushi Academy for Advanced Studies, Zhejiang University, Hangzhou, 310027, China
| | - Xiaoxiang Zheng
- Qiushi Academy for Advanced Studies, Zhejiang University, Hangzhou, 310027, China
| | - Gang Pan
- State Key Lab of CAD&CG, and College of Computer Science, Zhejiang University, Hangzhou, 310027, China
| | - Yiwen Wang
- Department of Electronic and Computer Engineering and Department of Chemical and Biological Engineering, Hong Kong University of Science and Technology, Kowloon, Hong Kong SAR, 999077, China
| |
Collapse
|
18
|
Extracting neuronal functional network dynamics via adaptive Granger causality analysis. Proc Natl Acad Sci U S A 2018; 115:E3869-E3878. [PMID: 29632213 DOI: 10.1073/pnas.1718154115] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Quantifying the functional relations between the nodes in a network based on local observations is a key challenge in studying complex systems. Most existing time series analysis techniques for this purpose provide static estimates of the network properties, pertain to stationary Gaussian data, or do not take into account the ubiquitous sparsity in the underlying functional networks. When applied to spike recordings from neuronal ensembles undergoing rapid task-dependent dynamics, they thus hinder a precise statistical characterization of the dynamic neuronal functional networks underlying adaptive behavior. We develop a dynamic estimation and inference paradigm for extracting functional neuronal network dynamics in the sense of Granger, by integrating techniques from adaptive filtering, compressed sensing, point process theory, and high-dimensional statistics. We demonstrate the utility of our proposed paradigm through theoretical analysis, algorithm development, and application to synthetic and real data. Application of our techniques to two-photon Ca2+ imaging experiments from the mouse auditory cortex reveals unique features of the functional neuronal network structures underlying spontaneous activity at unprecedented spatiotemporal resolution. Our analysis of simultaneous recordings from the ferret auditory and prefrontal cortical areas suggests evidence for the role of rapid top-down and bottom-up functional dynamics across these areas involved in robust attentive behavior.
Collapse
|
19
|
Karbasi A, Salavati AH, Vetterli M. Learning neural connectivity from firing activity: efficient algorithms with provable guarantees on topology. J Comput Neurosci 2018; 44:253-272. [PMID: 29464489 PMCID: PMC5851696 DOI: 10.1007/s10827-018-0678-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2017] [Revised: 01/06/2018] [Accepted: 01/22/2018] [Indexed: 11/18/2022]
Abstract
The connectivity of a neuronal network has a major effect on its functionality and role. It is generally believed that the complex network structure of the brain provides a physiological basis for information processing. Therefore, identifying the network's topology has received a lot of attentions in neuroscience and has been the center of many research initiatives such as Human Connectome Project. Nevertheless, direct and invasive approaches that slice and observe the neural tissue have proven to be time consuming, complex and costly. As a result, the inverse methods that utilize firing activity of neurons in order to identify the (functional) connections have gained momentum recently, especially in light of rapid advances in recording technologies; It will soon be possible to simultaneously monitor the activities of tens of thousands of neurons in real time. While there are a number of excellent approaches that aim to identify the functional connections from firing activities, the scalability of the proposed techniques plays a major challenge in applying them on large-scale datasets of recorded firing activities. In exceptional cases where scalability has not been an issue, the theoretical performance guarantees are usually limited to a specific family of neurons or the type of firing activities. In this paper, we formulate the neural network reconstruction as an instance of a graph learning problem, where we observe the behavior of nodes/neurons (i.e., firing activities) and aim to find the links/connections. We develop a scalable learning mechanism and derive the conditions under which the estimated graph for a network of Leaky Integrate and Fire (LIf) neurons matches the true underlying synaptic connections. We then validate the performance of the algorithm using artificially generated data (for benchmarking) and real data recorded from multiple hippocampal areas in rats.
Collapse
Affiliation(s)
- Amin Karbasi
- Inference, Information and Decision Systems Group, Yale Institute for Network Science, Yale University, New Haven, CT 06520 USA
| | - Amir Hesam Salavati
- Laboratory of Audiovisual Communications (LCAV), School of Computer and Communication Sciences, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
| | - Martin Vetterli
- Laboratory of Audiovisual Communications (LCAV), School of Computer and Communication Sciences, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
20
|
Magrans de Abril I, Yoshimoto J, Doya K. Connectivity inference from neural recording data: Challenges, mathematical bases and research directions. Neural Netw 2018; 102:120-137. [PMID: 29571122 DOI: 10.1016/j.neunet.2018.02.016] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2016] [Revised: 02/23/2018] [Accepted: 02/26/2018] [Indexed: 11/30/2022]
Abstract
This article presents a review of computational methods for connectivity inference from neural activity data derived from multi-electrode recordings or fluorescence imaging. We first identify biophysical and technical challenges in connectivity inference along the data processing pipeline. We then review connectivity inference methods based on two major mathematical foundations, namely, descriptive model-free approaches and generative model-based approaches. We investigate representative studies in both categories and clarify which challenges have been addressed by which method. We further identify critical open issues and possible research directions.
Collapse
Affiliation(s)
| | | | - Kenji Doya
- Okinawa Institute of Science and Technology, Graduate University, Japan
| |
Collapse
|
21
|
Hu M, Li W, Liang H. A Copula-Based Granger Causality Measure for the Analysis of Neural Spike Train Data. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:562-569. [PMID: 29610104 DOI: 10.1109/tcbb.2014.2388311] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
In systems neuroscience, it is becoming increasingly common to record the activity of hundreds of neurons simultaneously via electrode arrays. The ability to accurately measure the causal interactions among multiple neurons in the brain is crucial to understanding how neurons work in concert to generate specific brain functions. The development of new statistical methods for assessing causal influence between spike trains is still an active field of neuroscience research. Here, we suggest a copula-based Granger causality measure for the analysis of neural spike train data. This method is built upon our recent work on copula Granger causality for the analysis of continuous-valued time series by extending it to point-process neural spike train data. The proposed method is therefore able to reveal nonlinear and high-order causality in the spike trains while retaining all the computational advantages such as model-free, efficient estimation, and variability assessment of Granger causality. The performance of our algorithm can be further boosted with time-reversed data. Our method performed well on extensive simulations, and was then demonstrated on neural activity simultaneously recorded from primary visual cortex of a monkey performing a contour detection task.
Collapse
|
22
|
Wedgwood KCA, Satin LS. Six degrees of depolarization: Comment on "Network science of biological systems at different scales: A review" by Marko Gosak et al. Phys Life Rev 2018; 24:136-139. [PMID: 29395878 DOI: 10.1016/j.plrev.2018.01.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2017] [Accepted: 01/18/2018] [Indexed: 11/17/2022]
Affiliation(s)
| | - Leslie S Satin
- University of Michigan Medical School, Ann Arbor, MI, USA
| |
Collapse
|
23
|
Xiao X, Liu B, Zhang J, Xiao X, Pan Y. An Optimized Method for Bayesian Connectivity Change Point Model. J Comput Biol 2017; 25:337-347. [PMID: 29185805 DOI: 10.1089/cmb.2017.0154] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
The brain undergoes functional dynamic changes at all times. Investigating functional dynamics has been recently verified to be helpful for detecting psychological conditions and powerful for analyzing disease-related abnormalities of the brain. This article aims to detect functional dynamics. Specifically, we focus on how to effectively distinguish corresponding functional connectivity and change points from functional magnetic resonance imaging (fMRI) data. By combining Bayesian connectivity change point model (BCCPM), a modified genetic algorithm (GA) is presented to optimize the evolutionary procedure toward the most probable distributions of real change points in fMRI. We randomly initialize different binary indicator vectors to represent different distributions of change points. Each indicator vector represents an individual in GA, and together they form an initial population. Then we calculate Bayesian posterior probability and use it as the fitness of each individual. Finally, we evolve individuals of current generation toward the next higher fitness generation by a series of modified genetic operators. After several evolutionary procedures, individuals in the final generation may have outstanding fitness and the one with highest fitness can represent the most likely change point distribution in the corresponding fMRI data. Furthermore, the most probable change point distribution could be resolved. We test the optimized method for BCCPM on several synthesized data sets, and the experimental results verify that the proposed model produces higher accuracy results with lower time consumption. Also, we apply the new model to real block-designed task-based fMRI data set and excellent results are obtained.
Collapse
Affiliation(s)
- Xiuchun Xiao
- 1 College of Electronic and Information Engineering, Guangdong Ocean University , Zhanjiang, China .,2 Department of Computer Science, Georgia State University , Atlanta, Georgia
| | - Bing Liu
- 3 Department of Mathematics and Statistics, Georgia State University , Atlanta, Georgia
| | - Jing Zhang
- 3 Department of Mathematics and Statistics, Georgia State University , Atlanta, Georgia
| | - Xueli Xiao
- 2 Department of Computer Science, Georgia State University , Atlanta, Georgia
| | - Yi Pan
- 2 Department of Computer Science, Georgia State University , Atlanta, Georgia
| |
Collapse
|
24
|
Donner C, Obermayer K, Shimazaki H. Approximate Inference for Time-Varying Interactions and Macroscopic Dynamics of Neural Populations. PLoS Comput Biol 2017; 13:e1005309. [PMID: 28095421 PMCID: PMC5283755 DOI: 10.1371/journal.pcbi.1005309] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2016] [Revised: 01/31/2017] [Accepted: 12/12/2016] [Indexed: 11/29/2022] Open
Abstract
The models in statistical physics such as an Ising model offer a convenient way to characterize stationary activity of neural populations. Such stationary activity of neurons may be expected for recordings from in vitro slices or anesthetized animals. However, modeling activity of cortical circuitries of awake animals has been more challenging because both spike-rates and interactions can change according to sensory stimulation, behavior, or an internal state of the brain. Previous approaches modeling the dynamics of neural interactions suffer from computational cost; therefore, its application was limited to only a dozen neurons. Here by introducing multiple analytic approximation methods to a state-space model of neural population activity, we make it possible to estimate dynamic pairwise interactions of up to 60 neurons. More specifically, we applied the pseudolikelihood approximation to the state-space model, and combined it with the Bethe or TAP mean-field approximation to make the sequential Bayesian estimation of the model parameters possible. The large-scale analysis allows us to investigate dynamics of macroscopic properties of neural circuitries underlying stimulus processing and behavior. We show that the model accurately estimates dynamics of network properties such as sparseness, entropy, and heat capacity by simulated data, and demonstrate utilities of these measures by analyzing activity of monkey V4 neurons as well as a simulated balanced network of spiking neurons. Simultaneous analysis of large-scale neural populations is necessary to understand coding principles of neurons because they concertedly process information. Methods of thermodynamics and statistical mechanics are useful to understand collective phenomena of the interacting elements, and they have been successfully used to understand diverse activity of neurons. However, most analysis methods assume stationary data, in which activity rates of neurons and their correlations are constant over time. This assumption is easily violated in the data recorded from awake animals. Neural correlations likely organize dynamically during behavior and cognition, and this may be independent from the modulated activity rates of individual neurons. Recently several methods were proposed to simultaneously estimate dynamics of neural interactions. However, these methods are applicable to up to about 10 neurons. Here by combining multiple analytic approximation methods, we made it possible to estimate time-varying interactions of much larger neural populations. The method allows us to trace dynamic macroscopic properties of neural circuitries such as sparseness, entropy, and sensitivity. Using these statistics, researchers can now quantify to what extent neurons are correlated or de-correlated, and test if neural systems are susceptible within a specific behavioral period.
Collapse
Affiliation(s)
- Christian Donner
- Bernstein Center for Computational Neuroscience, Berlin, Germany
- Neural Information Processing Group, Department of Electrical Engineering and Computer Science, Technische Universität Berlin, Berlin, Germany
- Group for Methods of Artificial Intelligence, Department of Electrical Engineering and Computer Science, Technische Universität Berlin, Berlin, Germany
| | - Klaus Obermayer
- Bernstein Center for Computational Neuroscience, Berlin, Germany
- Neural Information Processing Group, Department of Electrical Engineering and Computer Science, Technische Universität Berlin, Berlin, Germany
| | | |
Collapse
|
25
|
Maren AJ. The Cluster Variation Method: A Primer for Neuroscientists. Brain Sci 2016; 6:E44. [PMID: 27706022 PMCID: PMC5187558 DOI: 10.3390/brainsci6040044] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Revised: 09/14/2016] [Accepted: 09/15/2016] [Indexed: 11/24/2022] Open
Abstract
Effective Brain-Computer Interfaces (BCIs) require that the time-varying activation patterns of 2-D neural ensembles be modelled. The cluster variation method (CVM) offers a means for the characterization of 2-D local pattern distributions. This paper provides neuroscientists and BCI researchers with a CVM tutorial that will help them to understand how the CVM statistical thermodynamics formulation can model 2-D pattern distributions expressing structural and functional dynamics in the brain. The premise is that local-in-time free energy minimization works alongside neural connectivity adaptation, supporting the development and stabilization of consistent stimulus-specific responsive activation patterns. The equilibrium distribution of local patterns, or configuration variables, is defined in terms of a single interaction enthalpy parameter (h) for the case of an equiprobable distribution of bistate (neural/neural ensemble) units. Thus, either one enthalpy parameter (or two, for the case of non-equiprobable distribution) yields equilibrium configuration variable values. Modeling 2-D neural activation distribution patterns with the representational layer of a computational engine, we can thus correlate variational free energy minimization with specific configuration variable distributions. The CVM triplet configuration variables also map well to the notion of a M = 3 functional motif. This paper addresses the special case of an equiprobable unit distribution, for which an analytic solution can be found.
Collapse
Affiliation(s)
- Alianna J Maren
- Northwestern University School of Professional Studies, Master of Science in Predictive Analytics Program, 405 Church St, Evanston, IL 60201, USA.
| |
Collapse
|
26
|
Consistent estimation of complete neuronal connectivity in large neuronal populations using sparse "shotgun" neuronal activity sampling. J Comput Neurosci 2016; 41:157-84. [PMID: 27515518 DOI: 10.1007/s10827-016-0611-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2015] [Revised: 06/09/2016] [Accepted: 06/13/2016] [Indexed: 11/27/2022]
Abstract
We investigate the properties of recently proposed "shotgun" sampling approach for the common inputs problem in the functional estimation of neuronal connectivity. We study the asymptotic correctness, the speed of convergence, and the data size requirements of such an approach. We show that the shotgun approach can be expected to allow the inference of complete connectivity matrix in large neuronal populations under some rather general conditions. However, we find that the posterior error of the shotgun connectivity estimator grows quickly with the size of unobserved neuronal populations, the square of average connectivity strength, and the square of observation sparseness. This implies that the shotgun connectivity estimation will require significantly larger amounts of neuronal activity data whenever the number of neurons in observed neuronal populations remains small. We present a numerical approach for solving the shotgun estimation problem in general settings and use it to demonstrate the shotgun connectivity inference in the examples of simulated synfire and weakly coupled cortical neuronal networks.
Collapse
|
27
|
Zhang C, Chai Y, Guo X, Gao M, Devilbiss D, Zhang Z. Statistical Learning of Neuronal Functional Connectivity. Technometrics 2016. [DOI: 10.1080/00401706.2016.1142904] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Chunming Zhang
- Department of Statistics, University of Wisconsin–Madison, Madison, WI 53706
| | | | - Xiao Guo
- Department of Statistics and Finance, University of Science and Technology of China, Hefei, Anhui, 230026, P.R. China
| | - Muhong Gao
- Department of Statistics, University of Wisconsin–Madison, Madison, WI 53706
| | - David Devilbiss
- Department of Psychology, University of Wisconsin–Madison, 53706, Madison, WI
| | - Zhengjun Zhang
- Department of Statistics, University of Wisconsin–Madison, 53706, Madison, WI
| |
Collapse
|
28
|
Oba S, Nakae K, Ikegaya Y, Aki S, Yoshimoto J, Ishii S. Empirical Bayesian significance measure of neuronal spike response. BMC Neurosci 2016; 17:27. [PMID: 27209433 PMCID: PMC4875706 DOI: 10.1186/s12868-016-0255-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2015] [Accepted: 05/10/2016] [Indexed: 12/01/2022] Open
Abstract
Background Functional connectivity analyses of multiple neurons provide a powerful bottom-up approach to reveal functions of local neuronal circuits by using simultaneous recording of neuronal activity. A statistical methodology, generalized linear modeling (GLM) of the spike response function, is one of the most promising methodologies to reduce false link discoveries arising from pseudo-correlation based on common inputs. Although recent advancement of fluorescent imaging techniques has increased the number of simultaneously recoded neurons up to the hundreds or thousands, the amount of information per pair of neurons has not correspondingly increased, partly because of the instruments’ limitations, and partly because the number of neuron pairs increase in a quadratic manner. Consequently, the estimation of GLM suffers from large statistical uncertainty caused by the shortage in effective information. Results In this study, we propose a new combination of GLM and empirical Bayesian testing for the estimation of spike response functions that enables both conservative false discovery control and powerful functional connectivity detection. We compared our proposed method’s performance with those of sparse estimation of GLM and classical Granger causality testing. Our method achieved high detection performance of functional connectivity with conservative estimation of false discovery rate and q values in case of information shortage due to short observation time. We also showed that empirical Bayesian testing on arbitrary statistics in place of likelihood-ratio statistics reduce the computational cost without decreasing the detection performance. When our proposed method was applied to a functional multi-neuron calcium imaging dataset from the rat hippocampal region, we found significant functional connections that are possibly mediated by AMPA and NMDA receptors. Conclusions The proposed empirical Bayesian testing framework with GLM is promising especially when the amount of information per a neuron pair is small because of growing size of observed network. Electronic supplementary material The online version of this article (doi:10.1186/s12868-016-0255-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Shigeyuki Oba
- Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Kyoto, Japan.
| | - Ken Nakae
- Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Kyoto, Japan
| | - Yuji Ikegaya
- Graduate School of Pharmaceutical Sciences, University of Tokyo, Tokyo, Japan
| | - Shunsuke Aki
- Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Kyoto, Japan
| | - Junichiro Yoshimoto
- Graduate School of Information Science, Nara Institute of Science and Technology, Ikoma, Japan
| | - Shin Ishii
- Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Kyoto, Japan
| |
Collapse
|
29
|
Reconstruction of recurrent synaptic connectivity of thousands of neurons from simulated spiking activity. J Comput Neurosci 2015; 39:77-103. [PMID: 26041729 PMCID: PMC4493949 DOI: 10.1007/s10827-015-0565-5] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2014] [Revised: 04/18/2015] [Accepted: 04/22/2015] [Indexed: 10/30/2022]
Abstract
Dynamics and function of neuronal networks are determined by their synaptic connectivity. Current experimental methods to analyze synaptic network structure on the cellular level, however, cover only small fractions of functional neuronal circuits, typically without a simultaneous record of neuronal spiking activity. Here we present a method for the reconstruction of large recurrent neuronal networks from thousands of parallel spike train recordings. We employ maximum likelihood estimation of a generalized linear model of the spiking activity in continuous time. For this model the point process likelihood is concave, such that a global optimum of the parameters can be obtained by gradient ascent. Previous methods, including those of the same class, did not allow recurrent networks of that order of magnitude to be reconstructed due to prohibitive computational cost and numerical instabilities. We describe a minimal model that is optimized for large networks and an efficient scheme for its parallelized numerical optimization on generic computing clusters. For a simulated balanced random network of 1000 neurons, synaptic connectivity is recovered with a misclassification error rate of less than 1 % under ideal conditions. We show that the error rate remains low in a series of example cases under progressively less ideal conditions. Finally, we successfully reconstruct the connectivity of a hidden synfire chain that is embedded in a random network, which requires clustering of the network connectivity to reveal the synfire groups. Our results demonstrate how synaptic connectivity could potentially be inferred from large-scale parallel spike train recordings.
Collapse
|
30
|
Song W, Cajigas I, Brown EN, Giszter SF. Adaptation to elastic loads and BMI robot controls during rat locomotion examined with point-process GLMs. Front Syst Neurosci 2015; 9:62. [PMID: 25972789 PMCID: PMC4411868 DOI: 10.3389/fnsys.2015.00062] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2015] [Accepted: 04/01/2015] [Indexed: 11/13/2022] Open
Abstract
Currently little is known about how a mechanically coupled BMI system's actions are integrated into ongoing body dynamics. We tested a locomotor task augmented with a BMI system driving a robot mechanically interacting with a rat under three conditions: control locomotion (BL), “simple elastic load” (E) and “BMI with elastic load” (BMI/E). The effect of the BMI was to allow compensation of the elastic load as a function of the neural drive. Neurons recorded here were close to one another in cortex, all within a 200 micron diameter horizontal distance of one another. The interactions of these close assemblies of neurons may differ from those among neurons at longer distances in BMI tasks and thus are important to explore. A point process generalized linear model (GLM), was used to examine connectivity at two different binning timescales (1 ms vs. 10 ms). We used GLM models to fit non-Poisson neural dynamics solely using other neurons' prior neural activity as covariates. Models at different timescales were compared based on Kolmogorov-Smirnov (KS) goodness-of-fit and parsimony. About 15% of cells with non-Poisson firing were well fitted with the neuron-to-neuron models alone. More such cells were fitted at the 1 ms binning than 10 ms. Positive connection parameters (“excitation” ~70%) exceeded negative parameters (“inhibition” ~30%). Significant connectivity changes in the GLM determined networks of well-fitted neurons occurred between the conditions. However, a common core of connections comprising at least ~15% of connections persisted between any two of the three conditions. Significantly almost twice as many connections were in common between the two load conditions (~27%), compared to between either load condition and the baseline. This local point process GLM identified neural correlation structure and the changes seen across task conditions in the rats in this neural subset may be intrinsic to cortex or due to feedback and input reorganization in adaptation.
Collapse
Affiliation(s)
- Weiguo Song
- Department of Neurobiology and Anatomy, Drexel University College of Medicine, Drexel University Philadelphia, PA, USA
| | - Iahn Cajigas
- Brain and Cognitive Sciences, Massachusetts Institute of Technology Cambridge, MA, USA
| | - Emery N Brown
- Brain and Cognitive Sciences, Massachusetts Institute of Technology Cambridge, MA, USA ; Institute for Medical Engineering and Science, Massachusetts Institute of Technology Cambridge, MA, USA ; Department of Anesthesia, Critical Care and Pain Medicine, Harvard Medical School, Massachusetts General Hospital Boston, MA, USA
| | - Simon F Giszter
- Department of Neurobiology and Anatomy, Drexel University College of Medicine, Drexel University Philadelphia, PA, USA ; School of Biomedical Engineering, Drexel University Philadelphia, PA, USA
| |
Collapse
|
31
|
Volgushev M, Ilin V, Stevenson IH. Identifying and tracking simulated synaptic inputs from neuronal firing: insights from in vitro experiments. PLoS Comput Biol 2015; 11:e1004167. [PMID: 25823000 PMCID: PMC4379067 DOI: 10.1371/journal.pcbi.1004167] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2014] [Accepted: 02/02/2015] [Indexed: 11/18/2022] Open
Abstract
Accurately describing synaptic interactions between neurons and how interactions change over time are key challenges for systems neuroscience. Although intracellular electrophysiology is a powerful tool for studying synaptic integration and plasticity, it is limited by the small number of neurons that can be recorded simultaneously in vitro and by the technical difficulty of intracellular recording in vivo. One way around these difficulties may be to use large-scale extracellular recording of spike trains and apply statistical methods to model and infer functional connections between neurons. These techniques have the potential to reveal large-scale connectivity structure based on the spike timing alone. However, the interpretation of functional connectivity is often approximate, since only a small fraction of presynaptic inputs are typically observed. Here we use in vitro current injection in layer 2/3 pyramidal neurons to validate methods for inferring functional connectivity in a setting where input to the neuron is controlled. In experiments with partially-defined input, we inject a single simulated input with known amplitude on a background of fluctuating noise. In a fully-defined input paradigm, we then control the synaptic weights and timing of many simulated presynaptic neurons. By analyzing the firing of neurons in response to these artificial inputs, we ask 1) How does functional connectivity inferred from spikes relate to simulated synaptic input? and 2) What are the limitations of connectivity inference? We find that individual current-based synaptic inputs are detectable over a broad range of amplitudes and conditions. Detectability depends on input amplitude and output firing rate, and excitatory inputs are detected more readily than inhibitory. Moreover, as we model increasing numbers of presynaptic inputs, we are able to estimate connection strengths more accurately and detect the presence of connections more quickly. These results illustrate the possibilities and outline the limits of inferring synaptic input from spikes.
Collapse
Affiliation(s)
- Maxim Volgushev
- Department of Psychology, University of Connecticut, Storrs, Connecticut, United States of America
| | - Vladimir Ilin
- Department of Psychology, University of Connecticut, Storrs, Connecticut, United States of America
| | - Ian H. Stevenson
- Department of Psychology, University of Connecticut, Storrs, Connecticut, United States of America
| |
Collapse
|
32
|
Wei J, Bai W, Liu T, Tian X. Functional connectivity changes during a working memory task in rat via NMF analysis. Front Behav Neurosci 2015; 9:2. [PMID: 25688192 PMCID: PMC4311635 DOI: 10.3389/fnbeh.2015.00002] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2014] [Accepted: 01/02/2015] [Indexed: 02/01/2023] Open
Abstract
Working memory (WM) is necessary in higher cognition. The brain as a complex network is formed by interconnections among neurons. Connectivity results in neural dynamics to support cognition. The first aim is to investigate connectivity dynamics in medial prefrontal cortex (mPFC) networks during WM. As brain neural activity is sparse, the second aim is to find the intrinsic connectivity property in a feature space. Using multi-channel electrode recording techniques, spikes were simultaneously obtained from mPFC of rats that performed a Y-maze WM task. Continuous time series converted from spikes were embedded in a low-dimensional space by non-negative matrix factorization (NMF). mPFC network in original space was constructed by measuring connections among neurons. And the same network in NMF space was constructed by computing connectivity values between the extracted NMF components. Causal density (Cd) and global efficiency (E) were estimated to present the network property. The results showed that Cd and E significantly peaked in the interval right before the maze choice point in correct trials. However, the increase did not emerge in error trials. Additionally, Cd and E in two spaces displayed similar trends in correct trials. The difference was that the measures in NMF space were significantly greater than those in original space. Our findings indicated that the anticipatory changes in mPFC networks may have an effect on future WM behavioral choices. Moreover, the NMF analysis achieves a better characterization for a brain network.
Collapse
Affiliation(s)
- Jing Wei
- School of Biomedical Engineering, Tianjin Medical University Tianjin, China
| | - Wenwen Bai
- School of Biomedical Engineering, Tianjin Medical University Tianjin, China
| | - Tiaotiao Liu
- School of Biomedical Engineering, Tianjin Medical University Tianjin, China
| | - Xin Tian
- School of Biomedical Engineering, Tianjin Medical University Tianjin, China ; Research Center of Basic Medicine, Tianjin Medical University Tianjin, China
| |
Collapse
|
33
|
Nakae K, Ikegaya Y, Ishikawa T, Oba S, Urakubo H, Koyama M, Ishii S. A statistical method of identifying interactions in neuron-glia systems based on functional multicell Ca2+ imaging. PLoS Comput Biol 2014; 10:e1003949. [PMID: 25393874 PMCID: PMC4230777 DOI: 10.1371/journal.pcbi.1003949] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2013] [Accepted: 09/29/2014] [Indexed: 11/18/2022] Open
Abstract
Crosstalk between neurons and glia may constitute a significant part of information processing in the brain. We present a novel method of statistically identifying interactions in a neuron-glia network. We attempted to identify neuron-glia interactions from neuronal and glial activities via maximum-a-posteriori (MAP)-based parameter estimation by developing a generalized linear model (GLM) of a neuron-glia network. The interactions in our interest included functional connectivity and response functions. We evaluated the cross-validated likelihood of GLMs that resulted from the addition or removal of connections to confirm the existence of specific neuron-to-glia or glia-to-neuron connections. We only accepted addition or removal when the modification improved the cross-validated likelihood. We applied the method to a high-throughput, multicellular in vitro Ca2+ imaging dataset obtained from the CA3 region of a rat hippocampus, and then evaluated the reliability of connectivity estimates using a statistical test based on a surrogate method. Our findings based on the estimated connectivity were in good agreement with currently available physiological knowledge, suggesting our method can elucidate undiscovered functions of neuron-glia systems.
Collapse
Affiliation(s)
- Ken Nakae
- Integrated Systems Biology Laboratory, Graduate School of Informatics, Kyoto University, Sakyo-ku, Kyoto, Japan
| | - Yuji Ikegaya
- Laboratory of Chemical Pharmacology, Graduate School of Pharmaceutical Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
- Center for Information and Neural Networks, Suita City, Osaka, Japan
- * E-mail: (YI); (SI)
| | - Tomoe Ishikawa
- Laboratory of Chemical Pharmacology, Graduate School of Pharmaceutical Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Shigeyuki Oba
- Integrated Systems Biology Laboratory, Graduate School of Informatics, Kyoto University, Sakyo-ku, Kyoto, Japan
| | - Hidetoshi Urakubo
- Integrated Systems Biology Laboratory, Graduate School of Informatics, Kyoto University, Sakyo-ku, Kyoto, Japan
| | - Masanori Koyama
- Integrated Systems Biology Laboratory, Graduate School of Informatics, Kyoto University, Sakyo-ku, Kyoto, Japan
| | - Shin Ishii
- Integrated Systems Biology Laboratory, Graduate School of Informatics, Kyoto University, Sakyo-ku, Kyoto, Japan
- * E-mail: (YI); (SI)
| |
Collapse
|
34
|
On Stationary Distributions of Stochastic Neural Networks. J Appl Probab 2014. [DOI: 10.1017/s0021900200011700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The paper deals with nonlinear Poisson neuron network models with bounded memory dynamics, which can include both Hebbian learning mechanisms and refractory periods. The state of the network is described by the times elapsed since its neurons fired within the post-synaptic transfer kernel memory span, and the current strengths of synaptic connections, the state spaces of our models being hierarchies of finite-dimensional components. We prove the ergodicity of the stochastic processes describing the behaviour of the networks, establish the existence of continuously differentiable stationary distribution densities (with respect to the Lebesgue measures of corresponding dimensionality) on the components of the state space, and find upper bounds for them. For the density components, we derive a system of differential equations that can be solved in a few simplest cases only. Approaches to approximate computation of the stationary density are discussed. One approach is to reduce the dimensionality of the problem by modifying the network so that each neuron cannot fire if the number of spikes it emitted within the post-synaptic transfer kernel memory span reaches a given threshold. We show that the stationary distribution of this ‘truncated’ network converges to that of the unrestricted network as the threshold increases, and that the convergence is at a superexponential rate. A complementary approach uses discrete Markov chain approximations to the network process.
Collapse
|
35
|
Borovkov K, Decrouez G, Gilson M. On Stationary Distributions of Stochastic Neural Networks. J Appl Probab 2014. [DOI: 10.1239/jap/1409932677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The paper deals with nonlinear Poisson neuron network models with bounded memory dynamics, which can include both Hebbian learning mechanisms and refractory periods. The state of the network is described by the times elapsed since its neurons fired within the post-synaptic transfer kernel memory span, and the current strengths of synaptic connections, the state spaces of our models being hierarchies of finite-dimensional components. We prove the ergodicity of the stochastic processes describing the behaviour of the networks, establish the existence of continuously differentiable stationary distribution densities (with respect to the Lebesgue measures of corresponding dimensionality) on the components of the state space, and find upper bounds for them. For the density components, we derive a system of differential equations that can be solved in a few simplest cases only. Approaches to approximate computation of the stationary density are discussed. One approach is to reduce the dimensionality of the problem by modifying the network so that each neuron cannot fire if the number of spikes it emitted within the post-synaptic transfer kernel memory span reaches a given threshold. We show that the stationary distribution of this ‘truncated’ network converges to that of the unrestricted network as the threshold increases, and that the convergence is at a superexponential rate. A complementary approach uses discrete Markov chain approximations to the network process.
Collapse
|
36
|
Chen Z. An overview of Bayesian methods for neural spike train analysis. COMPUTATIONAL INTELLIGENCE AND NEUROSCIENCE 2013; 2013:251905. [PMID: 24348527 PMCID: PMC3855941 DOI: 10.1155/2013/251905] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 05/16/2013] [Revised: 09/10/2013] [Accepted: 09/23/2013] [Indexed: 11/24/2022]
Abstract
Neural spike train analysis is an important task in computational neuroscience which aims to understand neural mechanisms and gain insights into neural circuits. With the advancement of multielectrode recording and imaging technologies, it has become increasingly demanding to develop statistical tools for analyzing large neuronal ensemble spike activity. Here we present a tutorial overview of Bayesian methods and their representative applications in neural spike train analysis, at both single neuron and population levels. On the theoretical side, we focus on various approximate Bayesian inference techniques as applied to latent state and parameter estimation. On the application side, the topics include spike sorting, tuning curve estimation, neural encoding and decoding, deconvolution of spike trains from calcium imaging signals, and inference of neuronal functional connectivity and synchrony. Some research challenges and opportunities for neural spike train analysis are discussed.
Collapse
Affiliation(s)
- Zhe Chen
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, 43 Vassar Street, Cambridge, MA 02139, USA
- Picower Institute for Learning and Memory, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| |
Collapse
|
37
|
Hamilton F, Berry T, Peixoto N, Sauer T. Real-time tracking of neuronal network structure using data assimilation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:052715. [PMID: 24329304 DOI: 10.1103/physreve.88.052715] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2013] [Revised: 09/18/2013] [Indexed: 05/09/2023]
Abstract
A nonlinear data assimilation technique is applied to determine and track effective connections between ensembles of cultured spinal cord neurons measured with multielectrode arrays. The method is statistical, depending only on confidence intervals, and requiring no form of arbitrary thresholding. In addition, the method updates connection strengths sequentially, enabling real-time tracking of nonstationary networks. The ensemble Kalman filter is used with a generic spiking neuron model to estimate connection strengths as well as other system parameters to deal with model mismatch. The method is validated on noisy synthetic data from Hodgkin-Huxley model neurons before being used to find network connections in the neural culture recordings.
Collapse
Affiliation(s)
- Franz Hamilton
- Electrical and Computer Engineering, George Mason University, Fairfax, Virginia 22030, USA
| | - Tyrus Berry
- Mathematical Sciences, George Mason University, Fairfax, Virginia 22030, USA
| | - Nathalia Peixoto
- Electrical and Computer Engineering, George Mason University, Fairfax, Virginia 22030, USA
| | - Timothy Sauer
- Mathematical Sciences, George Mason University, Fairfax, Virginia 22030, USA
| |
Collapse
|
38
|
Identification of sparse neural functional connectivity using penalized likelihood estimation and basis functions. J Comput Neurosci 2013; 35:335-57. [PMID: 23674048 DOI: 10.1007/s10827-013-0455-7] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2012] [Revised: 04/10/2013] [Accepted: 04/16/2013] [Indexed: 10/26/2022]
Abstract
One key problem in computational neuroscience and neural engineering is the identification and modeling of functional connectivity in the brain using spike train data. To reduce model complexity, alleviate overfitting, and thus facilitate model interpretation, sparse representation and estimation of functional connectivity is needed. Sparsities include global sparsity, which captures the sparse connectivities between neurons, and local sparsity, which reflects the active temporal ranges of the input-output dynamical interactions. In this paper, we formulate a generalized functional additive model (GFAM) and develop the associated penalized likelihood estimation methods for such a modeling problem. A GFAM consists of a set of basis functions convolving the input signals, and a link function generating the firing probability of the output neuron from the summation of the convolutions weighted by the sought model coefficients. Model sparsities are achieved by using various penalized likelihood estimations and basis functions. Specifically, we introduce two variations of the GFAM using a global basis (e.g., Laguerre basis) and group LASSO estimation, and a local basis (e.g., B-spline basis) and group bridge estimation, respectively. We further develop an optimization method based on quadratic approximation of the likelihood function for the estimation of these models. Simulation and experimental results show that both group-LASSO-Laguerre and group-bridge-B-spline can capture faithfully the global sparsities, while the latter can replicate accurately and simultaneously both global and local sparsities. The sparse models outperform the full models estimated with the standard maximum likelihood method in out-of-sample predictions.
Collapse
|
39
|
Pillow JW, Shlens J, Chichilnisky EJ, Simoncelli EP. A model-based spike sorting algorithm for removing correlation artifacts in multi-neuron recordings. PLoS One 2013; 8:e62123. [PMID: 23671583 PMCID: PMC3643981 DOI: 10.1371/journal.pone.0062123] [Citation(s) in RCA: 81] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2012] [Accepted: 03/19/2013] [Indexed: 12/05/2022] Open
Abstract
We examine the problem of estimating the spike trains of multiple neurons from voltage traces recorded on one or more extracellular electrodes. Traditional spike-sorting methods rely on thresholding or clustering of recorded signals to identify spikes. While these methods can detect a large fraction of the spikes from a recording, they generally fail to identify synchronous or near-synchronous spikes: cases in which multiple spikes overlap. Here we investigate the geometry of failures in traditional sorting algorithms, and document the prevalence of such errors in multi-electrode recordings from primate retina. We then develop a method for multi-neuron spike sorting using a model that explicitly accounts for the superposition of spike waveforms. We model the recorded voltage traces as a linear combination of spike waveforms plus a stochastic background component of correlated Gaussian noise. Combining this measurement model with a Bernoulli prior over binary spike trains yields a posterior distribution for spikes given the recorded data. We introduce a greedy algorithm to maximize this posterior that we call "binary pursuit". The algorithm allows modest variability in spike waveforms and recovers spike times with higher precision than the voltage sampling rate. This method substantially corrects cross-correlation artifacts that arise with conventional methods, and substantially outperforms clustering methods on both real and simulated data. Finally, we develop diagnostic tools that can be used to assess errors in spike sorting in the absence of ground truth.
Collapse
Affiliation(s)
- Jonathan W Pillow
- Center for Perceptual Systems, Department of Psychology and Section of Neurobiology, The University of Texas at Austin, Austin, Texas, USA.
| | | | | | | |
Collapse
|
40
|
Functional connectivity and tuning curves in populations of simultaneously recorded neurons. PLoS Comput Biol 2012; 8:e1002775. [PMID: 23166484 PMCID: PMC3499254 DOI: 10.1371/journal.pcbi.1002775] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2012] [Accepted: 09/25/2012] [Indexed: 11/25/2022] Open
Abstract
How interactions between neurons relate to tuned neural responses is a longstanding question in systems neuroscience. Here we use statistical modeling and simultaneous multi-electrode recordings to explore the relationship between these interactions and tuning curves in six different brain areas. We find that, in most cases, functional interactions between neurons provide an explanation of spiking that complements and, in some cases, surpasses the influence of canonical tuning curves. Modeling functional interactions improves both encoding and decoding accuracy by accounting for noise correlations and features of the external world that tuning curves fail to capture. In cortex, modeling coupling alone allows spikes to be predicted more accurately than tuning curve models based on external variables. These results suggest that statistical models of functional interactions between even relatively small numbers of neurons may provide a useful framework for examining neural coding. The number of simultaneous neurons that electrophysiologists can record is growing rapidly, and a central goal of computational neuroscience is to develop statistical methods that can make sense of this growing data. Here we present a unified statistical analysis of 10 different datasets recorded from several different species and brain areas. We show how functional interactions between neurons may be used to predict spiking in each of these different areas, and find that, in many cases, modeling interactions between a small number of neurons yields better spike predictions than modeling each neuron's relationship to the outside world using tuning curves. Although these statistical results cannot be linked to specific network architectures, since the measured interactions between neurons are purely functional rather than anatomical, they suggest that modeling interactions between neurons will be a useful approach to understanding neural coding as electrophysiologists record from increasing numbers of neurons.
Collapse
|
41
|
Devilbiss DM, Jenison RL, Berridge CW. Stress-induced impairment of a working memory task: role of spiking rate and spiking history predicted discharge. PLoS Comput Biol 2012; 8:e1002681. [PMID: 23028279 PMCID: PMC3441423 DOI: 10.1371/journal.pcbi.1002681] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2012] [Accepted: 07/19/2012] [Indexed: 12/19/2022] Open
Abstract
Stress, pervasive in society, contributes to over half of all work place accidents a year and over time can contribute to a variety of psychiatric disorders including depression, schizophrenia, and post-traumatic stress disorder. Stress impairs higher cognitive processes, dependent on the prefrontal cortex (PFC) and that involve maintenance and integration of information over extended periods, including working memory and attention. Substantial evidence has demonstrated a relationship between patterns of PFC neuron spiking activity (action-potential discharge) and components of delayed-response tasks used to probe PFC-dependent cognitive function in rats and monkeys. During delay periods of these tasks, persistent spiking activity is posited to be essential for the maintenance of information for working memory and attention. However, the degree to which stress-induced impairment in PFC-dependent cognition involves changes in task-related spiking rates or the ability for PFC neurons to retain information over time remains unknown. In the current study, spiking activity was recorded from the medial PFC of rats performing a delayed-response task of working memory during acute noise stress (93 db). Spike history-predicted discharge (SHPD) for PFC neurons was quantified as a measure of the degree to which ongoing neuronal discharge can be predicted by past spiking activity and reflects the degree to which past information is retained by these neurons over time. We found that PFC neuron discharge is predicted by their past spiking patterns for nearly one second. Acute stress impaired SHPD, selectively during delay intervals of the task, and simultaneously impaired task performance. Despite the reduction in delay-related SHPD, stress increased delay-related spiking rates. These findings suggest that neural codes utilizing SHPD within PFC networks likely reflects an additional important neurophysiological mechanism for maintenance of past information over time. Stress-related impairment of this mechanism is posited to contribute to the cognition-impairing actions of stress.
Collapse
|
42
|
Kelly RC, Kass RE. A framework for evaluating pairwise and multiway synchrony among stimulus-driven neurons. Neural Comput 2012; 24:2007-32. [PMID: 22509967 PMCID: PMC3374919 DOI: 10.1162/neco_a_00307] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
Several authors have previously discussed the use of log-linear models, often called maximum entropy models, for analyzing spike train data to detect synchrony. The usual log-linear modeling techniques, however, do not allow time-varying firing rates that typically appear in stimulus-driven (or action-driven) neurons, nor do they incorporate non-Poisson history effects or covariate effects. We generalize the usual approach, combining point-process regression models of individual neuron activity with log-linear models of multiway synchronous interaction. The methods are illustrated with results found in spike trains recorded simultaneously from primary visual cortex. We then assess the amount of data needed to reliably detect multiway spiking.
Collapse
Affiliation(s)
- Ryan C Kelly
- Department of Statistics and Center for the Neural Basis of Cognition, Carnegie Mellon University, Pittsburgh, PA 15213, U.S.A.
| | | |
Collapse
|
43
|
Berry T, Hamilton F, Peixoto N, Sauer T. Detecting connectivity changes in neuronal networks. J Neurosci Methods 2012; 209:388-97. [PMID: 22771714 DOI: 10.1016/j.jneumeth.2012.06.021] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2012] [Revised: 06/19/2012] [Accepted: 06/20/2012] [Indexed: 10/28/2022]
Abstract
We develop a method from semiparametric statistics (Cox, 1972) for the purpose of tracking links and connection strengths over time in a neuronal network from spike train data. We consider application of the method as implemented in Masud and Borisyuk (2011), and evaluate its use on data generated independently of the Cox model hypothesis, in particular from a spiking model of Izhikevich in four different dynamical regimes. Then, we show how the Cox method can be used to determine statistically significant changes in network connectivity over time. Our methodology is demonstrated using spike trains from multi-electrode array measurements of networks of cultured mammalian spinal cord cells.
Collapse
Affiliation(s)
- Tyrus Berry
- Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA
| | | | | | | |
Collapse
|
44
|
Lawhern V, Hatsopoulos NG, Wu W. Coupling Time Decoding and Trajectory Decoding using a Target-Included Model in the Motor Cortex. Neurocomputing 2012; 82:117-126. [PMID: 22379284 PMCID: PMC3286625 DOI: 10.1016/j.neucom.2011.10.030] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
Significant progress has been made within the last decade in motor cortical decoding that predicts movement behaviors from population neuronal activity in the motor cortex. A majority of these decoding methods have focused on estimating a subject's hand trajectory in a continuous movement. We recently proposed a time identification decoding approach and showed that if a stereotyped movement is well represented by a sequence of targets (or landmarks), then the main structure of the movement can be reconstructed by detecting the reaching times at those targets. Both trajectory decoding and landmark-time decoding have their particular advantages, whereas a coupling of these two different strategies has not been examined. In this article we propose a synergy that comes from combining these two approaches for a stereotyped movement under a linear state-space framework. We develop a new decoding procedure based on a forward-backward propagation where the target is used in the initial stage in the backward step. Experimental results show that the new method significantly improves decoding accuracy over the non-target-included models. Furthermore, the coupling based on the new target-included method effectively combines the time decoding and trajectory decoding and further improves the decoding accuracy.
Collapse
Affiliation(s)
- Vernon Lawhern
- Department of Statistics, Florida State University, Tallahassee, FL 32306-4330, USA
| | - Nicholas G. Hatsopoulos
- Department of Organismal Biology and Anatomy, Committees on Computational Neuroscience and Neurobiology, University of Chicago, Chicago, IL 60637, USA
| | - Wei Wu
- Department of Statistics, Florida State University, Tallahassee, FL 32306-4330, USA
| |
Collapse
|
45
|
Mishchenko Y, Paninski L. A Bayesian compressed-sensing approach for reconstructing neural connectivity from subsampled anatomical data. J Comput Neurosci 2012; 33:371-88. [DOI: 10.1007/s10827-012-0390-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2011] [Revised: 02/09/2012] [Accepted: 03/05/2012] [Indexed: 10/28/2022]
|
46
|
State-space analysis of time-varying higher-order spike correlation for multiple neural spike train data. PLoS Comput Biol 2012; 8:e1002385. [PMID: 22412358 PMCID: PMC3297562 DOI: 10.1371/journal.pcbi.1002385] [Citation(s) in RCA: 96] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2011] [Accepted: 12/28/2011] [Indexed: 11/23/2022] Open
Abstract
Precise spike coordination between the spiking activities of multiple neurons is suggested as an indication of coordinated network activity in active cell assemblies. Spike correlation analysis aims to identify such cooperative network activity by detecting excess spike synchrony in simultaneously recorded multiple neural spike sequences. Cooperative activity is expected to organize dynamically during behavior and cognition; therefore currently available analysis techniques must be extended to enable the estimation of multiple time-varying spike interactions between neurons simultaneously. In particular, new methods must take advantage of the simultaneous observations of multiple neurons by addressing their higher-order dependencies, which cannot be revealed by pairwise analyses alone. In this paper, we develop a method for estimating time-varying spike interactions by means of a state-space analysis. Discretized parallel spike sequences are modeled as multi-variate binary processes using a log-linear model that provides a well-defined measure of higher-order spike correlation in an information geometry framework. We construct a recursive Bayesian filter/smoother for the extraction of spike interaction parameters. This method can simultaneously estimate the dynamic pairwise spike interactions of multiple single neurons, thereby extending the Ising/spin-glass model analysis of multiple neural spike train data to a nonstationary analysis. Furthermore, the method can estimate dynamic higher-order spike interactions. To validate the inclusion of the higher-order terms in the model, we construct an approximation method to assess the goodness-of-fit to spike data. In addition, we formulate a test method for the presence of higher-order spike correlation even in nonstationary spike data, e.g., data from awake behaving animals. The utility of the proposed methods is tested using simulated spike data with known underlying correlation dynamics. Finally, we apply the methods to neural spike data simultaneously recorded from the motor cortex of an awake monkey and demonstrate that the higher-order spike correlation organizes dynamically in relation to a behavioral demand. Nearly half a century ago, the Canadian psychologist D. O. Hebb postulated the formation of assemblies of tightly connected cells in cortical recurrent networks because of changes in synaptic weight (Hebb's learning rule) by repetitive sensory stimulation of the network. Consequently, the activation of such an assembly for processing sensory or behavioral information is likely to be expressed by precisely coordinated spiking activities of the participating neurons. However, the available analysis techniques for multiple parallel neural spike data do not allow us to reveal the detailed structure of transiently active assemblies as indicated by their dynamical pairwise and higher-order spike correlations. Here, we construct a state-space model of dynamic spike interactions, and present a recursive Bayesian method that makes it possible to trace multiple neurons exhibiting such precisely coordinated spiking activities in a time-varying manner. We also formulate a hypothesis test of the underlying dynamic spike correlation, which enables us to detect the assemblies activated in association with behavioral events. Therefore, the proposed method can serve as a useful tool to test Hebb's cell assembly hypothesis.
Collapse
|
47
|
Vidne M, Ahmadian Y, Shlens J, Pillow JW, Kulkarni J, Litke AM, Chichilnisky EJ, Simoncelli E, Paninski L. Modeling the impact of common noise inputs on the network activity of retinal ganglion cells. J Comput Neurosci 2011; 33:97-121. [PMID: 22203465 DOI: 10.1007/s10827-011-0376-2] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2011] [Revised: 12/04/2011] [Accepted: 12/09/2011] [Indexed: 10/14/2022]
Abstract
Synchronized spontaneous firing among retinal ganglion cells (RGCs), on timescales faster than visual responses, has been reported in many studies. Two candidate mechanisms of synchronized firing include direct coupling and shared noisy inputs. In neighboring parasol cells of primate retina, which exhibit rapid synchronized firing that has been studied extensively, recent experimental work indicates that direct electrical or synaptic coupling is weak, but shared synaptic input in the absence of modulated stimuli is strong. However, previous modeling efforts have not accounted for this aspect of firing in the parasol cell population. Here we develop a new model that incorporates the effects of common noise, and apply it to analyze the light responses and synchronized firing of a large, densely-sampled network of over 250 simultaneously recorded parasol cells. We use a generalized linear model in which the spike rate in each cell is determined by the linear combination of the spatio-temporally filtered visual input, the temporally filtered prior spikes of that cell, and unobserved sources representing common noise. The model accurately captures the statistical structure of the spike trains and the encoding of the visual stimulus, without the direct coupling assumption present in previous modeling work. Finally, we examined the problem of decoding the visual stimulus from the spike train given the estimated parameters. The common-noise model produces Bayesian decoding performance as accurate as that of a model with direct coupling, but with significantly more robustness to spike timing perturbations.
Collapse
Affiliation(s)
- Michael Vidne
- Department of Applied Physics & Applied Mathematics, Center for Theoretical Neuroscience, Columbia University, New York, NY, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|
48
|
Park M, Pillow JW. Receptive field inference with localized priors. PLoS Comput Biol 2011; 7:e1002219. [PMID: 22046110 PMCID: PMC3203052 DOI: 10.1371/journal.pcbi.1002219] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2011] [Accepted: 08/19/2011] [Indexed: 11/18/2022] Open
Abstract
The linear receptive field describes a mapping from sensory stimuli to a one-dimensional variable governing a neuron's spike response. However, traditional receptive field estimators such as the spike-triggered average converge slowly and often require large amounts of data. Bayesian methods seek to overcome this problem by biasing estimates towards solutions that are more likely a priori, typically those with small, smooth, or sparse coefficients. Here we introduce a novel Bayesian receptive field estimator designed to incorporate locality, a powerful form of prior information about receptive field structure. The key to our approach is a hierarchical receptive field model that flexibly adapts to localized structure in both spacetime and spatiotemporal frequency, using an inference method known as empirical Bayes. We refer to our method as automatic locality determination (ALD), and show that it can accurately recover various types of smooth, sparse, and localized receptive fields. We apply ALD to neural data from retinal ganglion cells and V1 simple cells, and find it achieves error rates several times lower than standard estimators. Thus, estimates of comparable accuracy can be achieved with substantially less data. Finally, we introduce a computationally efficient Markov Chain Monte Carlo (MCMC) algorithm for fully Bayesian inference under the ALD prior, yielding accurate Bayesian confidence intervals for small or noisy datasets.
Collapse
Affiliation(s)
- Mijung Park
- Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin, Texas, United States of America
| | - Jonathan W. Pillow
- Center for Perceptual Systems, Department of Psychology and Section of Neurobiology, The University of Texas at Austin, Austin, Texas, United States of America
| |
Collapse
|
49
|
Mishchenko Y, Paninski L. Efficient methods for sampling spike trains in networks of coupled neurons. Ann Appl Stat 2011. [DOI: 10.1214/11-aoas467] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
50
|
Mishchenko Y, Vogelstein JT, Paninski L. A Bayesian approach for inferring neuronal connectivity from calcium fluorescent imaging data. Ann Appl Stat 2011. [DOI: 10.1214/09-aoas303] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|