Information in small neuronal ensemble activity in the hippocampal CA1 during delayed non-matching to sample performance in rats
© Takahashi and Sakurai; licensee BioMed Central Ltd. 2009
Received: 25 November 2008
Accepted: 15 September 2009
Published: 15 September 2009
The matrix-like organization of the hippocampus, with its several inputs and outputs, has given rise to several theories related to hippocampal information processing. Single-cell electrophysiological studies and studies of lesions or genetically altered animals using recognition memory tasks such as delayed non-matching-to-sample (DNMS) tasks support the theories. However, a complete understanding of hippocampal function necessitates knowledge of the encoding of information by multiple neurons in a single trial. The role of neuronal ensembles in the hippocampal CA1 for a DNMS task was assessed quantitatively in this study using multi-neuronal recordings and an artificial neural network classifier as a decoder.
The activity of small neuronal ensembles (6-18 cells) over brief time intervals (2-50 ms) contains accurate information specifically related to the matching/non-matching of continuously presented stimuli (stimulus comparison). The accuracy of the combination of neurons pooled over all the ensembles was markedly lower than those of the ensembles over all examined time intervals.
The results show that the spatiotemporal patterns of spiking activity among cells in the small neuronal ensemble contain much information that is specifically useful for the stimulus comparison. Small neuronal networks in the hippocampal CA1 might therefore act as a comparator during recognition memory tasks.
Hippocampal formation has been identified as an important substrate for declarative memory for a broad range of materials in humans [1, 2]. In contrast, in rodent studies, two views respectively hold that the hippocampus is dedicated to spatial memory processing  and that it associates general memory items . Regarding the spatial view, results of several analyses of the stability of place cells have shown that pattern separation and pattern completion are apparent in neuronal ensembles of the hippocampus [5, 6]. In support of the general view, several experiments directly showed activity related to match/non-match conditions using recognition memory tasks, such as a delayed non-matching-to-sample (DNMS) task [7–10]. Spatial information can be regarded as multiple items that are mutually associated according to temporal relations . Therefore, place cells might code multiple events constructing a place experienced in the past . From this viewpoint, the lines of evidence for the general and spatial views are not contradictory and are consistent with the view of the hippocampus as being capable of auto-associative functions to retrieve entire episodes . On the other hand, the matrix-like organization of the hippocampus with several inputs and outputs has inspired some researchers to propose the hippocampal comparator theory [13–15], which suggests that the hippocampus supports comparison, which might be one element in the match/non-match judgements. In this context, several lines of evidence obtained from analyses of single neurons in spatial and non-spatial behavioral tasks suggest that the hippocampus is critical for episodic-like representations. Nevertheless, our knowledge related to neuronal computations in the working brain of behaving animals is limited; most of it has been inferred exclusively from changes in the firing rates of individual cells accumulated through many trials [16, 17]. Therefore, to understand the hippocampal function completely, the actual encoding of information by the hippocampal neuronal networks of multiple neurons in a single trial during memory tasks must be elucidated. Recently, some researchers have attempted to elucidate the functioning of neuronal networks of the cerebral cortex using an artificial neuronal network classifier as a decoder that enables us to analyze spatiotemporal firing patterns among all observed cells in a single trial [18–20]. Consequently, using multi-neuronal recording and an artificial neural network classifier as a decoder, we analyzed spatiotemporal firing patterns among cells in the hippocampal CA1 of rats. We report its neuronal ensemble code in a single trial of a DNMS task.
Classification performance in neuronal ensembles and combination
Time resolution of code for stimulus comparison
By investigating the degree to which classification performance depended on the bin size of the spike count, we examined the temporal resolution of the ensemble code. Bin sizes of 2-50 ms yielded better performance than larger bin sizes (Figure 3G). A 2-ms bin typically contained zero or one spike. Consequently a few spikes from a small number of neurons are sufficient to encode the stimulus comparison in the hippocampal ensembles.
In this study, the activity of small neuronal ensembles (6-18 cells) in the hippocampal CA1 was used to classify the stimulus comparison with good accuracy over brief time intervals (2-50 ms). They were not useful to classify the stimulus perception, stimulus retention, or motor selection. The accuracy of the ensembles for the stimulus comparison was markedly higher than that of the combination (36 cells).
Small neuronal ensembles in the hippocampal CA1 specifically code stimulus comparison
During test periods in which the rats prepare for go/no-go responses, we examined the stimulus comparison as well as the motor selection. The motor selection of go/no-go responses contains the stimulus comparison of match/non-match tones because, in correct trials, go and no-go responses respectively correspond to non-match and match tones. It is possible that the low accuracy of the ensembles for the motor selection suggest that the activities of the neuronal ensembles are not useful to classify the stimulus comparison. To verify the stimulus comparison without the influence of the motor functions, we examined ensemble activities during match and non-match tones to which identical behavioral responses (go responses) were conducted. Moreover, to exclude the influences of correct and erroneous responses and the particular stimuli being presented, we examined the classification performance of the ensemble activity in relation to the correct-error difference and the particular-stimulus presentation. The results show that, using a linear classifier that can be realized easily in downstream neurons by summating appropriately weighted inputs, we can clearly characterize that at the ensemble level, the most available information in the hippocampal CA1 in each single trial is specifically the stimulus comparison: it is neither the correct-error difference nor the particular-stimulus presentation.
Based on the results of ensemble analyses on a trial-to-trial basis, Deadwyler and colleagues inferred multiple representation and conjunctive encoding in activities of small hippocampal ensembles consisting of 10 neurons . In contrast with the present results, they showed that the stimulus comparison in the hippocampal ensembles was highly correlated with correct and erroneous responses. The reason for the inconsistent results might be differences in the tasks (spatial vs non-spatial versions of DNMS task) and differences in the ensemble analysis tools (canonical discriminant analysis vs. SVM). Furthermore, because we examined neuronal ensemble activities only in the hippocampal CA1, the activity of the hippocampal CA3 analyzed in the Deadwyler's study might explain the inconsistency of the results. The results of the present study illustrate that small neuronal ensembles in the hippocampal CA1 are specifically dedicated to stimulus comparison. The results support the notion that the hippocampus contributes to memory by identifying consistencies across experiences that constitute new associations, as reported from some previous studies [4, 7, 10, 15].
The classification performance for the stimulus comparison of neuronal ensembles is better than that of their combination
Several lines of evidence [21–24] suggested that firing rate modulations of individual neurons and spike timing among neurons contain information in the neuronal ensembles. Given that the small neuronal ensembles in the hippocampal CA1 contains information related only to the firing rate modulations of individual neurons, the classification performance of the combination of all neurons in all the ensembles recorded from different sessions is expected to be the sum of the performance of the ensembles. However, in the present study, almost all classification performances of small neuronal ensembles in the hippocampal CA1 for the stimulus comparison were higher than that of the combination. Moreover, the ensembles yielded better performance in a 2-50 ms window. From this perspective, we presume that the timing of spikes of 2-50 ms resolution in the ensembles recorded in the identical session, such as second or higher order synchronies, contain much information for the stimulus comparison in the hippocampal CA1 at the ensemble level.
Information of the stimulus comparison may be distributed throughout the hippocampal CA1
To elucidate the decoding ability of small neuronal ensembles in the hippocampal CA1, we examined the spatiotemporal firing patterns of small neuronal ensembles consisting of a maximum of 18 cells. Results of a previous study suggested that each of the discreetly located clusters of neurons is dedicated to a different aspect of the spatial version of the DNMS task . Therefore, not all represented information in the hippocampus might be induced from activity of the small neuronal ensembles detected in the present study. That viewpoint, together with the results of analyses of firing rates of individual single cells [7–10], suggests that larger neuronal ensembles in the hippocampus are associated with all examined events in addition to the stimulus comparison in the present study. Nevertheless, our results demonstrate that, in spite of using randomly selected neurons, all small neuronal ensembles in all the events we examined specifically classified the stimulus comparison. For that reason, we infer that the information of stimulus comparison is not restricted, but is instead distributed throughout the hippocampal CA1. It can be retrieved robustly by the ensemble activity of any subcombination in it.
On the other hand, in this study, the SVM classifiers behave as the target of the hippocampal CA1, such as the entorhinal cortex and subiculum. Those targets receive spikes not from the entire hippocampal CA1 but from a subset of it . Consequently, the activity of a small neuronal ensemble in the hippocampal CA1 over brief time intervals during a single trial demonstrates the possibility that the targets of hippocampal CA1 work to perform, accurately, tasks that requires a stimulus comparison, as described in a recent study .
The results show that a neuronal ensemble in the hippocampal CA1 acts as a comparator during a recognition memory task. Our approaches and findings revealed the cognitive functioning of the hippocampal CA1 from the neuronal ensemble activity in a single trial, supporting the development of new solutions for reading cognitive functions from the brain such as brain-machine interfaces [18, 28]. Further assessment of the application of neuronal ensemble activity in the cognitive version of brain-machine interfaces can be achieved through experiments incorporating the sensory feedback of brain-controlled actuators in real-time.
Three male Wistar rats were extensively handled; then food was deprived to approximately 80-90% of their ad libitum body weight. Lights were left on in the colony room between 8 A.M. and 9 P.M Experiments were conducted between 9 A.M. and 6 P.M. The rats were trained to perform an auditory-guided continuous DNMS task [8, 9]. A training session consisted of approximately 200 trials performed for approximately 1 h. The criterion for performance was 80% correct trials per session. In each trial (Figure 2A), one of two tones (high tone: 10 kHz, 85 dB SPL; low tone: 2 kHz, 85 dB SPL) was randomly presented for 15 s following a 5 s delay period. One second after the tone onset, a guillotine door was opened for 3 s to show an illuminated response panel immediately behind the door. Trials were continuously performed with intervening delays (Figure 2B). A food reward was delivered immediately after the go response in a non-match trial, in which the presented tone (high/low) differed from that in the preceding trial. During each trial, a rat was required to make a go response in non-match trials and a no-go response in match trials (Figure 2C). During the delay period, the rat had to remember which stimulus (high/low tones) had been presented most recently. This task design enabled us to dissociate hippocampal neuronal activities associated with a) stimulus perception (difference in activity during sample periods (Figure 2A) between high and low tones) (Figure 2D), b) stimulus retention (difference in activity during delay periods (Figure 2A) between preceding high and low tones) (Figure 2E), c) motor selection (difference in activity during test periods (Figure 2A) between go and no-go responses) (Figure 2F), d) stimulus comparison (difference in activity during test periods between match tone (erroneous go response) and non-match tone (correct go response)) (Figure 2G), e) correct-error comparison (difference in activity during test periods between correct and erroneous responses when match tones was presented) (Figure 2H), and f) particular-stimulus comparison (difference in activity between low test tone periods preceded by high sample tones (low-after-high) and high test tone periods preceded by low sample tones (high-after-low)) (Figure 2I). All behavioral events were controlled using custom-written software running with another software program (Labview; National Instruments Corp., Austin, TX).
Animal preparation and recordings
Neuronal ensemble analysis
We explored codes by counting spikes in successive bins of size w within the interval starting i s s after the test-tone onset and ending i e s after the test-tone onset. We describe the results for different values of these parameters. The default condition was w = 25 ms. Parameter w controls the time resolution of the code: we used w = 2 ms, 4 ms, 8 ms, 10 ms, 25 ms, 50 ms, 100 ms, and 200 ms. The interval parameters were set as i s = -10 and i e = -5 for the stimulus perception because, during this period, a sensory stimulus (high/low tone) is obtained and a rat does not make a go/no-go response. The parameters were i s = -5 and i e = 0 for the stimulus retention because, during this period, a sensory stimulus (high/low tone) is not obtained and must be retained in memory. The parameters were i s = 0 and i e = 1 for the motor selection because, during this period, a rat prepares a motor response (go/no-go response). The parameters were i s = 0 and i e = 1 for the stimulus comparison, correct-error comparison, and particular-stimulus comparison because, during this period, a rat compares a sample tone with a test tone before it makes a go/no-go response.
This vector was used as the input to the decoding classifier. When considering the activity of multiple neurons, we concatenated the corresponding activity vectors and used the concatenated vector as the input to the classifier. The dimensionality of the input is therefore (n+1)N, where N represents the number of neurons in an ensemble or a combination. For a neuronal combination (neurons pooled over all three ensembles), this concatenation step assumes independence among different neurons. For a neuronal ensemble, correlations between simultaneously monitored neurons in an ensemble must contain additional information and must reveal additional aspects of the neural codes.
We used a leave-one-out cross-validation method for training and testing the data. Data were always divided in all cases into a training set and a test set. The training set comprised T-1 trials of each event, whereas the test set included the remaining one trial; T represents the number of trials in one session. Therefore, we used T pairs of training and test sets.
We specifically investigated six functional roles: stimulus perception, stimulus retention, motor selection, stimulus comparison, correct-error comparison and particular-stimulus comparison. For stimulus perception, the labels represent which tones were presented during the sample periods (high or low tone) (Figure 2D). For stimulus retention, the labels denote which tones were presented during the delay periods (high or low tone) (Figure 2E). In correct trials, the motor selection of go/no-go responses can be examined. For motor selection, the labels denote which behavioral responses a rat made during the test periods in correct trials (go or no-go response) (Figure 2F). The go responses can be divided into correct and erroneous go responses. Comparing activities between correct go responses during non-match tones and erroneous go responses during match tones, the stimulus comparison: match tone vs. non-match tone on the identical behavior (go response) can be examined. For stimulus comparison, the labels denote whether the sample tone and test tone differed in trials in which a rat made go responses during the test periods (match tone (erroneous go response) or non-match tone (correct go response)) (Figure 2G). For correct-error comparison, the labels denote whether a rat made correct or erroneous responses during the test periods (correct no-go response or erroneous go response during match tones) (Figure 2H). For particular-stimulus comparison, the labels show which pair of non-match tones was presented during the test periods (high-after-low or low-after-high tone) (Figure 2I). For all three rats examined, the number of trials in which the rats made erroneous no-go responses during non-match tones (< 3) was insufficient for statistical tests and analyses; we did not examine binary labels that include an erroneous no-go response during non-match tones for the stimulus comparison and correct-error comparison.
We trained one binary classifier for each role. The classification performance, as portrayed in all plots, represents the share of correct decoding for test data (i.e., data not used by the classifier during training).
We compared the performance of different statistical classifiers including Perceptron, Fisher's linear discriminant classifier, and Support Vector Machine (SVM)  with a linear kernel using sample ensemble activity. The SVM classifiers yielded the best performance. In addition, the SVM architecture can be easily realized in the hippocampal as a threshold sum of weighted synaptic inputs. Therefore, we used SVM implemented in the software package 'libsvm'http://www.csie.ntu.edu.tw/~cjlin/libsvm.
The linear classifier is expressed as y = sign [g(t)]. The calculation of g(t) is a linear function of the form of , where t i is single-cell activity, w i is the weight of i- th cell, and b is the bias. In fact, SVM finds the optimal weights and bias. We initially tested the performance of the classifier on a small subsample of the data, exploring a large set of parameters. Then we used the optimized parameters for analysis of the complete dataset. The parameters for libsvm were C = 10, S = 0, and T = 0.
The graphs presented in the text show that the N neurons used as the input to the classifier were all possible combinations of N neurons if the number of all possible combinations was < 1,000 or were 1,000 combinations chosen randomly from among all possible combinations of N neurons. We report the average obtained from all random neuronal subensembles and subcombinations.
Test of significance for decoding performance
Using this binomial distribution, we were able to establish the statistical significance of the observed values for the decoding performance.
All analyses were performed using custom-written software (C++) and another software program (MATLAB; The MathWorks Inc., Natick, MA).
This work was supported by the JST PRESTO Program (S.T.), and Grants-in-Aid for Scientific Research (18047012 and 16330140) and the JST CREST Program (Y.S.).
- Cohen NJ, Squire LR: Preserved learning and retention of pattern-analyzing skill in amnesia: dissociation of knowing how and knowing that. Science. 1980, 210 (4466): 207-210. 10.1126/science.7414331.View ArticlePubMedGoogle Scholar
- Scoville WB, Milner B: Loss of recent memory after bilateral hippocampal lesions. J Neurol Neurosurg Psychiatry. 1957, 20 (1): 11-21. 10.1136/jnnp.20.1.11.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Keefe J, Dostrovsky J: The hippocampus as a spatial map. Preliminary evidence from unit activity in the freely-moving rat. Brain Res. 1971, 34 (1): 171-175. 10.1016/0006-8993(71)90358-1.View ArticlePubMedGoogle Scholar
- Eichenbaum H, Dudchenko P, Wood E, Shapiro M, Tanila H: The hippocampus, memory, and place cells: is it spatial memory or a memory space?. Neuron. 1999, 23 (2): 209-226. 10.1016/S0896-6273(00)80773-4.View ArticlePubMedGoogle Scholar
- Leutgeb S, Leutgeb JK, Treves A, Moser MB, Moser EI: Distinct ensemble codes in hippocampal areas CA3 and CA1. Science. 2004, 305 (5688): 1295-1298. 10.1126/science.1100265.View ArticlePubMedGoogle Scholar
- Nakazawa K, McHugh TJ, Wilson MA, Tonegawa S: NMDA receptors, place cells and hippocampal spatial memory. Nat Rev Neurosci. 2004, 5 (5): 361-372. 10.1038/nrn1385.View ArticlePubMedGoogle Scholar
- Deadwyler SA, Bunn T, Hampson RE: Hippocampal ensemble activity during spatial delayed-nonmatch-to-sample performance in rats. J Neurosci. 1996, 16 (1): 354-372.PubMedGoogle Scholar
- Sakurai Y: Hippocampal cells have behavioral correlates during the performance of an auditory working memory task in the rat. Behav Neurosci. 1990, 104 (2): 253-263. 10.1037/0735-7044.104.2.253.View ArticlePubMedGoogle Scholar
- Sakurai Y: Involvement of auditory cortical and hippocampal-neurons in auditory Working-memory and reference memory in the rat. J Neurosci. 1994, 14 (5): 2606-2623.PubMedGoogle Scholar
- Wood ER, Dudchenko PA, Eichenbaum H: The global record of memory in hippocampal neuronal activity. Nature. 1999, 397 (6720): 613-616. 10.1038/17605.View ArticlePubMedGoogle Scholar
- Muller RU, Stead M, Pach J: The hippocampus as a cognitive graph. J Gen Physiol. 1996, 107 (6): 663-694. 10.1085/jgp.107.6.663.PubMed CentralView ArticlePubMedGoogle Scholar
- Lisman JE: Relating hippocampal circuitry to function: recall of memory sequences by reciprocal dentate-CA3 interactions. Neuron. 1999, 22 (2): 233-242. 10.1016/S0896-6273(00)81085-5.View ArticlePubMedGoogle Scholar
- Gray JA: The neuropsychology of anxiety. 1982, Great Britain: Oxford University PressGoogle Scholar
- Olton DS: Characteristics of spatial memory. 1978, Hillsdale, NJ: ErlbaumGoogle Scholar
- Vinogradova OS: Hippocampus as comparator: role of the two input and two output systems of the hippocampus in selection and registration of information. Hippocampus. 2001, 11 (5): 578-598. 10.1002/hipo.1073.View ArticlePubMedGoogle Scholar
- Georgopoulos A, Schwartz A, Kettner R: Neuronal population coding of movement direction. Science. 1986, 233 (4771): 1416-1419. 10.1126/science.3749885.View ArticlePubMedGoogle Scholar
- Hubel DH, Wiesel TN: Ferrier lecture. Functional architecture of macaque monkey visual cortex. Proc R Soc Lond B Biol Sci. 1977, 198 (1130): 1-59. 10.1098/rspb.1977.0085.View ArticlePubMedGoogle Scholar
- Chapin J, Moxon K, Markowitz R, Nicolelis M: Real-time control of a robot arm using simultaneously recorded neurons in the motor cortex. Nature Neurosci. 1999, 2 (7): 664-670. 10.1038/10223.View ArticlePubMedGoogle Scholar
- Hung CP, Kreiman G, Poggio T, DiCarlo JJ: Fast readout of object identity from macaque inferior temporal cortex. Science. 2005, 310 (5749): 863-866. 10.1126/science.1117593.View ArticlePubMedGoogle Scholar
- Mehring C, Rickert J, Vaadia E, Cardosa de Oliveira S, Aertsen A, Rotter S: Inference of hand movements from local field potentials in monkey motor cortex. Nat Neurosci. 2003, 6 (12): 1253-1254. 10.1038/nn1158.View ArticlePubMedGoogle Scholar
- Abeles M: Corticonics. 1991, Cambridge: Cambridge Univ. PressView ArticleGoogle Scholar
- Abeles M, Bergman H, Margalit E, Vaadia E: Spatiotemporal firing patterns in the frontal cortex of behaving monkeys. J Neurophysiol. 1993, 70 (4): 1629-1638.PubMedGoogle Scholar
- Harris K, Csicsvari J, Hirase H, Dragoi G, Buzsaki G: Organization of cell assemblies in the hippocampus. Nature. 2003, 424 (6948): 552-556. 10.1038/nature01834.View ArticlePubMedGoogle Scholar
- Harris KD, Hirase H, Leinekugel X, Henze DA, Buzsaki G: Temporal interaction between single spikes and complex spike bursts in hippocampal pyramidal cells. Neuron. 2001, 32 (1): 141-149. 10.1016/S0896-6273(01)00447-0.View ArticlePubMedGoogle Scholar
- Hampson RE, Simeral JD, Deadwyler SA: Distribution of spatial and nonspatial information in dorsal hippocampus. Nature. 1999, 402 (6762): 610-614. 10.1038/45154.View ArticlePubMedGoogle Scholar
- Freund T, Buzsaki G: Interneurons of the hippocampus. Hippocampus. 1996, 6 (4): 347-470. 10.1002/(SICI)1098-1063(1996)6:4<347::AID-HIPO1>3.0.CO;2-I.View ArticlePubMedGoogle Scholar
- O'Maria SM, Sanchez-Vives V, Brotons-Mas JR, O'Hare E: Roles for the subiculum in spatial information processing, memory, motivation and the temporal control of behaviour. Progress in Neuro-Psychopharmacology and Biological Psychiatry. 2009, 33 (5): 782-790. 10.1016/j.pnpbp.2009.03.040.View ArticleGoogle Scholar
- Nicolelis M: Actions from thoughts. Nature. 2001, 409 (6818): 403-407. 10.1038/35053191.View ArticlePubMedGoogle Scholar
- Sakurai Y, Takahashi S: Dynamic synchrony of firing in the monkey prefrontal cortex during working-memory tasks. J Neurosci. 2006, 26 (40): 10141-10153. 10.1523/JNEUROSCI.2423-06.2006.View ArticlePubMedGoogle Scholar
- Takahashi S, Anzai Y, Sakurai Y: A new approach to spike sorting for multi-neuronal activities recorded with a tetrode - how ICA can be practical. Neurosci Research. 2003, 46 (3): 265-272. 10.1016/S0168-0102(03)00103-2.View ArticleGoogle Scholar
- Takahashi S, Anzai Y, Sakurai Y: Automatic sorting for multi-neuronal activity recorded with tetrodes in the presence of overlapping spikes. J Neurophysiol. 2003, 89 (4): 2245-2258. 10.1152/jn.00827.2002.View ArticlePubMedGoogle Scholar
- Takahashi S, Sakurai Y: Real-time and automatic sorting of multi-neuronal activity for sub-millisecond interactions in vivo. Neuroscience. 2005, 134 (1): 301-315. 10.1016/j.neuroscience.2005.03.031.View ArticlePubMedGoogle Scholar
- Takahashi S, Sakurai Y: Coding of spatial information by soma and dendrite of pyramidal cells in the hippocampal CA1 of behaving rats. Eur J Neurosci. 2007, 26 (7): 2033-2045. 10.1111/j.1460-9568.2007.05827.x.View ArticlePubMedGoogle Scholar
- Csicsvari J, Hirase H, Czurko A, Buzsaki G: Reliability and state dependence of pyramidal cell-interneuron synapses in the hippocampus: an ensemble approach in the behaving rat. Neuron. 1998, 21 (1): 179-189. 10.1016/S0896-6273(00)80525-5.View ArticlePubMedGoogle Scholar
- Ranck JB: Studies on Single Neurons in Dorsal Hippocampal Formation and Septum in Unrestrained Rats.1. Behavioral Correlates and Firing Repertoires. Exp Neurol. 1973, 41 (2): 462-531. 10.1016/0014-4886(73)90290-2.View ArticleGoogle Scholar
- Schmitzer-Torbert N, Jackson J, Henze D, Harris K, Redish AD: Quantitative measures of cluster quality for use in extracellular recordings. Neuroscience. 2005, 131 (1): 1-11. 10.1016/j.neuroscience.2004.09.066.View ArticlePubMedGoogle Scholar
- Maurer AP, Cowen SL, Burke SN, Barnes CA, McNaughton BL: Organization of hippocampal cell assemblies based on theta phase precession. Hippocampus. 2006, 16 (9): 785-794. 10.1002/hipo.20202.View ArticlePubMedGoogle Scholar
- Vapnik VN: Statistical learning theory. 1998, New York: WileyGoogle Scholar