|
|
||||||||
1Vanderbilt Brain Institute, 2Department of Psychology, 3Center for Integrative and Cognitive Neuroscience, Vanderbilt University, Nashville, Tennessee, and 4Vanderbilt Vision Research Center, Nashville, Tennessee
Submitted 9 May 2007; accepted in final form 6 September 2007
| ABSTRACT |
|---|
|
|
|---|
100 ms before saccade onset. The minimum in response variability occurred
20 ms earlier for the easier search task compared with the harder one. This difference is positively correlated with the difference in saccade reaction times for the two tasks. These findings show that a multivariate analysis can provide a measure of neuronal interactions and characterize the spiking activity of FEF neurons in the context of a population of neurons. | INTRODUCTION |
|---|
|
|
|---|
Recent developments in multivariate point process modeling of neural responses have provided an analytical framework to characterize neural activity in the context of interactions between simultaneously recorded neurons (Brown et al. 2002
; Okatan et al. 2005
; Truccolo et al. 2005
). This analysis models the instantaneous firing rate of a neuron using its own spiking history and that of other simultaneously recorded neurons as covariates in the model. The significance of each covariate provides an estimate of its contribution to the response of the modeled neuron. If the multivariate model's estimate of instantaneous firing rate is significantly improved by including covariates representing the activity of other neurons, it is evidence of interactions among neurons in the recorded ensemble. Such interactions may be direct synaptic connections between neurons or may be mediated polysynaptically or by shared input. Compared with conventional univariate estimates of neuronal activity, such as the peristimulus time histogram (PSTH), the multivariate approach can distinguish between a neuron's response and its response in the context of interactions in a population of neurons. Furthermore, point process modeling of neural activity preserves spike timing information which is distorted by measures that average over time-interval windows (e.g., PSTH). Traditional approaches for analysis of interactions between neurons, such as the covariogram and joint PSTH (Aertsen et al. 1989
; Brody 1999a
,b
; Constantinidis et al. 2001
), are limited to pairwise comparisons and do not provide adequate measures of ensemble interactions. The multivariate point process model estimates the instantaneous firing rate of a neuron in real-time without limits on ensemble size. Nonetheless, it can extract the same interactions as pairwise measures.
We show that the multivariate approach can accurately model spiking activity in the FEF and characterize interactions among simultaneously recorded neurons while monkeys perform a visual search task. To better understand how neuronal interactions in the FEF contribute to target selection, we addressed the following questions. 1) Do interactions among FEF neurons differ for hard and easy search? 2) Are there differences in interactions between the three classes of FEF neurons? 3) Are interactions of movement-related neurons associated with saccade onset times?
| METHODS |
|---|
|
|
|---|
Activity of FEF neurons was recorded in macaques performing a visual search task in which they were required to saccade to a singleton target defined by color (Sato et al. 2001
). Each trial began with the monkey fixating a central spot for
600 ms. A target was presented at one of eight isoeccentric locations equally spaced around the fixation spot (Fig. 1A, inset). The other seven locations contained distractor stimuli. Monkeys were given a juice reward for making a saccade to the target location and holding their gaze on the target for
400 ms. There were two levels of task difficulty, hard and easy, determined by the degree of target-distractor similarity. The hard task contained a green target among yellow-green distractors. The easy task contained a green target among red distractors. Recordings were made simultaneously from two to four tungsten electrodes placed in the rostral bank of the arcuate sulcus. A neural ensemble was defined as a set of simultaneously recorded neurons with overlapping receptive fields. Our data set consisted of 91 neurons in 29 ensembles from one monkey (Macaca mulatta) and 21 neurons in 7 ensembles from a second monkey (Macaca radiata). Of the 36 ensembles, 15 contained two neurons, 10 contained three neurons, 6 contained four neurons, 4 contained five neurons, and 1 contained eight neurons. Spikes were sorted off-line using principal components analysis (Plexon). Because data from each monkey were collected during the same behavioral tasks and were similar, we pooled data from both monkeys.
|
Each monkey was surgically implanted with a head post, a scleral eye coil, and a recording chamber. Surgery was conducted under aseptic conditions with isofluorane anesthesia (Schall et al. 1995a
). Antibiotics and analgesics were administered postoperatively. All experimental procedures were performed in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and approved by the Vanderbilt Institutional Animal Care and Use Committee.
Data analysis
To assess simultaneously the interactions of several neurons, we used a point process multivariate analysis (Okatan et al. 2005
; Truccolo et al. 2005
). The point process framework can provide for comparisons of arbitrarily large ensembles of simultaneously recorded neurons that pairwise measures cannot. We constructed a statistical model of the firing rate of a neuron by incorporating its firing history and the firing history of other neurons in its ensemble from stimulus onset to saccade onset (5,315 correct hard search trials, 7,414 correct easy search trials). We used a modified version of the generalized linear model (GLM) approach recently applied by Truccolo et al. (2005)
. The modification was necessary because the fitted GLMs for the hard and easy tasks may be different. Therefore unless they are nested, comparing such models statistically is intractable. We modified the GLM in the following way so that the data for both the hard and the easy tasks (for a particular neuron) were combined in a single GLM.
Using the theory of point processes, we represented recorded spike trains as sets of discrete event times. We modeled the instantaneous firing rate (conditional intensity function) of a neuron as a combination of terms of covariates (Truccolo et al. 2005
). The conditional intensity function (
t) is more informative of the instantaneous firing rate than univariate measures (e.g., PSTH) because its estimate is derived in the context of interacting neurons. We predicted the firing of a neuron using its firing history (autoregressive process) and the firing history of other neurons recorded simultaneously (crossregressive process). We concatenated data from all correct trials for each neuron and used a GLM to predict the firing rate of a neuron as
![]() |
t is the firing rate at time t, µ is a baseline term associated with the hard search condition,
is a baseline term associated with the easy search condition, {
i} is the set of hard search autoregressive (AR) parameters, {
i} is the set of easy search AR parameters, {
j} is the set of hard search crossregressive (CR) parameters (1 for each neuron in the ensemble at each lag), {
j} is the set of easy search CR parameters, and
Nt-k(c) is the spike count in the k th ms before the current time t, for neuron c in the ensemble. Q and R are the AR and CR lags of the model, respectively. The indicator variables x1 and x2 combine the parameters associated with each task into a single model. x1 is 0 for easy search trials and 1 for hard search trials. x2 is 1 for easy search trials and 0 for hard search trials. Because of constraints on the length of the recordings (relative to the firing rate) we set Q = R = 30. This constraint was not of consequence to our analysis because the parameter fits stabilized well before lags of 30 ms. In the GLM above, the AR parameters describe the timing of the modeled neuron's dependence on its firing history and the CR ones describe the timing of interactions between neurons. To compare models of hard and easy visual search trials, x1 and x2 terms were merged to create a separate GLM for each task.
We fit the GLM using an iteratively reweighed least squares algorithm (McCullagh and Nelder 1989
). This algorithm provides a robust maximum likelihood estimate of model parameters. If the assumptions of the GLM are met, the fitted model's residuals should have a normal distribution around mean 0 and constant variance with no autocorrelations (McCullagh and Nelder 1989
; Truccolo et al. 2005
). Thus we examined the residuals of each of our model fits.
To compare the fit of nested models, we used likelihood ratio tests. For each neuron, we compared the model deviance (D = –2logL) for AR-only models and for AR-CR models for hard and easy search trials separately. This deviance comes from a
2 distribution. Thus we can test the hypothesis that adding CR terms to an AR model does not improve the GLM fits (McCullagh and Nelder 1989
). If the likelihood ratio is large, the modeled neuron's response depends heavily on the ensemble neuron's responses. We measured variability in
t across neurons using a standardized measure of variability, the coefficient of variation (CV) (SD/mean). All analyses were performed in R (R Development Core Team 2006
).
Model fits
To show that the GLM accurately accounted for the firing of FEF neurons, we compared the conditional intensity functions (
t, instantaneous firing rates) against the observed spike trains. Figure 2A shows the modeled intensity (black curve) and the observed spike train (gray vertical lines) for one example neuron over the course of 300 ms. The model covariates are the neuron's spiking history and the history of a second neuron in the ensemble. The modeled rate closely follows the observed spike times. Note that the magnitude of the intensity increases with the frequency of spikes. The brief (1 ms) decrease in the intensity after a spike likely corresponds to the neuron's absolute refractory period (Truccolo et al. 2005
). This decrease is evident in Fig. 2B, which depicts the spike-triggered average intensity for the neuron in Fig. 2A. We also examined the Pearson residuals for each model. If the Pearson residuals were distributed normally with mean 0 and variance 1, the firing left unexplained by the model was insignificant (McCullagh and Nelder 1989
). For the model in Fig. 2, which was representative of the model fits of FEF neurons in our sample, the mean of the Pearson residuals was 6.08 x 10–4, the variance was 0.892, and they were distributed normally. Models for all neurons had normally distributed Pearson residuals (Shapiro-Wilk test, P < 0.05).
|
A simulation was performed to test whether the GLM analysis extracted the same interactions that standard measures of pairwise interactions do. Figure 3, A and B, shows a comparison of the GLM fitting results for a simulated pair of spike trains with a standard measure of interactions between pairs of neurons, the covariogram (Brody 1999a
). The CR parameter values in Fig. 3B match the lags of high correlation between the pair of simulated spike trains in Fig. 3A. Significant parameters (by Wald Z test, P = 0.05) corresponded to interactions between the covariate and the target neuron at that lag, although not necessarily "monosynaptically." Positive parameter values corresponded to high probability that the covariate neuron excited the modeled neuron. Negative values corresponded to high probability of inhibition. Nonsignificant parameters (gray points) corresponded to low probability that the covariate interacted with the modeled neuron at those lags.
|
10 ms and was inhibitory (negative parameter values). This likely corresponded to effects of absolute and relative refractory periods (Truccolo et al. 2005There were no significant differences between AR or CR parameters between hard and easy tasks. Figure 3, C and D, displays the overall pattern of history dependence on a neuron's firing.
| RESULTS |
|---|
|
|
|---|
Because both monkeys showed similar differences in performance for the easy and hard search tasks, these data were pooled. The easy and hard search tasks resulted in significant differences in percent correct (hard task = 71.6%; easy task = 96.6%; Wilcoxon rank sum test, P < 10–15) and mean saccade reaction time (hard task = 237 ms; easy task = 196 ms; P < 10–15). Figure 1A shows the densities of saccade reaction times (RTs) for the hard and easy search tasks. In addition to the significant difference in the mean RT (44.3 ms) for the two tasks, the variability of saccade RT was greater for the hard search task (SD = 65.3 ms) than the easy search task (SD = 41.2 ms). The mean difference in reaction time (hard minus easy) within each session was 44.4 ± 3.6 (SE) ms, with a minimum of 15.4 ms (Fig. 1B).
Do interactions among FEF neurons differ for hard and easy search?
We asked whether the amount of interactions among neurons in the FEF was affected by the difficulty of the visual search task. To measure interactions between neurons, we used a point process multivariate analysis to model the instantaneous firing rate of each neuron for each task (hard and easy) taking into account the firing history of all neurons in the recorded ensemble (see METHODS). The model estimated the conditional intensity function (instantaneous firing rate) of the neuron with respect to the covariates (AR and CR) when a saccade was made to a target located within the receptive field. We computed the likelihood ratios of models for each neuron for each task by subtracting the deviance of the AR-CR model from the deviance of the AR model. Large likelihood ratios indicated that including ensemble neurons in the model of a neuron greatly improved the prediction of the firing of that neuron. Because the likelihood ratio is distributed as
2, it is an ideal measure of the degree to which addition of CR covariates improved the estimate of the firing rate (McCullagh and Nelder 1989
).
We found that for 63.7% of easy task models and 51.6% of hard task models,
2 P values were <0.05, indicating that the inclusion of CR covariates improved the prediction of firing rate in the majority of models. We compared the likelihood ratios of the two pairs of models (hard vs. easy) to determine under which task the addition of ensemble responses (CR covariates) improved the model the most. The model that benefited the most from addition of CR covariates was judged to convey more interactions. Figure 4A is a histogram of the difference between likelihood ratios for each neuron (easy minus hard). The histogram is shifted significantly to the right of zero (paired Wilcoxon rank sum test, P = 3.05 x 10–3). Thus we could better predict the firing of neurons in the easy task than the hard task when including the firing history of other neurons recorded simultaneously, indicating that neurons interacted more with each other during easy visual search. Simulations revealed that, regardless of visual search task, adding randomly firing simulated neurons to models of FEF neurons did not improve those model fits. Therefore improvements to models of FEF neurons by adding CR covariates were caused by neuronal interactions, whether monosynaptic, polysynaptic, or through shared input. There was no significant correlation (P = 0.490) between RT difference (hard minus easy) and likelihood ratio difference (easy minus hard), possibly because of low variability in RTs.
|
Are there differences in interactions between the three classes of FEF neurons?
We classified each neuron as visual-related, movement-related, or visual-and-movement–related. There were 46 visual-related neurons, 16 movement-related neurons, 49 visual-and-movement–related neurons, and 1 that was unclassified. Figure 5 shows the responses of representative neurons of each class from our sample of recordings. The model for each showed improved fit on addition of CR covariates. Each PSTH in this figure shows the average firing rate of the neuron when the target was presented in the neuron's receptive field. Visual-related neurons had clear responses
50 ms after target onset (Fig. 5A), visual-and-movement–related neurons had similar visual latencies and increased firing leading up to saccades (Fig. 5B), and movement-related neurons fired at baseline until
60 ms before saccades (Fig. 5C). The difference in visual-and-movement–and movement-related responses during the hard and easy tasks correlated with the difference in mean RT between the hard (gray arrowhead) and easy (black arrowhead) tasks.
|
Are interactions among FEF neurons time-locked to saccade onset?
Because movement-related neurons in the FEF characteristically increase their firing leading up to saccades (Hanes and Schall 1996
) (see Fig. 5C), we analyzed the conditional intensity functions (
t, instantaneous firing rates) for movement-related neurons to determine whether the timing of interactions reflected the decision of the monkeys to move their eyes. To observe the neuronal responses around saccades, we modeled the intensity of these neurons from target onset to 50 ms after the saccade. Figure 6 shows the mean intensity (
t) for the 16 movement-related neurons during the hard task (Fig. 6A) and during the easy task (Fig. 6B). The gray curves show the mean intensities of AR models during the 100 ms leading up to saccades for the hard task; the black curves are the mean intensities of AR-CR models.
|
t leading up to saccades, relative to
t for AR models. Thus interactions between neurons statistically accounted for a significant portion of the presaccadic activity in movement-related neurons. Figure 6, E and F, shows the normalized mean PSTH (gray curve) and intensity function (black curve) for hard trials (E) and easy trials (F). A comparison of the intensity functions and PSTHs shows that the AR-CR models accurately describe both the magnitude and the dynamics of the PSTHs for the two search tasks. The intensity functions are shifted
30 ms to the right of the PSTHs, reflecting the integration of 30 ms of firing history into the estimate of the intensities.
To further explore the effect of ensemble interactions on the activity of movement-related neurons, we compared the variability of the conditional intensity functions for each search task. The CV (SD/mean) is a measure of standardized variability frequently used in neuroscience and may be interpreted as a "noise-to-signal" ratio (de Ruyter van Steveninck et al. 1997
; Feng and Brown 1999
; Stein and Matthews 1965
; Stevens and Zador 1998
). We used the CV to measure changes in instantaneous firing rate of movement-related neurons that occurred just before saccades. Decreases in the CV over time reflected less variability (or noise) in the system.
Figure 7 shows the CV of the mean intensity function for the movement-related neurons for the hard (gray curves) and easy (black curves) tasks for AR-CR models (Fig. 7A) and AR models (Fig. 7B). For both hard and easy tasks, the CV decreased leading up to saccades in the AR-CR models until just before saccades, at which point the CV increased. Thus addition of CR covariates decreased standardized variability (CV) leading up to saccades.
|
20 ms earlier in the easy task, despite the fact that we aligned intensities on saccade times. The difference in minimum CV times between easy and hard search was closer for AR-CR models (Fig. 7A) than for AR models (Fig. 7B). Presumably, this was because of poorer fits of AR models. It is important to note that, although intensities were shifted to the right of PSTHs, the timing results in Fig. 7 compare identical measures to each other (i.e., intensity to intensity and PSTH to PSTH). Thus the differences in time course of the CV between the hard and easy tasks reflect the dynamics of the conditional intensity functions. We compared the results using the intensity function to a standard model of the firing rate, the PSTH. Figure 7C shows the CV of the PSTHs of movement-related neurons. Similar to the results obtained from the intensity functions (Fig. 7, A and B), the PSTH CV decreased leading up to saccades and increased 5–10 ms before saccades. However, there was no difference between time of increase for hard and easy tasks. It is important to keep in mind that comparing changes in the PSTH between hard and easy tasks cannot reveal the same information as the intensity can, because the latter measures changes in the system in the context of interacting neurons. Thus the intensity function is a better measure of instantaneous firing rate than the PSTH to measure time of response variability in the FEF.
| DISCUSSION |
|---|
|
|
|---|
These results indicate that easier visual search tasks are associated with greater interactions among populations of FEF neurons. This is seemingly counterintuitive. After all, why should neurons interact more during a task that seemingly requires fewer resources to solve? The average firing rates are not significantly different between the two conditions, so there must be a difference in the structure of the spike trains. There are at least two potential explanations. 1) During hard tasks, FEF neurons fire less randomly, requiring fewer interactions to accomplish the same goals (target selection and saccade preparation). 2) FEF neurons fire more randomly during hard tasks, requiring more information from other neurons. We measured the CV of the interspike intervals during hard and easy tasks. We found no significant difference between the two (P = 0.255), although there was a trend toward larger CV during the easy task. Further studies are required to determine whether FEF neurons' firing patterns differ between tasks.
Hanes and Schall (1996)
determined that single movement-related neurons in the FEF could reliably predict saccade reaction time. This finding leads to the question of how movement-related neurons integrate visual information. It is possible that interactions between the three neuron classes in the FEF fully predict the timing of saccades. Simultaneous recordings from the FEF and other prefrontal areas (e.g., area 46) or parietal areas (e.g., lateral intraparietal cortex) would address this issue. In concert with a point process analysis, such an experiment would describe the timing of interactions between areas and determine how neurons in other cortical areas interact with FEF neurons to decide to move the eyes.
We propose that the timing of changes in the intensity function of movement-related neurons reveals the time at which networks of neurons decide to initiate a saccade. For models that include ensemble activity, the CV decreases until just before saccades, which coincides with the period of time when information about the decision to saccade accumulates. The CV is a second-order measure. As such, it reflects changes in noise in the system. Thus it provides an estimate of the time between the saccade decision time of FEF movement-related neurons and execution of the saccade. Because the CV reflects noise in the system, remarkably, this also suggests that the noise in the system actually decreases as the firing rate of these neurons increases leading up to saccades (until the increase in CV just before saccades). Therefore we show not only that movement-related neurons have activity sufficient to trigger a saccade (Hanes and Schall 1996
), but that their presaccadic activity reflects a decrease in firing variability.
The difference in time of increase of CV between hard and easy tasks has important implications for when the decision to make a saccade occurs. In our experiments, the neurons spent 20 ms longer reaching their minimum in response variability in the hard task than in the easy one. This may correspond to earlier decision times in the easy task than in the hard one. Therefore it seems that the decision to move the eyes is mediated by the movement-related neurons. The larger likelihood ratios for movement-related neurons versus visual-related neurons suggest that the movement-related neurons receive greater interactions from other FEF neurons. This agrees with models of visual search in the FEF, which assert that projections from visual-related neurons to movement-related neurons transform information about visual stimuli into a saccade execution plan (Thompson et al. 1996
).
How does the estimate of 20-ms difference in time of minimum response variability correspond to differences in saccade RT? The mean RT difference between the hard and easy tasks was 44.3 ms. The CV during the 100 ms after target onset did not decrease for visual-related neurons (data not shown). This may be because of the location of presynaptic neurons connecting with these visual-related neurons. If, as the anatomy suggests (Schall et al. 1995b
), visual-related neurons receive synaptic connections from visually responsive neurons in the parietal, temporal, and occipital cortices, we would not have observed visual-related FEF neurons in the context of those interactions. Simultaneous recording from the FEF and cortical areas that project to the FEF would address this issue. Such an experiment would address whether sensorimotor integration occurs in a single bottom-up volley or is the result of continued flow of information between neurons with sensory responses and those with motor responses (Riehle et al. 1997
; Woodman et al. 2007
).
We showed that accounting for ensemble activity is a powerful method of modeling the firing rate of a neuron. There are several ways to measure activity of single neurons, however. A continuous function best represents the firing of a neuron if it accounts for the synaptic input to the neuron (a network property) and approximates the instantaneous firing rate of the neuron. The intensity function has two advantages over other measures of neuronal firing, such as the PSTH. First, it accounts for network activity (the CR covariates) and second, it approximates instantaneous firing rate better than the PSTH (Truccolo et al. 2005
). We have shown that using the intensity function to describe the firing of FEF neurons reveals the time-course of activity leading to the decision to saccade more sensitively than the PSTH.
Although the point process model can account for a large number of influences on a neuron's firing, it is not a mechanistic model. It cannot, therefore, distinguish between neurons that are synaptically connected and neurons that share common input. This limitation is shared by other common techniques such as the covariogram and joint PSTH. A challenge for future studies of neuronal interactions is to include knowledge about anatomical connections and biophysical properties of neurons in mechanistic models of networks of neurons.
We postulate that the distinction between easy and hard tasks, as detected by our results (e.g., difference in lag between minimum CV and saccade) is continuous, rather than dichotomous. Theoretically, there must exist a minimal set of neurons that are required to complete a given visual search task and a latest time before saccade that the decision is made. Our results show that easier visual search tasks are associated with greater interactions among populations of FEF neurons and may result in earlier saccade decision times.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: A. F. Rossi, 301 Wilson Hall, 111 21st Ave. S., Vanderbilt Univ., Nashville, TN 37240 (E-mail: andrew.rossi{at}vanderbilt.edu)
| REFERENCES |
|---|
|
|
|---|
Bichot NP, Schall JD. Effects of similarity and history on neural mechanisms of visual selection. Nat Neurosci 2: 549–554, 1999.[CrossRef][ISI][Medline]
Brody CD. Correlations without synchrony. Neural Comput 11: 1537–1551, 1999a.
Brody CD. Disambiguating different covariation types. Neural Comput 11: 1527–1535, 1999b.
Brown EN, Barbieri R, Ventura V, Kass RE, Frank LM. The time-rescaling theorem and its application to neural spike train data analysis. Neural Comput 14: 325–346, 2002.
Bruce CJ, Goldberg ME. Primate frontal eye fields. I. Single neurons discharging before saccades. J Neurophysiol 53: 603–635, 1985.
Constantinidis C, Franowicz MN, Goldman-Rakic PS. Coding specificity in cortical microcircuits: a multiple-electrode analysis of primate prefrontal cortex. J Neurosci 21: 3646–3655, 2001.
de Ruyter van Steveninck RR, Lewen GD, Strong SP, Koberle R, Bialek W. Reproducibility and variability in neural spike trains. Science 275: 1805–1808, 1997.
Duncan J, Humphreys GW. Visual search and stimulus similarity. Psychol Rev 96: 433–458, 1989.[CrossRef][ISI][Medline]
Feng JF, Brown D. Coefficient of variation of interspike intervals greater than 0.5. How and when? Biol Cybern 80: 291–297, 1999.[CrossRef][ISI][Medline]
Hanes DP, Schall JD. Neural control of voluntary movement initiation. Science 274: 427–430, 1996.
Helminski JO, Segraves MA. Macaque frontal eye field input to saccade-related neurons in the superior colliculus. J Neurophysiol 90: 1046–1062, 2003.
Hikosaka O, Wurtz RH. Visual and oculomotor functions of monkey substantia nigra pars reticulata. III. Memory-contingent visual and saccade responses. J Neurophysiol 49: 1268–1284, 1983.
McCullagh P, Nelder JA. Generalized Linear Models (2nd ed.). Boca Raton, FL: Chapman and Hall/CRC, 1989.
Murthy A, Thompson KG, Schall JD. Dynamic dissociation of visual selection from saccade programming in frontal eye field. J Neurophysiol 86: 2634–2637, 2001.
Okatan M, Wilson MA, Brown EN. Analyzing functional connectivity using a network likelihood model of ensemble neural spiking activity. Neural Comput 17: 1927–1961, 2005.
R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria, 2006.
Riehle A, Kornblum S, Requin J. Neuronal correlates of sensorimotor association in stimulus-response compatibility. J Exp Psychol Hum Percept Perform 23: 1708–1726, 1997.[CrossRef][ISI][Medline]
Sato T, Murthy A, Thompson KG, Schall JD. Search efficiency but not response interference affects visual selection in frontal eye field. Neuron 30: 583–591, 2001.[CrossRef][ISI][Medline]
Schall JD. Neuronal activity related to visually guided saccades in the frontal eye fields of rhesus monkeys: comparison with supplementary eye fields. J Neurophysiol 66: 559–579, 1991.
Schall JD, Hanes DP. Neural basis of saccade target selection in frontal eye field during visual search. Nature 366: 467–469, 1993.[CrossRef][Medline]
Schall JD, Hanes DP, Thompson KG, King DJ. Saccade target selection in frontal eye field of macaque. I. Visual and premovement activation. J Neurosci 15: 6905–6918, 1995a.
Schall JD, Morel A, King DJ, Bullier J. Topography of visual cortex connections with frontal eye field in macaque: convergence and segregation of processing streams. J Neurosci 15: 4464–4487, 1995b.[Abstract]
Schall JD, Thompson KG. Neural selection and control of visually guided eye movements. Annu Rev Neurosci 22: 241–259, 1999.[CrossRef][ISI][Medline]
Segraves MA, Goldberg ME. Functional properties of corticotectal neurons in the monkey's frontal eye field. J Neurophysiol 58: 1387–1419, 1987.
Sommer MA, Wurtz RH. Frontal eye field neurons orthodromically activated from the superior colliculus. J Neurophysiol 80: 3331–3335, 1998.
Sommer MA, Wurtz, RH. Frontal eye field sends delay activity related to movement, memory, and vision to the superior colliculus. J Neurophysiol 85: 1673–1685, 2001.
Stein RB, Matthews PBC. Differences in variability of discharge frequency between primary and secondary muscle spindle afferent endings of the cat. Nature 208: 1217–1218, 1965.[CrossRef][Medline]
Stevens CF, Zador AM. Input synchrony and the irregular firing of cortical neurons. Nat Neurosci 1: 210–217, 1998.[CrossRef][ISI][Medline]
Thompson KG, Bichot NP, Schall JD. Dissociation of visual discrimination from saccade programming in macaque frontal eye field. J Neurophysiol 77: 1046–1050, 1997.
Thompson KG, Hanes DP, Bichot NP, Schall JD. Perceptual and motor processing stages identified in the activity of macaque frontal eye field neurons during visual search. J Neurophysiol 76: 4040–4055, 1996.
Truccolo W, Eden UT, Fellows MR, Donoghue JP, Brown EN. A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. J Neurophysiol 93: 1074–1089, 2005.
Woodman GF, Kang MS, Thompson KG, Schall JD. Visual search efficiency modulates the beginning of response preparation: neurophysiological evidence for discrete flow. Psych Sci In press.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |