Real and Modeled Spike Trains: Where Do They Meet?
|
|
- Clemence Johnston
- 5 years ago
- Views:
Transcription
1 Real and Modeled Spike Trains: Where Do They Meet? Vasile V. Moca 1, Danko Nikolić,3, and Raul C. Mureşan 1, 1 Center for Cognitive and Neural Studies (Coneural), Str. Cireşilor nr. 9, 4487 Cluj-Napoca, Romania moca@coneural.org, muresan@coneural.org, WWW home page: Max Planck Institute for Brain Research, Deutschordenstrasse 46, 658 Frankfurt am Main, Germany 3 Frankfurt Institute for Advanced Studies, Ruth-Moufang-Str. 1, 6438 Frankfurt am Main, Germany Abstract. Spike train models are important for the development and calibration of data analysis methods and for the quantification of certain properties of the data. We study here the properties of a spike train model that can produce both oscillatory and non-oscillatory spike trains, faithfully reproducing the firing statistics of the original spiking data being modeled. Furthermore, using data recorded from cat visual cortex, we show that despite the fact that firing statistics are reproduced, the dynamics of the modeled spike trains are significantly different from their biological counterparts. We conclude that spike train models are difficult to use when studying collective dynamics of neurons and that there is no universal recipe for modeling cortical firing, as the latter can be both very complex and highly variable. 1 Introduction The activity of cortical neurons arises from complex firing patterns that are determined by the intricate brain architecture, the external stimuli and the interaction with subcortical structures [1]. Each cortical neuron receives on the order of inputs from other neurons, and hence it is prone to being exposed to a high input bombardment []. Indeed, it has been suggested that such a bombardment exists, and, in addition to a balanced excitation-inhibition state, it keeps the neuron in a so-called high-conductance state [3]. Moreover, under such heavy input, cortical neurons have a tendency to fire highly irregularly, such that the distribution of their inter-spike intervals () takes an exponential form [4]. As a result, numerous models of spike trains have been proposed [5, 6], which assume an underlying homogenous (constant firing probability over time) or inhomogeneous (the instantaneous firing probability can fluctuate over time) Poisson process [7]. Recently however, the heavy input bombardment hypothesis has been challenged [8], new [9] and old [1] data suggesting that the assumption
2 of neurons being simple Poisson spike generators is an extreme oversimplification. The firing properties of neurons depend critically on neuron type [11], brain structure [1], brain state [13], arousal [14], and other factors. A particular case, where firing is non-poissonian, is represented by the oscillatory discharge of cortical neurons. In such cases, neurons can engage into various rhythms in multiple frequency bands [15]. The oscillatory firing properties of neurons can be characterized by computing auto-correlation histograms (), and we have previously introduced a measure, called Oscillation Score [16] that is useful in determining the degree to which a neuron oscillates. In order to precisely quantify the oscillation strength we needed to develop a model of oscillatory spike trains, in which the oscillatory behavior and the firing rate could be independently controlled [16]. Here, we explicitly study the properties of this model: the precision with which firing rate and oscillation strength can be controlled and the ability to independently express the two properties in the spike trains. In addition, we are interested in the degree to which artificial spike trains, generated from statistical parameter distributions of recorded neuronal data, reflect the temporal structure of the latter. The model The model produces artificial spike trains that retain basic properties of a recorded data set. The considered properties are: the firing rates, burst probability, oscillation strength and frequency, spike counts in bursts, refractory periods, and intra-burst inter-spike intervals, all of which are quantitatively determined from biological data. Two processes are used to generate the spike trains: at a coarse scale a discharge probability function p s (t) (Fig. 1A, B), and at a finer scale another process (Fig. 1C) that controls the exact spike timings (refractoriness, burst properties, etc). The spike discharge probability p s (t) should have the following properties: first, it should allow the spike train to exhibit a preferred oscillation frequency for transient periods of time; second, it should allow control over the strength and stability of the desired oscillation; and third, it should enable the control over the firing rates. To control the amount of oscillations, two discharge probabilities p o (t) and p b (t) corresponding respectively to an oscillatory and a background process are intermixed with a factor o (oscillation strength) (1). To obtain the transient oscillatory behavior we modulate the frequency of a sine probability function p o (t) () by a random process f o (t) (4). This random process takes into account past values and thus it has memory. The history dependence is given by a decay constant, τ, while another factor, m, controls the amount of noise added to the random process. The interplay between τ and m controls how the oscillation frequency changes over time. The frequency range of the modulatory function, f o (t), can be bounded to increase the stability of the oscillation. The background probability p b (t) (3) is generated in the same way except that in this case, the frequency of the process f b (t) (which has the same form as f o (t)) varies in a much broader frequency range (Fig. 1B). For a spike to be generated, p s (t)
3 5Hz fo( t) fb( t) p t s( ) offset 1Hz Hz Hz Time [ms] Time [s] 1 A p t b( ) p t o( ) B 1 C f t b( ) p t s( ) r b r b SpkCnt f t o( ) Fig. 1. Model for simulation of oscillatory spiking activity. A: Spike discharge probability. The spike discharge probability p s(t) is obtained by mixing two time-varying processes p b (t) and p o (t) with frequencies f b (t) and f o (t), respectively and with a mixing factor o =.3 (see Eq. 1). The amplitude and offset of the two processes are fixed for one run of the algorithm. A spike can be generated only where p s(t) is positive (grey band). B: Fluctuations of the modulation functions f b (t) and f o (t). While f b (t) varies between and 1 Hz (gray) f o(t) varies slowly around 5 Hz due to its strong history dependence and boundaries ( to 3 Hz). C: Spiking. A burst of spikes is represented by grey vertical bars, while a tonic spike is depicted in black. The spikes are generated taking into account the refractory period after a burst or a tonic spike r, the intra-burst inter-spike interval b and the burst spike count b SpkCnt. must be positive. The greater p s (t), the greater is the chance that a spike will be generated and thus, most of the spikes will be concentrated at the peaks of the probability functions (Fig. 1C). By manipulating the of f set (1) we can control how well the spikes are aligned to the desired oscillation, and by manipulating the spike-probability-positive-integral (SP P I) (5) trough the amplitude parameter, A (1), we control the firing rates (Fig. 1A and Fig. - left column) for the duration T of the spike train. Thus, the value of p s (t) controls the periodicity and firing rates of spike-trains in a manner that realistically mimics the oscillatory behavior of the recorded neurons. p s (t) =A [o p o (t) + (1 o) p b (t) offset] (1) p o (t) = t ) ( sin π f o (x)dx () f o (t) = lim ɛ p b (t) = sin ( π 3τ ɛ e x τ f o (t x)dx 3τ ɛ e x τ dx t ) f b (x)dx (3) + m rand [ 1,1] (t) (4)
4 T SP P I = 1 p + s (t)dt (5) T { where: p + ps (t), p s (t) = s (t) (6), p s (t) < At smaller time scales, the model controls the timings between the spikes and the burstiness of the spike trains. The burst probability, p burst, is modeled as a constant ratio between the count of bursts relative to the number of tonic spikes and burst occurrences altogether. These values can be extracted from a recorded data set. Once a discharge is initiated based on, p s (t), p burst determines whether that discharge will be a burst or tonic spike. If a burst is generated, the number of spikes in the burst b SpkCnt and the spacing between the spikes within the burst b (Fig. 1C) are set according to probabilities measured from real data. After a tonic spike a refractory period, r, prevents the occurrences of other discharges in a given period of time. The model makes a clear distinction between the tonic spikes and bursts. The bursts are defined as groups of spikes with successive s smaller or equal to 8 ms [17]. Thus, the refractory period, r, is set to 8 ms such that tonic spikes can occur with an of at least 9 ms. The control over these timings at a small scale produces a realistic center-peak shape in the (Fig. 3) of artificial spike trains which model real spike trains. 3 Results We wanted to address here two important aspects related to the model: the relationship between firing rate and oscillation strength, on one hand, and the degree to which the model can reproduce the fine temporal structure of recorded spike trains, on the other. 3.1 Firing rate and oscillation strength We computed the relationship between firing rate and oscillation strength by producing artificial spike trains, with a length of 3 s, and controlling independently the two parameters for rate and oscillation strength. The rate parameter (SP P I; see Methods) and the oscillation strength parameter (o) were varied in 4 steps. For each combination, two spike trains were produced, yielding a total of 115 spike trains. In addition, we generated such spike trains for each oscillation frequency band separately: theta (4-8 Hz), alpha (8-1 Hz), beta-low (1- Hz), beta-high (-3 Hz) and gamma (3-8 Hz) [16]. For the oscillation frequency in the model (f o ), we took the central (middle) frequency for each frequency band. The firing rate and the oscillation score [16] (the measured strength of oscillation) were measured for each spike train. We found that, the firing rate could be more precisely controlled for higher than for lower frequency bands (Fig., left column, note the lower variance for higher frequency bands) while the oscillation score seemed to be controlled equally well across all frequency bands (Fig., middle column). Moreover, the firing rate scaled as
5 a power function (R >.93) with respect to the SP P I while the oscillation score scaled as a sigmoid with respect to the oscillation strength parameter. The noisy clouds in Fig. (middle column), that deviate from the sigmoid shape are due to trials with very small firing rates, for which the oscillation score estimate becomes imprecise [16]. Finally, there was no correlation between the firing rate and oscillation score, as shown in Fig., right column (R <.8). Thus, the resulting firing rate and oscillation strength can be independently controlled by the two corresponding parameters in the model. 3. Temporal structure of spike trains Although our model is quite detailed, including statistics of bursting, oscillatory modulation, and other spike train properties, we wanted to estimate how closely some model spike trains resemble their corresponding, real, spike trains recorded from cat visual cortex. We considered data recorded from an anaesthetized cat on two experimental sessions: one without or poor oscillatory responses (col11b44; see Fig. 3) and one with strong oscillations (col11b68; see Fig. 3) in the betahigh frequency band. Since the experimental trials were rather short (6.5 s) we selected, from one single electrode, multi-unit activity (MUA) such as to yield sufficient spikes when computing single-trial statistics (see next). We also chose a stimulation condition (center-surround stimulus, with a small sinusoidal grating placed in the center and a larger one in the surround) for which the MUA showed strong oscillatory behavior in col11b68 session but not in col11b44. For this stimulation condition, we obtained experimental trials in both sessions. Next, for each experimental trial, one corresponding model trial (spike train) was generated. We computed the distributions of spikes within bursts (burst ), the burst spike-count distributions, and the firing rate on the experimental trial, and plugged these parameters into the model, producing an artificial spike train having the same statistics (see Methods). Importantly, we wanted to model each trial independently such that the temporal structure of firing is maximally similar to the corresponding recorded trial. Statistics are thus computed separately for each trial. For oscillatory spike trains, the oscillation frequency was computed on the whole stimulation condition and plugged into each model trial (because the oscillation frequency was very stable across trials in the analyzed dataset). Next, the autocorrelation histograms () per stimulation condition were computed by averaging the s computed on each of the trials. The inter-spike interval () distributions were also computed. In Figs. 3A and 3B we show that the s and distributions of the real and model spike trains were remarkably similar, both for the non-oscillatory and oscillatory case. This indicates that statistically, the real and modeled spike trains had similar local structure. At a first look, one could say that the model spike train can successfully replace the real one. To get further insight into this hypothesis, we computed the cross-correlation histograms (CCH) between the model and their corresponding, real, spike trains. Ideally, if the model successfully replaces the real spike train, the CCH should be very similar to the s of both the real and modeled spike trains. This was, however, not the case. For
6 Theta (4-8Hz) Alpha (8-1Hz) Beta-low (1-Hz) Beta-high (-3Hz) Gamma (3-8Hz) 1 5 Pow. fit R = SPPI Pow. fit R = SPPI Pow. fit R = SPPI Pow. fit R = SPPI Pow. fit R = SPPI Osc. Strength Osc. Strength Osc. Strength Osc. Strength Osc. Strength 1 5 Lin. fit R = Lin. fit R = Lin. fit R = Lin. fit R = Lin. fit R = Fig.. Control of rate and oscillation strength. On the left column, scatter plots show the dependence of firing on the spike-probability-positive-integral: SP P I (5). The black line shows a power fit. The center column depicts the dependency of the measured oscillation score on the oscillation strength parameter, o, in the model (1). In the right column, scatter plots of the rate and oscillation score show no dependency between these two, as indicated by linear fits (black lines).
7 A col11b44 MU N1 C col11b44 MU N1 C col11b44- CCH - B col11b68mu N1 C col11b68 MU N1 C #/Trial col11b68- CCH - C col11b68 MU N1 C Stim. On Time [s] Stim. Off Time [s] col11b68mu N1 C9 T #/Trial #/Trial col11b68mu N1 C9 T
8 Fig. 3. Comparison between recorded and spike trains. A: Non oscillatory data. The first row shows the (left panel) and the distribution (center panel) for non-oscillatory data. On the second row, the (left panel) and for data (center panel) are shown. Note that the and plots are averages across trials of the same condition. On the right panel the average cross-correlation between each recorded and artificial spike train is shown. B: The same plots as in A are shown for oscillatory data. The similarity of s and distributions between the recorded and data is again clear. C: Differences between the recorded and spike trains. The spikes are shown on the left: recorded data (upper traces, consecutive trials) and data (lower traces). For one recorded spike train and for its artificial model (both depicted in black) the s (right panels) and distributions (bottom panels) are shown. For each plot depicting recorded data the dataset, neuron, and condition labels are shown above the graph; for data simulation is shown. non-oscillatory spike trains (Fig. 3A) the CCH was flat, without a peak, indicating that there was no consistent relationship (fine temporal correlation) between real and modeled spike timings. The oscillatory case was somewhat better, such that weak oscillatory modulation could be seen in the CCH (Fig. 3B), but the central peak was missing here as well. These results suggest that, although the statistical properties of the original spike trains are faithfully reproduced, the temporal dynamics of the trials are not reproduced. We show in Fig. 3C that, for a given spike train and its model, the s and distributions had reasonably close structure, but the location of spikes within the real and model trials was strikingly different. In the real spike train there was a clear modulation of the firing by the drifting sinusoidal grating stimulus. In the modeled spike train, the firing events were more or less arbitrarily located, despite the fact that global statistics and s structure were accurately reproduced (Fig. 3C). 4 Discussion We have shown that the proposed spike train model allows one to independently control firing rate and oscillation strength. Thus, the spike train model can isolate reasonably well the two processes such that one can use it to calibrate data analysis methods which need to separate the effects of the firing rate from those of the oscillation strength [16]. Moreover, our results indicate that the model can reproduce the statistical properties of the real data quite faithfully. The s and the distributions of the model were strikingly similar to those of the real data and thus, we validated the model as being appropriate for studying properties of spike trains in terms of their s. On the other hand, we have found that reproducing the statistics of neuronal firing can not account for the dynamics of the spike trains. There are multiple conclusions that stem from the above. First, given the complexity required by our model in order to reproduce realistic spike trains, we infer that simple, homogenous or inhomogeneous Poisson processes are crude and largely
9 inappropriate approximations of cortical firing. These two types of processes need to be complemented by other, more complex ones, as real spike trains are very hard to model. One needs to add many constraints to the model and these constraints vary across neurons, stimuli, cortical states, and so forth. There is no universal recipe for generating spike-trains, since neuronal behavior can be very rich and highly variable. Second, modeling the fine temporal structure of spike trains, including their temporal dynamics can become very difficult because one needs to know the underlying drive received by neurons from the sensory (thalamic) and cortical inputs. An important implication of this fact is that collective coding strategies (where multiple, simultaneous spike trains are analyzed) cannot be easily studied with model spike trains. The reason is that, for collective codes, the relative dynamics of different neurons plays a crucial role. Reproducing the statistical properties of each individual spike train is obviously not enough. One could use a hybrid approach and also measure some instantaneous firing probability over time, but it is doubtful that the exact spike times can be accurately reproduced. Some have been able to achieve accurate modeling of firing processes, but only when the visual stimulus is known and mostly for early processing stages, such as the retina or LGN [18]. Finally, we want to emphasize that it is always important to properly assess the usefulness of a given spike train model. Models can prove useful for developing and calibrating data analysis methods [16], or for studying and quantifying given properties of the data [18]. In general however, the spiking behavior of neurons is both complex and variable. One needs to judge carefully which model captures the interesting properties relevant for a given scientific question. Furthermore, one has to be aware that, even if statistical properties of neuronal firing are precisely reproduced, the exact spiking dynamics stem from the complex interactions with other neurons and external stimuli. When one needs to generate a spike train that is very close to the original, the best model for the spike train is probably the spike train itself. Appendix To generate the data for Fig. the following parameters were used: for f o (t), τ =.5 1, m =.5 and f o (t) was allowed to deviate from the target oscillation frequency with at most 1 Hz; for f b (t), τ =.1, m = 3 and f b (t) was bounded between and 1 Hz. Note that τ is given relative to the period of the target oscillation frequency. The oscillation strength, o, took values from. to.8 in 4 steps, SP P I was varied between.4 to.1 in 4 steps and the offset was set to.5. To generate the artificial data presented in Fig. 3, some of the parameters given above had to be adapted to each recorded spike train (trial) that was modeled.
10 Acknowledgments The authors acknowledge three grants of the Romanian Government (Human Resources Program RP-5/7 contract 1/1.1.7 and Ideas Program ID 48/7 contract 4/1.1.7 both financed by MECT/UEFISCSU, and Partnerships Program contract 1139/ financed by MECT/ANCS), a grant for the Max Planck - Coneural Partner Group, and a Deutsche Forschungsgemeinschaft grant (NI 78/-1). References 1. Braitenberg, V., Schz, A.: Cortex: Statistics and Geometry of Neuronal Connectivity. nd edn, Springer, Berlin (1998). Churchland, Patricia S., Sejnowski, T.J.: The Computational Brain. MIT Press, Cambridge MA (1999) 3. Destexhe, A., Rudolph, M., Par, D.: The high-conductance state of neocortical neurons in vivo. Nat. Rev. Neurosci. 4 (3) Tuckwell, H.C.: Stochastic Processes in the Neurosciences. Society for Industrial and Applied Mathematics, Philadelphia PA (1989) 5. Rieke, F., Warland, D., de Ruyter van Steveninck, R. and Bialek, W.: Spikes: Exploring the Neural Code. MIT Press, Cambridge MA (1997) 6. Dayan, P., Abbott, L.F.: Theoretical Neuroscience. MIT Press, Cambridge MA (1) 7. Heeger, D.: Poisson model of spike generation. david/handouts /poisson.pdf () 8. Waters, J., Helmchen F.: Background synaptic activity is sparse in neocortex. J. Neurosci. 6(3) (6) Tiesinga, P., Fellous, J.M., Sejnowski, T.J.: Regulation of spike timing in visual cortical circuits. Nat. Rev. Neurosci. 9 (8) Evarts, E.V.: Temporal patterns of discharge of pyramidal tract neurons during sleep and waking in the monkey. J. Neurophysiol. 7 (1964) Soltesz, I.: Diversity in the neuronal machine. Oxford University Press, New York (5) 1. Kara, P., Reinagel, P., Reid, R.C.: Low response variability in simultaneously recorded retinal, thalamic, and cortical neurons. Neuron 7 () Steriade, M., McCormick, D.A., Sejnowski, T.J.: Thalamocortical oscillations in the sleeping and aroused brain. Science 6 (1993) Steriade, M., Timofeev, I., Drmller, N., Grenier, F.: Dynamic properties of cortico-thalamic neurons and local cortical interneurons generating fast rhythmic (3-4 Hz) spike-bursts. J. Neurophysiol. 79 (1998) Buzsáki, G.: Rhythms of the brain. Oxford University Press, Oxford (6) 16. Mureşan, R.C., Jurjuţ, O.F., Moca, V.V., Singer, W., Nikolić, D.: The Oscillation Score: An Efficient Method for Estimating Oscillation Strength in Neuronal Activity. J. Neurophysiol. 99 (8) DeBusk, B.C., DeBruyn, E.J., Snider, R.K., Kabara, J.F., Bonds, A.B.: Stimulus-dependent modulation of spike burst length in cat striate cortical cells. J. Neurophysiol. 78 (1997) Keat, J., Reinagel, P., Reid, R.C., Meister, M.: Predicting Every Spike: A Model for the Responses of Visual Neurons. Neuron 3 (1)
Consider the following spike trains from two different neurons N1 and N2:
About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in
More informationThe Spike Response Model: A Framework to Predict Neuronal Spike Trains
The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology
More informationThe homogeneous Poisson process
The homogeneous Poisson process during very short time interval Δt there is a fixed probability of an event (spike) occurring independent of what happened previously if r is the rate of the Poisson process,
More informationSpike-Frequency Adaptation: Phenomenological Model and Experimental Tests
Spike-Frequency Adaptation: Phenomenological Model and Experimental Tests J. Benda, M. Bethge, M. Hennig, K. Pawelzik & A.V.M. Herz February, 7 Abstract Spike-frequency adaptation is a common feature of
More informationSPATIOTEMPORAL ANALYSIS OF SYNCHRONIZATION OF NEURAL ENSEMBLES FOR SPATIAL DISCRIMINATIONS IN CAT STRIATE CORTEX
SPATIOTEMPORAL ANALYSIS OF SYNCHRONIZATION OF NEURAL ENSEMBLES FOR SPATIAL DISCRIMINATIONS IN CAT STRIATE CORTEX Jason M Samonds Department of Biomedical Engineering Vanderbilt University The auto-content
More informationExercises. Chapter 1. of τ approx that produces the most accurate estimate for this firing pattern.
1 Exercises Chapter 1 1. Generate spike sequences with a constant firing rate r 0 using a Poisson spike generator. Then, add a refractory period to the model by allowing the firing rate r(t) to depend
More informationHow to read a burst duration code
Neurocomputing 58 60 (2004) 1 6 www.elsevier.com/locate/neucom How to read a burst duration code Adam Kepecs a;, John Lisman b a Cold Spring Harbor Laboratory, Marks Building, 1 Bungtown Road, Cold Spring
More informationModel neurons!!poisson neurons!
Model neurons!!poisson neurons! Suggested reading:! Chapter 1.4 in Dayan, P. & Abbott, L., heoretical Neuroscience, MI Press, 2001.! Model neurons: Poisson neurons! Contents: Probability of a spike sequence
More informationA Mass Model of Interconnected Thalamic Populations including both Tonic and Burst Firing Mechanisms
International Journal of Bioelectromagnetism Vol. 12, No. 1, pp. 26-31, 2010 www.ijbem.org A Mass Model of Interconnected Thalamic Populations including both Tonic and Burst Firing Mechanisms Pirini M.
More informationencoding and estimation bottleneck and limits to visual fidelity
Retina Light Optic Nerve photoreceptors encoding and estimation bottleneck and limits to visual fidelity interneurons ganglion cells light The Neural Coding Problem s(t) {t i } Central goals for today:
More informationSean Escola. Center for Theoretical Neuroscience
Employing hidden Markov models of neural spike-trains toward the improved estimation of linear receptive fields and the decoding of multiple firing regimes Sean Escola Center for Theoretical Neuroscience
More informationNature Neuroscience: doi: /nn.2283
Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S. Overall distribution of
More informationCORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE. by Pamela Reitsma. B.S., University of Maine, 2007
CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE by Pamela Reitsma B.S., University of Maine, 27 Submitted to the Graduate Faculty of the Department of Mathematics in partial
More informationEmpirical Bayes interpretations of random point events
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 (25) L531 L537 doi:1.188/35-447/38/29/l4 LETTER TO THE EDITOR Empirical Bayes interpretations of
More informationDesign principles for contrast gain control from an information theoretic perspective
Design principles for contrast gain control from an information theoretic perspective Yuguo Yu Brian Potetz Tai Sing Lee Center for the Neural Computer Science Computer Science Basis of Cognition Department
More informationMid Year Project Report: Statistical models of visual neurons
Mid Year Project Report: Statistical models of visual neurons Anna Sotnikova asotniko@math.umd.edu Project Advisor: Prof. Daniel A. Butts dab@umd.edu Department of Biology Abstract Studying visual neurons
More informationAn Introductory Course in Computational Neuroscience
An Introductory Course in Computational Neuroscience Contents Series Foreword Acknowledgments Preface 1 Preliminary Material 1.1. Introduction 1.1.1 The Cell, the Circuit, and the Brain 1.1.2 Physics of
More information1/12/2017. Computational neuroscience. Neurotechnology.
Computational neuroscience Neurotechnology https://devblogs.nvidia.com/parallelforall/deep-learning-nutshell-core-concepts/ 1 Neurotechnology http://www.lce.hut.fi/research/cogntech/neurophysiology Recording
More informationA MODEL OF A NEURONAL STRUCTURE ABLE TO GENERATE COLLECTIVE OSCILLATIONS SIMILAR TO HUMAN PARIETAL ALPHA RHYTHM
MODEL OF NEURONL STRUCTURE LE TO GENERTE COLLECTIVE OSCILLTIONS SIMILR TO HUMN PRIETL LPH RHYTHM Otilia Păduraru *, Hariton Costin *+ * Institute for Theoretical Computer Science, Romanian cademy, Iaşi
More informationSpike Count Correlation Increases with Length of Time Interval in the Presence of Trial-to-Trial Variation
NOTE Communicated by Jonathan Victor Spike Count Correlation Increases with Length of Time Interval in the Presence of Trial-to-Trial Variation Robert E. Kass kass@stat.cmu.edu Valérie Ventura vventura@stat.cmu.edu
More informationNeuronal Dynamics: Computational Neuroscience of Single Neurons
Week 5 part 3a :Three definitions of rate code Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 5 Variability and Noise: The question of the neural code Wulfram Gerstner EPFL, Lausanne,
More informationComparing integrate-and-fire models estimated using intracellular and extracellular data 1
Comparing integrate-and-fire models estimated using intracellular and extracellular data 1 Liam Paninski a,b,2 Jonathan Pillow b Eero Simoncelli b a Gatsby Computational Neuroscience Unit, University College
More information!) + log(t) # n i. The last two terms on the right hand side (RHS) are clearly independent of θ and can be
Supplementary Materials General case: computing log likelihood We first describe the general case of computing the log likelihood of a sensory parameter θ that is encoded by the activity of neurons. Each
More informationEmergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity
Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity Jorge F. Mejias 1,2 and Joaquín J. Torres 2 1 Department of Physics and Center for
More informationLateral organization & computation
Lateral organization & computation review Population encoding & decoding lateral organization Efficient representations that reduce or exploit redundancy Fixation task 1rst order Retinotopic maps Log-polar
More informationCHARACTERIZATION OF NONLINEAR NEURON RESPONSES
CHARACTERIZATION OF NONLINEAR NEURON RESPONSES Matt Whiteway whit8022@umd.edu Dr. Daniel A. Butts dab@umd.edu Neuroscience and Cognitive Science (NACS) Applied Mathematics and Scientific Computation (AMSC)
More informationEffects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell
Effects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell 1. Abstract Matthew Dunlevie Clement Lee Indrani Mikkilineni mdunlevi@ucsd.edu cll008@ucsd.edu imikkili@ucsd.edu Isolated
More informationarxiv: v1 [q-bio.nc] 9 Mar 2016
Temporal code versus rate code for binary Information Sources Agnieszka Pregowska a, Janusz Szczepanski a,, Eligiusz Wajnryb a arxiv:1603.02798v1 [q-bio.nc] 9 Mar 2016 a Institute of Fundamental Technological
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationAdaptation in the Neural Code of the Retina
Adaptation in the Neural Code of the Retina Lens Retina Fovea Optic Nerve Optic Nerve Bottleneck Neurons Information Receptors: 108 95% Optic Nerve 106 5% After Polyak 1941 Visual Cortex ~1010 Mean Intensity
More informationSUPPLEMENTARY INFORMATION
Supplementary discussion 1: Most excitatory and suppressive stimuli for model neurons The model allows us to determine, for each model neuron, the set of most excitatory and suppresive features. First,
More informationBayesian probability theory and generative models
Bayesian probability theory and generative models Bruno A. Olshausen November 8, 2006 Abstract Bayesian probability theory provides a mathematical framework for peforming inference, or reasoning, using
More informationWhat do V1 neurons like about natural signals
What do V1 neurons like about natural signals Yuguo Yu Richard Romero Tai Sing Lee Center for the Neural Computer Science Computer Science Basis of Cognition Department Department Carnegie Mellon University
More informationNeural variability and Poisson statistics
Neural variability and Poisson statistics January 15, 2014 1 Introduction We are in the process of deriving the Hodgkin-Huxley model. That model describes how an action potential is generated by ion specic
More informationNonlinear reverse-correlation with synthesized naturalistic noise
Cognitive Science Online, Vol1, pp1 7, 2003 http://cogsci-onlineucsdedu Nonlinear reverse-correlation with synthesized naturalistic noise Hsin-Hao Yu Department of Cognitive Science University of California
More informationFinding a Basis for the Neural State
Finding a Basis for the Neural State Chris Cueva ccueva@stanford.edu I. INTRODUCTION How is information represented in the brain? For example, consider arm movement. Neurons in dorsal premotor cortex (PMd)
More informationEstimation of information-theoretic quantities
Estimation of information-theoretic quantities Liam Paninski Gatsby Computational Neuroscience Unit University College London http://www.gatsby.ucl.ac.uk/ liam liam@gatsby.ucl.ac.uk November 16, 2004 Some
More informationEffects of synaptic conductance on the voltage distribution and firing rate of spiking neurons
PHYSICAL REVIEW E 69, 051918 (2004) Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons Magnus J. E. Richardson* Laboratory of Computational Neuroscience, Brain
More informationWhen do Correlations Increase with Firing Rates? Abstract. Author Summary. Andrea K. Barreiro 1* and Cheng Ly 2
When do Correlations Increase with Firing Rates? Andrea K. Barreiro 1* and Cheng Ly 2 1 Department of Mathematics, Southern Methodist University, Dallas, TX 75275 U.S.A. 2 Department of Statistical Sciences
More informationInformation Theory. Mark van Rossum. January 24, School of Informatics, University of Edinburgh 1 / 35
1 / 35 Information Theory Mark van Rossum School of Informatics, University of Edinburgh January 24, 2018 0 Version: January 24, 2018 Why information theory 2 / 35 Understanding the neural code. Encoding
More informationDivisive Inhibition in Recurrent Networks
Divisive Inhibition in Recurrent Networks Frances S. Chance and L. F. Abbott Volen Center for Complex Systems and Department of Biology Brandeis University Waltham MA 2454-911 Abstract Models of visual
More informationNeural Networks 1 Synchronization in Spiking Neural Networks
CS 790R Seminar Modeling & Simulation Neural Networks 1 Synchronization in Spiking Neural Networks René Doursat Department of Computer Science & Engineering University of Nevada, Reno Spring 2006 Synchronization
More informationSynchrony in Neural Systems: a very brief, biased, basic view
Synchrony in Neural Systems: a very brief, biased, basic view Tim Lewis UC Davis NIMBIOS Workshop on Synchrony April 11, 2011 components of neuronal networks neurons synapses connectivity cell type - intrinsic
More informationCausality and communities in neural networks
Causality and communities in neural networks Leonardo Angelini, Daniele Marinazzo, Mario Pellicoro, Sebastiano Stramaglia TIRES-Center for Signal Detection and Processing - Università di Bari, Bari, Italy
More informationLinearization of F-I Curves by Adaptation
LETTER Communicated by Laurence Abbott Linearization of F-I Curves by Adaptation Bard Ermentrout Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A. We show that negative
More informationHierarchy. Will Penny. 24th March Hierarchy. Will Penny. Linear Models. Convergence. Nonlinear Models. References
24th March 2011 Update Hierarchical Model Rao and Ballard (1999) presented a hierarchical model of visual cortex to show how classical and extra-classical Receptive Field (RF) effects could be explained
More informationMembrane equation. VCl. dv dt + V = V Na G Na + V K G K + V Cl G Cl. G total. C m. G total = G Na + G K + G Cl
Spiking neurons Membrane equation V GNa GK GCl Cm VNa VK VCl dv dt + V = V Na G Na + V K G K + V Cl G Cl G total G total = G Na + G K + G Cl = C m G total Membrane with synaptic inputs V Gleak GNa GK
More informationHigh-conductance states in a mean-eld cortical network model
Neurocomputing 58 60 (2004) 935 940 www.elsevier.com/locate/neucom High-conductance states in a mean-eld cortical network model Alexander Lerchner a;, Mandana Ahmadi b, John Hertz b a Oersted-DTU, Technical
More informationBursts generate a non-reducible spike-pattern code
FOCUSED REVIEW published: 1 May 2009 doi: 10.3389/neuro.01.002.2009 Bursts generate a non-reducible spike-pattern code Hugo G. Eyherabide 1,2, Ariel Rokem 1,3, Andreas V. M. Herz 1,4 and Inés Samengo 1,2
More informationTHE retina in general consists of three layers: photoreceptors
CS229 MACHINE LEARNING, STANFORD UNIVERSITY, DECEMBER 2016 1 Models of Neuron Coding in Retinal Ganglion Cells and Clustering by Receptive Field Kevin Fegelis, SUID: 005996192, Claire Hebert, SUID: 006122438,
More informationNeural Encoding: Firing Rates and Spike Statistics
Neural Encoding: Firing Rates and Spike Statistics Dayan and Abbott (21) Chapter 1 Instructor: Yoonsuck Choe; CPSC 644 Cortical Networks Background: Dirac δ Function Dirac δ function has the following
More informationVisual motion processing and perceptual decision making
Visual motion processing and perceptual decision making Aziz Hurzook (ahurzook@uwaterloo.ca) Oliver Trujillo (otrujill@uwaterloo.ca) Chris Eliasmith (celiasmith@uwaterloo.ca) Centre for Theoretical Neuroscience,
More informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More informationAdaptive contrast gain control and information maximization $
Neurocomputing 65 66 (2005) 6 www.elsevier.com/locate/neucom Adaptive contrast gain control and information maximization $ Yuguo Yu a,, Tai Sing Lee b a Center for the Neural Basis of Cognition, Carnegie
More informationThe Bayesian Brain. Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester. May 11, 2017
The Bayesian Brain Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester May 11, 2017 Bayesian Brain How do neurons represent the states of the world? How do neurons represent
More informationarxiv:physics/ v1 [physics.data-an] 7 Jun 2003
Entropy and information in neural spike trains: Progress on the sampling problem arxiv:physics/0306063v1 [physics.data-an] 7 Jun 2003 Ilya Nemenman, 1 William Bialek, 2 and Rob de Ruyter van Steveninck
More informationSelf-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire Neuron Model Based on Small World Networks
Commun. Theor. Phys. (Beijing, China) 43 (2005) pp. 466 470 c International Academic Publishers Vol. 43, No. 3, March 15, 2005 Self-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire
More informationBroadband coding with dynamic synapses
Broadband coding with dynamic synapses Benjamin Lindner Max-Planck-Institute for the Physics of Complex Systems, Nöthnitzer Str. 38 1187 Dresden, Germany André Longtin Department of Physics and Center
More informationNeural characterization in partially observed populations of spiking neurons
Presented at NIPS 2007 To appear in Adv Neural Information Processing Systems 20, Jun 2008 Neural characterization in partially observed populations of spiking neurons Jonathan W. Pillow Peter Latham Gatsby
More informationOn similarity measures for spike trains
On similarity measures for spike trains Justin Dauwels, François Vialatte, Theophane Weber, and Andrzej Cichocki RIKEN Brain Science Institute, Saitama, Japan Massachusetts Institute of Technology, Cambridge,
More informationStatistical models for neural encoding
Statistical models for neural encoding Part 1: discrete-time models Liam Paninski Gatsby Computational Neuroscience Unit University College London http://www.gatsby.ucl.ac.uk/ liam liam@gatsby.ucl.ac.uk
More informationWhen is an Integrate-and-fire Neuron like a Poisson Neuron?
When is an Integrate-and-fire Neuron like a Poisson Neuron? Charles F. Stevens Salk Institute MNL/S La Jolla, CA 92037 cfs@salk.edu Anthony Zador Salk Institute MNL/S La Jolla, CA 92037 zador@salk.edu
More informationFitting a Stochastic Neural Network Model to Real Data
Fitting a Stochastic Neural Network Model to Real Data Christophe Pouzat, Ludmila Brochini, Pierre Hodara and Guilherme Ost MAP5 Univ. Paris-Descartes and CNRS Neuromat, USP christophe.pouzat@parisdescartes.fr
More informationPatterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops
Math. Model. Nat. Phenom. Vol. 5, No. 2, 2010, pp. 67-99 DOI: 10.1051/mmnp/20105203 Patterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops J. Ma 1 and J. Wu 2 1 Department of
More informationGamma Oscillation Maintains Stimulus Structure-Dependent Synchronization in Cat Visual Cortex
J Neurophysiol 93: 223 236, 2005. First published July 28, 2004; doi:10.1152/jn.00548.2004. Gamma Oscillation Maintains Stimulus Structure-Dependent Synchronization in Cat Visual Cortex Jason M. Samonds
More informationTransformation of stimulus correlations by the retina
Transformation of stimulus correlations by the retina Kristina Simmons (University of Pennsylvania) and Jason Prentice, (now Princeton University) with Gasper Tkacik (IST Austria) Jan Homann (now Princeton
More informationNeural Encoding I: Firing Rates and Spike Statistics
Chapter 1 Neural Encoding I: Firing Rates and Spike Statistics 1.1 Introduction Neurons are remarkable among the cells of the body in their ability to propagate signals rapidly over large distances. They
More informationThe Mixed States of Associative Memories Realize Unimodal Distribution of Dominance Durations in Multistable Perception
The Mixed States of Associative Memories Realize Unimodal Distribution of Dominance Durations in Multistable Perception Takashi Kanamaru Department of Mechanical Science and ngineering, School of Advanced
More informationIntroduction to neural spike train data for phase-amplitude analysis
Electronic Journal of Statistics Vol. 8 (24) 759 768 ISSN: 935-7524 DOI:.24/4-EJS865 Introduction to neural spike train data for phase-amplitude analysis Wei Wu Department of Statistics, Florida State
More informationCoarse-grained event tree analysis for quantifying Hodgkin-Huxley neuronal network dynamics
J Comput Neurosci (212) 32:55 72 DOI 1.17/s1827-11-339-7 Coarse-grained event tree analysis for quantifying Hodgkin-Huxley neuronal network dynamics Yi Sun Aaditya V. Rangan Douglas Zhou David Cai Received:
More informationSubthreshold cross-correlations between cortical neurons: Areference model with static synapses
Neurocomputing 65 66 (25) 685 69 www.elsevier.com/locate/neucom Subthreshold cross-correlations between cortical neurons: Areference model with static synapses Ofer Melamed a,b, Gilad Silberberg b, Henry
More informationCharacterizing the firing properties of an adaptive
Characterizing the firing properties of an adaptive analog VLSI neuron Daniel Ben Dayan Rubin 1,2, Elisabetta Chicca 2, and Giacomo Indiveri 2 1 Department of Bioengineering, Politecnico di Milano, P.zza
More informationSPIKE TRIGGERED APPROACHES. Odelia Schwartz Computational Neuroscience Course 2017
SPIKE TRIGGERED APPROACHES Odelia Schwartz Computational Neuroscience Course 2017 LINEAR NONLINEAR MODELS Linear Nonlinear o Often constrain to some form of Linear, Nonlinear computations, e.g. visual
More informationAdaptive changes of rhythmic EEG oscillations in space implications for brain-machine interface applications.
ESA 29 Adaptive changes of rhythmic EEG oscillations in space implications for brain-machine interface applications. G. Cheron, A.M. Cebolla, M. Petieau, A. Bengoetxea, E. Palmero-Soler, A. Leroy, B. Dan
More informationMaximum Likelihood Estimation of a Stochastic Integrate-and-Fire Neural Model
Maximum Likelihood Estimation of a Stochastic Integrate-and-Fire Neural Model Jonathan W. Pillow, Liam Paninski, and Eero P. Simoncelli Howard Hughes Medical Institute Center for Neural Science New York
More informationConductance-Based Integrate-and-Fire Models
NOTE Communicated by Michael Hines Conductance-Based Integrate-and-Fire Models Alain Destexhe Department of Physiology, Laval University School of Medicine, Québec, G1K 7P4, Canada A conductance-based
More informationSpike-Frequency Adaptation of a Generalized Leaky Integrate-and-Fire Model Neuron
Journal of Computational Neuroscience 10, 25 45, 2001 c 2001 Kluwer Academic Publishers. Manufactured in The Netherlands. Spike-Frequency Adaptation of a Generalized Leaky Integrate-and-Fire Model Neuron
More informationProbabilistic Models in Theoretical Neuroscience
Probabilistic Models in Theoretical Neuroscience visible unit Boltzmann machine semi-restricted Boltzmann machine restricted Boltzmann machine hidden unit Neural models of probabilistic sampling: introduction
More informationActivity Driven Adaptive Stochastic. Resonance. Gregor Wenning and Klaus Obermayer. Technical University of Berlin.
Activity Driven Adaptive Stochastic Resonance Gregor Wenning and Klaus Obermayer Department of Electrical Engineering and Computer Science Technical University of Berlin Franklinstr. 8/9, 187 Berlin fgrewe,obyg@cs.tu-berlin.de
More informationDimensionality reduction in neural models: an information-theoretic generalization of spiketriggered average and covariance analysis
to appear: Journal of Vision, 26 Dimensionality reduction in neural models: an information-theoretic generalization of spiketriggered average and covariance analysis Jonathan W. Pillow 1 and Eero P. Simoncelli
More informationSeveral ways to solve the MSO problem
Several ways to solve the MSO problem J. J. Steil - Bielefeld University - Neuroinformatics Group P.O.-Box 0 0 3, D-3350 Bielefeld - Germany Abstract. The so called MSO-problem, a simple superposition
More informationPatterns of Synchrony in Neural Networks with Spike Adaptation
Patterns of Synchrony in Neural Networks with Spike Adaptation C. van Vreeswijky and D. Hanselz y yracah Institute of Physics and Center for Neural Computation, Hebrew University, Jerusalem, 9194 Israel
More informationSpatiotemporal Structure in Large Neuronal Networks Detected from Cross-Correlation
LETTER Communicated by Peter König Spatiotemporal Structure in Large Neuronal Networks Detected from Cross-Correlation Gaby Schneider gaby.schneider@math.uni-frankfurt.de Department of Computer Science
More informationNeural Coding: Integrate-and-Fire Models of Single and Multi-Neuron Responses
Neural Coding: Integrate-and-Fire Models of Single and Multi-Neuron Responses Jonathan Pillow HHMI and NYU http://www.cns.nyu.edu/~pillow Oct 5, Course lecture: Computational Modeling of Neuronal Systems
More informationSUPPLEMENTARY INFORMATION
Spatio-temporal correlations and visual signaling in a complete neuronal population Jonathan W. Pillow 1, Jonathon Shlens 2, Liam Paninski 3, Alexander Sher 4, Alan M. Litke 4,E.J.Chichilnisky 2, Eero
More informationSynaptic Input. Linear Model of Synaptic Transmission. Professor David Heeger. September 5, 2000
Synaptic Input Professor David Heeger September 5, 2000 The purpose of this handout is to go a bit beyond the discussion in Ch. 6 of The Book of Genesis on synaptic input, and give some examples of how
More informationT CHANNEL DYNAMICS IN A SILICON LGN. Kai Michael Hynnä A DISSERTATION. Bioengineering
T CHANNEL DYNAMICS IN A SILICON LGN Kai Michael Hynnä A DISSERTATION in Bioengineering Presented to the Faculties of the University of Pennsylvania in Partial Fulfillment of the Requirements for the Degree
More informationA Three-dimensional Physiologically Realistic Model of the Retina
A Three-dimensional Physiologically Realistic Model of the Retina Michael Tadross, Cameron Whitehouse, Melissa Hornstein, Vicky Eng and Evangelia Micheli-Tzanakou Department of Biomedical Engineering 617
More informationA Model for Real-Time Computation in Generic Neural Microcircuits
A Model for Real-Time Computation in Generic Neural Microcircuits Wolfgang Maass, Thomas Natschläger Institute for Theoretical Computer Science Technische Universitaet Graz A-81 Graz, Austria maass, tnatschl
More informationTwo dimensional synaptically generated traveling waves in a theta-neuron neuralnetwork
Neurocomputing 38}40 (2001) 789}795 Two dimensional synaptically generated traveling waves in a theta-neuron neuralnetwork Remus Osan*, Bard Ermentrout Department of Mathematics, University of Pittsburgh,
More informationAnalyzing Neuroscience Signals using Information Theory and Complexity
12th INCF Workshop on Node Communication and Collaborative Neuroinformatics Warsaw, April 16-17, 2015 Co-Authors: Analyzing Neuroscience Signals using Information Theory and Complexity Shannon Communication
More informationProcessing of Time Series by Neural Circuits with Biologically Realistic Synaptic Dynamics
Processing of Time Series by Neural Circuits with iologically Realistic Synaptic Dynamics Thomas Natschläger & Wolfgang Maass Institute for Theoretical Computer Science Technische Universität Graz, ustria
More informationMethods for Estimating the Computational Power and Generalization Capability of Neural Microcircuits
Methods for Estimating the Computational Power and Generalization Capability of Neural Microcircuits Wolfgang Maass, Robert Legenstein, Nils Bertschinger Institute for Theoretical Computer Science Technische
More informationAuto-correlation of retinal ganglion cell mosaics shows hexagonal structure
Supplementary Discussion Auto-correlation of retinal ganglion cell mosaics shows hexagonal structure Wässle and colleagues first observed that the local structure of cell mosaics was approximately hexagonal
More informationSearching for Nested Oscillations in Frequency and Sensor Space. Will Penny. Wellcome Trust Centre for Neuroimaging. University College London.
in Frequency and Sensor Space Oscillation Wellcome Trust Centre for Neuroimaging. University College London. Non- Workshop on Non-Invasive Imaging of Nonlinear Interactions. 20th Annual Computational Neuroscience
More informationLimulus. The Neural Code. Response of Visual Neurons 9/21/2011
Crab cam (Barlow et al., 2001) self inhibition recurrent inhibition lateral inhibition - L16. Neural processing in Linear Systems: Temporal and Spatial Filtering C. D. Hopkins Sept. 21, 2011 The Neural
More informationCHARACTERIZATION OF NONLINEAR NEURON RESPONSES
CHARACTERIZATION OF NONLINEAR NEURON RESPONSES Matt Whiteway whit8022@umd.edu Dr. Daniel A. Butts dab@umd.edu Neuroscience and Cognitive Science (NACS) Applied Mathematics and Scientific Computation (AMSC)
More informationAT2 Neuromodeling: Problem set #3 SPIKE TRAINS
AT2 Neuromodeling: Problem set #3 SPIKE TRAINS Younesse Kaddar PROBLEM 1: Poisson spike trains Link of the ipython notebook for the code Brain neuron emit spikes seemingly randomly: we will aim to model
More informationTemporal whitening by power-law adaptation in neocortical neurons
Temporal whitening by power-law adaptation in neocortical neurons Christian Pozzorini, Richard Naud, Skander Mensi and Wulfram Gerstner School of Computer and Communication Sciences and School of Life
More informationHow Complex Cells Are Made in a Simple Cell Network
How Complex Cells Are Made in a Simple Cell Network Louis Tao Courant Institute, New York University Collaborators: Michael Shelley (Courant, NYU) Robert Shapley (CNS, NYU) David McLaughlin (Courant, NYU)
More information