JN AJP: Gastrointestinal and Liver Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 88: 1097-1118, 2002;
0022-3077/02 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (118)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Peterka, R. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Peterka, R. J.

The Journal of Neurophysiology Vol. 88 No. 3 September 2002, pp. 1097-1118
Copyright ©2002 by the American Physiological Society

Sensorimotor Integration in Human Postural Control

R. J. Peterka

Neurological Sciences Institute, Oregon Health & Science University, Portland, Oregon 97006


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Peterka, R. J.. Sensorimotor Integration in Human Postural Control. J. Neurophysiol. 88: 1097-1118, 2002. It is generally accepted that human bipedal upright stance is achieved by feedback mechanisms that generate an appropriate corrective torque based on body-sway motion detected primarily by visual, vestibular, and proprioceptive sensory systems. Because orientation information from the various senses is not always available (eyes closed) or accurate (compliant support surface), the postural control system must somehow adjust to maintain stance in a wide variety of environmental conditions. This is the sensorimotor integration problem that we investigated by evoking anterior-posterior (AP) body sway using pseudorandom rotation of the visual surround and/or support surface (amplitudes 0.5-8°) in both normal subjects and subjects with severe bilateral vestibular loss (VL). AP rotation of body center-of-mass (COM) was measured in response to six conditions offering different combinations of available sensory information. Stimulus-response data were analyzed using spectral analysis to compute transfer functions and coherence functions over a frequency range from 0.017 to 2.23 Hz. Stimulus-response data were quite linear for any given condition and amplitude. However, overall behavior in normal subjects was nonlinear because gain decreased and phase functions sometimes changed with increasing stimulus amplitude. "Sensory channel reweighting" could account for this nonlinear behavior with subjects showing increasing reliance on vestibular cues as stimulus amplitudes increased. VL subjects could not perform this reweighting, and their stimulus-response behavior remained quite linear. Transfer function curve fits based on a simple feedback control model provided estimates of postural stiffness, damping, and feedback time delay. There were only small changes in these parameters with increasing visual stimulus amplitude. However, stiffness increased as much as 60% with increasing support surface amplitude. To maintain postural stability and avoid resonant behavior, an increase in stiffness should be accompanied by a corresponding increase in damping. Increased damping was achieved primarily by decreasing the apparent time delay of feedback control rather than by changing the damping coefficient (i.e., corrective torque related to body-sway velocity). In normal subjects, stiffness and damping were highly correlated with body mass and moment of inertia, with stiffness always about 1/3 larger than necessary to resist the destabilizing torque due to gravity. The stiffness parameter in some VL subjects was larger compared with normal subjects, suggesting that they may use increased stiffness to help compensate for their loss. Overall results show that the simple act of standing quietly depends on a remarkably complex sensorimotor control system.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Bipedal upright stance is inherently unstable. A small sway deviation from a perfect upright position results in a torque due to gravity that accelerates the body further away from the upright position. To maintain upright stance, the destabilizing torque due to gravity must be countered by a corrective torque exerted by the feet against the support surface. A widely held view is that the corrective torque is generated through the action of a feedback control system (see reviews by Horak and Macpherson 1996; Johansson and Magnusson 1991). We will refer to this corrective torque, which necessarily involves a time delay due to sensory transduction, transmission, processing, and muscle activation, as "active" torque. However, controversy remains (Morasso and Schieppati 1999) because another view holds that the corrective torque is generated by muscle "tone" that acts without time delay (Winter et al. 1998, 2001). In this paper, we will refer to corrective torque that acts without time delay as "passive" torque. Finally, a third view recognizes that feedback mechanisms contribute to postural stabilization but states that feedback alone is insufficient and that feedforward predictive mechanisms are required to explain postural control behavior (Fitzpatrick et al. 1996). Our results support the view that active torque generated by feedback control mechanisms is the dominant contributor to quiet stance control.

Visual, proprioceptive, and vestibular systems clearly contribute to postural control because numerous studies have shown that stimulation of visual (Berthoz et al. 1979; Bronstein 1986; Dijkstra et al. 1994a; Lee and Lishman 1975; Lestienne et al. 1977; van Asten et al. 1988a), proprioceptive (Allum 1983; Jeka et al. 1997; Johansson et al. 1988; Kavounoudias et al. 1999), or vestibular systems (Day et al. 1997; Hlavacka and Nijiokiktjien 1985; Johansson et al. 1995; Nashner and Wolfson 1974) evoke body sway. However, little is known about how information from these senses is processed and combined to generate appropriate corrective torque when there is conflicting or inaccurate orientation information from different sensory systems. One possibility is that sensory cues are combined in an essentially linear manner. That is, each sensory system detects an "error" indicating deviation of body orientation from some reference position. Vestibular sensory cues detect deviations of head orientation from earth-vertical (gravity), visual sensors detect head orientation relative to the visual world, and proprioceptors detect leg orientation relative to the support surface. The individual error signals are summed, and appropriate corrective torque is generated as a function of this summed signal. Note that in this paper, for modeling purposes, we use a restricted definition of proprioceptive cues as only those sensory cues signaling body motion relative to the support surface. Additionally, we assume that appropriate neural transformations are performed on the various sensory cues so that the nervous system has information on body center-of-mass (COM) motion relative to each sensory reference (i.e., the direction of gravity for vestibular cues, visual world orientation for visual cues, and support surface orientation for proprioceptive cues). Psychophysical studies support the fact that such transformations can occur (Mergner et al. 1991, 1997).

Previous experimental results, where body sway was evoked by manipulation of individual and combined sensory cues, appear to be consistent with an essentially linear model (Fitzpatrick et al. 1996; Hajos and Kirchner 1984; Jeka et al. 1998, 2000; Johansson et al. 1988; Maki et al. 1987; Schöner 1991; van Asten et al. 1988b). Many of these earlier studies developed linear models that assumed that the postural control system was inherently stable, with experimental stimuli merely perturbing this inherently stable system. A more complete understanding of postural control must explain how this apparent inherent stability is actually achieved.

Most studies of human postural control have employed transient stimuli (e.g., sudden support surface motions) to evoke characteristic postural responses (Allum 1983; Diener et al. 1984b; Horak and Nashner 1986; Nashner 1977) or methods that artificially stimulate individual sensory receptors [i.e., muscle or tendon vibration (Kavounoudias et al. 1999) and galvanic vestibular stimulation (Nashner and Wolfson 1974; Watson and Colebatch 1997)]. We chose to investigate postural control using motion stimuli (tilts of the support surface and/or visual surround) that continuously perturb the system. Continuously varying stimuli (often sinusoidal stimuli over a range of frequencies, or more complex random or pseudorandom time series) evoke responses that eventually achieve a steady state. Steady-state stimulus-response data can be used to obtain transfer functions that characterize the dynamic properties of the system (Bendat and Piersol 2000). These techniques have been used previously to investigate postural control in humans and animals (Dijkstra et al. 1994b; Fitzpatrick et al. 1996; Hajos and Kirchner 1984; Ishida and Imai 1983; Jeka et al. 1998, 2000; Johansson et al. 1988; Maki et al. 1987; Peterka and Benolken 1995; Talbott 1980) but have not been systematically applied to investigate dynamic behavior over a wide range of conditions. The use of continuously applied perturbations seems appropriate to study quiet stance behavior, which itself is a continuously active process. This is in contrast to transient stimuli that may trigger specific and equally transient motor programs, which may not be directly related to the continuous regulation of balance.

Our results show that sensory integration and postural regulation do appear to be essentially linear processes for a specific sensory condition and a given stimulus amplitude. However, as stimulus conditions change, nonlinearities become apparent. The major nonlinearity occurs with changing stimulus amplitudes where there is an apparent graded shift in the source of sensory information contributing to postural control, with increasing utilization of vestibular cues as visual and proprioceptive perturbations increase. In subjects with absent vestibular function, this shift cannot occur, and their overall behavior remains quite linear independent of stimulus amplitude.

Such context-dependent changes in sensory utilization are in general agreement with previous views of postural behavior (Forssberg and Nashner 1982; Horak and Macpherson 1996; Nashner et al. 1982), experimental findings using galvanic vestibular stimulation (Britton et al. 1993; Fitzpatrick et al. 1994), and motor control in general (Hultborn 2001; Prochazka 1989). Our results provide quantitative measures of the stimulus-dependent changes in sensory contributions to postural control. In addition, our results provide estimates of important postural control parameters (stiffness, damping, time delay) and demonstrate how these parameters change in different sensory environments and stimulus conditions.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Subjects

The experimental protocols were approved by the Institutional Review Board at Oregon Health & Science University and were performed in accordance with the 1964 Helsinki Declaration. Prior to testing, all subjects gave their informed consent. Twelve subjects participated in this study. Eight were adults who had normal results on clinical sensory organization tests of postural control (Peterka and Black 1990) and had no known history of balance impairment or dizziness. The other four subjects had profound bilateral vestibular loss (VL subjects) as confirmed by clinical rotation testing and results from sensory organization tests of postural control (Nashner 1993a). The causes and durations of vestibular loss in these subjects are given in Table 1. Table 1 also shows the gain of horizontal vestibuloocular reflex (HVOR) for 0.05- and 0.2-Hz yaw rotations for the VL subjects. The HVOR gain for all VL subjects was well below the 95th percentile for the normal population (Peterka et al. 1990). The age range of normal subjects was 24-46 yr, while the VL subjects ranged in age from 45 to 58 yr. Although there was an age difference between these two groups, previous research has identified only minor changes in postural control in subjects in these age ranges (Peterka and Black 1990).


                              
View this table:
[in this window]
[in a new window]
 
Table 1. Vestibular loss subject demographics

Experimental setup

All experiments were performed on a custom balance-testing device that included a motor-driven support surface and visual surround (Fig. 1). Position servo-controlled motors produced anterior/posterior (AP) tilts of the support surface and visual surround with the rotation axes collinear with the subject's ankle joints. Vertical force sensors in the support surface were used to measure center-of-pressure (COP) data. The visual surround had a half-cylinder shape (70-cm radius) and was lined with a complex checkerboard pattern consisting of white, black, and three gray levels (see Fig. 1). During testing, the room lights were off, and the visual surround was illuminated by fluorescent lights attached to the right and left edges of the surround.



View larger version (121K):
[in this window]
[in a new window]
 
Fig. 1. Balance test device. The subject stands on a support surface and views a high contrast visual surround. Both the support surface and visual surround can rotate in the anterior/posterior (AP) direction about the ankle joint axis. Single-link inverted pendulum dynamics are ensured by use of a backboard assembly. During freestanding trials, lightweight rods attached to potentiometers, located to the right and behind the subject, are used to measure AP body sway at hip and shoulder levels.

During most experiments, the backboard assembly shown in Fig. 1 was used to constrain body motion to that of a single-link inverted pendulum with rotational motion occurring only in the AP direction. The subjects were secured to the backboard with a padded head rest and straps at the level of the subject's knees, hips, and shoulders. The backboard was attached at its base to a pair of bearings aligned with the subject's ankle joint axis. Therefore the subjects were not required to support the full weight of the backboard. However, the influence of backboard mass (12.6 kg) and moment of inertia (11.4 kg m2) was accounted for in the data analysis. The body and backboard together acted as a single-link inverted pendulum. Body/backboard AP COM sway motion was provided by measures of backboard angular position, using a potentiometer, and angular velocity, using a rate sensor (Watson Industries, Eau Claire, WI). We considered these measures to be the output variables of interest.

During some experiments, the backboard assembly was not used and subjects were freestanding. In these experiments, AP body motion was measured by two horizontal rods attached to the subject at hip and shoulder levels (seen in Fig. 1, to the right of the subject). Rotational motion of the sway rods was recorded by potentiometers. Appropriate trigonometric conversions were made later to determine AP body displacement at hip and shoulder levels. A 120-s calibration trial was performed where subjects slowly leaned forward and backward using different combinations of leg and trunk rotations and minimizing knee flexion. A least-squared error curve fit of the following equation was used to obtain estimates of the coefficients ah, as, and b
<IT>x</IT><SUB><IT>cop</IT></SUB>(<IT>t</IT>)<IT>=</IT><IT>b</IT><IT>+</IT><IT>a</IT><SUB><IT>h</IT></SUB><IT>x</IT><SUB><IT>h</IT></SUB>(<IT>t</IT>)<IT>+</IT><IT>a</IT><SUB><IT>s</IT></SUB><IT>x</IT><SUB><IT>s</IT></SUB>(<IT>t</IT>) (1)
where xcop is AP COP displacement, xh is AP body displacement at hip level, xs is AP body displacement at shoulder level, and t is time. Because body movements were very slow, xcop is essentially equal to AP COM displacement (except for small rapid oscillations about the local COM position indicative of AP body acceleration) (Brenière 1996; Winter et al. 1998). In subsequent trials, Eq. 1 was used to calculate AP COM displacement from measures of xh(t) and xs(t). Then an estimate of the subject's COM height (based on anthropometric measures) above the ankle joint was used to calculate the COM rotation angle, which we considered to be the final output variable of interest.

Stimulus delivery and data sampling were computer controlled at a rate of 100/s. Sampled data included: visual surround and support surface angular position, four vertical forces from sensors at the corners of the support surface, rotational position of the hip and shoulder sway rods, and rotational position and velocity of the backboard assembly.

Pseudorandom stimuli

Rotational motion stimuli were based on a pseudorandom ternary sequence (PRTS) of numbers (Davies 1970). The method used to generate the pseudorandom stimulus waveforms is shown in Fig. 2. A stimulus was created from a 242-length PRTS sequence by assigning a rotational velocity waveform a fixed value of +v, 0, or -v°/s according to the PRTS sequence for a duration of Delta t = 0.25 s (Fig. 2B, top). The duration of each stimulus cycle was 60.5 s. The mathematical integration of this PRTS velocity waveform gave a position waveform (Fig. 2B, bottom) that was delivered to the position servo motors to drive visual surround and/or support surface rotation. The PRTS stimulus has a wide spectral bandwidth (Fig. 2C) with the velocity waveform having spectral and statistical properties approximating a white noise stimulus (Davies 1970). As such, this stimulus appeared to be unpredictable to the test subject and thus likely limited any predictive contributions to postural responses that are known to occur (McIlroy and Maki 1994).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 2. Generation of a pseudorandom ternary sequence (PRTS) stimulus. A: a shift register with feedback is used to generate the PRTS. At each time increment, Delta t, the value of each register is shifted to the right. A new value is entered into the left register based on modulo-3 addition of values in registers 3-5, and this new value is taken as the output. With the initial values shown in the 5 shift registers, the first 12 output values of the 242 value periodic sequence are shown below the shift register. This sequence is transformed into a velocity command sequence with values of +v, -v, or 0°/s. B: the velocity sequence is transformed into a velocity waveform by holding each velocity command for Delta t = 0.25 s. The position waveform is the time integral of the velocity waveform. With Delta t = 0.25 s, the PRTS waveform has a period of 60.5 s. C: the spectral composition of the complete velocity and position waveforms for a 1° peak-to-peak PRTS stimulus. Only odd-numbered spectral components have nonzero energy.

Protocol

The PRTS stimulus position waveform was scaled to provide five different stimulus amplitudes (0.5, 1, 2, 4, and 8° peak-to-peak). On most trials, six complete stimulus cycles were presented. On all 0.5° trials, and some 1° trials, eight cycles of the stimulus were presented. Table 2 provides a summary of all six test conditions. Appropriate control trials with duration equivalent to a six-cycle stimulus were also given. The starting value of the PRTS stimulus was selected so that about 80% of the rotational tilt occurred in the positive direction (i.e., toe down support surface tilt and visual surround tilt directed forward and away from the subject) because subjects are able to tolerate larger angle forward body tilts.


                              
View this table:
[in this window]
[in a new window]
 
Table 2. Description of test conditions and performance of vestibular loss subjects

Two of the test conditions used "sway-referencing" of the support surface or visual surround to manipulate the orientation cues available for postural control (Nashner and Berthoz 1978). Sway-referencing was performed by commanding the support surface or visual surround servo systems to continuously track the subject's AP body-sway angle. When the backboard assembly was used, sway-referenced rotations of the visual surround or support surface tracked the rotational movement of the backboard. For freestanding trials, sway-referenced motions of the support surface or the visual surround were in direct proportion to body-sway angles determined from sway-rod measures at the level of the hip or the shoulder, respectively. Sway-referencing alters the normal relationship between body sway and proprioceptive cues (during support surface sway-referencing) or visual cues (during visual surround sway-referencing) and presumably greatly reduces the contribution of these sensory orientation cues.

For normal subjects, all five stimulus amplitudes were presented for each of the six test conditions, and the stimulus and control trials were presented in a randomized order. Trials with sway-referencing included a 10-s delay prior to stimulus onset during which sway-referencing was active. All eight normal subjects completed a full set of trials with the backboard assembly. Four of these subjects also completed the full set of trials without the backboard assembly (free-standing). Data collection for backboard and freestanding conditions was completed in five testing sessions. Each session lasted about 2 h with one session performed per day. A 5-min break was given between trials.

For VL subjects, testing time was limited. Only backboard trials were performed, and the duration of control trials was reduced (to 130 s) compared with those used for normal subjects. Results from the first VL subject showed that 8° amplitude trials were extremely difficult to complete without falling. Therefore all other VL subjects were presented with 1, 2, and 4° PRTS stimuli. Some VL subjects were unable to complete a 4° trial for a given test condition (Table 2). If this occurred, a 0.5° stimulus was substituted so that results were obtained for each test condition at a minimum of three different stimulus amplitudes. There were three testing sessions for each VL subject lasting about 2.5 h each. Two sessions were performed on the same day (morning and afternoon), and a third session was completed on another day. A 5-min break was given between trials.

All subjects were instructed to maintain a relaxed upright stance position. If a subject fell during a trial, the trial was immediately repeated once or twice. Subjects wore headphones and listened to audio tapes of novels and short stories to mask equipment sounds, maintain alertness, and distract them from concentrating on their balance control.

Transfer function analysis

Results were analyzed primarily by calculating transfer function and coherence function estimates from stimulus-response data from each test trial. A transfer function characterizes the dynamic behavior of a system by showing how response sensitivity (gain) and timing (phase) change as a function of stimulus frequency. At each frequency, the transfer function gain gives the ratio of the amplitude of the response (COM body-sway angle) to the stimulus amplitude (support surface and/or visual surround tilt angle) at that frequency. The transfer function phase (in degrees) expresses the relative timing of the response compared with the stimulus at each frequency.

Transfer functions were computed using a discrete Fourier transform (DFT) to decompose the PRTS stimulus and response waveforms into sinusoidal component parts (Bendat and Piersol 2000). The DFT was applied to each 60.5 s cycle of each trial's stimulus and response waveforms. The DFT was calculated at 150 frequencies ranging from f = 1/60.5 = 0.0165 Hz to f = 150/60.5 = 2.48 Hz. A property of the PRTS stimulus is that all even frequency components have zero amplitude (Davies 1970). These even frequency points were discarded leaving 75 frequency samples.

Various power spectra were computed and averaged over the N cycles (N = 6 or 8)
<IT>G<SUB>x</SUB></IT>(<IT>&ohgr;</IT>)<IT>=</IT><FR><NU><IT>1</IT></NU><DE><IT>N</IT></DE></FR> <LIM><OP>∑</OP><LL><IT>i</IT><IT>=1</IT></LL><UL><IT>N</IT></UL></LIM> <IT>X<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>*·</IT><IT>X<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>) (2)

<IT>G<SUB>y</SUB></IT>(<IT>&ohgr;</IT>)<IT>=</IT><FR><NU><IT>1</IT></NU><DE><IT>N</IT></DE></FR> <LIM><OP>∑</OP><LL><IT>i</IT><IT>=1</IT></LL><UL><IT>N</IT></UL></LIM> <IT>Y<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>*·</IT><IT>Y<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>) (3)

<IT>G<SUB>xy</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>=</IT><FR><NU><IT>1</IT></NU><DE><IT>N</IT></DE></FR> <LIM><OP>∑</OP><LL><IT>i</IT><IT>=1</IT></LL><UL><IT>N</IT></UL></LIM> <IT>X<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>*·</IT><IT>Y<SUB>i</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>) (4)
where Xi(jomega ) and Yi(jomega ) are the DFTs of the ith stimulus and response cycles, respectively, omega  = 2pi f, j is <RAD><RCD>−1</RCD></RAD>, and * indicates a complex conjugate. Higher-frequency portions of these spectra were further smoothed by averaging adjacent spectral points to produce the final smoothed spectra, Gxs(omega ), Gys(omega ), and Gxys(jomega ) at 17 frequencies ranging from 0.0165 to 2.23 Hz. (The number of adjacent points averaged increased with increasing frequency such that 15 points contributed to the highest frequency spectral estimate).

The transfer function was then estimated
<IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>=</IT><FR><NU><IT>G<SUB>xys</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)</NU><DE><IT>G<SUB>xs</SUB></IT>(<IT>&ohgr;</IT>)</DE></FR> (5)
and the gain function (ratio of response amplitude to stimulus amplitude) and phase function (in degrees) were computed by
‖<IT>H</IT>(<IT>&ohgr;</IT>)<IT>‖=</IT><RAD><RCD><IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>*·</IT><IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT>)</RCD></RAD> (6)

∠<IT>H</IT>(<IT>&ohgr;</IT>)<IT>=</IT><FR><NU><IT>180°</IT></NU><DE><IT>&pgr;</IT></DE></FR> tan<SUP><IT>−1</IT></SUP>(Im(<IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT>))<IT>/</IT>Re(<IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT>))) (7)
where Im(H(jomega )) and Re(H(jomega )) are the imaginary and real portions, respectively, of the complex numbers representing transfer function estimates. The phase function was computed using the Matlab (The MathWorks, Natick, MA, version 5.3) function "phase" from the Systems Identification Toolbox. This function "unwraps" the phase values, meaning that phase measures more negative than -180° could be obtained in cases where phase lags were increasing with increasing frequency.

COHERENCE FUNCTION ESTIMATES. The smoothed power spectra were used to estimate a "coherence function" given by
&ggr;<SUP>2</SUP>(&ohgr;)=<FR><NU>‖<IT>G<SUB>xys</SUB></IT>(<IT>j</IT><IT>&ohgr;</IT>)<IT>‖<SUP>2</SUP></IT></NU><DE><IT>G<SUB>xs</SUB></IT>(<IT>&ohgr;</IT>)<IT>·</IT><IT>G<SUB>ys</SUB></IT>(<IT>&ohgr;</IT>)</DE></FR> (8)
Values of the coherence function can vary from 0 to 1, with 0 indicating that there is no linear correlation between the stimulus and response, and 1 indicating a perfect linear correlation with no noise. Values less than 1 occur in practice either because there is noise in the system or there is a nonlinear relation between stimulus and response (Bendat and Piersol 2000). Coherence function estimates were also used in the calculation of 95% confidence limits on the transfer function gain and phase data (Otnes and Enochson 1972).

CURVE FITTING. Transfer function curve fits were made to the experimentally determined transfer function estimates. These curve fits were based on a theoretical model of the postural control system (Fig. 9, Eqs. 14-17) and were used to derive estimates of model parameters. The curve fits were made using a constrained nonlinear optimization algorithm ("constr" from the Matlab Optimization Toolbox) that adjusted model transfer function parameters to minimize the error value E given by
<IT>E</IT><IT>=</IT><LIM><OP>∑</OP><LL><IT>i</IT><IT>=1</IT></LL><UL><IT>P</IT></UL></LIM> <FR><NU><IT>‖</IT><IT>M</IT>(<IT>j</IT><IT>&ohgr;</IT><SUB><IT>i</IT></SUB>)<IT>−</IT><IT>H</IT>(<IT>j</IT><IT>&ohgr;</IT><SUB><IT>i</IT></SUB>)<IT>‖</IT></NU><DE><IT>‖</IT><IT>M</IT>(<IT>j</IT><IT>&ohgr;</IT><SUB><IT>i</IT></SUB>)<IT>‖</IT></DE></FR> (9)
where M(jomega i) and H(jomega i) represent the complex value of the model and experimental transfer function, respectively, at the ith frequency point. The error magnitude at each frequency point was summed over the P frequencies included in the fit. The error was normalized by the magnitude of the model to account for the fact that estimation errors were lower for the higher-frequency lower-gain data due to averaging of spectral data. Although 17 frequency points were calculated for experimental transfer functions, P was typically less than 17 (usually 12). The lowest frequency and highest four frequencies were excluded from the fit procedure because the model transfer function was often unable to account for these data. (Large systematic errors between experimental data and fits occurred if all points were included.) Extensive simulations (using the Fig. 9 model and including variability to represent spontaneous body sway) were performed to validate our data analysis and parameter estimation procedures.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Stimulus-evoked body sway

Examples of COM body sway rotational responses to five different amplitudes of the PRTS stimulus for one representative normal subject and one VL subject are shown in Fig. 3. In the example in Fig. 3B, the subjects were standing with eyes closed, and the PRTS stimulus was applied to the support surface. The COM sway responses clearly followed the general PRTS stimulus waveform (Fig. 3A), indicating that both normal and VL subjects tended to orient to the moving support surface. At the three lowest stimulus amplitudes (0.5, 1, and 2°) for this test condition, the amplitude of body sway was noticeably larger than the stimulus amplitude in both the normal and VL subjects. In the normal subject, the sway amplitude appeared to saturate as the stimulus amplitude increased, such that the body-sway responses to the two highest stimuli (4 and 8°) were clearly smaller than the stimulus. However, in the VL subject, body sway continued to increase with increasing stimulus amplitude, and body-sway amplitude remained noticeably larger than the stimulus at the highest stimulus amplitude tested (4°). A similar pattern of body-sway responses was seen with the PRTS stimulus applied to the visual surround while subjects stood on a sway-referenced support surface (Fig. 3C).



View larger version (33K):
[in this window]
[in a new window]
 
Fig. 3. Example PRTS stimulus and body-sway responses for normal and vestibular loss subjects. A: one cycle of the PRTS stimulus waveform with peak-to-peak amplitude ranging from 0.5 to 8°. B: corresponding body-sway responses of a normal and vestibular loss subject to support surface rotations with eyes closed (---, mean of 6 cycles, grayed area represents the 95% confidence interval about the mean). C: body-sway responses of a normal and vestibular loss subject to visual surround rotations with the support surface sway-referenced.

To illustrate the general pattern of responses to all six test conditions, Fig. 4 plots the mean root-mean-square (rms) COM body sway for normal subjects and individual results for the VL subjects. The plots include results for all six test conditions as a function of the peak-to-peak PRTS stimulus amplitude. The rms amplitude of the PRTS stimulus is shown for reference. In the two test conditions where two or more sensory systems were providing veridical orientation cues (i.e., fixed support surface with visual stimulation and fixed visual surround with support surface stimulation), normal subjects showed reduced rms body sway compared with the other test conditions, where fewer sensory systems were providing veridical orientation cues. In the fixed-support-surface condition with visual stimulation, there was no significant difference among mean rms sway for 0.5, 1, 2, and 4° amplitudes for normal subjects, although rms sway was greater for the 8° stimulus compared with the 0.5 or 1° stimuli (Tukey-Kramer multiple comparisons test, P = 0.05 significance level). In the fixed-visual-surround condition with support surface stimulation, the mean body-sway response increased in proportion to the stimulus for the 0.5 and 1° amplitudes but then showed no further increase with increasing stimulus amplitude (no significant differences among mean rms sway for 1-8° stimulus amplitudes, Tukey-Kramer multiple comparisons test, P > 0.05).



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 4. Stimulus-response relationships for all 6 test conditions. The root-mean-square (rms) center-of-mass (COM) body sway for normal subjects (, ---, mean ± SD), and rms COM body sway of individual vestibular loss subjects (open circle , · · · ) are shown as a function of the peak-to-peak stimulus amplitude. The rms amplitude of the PRTS stimulus is shown for reference (, ---).

In the four conditions where veridical orientation cues were provided by fewer sensory systems, the mean rms sway amplitude was always greater than the rms stimulus amplitude for 0.5 and 1° stimuli. For stimuli greater than 1°, the rms body-sway response did not increase in proportion to the stimulus. For the 2° stimuli, the mean rms sway amplitude was approximately equal to the rms stimulus amplitude. For the 4 and 8° PRTS stimuli, the mean rms sway amplitude was always less than the rms stimulus amplitude. In most conditions, the general pattern of results suggests that the body-sway response saturated, with no further increase in response occurring with PRTS stimuli at or above about 4°. Pairwise multiple comparisons showed no significant difference between 4 and 8° responses in normal subjects (P > 0.05, Tukey-Kramer multiple comparisons test).

In contrast, responses of VL subjects generally showed greater levels of stimulus-evoked body sway for each test condition and at each stimulus amplitude, and some subjects were unable to complete all trials (Table 2). Figure 4 plots the rms response amplitudes of individual VL subjects. With the exception of only two individual trials, the individual rms responses of VL subjects were greater than the mean normal responses on all test conditions and at all stimulus amplitudes. In the four test conditions where the number of sensory systems providing veridical orientation cues was reduced, the rms responses of individual VL subjects were always greater than the rms stimulus amplitude, even at the 4° stimulus amplitude where normal subjects showed rms responses that were always less than the stimulus. Overall, the VL subjects' responses increased with increasing stimulus amplitude and showed little evidence for the response saturation seen in normal subjects.

Adaptation and habituation

There was no evidence for adaptation or habituation in normal subjects in their responses to the PRTS stimuli over the six or eight cycles presented in the different test conditions. When cycle by cycle rms sway data from all six test conditions were grouped according to the stimulus amplitude, rather than a decrement in the response due to habituation, the only trends were small increases in the rms response amplitude over the course of the stimulus. Linear regression analyses showed slopes ranging from 0.002°/cycle for the 1° PRTS data to 0.021°/cycle for the 8° PRTS data. The linear regression slope was significantly different from zero only for the 8° PRTS stimulus data (P < 0.046). VL subjects also did not appear to habituate to the various stimuli.

Postural control dynamics

Body-sway responses to PRTS stimuli were analyzed using Fourier methods to compute stimulus-response transfer functions (gain and phase as a function of stimulus frequency) and coherence functions. In these transfer functions, a gain measure of one and phase of zero indicates that the body orientation remained perfectly aligned (in amplitude and timing) with the support surface (for tests with support surface stimulus), or with the visual surround (for tests with visual surround stimulus). A gain greater (less) than one at a particular component stimulus frequency indicates that the body sway amplitude was larger (smaller) than the stimulus amplitude at that frequency. A gain of zero indicates that the body orientation was not influenced by the stimulus and remain perfectly aligned with earth-vertical.

EXAMPLE TRANSFER FUNCTIONS. Figure 5 shows individual transfer functions and their associated coherence functions from four different normal subjects and four different test conditions. All of these transfer functions are from trials where the subjects were restrained to sway as a single-link inverted pendulum by the backboard assembly. The general pattern of gain and phase changes, as a function of stimulus frequency, was similar for all test conditions and stimulus amplitudes. The gain was largest in the 0.1-0.8 Hz frequency range and was often greater than unity in this range. At frequencies less than about 0.1 Hz, the gain gradually decreased with decreasing frequency. At higher frequencies, the gain typically showed a steep decline with increasing frequency. A prominent phase lead was typically present for stimulus frequencies less than about 0.1 Hz. The phase crossed zero in the 0.1-0.2 Hz region and then showed increasing phase lags with increasing frequency. Phase lags of as much as -400° were seen at the highest frequency (2.23 Hz) on some tests (Fig. 5A), but on other tests, the phase lag at the highest frequency was much less (about -200°, Fig. 5B). Coherence functions typically had values in the 0.6-0.8 range at the lower stimulus frequencies, and the values declined with increasing frequency.



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 5. Example transfer functions and coherence functions from 4 different stimulus conditions in 4 different normal subjects. Transfer function gain data (log scale) and phase data (linear scale) are plotted against stimulus frequency (log scale). Error bars indicate 95% confidence intervals around the gain and phase estimates ( and open circle ). Unity gain and 0 phase responses ( · · · ) signify the result expected if subjects were able to maintain perfect body alignment to the moving visual surround and/or support surface stimulus. ---, fits of Eq. 14 to the transfer function data. Only transfer function data points indicated by  were used for the curve fit procedure. Coherence function estimates as a function of stimulus frequency are shown for each test condition. Results are from backboard supported trials.

The general form of the example transfer functions shown in Fig. 5, A and D, are representative of the majority of experimentally obtained transfer functions. The two transfer functions in Fig. 5, B and C, represent the most extreme deviations from the more common transfer function form. Specifically, these two transfer functions show resonant peaks in their gain functions. In one case, a prominent resonant peak was located at about 0.2 Hz. In the other, a smaller resonant peak was located near 0.9 Hz. Of the 240 tests performed by the eight normal subjects, only about 20 tests showed resonant behavior qualitatively similar to the examples shown in Fig. 5, B and C. Four of the normal subjects did not show resonant behavior on any test. Resonant behavior was most common on test conditions with visual stimulation on a sway-referenced support surface, or with support surface stimulation with eyes closed.

MEAN TRANSFER FUNCTION. Backboard support. The mean transfer functions from eight normal subjects are shown in Fig. 6 for the six different test conditions at each of the five different stimulus amplitudes. These transfer function families illustrate that mean gains were greater than unity for all test conditions at the 0.5 and 1° stimulus amplitudes in the 0.1-0.8 Hz frequency range. As the stimulus amplitude increased, the gain generally decreased. For the two visual stimulus conditions with a fixed or sway-referenced support surface, the decrease in gain with increasing stimulus amplitude was more uniform across the bandwidth of test frequencies compared with the other four conditions with support surface stimulation or combined visual and support surface stimulation. In these four support surface stimulus conditions, there was little or no decrease in gain at the highest test frequencies. However, for frequencies less than about 1 Hz, the decline in gain with increasing stimulus amplitude was similar to that for the two visual stimulus conditions. As a result, the gain functions corresponding to the five stimulus amplitudes for all four support surface stimulus conditions converged at about 2 Hz.



View larger version (37K):
[in this window]
[in a new window]
 
Fig. 6. Mean transfer function gain and phase data, and coherence function data plotted vs. stimulus frequency for the 6 different test conditions and 5 different stimulus amplitudes (mean of data from 8 normal subjects on backboard supported trials).

Of all test conditions and amplitudes, the largest mean gain value occurred on the 0.5° visual stimulus condition with a sway-referenced support surface. This gain value was 4.1 at a frequency of 0.36 Hz. Among the other tests performed using the 0.5° stimulus amplitude, maximum gains of about 3.7 were obtained on three conditions (support surface stimulation with eyes closed and with sway-referenced vision, and combined support surface and visual stimulation). A lower maximum gain of about 2.2 was obtained with the 0.5° support surface stimulation with fixed vision, and the lowest maximum gain of about 1.8 was obtained with visual stimulation on the fixed support surface.

Phase functions also showed characteristic differences between the two conditions with only visual stimulation and the four conditions with support surface or combined support surface and visual stimulation. For the two conditions with only visual stimulation, the phase functions for all five stimulus amplitudes were indistinguishable from one another for frequencies greater than about 0.1 Hz. Large phase lags ranging from -342 to -382° were measured in these two conditions at 2.23 Hz. For frequencies less than 0.1 Hz, there was some divergence of these phase functions. At the lowest frequency (0.017 Hz) on the fixed support surface with visual stimulus condition, there was a systematic relationship between the phase and stimulus amplitude with the largest phase advance associated with the lowest stimulus amplitude, and the least phase advance occurring with the largest stimulus amplitude.

The phase functions for all conditions with support surface stimulation showed a different pattern. For three of the four support surface stimulus conditions, the phase functions from different stimulus amplitudes were indistinguishable from one another for frequencies less than about 0.1 Hz. (The exception was the fixed vision condition were there was some divergence of phases at lower frequencies, although no systematic relationship between phase and stimulus amplitude was present.) At frequencies greater than about 0.2 Hz, the phase functions associated with different stimulus amplitudes diverged from one another, and this divergence increased with increasing frequency. At a given frequency greater than 0.2 Hz, there was less phase lag associated with larger stimulus amplitudes. The largest difference between mean phases was seen on the combined support surface and visual stimulus condition where the mean phase from the 8° stimulus was -212°, and the mean phase from the 0.5° stimulus was -386° at 2.23 Hz.

Conditions that evoked body sway, using support surface rotation with either eyes closed or sway-referenced visual surround conditions, resulted in very similar responses in normal subjects (Figs. 4 and 6). This suggests that eye closure and visual surround sway-referencing were both equally effective in eliminating the contribution of visual sensory cues to postural control.

Coherence functions were generally largest at the lowest test frequencies and gradually declined with increasing frequency. Within each of the six test conditions, the coherence function from the 0.5° stimulus typically had lower values than the coherence functions associated with the larger stimulus amplitudes. This is consistent with a lower signal-to-noise ratio for responses evoked by this very low-amplitude stimulus. For 1-8° stimuli, there was no systematic change in coherence with stimulus amplitude for stimulus frequencies less than about 1 Hz. This indicates that each of the responses to the different amplitude stimuli showed a similar degree of linearity even though the overall system responses were clearly nonlinear since gains declined with increasing stimulus amplitude. This suggests that the overall gain decline may be caused by a resetting of system parameters for different amplitudes rather than by a specific nonlinear processing of sensory or motor signals (see DISCUSSION).

One might question the validity of averaging transfer functions across subjects, based on the concern that important individual variations might be obscured or biases introduced. However, examination of each individual's results showed that they were consistent with the population mean. In addition, quantitative estimates of the stiffness and damping properties of each subject's control behavior (see Model interpretation section) showed that each subject effectively normalized his or her postural transfer function dynamics by systematically setting stiffness and damping properties in proportion to individual body mass and height parameters. Similar self-normalizing behavior has been reported previously in a study of head stabilization (Keshner et al. 1999).

Freestanding. The major features of transfer and coherence functions obtained with subjects restrained to sway as a single-link inverted pendulum by the backboard support were also observed in freestanding subjects. As an example, the mean transfer and coherence functions from the four subjects who completed a full set of freestanding tests are shown in Fig. 7 for the combined support surface and visual stimulation condition. Both backboard and freestanding gain functions show a decline with increasing stimulus amplitude, and both show the characteristic convergence of the gain functions from different stimulus amplitudes at higher stimulus frequencies.



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 7. Comparison of transfer functions and coherence functions obtained in backboard supported (left) and freestanding (right) conditions. Results are from the combined visual and support surface stimulus condition at 5 different stimulus amplitudes. Each curve represents the mean data from 4 normal subjects who performed both backboard and freestanding trials.

Both the backboard and freestanding phase functions show a phase lead at frequencies less than 0.1 Hz. The phase functions from different amplitude stimuli were not distinguishable from one another at these lower frequencies. However, as frequency increased greater than about 0.2 Hz, both the backboard and freestanding phase functions diverged from one another, with phase functions showing less phase lag with increasing stimulus amplitude. Coherence functions were also similar for the backboard and freestanding tests.

There was a similar correspondence between backboard-restrained and freestanding results in the other five test conditions. All of the major features of transfer functions for the six test conditions obtained from backboard-restrained subjects (shown in Fig. 6) were also present in the freestanding results.

VL TRANSFER FUNCTIONS. Unlike the normal subject transfer functions, transfer functions derived from VL subjects showed specific characteristic variations between the subjects. These variations precluded any averaging of transfer function results across VL subjects. However, there were some consistent findings across all VL subjects. The most obvious and consistent difference between results from normal and VL subjects was that the transfer function gain functions showed little change with increasing PRTS stimulus amplitude. Figure 8A overlays four transfer functions obtained from one representative VL subject in response to support surface stimulation with eyes closed and PRTS stimulus amplitude varying from 0.5 to 4°. A small decrease in gain can be appreciated by comparing the 0.5° with the 4° data in the lower and mid-frequency range, but this decrease is small compared with the approximate fourfold gain decrease in normal subjects on the same test conditions (see Fig. 6). The phase functions in Fig. 8A also show that there was little change in phase with increasing stimulus amplitude. This is in contrast to the phase data from normal subjects who showed less phase lag with increasing stimulus amplitude in this test condition (Fig. 6).



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 8. Example transfer function gain and phase data and coherence function data vs. stimulus frequency from 3 VL subjects. A: results from subject VL3 showing little change in gain and phase functions on tests performed with eyes closed and support surface stimulus amplitudes varying from 0.5 to 4°. B: results from subject VL3 showing a decrease in gain and some change in phase functions when the support surface stimulus amplitude increased from 1 to 4° on tests performed with a fixed visual surround. C: results from subject VL1 performing the same tests as in B but showing little change in gain and phase functions as the stimulus amplitudes varied from 1 to 4°. D: results from subject VL2 showing large mid-frequency gains in response to a 1° support surface stimulus with a sway-referenced visual surround. ---, from a fit of Eq. 14 to the transfer function data. Error bars indicate 95% confidence intervals. Only transfer function data points indicated () were used for the curve fit procedure. All results are from backboard supported trials.

Results similar to those in Fig. 8A (showing little or no gain reduction or phase changes with increasing stimulus amplitude) were obtained in all four VL subjects with support surface stimulation with eyes closed or sway-referenced vision, with combined support surface and visual stimulation, and with visual stimulation on a sway-referenced support surface. In contrast to this pattern of results, different VL subjects showed varying amounts of reduction in gain and changes in phase with increasing stimulus amplitude during support surface stimulation with a fixed visual surround. The data from the subject showing the most change are shown in Fig. 8B. For comparison, the data from the VL subjects who showed the least change in this same test condition are shown in Fig. 8C.

The VL subject's transfer function data shown in Fig. 8C also illustrate a property of this subject's responses that distinguished her from the other three VL subjects. Note that all three of the transfer functions from this subject in Fig. 8C are very similar to the transfer function of the other VL subject (Fig. 8B) obtained in response to the 4° PRTS stimulus. That is, the gain functions in Fig. 8C are relatively flat in the low- and mid-frequencies before beginning to decrease at higher frequencies. Additionally, the maximum phase lag in Fig. 8C at the highest frequency is about -280°, which is approximately equal to the phase lag in Fig. 8B at the highest frequency for the 4° stimulus. The pattern of results for the VL subject shown in Fig. 8C (subject VL1) will later be shown to be consistent with this subject maintaining a higher level of postural stiffness than other VL subjects and other normal subjects. This subject's high stiffness strategy was evident in all test conditions.

In the test condition with visual stimulation on a fixed support surface, one of the VL subjects (VL2) showed some reduction in gain with increasing stimulus amplitude, whereas the other three VL subjects did not.

VL subjects generally had more difficulty maintaining balance in the test condition where the visual surround was sway-referenced than in the eyes closed condition. Only one of the four subjects was able to complete the 4° support surface stimulus with a sway-referenced visual surround while three of the four were able to complete the 4° stimulus with eyes closed. Some of the highest gains were obtained on tests with support surface stimulation with sway-referenced vision. An example is shown in Fig. 8D (subject VL2). Subject VL4 was unable to complete any of the tests with the visual surround sway-referenced, including the control test with a fixed support surface, although her performance on all five other test conditions was not distinguishable from the other VL subjects. This poor performance on tests with sway-referenced vision is in contrast to the results from normal subjects where performance on tests with eyes closed was indistinguishable from tests with sway-referenced vision (Fig. 6).

Model interpretation

INDEPENDENT CHANNEL MODEL. To achieve further insight into the postural control behavior revealed in the experimental results, we used a simple control system model to parameterize our transfer function results. This "independent channel" model (shown in Fig. 9) uses negative feedback control from sensory systems to generate an active torque, Ta, and muscle/tendon stretch caused by body sway relative to the support surface to generate a passive torque, Tp. Ta and Tp sum to produce the total corrective torque, Tc, acting about the ankle joint. This model assumes that three types of sensory information contribute to the generation of Ta. The sensory information is provided by the visual system, detecting body orientation with respect to the visual environment, proprioceptive system, detecting body orientation with respect to the support surface, and graviceptive system, detecting body orientation with respect to earth vertical. All sensory channels act separately, and the "goal" of each channel is to minimize deviation from its individual internal reference.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 9. "Independent channel" model of sensory integration in postural control showing a weighted addition of contributions from visual, proprioceptive, and graviceptive systems to the generation of an active corrective torque, Ta, as well as a passive torque contribution, Tp, related to body movement relative to the feet.

The body rotation produced by Tc is described by the differential equation of a single-link inverted pendulum
<IT>J</IT> <FR><NU><IT>d</IT><SUP><IT>2</IT></SUP><IT>BS</IT></NU><DE><IT>dt</IT><SUP><IT>2</IT></SUP></DE></FR><IT>=</IT><IT>mgh</IT><IT> sin</IT>(<IT>BS</IT>)<IT>+</IT><IT>T</IT><SUB><IT>c</IT></SUB> (10)
where J is the body's moment of inertia about the ankle joint axis, m is body mass (not including the feet), g is acceleration due to gravity, h is the height of the COM above the ankle joint, and BS is the angular position of the body with respect to earth vertical. Both BS and Tc are time varying. Expressing this equation as a Laplace transform gives
<FR><NU>BS(<IT>s</IT>)</NU><DE><IT>T</IT><SUB><IT>c</IT></SUB>(<IT>s</IT>)</DE></FR>=<FR><NU>1</NU><DE><IT>Js</IT><SUP><IT>2</IT></SUP><IT>−</IT><IT>mgh</IT></DE></FR> (11)
where s is the Laplace variable and the small angle approximation BS approx  sin(BS) was used to simplify the equation.

An inverted pendulum is inherently unstable because a small deviation of body position from a perfect upright position produces a torque due to gravity, mgh sin(BS), that accelerates the body further from upright. To counteract this gravitational torque, Tc must be generated with a sign opposite to the gravitational torque. The stability and dynamic behavior of the inverted pendulum control depends on the time course of Tc. If we assume initially that the passive torque contribution to Tc is negligible, then all of the stabilizing torque must be derived from the sensory error signal, e, given by
<IT>e</IT>(<IT>t</IT>)<IT>=</IT>−<IT>w</IT><SUB><IT>g</IT></SUB><IT>BS</IT>(<IT>t</IT>)<IT>−</IT><IT>w</IT><SUB><IT>v</IT></SUB>[<IT>BS</IT>(<IT>t</IT>)<IT>−VS</IT>(<IT>t</IT>)]<IT>−</IT><IT>w</IT><SUB><IT>p</IT></SUB>[<IT>BS</IT>(<IT>t</IT>)<IT>−FS</IT>(<IT>t</IT>)] (12)
It is known from control system theory that stabilization of an inverted pendulum requires that Ta minimally contain two components: one proportional to e (a "stiffness" factor) and a second proportional to the time derivative of e (a "damping" factor). An optional third component contributes a torque proportional to the mathematical time integral of e. This integral component adds low frequency error correcting properties to the overall control system but is not necessary for stability. The equation for a controller that contains all three of these components is
<IT>T</IT><SUB><IT>a</IT></SUB>(<IT>t</IT>)<IT>=</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>e</IT><IT>+</IT><IT>K</IT><SUB><IT>D</IT></SUB> <FR><NU><IT>de</IT></NU><DE><IT>dt</IT></DE></FR><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB> <LIM><OP>∫</OP></LIM> <IT>e dt</IT> (13)
where t is time, and KP, KD, and KI are gain constants that determine the magnitude of the position, velocity, and integral components, respectively. This type of controller, commonly used in man-made control systems, is referred to as a PID controller (proportional, integral, derivative control).

If all of the sensory systems are assumed to have no dynamic behavior over the bandwidth of movement associated with body sway, then the overall transfer function relating body sway evoked by either visual surround and/or support surface motion is given by
<FR><NU>BS(<IT>s</IT>)</NU><DE><IT>VS</IT>(<IT>s</IT>)</DE></FR><IT> or </IT><FR><NU><IT>BS</IT>(<IT>s</IT>)</NU><DE><IT>FS</IT>(<IT>s</IT>)</DE></FR><IT>=</IT><FR><NU><IT>W</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></NU><DE><IT>Js</IT><SUP><IT>3</IT></SUP><IT>−</IT><IT>mghs</IT><IT>+</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></DE></FR> (14)
where VS is the rotational position of the visual surround in space, FS is the foot-in-space rotational position (equal to the support surface rotational angle), and tau d is a time delay that includes sensory transduction, neural processing, transmission, and muscle activation delays.

Passive torque may not be negligible. Tp could contribute to Tc in test conditions where there is body movement relative to the feet (i.e., where BF, defined in Fig. 9, varies over time). This passive torque is assumed to contain both a stiffness component, K, and a damping component, B, and it acts without a time delay. The overall transfer function for conditions with support surface stimulation or combined visual and support surface stimulation is given by
<FR><NU>BS(<IT>s</IT>)</NU><DE><IT>FS</IT>(<IT>s</IT>)</DE></FR><IT>=</IT><FR><NU><IT>Bs</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>Ks</IT><IT>+</IT><IT>W</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></NU><DE><IT>Js</IT><SUP><IT>3</IT></SUP><IT>+</IT><IT>Bs</IT><SUP><IT>2</IT></SUP><IT>+</IT>(<IT>K</IT><IT>−</IT><IT>mgh</IT>)<IT>s</IT><IT>+</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></DE></FR> (15)
and for the condition with visual stimulation on a fixed support surface
<FR><NU>BS(<IT>s</IT>)</NU><DE><IT>VS</IT>(<IT>s</IT>)</DE></FR><IT>=</IT><FR><NU><IT>W</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></NU><DE><IT>Js</IT><SUP><IT>3</IT></SUP><IT>+</IT><IT>Bs</IT><SUP><IT>2</IT></SUP><IT>+</IT>(<IT>K</IT><IT>−</IT><IT>mgh</IT>)<IT>s</IT><IT>+</IT>(<IT>K</IT><SUB><IT>D</IT></SUB><IT>s</IT><SUP><IT>2</IT></SUP><IT>+</IT><IT>K</IT><SUB><IT>P</IT></SUB><IT>s</IT><IT>+</IT><IT>K</IT><SUB><IT>I</IT></SUB>)<IT>e</IT><SUP><IT>−&tgr;<SUB>d</SUB>s</IT></SUP></DE></FR> (16)
In Eqs. 14-16, the parameter W is a factor that influences the overall gain of the transfer function but does not influence the phase. We refer to W as the "sensory weighting factor" because W is a function of the sensory channel weights in the different test conditions. Specifically, in the condition with support surface stimulation with eyes closed or sway-referenced vision, W = wp = 1 - wg, where wp and wg are the proprioceptive and graviceptive channel weights, respectively. (We assume that the sum of all sensory channels that are contributing to balance control is unity in steady state conditions, in this case wp + wg = 1.) For the test condition with visual stimulation on a sway-referenced support surface, W = wv = 1 - wg. For visual stimulation on a fixed support surface, W = wv = 1 - wpwg. For support surface stimulation with fixed visual surround, W = wp = 1 - wvwg. Finally, for combined support surface and visual stimulation, W = wv + wp = 1 - wg. Therefore, estimates of W obtained from transfer function data can be used to derive estimates of the relative contributions made by different sensory systems to balance control in different test conditions with different stimulus amplitudes.

Although the Fig. 9 model represents the output of the sensory systems as position signals, this is not meant to imply that the postural control system makes use of only position information. This model could be equivalently drawn to include both position and velocity signals (as well as other functions of body motion) derived from the sensory systems, and the active corrective torque would then be generated by appropriately scaling and summing the individual motion signals.

Despite the various simplifying assumptions, the independent channel model can explain experimental transfer function results over a wide range of test frequencies, as illustrated by the curve fits of Eq. 14 to transfer function data shown in Figs. 5 and 8D. These curve fits show that the experimental gain and phase data over a frequency range of about 0.05-1.2 Hz are consistent with Eq. 14, including results from some tests that showed prominent resonant properties. Limitations of this model were evident in that gain and phase results below and above this frequency range were not always well fit by Eq. 14. These limitations were not overcome by inclusion of passive properties (Eqs. 15 or 16).

Curve fits were made to all individual transfer function curves of normal and VL subjects and to the average transfer function data from normal subjects using Eq. 14 (no passive torque) and Eqs. 15 or 16 (including passive torque). To limit the number of unconstrained fit parameters, we used anthropometric estimates of J and h (Winter 1990) and direct measurement of m (less 2.6% of total m to account for mass of the feet). For curve fits to transfer function data obtained using the backboard restraint, the values of J, h, and m included the backboard contribution.

Curve fits using Eq. 14 reliably converged to single solutions. Results from these curve fits are discussed in the following text. Later we describe estimates of passive stiffness, K, and