Analysis of higher-order neuronal interactions based on conditional inference

 Higher-order neural interactions, i.e., interactions that cannot be reduced to interactions between pairs of cells, have received increasing attention in the context of recent attempts to understand the cooperative dynamics in cortical neural networks. Typically, likelihood-ratio tests of log-linear models are being employed for statistical inference. The parameter estimation of these models for simultaneously recorded single-neuron spiking activities is a crucial ingredient of this approach. Extending a previous investigation of a two-neuron system, we present here the general formulation of an exact test suited for the detection of positive higher-order interactions between m neurons. This procedure does not require the estimation of any interaction parameters and additionally optimizes the test power of the statistical inference. We apply the approach to a three-neuron system and show how second-order and third-order interactions can be reliably distinguished. We study the performance of the method as a function of the interaction strength. Received: 18 January 2002 / Accepted in revised form: 26 November 2002 / Published online: 13 March 2003 RID="*" ID="*" Present address: Institute for Theoretical Biology, Humboldt University, 10115 Berlin, Germany Correspondence to: R. Gütig (e-mail: r.guetig@biologie.hu-berlin.de, Tel.: +49 30 2093 9112, Fax: +49 30 2093 8801) Acknowledgements. We thank Shun-ichi Amari and Hiro Nakahara for valuable discussions on the information geometry of the exponential family of probability distributions underlying the present approach. Supported in part by the Studienstiftung des deutschen Volkes, the German-Israeli Foundation for Scientific Research and Development (GIF), the Deutsche Forschungsgemeinschaft (DFG), and the Institut für Grenzgebiete der Psychologie, Freiburg.