Abstract
Understanding how neurons encode multiple simultaneous stimuli is a fundamental question in neuroscience. We have previously introduced a novel theory of stochastic encoding patterns wherein a neuron's spiking activity dynamically switches among its constituent single-stimulus activity patterns when presented with multiple stimuli (Groh et al., 2024). Here, we present an enhanced, comprehensive statistical testing framework for such "multiplexing" or "code juggling". Our new approach evaluates whether dual-stimulus responses can be accounted for as mixtures of Poissons either anchored to or bounded by single-stimulus benchmarks. Our enhanced framework improves upon previous methods in two key ways. First, it introduces a stronger set of foils for multiplexing, including an "overreaching" category that captures overdispersed activity patterns unrelated to the single-stimulus benchmarks, reducing false detection of multiplexing/code-juggling. Second, it detects faster fluctuations - i.e. at sub-trial timescales - that would have been overlooked before. We utilize a Bayesian inference framework, considering the hypothesis with the highest posterior probability as the winner, and employ predictive recursion marginal likelihood method for the involving nonparametric density estimation. Reanalysis of previous findings confirms the general observation of "code juggling" and indicates that such juggling may well occur on faster timescales than previously suggested. We further confirm that juggling is more prevalent in (a) the inferotemporal face patch system for combinations of face stimuli than for faces and non-face objects; and (b) the primary visual cortex for distinct vs fused objects.