Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Twentieth Annual Computational Neuroscience Meeting: CNS*2011

Open Access Poster presentation

Prediction of spatiotemporal patterns of neural activity using a higher-order Markov representation of instantaneous pairwise maximum entropy model

Mehdi Aghagolzadeh1 and Karim Oweiss12*

Author Affiliations

1 Department of Electrical and Computer Engineering, Michigan State University, East Lansing, MI, 48823, USA

2 Neuroscience Program, Michigan State University, East Lansing, MI, 48823, USA

For all author emails, please log on.

BMC Neuroscience 2011, 12(Suppl 1):P245  doi:10.1186/1471-2202-12-S1-P245

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2202/12/S1/P245


Published:18 July 2011

© 2011 Aghagolzadeh and Oweiss; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Poster presentation

Maximum entropy (MaxEnt) models have been suggested to fit the spatiotemporal patterns of simultaneously observed neuronal spike trains. To build such models, the distribution, px(X), needs to be determined, where X is a sequence of neural patterns formed from the vector X1,…,Xn for n time samples, where Xi is a sample of the spike train of a population at time i . Each neural pattern itself is a binary sequence, for which one represents a spike, and zero otherwise. MaxEnt models are constructed by matching the moments between neuronal constituents with the analogous empirical estimates from the distribution px(X).

MaxEnt models that incorporate the first and second order interactions among neurons, also referred to as the instantaneous pairwise MaxEnt model or the Ising model, are expressed as

where the individual, hi, and coupling, Jij, terms are parameters associated with the first and second order interactions, respectively, and Z is a normalization factor to insure . This model was shown to be effective in representing almost 90 percent of the information in small populations (less than 20 neurons) [1,2]. The main disadvantage of this model, however, is that it ignores the temporal correlation that represents the effect of the history of spiking activity on the current population state. Extending the width of neural patterns has been proposed as a method to incorporate spiking history. It destroys, however, the information in precise spike timings [3].

To incorporate the history term without compromising information in spike times, a higher-order Markov model is constructed in which patterns from different time samples are added. An example first-order Markov representation is expressed as

where and are the individual, and , and are the coupling terms among the neuronal elements. This spatiotemporal model is more general than the representation introduced in Marre et al [4] to account for history terms , but at the expense of increasing the number of parameters.

We evaluated the performance of the extended MaxEnt model in predicting the activity pattern of cat V1 neurons in response to drifting grating stimuli. A population of 21 neurons was recorded with 5 tetrodes. We compared the prediction power of this model with the instantaneous pairwise MaxEnt model, as well as the independent model. We show that the iterative scaling algorithm makes the extended model converge faster compared to [4], and also reduces the Kullback-Leibler divergence between the estimated and true distributions by 20 percent. Taken together, these results suggest the importance of accounting for temporal correlations in predicting the spatiotemporal patterns.

References

  1. Schneidman E, Berry MJ, Segev R, Bialek W: Weak pairwise correlations imply strongly correlated network states in a neural population.

    Nature 2006, 440:1012-1007. Publisher Full Text OpenURL

  2. Tang A, Jackson D, Hobbs J, Chen W, Smith JL, Patel H, Prieto A, Petrusca D, Grivich MI, Sher A: A maximum entropy model applied to spatial and temporal correlations from cortical networks in vitro.

    Journal of Neuroscience 2008, 28:505. PubMed Abstract | Publisher Full Text OpenURL

  3. Truccolo W, Hochberg LR, Donoghue JP: Collective dynamics in human and monkey sensorimotor cortex: predicting single neuron spikes.

    Nature neuroscience 2009. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Marre O, Boustani El S, Fregnac Y, Destexhe A: Prediction of spatiotemporal patterns of neural activity from pairwise correlations.

    Physical review letters 2009, 102:138101. PubMed Abstract | Publisher Full Text OpenURL