| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1Brain Sciences Center, Veterans Affairs Medical Center, Minneapolis Minnesota; and Departments of 2Neuroscience, 3Neurology, and 4Psychiatry and 5Cognitive Sciences Center, University of Minnesota, Minneapolis Minnesota
Submitted 4 July 2004; accepted in final form 17 July 2004
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
5 Hz. Increasing the width of the transfer function, up to a width of 700–800 ms, improved the fit of the model. Furthermore, the Rsqr of the linear model improved monotonically with the number of cells in the ensemble, saturating at 60–80% for a filter width of 700 ms. Finally, it was found that including an interaction term, which allowed the transfer function to shift with the eye position, did not improve the fit of the model. Thus ensemble neural responses in superior parietal cortex provide a high-fidelity, linear representation of hand kinematics within our task. |
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
; Sakata et al. 1973), following its known anatomical connections with primary sensory cortex (Jones and Powell 1970). Later experiments explored its relation to movements (Kalaska et al. 1983; Mountcastle et al. 1975), and arm position (Georgopoulos and Massey 1985; Georgopoulos et al. 1984) using more rigorously controlled paradigms. These experiments as well as others (Chapman et al. 1984; Seal and Commenges 1985) established the presence of both a sensory representation in SPC, which could be eliminated by cutting the dorsal roots (Seal et al. 1982), and a motor representation, which could precede the earliest electromyographic (EMG) activity (Kalaska et al. 1983).
More recent experiments in the SPC have extended these findings by characterizing the neural activity with respect to task factors including delay periods (Crammond and Kalaska 1989), dynamics (Kalaska et al. 1990), different periods of a delayed reach task (Johnson et al. 1996), and movements made in a similar direction but within different parts of space (Ferraina and Bianchi 1994; Lacquaniti et al. 1995). Several of these experiments have shown that neural activity in the SPC can be modulated by task factors other than movements, including the planned direction of a movement (Crammond and Kalaska 1989) and the onset of a target that indicates the direction of an upcoming movement (Johnson et al. 1996). Finally, work on the coordinate system in which reaches are represented in parietal cortex has shown that the preferred reach directions of cells in medial area 5 can be dependent on eye position (Buneo et al. 2002).
We trained monkeys to copy triangles, squares, trapezoids, and inverted triangles (Averbeck et al. 2003) and recorded the activity of small ensembles of neurons while the monkeys performed the task. Previous experiments have shown that the neural activity in primary motor cortex can be used to predict the hand velocity during the production of complex trajectories in a tracing task (Schwartz 1992–1994; Schwartz and Moran 1999). Our analyses extend these results to parietal cortex. Although the studies cited in the preceding text have shown that there is a sensory-motor representation of the arm in parietal cortex, the accuracy of that representation has not been examined. The use of simultaneously recorded neurons is important because noise correlations measured in most systems may decrease the amount of information contained in neural ensembles (Sompolinsky et al. 2001). Thus if sequentially recorded neurons are used to estimate hand velocity, the accuracy of the reconstruction will be overestimated. Although simpler tasks could have been used to study the parietal representation of hand velocity, the complex response properties of cortical neurons (Muir and Lemon 1983) would not guarantee that the representation of movements in a simple task would correlate with the representation of movements in a more complex task. Finally, the diversity of movements generated in the execution of the trajectories in our task generated a rich dataset on which to test our decoding models.
For our analyses, we applied both a linear model and a nonlinear model to predict the hand velocity in 25-ms bins, using the neural responses of the small ensembles. These analyses allowed us to establish how well the hand movements were represented in small ensembles of parietal cortex neurons, an issue relevant to the neural control of prosthetic devices (Taylor et al. 2002; Wessberg et al. 2000) as well as the neural coding of complex arm movements. Furthermore, by comparing linear and nonlinear models, we were able to examine the relative linearity of the representation in parietal cortex. If the hand velocity was a strongly nonlinear function of the neural responses, a nonlinear model should outperform a linear model. We found that the linear model generally outperformed the nonlinear model. This suggests a relatively linear relationship between neural activity and hand velocity. The estimate of hand velocity by the linear model was highly accurate, accounting for 60–80% of the variance of the hand velocity, with ensembles of 6–18 cells. We also found that the accuracy (Rsqr) of the reconstruction scaled monotonically with the number of cells in the ensemble, saturating for larger ensembles. In the frequency domain, the average transfer function was a low-pass filter, reflecting the fact that the movement was band limited to frequencies below 5 Hz, whereas the neural activity contained power at higher frequencies. Finally, we found that there was no effect of eye position, such that an interaction term between the neural activity and the eye position did not increase the accuracy of the prediction of the hand velocity, suggesting that in the part of the cortex from which we were recording, the representation was in a hand-centered coordinate system.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
). Animals and task
Two male rhesus macaques, Macaca mulatta (M157 and M555, 8–10 kg body wt) were used in the experiments. Care and treatment of the animals during all stages of the experiments conformed to the Principles of Laboratory Animal Care (National Institutes of Health Publication No. 86–23, revised 1995). All experimental protocols were approved by the appropriate institutional review boards.
During the experiments, monkeys were seated in a primate chair facing a rear projection screen. They used a joystick (model 541 FP, Measurement Systems, Norwalk, CT) to control a cursor on the screen. A 26-mm excursion of the joystick (1 side of the square) resulted in a 113-mm (13.4° of visual angle) excursion of the cursor on the screen. The monkeys began a trial by moving the cursor into a start hold circle presented on the left half of the screen. This initiated a wait time (WT, 1 s for M157, 2 s for M555). At the end of the WT, a shape appeared on the right half of the screen, and the monkey drew the shape presented in a copy space on the left half of the screen. Monkey M555 drew triangles and squares and monkey M157 drew triangles, squares, trapezoids, and inverted triangles. To complete a trial, the monkey had to execute a trajectory while remaining within invisible corridors. The corridors were small channels, essentially outline versions of the shapes, that constrained the trajectory appropriately. A trajectory was completed when the monkey returned to the start hold circle. If at any time during the trial the cursor went beyond the allowed corridors, the trial was terminated, and the monkey was not given a reward. Animals were allowed as much time as they needed to complete the trajectory, although they normally drew continuously and rapidly. M157 drew 30 correct trials of each shape in five blocks of six trials each. The average percent correct performance for M157 was 72.4 ± 0.12, 63.2 ± 0.14, 56.4 ± 0.10, and 62.6 ± 0.13% (means ± SD) for the triangle, square, trapezoid, and inverted triangle, respectively. M555 was given 60 trials of each shape, in one block of 60 trials each. The average percent correct performance for M555 was 55.6 ± 0.19 and 34.1 ± 0.15% for the triangle and the square, respectively.
Data collection and preprocessing
Eye position was monitored using a scleral search coil technique (CNC Engineering, Seattle, WA) (Fuchs and Robinson 1966; Judge et al. 1980). The coil and a halo for holding the head (Nakasawa Works, Tokyo, Japan) were implanted during an aseptic surgery performed under general gas (halothane) anesthesia. We recorded the activity of multiple single cells simultaneously using a 16-microelectrode recording matrix (UWE THOMAS Recording, Marburg, Germany) (see Lee et al. 1998; Mountcastle et al. 1991). The joystick and eye position were sampled at 200 Hz. The hand position, eye position, and spike times were prefiltered, using an FIR low-pass filter, with a cutoff frequency of 15 Hz (Oppenheim and Schafer 1989), and sub-sampled at 40 Hz (similar to using a binwidth of 25 ms) to reduce the computational load, and eliminate noise. The FIR filter order was sufficient to minimize aliasing.
Frequency domain analyses presented in the results were carried out using the fast Fourier transform (FFT) algorithm of Matlab (The Mathworks, Natick, MA). The periodogram for the neural signal (see Fig. 5) was estimated by averaging across trials and across the cells from the ensemble, using the spike times for each cell before they were prefiltered. The periodogram for the transfer function was found by computing the periodogram for each cell's transfer function, and then averaging these across cells. Finally, the periodogram for the predicted and measured velocity were found by averaging across trials, and across the x and y components of the velocity.
|
Nonlinear model
We carried out a nonlinear regression analysis (Chan and Tong 2001) by numerically computing the expected value of the velocity given the neural activity. Theoretically such a model represents the best estimate of a variable in a least squares sense (Papoulis 1991). The expected value can be written as
 | (1) |
(t) is the velocity to be predicted at time t, and s(t) is the neural activity vector at time t. This integral can be estimated numerically, using the formula
 | (2) |
(t) is a neural activity vector defined below and b is a bandwidth parameter, which was optimized separately for each model. The vertical bars indicate Euclidean distance.
Estimates were constructed using Eq. 2 by first splitting the dataset in half. The current neural activity vector, (t) was then compared with all the neural activity vectors (t) in the other half of the dataset, i.e., the sum over T is taken across the half of the dataset, which did not contain (t). Thus to predict the hand velocity at a point in time, the Euclidean distance was calculated, between the neural activity vector (t), and all the neural activity vectors in the other half of the data, denoted by (t). These distances were then scaled by the bandwidth parameter b and used as arguments to f, which was a zero mean, unit variance Gaussian distribution. The velocity, (t), which corresponded to each neural activity vector (t), was then weighted by f. This regression is similar to a k-nearest neighbors analysis, in which the estimated velocity would be the average of the k nearest neighbors, except in our estimate all the data points were used, and our average was weighted by the distance of each data point. In practice, however, b was quite small, such that the mean was heavily weighted toward the nearby data points, making it a local average.
Because we were analyzing data from sensory cortex, which is also known to have neural activity preceding (i.e., motor) and following (i.e., sensory) the movement in time, the window of neural responses was centered on the hand velocity. Thus when predicting the velocity at time t, the corresponding neural activity vector was given by
 |
|
We also used a linear model to predict the hand velocity from the neural activity (Ljung 1999). For the linear model, we estimated the optimal transfer function h, which, when convolved with a spike train, s(t) minimized the squared error of the estimate of hand velocity (t). In our case h(k) and (t) are two-dimensional vectors, one for the x and one for the y dimension. The convolution formula for one dimension is
 | (3) |
is a constant, w is the width of the filter in bins, n is the number of neurons, x is the x dimension of the velocity, the hxi(k) are the coefficients of the filter for the x dimension of the hand velocity, for cell i, and si(t) is the neural response at time t. We fit models at a range of different widths W, and lags, 
. For the analyses presented in Figs. 2 and 10, w was eight bins. We used a smaller window in Fig. 2 to allow direct comparison with the nonlinear model and in Fig. 10 to control the complexity of the model. For the other cases, (except Fig. 8A), we present results with 28 bins (700 ms) because filter widths between 700 and 800 ms extracted the most information. For the analyses presented in Figs. 2 and 10, was 4 bins. For the other noncausal filters, was 8 bins, and for the causal filters, was 0.
|
|
To estimate the effect of eye position on the motor representation, we fit models that included interaction terms between the eye position and the neural response. Eye position was correlated with hand velocity in our task, therefore we had to test for the significance of an interaction between neural activity and eye position, in the presence of linear terms for both members of the interaction. Therefore the first model was
 | (4) |
 | (5) |
 | (6) |
Model performance
For both the linear and the nonlinear analyses, a single model, fit across the data from all the shapes, was used to predict the hand velocity from movement onset until the end of the trial. Model performance was assessed by calculating an Rsqr statistic, in conjunction with cross validation (CV). However, CV, especially in the form of leave-one-out cross validation, tends to overfit models (Larsen and Goutte 1999). To minimize the effect of overfitting, we used twofold CV. The R2 results presented in the paper were estimated as follows. The dataset, which consisted of all correct trials for all the shapes, was split into two halves. The split was carried out such that both halves of the dataset contained half of the trials for each shape. The model parameters were estimated on the first half, and the sum of the squared error (sse) was estimated on the second half, according to
 | (7) |
 | (8) |
x is the mean velocity in the x dimension and |y is the mean velocity in the y dimension. The Rsqr was then calculated as
 | (9) |
). Unless it is assumed that the correlations between the x and y components will change for different models, however, using the determinant should not produce different results, and it would be less easy to interpret. Anatomical location of recordings
Figure 1 shows the electrode penetrations. All penetrations are shown, but data were only analyzed for penetrations anterior to the intraparietal sulcus. Most cells were recorded on the crown and upper bank of the gyrus. Although the data for M157 likely came from area 2 as well as area 5, we were unable to find significant effects of the x-y position of the recording location on the Rsqr, and therefore all data were pooled.
|
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
In Fig. 3, we plot histograms of the spike rate on the mean trajectory, as well as the corresponding single trial rasters, for a single session for four single cells. The width of the line which plots the trajectory is proportional to the mean spike rate in the corresponding time bin. These cells show activity typical of our population in that they respond similarly for segments with similar directions across the shapes. These four cells belong to an ensemble of 13 cells that we explore in detail in the following figures. The cell number in the triangle histogram plot corresponds to the cell number in Fig. 4, which shows the filter coefficients, h (see METHODS, Eq. 3), for all 13 cells from this ensemble. Figure 4A shows the filter coefficients for the cells, when considered as part of the ensemble, and B shows the filter coefficients estimated separately for each cell. When the coefficients for each cell were estimated as part of an ensemble, Eq. 3 was estimated once for all cells simultaneously. However, when the filter coefficients were estimated for each cell separately, Eq. 3 was fit separately for each cell. The coefficients in the two cases are generally different because the filter coefficients estimated by the linear analysis are modified to account for correlations between the neurons to optimize the prediction of the hand velocity (Zhang et al. 1998). This is similar to the difference between the population vector and the optimal linear estimator (Salinas and Abbott 1994). Therefore the transfer functions estimated when the cells are treated as an ensemble depend arbitrarily on the other cells that were simultaneously recorded. The main features of the transfer functions, especially of the transfer functions which have the largest amplitudes (for example, cells 1, 8, and 13), are relatively unaffected by estimating them as part of the ensemble. However, as cell 8 shows, even small differences in the velocity transfer functions can lead to position representations that look much different. These transfer functions are also closely related to the commonly estimated spike-triggered averages, in our case of hand velocity. However, our transfer functions are controlled for the autocorrelation of each neuron's response.
|
|
1 Hz. This peak, however, is clearly due to the peak in the signal. Therefore the model appears to be able to reconstruct the power in the hand velocity relatively well, at all relevant frequencies. For comparison, Fig. 6, A and B, plots the population average of the measured and predicted hand velocity as well as the noise and the SNR. Figure 6A shows the average across the entire population; Fig. 6B shows the average using only ensembles with Rsqr >0.5. The SNR for the total population average has a peak <1 Hz, whereas the ensembles with better prediction performance have a peak SNR that is more similar to the measured hand velocity and more similar to our example ensemble given in Fig. 5. Thus Fig. 6B likely provides a better estimate of the noise in our analysis because the results in A are due to the lack of signal in the smaller ensembles.
|
|
= 0 in Eq. 3), whereas the noncausal filter included neural activity that followed the movement in time by 200 ms ( = 8 in Eq. 3). This plot suggests that sensory information contributes to the representation in parietal cortex because the noncausal filter is more accurate than the causal filter. Most of the information that precedes the movement is present in the interval from 500 ms to the time of the movement. Extending the filter past 500 ms before movement onset adds only a small amount of additional predictive capability. This has to be interpreted carefully with respect to the maximum lead at which neural activity predicts hand velocity because the hand velocity is also autocorrelated, and thus a spike that predicts a feature of the trajectory 200 ms in the future will also predict features of the trajectory further into the future because the trajectory also predicts features of the trajectory into the future. The average differences in the Rsqr between the causal and noncausal filters are rather small, being <0.04 at their maximum. Figure 8B shows a plot of the Rsqr of the 700-ms noncausal filter versus the number of cells in the ensemble. This plot shows that for small ensemble sizes, the Rsqr increases approximately linearly with the number of neurons in the ensemble. However, for larger ensembles, the increase in the Rsqr saturates.
In Fig. 9, we show a plot of the eye position overlaid on the evolving copy trajectory for a single representative trial for each shape. Each column shows a snapshot of the eye-position trace superimposed on the drawing at a different point in the trial with the trial time increasing from left to right. The monkey normally made a saccade to the template, which is not shown, when it appeared, then made a saccade back to the drawing area as the trial began, and then made a sequence of small saccades to follow the progress of the drawing. Thus the eye position changes considerably through the course of a trial. Although there is a stereotyped relation between the eye and hand positions, multiple regression analyses showed that only 20% of the variance in either the x or the y hand position could be accounted for by the combined effects of the x and y components of the eye position and even less of the hand velocity could be predicted by the eye position. Thus the variability in the eye position was only weekly correlated with the hand kinematics. This regression was, however, carried out at a lag of 1 (i.e., 25 ms) between the eye and the hand position because this was the lag used in our decoding analyses (see following text and Eqs. 5 and 6). It is possible that other lags could have revealed a much larger correlation between the hand and the eye (Ariff et al. 2002), although we did not explore this possibility systematically.
|
) to test for the significance of the interaction term, we wanted to test for its effect in the presence of both linear terms. Therefore both models contain linear terms for eye position and neural activity (compare Eqs. 4 and 5). The interaction model includes an additional term for the interaction between the eye position and the neural activity. Figure 10A shows a comparison between models with and without interaction terms with both models including a linear term for eye position. Across the population, there appears to be little effect of eye position on the tuning functions. This result could be due to the fact that only a subset of the population shows the effect, and therefore carrying out this analysis at the ensemble level washes out the effect as well as the fact that the model becomes quite complex at the ensemble level. Therefore we also carried out the analysis at the single-cell level. The results are plotted in Fig. 10B. Although a few cells seem to show some effect, across the population, there is relatively little change in the variance explained. |
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
The SPC is richly interconnected within the cortical network that subserves reaching movements (Caminiti et al. 1996). On the sensory side, it receives somatosensory afferents from area 2 of sensory cortex (Jones et al. 1978), which contain information about arm position and velocity (Mountcastle and Powell 1959), and visual afferents from area MIP in the anterior bank of the intraparietal sulcus (Caminiti et al. 1996), which gets visual input from area 7m and area V6a. The SPC is thus capable of integrating information about hand position and potential targets for reaching movements. On the motor side, the SPC has connections with primary motor cortex (Caminiti et al. 1996; Jones and Powell 1970) and the supplementary motor area (Pandya and Seltzer 1982) as well as a weak, though consistently reported, connection with the caudal periprincipalis area of prefrontal cortex (Pandya and Kuypers 1969; Petrides and Pandya 1984; Yeterian and Pandya 1985). Subcortically, area 5 is primarily connected with the lateral posterior and anterior pulvinar nuclei of the thalamus; these connections being the characteristic that distinguishes area 5 from area 2, which is connected with the ventral basal nuclei (Jones et al. 1979; Yeterian and Pandya 1985). There is relatively little known about the thalamic nuclei with which area 5 connects (Steriade et al. 1997). Along with its position in the parietal frontal network involved in reaching movements, area 5 also has direct projections to the spinal cord (Coulter and Jones 1977; Galea and Darian-Smith 1994; Murray and Coulter 1981), similar to the other somatomotor areas (Dum and Strick 1996; Kuypers 1981). The cortical spinal projections from the parietal cortex differ from those from the frontal cortex in that they tend to synapse on the dorsal horn in the spinal cord gray matter as opposed to the projections from the frontal cortex, which synapse on the ventral horn (Coulter and Jones 1977). This cortical-spinal projection gives the SPC the ability to directly modulate sensory responses and indirectly modulate motor output at the level of the spinal cord.
The anatomical projections, as well as previous neurophysiological results (Ashe and Georgopoulos 1994; Georgopoulos and Massey 1985; Mountcastle and Powell 1959; Prud'homme and Kalaska 1994; Tillery et al. 1996), suggest that the SPC contains both position and velocity information. In our analyses, we estimated hand velocity from the neural responses. We do not consider the problem of estimating the relative contribution of position and velocity to parietal cortex neural responses (Paninski et al. 2004) because we are principally interested in the accuracy of the representation. Given that SPC is highly interconnected with the cortical and subcortical reaching network, it is not surprising that it contains an accurate representation of hand velocity. We found that the linearly filtered activity of small ensembles of 6–18 neurons accounted for 60–80% of the variance in hand trajectories. It seems likely that larger ensembles would allow a better prediction of the hand velocity, although our performance may have been asymptoting with only 16–18 cells. Some theoretical analyses have suggested that nonlinear models would be required to extract all of the available information from larger ensembles (Shamir and Sompolinsky 2004). Further analyses and the actual recording of larger ensembles will be required to establish the limits of the effectiveness of linear or other models. We also fit both causal and noncausal models to predict the hand velocity. The causal models were less effective than the noncausal models, presumably because they were unable to utilize the sensory information in SPC. The average decrease in Rsqr, however, was <0.04 at a filter width of 700 ms. These results suggest that neural activity in parietal cortical areas 2 and 5 could be useful in driving prosthetics devices (Pesaran et al. 2002), assuming it is not strictly sensory in origin. It is also interesting to note that task factors other than position and velocity, for example, the shape being drawn (Averbeck et al. 2003), may be affecting neural responses in area 5, and these influences will also decrease our ability to decode velocity.
The effectiveness of any decoding algorithm is dependent on two factors, how well the variable of interest is represented in the neuronal population and how effective the algorithm is at extracting the information from the neuronal signal. Linear models for decoding neural activity have proved successful in a number of different situations, including motion-sensitive neurons in the fly (Bialek et al. 1991), responses of LGN neurons to movies (Stanley et al. 1999), V1 neural responses to drifting gratings (Jones and Palmer 1987), the representation of reaching movements in motor cortex (Paninski et al. 2004; Wessberg et al. 2000), and the transfer of information from synaptic input to output (Eckhorn et al. 1976). Eckhorn et al. (1976) further showed that linear decoding was able to extract between 90 and 99% of the total information in the input signal to a cell. Furthermore, Bialek and his colleagues (1993) have suggested, based on empirical observations and a model of neural firing rates, that if mean neural firing rates are on the scale of the autocorrelation interval of the signal to be estimated, it will not be possible to do better than a linear model. Analyses of the position representation in S1 have shown that hand position is often encoded nonlinearly (Tillery et al. 1996). However, the movements made in our task were relatively small (a few centimeters). Thus even if this system is generally nonlinear, the linear approximation may be valid over the range of movements tested.
We have also explored the effect of eye position on the transfer function of our neurons in the SPC. We found that allowing the transfer functions to shift with eye position, in general, did not improve the ability of our model to predict the hand velocity. This differs from results reported in the parietal reach region (Batista et al. 1999), as well as area 5 (Buneo et al. 2002). There are several reasons that can account for this difference, including differences in the task and the analytical techniques employed as well as the exact areas from which the respective studies recorded. Perhaps the largest difference is that our recordings were located relatively anterior and lateral, and published results (Buneo et al. 2002) suggest that the effect of the eye position on the movement field of a neuron decreases at more anterior locations in the cortex.
Thus we have shown that a simple linear model can capture much of the useful signal in simultaneously recorded neural responses. In some respects, the linear model can be considered an extension of a rate code because instead of counting spikes within a window, which amounts to giving each spike in the window the same weight, we simply compute a weighted sum of the spikes within a window to get our estimate. Furthermore, this model was able to extract up 80% of the variance of a complex hand movement, using small ensembles of neurons. Thus parietal cortex contains a robust neural representation of hand velocity.
|
GRANTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: A. P. Georgopoulos, Brain Sciences Center (11B), Veterans Affairs Medical Center, One Veterans Drive, Minneapolis, MN 55417 (E-mail: omega{at}umn.edu)
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
Ashe J and Georgopoulos AP. Movement parameters and neural activity in motor cortex and area 5. Cereb Cortex 4: 590–600, 1994.
Averbeck BB, Chafee MV, Crowe DA, and Georgopoulos AP. Parallel processing of serial movements in prefrontal cortex. Proc Natl Acad Sci USA 99: 13172–13177, 2002.
Averbeck BB, Chafee MV, Crowe DA, and Georgopoulos AP. Neural activity in prefrontal cortex during copying geometrical shapes. I. Single cells encode shape, sequence, and metric parameters. Exp Brain Res 150: 127–141, 2003.[ISI][Medline]
Batista AP, Buneo CA, Snyder LH, and Andersen RA. Reach plans in eye-centered coordinates. Science 285: 257–260, 1999.
Bialek W, DeWeese M, Rieke F, and Warland D. Bits and Brains: Information flow in the nervous system. Physica A 200: 581–593, 1993.[CrossRef]
Bialek W, Rieke F, de Ruyter van Steveninck RR, and Warland D. Reading a neural code. Science 252: 1854–1857, 1991.
Buneo CA, Jarvis MR, Batista AP, and Andersen RA. Direct visuomotor transformations for reaching. Nature 416: 632–636, 2002.[CrossRef][Medline]
Caminiti R, Ferraina S, and Johnson PB. The sources of visual information to the primate frontal lobe: a novel role for the superior parietal lobule. Cereb Cortex 6: 319–328, 1996.
Chan K-S and Tong H. Chaos: A Statistical Perspective. New York: Springer, 2001.
Chapman CE, Spidalieri G, and Lamarre Y. Discharge properties of area 5 neurons during arm movements triggered by sensory stimuli in the monkey. Brain Res 309: 63–77, 1984.[CrossRef][ISI][Medline]
Coulter JD and Jones EG. Differential distribution of corticospinal projections from individual cytoarchitectonic fields in the monkey. Brain Res 129: 335–340, 1977.[CrossRef][ISI][Medline]
Crammond DJ and Kalaska JF. Neuronal activity in primate parietal cortex area 5 varies with intended movement direction during an instructed-delay period. Exp Brain Res 76: 458–462, 1989.[ISI][Medline]
Draper NR and Smith H. Applied Regression Analysis. New York: Wiley, 1998.
Duffy FH and Burchfiel JL. Somatosensory system: organizational hierarchy from single units in monkey area 5. Science 172: 273–275, 1971.
Dum RP and Strick PL. The corticospinal system: a structural framework for the central control of movement. In: Handbook of Physiology. Exercise: Regulation and Integration of Multiple Systems, edited by Rowell LB and Shepherd JT. New York: Oxford, 1996, sect. 12, p. 217–254.
Eckhorn R, Grusser OJ, Kroller J, Pellnitz K, and Popel B. Efficiency of different neuronal codes: information transfer calculations for three different neuronal systems. Biol Cybern 22: 49–60, 1976.[CrossRef][ISI][Medline]
Ferraina S and Bianchi L. Posterior parietal cortex: functional properties of neurons in area 5 during an instructed-delay reaching task within different parts of space. Exp Brain Res 99: 175–178, 1994.[ISI][Medline]
Fuchs AF and Robinson DA. A method for measuring horizontal and vertical eye movement chronically in the monkey. J Appl Physiol 21: 1068–1070, 1966.
Galea MP and Darian-Smith I. Multiple corticospinal neuron populations in the macaque monkey are specified by their unique cortical origins, spinal terminations, and connections. Cereb Cortex 4: 166–194, 1994.
Georgopoulos AP, Caminiti R, and Kalaska JF. Static spatial effects in motor cortex and area 5: quantitative relations in a two-dimensional space. Exp Brain Res 54: 446–454, 1984.[ISI][Medline]
Georgopoulos AP and Massey JT. Static versus dynamic effects in motor cortex and area 5: comparison during movement time. Behav Brain Res 18: 159–166, 1985.[CrossRef][ISI][Medline]
Johnson PB, Ferraina S, Bianchi L, and Caminiti R. Cortical networks for visual reaching: physiological and anatomical organization of frontal and parietal lobe arm regions. Cereb Cortex 6: 102–119, 1996.
Johnson RA and Wichern DW. Applied Multivariate Statistical Analysis. Saddle River, NJ: Prentice Hall, 1998.
Jones EG, Coulter JD, and Hendry SH. Intracortical connectivity of architectonic fields in the somatic sensory, motor and parietal cortex of monkeys. J Comp Neurol 181: 291–347, 1978.[CrossRef][ISI][Medline]
Jones EG and Powell TP. An anatomical study of converging sensory pathways within the cerebral cortex of the monkey. Brain 93: 793–820, 1970.
Jones EG, Wise SP, and Coulter JD. Differential thalamic relationships of sensory-motor and parietal cortical fields in monkeys. J Comp Neurol 183: 833–882, 1979.[CrossRef][ISI][Medline]
Jones JP and Palmer LA. The two-dimensional spatial structure of simple receptive fields in cat striate cortex. J Neurophysiol 58: 1187–1211, 1987.
Judge SJ, Richmond BJ, and Chu FC. Implantation of magnetic search coils for measurement of eye position: an improved method. Vision Res 20: 535–538, 1980.[CrossRef][ISI][Medline]
Kalaska JF, Caminiti R, and Georgopoulos AP. Cortical mechanisms related to the direction of two-dimensional arm movements: relations in parietal area 5 and comparison with motor cortex. Exp Brain Res 51: 247–260, 1983.[ISI][Medline]
Kalaska JF, Cohen DA, Prud'homme M, and Hyde ML. Parietal area 5 neuronal activity encodes movement kinematics, not movement dynamics. Exp Brain Res 80: 351–364, 1990.[ISI][Medline]
Kuypers HGJM. Anatomy of the descending pathways. In: Handbook of Physiology. The Nervous System. Motor Control. Bethesda MD: Am. Physiol. Soc., 1981, sect. 1, vol. II, p. 597–666.
Lacquaniti F, Guigon E, Bianchi L, Ferraina S, and Caminiti R. Representing spatial information for limb movement: role of area 5 in the monkey. Cereb Cortex 5: 391–409, 1995.
Larsen J and Goutte C. On optimal data split for generalization estimation and model selection. In: Proceedings of the IEEE Workshop on Neural Networks for Signal Processing IX, edited by Hu YH, Larsen J, Wilson E, and Douglas S. Piscataway, NJ: IEEE, 1999, p. 225–234.
Lee D, Port NL, Kruse W, and Georgopoulos AP. Neuronal population coding: multielectrode recordings in primate cerebral cortex. In: Neuronal Ensembles: Strategies for Recording and Decoding. New York: Wiley, 1998, p. 117–136.
Ljung L. System Identification, Theory for the User. Upper Saddle River, NJ: Prentice Hall, 1999.
Mountcastle VB, Lynch JC, Georgopoulos A, Sakata H, and Acuna C. Posterior parietal association cortex of the monkey: command functions for operations within extrapersonal space. J Neurophysiol 38: 871–908, 1975.
Mountcastle VB and Powell TPS. Central nervous mechanisms subserving position sense and kinesthesis. Bull Johns Hopkins Hosp 105: 173–200, 1959.[ISI][Medline]
Mountcastle VB, Reitboeck HJ, Poggio GF, and Steinmetz MA. Adaptation of the Reitboeck method of multiple microelectrode recording to the neocortex of the waking monkey. J Neurosci Methods 36: 77–84, 1991.[CrossRef][ISI][Medline]
Muir RB and Lemon RN. Corticospinal neurons with a special role in precision grip. Brain Res 261: 312–316, 1983.[CrossRef][ISI][Medline]
Murray EA and Coulter JD. Organization of corticospinal neurons in the monkey. J Comp Neurol 195: 339–365, 1981.[CrossRef][ISI][Medline]
Oppenheim A and Schafer R. Discrete-Time Signal Processing. Englewood Cliffs, NJ: Prentice Hall, 1989.
Pandya DN and Kuypers HG. Cortico-cortical connections in the rhesus monkey. Brain Res 13: 13–36, 1969.[CrossRef][ISI][Medline]
Pandya DN and Seltzer B. Intrinsic connections and architectonics of posterior parietal cortex in the rhesus monkey. J Comp Neurol 204: 196–210, 1982.[CrossRef][ISI][Medline]
Paninski L, Fellows MR, Hatsopoulos NG, and Donoghue JP. Spatiotemporal tuning of motor cortical neurons for hand position and velocity. J Neurophysiol 91: 515–532, 2004.
Papoulis A. Probability, Random Variables and Stochastic Processes. New York: McGraw-Hill Higher Education, 1991.
Pesaran B, Pezaris JS, Sahani M, Mitra PP, and Andersen RA. Temporal structure in neuronal activity during working memory in macaque parietal cortex. Nat Neurosci 5: 805–811, 2002.[CrossRef][ISI][Medline]
Petrides M and Pandya DN. Projections to the frontal cortex from the posterior parietal region in the rhesus monkey. J Comp Neurol 228: 105–116, 1984.[CrossRef][ISI][Medline]
Prud'homme MJ, and Kalaska JF. Proprioceptive activity in primate primary somatosensory cortex during active arm reaching movements. J Neurophysiol 72: 2280–2301, 1994.
Sakata H, Takaoka Y, Kawarasaki A, and Shibutani H. Somatosensory properties of neurons in the superior parietal cortex (area 5) of the rhesus monkey. Brain Res 64: 85–102, 1973.[CrossRef][ISI][Medline]
Salinas E and Abbott LF. Vector reconstruction from firing rates. J Comput Neurosci 1: 89–107, 1994.[CrossRef][Medline]
Schwartz AB. Motor cortical activity during drawing movements: single-unit activity during sinusoid tracing. J Neurophysiol 68: 528–541, 1992.
Schwartz AB. Motor cortical activity during drawing movements: population representation during sinusoid tracing. J Neurophysiol 70: 28–36, 1993.
Schwartz AB. Direct cortical representation of drawing. Science 265: 540–542, 1994.
