Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2010; 5(9): e12574.
Published online 2010 September 3. doi:  10.1371/journal.pone.0012574
PMCID: PMC2933244

Incorporating Prediction in Models for Two-Dimensional Smooth Pursuit

Paul L. Gribble, Editor


A predictive component can contribute to the command signal for smooth pursuit. This is readily demonstrated by the fact that low frequency sinusoidal target motion can be tracked with zero time delay or even with a small lead. The objective of this study was to characterize the predictive contributions to pursuit tracking more precisely by developing analytical models for predictive smooth pursuit. Subjects tracked a small target moving in two dimensions. In the simplest case, the periodic target motion was composed of the sums of two sinusoidal motions (SS), along both the horizontal and the vertical axes. Motions following the same or similar paths, but having a richer spectral composition, were produced by having the target follow the same path but at a constant speed (CS), and by combining the horizontal SS velocity with the vertical CS velocity and vice versa. Several different quantitative models were evaluated. The predictive contribution to the eye tracking command signal could be modeled as a low-pass filtered target acceleration signal with a time delay. This predictive signal, when combined with retinal image velocity at the same time delay, as in classical models for the initiation of pursuit, gave a good fit to the data. The weighting of the predictive acceleration component was different in different experimental conditions, being largest when target motion was simplest, following the SS velocity profiles.


Interacting with a moving target involves prediction, whether the target be pursued by means of eye movements or intercepted by a movement of the hand. Thus, for manual interception, the initial direction of the hand is in advance of the target's location [1][5]. Similarly, when the direction and the time of the onset of target motion is predictable, smooth pursuit eye movements anticipate the target's motion [6][7] in contrast to the latency of 100 ms or more when the motion onset is not predictable [8]. The fact that pursuit is maintained, albeit at reduced gain, even when the target disappears transiently [9][11] also provides evidence for predictive mechanisms.

In the examples cited above, the target speed for smooth pursuit was always constant. However, anticipatory effects in pursuit can also be demonstrated when target motion varies predictively, for example when it is generated by a combination of sine waves of different frequencies [12][15]. In that instance, eye velocity leads target velocity at low frequencies (below 0.5 Hz), whereas it lags target velocity at higher frequencies (as would be expected if eye velocity were delayed with respect to target velocity). Furthermore, the amount of phase lead depends on the spectral composition of the target signal. For example, the addition of a harmonic component to a fundamental frequency alters the amplitude and timing of the response to the fundamental component [14]. This suggests that the nature or the extent of predictive mechanisms depends on the precise characteristics of the stimulus.

Thus the aim of the present experiments was to characterize more precisely the predictive component of pursuit eye movements by developing a quantitative model. Quantitative models have been successful in predicting the initiation of eye movements in response to a step of target velocity [16][18]. In all of these models, pursuit eye velocity is driven in part by retinal image velocity (the difference between target and eye velocity) with a time delay. This signal is modified in different ways, either by means of internal feedback loops of eye velocity [18] or acceleration [17] or by a sensitivity to image acceleration [16]. The latter model was also able to account for the effect of high frequency (2–10 Hz) perturbations during maintained pursuit [19]. However, these models do not account for anticipatory properties of pursuit eye movements. Here we show that modifying the Krauzlis and Lisberger [16] model by adding a low-pass filtered target acceleration signal to image velocity is able to reproduce the major aspects of two dimensional smooth pursuit of periodic target motion, and that the gain of this predictive signal depends on the characteristics of the target motion.


Six subjects participated in this experiment. All had normal or corrected to normal vision and gave informed written consent to procedures approved by the Institutional Review Board of the University of Minnesota.

Target motions

Subjects tracked a cyan, circular target 0.5° in diameter that was displayed on a computer monitor. The target underwent periodic motion in two dimensions with a period of 4.5 seconds. Target trajectories were constructed by a sum of sines (fundamental and 2nd or 3rd harmonic) in the x (horizontal) and y (vertical) directions. Five such target paths were constructed, illustrated in Fig. 1. These base trajectories (SS – sum of sines) were modified such that the target followed the same path with the same period (4.5 s) but at a constant speed (CS) [4]. This transformation added additional harmonic frequency components to the x and y components of the target velocities, in that sense making target motion less predictable. Finally, we constructed two additional target trajectories by combining the x SS velocity with the y CS velocity and vice versa. These combinations resulted in target paths that were modified slightly from the original paths, as shown for one example (path 5) in the lower right panel of Fig. 1. However, on the assumption that smooth pursuit eye velocity depends primarily on target velocity and its derivatives and not on position, these additional trajectories permitted us to assess the extent to which the x – and y – components of pursuit velocity are independent of each other. If they are, then the horizontal and vertical eye velocity components of these combined stimuli should be identical to their respective counterparts in the SS and CS trajectories.

Figure 1
Target motion paths.

Experimental procedures and data analysis

The experimental procedures have been described in detail previously [1], [3]. Subjects sat with their eyes 40 cm from a computer monitor with the head stabilized by a chin rest. At the beginning of each trial the target appeared at the starting location of the particular trajectory, and after a brief interval began moving with the prescribed motion profile. As described above, there were 5 different paths and 4 motion profiles for each path, resulting in 20 combinations. Each was presented 10 times, in random order, for a total of 200 trials. Each trial lasted 5.7 seconds, i.e., 1.25 times the period of the motion. Trials in which pursuit was not maintained or in which there was an eye blink were rejected during the experimental session and subsequently repeated. Subjects could take breaks during the course of the experiment and each experiment lasted about 1.5 hours.

Eye movements were recorded using head-mounted infrared cameras (EyeLink, SR Research, Mississauga, Ontario) at 250 Hz. To increase the signal to noise ratio, we generally combined recordings from the left and right eyes. Position data were first smoothed with a double-sided exponential filter (time constant 4 ms) and differentiated numerically. Saccades were identified and removed by interpolation with a cubic spline [1]. Desaccaded velocity traces for the 10 trials for each experimental condition were then averaged and all subsequent analysis was performed on these averaged data.

Eye velocities were first analyzed in the frequency domain. We restricted our analysis to the steady state response, neglecting eye movements during the first second of target motion. For this purpose, we resampled eye or target velocities to generate 512 equally spaced points in the interval from 1.0 to 5.5 seconds and transformed the data into the frequency domain using the fast Fourier transform (fft). From this analysis, we computed gains and phases of the response at various frequencies. Average responses from the six subjects were computed, using circular statistics for the phase data [20] and the Rayleigh test for significance with p<0.01.

We also tested several different models relating smooth pursuit eye velocity to target motion. These models will be described in detail in Results. For each model, we solved the differential equations relating input and output using the Runge-Kutta method with adaptive step size [21] and found the parameters providing the best fit of the model to the data by minimizing the square error between model velocity and eye velocity in the interval from 1.0 to 5.5 seconds, i.e. neglecting the onset of pursuit. For this purpose, we used the simplex algorithm of Nelder and Mead [22].


Figure 2 illustrates representative results from one subject (4) and one path (2), each of the 4 panels corresponding to one of the motion profiles: A is sum of sines (SS), B is constant speed (CS), and C and D are combinations of SS and CS for the x – and y – components of the target motion, generated from the same component velocity profiles as in A and B. Each panel shows from bottom to top the x– and y – velocities, the speed, and the direction of motion, positive being up and to the right and direction being measured in the counterclockwise direction from the right horizontal. Target motion is shown by the heavy black trace and the blue traces show the mean (±1 SE) of the ocular response.

Figure 2
Representative results for one target path (2) and one subject (4).

From inspection of the traces, it is apparent that pursuit was initiated at a latency of about 100 ms. Thereafter, for the sinusoidal motion (Fig. 2A) pursuit velocity matched target velocity with a gain that was close to unity, but with a slight delay. These delays are more apparent for the CS profile in Fig. 2B. Also, note that transforming the trajectory so that it moved at constant speed introduced higher frequency components to the x – and y – target motion profiles and that the ocular response appears to reflect a low pass filtering of the target motion. Inspection of Fig. 2B may suggest that the amount by which eye velocity lagged the target's velocity was not constant, but varied over time (see for example the x velocity at ~3.5 s). However, the results of modeling to be presented in a subsequent section do not support this supposition. Finally, even though the target's speed was constant in Fig. 2B, eye speed was modulated substantially, decreasing when the direction of target motion changed abruptly (i.e. also at ~3.5 s).

The responses in x – and y – velocity to combined (SS-CS) stimuli were largely similar to the component responses to SS and CS target motions, as can be ascertained for example by comparing the x – velocities in Fig. 2C and 2A, and the y – velocities in Fig. 2C and 2B. Nevertheless, there were instances where there were clear differences, for example in the y – velocity component around 2.0 s in Fig. 2D and 2A, and the x – velocity component at about 3.5 s in Fig. 2D and 2B, suggesting an interaction between the horizontal and vertical components of pursuit. As we will show in a following section, some of these differences, although small, were found consistently in the responses of the six subjects.

Frequency response of pursuit eye velocity

When the target velocity consisted of a sum of sines (a fundamental component at 0.22 Hz and a 2nd or 3rd harmonic), smooth pursuit velocity led target velocity at the fundamental frequency, but lagged it at higher frequencies (Table 1). The gain of the harmonic responses was greater than the gain of the fundamental and the horizontal gain was greater than the gain in the vertical direction. All of these results are in agreement with previous observations [1], [12], [14], [23]. An ANOVA on the gain and phase at each frequency using speed profile and path as factors showed that, with the exception of the phase of the y-velocity response at 0.22 Hz (see Table 1), neither the phase nor the gain depended on the speed profile (SS or CS) of the other directional component (p>0.05). In several instances, the gain and/or the phase of the smooth pursuit did depend on the path (F4, 50>4.36, p<0.01).

Table 1
Frequency Response to Sum of Sines Stimuli.

There was power at frequency components not contained in the target signals, but this was generally small, never exceeding 3.4°/s, compared to peak target velocities that ranged from 30 to 40°/s (see Fig. 2). The largest non-target related frequency responses were at the 4th and 5th harmonics. At these frequencies, pursuit response amplitude averaged 4.3% (0.89 Hz) and 6.1% (1.11 Hz) of the amplitude of the largest target-related frequency component. Although small, some of these responses could not be attributed to random noise because sometimes the phase of the response for a given frequency and path was not randomly distributed for the six subjects (Rayleigh test with p<0.01). This was the case for 75% of the instances at 1.11 Hz, but much less common at all other frequencies in the range of 0.89 to 3.33 Hz (maximum 40% at 1.33 Hz, average 20%). Furthermore, when the phases were not randomly distributed they generally (in 75% of the instances at 1.11 Hz) did not differ significantly when the velocity along the other coordinate (x or y) was generated from a sum of spines (SS) or a constant speed (CS) trajectory. Thus it is unlikely that these distortions are attributable to interactions between the x- and y- components of pursuit velocity.

When the target velocity was generated according to the CS criterion, its power was distributed throughout the frequency spectrum, being highest at the frequency components for the corresponding SS trajectory. This is illustrated in Fig. 3 for two of the paths (3 and 4). In each panel, the black trace in the lower plot depicts the amplitude at each of the frequency components up to 3.3 Hz on a logarithmic scale. For example, for the x – component of Path 3, the SS trajectory was generated from the fundamental (0.22 Hz) and the 3rd harmonic (0.67 Hz), whereas the y – component contained power at the fundamental and the 2nd harmonic. Generally, power in the target signal decreased with frequency, albeit in a non-monotonic fashion. For those frequencies in which there was consistent response (Rayleigh test on the phase, p<0.01), the plots also depict the amplitude of the pursuit eye velocity, its gain and its phase. For each panel, two traces are shown, coded according to whether the target velocity of the other component corresponded to the SS trajectory (cyan) or the CS trajectory (red).

Figure 3
Frequency analysis of response to CS stimuli.

Comparing the results for the two paths and the two directions (x – and y – axes), one notes that the phases are quite consistent, independent of path and component, whereas the gains are considerably more variable, within a particular path and direction (see error bars denoting the standard deviation across subjects) and also between conditions. Specifically, in the frequency range from 0.22 to 3.33 Hz, the phase decreased monotonically from a slight phase lead at 0.22 Hz to a lag of ~180° at 3.33 Hz. Furthermore, the phase of the response was quite similar, irrespective of whether the other component of the target's velocity followed the SS or the CS trajectory. The only exception to this in Fig. 3 is for the y – component for Path 3 (top left panel), where the phase at 1.56 Hz for the SS trajectory (cyan) was close to zero. (The phase for the CS trajectory at this frequency is not shown, because this value was not consistent from subject to subject, p>0.01). The results for the two paths shown in Fig. 3 are representative of the results for all 5 paths (see Table 2). A statistical comparison of the phases of the response showed that in only a small fraction of the cases (7%, 9/128, paired t-test for each path and frequency, p<0.01) did this value depend on the target motion of the other coordinate.

Table 2
Frequency Response to Constant Speed Stimuli.

As noted before, the gains of the responses were more variable. In general, the gain declined gradually with frequency (for example, the x – direction for path 3 and the y – direction for path 4), but there were also instances where the gain remained elevated or increased with frequency (y – direction, path 3). On average the gain, which was 0.8 to 0.9 at the lowest frequencies, decreased to a value of about 0.5 at 3 Hz. At the lowest 3 frequencies, the gains for the CS condition tended to decrease slightly in contrast to the increase found when the target followed a sum of sines (SS) trajectory and the phase lag at 0.67 Hz tended to be greater (−21.3° on average) compared to the phase lag at that frequency for the SS trajectories (−11.3° on average).

Interactions between horizontal and vertical components of pursuit

As described above, the analysis of the results in the frequency domain did not provide strong evidence for interactions between pursuit along the two coordinate axes in the sense that the phase of the responses generally did not depend on the target motion in the other dimension. To further examine this issue, we also analyzed the data in the time domain, computing the difference in pursuit velocity along one dimension (e.g. x-axis) for the two different target motions along the other axis (e.g. y-axis). We first eliminated any effects arising from potential differences in the gain of pursuit in the two conditions by scaling. This effect was modest, the slopes of the regressions between the two conditions averaging 0.93±0.05 for SS and 0.98±0.33 for CS. Since the coefficients of determination were also uniformly high (r2>0.944 in all cases, averaging 0.982), differences in pursuit eye velocity along one axis arising from different target velocities along the other axis can be expected to be minor.

Nevertheless, as shown in Figs. 4 and and5,5, occasionally there were differences in pursuit velocity that were consistent from subject to subject. Fig. 4 shows the difference in pursuit velocities for path 5, results for the x-direction being shown in the left column and for the y-direction in the right column. The two target velocities along the other coordinate direction are shown in the overlay in the middle or lower trace of each panel, and the top traces show the difference in pursuit eye velocity for the two conditions for each of the six subjects. Note that this difference plot has been scaled by a factor of 5 compared to the target velocity.

Figure 4
Degree of independence of pursuit motion along the horizontal and vertical directions.
Figure 5
Interactions between horizontal and vertical eye velocities.

While the differences are generally small and variable from subject to subject, occasionally, such as in the interval from 1 to 2 s for the SS Vx (top left panel) and for the CS Vy at 1 s (bottom right panel), there were consistent differences in the results. This can also be seen in Fig. 5 which shows the mean (± SE) for the results for a different path (path 2) plotted in the same format. Time intervals in which the difference in velocities was significant (t-test, p<0.01) are indicated above the traces of the mean velocities. For this path, these instances were sparse and this was the case for all 5 paths; overall, the difference in velocities was significantly different from zero 8.6% of the time at the p<0.01 level and 2.8% of the time at the p<0.001 level.

We were unable to relate the differences in the response to the differences in the time course of target motion along the other axis, such as differences in target velocity or acceleration. A frequency analysis of the differences in pursuit velocity showed that the power in the difference was uniformly distributed in the interval from 0.22 to 2.0 Hz, with an amplitude of ~0.75°/s, decreasing at frequencies above 2 Hz. Accordingly at the lowest frequencies (<0.67 Hz), the difference signal was small compared to the response to the target motion (<8%), whereas at higher frequencies the ratio was much greater (e.g. 48% at 2 Hz).

Quantitative models of prediction in smooth pursuit

A flat gain and a phase lag that increases with frequency implies a pure time delay between the input and the output. Thus, the frequency responses described in Fig. 3 and Tables 1 and and22 suggest the presence of a time delay between target motion and eye velocity. This conclusion is hardly surprising because it is consistent with a large body of evidence. However, at the lowest frequency (0.22 Hz), eye velocity actually led target velocity irrespective of whether target motion only consisted of a few frequency components (SS) or whether the spectral composition of the motion was more complex. This observation suggests the presence of a predictive component in the control of a smooth pursuit and we tested a variety of models in an attempt to characterize this component more precisely.

We began with a simple model in which eye acceleration is proportional to the difference of target velocity (vt) and eye velocity (ve) , each delayed in time:

equation image

where τt and τe are the time delays of eye acceleration dve/dt relative to the target and eye velocity, respectively. The parameter g (gain) is included because the gain of smooth pursuit was typically less than unity. Finally, the parameter a defines the time constant of the response, having units of s−1. This model is a simplified version of one developed by Krauzlis and Lisberger [16]. Their model also included a term proportional to retinal image acceleration, which was considerably smaller, and a transient response triggered at the onset of target motion. Since we were interested primarily in the steady-state response, this second term could be neglected. If the two time delays τt and τe are equal and g = 1, the model implies that eye acceleration is proportional to the retinal image velocity. However, to allow the model more flexibility, we permitted the two time delays to differ.

We fitted the model to the data for the constant speed (CS) condition, fitting the x- and y-components of the velocity separately. For each subject (6) and path (5), we identified the 4 parameters that minimized the square error between the model and the actual eye velocity over the interval from 1 to 5.5 seconds. That is, we neglected the first second of the response in the fitting. Overall, the model gave a good fit to the data, the variance not accounted for (VNAF) averaging 9.0±3.0% for the fits to the 30 trial averages. Moreover, the parameters that gave the best fit were also quite consistent for all subjects and paths. The parameter a averaged 6.2±1.3 s−1, corresponding to a time constant of 160 ms. This is somewhat larger than the value of ~10 s−1 obtained by Krauzlis and Lisberger [16]. The estimates for the time delays of target and eye velocity were also consistent, but differed significantly, with an estimate of 120±30 ms for the delay in eye velocity (τe) compared to a value of 20±20 ms for target velocity (τt). While the value for τe is consistent with the latency for the initiation for smooth pursuit, the value for τt clearly is much smaller. Finally, the average gain (g) in the x-direction (0.73±0.11) was slightly larger than the value (0.62±0.10) in the y-direction (see Table 2).

Not surprisingly, this model gave an even better fit to the results for the SS condition with an average VNAF of 2.9%. The estimates for the two time delays (100±20 ms for τe and 40±20 ms for τt) were comparable to the values obtained for the CS target trajectories, as were the values for the gain g and the time constant a.

The fact that the time delay for target velocity was much smaller than the time delay for eye velocity and also much less than the latency for the initiation of smooth pursuit suggests that there was an additional, predictive signal providing for a phase advance of the target velocity. One means of achieving this would be to add a component proportional to target acceleration, since by Taylor's theorem,

equation image

Moreover, the frequency analysis suggested that that such a predictive component would be most important at low frequencies (Fig. 3 and Table 2). One means of achieving this would be to low-pass filter target acceleration:

equation image

where An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e004.jpg is the target acceleration and An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e005.jpg is a predictive version of this signal, obtained by low-pass filtering with a time constant equal to 1/b. Thus, we added a predictive target acceleration component to the original model (eq. 1), and constrained the time delays for target and eye velocity to be the same:

equation image

The tangential acceleration (An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e007.jpg) is equal to the rate of change in speed whereas the normal component (An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e008.jpg) is proportional to the curvature or the rate of change in direction. Since visual sensitivity to these two components may differ, we permitted the weightings for the two predictive acceleration components to differ [24]. Thus, this model had one more free parameter than did the simpler one (eq. 1).

We first fitted the new model to all of the data (5 paths and 4 conditions) for each of the subjects, assuming that that the gains in the x- and y-direction (gx and gy) could differ, minimizing the squared error between the model and the actual data over one cycle (from 1.0 to 5.5 s). This procedure gave a good fit to the data, with an average variance not accounted for equal to 7.43%. However, the estimates for the time delay τ were very variable (ranging from 0 to 70 ms), as were the estimates for the time constant b (eq. 3), which ranged from 3.4 to 7.8 s−1. This was not unexpected; according to eq. 2 changes in time delay can be compensated for by changing the weighting (and filtering) of the acceleration term. We explored the effect of fixing the time delay (τ) and found that fixing it at 80 ms, a plausible value given neural time delays in the feedback loop, gave only a modest 4% increase in the error (from 7.43 to 7.72%). Therefore, we fixed the time delay, and found the other parameter values that gave a best fit to the data (see Table 3). With this restriction, the estimates for the other parameters were quite consistent across the 6 subjects. The average value of the parameter b was 3.47, corresponding to a time constant of ~300 ms for the predictive acceleration term and the average value for a was 7.12 , corresponding to a main time constant of 140 ms. This latter value is comparable to the value of 160 ms found using the simpler model. The predictive feedback weighting terms c1 and c2 were similar to one another, 0.27 and 0.29 respectively. Finally, the gain in the x direction (gx = 0.53) was somewhat larger than the value for the y direction (gy = 0.43) as in the simpler model. However, both values were smaller than the values found with the simpler model but this is to be expected, since the second model had additional acceleration terms that added to the target motion signal.

Table 3
Best Fit Parameter Values.

While one might expect that the two time constants (a and b) would be consistent from trial to trial for a given subject, one might also expect that the gains of the feedback of target velocity (g) and the amount of the predictive contribution (c) could depend on the speed profile, in the sense that a simpler velocity profile (SS) would be more predictable. Therefore, we repeated the modeling, fixing the two time constants and the time delay τ at the values in Table 3 and obtaining the values for the other 4 parameters that gave the best fit for each of the 20 combinations of path and speed profile. This gave a substantially better fit to the data, the average VNAF decreasing by 23% from 7.72% to 5.96%. Some examples of this fitting procedure are shown in Figs. 6 and and77 for two different subjects (1 in Fig. 6 and 6 in Fig. 7) and two different paths, 2 and 5 respectively. In each panel, the x- and y- components of the eye velocity are shown in cyan and the predictions of the model are superimposed in red. For these two examples, the average VNAF was 5.3% (Fig. 6) and 4.6% (Fig.7), the worst fit being for the y-velocity components in Fig. 6B (10.8%) and Fig. 7C (8.0%). Note that the model was able to reproduce well most of the fluctuations in pursuit eye velocity, some exceptions occurring for the y-velocity (Fig. 6B and 6C, interval 1.5 to 2.0 s) and the x-velocity in Fig. 7B (at about 1.5 s). Recall that we did not include the first 1.0 s after motion onset in the fitting procedure, and occasionally this initial interval was also not well-fit by the model.

Figure 6
Fit of the predictive model to the experimental data.
Figure 7
Fit of the predictive model to experimental data for another path (5) and subject (6).

In two instances, both involving the CS velocity profiles, there was a discrepancy between target motion and eye velocity in the sense that eye velocity decayed to zero while target motion increased. In those instances, the model gave a poor fit to the data. One is shown in Fig. 8, the two arrows denoting times where the model's prediction deviated substantially from the measured responses. These results are for path 1, and in both instances, the y-component of target velocity was large and negative, whereas the y-component of pursuit velocity decayed to zero. These intervals correspond to times at which the target approached the lower boundary of the monitor at constant speed (see Fig. 1), before reversing. A similar discrepancy was also observed for the x-velocity for path 3 as the target approached the right-hand border (data not shown). These discrepancies most likely reflect the effect of cognitive influences on tracking behavior, subjects expecting with some confidence that the target would reverse direction as it approached the edge of the screen. The model was not intended to capture such effects, which have also been observed previously [25].

Figure 8
Instance in which the model gave a poor fit to the data.

The two weighting coefficients for the predictive acceleration component (c1 and c2) did depend on the speed profile (ANOVA, F3,116>12.4, p<0.001), but they did not differ from each other (paired t-test, p = 0.18). The coefficients were largest for the SS condition (averaging 0.39) and they were smallest for the CS condition (averaging 0.21), with intermediate values for the two combined CS-SS conditions (0.29 on average). The values for the two combined conditions (SSx-CSy and CSx-SSy) did not differ from each other. Thus, as was expected, the weighting of the predictive acceleration component in eq. 4 became greater when the target motion was more predictable, in the sense that its power spectrum was more confined. However, the components of acceleration tangential and normal to the direction of motion were equal.

Fig. 9 illustrates the effect of the predictive acceleration component in the frequency domain, with three different values for c (ranging from 0.0 to 0.5). (The gain was normalized in these plots so that the maximum value was always equal to unity.) Note that the effect of this predictive component is most prominent at low frequencies, introducing a phase advance but decreasing the relative gain at frequencies below 1.0 Hz. At higher frequencies, the effect of this predictive component is negligible.

Figure 9
Effect of predictive target acceleration on the frequency response of pursuit eye movements.

The fact that the weighting of the predictive acceleration component depended on the target's velocity profile (SS or CS) could potentially account for the differences pursuit velocity illustrated in Figs. 4 and and5.5. We tested this possibility by computing the regression between the differences in the fits of the model to the data (for the two conditions in which the motion along the other coordinate differed) and the experimentally observed differences at intervals in which the latter differed significantly from 0 (see Fig. 5). The correlation was significant (p<0.001), but the regression accounted for only a small amount of the variance (r2 = 0.195) and the slope was considerably less than unity (0.18). Furthermore, while the model predicted low-frequency variations in the differences (as expected from Fig. 9), it failed the account for the higher frequency fluctuations observed in Figs. 4 and and55.

Alternative models of prediction in smooth pursuit

We also tested a different set of models in which the speed (ve) and direction (θe) of pursuit were the controlled variables, rather than the x- and y-velocities. A formulation in terms of speed and direction could potentially account for the lack of independence of the horizontal and vertical components of velocity described in Figs. 4 and and5.5. Furthermore, this formulation was suggested by modeling studies of manual tracking [26] and manual interception [4] where the speed and direction of hand movements appeared to be the controlled parameters.

We again assumed a predictive acceleration component An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e009.jpg achieved by means of low-pass filtering as defined by eq. 3. We then assumed that the rate of change in speed and direction would be proportional to error signals proportional to the difference between eye speed (ve) and predicted target speed and between the direction of pursuit (θe) and the predicted target direction (θt):

equation image
equation image
equation image

In eq. 5, An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e013.jpg is a unit vector in the direction of eye velocity ve at time t and g is a gain factor for target velocity vt. As was the case for the model given by eq. 4, we assumed that the gains in the x- and y-directions could differ. The model also assumed that the weighting coefficients c of the predictive acceleration component for direction and speed could differ. The parameters a (eq. 5) and d (eq. 7) represent the time constants for speed and direction, respectively.

As was the case for the previous model, we fixed the time delay at 80 ms and then identified the 7 parameters (a, b, c1, c2, d, gx, gy) that gave the best fit to all of the data (5 paths×4 speed profiles) for each of the 6 subjects. We then fitted the data separately for each path and speed condition, assuming a fixed value for the three time constants (a, b, d) but letting the gains and the weighting coefficients for the acceleration vary from condition to condition. Even though the model given by eq. 5–7 had one more free parameter than did the previous model (eq. 4), it gave a poorer fit to the data, the VNAF increasing by 8.5%. A paired samples t-test confirmed that the previous model gave a significantly better fit to the data (t119 = 3.673, p<0.001).

We also tested a variant of the model in eq. 5–7 in which the normal acceleration An external file that holds a picture, illustration, etc.
Object name is pone.0012574.e014.jpg was proportional to the directional error, rather than the angular velocity of the eye in eq. 7. However, this variant gave an even poorer fit, the VNAF increasing by an additional 22%. Therefore, we could not find any support for the hypothesis that speed and direction, rather than eye velocity, were the variables that were controlled in smooth pursuit.


We presented subjects with a target that moved in two dimensions along trajectories that varied in their predictability. In the simplest case, we generated the trajectory from a sum of two sinusoids (SS): the fundamental frequency and the second or the third harmonic. With this formulation, the target's velocity, speed and direction changed smoothly throughout the trial and speed and the rate of change of direction were inversely correlated, i.e. they approximately followed the power law relation between speed and curvature [27], [28]. It has been shown that tracking errors are smaller for motions obeying this relation [29] and that the motion appears perceptually to be more uniform [30].

In a second condition, the target followed the same path but with a different time course, namely at a constant speed (CS). However, even though the speed was entirely predictable, target direction could change abruptly, as could the horizontal and vertical components of target velocity (see Fig. 2). We reasoned that, due to such abrupt changes, the target motion would be less predictable. In fact, the speed of pursuit in this condition was not constant but it was modulated in a consistent fashion (Fig. 2B). We also examined two more conditions, in which we combined smooth target motion (SS) along one axis with the target's velocity for the CS condition along the other axis. This set of stimuli also permitted us to assess the extent to which the horizontal and vertical components of smooth pursuit were independent of each other [14].

We tested several different quantitative models with the aim of identifying the signals that governed smooth pursuit under steady state conditions, and the influence of the motion's predictability on the relative contribution of these signals. We began with a simplified version of the model introduced by Krauzlis and Lisberger [16], namely one in which eye acceleration is proportional to retinal image velocity. This model gave a reasonable fit to the data, with a time constant (160 ms) that is close to the time constant that can be estimated from the response to a step in target velocity [16], [18], [25]. However, in this model we permitted the time delays for target and eye velocity to differ, and the time delay for target velocity that gave the best fit (20 ms) was much less than the time delay for eye velocity (120 ms) and also much less than the latency for the initiation of smooth pursuit (~100 ms). This observation suggested that a signal predicting future target motion was added to the retinal image motion signal. We modeled this as low-pass filtered version of target acceleration, reasoning that high frequency changes in acceleration would be essentially unpredictable. This model gave a reasonable fit to the data, and the relative importance of this predictive signal was greater for the more predictable set of target motions.

The model we have proposed is physiologically plausible. Target motion, rather than retinal image motion is encoded in the activity of MST neurons [31][34]. An acceleration signal is represented, albeit weakly in the activity of MT neurons [35], [36]. Whether or not this parameter is encoded by activity of neurons in MST is not known. However, a predictive acceleration signal, as contemplated by our model, could also be derived from a velocity signal by intrinsic mechanisms such as short-term synaptic depression and spike-frequency adaptation [37].

In our model, we did not include a term proportional to image acceleration, as was done by Krauzlis and Lisberger [16]. However, in the range of target accelerations experienced in the present experiment, they found the gain for acceleration to be <10% of the velocity contribution and thus smaller than the predictive acceleration term in our model. Recently, two groups [38], [39] have probed the dynamics of smooth pursuit using random perturbations. Both groups estimated the time constants of the response to such perturbations in the range from 40 to 60 ms, i.e. much faster than the time constants estimated in our experimental conditions. However, the gain of the response was found to decrease rapidly with the amplitude of the stimulus, decreasing to a value of about 0.2 when the standard deviation of the noise stimulus was 8°/s. Conceivably, the addition of a low-gain, faster response to our model could have improved the fit of the model to our data, especially in the high-frequency domain.

We assumed that the relative contribution of the predictive acceleration component was constant throughout one cycle of target motion, after steady-state conditions had been achieved (i.e. 1 s after motion onset). Prior to this time, the contribution of retinal image acceleration or target acceleration to smooth pursuit is modest, at best [16], [25]. Target acceleration also generally is not perceived directly [40] nor does it contribute to direct the arm during manual interception tasks [4] for targets following trajectories similar to the ones studies here. Under other conditions, such as the interception of objects accelerated by gravity, subjects can accurately time an interception movement [41], but this behavior is highly context dependent [5] and appears to reflect the implementation of an internal model of motion in a gravitational field rather than the use of a predictive extrapolation as described by our model. (Note that in most manual interception tasks, target motion is viewed for less than one second).

Thus the gain of the predictive acceleration component in our model should be expected to increase gradually over the course of the first second or so. In fact, in a situation where the target was subsequently occluded, Bennett et al. [42] have reported that only when an accelerating target was in view for more than 500 ms was the acceleration reflected in pursuit eye movements. In our model, we assumed that the gain was constant after the first second. This assumption is probably incorrect and the gain of this component most likely fluctuates over time. For example, a Kalman filter would generate a predictive component whose importance could fluctuate [43][45]. However, since the experiments were not specifically designed to test for this possibility, we did not incorporate it into our model. Finally, we examined only linear models and some of the discrepancies between the model's behavior and the experimental data could have arisen from amplitude-dependent nonlinearities in the response.

One aim of our experiments was to assess the extent to which horizontal and vertical components of pursuit are independent of each other. We found that, for the most part they are, but we also found instances in which horizontal eye velocity was affected by the vertical component of the target signal and vice versa (see Figs. 3 and and4).4). The effect was generally small (<10%) and infrequent (<10% of the time) and was more pronounced for the higher frequency components of the response. The interdependence of the horizontal and vertical components of pursuit eye velocity could be accounted for in part by allowing the gain of the predictive acceleration component to depend on the overall target motion. However, this factor accounted for only a small percentage of the variance and it accounted mostly for the low frequency components of the difference, i.e. the region where the predictive feedback component is most effective (see Fig. 9). However, we are unable to account for the higher frequency components of this phenomenon.


We thank Dr. Martha Flanders for helpful suggestions during the course of this work.


Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by National Institutes of Health Grant NS-15018 ( The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Mrotek LA, Soechting JF. Target interception: hand-eye coordination and strategies. J Neurosci. 2007;27:7297–7309. [PubMed]
2. Port NL, Lee D, Dassonville P, Georgopoulos AP. Manual interception of moving targets. I. Performance and movement initiation. Exp Brain Res. 1997;116:406–420. [PubMed]
3. Soechting JF, Flanders M. Extrapolation of visual motion for manual interception. J Neurophysiol. 2008;99:2956–2967. [PubMed]
4. Soechting JF, Juveli JZ, Rao HM. Models for the extrapolation of target motion for manual interception. J Neurophysiol. 2009;102:1491–1502. [PubMed]
5. Zago M, McIntyre J, Senot P, Lacquaniti F. Visuo-motor coordination and internal models for object interception. Exp Brain Res. 2009;192:571–604. [PubMed]
6. Kowler E, Martins AJ, Pavel M. The effect of expectations on slow oculomotor control. IV. Anticipatory smooth eye movements depend on prior target motions. Vision Res. 1984;24:197–210. [PubMed]
7. Kowler E. Cognitive expectations, not habits, control anticipatory smooth oculomotor pursuit. Vision Res. 1989;29:1049–1057. [PubMed]
8. Lisberger S, Morris EJ, Tychsen L. Visual motion processing and sensory-motor integration for smooth pursuit eye movements. Ann Rev Neurosci. 1987;10:97–129. [PubMed]
9. Becker W, Fuchs AF. Prediction in the oculomotor system: smooth pursuit during transient disappearance of a visual target. Exp Brain Res. 1985;57:562–575. [PubMed]
10. Bennett SJ, Barnes GR. Human ocular pursuit during the transient disappearance of a visual target. J Neurophysiol. 2003;90:2504–2520. [PubMed]
11. Bennett SJ, Barnes GR. Predictive smooth ocular pursuit during the transient disappearance of a visual target. J Neurophysiol. 2004;92:578–590. [PubMed]
12. Collewijn H, Tamminga EP. Human smooth and saccadic eye movements during voluntary pursuit of different target motions on different backgrounds. J Physiol. 1984;351:217–250. [PubMed]
13. Deno DC, Crandall WF, Sherman K, Keller EL. Characterization of prediction in the primate visual smooth pursuit system. BioSystems. 1995;34:107–128. [PubMed]
14. Kettner RE, Leung HC, Peterson BW. Predictive smooth pursuit of complex two-dimensional trajectories in monkey: component interactions. Exp Brain Res. 1996;108:221–235. [PubMed]
15. Yasui S, Young LR. On the predictive control of foveal eye tracking and slow phases of optokinetic and vestibular nystagmus. J Physiol. 1984;347:17–33. [PubMed]
16. Krauzlis RJ, Lisberger SG. A model of visually-guided smooth pursuit eye movements based on behavioral observations. J Comput Neurosci. 1994;1:265–283. [PubMed]
17. Ringach DL. A tachometer feedback model of smooth pursuit eye movements. Biol Cybern. 1995;73:561–591. [PubMed]
18. Robinson DA, Gordon JL, Gordon SE. A model of the smooth pursuit eye movement system. Biol Cybern. 1986;55:43–57. [PubMed]
19. Churchland MM, Lisberger SG. Experimental and computational analysis of monkey smooth pursuit eye movements. J Neurophysiol. 2001;86:741–759. [PubMed]
20. Batschelet E. Circular statistics in biology. London: Academic Press; 1981.
21. Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in C. 1992. Cambridge University Press, New York, NY.
22. Nelder JA, Mead R. A simplex method for function minimization. Comp J. 1964;7:308–313.
23. Barnes GR, Donnelly SF, Eason RD. Predictive velocity estimation in the pursuit reflex response to pseudo-random and step displacement stimuli in man. J Physiol. 1987;389:111–136. [PubMed]
24. Mrotek LA, Flanders M, Soechting JF. Oculomotor responses to gradual changes in target direction. Exp Brain Res. 2006;172:175–192. [PubMed]
25. Soechting JF, Mrotek LA, Flanders M. Smooth pursuit tracking of an abrupt change in target direction: vector superposition of discrete responses. Exp Brain Res. 2005;160:245–258. [PubMed]
26. Engel KC, Soechting JF. Manual tracking in two dimensions. J Neurophysiol. 2000;83:3483–3496. [PubMed]
27. Lacquaniti F, Terzuolo C, Viviani P. The law relating the kinematic and figural aspects of drawing movements. Acta psychol. 1983;54:115–130. [PubMed]
28. Viviani P, Cenzato M. Segmentation and coupling in complex movements. J Exp Psychol Hum Perc Perf. 1985;11:828–845. [PubMed]
29. de'Sperati C, Viviani P. The relationship between curvature and velocity in two-dimensional smooth pursuit eye movements. J Neurosci. 1997;17:3932–3945. [PubMed]
30. Viviani P, Stucchi N. Biological movements look uniform: evidence of motor-perceptual interactions. J Exp Psychol Hum Percept Perform. 1992;18:603–623. [PubMed]
31. Chukoskie L, Movshon JA. Modulation of visual signals in macaque MT and MST neurons during pursuit eye movement. J Neurophysiol. 2009;102:3225–3233. [PubMed]
32. Duffy CH, Wurtz RH. Sensitivity of MST neurons to optic flow stimuli. I. A continuum of response selectivity to large field stimuli. J Neurophysiol. 1991;65:1329–1345. [PubMed]
33. Inaba N, Shinomoto S, Yamane S, Takemura A, Kawano K. MST neurons code for visual motion in space independent of pursuit eye movements. J Neurophysiol. 2007;97:3473–3483. [PubMed]
34. Newsome WT, Wurtz RH, Komatsu H. Relation of cortical areas MT and MST to pursuit eye movements. II. Differentiation of retinal from extraretinal inputs. J Neurophysiol. 1988;60:604–620. [PubMed]
35. Lisberger SG, Movshon JA. Visual motion analysis for pursuit eye movements in area MT of macaque monkeys. J Neurosci. 1999;19:2224–2246. [PubMed]
36. Schlack A, Krekelberg B, Albright TD. Recent history of stimulus speeds affects the speed tuning of neurons in area MT. J Neurosci. 2007;27:11009–11018. [PubMed]
37. Puccini GD, Sanchez-Vives MV, Compte A. Integrated mechanisms of anticipation and rate-of-change computations in cortical circuits. PLoS Comput Biol. 2007;3(5):e82. [PubMed]
38. Osborne LC, Lisberger SG. Spatial and temporal integration of visual motion signals for smooth pursuit eye movements in monkeys. J Neurophysiol. 2009;102:2013–2025. [PubMed]
39. Tavassoli A, Ringach DL. Dynamics of smooth pursuit maintenance. J Neurophysiol. 2009;102:10–118. [PubMed]
40. Brouwer AM, Brenner E, Smeets JB. Perception of acceleration with short presentation times: can acceleration be used in interception? Percept Psychophys. 2002;64:1160–1168. [PubMed]
41. Zago M, Bosco G, Maffei V, Iosa M, Ivanenko YP, Lacquaniti F. Internal models of target motion: expected dynamics overrides measured kinematics in timing manual interceptions. J Neurophysiol. 2004;91:1620–1634. [PubMed]
42. Bennett SJ, Orban de Xivry JJ, Barnes GR, Lefèvre P. Target acceleration can be extracted and represented within the predictive drive to ocular pursuit. J Neurophysiol. 2007;98:1405–1414. [PubMed]
43. Kramer KA, Stubberud SC. An analysis and implementation of a neural extended Kalman filter for target tracking. Int J Neur Sys. 2006;16:1–13. [PubMed]
44. Welch G, Bishop G. An introduction to the Kalman filter. 2003. (Tech. Rep. 95-041) Chapel Hill: Univ North Carolina. (
45. Wu W, Gao Y, Bienenstock E, Donoghue JP, Black MJ. Bayesian population decoding of motor cortical activity using a Kalman filter. Neur Comp. 2006;18:80–118. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science