A Test of Torque-Control and Equilibrium-Point Models of Motor Control
ABSTRACT
Based upon predictions of equilibrium point models of motor control, a recent article by Jaric et al. (1999), reported that the peak velocity of an inertially loaded single degree-of-freedom elbow movement was not altered to a statistically distinguishable degree by the accuracy of the subjects expectation about the size of the load. They concluded that their results confirmed their equilibrium point models and were incompatible with what they described as torque control models. A review of the literature and analysis of a more extensive set of data does not support that conclusion. To the contrary, if the actual load is heavier (or lighter) than expected, it is moved more slowly (or faster) than it would have been, had expectation matched reality. Torque-control models predict those consequences, equilibrium point models do not.
INTRODUCTION
A recent paper by Jaric et al. (1999) describes a series of experiments on elbow flexion movements in which inertial loads are intermittently changed without the subjects' knowledge. It is generally found that movements are slower with heavier than with lighter loads. At issue is whether subject's expectation about the size of the load is also a factor. They analyze the movements in terms of two versions of the equilibrium point (EP) hypotheses (Latash 1993; Feldman and Levin 1995). According to Jaric et al. (1999), whether a subject expects a light load or a heavy load, "movement velocity is supposed to reach approximately the same level in both cases" (page 53) in Feldman's version of the EP hypothesis. In Latash's version "peak velocity may be expected to drop slightly [if the subject encounters a heavier than expected load] the difference may be expected to be about 10%" (page 55). They contrast this prediction with what they term "torque-control hypotheses" (Atkeson 1989; Bock 1990; Gottlieb 1996) that they suggest predict substantial (as much as 30-50%) decreases in peak velocity with loads that are heavier than expected. Jaric et al. (1999) find, for six subjects, no statistically significant effects of expectation on peak velocity. They conclude (page 64) that "Gottlieb's [torque-control] hypothesis is wrong in 2 of four predictions [while] Both versions of the EP-hypothesis make qualitatively similar predictions which all happen to be correct."
Given such clearly stated differences, it seems appropriate to examine our own model of motor control as well as those that seem to be confirmed by Jaric et al. (1999). A clear, unambiguously negative finding for any theory's predictions provides a strong rational basis for rejecting it. We will only examine the effects of expectation on peak velocity because it is a variable that can be measured unambiguously and is also the variable of clearest disagreement.
METHODS
The experiments performed for this study are a subset of some that are described in greater detail in (Gottlieb 1996). Neither the data nor its mode of analysis reported here were part of the earlier paper.
Subjects made 54° elbow flexions, as "fast and accurately" as possible, using a manipulandum that allowed rotation about the elbow in the horizontal plane. The right upper arm was abducted 90° and the elbow started from 60° of flexion from a straight arm. A digital computer sampled joint angle and tangential acceleration at 1000/s using a 12-bit A/D converter. A torque motor that could vary its torque using acceleration feedback controlled the effective inertia of the manipulandum. All subjects gave informed consent according to Rush Medical Center procedures before participating.
Subjects made movements under two conditions. Under known conditions they practiced with a fixed inertial load and then 11-15 movements were recorded. Under uncertain conditions, they made 30-60 movements. For most movements, the actual moment of inertia of a typical arm plus manipulandum ( about 0.18 kg m2) was the total load. Infrequently (1/3 to 1/6 of the time) the torque motor added an additional load. It could take 1 to 3 different values over a range of 0.13 kg m2 to +0.9 kg m2. Loads were intermixed in a pseudo-randomized order. We assumed that under uncertain conditions, the subjects expected the most likely load and planned accordingly. Hence, they experienced actual loads that were either equal to expected (A=X), heavier than expected (A>X) or lighter than expected (A<X). In our analysis, we compare peak movement velocity under known conditions with the velocity under uncertain conditions for the same load.
Peak Velocity under uncertain load conditions is plotted against Peak Velocity under known load conditions. + symbols : The actual and expected load are the same. Solid circles: The actual load is greater than the expected load. Solid boxes: The actual load is less than the expected load. Individual values represent the average of 11-40 movements by a subject with a given load and expectation. The dashed lines have slopes or 1.1, 1.0 and 0.9. The open symbols show the data from Table 2 of
Jaric et al. (1999) with each point representing the average of their 6 subjects. Solid lines are linear regressions through the origin.RESULTS
Figure 1 shows peak movement velocities of 7 subjects. Each point represents the velocity for a particular load, moved under known conditions (on the abscissa) and the velocity with the same load, moved under uncertain conditions (on the ordinate). Different symbols distinguish the three groups described above. A dashed line of unity slope through the origin defines the null-hypothesis. Two dotted lines bound ±10% changes from unity. Solid lines show linear regressions through the origin for each of the three conditions. The null hypothesis is that expectation of load size has no effect on movement velocity. According to Jaric et al. (1999), Feldman's version of the EP-hypothesis predicts that the points from all three groups should lie along the line unity slope while Latash's predicts that they lie within the ±10% lines. The Table presents the results of analysis comparing each uncertain group to the known group using Student's t-test.
|
Hypothesis |
DOF |
P |
|
PVA=X=PVKnown |
11 |
.90 |
|
PVA>X=PVKnown |
9 |
.0018 |
|
PVA<X=PVKnown |
11 |
<.0001 |
TABLE 1: A test of the hypothesis that Peak Velocity (PV) under uncertain load conditions is the same under known load conditions using Student's t-test. In row 1 the actual load was equal to the expected (i.e. the most likely) load. In row 2 the actual load was greater than expected. In row 3 the actual load was less than expected.
These results show a three-fold variation in peak velocity with load and across subjects. Subjects moved at the same velocity under known conditions as they did under uncertain conditions if the actual load equaled the expected one (A=X). The slope of linear regression through the origin is 1.003 (1.041-0.965 are the 95% confidence limits), r2=0.97. Similar velocities were not found if the actual load was greater or less than expected. For any actual load, subjects moved significantly more slowly if it was heavier than expected and significantly faster if it was lighter. The slope of linear regression for A<X is 1.171 (1.262-1.080), r2=0.83 and for A>X is 0.807 (0.859-0.766) r2=0.81.
DISCUSSION
Under known load conditions, it seems reasonable to assume that subjects expect and can plan for the actual load. Under uncertain load conditions, what do subjects do? Jaric et al. (1999) assume that they plan for the most recent load. This seems reasonable for the load sequences they used. In our more randomized presentation we assume that subjects plan for the most likely load. The reasonableness of this assumption is supported by the fact that peak velocities under known conditions are indistinguishable from those seen if the most likely load was encountered under uncertain conditions. An alternative might be for subjects to assume a weighted, probabilistic estimate based on the frequency of the different loads. This would seem a better scheme for experiments in which two different loads are equally likely (Latash 1994). In our experiment it would produce only slightly different expectations because our load changes were presented infrequently.
Under most conditions, the fact that subjects move faster with lighter loads and slower with heavier loads is not obligatory. It is a consequence of how the movement is planned, based on experience and incorporated into some form of internal model (Bhushan and Shadmehr 1999). If subjects plan for the most likely load, what is it that they plan? In (Gottlieb et al. 1989) we postulated that "Strategies [for the control of voluntary movement] consist of sets of rules that determine the patterns of muscle activation [which] lead to patterns of muscle torques and EMGs." This may be informally called a "torque-control" hypothesis.
[Even before Hill (1938) it was widely understood that the external force generated by a muscle depends upon its length and rate of shortening as well as its level of excitation. Those compliant properties are ever present and important to any control scheme (Hogan 1984; Hasan 1986; Zajac 1989). In fact, "control of torque" usually comes down to control of muscle excitation patterns (Gottlieb et al. 1989; Hoffman and Strick 1993)). Therefore, "compliant torque-control" is a better term but the simpler one is commonly used and is acceptable as long as one does not take it literally to mean that the CNS can control the instantaneous value of muscle torque regardless of the load.] In such a scheme, a significant contribution from reflexes only occurs if the planned torque is not well matched to the load and as a result, the actual trajectory deviates from the expected trajectory (Smeets et al. 1990; Gottlieb 1996). Because muscle activation patterns are planned, in part, on the expected size of the load, false expectations will change muscle forces and thereby movement velocities. EP-hypothesis planning is done in terms of an equilibrium trajectory, albeit one that varies in shape according to the version of the hypothesis. Jaric et al. (1999) develops the basis for their EP predictions in considerable detail that need not be recapitulated here. The Feldman/Latash EP hypotheses are based upon the principal that reflex pathways make significant contributions to muscle activation patterns under all movement conditions because it is the deviation between the actual and the equilibrium trajectory that drives motor neuron activation.[EP hypotheses actually require at least three commands . Jaric et al. discuss the equilibrium trajectory command (sometimes called a reciprocal command, R) but have little to say about the coactivation command (C) or the damping coefficient m . These two commands act as "gain controls" on the effects of R and on reflex inputs.].Jaric et al. (1999) described the predictions made by EP and torque control approaches and summarized them in their Table 1. According to the EP-hypotheses, unexpected changes in inertial load should produce peak velocities that differ from those with equal but known loads by only 0-10%. They predict those differences to be as much as 30-50% according to torque-control hypotheses. We find that average differences are +17±9% and -19±5%. The extrema of the data are +43% and -38%. With expected loads, 12 of 12 data points lie within the ±10% deviation lines. With unexpected loads, 17 of the 21 points lie outside those lines.
Furthermore, figure 2 of (Latash 1994), figure 2 of (Gottlieb 1994) and figures 4,5 and 10A of (Gottlieb 1996) all show movements with the same load and different expectations that visibly differ in peak velocity. In all cases, if X>A, peak velocity is larger and if X<A, it is smaller. The data in table 2 of Jaric et al. (1999) also show this directional consistency and the six open symbols in figure 1 plot those data [They did not conduct experimens under known conditions, so the (A=X) condition is used on the abscissa.]. Although their results did not reach statistical significance, all data lie on the sides of the unity-slope line that are predicted by torque-control hypotheses. These results are consistent with the predictions of the torque-control hypothesis and with neither version of the EP hypothesis.
Immediately following this paper in the same issue is a reply by Jaric & Latash. They were sent a copy of my manuscript by the editors. To the best of my knowledge, their reply was accepted without external review. My request to publish a further reply was refused. That rejected reply, in extended form, is found here.
REFERENCES CITED
Atkeson, C.G., 1989 Learning arm kinematics and dynamics. Annual Review of Neuroscience 12: 157-183
Bhushan, N, Shadmehr, R,. 1999 Computational nature of human adaptive control during learning of reaching movements in force fields. Biol Cybern 81: 39-60
Bock, O,. 1990 Load compensation in human goal-directed arm movements. Behavioral Brain Research 41: 167-177
Feldman, A.G., Levin, M.F., 1995 The origin and use of positional frames of
reference in motor control. Behavioral and Brain Sciences 18: 723-806Gottlieb, G.L., 1994 The generation of the efferent command and the importance of joint compliance in fast elbow movements. Experimental Brain Research 97: 545-550
Gottlieb,G.L., 1996 On the voluntary movement of compliant inertial-viscoelastic loads by parcellated control mechanisms. Journal of Neurophysiology 76: 3207-3229
Gottlieb, G.L., Corcos DM, Agarwal GC 1989 Organizing principles for single joint movements: I - A speed-insensitive strategy. Journal of Neurophysiology 62: 342-357
Hasan, Z., 1986 Optimized movement trajectories and joint stiffness in unperturbed, inertially loaded movements. Biol. Cybernetics 53: 373-382
Hill, A.V., 1938 The heat of shortening and the dynamic constants of muscle. Proceedings of the Royal Society London B126: 136-195
Hoffman, D.S., Strick PL 1993 Step-tracking movements of the wrist: III. Influence of changes in load on patterns of muscle activity. Journal of Neuroscience 13: 5212-5227
Hogan, N., 1984 Adaptive control of mechanical impedance by coactivation of antagonist muscles. IEEE Trans. Auto. Control AC-29: 681-690
Jaric, S, Milanovic, S, Blesic, S, Latash, M.L, 1999 Changes in movement kinematics during single-joint movements against expectedly and unexpectedly changed inertial loads. Human Movement Sciences 18: 49-66
Latash, M.L., 1993 Control of Human Movement. Human Kinetics, Urbana, IL
Latash, M.L., 1994 Control of fast elbow movement: A study of electromyographic patterns during movements against unexpectedly decreased inertial load. Experimental Brain Research 98: 145-152
Smeets, J.B.J., Erkelens, C.J., Denier van der Gon, J.J., 1990 Adjustments of fast goal-directed movements in response to an unexpected inertial load. Experimental Brain Research 81: 303-312
Zajac, F.E., 1989 Muscle and Tendon: Properties, Models, Scaling and Application to Biomechanics and Motor Control. CRC Crit. Rev. in Biomed. Eng. 17: 359-411