Skip to main content

Is there a trade-off between economy and task goal variability in transfemoral amputee gait?

Abstract

Background

Energy cost minimization has been widely accepted to regulate gait. Optimization principles have been frequently used to explain how individuals adapt their gait pattern. However, there have been rare attempts to account for the role of variability in this optimization process. Motor redundancy can enable individuals to perform tasks reliably while achieving energy optimization. However, we do not know how the non-goal-equivalent and goal-equivalent variability is regulated. In this study, we investigated how unilateral transfemoral amputees regulate step and stride variability based on the task to achieve energy economy.

Methods

Nine individuals with unilateral transfemoral amputation walked on a treadmill at speeds of 0.6, 0.8, 1.0, 1.2 and 1.4 m/s using their prescribed passive prostheses. We calculated the step-to-step and stride-to-stride variability and applied goal equivalent manifold (GEM) based control to decompose goal-equivalent and non-goal-equivalent manifold. To quantify the energy economy, the energy recovery rate (R) was calculated based on potential energy and kinetic energy. Comparisons were made between GEM variabilities and commonly used standard deviation measurements. A linear regression model was used to investigate the trade-off between R and GEM variabilities.

Results

Our analysis shows greater variability along the goal-equivalent manifold compared to the non-goal-equivalent manifold (p < 0.001). Moreover, our analysis shows lower energy recovery rate for amputee gait compared to nonamputee gait (at least 20% less at faster walking speed). We found a negative relationship between energy recovery rate and non-goal-equivalent variability. Compared to the standard deviation measurements, the variability decomposed using GEM reflected the preferred walking speed and the limitation of the passive prosthetic device.

Conclusion

Individuals with amputation cleverly leverage task redundancy, regulating step and stride variability to the GEM. This result suggests that task redundancy enables unilateral amputees to benefit from motor variability in terms of energy economy. The differences observed between prosthetic step and intact step support the development of prosthetic limbs capable of enhancing positive work during the double support phase and of powered prosthesis controllers that allow for variability along the task space while minimizing variability that interferes with the task goal. This study provides a different perspective on amputee gait analysis and challenges the field to think differently about the role of variability.

Introduction

Energy cost minimization has been widely accepted to regulate gait. Experiments [1, 2] and computational models [3, 4] suggest that individuals select their preferred step length [4, 5], step width [6], step frequency [3], and walking speed [1, 2] to minimize the energy cost of walking. Energy optimization characterizes both normal and pathological gait, acting as an important factor in clinical applications [7, 8]. Individuals with gait impairments operate close to their maximum level of effort and are prone to fatigue even at low walking speeds [9, 10]. Hence, reducing the energy cost of walking is an important goal for rehabilitation [8, 11, 12] and a key metric in assessing the effectiveness of assistive devices such as prostheses, orthoses, and exoskeletons [13, 14].

Optimization principles have been frequently used to describe how individuals adapt their gait pattern. However, there have been rare attempts to assess variability in this gait optimization process. Variability is ubiquitous in motor performance [15, 16]. Yet, a deviation from the average gait pattern might cause a divergence from energy optimality [17, 18]. The fact that variation is inevitable raises a question: Is human variability regulated to assist the energy optimization process? If yes, then how?

Traditionally, movement variability has been linked to noise or error, and a large variability is considered “bad motor control”. This concept might be influenced by traditional statistical and assessment methods of movement variability (e.g., standard deviation) that assume randomness and independence of observations [16, 19]. Opposite to this traditional perspective, variability can be seen to give rise to equifinality, allowing infinite number of ways to perform the same task [20, 21]. Thus, individuals can tolerate movement variability that does not interfere with the task goal without requiring extra effort [22, 23]. These observations suggest that variability and movement redundancy can enable individuals to perform complex tasks reliably and repeatedly, increasing adaptability in motor performance [20, 24].

Researchers have analyzed movement redundancy and variability with different approaches, including uncontrolled manifold (UCM) analysis [24, 25], minimum intervention principle (MIP) [23, 26], and goal equivalent manifold (GEM) [21, 27]. Despite methodological differences, all these approaches decompose the total variability at the body-level into goal-equivalent variability and non-goal-equivalent variability. The goal-equivalent variability is the variability component tangent to the manifold, which does not affect the task goal. The non-goal-equivalent variability is the variability component perpendicular to the manifold, which causes deviation from the task goal. Experiments focusing on quiet standing [24, 28], walking [29, 30], reaching [21, 23], and aiming [27] show that individuals preferentially constrain the non-goal-equivalent variability rather than the goal-equivalent variability, most likely because the latter is not detrimental to the achievement of the task goal. If we assume that non-goal-equivalent variability interferes with the task goal and that this interference requires extra effort to correct, then higher non-goal-equivalent variability should be linked to increased energy cost. Based on this rationale, there may be a trade-off between energy optimization and regulation of task variability.

Understanding energy optimization requires us to determine how the features between the non-goal-equivalent and goal-equivalent variability emerge to reconcile with minimal energy cost. This problem is especially interesting in individuals with unilateral lower-limb amputations, who only have indirect control of the motion of the prostheses. Lower-limb amputees spend more metabolic energy than nonamputee individuals during walking and more proximal amputations are associated with greater metabolic energy cost than distal amputations [12, 31]. The increased energy cost of walking might be related to a lack of ankle push-off power from the prosthesis. Individuals with amputations compensate for the lack of push off power by pulling the thigh forward at the end of stance. This compensatory movement allows for continuous forward propulsion but is highly inefficient [14, 31, 32]. In addition, walking with a passive prosthesis requires greater mechanical work for the step-to-step transition from prosthetic to intact limb, increasing the overall energy cost of walking.

Individuals with unilateral amputation need to adapt at every step to respond to external (i.e., lack of push-off) and internal (i.e., compensatory hip movements) perturbations. Thus, knowing how individuals control the step/stride variability would be beneficial to understand how they achieve optimal energy cost. The purpose of this study is to investigate how individuals with unilateral transfemoral amputation walking with a passive prosthesis regulate step and stride variability. Specifically, we aim to understand how goal-equivalent and non-goal-equivalent variability are regulated with respect to energy economy.

We hypothesize that if individuals with amputations regulate their step/stride variability explicitly to minimize energy cost, then goal-equivalent variability will be larger than non-goal-equivalent variability given that these movement variations do not interfere with task performance. Moreover, we hypothesize that if individuals with amputations regulate their step/stride variability explicitly to minimize energy cost and the variability that interferes with task performance requires extra effort to correct, there will be a trade-off between non-goal-equivalent variability and the energy recovery rate across different walking speeds.

Methods

Experimental protocol and data collection

The data collection protocol has been carefully documented in Hood et al. [33], and data has been shared on Springer Nature, scientific data. The present study examined unilateral amputees walking on the treadmill who comfortably walked at or above 0.8 m/s and were classified as full community ambulators (K3) [33]. Seven male and two female transfemoral amputees (39.56 ± 12.5 years old) were recruited in this study. Participant inclusion criteria included the following: age between 18 and 85 years, at least 6 months after amputation, daily use of their prescribed prosthesis, and ability to walk on a treadmill. Exclusion criteria included serious comorbidities (including musculoskeletal, cardiac, neuromuscular, skin, or vascular conditions) and inability to communicate and/or be understood by investigators. The institutional review board at the University of Utah (IRB00099066) and the US Human Research Protection Office (HRPO) of the US Army Medical Research and Development Command (HSRRB log number: A-19840) approved the study protocol. Participants provided written informed consent to participate in the study. Participants also provided written consent for the use of pictures and videos of the experiments. The details of the demographic data are provided in Table 1. Upon enrollment, subjects were asked about their previous experience with walking on a treadmill. Subjects who reported using the treadmill regularly were asked to report their maximum comfortable treadmill speed and reliance on handrails during treadmill use. If a subject reported little to no experience walking on a treadmill with a prosthesis, they were provided training. Training consisted of the subject walking on the treadmill for 2–5 min intervals with periods of rest in between each interval. During each training session, the experimenter started at the slowest speed of the treadmill (i.e., 0.2 m/s) and slowly increased speed until the subject reported that the speed was their maximum comfortable treadmill speed. Each training session lasted less than 2 h. Weekly training sessions were conducted until no further improvements in comfortable walking speed were observed. Before the start of the data acquisition session, subjects were given time to familiarize with all tested speeds.

Table 1 Subject information

Each participant was asked to walk on a treadmill with their daily passive prostheses at speeds of 0.6, 0.8, 1.0, 1.2 and 1.4 m/s. During acceleration and deceleration of the treadmill, participants were instructed to hold on to the treadmill handrails. After the treadmill reached the constant speed, the participants were encouraged to walk without holding the handrails and then recording started. For each treadmill speed, five trials contained around 10 continuous strides except for TF01 having four trials at 0.6, 0.8 and 1.0 m/s and TF02 having 4 trials at speed 0.6 m/s. Only TF05 and TF06 held the handrails all the times during the 1.4 m/s condition. During walking, a 10-camera motion capture system (Vicon, Oxford, UK) was used to record the reflective markers position at 200 Hz and a split-belt Bertec fully instrumented treadmill (Bertec Co; Columbus, OH, USA) was used to record the bilateral ground reaction forces at 1000 Hz. Reflective markers were placed on the participants following a modified Plug-in-Gait Model [34, 35]. Marker trajectories and ground reaction force data were synchronized, recorded, and pre-processed using Vicon Nexus 2 software. A low-pass Butterworth filter with a cut-off frequency of 6 Hz was applied for the marker trajectories. Inverse dynamic of rotational, translational, and potential energy of 15 segments (head, thorax, pelvis, and left and right forearms, upper arms, hands, thighs, shanks, feet) were calculated using Visual 3D. Anthropometric measurements included body weight, height, and segment lengths were taken from each participant to accurately reconstruct the representative model. For each segment, the COM location, the radius of head, and the base radii of the other segments were estimated based on the anthropometric dimension of the 50 percentile composite subjects from Hanavan et al. [36]. The segment mass for each participant was calculated as a percentage of whole-body mass based on de Leva [37], and the mass of the prosthetic ankle and foot was the actual measured mass.

Data analysis

We applied GEM-based control to decompose the goal-equivalent and non-goal-equivalent variability [30, 38] (details see Fig. 1 and “Methods”). To quantify the energy economy, we calculated the energy recovery rate for each step based on potential energy (\({E}_{p}\)) and kinetic energy (\({E}_{k}\)) [39]. To optimize the recovery of mechanical energy, the \({E}_{p}\) and \({E}_{k}\) curves must have the same shape, be equal in amplitude, and be opposite in phase, as in a pendulum. A higher energy recovery rate indicates more efficient walking (see Fig. 2B and “Methods” for details).

Fig. 1
figure 1

A Schematic depicting the goal equivalent manifold (GEM) for maintaining constant walking speed (v) using data from TF01. The red dot is the preferred mean operating point. Each blue dot represents the combination of stride time and stride length for an amputee participant with walking speed at 1.2 m/s. Dots that lie exactly on the red diagonal line (GEM) achieved the same speed and satisfied the goal. B The time series of δT and δp deviations for the data set shown in A

Fig. 2
figure 2

A Example of \({\mathrm{E}}_{p}\) and \({\mathrm{E}}_{k}\) oscillations are plotted as a function of gait cycle duration. The shaded area indicates the periods of double support. The plot shows an out-of-phase oscillation of kinetic and potential energy, allowing energy exchange to occur between \({\mathrm{E}}_{k}\) and \({\mathrm{E}}_{p}\). B Demonstration the stride and step (prothesis/intact step) period determined in this study. HC is heel contact of intact limb; PHC is heel contact of prosthetic limb. Stride period is determined from HC to the next HC. Prosthetic step is determined from PHC to HC where prosthetic limb leads, intact limb trails, and single support is on the prosthetic limb. Intact step is determined from HC to the PHC where intact limb leads, prosthetic limb trails, and single support is on the intact limb

Stride-to-stride and step-to-step analysis were both performed in this study. For each stride, n, stride time (Tn) was taken as the time between consecutive right heel strikes. Stride length (Ln) was taken as the anterior–posterior distance travelled during each stride using the position of right and left heel markers. Stride speed (Sn) was calculated as Sn = Ln/Tn. Step analysis was calculated using the same methods and taken the time and length between heel contact to the contralateral side heel contact. Prosthetic step was determined from prosthetic heel contact to intact heel contact where the intact limb provided push-off power during the initial double support phase and then transferred to the prosthetic single support phase. Intact step was determined from intact heel contact to the prosthetic heel contact where the prosthesis leg provided push-off power during the initial double support phase and then transferred to the intact leg single support phase. (see Fig. 2B).

Goal equivalent manifold (GEM)

We applied a GEM based method to decompose the goal-equivalent and non-goal-equivalent variabilities on treadmill walking as originally proposed by Dingwell [38]. The primary goal for treadmill walking with speed v is to not walk off the treadmill. Thus, we can define the task goal specifically considered in the sagittal plane. The treadmill walking task is specified by an inequality constraint of the form

$$\left|\sum_{n=\mathrm{0,1},2\dots m}^{m}\left({L}_{n}-v{T}_{n}\right)\right| \le K is true for any m < N,$$
(1)

where Ln and Tn are the stride length and time, respectively, at stride nv is the treadmill speed; K is nominally half the treadmill’s length; m is any number of stride smaller than N; and N is the total number of strides. The simplest strategy to satisfy Eq. (1) is to keep v constant at each step, which was formulated using the goal function:

$${L}_{n}-v{T}_{n}=0 \to \frac{{L}_{n}}{{T}_{n}}=v.$$
(2)

This equation is equivalent to a strategy of matching the treadmill speed at each stride. Due to the body redundancy, there are many possible gait strategies that one could use to satisfy this requirement, including a variety of “drunken” or “silly” walks.

Hence, all [\({L}_{n}\), \({T}_{n}\)] pairs that satisfied Eq. (2) defines the GEM, which was a red solid line in the \({L}_{n}\) versus \({T}_{n}\) plane (see Fig. 1A).

To analyze walking dynamics relative to the GEM, we first normalized each \({L}_{n}\) and \({T}_{n}\) to unit variance (\(\sigma =1\)) to provide an intuitive reference for comparison by dividing its own standard deviation:

$$\stackrel{\sim }{{T}_{n}}={T}_{n} /\sigma \left({T}_{n}\right) and \stackrel{\sim }{{L}_{n}}={L}_{n} /\sigma ({L}_{n}).$$
(3)

This yielded a GEM defined by the dimensionless walking speed: \(\tilde{v }=mean(\stackrel{\sim }{{L}_{n}} /\stackrel{\sim }{{T}_{n}})\). We then define the new coordinates centered at a mean “preferred operating point”: [\(\stackrel{\sim }{{T}_{n}}\)*, \(\stackrel{\sim }{{L}_{n}}\)*] = [mean(\(\stackrel{\sim }{{T}_{n}}\)), \(\tilde{v }\stackrel{\sim }{{T}_{n}}\)*] and re-expressed as \({\stackrel{\sim }{{T}_{n}}}^{^{\prime}}= \stackrel{\sim }{{T}_{n}}\) − \(\stackrel{\sim }{{T}_{n}}\)* and \({\stackrel{\sim }{{L}_{n}}}^{^{\prime}}= \stackrel{\sim }{{L}_{n}}\) − \(\stackrel{\sim }{{L}_{n}}\)*.

Finally, the goal-equivalent deviation along the GEM (δT is the goal-equivalent deviation), and non-goal-equivalent deviations perpendicular to the GEM (δp is the non-goal-equivalent deviation) were calculated as (see Additional file 1: Appendix for the derive process):

$$\left[\begin{array}{c}{\delta }_{T}\\ {\delta }_{p}\end{array}\right]= \frac{1}{\sqrt{1+{v}^{2}}} \left[\begin{array}{cc}1& v\\ -v& 1\end{array}\right] \left[\begin{array}{c}{\stackrel{\sim }{{T}_{n}}}^{^{\prime}}\\ {\stackrel{\sim }{{L}_{n}}}^{^{\prime}}\end{array}\right].$$
(4)

To quantify how variability was distributed relative to the GEM, the standard deviation of each new time series (δT and δp) was calculated (see Fig. 1B).

Assuming that the participant intent is to walk in a funky gait pattern, this intention is nested within the defined goal function and the modification of stride time and stride length due to the deviation of the position or the strategy will be reflected on the quantified goal-relevant and goal-irrelevant variabilities.

Energy recovery rate

The potential and kinetic energy of the whole body, divided unto s segments of mass m, can be measured from the gravitational and the kinetic energy of each segment calculated at each instant of time (t) relative to the frame of reference (Fig. 2):

The whole-body potential energy was calculated as:

$${\mathrm{E}}_{p}(\mathrm{t})= \sum_{i=1}^{s}\left({m}_{i}g{h}_{i}\left(t\right)\right),$$
(5)

where hi and mi, respectively, are the vertical distance of center of mass to the ground and the mass of the ith segment, relative to the frame of global reference; \(g\) is the gravitational acceleration (9.81 m/s).

The whole-body kinetic energy was calculated as:

$${E}_{k}\left(t\right)= \sum_{i=1}^{n}(\frac{1}{2}{m}_{i} {[v}_{xi}^{2}\left(t\right)+{v}_{yi}^{2}\left(t\right)+{v}_{zi}^{2}\left(t\right)]+ \frac{1}{2}{m}_{i} {[k}_{xi}^{2}{\omega }_{xi}^{2}(t)+{k}_{yi}^{2}{\omega }_{yi}^{2}\left(t\right)+{k}_{zi}^{2}{\omega }_{zi}^{2}(t)]),$$
(6)

where \({v}_{i}\), is the linear velocity relative to the frame of the global references in x, y and z axis; \({k}_{i}\) and \({\omega }_{i}\), respectively, are the radius of gyration about the axis through the center of mass and angular velocity of the ith segment (frame of local reference).

Thus, the total energy can be computed as:

$${\mathrm{E}}_{\mathrm{total}}\left(\mathrm{t}\right)= {\mathrm{E}}_{p}\left(\mathrm{t}\right)+{\mathrm{E}}_{k}\left(\mathrm{t}\right).$$
(7)

Over the stride and step period, the energy exchange over the stride and step period was computed according to Winter [39]. The external work (\({w}_{ext}\)) on the COM during the N sample period was computed as:

$${w}_{ext}=\sum_{i=1}^{N}\left(\left|\Delta {E}_{p}+\Delta {E}_{k}\right|\right).$$
(8)

Assuming no energy exchanges between \({E}_{k}\) and \({E}_{p}\), the work (\({w}_{net}\)) done by a segment during N sample periods is:

$${w}_{net}= \sum_{i=1}^{N}\left(\left|\Delta {E}_{p}\right|+\left|\Delta {E}_{k}\right|\right),$$
(9)

where the \(\Delta {E}_{p}\) and \(\Delta {E}_{k}\) is the energy difference between two successive time points.

The energy recovery rate (R) represents the percentage of mechanical energy recovery via exchange between kinetic and potential energy in the COM movement. This is computed as:

$$R=100 \cdot \frac{({{w}_{net}-w}_{ext}) }{{w}_{net}}.$$
(10)

In an ideal energy recovery mechanism, the work associated with changes of potential energy is exactly the same as the work associated with kinetic energy changes, but with different sign: Wp = − Wk. That means that work produced to increase the potential energy can be obtained by reducing the kinetic energy and can again be returned to increase the kinetic energy at the next step-to-step transition. Hence, the larger value of R indicates better energy economy.

Statistics

We first performed one-way repeated ANOVA to examine the commonly used gait features across walking speeds (mean and standard deviation (SD) of stride/step length and time). The trend of the effect on speed was tested using within-subject contrast on a linear, quadratic, cubic, and order 4. To answer the question of the first hypothesis, the differences between the two types of decomposed variabilities (δT and δp) across five walking speeds were tested using 2-way repeated ANOVA (2 variabilities × 5 speeds) in stride, intact step, and prosthesis step. To answer the second question, we performed a simple linear regression of energy recovery rate (R) as the dependent variable and δp and walking speed as the independent variable for each stride, prosthetic step, and intact step (model: Ri = β0 + β1 δpi + β2 speedi + εi). If the coefficient of δp is negative and reaches a significant level, it indicates a trade-off between R and δp. To understand the role of push-off for prosthetic and intact leg in respect to the energy exchange, a follow-up analysis was also conducted. We performed 2-way repeated ANOVA (2 (intact & prosthesis step) × 5 (Speeds)) to compare the energy recovery rate between step types and speeds. The significant level was set at α = 0.05. For the main effects that reached the significance level, the Bonferroni test was used for the post hoc comparisons. The Shapiro–Wilk Test was used to test the normality prior applying the ANOVAs and all the tests did not reach significant difference.

Results

Primary gait features across walking speeds

Figure 3A, B depict the standard deviations (SD) of the length and time of prosthetic step, intact step, and stride, respectively. For the length SD of prosthetic step, intact step, and stride linearly decreased with faster walking speed (Fig. 3C) (prosthetic step: F(4,32) = 11.629, p < 0.001, \(\eta_{p}^{2}\) = 0.592, linear effect F(1,8) = 24.876, p < 0.001, \(\eta_{p}^{2}\) = 0.757; intact step: F(4,32) = 6.826, p < 0.001, \(\eta_{p}^{2}\) = 0.460, linear effect F(1,8) = 22.238, p < 0.001, \(\eta_{p}^{2}\) = 0.735; stride: F(4,32) = 6.646, p < 0.001, \(\eta_{p}^{2}\) = 0.454, linear effect F(1,8) = 14.889, p < 0.01, \(\eta_{p}^{2}\) = 0.650). The time variability of prosthetic step, intact step, and stride also linearly decreased with faster walking speed (Fig. 3D) (prosthetic step: F(4,32) = 51.593, p < 0.001, \(\eta_{p}^{2}\) = 0.866, linear effect F(1,8) = 158.259, p < 0.001, \(\eta_{p}^{2}\) = 0.9527; intact step: F(4,32) = 29.234, p < 0.001, \(\eta_{p}^{2}\) = 0.785, linear effect F(1,8) = 51.545, p < 0.001, \(\eta_{p}^{2}\) = 0.866; stride: F(4,32) = 53.982, p < 0.001, \(\eta_{p}^{2}\) = 0.871, linear effect F(1,8) = 117.232, p < 0.01, \(\eta_{p}^{2}\) = 0.936).

Fig. 3
figure 3

Primary gait parameters. Standard deviations (A, B) for stride and step length (Ln) and time (Tn) as a function of walking speed from 0.6 m/s to 1.4 m/s. Error bars indicate between subject 95% confidence intervals at each speed. The SD of length and the SD of time in stride and step length linearly decreased along with speed

Gait variability decomposed based on GEM

The first hypothesis was whether amputees exhibited greater variability along the GEM (task goal) rather than perpendicular to it. Figure 4A–C, respectively demonstrate that the task relevant (δT) variability was significantly larger than task irrelevant (δP) variability in stride, prosthetic step, and intact step (stride: F(1,8) = 47.885, p < 0.001, \(\eta_{p}^{2}\) = 0.857; prosthetic step: F(1,8) = 158.198, p < 0.001, \(\eta_{p}^{2}\) = 0.952; intact step: F(1,8) = 172.667, p < 0.001, \(\eta_{p}^{2}\) = 0.956). For the interaction effect between speed and variability, there was a significant difference in prosthetic step (F(4,32) = 6.827, p < 0.001, \(\eta_{p}^{2}\) = 0.46) and a marginal level on stride (F(4,32) = 2.670, p = 0.051, \(\eta_{p}^{2}\) = 0.26). The post hoc analyses for δT in prosthetic step showed that at the speed of 1.0 m/s, 1.2 m/s and 1.4 m/s were equal but were significantly higher than the speed at 0.6 and 0.8 m/s. δp in prosthetic step, on the other hand, showed a trend as a U shape that walking speed at 1.0 m/s was significantly smaller than 0.6 m/s and 0.8 m/s and significantly smaller than 1.2 m/s and 1.4 m/s. The post hoc analyses for δT in stride showed that 1.2 m/s was significantly smaller than 1.0 m/s and 1.4 m/s, and δp showed 1.2 m/s was significantly larger than 1.0 m/s and 1.4 m/s.

Energy recovery rate

Figure 5 depicts the energy recovery rate of prosthetic step and intact step. There was a significant interaction between step types and speed (F (4,32) = 5.552, p = 0.002, \(\eta_{p}^{2}\) = 0.40), and the post hoc analysis indicated only in prosthetic step at speeds 1.2 m/s and 1.4 m/s were significantly larger than other speeds. The second hypothesis was whether there is a trade-off between δp and energy recovery rate (negative relationship). We performed a linear regression between δp, speed, and energy recovery rate of prosthetic step and intact step. The results showed the coefficient in prosthetic step and intact step reached significant differences with a negative coefficient (see Table 2).

Fig. 4
figure 4

GEM-Based Decomposition of Gait Variability. A Standard deviations for all δT and δp time series at all 5 walking speeds calculated from each stride period. B Standard deviations for all δT and δp time series at all 5 walking speeds calculated from each prosthetic step. C Standard deviations for all δT and δp time series at all 5 walking speeds calculated from each intact step

Fig. 5
figure 5

reproduced from Willems et al., 1995

Energy recovery rate (R) are plotted as a function of walking speed for healthy adults (hexagram), Intact Step (circles) and Prosthetic Step (square). Adult data were

Table 2 The results of the linear regressions of for energy recovery rate of prosthetic step and intact step as a function of δp and speed (model: Ri = β0 + β1δpi + β2 speedi + εi)

Discussion

The purpose of this study was to investigate how individuals with unilateral transfemoral amputation regulate step and stride variability on the goal-equivalent and non-goal-equivalent manifold to achieve energy economy. We examined the hypothesis that if amputees are regulating variability in terms of task goal, a large portion of goal-equivalent variability would be observed. In addition, we examined whether there would be a trade-off between non-goal-equivalent variability and energy recovery rate across different walking speeds under the assumption that individuals require more energy to correct for the non-goal-equivalent variability.

Looking at the main results of primary gait features, the standard deviation of stride/step length and time linearly decreased with faster walking speeds (Fig. 3A, B). This trend was consistent for both prosthesis steps and intact steps. While the standard deviation quantified the magnitude of deviation from the mean of all strides/steps, the time series of goal-equivalent (δT) and non-goal-equivalent (δp) variabilities exhibited the step-to-step change in respect to the task goal with the temporal order. δT, qualitatively, showed larger amplitudes than the δp (Fig. 1B) and the variabilities of δT were significantly higher than the δp across all speeds (Fig. 4). These results support our first hypothesis and indicate that amputees explicitly regulate step/stride variability tuned to the goal equivalent manifold (GEM). Our findings on unilateral amputees followed the same trend as healthy young adults [38, 40] and elderly [41, 42]. The ability to leverage the task redundancy could be critical when tasks become more demanding. For example, walking on uneven terrain requires individuals to adjust each step rather than change their average gait pattern. In addition, this concept can be used as a guiding principle for designing fault-tolerant controllers for powered prostheses. To achieve this goal, a powered prosthesis controller should minimize the δp and only correct errors distributed on the δp because these errors interfere with task performance. In contrast, errors that vary along with the δT should be tolerated as they do not interfere with the task performance. A prosthesis controller following this principle may improve safety while reducing control effort.

The goal-equivalent (δT) and non-goal-equivalent (δp) variabilities demonstrated a different speed effect compared to the standard deviation of stride/step length and time, and such effects reflect the unique characteristics of amputee gait (Fig. 4). During the prosthetic step, the intact leg regulates the push-off power, which propels the body forward, and the swing movement. Thus, the intact leg largely determines the step time and step length. Our results show that the variability of δT and δp changed with walking speed forming an inverted U shape and U shape, respectively (Fig. 4B). A similar trend was observed with young and elderly healthy adults [38, 40, 43]. This trend suggests that the intact leg can adjust the δT and δp variabilities that tuned to the preferred walking speed for unilateral amputees (around 1 m/s) [44]. Interestingly, the intact step did not show the same adaptability to walking speed. During the intact step, the passive prosthetic leg provides limited push-off power [32, 45] and swing time adaptability. Thus, although microprocessor-controlled knee prostheses can adjust the joint damping to walking speed [46, 47], our results indicate that this adaptation mechanism is not enough to achieve similar task performance to the intact leg. Future work could test whether the U-shaped relationship observed in the intact leg occur in powered prostheses capable of imitating the intact-leg push-off power and swing time adaptability.

In this study, we quantified energy recovery rate at each step based on potential energy (\({E}_{p}\)) and kinetic energy (\({E}_{k}\)) [39]. Previous studies have assessed the effectiveness of the pendulum mechanism during human walking by measuring the fraction of the total mechanical energy that is recovered as a result of the transduction between \({E}_{p}\) and \({E}_{k}\) [48, 49]. To optimize the recovery of mechanical energy, which leads to the most efficient walking pattern, the \({E}_{p}\) and \({E}_{k}\) curves must have the same shape, be equal in amplitude, and be opposite in phase, as in a pendulum. A higher energy recovery rate indicates more efficient walking (see Fig. 2B and “Methods” for details). Previous studies showed that for healthy adults at preferred walking speed (~ 1.6 m/s), as much as 70% of the required external mechanical energy can be recovered due to this energy saving mechanism [48, 50, 51]. The other 30% of external mechanical energy is lost from the system and must be supplied by the muscles [52]. Our study shows that unilateral amputees have lower recovery rate than healthy adults. At 1.4 m/s, the average energy recovery rate during the prosthetic step was about 40%. Moreover, the recovery rate could be up to 60% depending on the amputee individual (see Fig. 5). This result explains why individuals with amputation using a passive prosthesis spend, on average, more energy during walking than nonamputee individuals.

As expected, the prosthetic step had a larger energy recovery rate than the intact step. Interestingly, the speed trends for the prosthetic and intact steps were different. The prosthetic step showed a linearly increase of the energy recovery rate with walking speed, whereas the intact step showed an inverted U-shaped profile. This difference in the speed trends could be due to the different ability of the intact and prosthetic leg to generate positive work during push-off. Previous studies in healthy adults found an inverted U-shaped relationship, although for a wider range of walking speeds than the one used in this study (i.e., 0.5–2.5 m/s vs. 0.8–1.4 m/s) [48, 50, 51]. In these studies, energy recovery rate gradually increased from 0.6 m/s to 1.6 m/s, which is similar to the walking speed range of our study. This result is in agreement with our finding in the prosthetic side.

Previous studies in nonamputee individuals show that, at low walking speeds, the potential energy changes are larger than the kinetic energy changes (see Fig. 5). When walking speed increases, the kinetic energy changes increase relative to the potential energy changes, exceeding the potential energy changes after the walking speed is above the preferred walking speed. Thus, the energy recovery rate reaches its maximum for the self-selected speed and it is lower than its maximum both above and below the self-selected speed [48]. Our study shows that when the intact leg is in charge of forward propulsion (i.e., prosthesis step), the amplitude of kinetic energy continuously increases along with walking speed. However, when the prosthetic leg is in charge of forward propulsion (i.e., intact-leg step), the capacity of increasing the amplitude of kinetic energy relative to the potential energy reaches a limit after the walking speed is above the amputee’s preferred walking speed. We speculate that this inability to increase the amplitude of kinetic energy might be due to the passive prosthesis generating less power than the intact leg during push-off, resulting in net negative work during the intact step [32, 45, 53]. The increase of energy dissipation during double support phase at fast speed might explain the inverted U-shaped relationship in a smaller range of walking speed. We speculate that, with proper control of propulsive torque, powered prostheses [14, 45] could enable amputees to achieve better energy recovery rate. Future studies should investigate the step-to-step change of energy recovery rate with powered prostheses to test this hypothesis.

The second hypothesis was tested by applying linear regression to the time series of δp and R calculated extracted from each stride and step. Using this regression, we found that adding the δp in the linear model significantly explains the energy recovery rate. Moreover, the negative coefficient of the linear regression supports the trade-off relationship between δp variability and recovery rate. Since the fitting considered the temporal order with δp and the energy recovery rate, this result demonstrates that individuals with unilateral amputations leverage movement redundancy to correct small deviations at each step to control the task variability while minimizing energy cost. This short-term control mechanism (step-to-step) exploits the inherent task redundancy enabling unilateral amputees to receive potential benefits from motor variability in achieving energy economy. This result could inform a new assessment method for the functionality of powered lower-limb prostheses. For example, better prosthesis designs would be reflected in the GEM analysis by showing a larger δT than δp as it indicates a better energy recovery rate. This assessment method could be the object of future work for rehabilitation engineers.

We recognize that the energy recovery rate might not accurately estimate the energy cost of walking, especially for fast walking speeds [32]. Although measuring the oxygen consumption can lead to a more accurate estimate of the energy cost of walking, it cannot provide step-to-step cost, which is needed for the proposed analysis of GEM variabilities. Using more elaborate musculoskeletal models might provide a better estimation of the energy cost of walking; however, inverse dynamics calculations have been shown not to solve this problem easily [54] and the current validity of more elaborate musculoskeletal models does not warrant success with such an approach either. In addition, previous studies have shown that the effect of speed on the oxygen consumption showing an inverted-U shape turned to the amputees’ preferred walking speed [10, 55]. Considering the average \({\delta }_{P}\) in Fig. 4, we postulated that under a longer time-scale, the trade-off between \({\delta }_{P}\) and oxygen consumption will still hold at least for the prosthetic side. We believe the interaction of task variability and energy recovery rate could provide some clarification on the issue. Further interpretation would be highly speculative and require more studies.

In this study, we defined the goal function under the rationale of treadmill walking with a constant speed. However, one of the limitations of this study is that the goal (maintaining walking speed) may not be what the participants intend to do; instead, it is the necessary criteria to perform the task. This characteristic provides us a simple way to capture the goal equivalent manifold as well as to decompose the variabilities. Assuming the participant has the intent to walk in a “funky” gait pattern, this intention is nested within the defined goal function and the modification of stride time and stride length due to the deviation of the position or the strategy will be reflected on the quantified goal-equivalent and non-goal-equivalent variabilities. However, we do not know if the task goal was indeed the goal adopted by the human subjects. It is an interesting question for a future study to investigate how the task goal and human goal intervene to regulate step and stride variabilities on the GEM. Moreover, treadmill walking is a well-defined task where it is easy to define the task goal. Considering the variety of tasks in clinical settings (e.g., walking, running, stair climbing, etc.), it is not clear how to systematically define these task goals to decompose the goal variabilities, and thus limited the application of the GEM.

In addition, this study recruited nine transfemoral amputees. Even though all of them are K3 level, their demographic data are widely spread (e.g., age, gender, proscribed prosthesis, etc.). Also, we did not record participants’ preferred walking speed and that prevents us to investigate this effect associated with energy cost and variabilities. Thus, future study could expand the sample size and consider the demographic analysis into the model fitting to improve the conclusion and take into account the preferred walking speed into consideration.

To the best of our knowledge, this is the first study to analyze the gait variability at the level of task goal in connection to energy optimization in individuals with lower-limb amputations. Given that human behavior is goal driven, instead of focusing on common gait features such as step mean and standard deviation, we took the task goal into consideration. We believe that the results of this study will provide a different perspective on amputee gait analysis and challenge the field to think differently about the role of variability.

Conclusion

This study investigated how individuals with unilateral transfemoral amputation walking with a passive prosthesis regulate step and stride variability on the goal-equivalent and non-goal-equivalent manifolds. Our results suggest that individuals with amputation cleverly leverage task redundancy, regulating step and stride variability to the goal equivalent manifold (GEM). This result suggests that task redundancy enables unilateral amputees to benefit from motor variability in terms of energy economy. The differences observed between prosthetic step and intact step support the development of prosthetic leg capable of enhancing positive work during the double support phase. Moreover, the results of this study motivate the development of powered prosthesis controllers that allow for variability along the task space while minimizing variability that interferes with the task goal. This study provides a different perspective on amputee gait analysis and challenge the field to think differently about the role of variability.

Availability of data and materials

The data has published at scientific data.

References

  1. Browning RC, Baker EA, Herron JA, Kram R. Effects of obesity and sex on the energetic cost and preferred speed of walking. J Appl Physiol. 2006;100(2):390–8.

    Article  PubMed  Google Scholar 

  2. Ralston HJ. Energy-speed relation and optimal speed during level walking. Int Z Angew Physiol Einschl Arbeitsphysiol. 1958;17(4):277–83.

    CAS  Google Scholar 

  3. Zarrugh M, Todd F, Ralston H. Optimization of energy expenditure during level walking. Eur J Appl Physiol. 1974;33(4):293–306.

    Article  CAS  Google Scholar 

  4. Minetti AE, Alexander RM. A theory of metabolic costs for bipedal gaits. J Theor Biol. 1997;186(4):467–76.

    Article  CAS  PubMed  Google Scholar 

  5. Sekiya N, Nagasaki H, Ito H, Furuna T. Optimal walking in terms of variability in step length. J Orthop Sports Phys Ther. 1997;26(5):266–72.

    Article  CAS  PubMed  Google Scholar 

  6. Donelan JM, Kram R, Kuo AD. Mechanical and metabolic determinants of the preferred step width in human walking. Proc Biol Sci. 2001;268(1480):1985–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Gordon NF, Gulanick M, Costa F, Fletcher G, Franklin BA, Roth EJ, et al. Physical activity and exercise recommendations for stroke survivors: an American Heart Association scientific statement from the Council on Clinical Cardiology, Subcommittee on Exercise, Cardiac Rehabilitation, and Prevention; the Council on Cardiovascular Nursing; the Council on Nutrition, Physical Activity, and Metabolism; and the Stroke Council. Circulation. 2004;109(16):2031–41.

    Article  PubMed  Google Scholar 

  8. Waters RL, Mulroy S. The energy expenditure of normal and pathologic gait. Gait Posture. 1999;9(3):207–31.

    Article  CAS  PubMed  Google Scholar 

  9. Stoquart G, Detrembleur C, Lejeune TM. The reasons why stroke patients expend so much energy to walk slowly. Gait Posture. 2012;36(3):409–13.

    Article  CAS  PubMed  Google Scholar 

  10. Genin JJ, Bastien GJ, Franck B, Detrembleur C, Willems PA. Effect of speed on the energy cost of walking in unilateral traumatic lower limb amputees. Eur J Appl Physiol. 2008;103(6):655.

    Article  PubMed  Google Scholar 

  11. da Cunha IT, Lim PA, Qureshy H, Henson H, Monga T, Protas EJ. Gait outcomes after acute stroke rehabilitation with supported treadmill ambulation training: a randomized controlled pilot study. Arch Phys Med Rehabil. 2002;83(9):1258–65.

    Article  PubMed  Google Scholar 

  12. Waters R, Perry J, Antonelli D, Hislop H. Energy cost of walking of amputees: the influence of level of amputation. J Bone Joint Surg Am. 1976;58(1):42–6.

    Article  CAS  PubMed  Google Scholar 

  13. Zhang J, Fiers P, Witte KA, Jackson RW, Poggensee KL, Atkeson CG, et al. Human-in-the-loop optimization of exoskeleton assistance during walking. Science. 2017;356(6344):1280–4.

    Article  CAS  PubMed  Google Scholar 

  14. Herr HM, Grabowski AM. Bionic ankle–foot prosthesis normalizes walking gait for persons with leg amputation. Proc R Soc B Biol Sci. 2012;279(1728):457–64.

    Article  Google Scholar 

  15. Bernstein NA. Dexterity and its development. UK: Psychology Press; 2014.

    Book  Google Scholar 

  16. Newell KM, Corcos DM. Variability and motor control. USA: Human Kinetics Publishers Champaign; 1993.

    Google Scholar 

  17. Wong JD, Selinger JC, Donelan JM. Is natural variability in gait sufficient to initiate spontaneous energy optimization in human walking? J Neurophysiol. 2019;121(5):1848–55.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Malatesta D, Simar D, Dauvilliers Y, Candau R, Borrani F, Préfaut C, et al. Energy cost of walking and gait instability in healthy 65-and 80-yr-olds. J Appl Physiol. 2003;95(6):2248–56.

    Article  PubMed  Google Scholar 

  19. Stergiou N, Decker LM. Human movement variability, nonlinear dynamics, and pathology: is there a connection? Hum Mov Sci. 2011;30(5):869–88.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Latash ML, Scholz JP, Schöner G. Motor control strategies revealed in the structure of motor variability. Exerc Sport Sci Rev. 2002;30(1):26–31.

    Article  PubMed  Google Scholar 

  21. Cusumano JP, Dingwell JB. Movement variability near goal equivalent manifolds: fluctuations, control, and model-based analysis. Hum Mov Sci. 2013;32(5):899–923.

    Article  PubMed  Google Scholar 

  22. Todorov E. Optimality principles in sensorimotor control. Nat Neurosci. 2004;7(9):907–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Todorov E, Jordan MI. Optimal feedback control as a theory of motor coordination. Nat Neurosci. 2002;5(11):1226–35.

    Article  CAS  PubMed  Google Scholar 

  24. Hsu W-L, Scholz JP, Schöner G, Jeka JJ, Kiemel T. Control and estimation of posture during quiet stance depends on multijoint coordination. J Neurophysiol. 2007;97(4):3024–35.

    Article  PubMed  Google Scholar 

  25. Freitas SM, Scholz JP, Latash ML. Analyses of joint variance related to voluntary whole-body movements performed in standing. J Neurosci Methods. 2010;188(1):89–96.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Liu D, Todorov E. Evidence for the flexible sensorimotor strategies predicted by optimal feedback control. J Neurosci. 2007;27(35):9354–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Cusumano JP, Cesari P. Body-goal variability mapping in an aiming task. Biol Cybern. 2006;94(5):367–79.

    Article  PubMed  Google Scholar 

  28. Hsu W-L, Scholz JP. Motor abundance supports multitasking while standing. Hum Mov Sci. 2012;31(4):844–62.

    Article  PubMed  Google Scholar 

  29. Black DP, Smith BA, Wu J, Ulrich BD. Uncontrolled manifold analysis of segmental angle variability during walking: preadolescents with and without Down syndrome. Exp Brain Res. 2007;183(4):511–21.

    Article  PubMed  Google Scholar 

  30. Dingwell J, Cusumano JP, Cavanagh P, Sternad D. Local dynamic stability versus kinematic variability of continuous overground and treadmill walking. J Biomech Eng. 2001;123(1):27–32.

    Article  CAS  PubMed  Google Scholar 

  31. Quesada RE, Caputo JM, Collins SH. Increasing ankle push-off work with a powered prosthesis does not necessarily reduce metabolic rate for transtibial amputees. J Biomech. 2016;49(14):3452–9.

    Article  PubMed  Google Scholar 

  32. Houdijk H, Pollmann E, Groenewold M, Wiggerts H, Polomski W. The energy cost for the step-to-step transition in amputee walking. Gait Posture. 2009;30(1):35–40.

    Article  PubMed  Google Scholar 

  33. Hood S, Ishmael MK, Gunnell A, Foreman KB, Lenzi T. A kinematic and kinetic dataset of 18 above-knee amputees walking at various speeds. Sci Data. 2020;7(1):150.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Limited VMS. Plug-in gait reference guide. Vicon Documentation. 2016.

  35. Phinyomark A, Petri G, Ibáñez-Marcelo E, Osis ST, Ferber R. Analysis of big data in gait biomechanics: current trends and future directions. J Med Biol Eng. 2018;38(2):244–60.

    Article  PubMed  Google Scholar 

  36. Hanavan Jr EP. A mathematical model of the human body. Air Force Aerospace Medical Research Lab Wright-patterson AFB OH; 1964.

  37. de Leva P. Adjustments to Zatsiorsky-Seluyanov’s segment inertia parameters. J Biomech. 1996;29(9):1223–30.

    Article  PubMed  Google Scholar 

  38. Dingwell JB, John J, Cusumano JP. Do humans optimally exploit redundancy to control step variability in walking? PLoS Comput Biol. 2010;6(7):e1000856.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Winter DA. Biomechanics and motor control of human movement. USA: Wiley; 2009.

    Book  Google Scholar 

  40. Roerdink M, de Jonge CP, Smid LM, Daffertshofer A. Tightening up the control of treadmill walking: effects of maneuverability range and acoustic pacing on stride-to-stride fluctuations. Front Physiol. 2019;10(257).

  41. Decker LM, Cignetti F, Potter JF, Studenski SA, Stergiou N. Use of motor abundance in young and older adults during dual-task treadmill walking. PLoS ONE. 2012;7(7):e41306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Dingwell JB, Salinas MM, Cusumano JP. Increased gait variability may not imply impaired stride-to-stride control of walking in healthy older adults: Winner: 2013 Gait and Clinical Movement Analysis Society Best Paper Award. Gait Posture. 2017;55:131–7.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Yentes JM, Hunt N, Schmid KK, Kaipust JP, McGrath D, Stergiou N. The appropriate use of approximate entropy and sample entropy with short data sets. Ann Biomed Eng. 2012;41(2):349–65.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Batten HR, McPhail SM, Mandrusiak AM, Varghese PN, Kuys SS. Gait speed as an indicator of prosthetic walking potential following lower limb amputation. Prosthet Orthot Int. 2019;43(2):196–203.

    Article  PubMed  Google Scholar 

  45. Fylstra BL, Lee IC, Huang S, Brandt A, Lewek MD, Huang H. Human-prosthesis coordination: a preliminary study exploring coordination with a powered ankle-foot prosthesis. Clin Biomech. 2020;80:105171.

    Article  Google Scholar 

  46. Michael J. Prosthetic suspensions and components. Atlas of amputations and limb deficiencies: Surgical, prosthetic, and rehabilitation principles 3rd ed Rosemont (IL): American Academy of Orthopaedic Surgeons. 2004:409–25.

  47. Gard SA. The influence of prosthetic knee joints on gait. In: Handbook of human motion. Berlin: Springer International Publishing; 2018. p. 1359–82.

    Chapter  Google Scholar 

  48. Cavagna G, Willems P, Heglund N. The role of gravity in human walking: pendular energy exchange, external work and optimal speed. J Physiol. 2000;528(3):657–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Cavagna GA, Thys H, Zamboni A. The sources of external work in level walking and running. J Physiol. 1976;262(3):639–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Willems PA, Cavagna GA, Heglund NC. External, internal and total work in human locomotion. J Exp Biol. 1995;198(Pt 2):379–93.

    Article  CAS  PubMed  Google Scholar 

  51. Cavagna G, Willems PA, Legramandi M, Heglund NC. Pendular energy transduction within the step in human walking. J Exp Biol. 2002;205(21):3413–22.

    Article  CAS  PubMed  Google Scholar 

  52. Cavagna GA, Heglund NC, Taylor CR. Mechanical work in terrestrial locomotion: two basic mechanisms for minimizing energy expenditure. Am J Physiol Regulat Integr Comparat Physiol. 1977;233(5):R243–61.

    Article  CAS  Google Scholar 

  53. Kuo AD. Energetics of actively powered locomotion using the simplest walking model. J Biomech Eng. 2002;124(1):113–20.

    Article  PubMed  Google Scholar 

  54. Gitter A, Czerniecki J, Weaver K. A reassessment of center-of-mass dynamics as a determinate of the metabolic inefficiency of above-knee amputee ambulation. Am J Phys Med Rehabil. 1995;74(5):332–8.

    Article  CAS  PubMed  Google Scholar 

  55. Hoffman MD, Sheldahl LM, Buley KJ, Sandford PR. Physiological comparison of walking among bilateral above-knee amputee and able-bodied subjects, and a model to account for the differences in metabolic cost. Arch Phys Med Rehabil. 1997;78(4):385–92.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This work was partly supported by the National Institutes of Health [NIH EB024570, R01HD098154] and by the Department of the Army [W81XWH-16-1-0701].

Author information

Authors and Affiliations

Authors

Contributions

IL and BLF processed the data. TL directed the data collection. ML and IL discussed the results. IL drafted the manuscript. HH, ML, BLF and TL revised the manuscript. All authors read and approved the final manuscript prior to publication.

Corresponding author

Correspondence to I-Chieh Lee.

Ethics declarations

Ethics approval and consent to participate

The experimental protocol for this study was approved by the University of Utah Institutional Review Board.

Consent for publication

All subjects provided written informed consent, including written permission to publish photos and videos of the experiment.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: 

Derivation of the Equation 4.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lee, IC., Fylstra, B.L., Liu, M. et al. Is there a trade-off between economy and task goal variability in transfemoral amputee gait?. J NeuroEngineering Rehabil 19, 29 (2022). https://doi.org/10.1186/s12984-022-01004-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12984-022-01004-8

Keywords