The Annals of Applied Statistics

Adjusted regularization in latent graphical models: Application to multiple-neuron spike count data

Giuseppe Vinci, Valérie Ventura, Matthew A. Smith, and Robert E. Kass

Full-text: Access denied (no subscription detected)

We're sorry, but we are unable to provide you with the full text of this article because we are not able to identify you as a subscriber. If you have a personal subscription to this journal, then please login. If you are already logged in, then you may need to update your profile to register your subscription. Read more about accessing full-text


A major challenge in contemporary neuroscience is to analyze data from large numbers of neurons recorded simultaneously across many experimental replications (trials), where the data are counts of neural firing events, and one of the basic problems is to characterize the dependence structure among such multivariate counts. Methods of estimating high-dimensional covariation based on $\ell_{1}$-regularization are most appropriate when there are a small number of relatively large partial correlations, but in neural data there are often large numbers of relatively small partial correlations. Furthermore, the variation across trials is often confounded by Poisson-like variation within trials. To overcome these problems we introduce a comprehensive methodology that imbeds a Gaussian graphical model into a hierarchical structure: the counts are assumed Poisson, conditionally on latent variables that follow a Gaussian graphical model, and the graphical model parameters, in turn, are assumed to depend on physiologically-motivated covariates, which can greatly improve correct detection of interactions (nonzero partial correlations). We develop a Bayesian approach to fitting this covariate-adjusted generalized graphical model and we demonstrate its success in simulation studies. We then apply it to data from an experiment on visual attention, where we assess functional interactions between neurons recorded from two brain areas.

Article information

Ann. Appl. Stat., Volume 12, Number 2 (2018), 1068-1095.

Received: December 2017
Revised: June 2018
First available in Project Euclid: 28 July 2018

Permanent link to this document

Digital Object Identifier

Bayesian inference Gaussian graphical models Gaussian scale mixture high dimensionality lasso latent variable models macaque prefrontal cortex macaque visual cortex Poisson-lognormal sparsity spike-counts


Vinci, Giuseppe; Ventura, Valérie; Smith, Matthew A.; Kass, Robert E. Adjusted regularization in latent graphical models: Application to multiple-neuron spike count data. Ann. Appl. Stat. 12 (2018), no. 2, 1068--1095. doi:10.1214/18-AOAS1190.

Export citation


  • Andrews, D. F. and Mallows, C. L. (1974). Scale mixtures of normal distributions. J. Roy. Statist. Soc. Ser. B 36 99–102.
  • Banerjee, S. and Ghosal, S. (2015). Bayesian structure learning in graphical models. J. Multivariate Anal. 136 147–162.
  • Behseta, S., Berdyyeva, T., Olson, C. R. and Kass, R. E. (2009). Bayesian correction for attenuation of correlation in multi-trial spike count data. J. Neurophysiol. 101 2186–2193.
  • Benjamini, Y. and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. Roy. Statist. Soc. Ser. B 57 289–300.
  • Benjamini, Y. and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Ann. Statist. 29 1165–1188.
  • Brillinger, D. R. (1988). Maximum likelihood analysis of spike trains of interacting nerve cells. Biol. Cybernet. 59 189–200.
  • Chandrasekaran, V., Parrilo, P. A. and Willsky, A. S. (2012). Latent variable graphical model selection via convex optimization. Ann. Statist. 40 1935–1967.
  • Churchland, A. K., Kiani, R., Chaudhuri, R., Wang, X. J., Pouget, A. and Shadlen, M. N. (2011). Variance as a signature of neural computations during decision making. Neuron 69 818–831.
  • Cohen, M. R. and Maunsell, J. H. (2009). Attention improves performance primarily by reducing interneuronal correlations. Nat. Neurosci. 12 1594–1600.
  • Cox, D. R. and Lewis, P. A. W. (1972). Multivariate point processes. 401–448.
  • Ecker, A. S., Berens, P., Cotton, R. J., Subramaniyan, M., Denfield, G. H., Cadwell, C. R., Smirnakis, S. M., Bethge, M. and Tolias, A. S. (2014). State dependence of noise correlations in macaque primary visual cortex. Neuron 82 235–248.
  • Efron, B., Tibshirani, R., Storey, J. D. and Tusher, V. (2001). Empirical Bayes analysis of a microarray experiment. J. Amer. Statist. Assoc. 96 1151–1160.
  • Fan, J., Feng, Y. and Wu, Y. (2009). Network exploration via the adaptive lasso and SCAD penalties. Ann. Appl. Stat. 3 521–541.
  • Fienberg, S. E. (1974). Stochastic models for single neuron firing trains: A survey. Biometrics 30 399–427.
  • Friedman, J., Hastie, T. and Tibshirani, R. (2008). Sparse inverse covariance estimation with the graphical lasso. Biostatistics 9 432–441.
  • George, E. I. and McCulloch, R. E. (1997). Approaches for Bayesian variable selection. Statist. Sinica 7 339–373.
  • Georgopoulos, A. P. and Ashe, J. (2000). One motor cortex, two different views. Nat. Neurosci. 3 964–963; author reply 965.
  • Giraud, C. and Tsybakov, A. (2012). Discussion: Latent variable graphical model selection via convex optimization [MR3059067]. Ann. Statist. 40 1984–1988.
  • Goris, R. L., Movshon, J. A. and Simoncelli, E. P. (2014). Partitioning neuronal variability. Nat. Neurosci. 17 858–865.
  • Hinne, M., Ambrogioni, L., Janssen, R. J., Heskes, T. and van Gerven, M. A. J. (2014). Structurally-informed Bayesian functional connectivity analysis. NeuroImage 86 294–305.
  • Hoerl, A. E. and Kennard, R. W. (1970). Ridge regression: Biased estimation for nonorthogonal problems. Technometrics 12 55–67.
  • Inouye, D. I., Yang, E., Allen, G. I. and Ravikumar, P. (2017). A review of multivariate distributions for count data derived from the Poisson distribution. Wiley Interdiscip. Rev.: Comput. Stat. 9 e1398, 25.
  • Kass, R. E., Eden, U. T. and Brown, E. N. (2014). Analysis of Neural Data. Springer, New York.
  • Kass, R. E., Ventura, V. and Brown, E. N. (2005). Statistical issues in the analysis of neuronal data. J. Neurophysiol. 94 8–25.
  • Kass, R. E., Amari, S. I., Arai, K., Brown, E. N., Diekman, C. O., Diesmann, M. …and Fukai, T. (2018). Computational neuroscience: Mathematical and statistical perspectives. Ann. Rev. Statist. Appl. 5 183–214.
  • Maunsell, J. H. (2015). Neuronal mechanisms of visual attention. Ann. Rev. Vision Sci. 1 373–391.
  • Mazumder, R. and Hastie, T. (2012). The graphical lasso: New insights and alternatives. Electron. J. Stat. 6 2125–2149.
  • Mitchell, J. F., Sundberg, K. A. and Reynolds, J. H. (2009). Spatial attention decorrelates intrinsic activity fluctuations in macaque area V4. Neuron 63 879–888.
  • Ng, B., Varoquaux, G., Poline, J. B. and Thirion, B. (2012). A novel sparse graphical approach for multimodal brain connectivity inference. In International Conference on Medical Image Computing and Computer-Assisted Intervention 707–714. Springer, Berlin, Heidelberg.
  • Perkel, D. H., Gerstein, G. L. and Moore, G. P. (1967a). Neuronal spike trains and stochastic point processes: I. The single spike train. Biophys. J. 7 391–418.
  • Perkel, D. H., Gerstein, G. L. and Moore, G. P. (1967b). Neuronal spike trains and stochastic point processes: II. Simultaneous spike trains. Biophys. J. 7 419–440.
  • Pineda-Pardo, J. A., Bruña, R., Woolrich, M., Marcos, A., Nobre, A. C., Maestú, F. and Vidaurre, D. (2014). Guiding functional connectivity estimation by structural connectivity in MEG: An application to discrimination of conditions of mild cognitive impairment. NeuroImage 101 765–777.
  • Polson, N. G., Scott, J. G. and Windle, J. (2013). Bayesian inference for logistic models using Pólya–Gamma latent variables. J. Amer. Statist. Assoc. 108 1339–1349.
  • Rabinowitz, N. C., Goris, R. L., Cohen, M. and Simoncelli, E. P. (2015). Attention stabilizes the shared gain of V4 populations. eLife 4 e08998.
  • Ravikumar, P., Wainwright, M. J., Raskutti, G. and Yu, B. (2011). High-dimensional covariance estimation by minimizing $\ell_{1}$-penalized log-determinant divergence. Electron. J. Stat. 5 935–980.
  • Rothman, A. J., Bickel, P. J., Levina, E. and Zhu, J. (2008). Sparse permutation invariant covariance estimation. Electron. J. Stat. 2 494–515.
  • Ruff, D. A. and Cohen, M. R. (2016). Stimulus dependence of correlated variability across cortical areas. J. Neurosci. 36 7546–7556.
  • Scott, J. G., Kelly, R. C., Smith, M. A., Zhou, P. and Kass, R. E. (2015). False discovery rate regression: An application to neural synchrony detection in primary visual cortex. J. Amer. Statist. Assoc. 110 459–471.
  • Smith, M. A. and Kohn, A. (2008). Spatial and temporal scales of neuronal correlation in primary visual cortex. J. Neurosci. 28 12591–12603.
  • Smith, M. A. and Sommer, M. A. (2013). Spatial and temporal scales of neuronal correlation in visual area V4. J. Neurosci. 33 5422–5432.
  • Snyder, A. C., Morais, M. J. and Smith, M. A. (2016). Dynamics of excitatory and inhibitory networks are differentially altered by selective attention. J. Neurophysiol. 116 1807–1820.
  • Tibshirani, R. (1996). Regression shrinkage and selection via the lasso. J. Roy. Statist. Soc. Ser. B 58 267–288.
  • Vinci, G., Ventura, V., Smith, M. A. and Kass, R. E. (2016). Separating spike count correlation from firing rate correlation. Neural Comput. 28 849–881.
  • Vinci, G., Ventura, V., Smith, M. A. and Kass, R. E. (2018a). Adjusted regularization of cortical covariance. J. Comput. Neurosci. To appear.
  • Vinci, G., Ventura, V., Smith, M. A. and Kass, R. E. (2018b). Supplement to “Adjusted regularization in latent graphical models: Application to multiple-neuron spike count data.” DOI:10.1214/18-AOAS1190SUPP.
  • Wang, H. (2012). Bayesian graphical lasso models and efficient posterior computation. Bayesian Anal. 7 867–886.
  • Wang, H. (2015). Scaling it up: Stochastic search structure learning in graphical models. Bayesian Anal. 10 351–377.
  • Wilson, H. R. and Cowan, J. D. (1972). Excitatory and inhibitory interactions in localized populations of model neurons. Biophys. J. 12 1–24.
  • Yatsenko, D., Josić, K., Ecker, A. S., Froudarakis, E., Cotton, R. J. and Tolias, A. S. (2015). Improved estimation and interpretation of correlations in neural circuits. PLoS Comput. Biol. 11 e1004083.
  • Yu, B. M., Cunningham, J. P., Santhanam, G., Ryu, S. I., Shenoy, K. V. and Sahani, M. (2009). Gaussian-process factor analysis for low-dimensional single-trial analysis of neural population activity. In Advances in Neural Information Processing Systems 1881–1888.
  • Yuan, M. (2012). Discussion: Latent variable graphical model selection via convex optimization [MR3059067]. Ann. Statist. 40 1968–1972.
  • Yuan, M. and Lin, Y. (2007). Model selection and estimation in the Gaussian graphical model. Biometrika 94 19–35.
  • Zou, H. and Hastie, T. (2005). Regularization and variable selection via the elastic net. J. R. Stat. Soc. Ser. B. Stat. Methodol. 67 301–320.

Supplemental materials

  • Supplement to “Adjusted regularization in latent graphical models: Application to multiple-neuron spike count data.”. Appendix containing: Appendix A additional details about adjusted regularization and simulations, Appendix B lemmas with proofs, Appendix C algorithms, and Appendix D additional data analyses.