- Poster presentation
- Open Access
Bayes factor analysis for detection of time-dependent higher-order spike correlations
BMC Neurosciencevolume 10, Article number: P99 (2009)
Precise spike coordination in the spiking activities of a neuronal population is discussed as an indication of coordinated network activity in form of a cell assembly relevant for information processing. Supportive evidence was provided by the existence of excess spike synchrony occurring dynamically in relation to behavioral context [1–3]. These findings are based on the measured dependence of multiple neurons against the null-hypothesis of full independence. However, neurons jointly involved in assemblies may express higher-order correlations (HOCs) in their spiking activities . By characterizing the spatio-temporal pattern of parallel spikes with the HOCs, one may elucidate assembly activities and possibly their behavioral relevance. To describe the HOCs in parallel spike trains, the log-linear model is an useful model because it provides a well-defined measure of correlation based on information geometry . Former studies on HOCs performed a regression analysis on parallel spike trains using either a full log-linear model  or a log-linear model containing up to pairwise interaction only (maximum entropy model) . The existing approaches, however, assume stationarity, a condition that is typically not fulfilled in neuronal spike data from awake behaving animals. Recently, we established a method for estimating the dynamics of HOCs by means of a state-space analysis  with a log-linear observation model to trace active assemblies [Abstracts in SAND4, NIPS08 Workshop, and Cosyne09, ]. However, presentation of the smoothed posterior estimates only may mislead neurophysiologists to presume the existence of the HOC at the moment when no or weak evidence is available. Furthermore, the method did not provide a statistic to detect an assembly in which cells are jointly connected through multiple correlations. In this contribution, we investigate the method of the Bayesian hypothesis testing to answer which compositions of parallel spikes exhibit the joint HOCs, and if they do, when those HOCs appear. We computed the Bayes factor  of temporally local spike observation to gain evidence of positive joint HOCs of a specific set of parallel spike sequences against negative HOCs by using filter and one-step prediction odds. The proposed method may be useful to detect the dynamic assembly activities, their composition and behavioral relevance when applied to simultaneous recordings of neuronal activity of behaving animal.
Vaadia E, Haalman I, Abeles M, Bergman H, Prut Y, Slovin H, Aertsen A: Dynamics of neuronal interactions in monkey cortex in relation to behavioral events. Nature. 1995, 373: 515-518.
Riehle A, Grün S, Diesmann M, Aertsen A: Spike synchronization and rate modulation differentially involved in motor cortical function. Science. 1997, 278: 1950-1953.
Fujisawa S, Amarasingham A, Harrison MT, Buzsáki G: Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex. Nat Neurosci. 2008, 11: 823-833.
Grün S, Abeles M, Diesmann M: Impact of higher-order correlations on coincidence distributions of massively parallel data. Lect Notes Comput Sci. 2008, 5286: 96-114.
Amari S: Information geometry on hierarchy of probability distributions. IEEE Trans Inf Theory. 2001, 47: 1701-1711.
Martignon L, Hasseln HV, Grün S, Aertsen A, Palm G: Detecting higher-order interactions among the spiking events in a group of neurons. Biol Cybern. 1995, 73: 69-81.
Schneidman E, Berry MJ, Segev R, Bialek W: Weak pairwise correlations imply strongly correlated network states in a neural population. Nature. 2006, 440: 1007-1012.
Smith AC, Brown EN: Estimating a state-space model from point process observations. Neural Comput. 2003, 15: 965-9919.
Shimazaki H, Amari S, Brown EN, Grün S: State-space analysis on time-varying correlations in parallel spike sequences. Proc IEEE ICASSP. 2009
Kass RE, Raftery AE: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795.
This work was supported in part by JSPS Research Fellowships for Young Scientists (HS), R01 MH59733/DP1 OD 003646 (EB), and RIKEN Strategic Programs for R&D (SG).