Signatures of criticality arise from random subsampling in simple population models

PLoS Comput Biol. 2017 Oct 3;13(10):e1005718. doi: 10.1371/journal.pcbi.1005718. eCollection 2017 Oct.

Abstract

The rise of large-scale recordings of neuronal activity has fueled the hope to gain new insights into the collective activity of neural ensembles. How can one link the statistics of neural population activity to underlying principles and theories? One attempt to interpret such data builds upon analogies to the behaviour of collective systems in statistical physics. Divergence of the specific heat-a measure of population statistics derived from thermodynamics-has been used to suggest that neural populations are optimized to operate at a "critical point". However, these findings have been challenged by theoretical studies which have shown that common inputs can lead to diverging specific heat. Here, we connect "signatures of criticality", and in particular the divergence of specific heat, back to statistics of neural population activity commonly studied in neural coding: firing rates and pairwise correlations. We show that the specific heat diverges whenever the average correlation strength does not depend on population size. This is necessarily true when data with correlations is randomly subsampled during the analysis process, irrespective of the detailed structure or origin of correlations. We also show how the characteristic shape of specific heat capacity curves depends on firing rates and correlations, using both analytically tractable models and numerical simulations of a canonical feed-forward population model. To analyze these simulations, we develop efficient methods for characterizing large-scale neural population activity with maximum entropy models. We find that, consistent with experimental findings, increases in firing rates and correlation directly lead to more pronounced signatures. Thus, previous reports of thermodynamical criticality in neural populations based on the analysis of specific heat can be explained by average firing rates and correlations, and are not indicative of an optimized coding strategy. We conclude that a reliable interpretation of statistical tests for theories of neural coding is possible only in reference to relevant ground-truth models.

MeSH terms

  • Animals
  • Cats
  • Computational Biology / methods*
  • Models, Neurological*
  • Neurons / physiology*
  • Retinal Ganglion Cells / cytology
  • Retinal Ganglion Cells / physiology
  • Thermodynamics

Grants and funding

Work was funded by the German Federal Ministry of Education and 568 Research (BMBF; FKZ: 01GQ1002, Bernstein Center Tübingen, FKZ 01GQ1601 to PB), the 569 German Research Foundation (BE 5601/1-1 to PB) the Max Planck Society and the caesar foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.