A method for constructing data-based models of spiking neurons using a dynamic linear-static nonlinear cascade

This paper describes a general nonlinear dynamical model for neural system identification. It describes an algorithm for fitting a simple form of the model to spike train data, and reports on this algorithm's performance in identifying the structure and parameters of simulated neurons. The central element of the model is a Wiener-Bose dynamic nonlinearity that ensures that the model is able to approximate the behaviour of an arbitrary nonlinear dynamical system. Nonlinearities associated with spike generation and transmission are treated by placing the Wiener-Bose system in cascade with pulse frequency modulators and demodulators, and the static nonlinearity at the output of the Wiener-Bose system is decomposed into a rectifier and a multinomial. This simplifies the model without reducing its generality for neuronal system identification. Model elements can be characterised using standard methods of dynamical systems analysis, and the model has a simple form that can be implemented and simulated efficiently. This model bears a structural resemblance to real neurons; it may be regarded as a connectionist neuron that has been generalized in a realistic way to enable it to mimic the behaviour of an arbitrary nonlinear system, or conversely as a general nonlinear model that has been constrained to make it easy to fit to spike train data. Tests with simulated data show that the identification algorithm can accurately estimate the structure and parameters of neuron-like nonlinear dynamical systems using data sets containing only a few hundred spikes.