Volume 15 Supplement 1

## Abstracts from the Twenty Third Annual Computational Neuroscience Meeting: CNS*2014

# Probability-based nonlinear modeling of neural dynamical systems with point-process inputs and outputs

- Roman Sandler
^{1}Email author, - Dong Song
^{1}, - Robert E Hampson
^{2}, - Sam A Deadwyler
^{2}, - Theodore Berger
^{1}and - Vasilis Marmarelis
^{1}

**15(Suppl 1)**:P102

https://doi.org/10.1186/1471-2202-15-S1-P102

© Sandler et al.; licensee BioMed Central Ltd. 2014

**Published: **21 July 2014

One of the great challenges of contemporary neuroscience is understanding in a quantitative manner how neurons and neural ensembles encode and process information in the form of action potentials. The advent of multi-electrode arrays which can record large amounts of data simultaneously from several neurons has made this task more urgent. This task, however, is made difficult by the inherent complexity of neural systems which are highly nonlinear, interconnected, dynamic, and subject to stochastic variations. Furthermore, while several methods exist which offer good predictive performance for spike train modeling, their adoption in the broader neuroscience community has been limited due to their mathematical complexity and lack of interpretability.

^{th}order PBV kernel, PBV

_{n}, is derived in two steps. First, we calculate the conditional probability of an output spike given n input spikes at various lags. Second, we subtract lower order effects to isolate the nth order nonlinearity. Thus, the first PBV kernel, PBV

_{1}(

*τ*), is the conditional probability of an output spike given a spike in the input spike train at time

*τ*, minus the probability that there will be a spike in the output, i.e.:

_{2}(

*τ*

_{ 1 }

*, τ*

_{ 2 }), which is the first nonlinear kernel, is the conditional probability of an output spike given a

*pair*of input spikes minus the conditional probability of an output spike given either one of those input spikes individually, i.e.:

This method may be extended to describe the contribution of n pairs spikes to the output in the form of the n^{th} PBV kernel. We show that the PBV kernels are equivalent to the Wiener kernels when the input is a Poisson process, thus placing the PBV kernels in the context of a well-established and rigorous mathematical theory [1].

The proposed PBV methodology was applied to synthetic systems where the ground truth of the model was available. The PBV kernels were found to both accurately estimate the ground truth kernels and to reproduce the given output, thus validating the method. Finally, the proposed PBV methodology was applied to real neural data derived from the CA3 and CA1 regions of the rodent hippocampus [2]. Although here ground truth was not available, the PBV kernels were able to reproduce the output as well as other models which have been validated both mathematically and in-vivo in the context of neural prosthetics [2].

## Authors’ Affiliations

## References

- Marmarelis VZ, Berger TW: General methodology for nonlinear modeling of neural systems with Poisson point-process inputs. Mathematical biosciences. 2005, 196 (1): 1-13. 10.1016/j.mbs.2005.04.002.View ArticlePubMedGoogle Scholar
- Berger , Theodore W, et al: A hippocampal cognitive prosthesis: multi-input, multi-output nonlinear modeling and VLSI implementation. Neural Systems and Rehabilitation Engineering, IEEE Transactions. 2012, 20 (2): 198-211.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.