Analysis of gait data in manifold space: The control of the centre of mass during gait.

Similar documents
Megan E. Krause, BSBSE, Young Hui Chang, Ph.D. Comparative Neuromechanics Laboratory. Georgia Institute of Technology

An investigation of kinematic and kinetic variables for the description of prosthetic gait using the ENOCH system

video Purpose Pathological Gait Objectives: Primary, Secondary and Compensatory Gait Deviations in CP AACPDM IC #3 1

The Influence of Load Carrying Modes on Gait variables of Healthy Indian Women

RUNNING SHOE STIFFNESS: THE EFFECT ON WALKING GAIT

Analysis of Foot Pressure Variation with Change in Stride Length

Purpose. Outline. Angle definition. Objectives:

Normal and Abnormal Gait

Assessments SIMPLY GAIT. Posture and Gait. Observing Posture and Gait. Postural Assessment. Postural Assessment 6/28/2016

The overarching aim of the work presented in this thesis was to assess and

Supplementary Figure S1

Impact Points and Their Effect on Trajectory in Soccer

Posture influences ground reaction force: implications for crouch gait

Complex movement patterns of a bipedal walk

Spasticity in gait. Wessex ACPIN Spasticity Presentation Alison Clarke

Biomechanical analysis of the medalists in the 10,000 metres at the 2007 World Championships in Athletics

Using GPOPS-II to optimize sum of squared torques of a double pendulum as a prosthesis leg. Abstract

Available online at ScienceDirect. Procedia Engineering 112 (2015 )

INTERACTION OF STEP LENGTH AND STEP RATE DURING SPRINT RUNNING

Dynamically stepping over large obstacle utilizing PSO optimization in the B4LC system

GROUND REACTION FORCE DOMINANT VERSUS NON-DOMINANT SINGLE LEG STEP OFF

A COMPARISON OF SELECTED BIOMECHANICAL PARAMETERS OF FRONT ROW SPIKE BETWEEN SHORT SET AND HIGH SET BALL

C-Brace Orthotronic Mobility System

COMPARISON STUDY BETWEEN THE EFFICIENY OF THE START TECHNIQUES IN THE ROMANIAN COMPETITIVE SWIMMING

Gait Analysis at Your Fingertips:

ASSESMENT Introduction REPORTS Running Reports Walking Reports Written Report

Investigation of Bio-Kinematic Elements of Three Point Shoot in Basketball

Toward a Human-like Biped Robot with Compliant Legs

Center of Mass Acceleration as a Surrogate for Force Production After Spinal Cord Injury Effects of Inclined Treadmill Walking

Impact of heel position on leg muscles during walking

Artifacts Due to Filtering Mismatch in Drop Landing Moment Data

Transformation of nonfunctional spinal circuits into functional states after the loss of brain input

Biomechanics and Models of Locomotion

Analysis of ankle kinetics and energy consumption with an advanced microprocessor controlled ankle foot prosthesis.

Brief Biomechanical Analysis on the Walking of Spinal Cord Injury Patients with a Lower Limb Exoskeleton Robot

WORKBOOK/MUSTANG. Featuring: The R82 Next Step Development Plan. mustang. R82 Education

CHAPTER IV FINITE ELEMENT ANALYSIS OF THE KNEE JOINT WITHOUT A MEDICAL IMPLANT

HPW Biomechanics

Gait Analyser. Description of Walking Performance

Gait Analysis of a Little Biped Robot. Received May 2015; accepted July 2015

The effect of different backpack loading systems on trunk forward lean angle during walking among college students

THE ANKLE-HIP TRANSVERSE PLANE COUPLING DURING THE STANCE PHASE OF NORMAL WALKING

Sample Solution for Problem 1.a

The Starting Point. Prosthetic Alignment in the Transtibial Amputee. Outline. COM Motion in the Coronal Plane

Body Stabilization of PDW toward Humanoid Walking

A bit of background. Session Schedule 3:00-3:10: Introduction & session overview. Overarching research theme: CPTA

Ball impact dynamics of knuckling shot in soccer

Sensitivity of toe clearance to leg joint angles during extensive practice of obstacle crossing: Effects of vision and task goal

In memory of Dr. Kevin P. Granata, my graduate advisor, who was killed protecting others on the morning of April 16, 2007.

COMPARISON OF BIOMECHANICAL DATA OF A SPRINT CYCLIST IN THE VELODROME AND IN THE LABORATORY

The Effect of Military Load Carriage on Ground Reaction Forces. Technical Note. Stewart A Birrell 1 Robin H Hooper 1 Roger A Haslam 1

SPRINTING CHARACTERISTICS OF WOMEN S 100 METER FINALS AT THE IAAF WORLD CHAMPIONSHOPS DAEGU 2011

by Michael Young Human Performance Consulting

Figure 1 betois (bending torsion insole system) system with five measuring points and A/D- converter.

APPLICATION OF THREE DIMENSIONAL ACCELEROMETRY TO HUMAN MOTION ANALYSIS

INCLINOMETER DEVICE FOR SHIP STABILITY EVALUATION

+ t1 t2 moment-time curves

KASAMATSU VERSUS TSUKAHARA VAULT

Steffen Willwacher, Katina Fischer, Gert Peter Brüggemann Institute of Biomechanics and Orthopaedics, German Sport University, Cologne, Germany

Mutual and asynchronous anticipation and action in sports as globally competitive

Legendre et al Appendices and Supplements, p. 1

The importance of physical activity throughout an individual's life is indisputable. As healthcare

Gait. Kinesiology RHS 341 Lecture 12 Dr. Einas Al-Eisa

Palacký Univerzity in Olomouc Faculty of Physical Culture

The Examination of Upper Limb Ambidexterity in Wrestling Snap Down Technique

Inter and intra-limb processes of gait control.

Ankle biomechanics demonstrates excessive and prolonged time to peak rearfoot eversion (see Foot Complex graph). We would not necessarily expect

The technique of reciprocal walking using the hip guidance orthosis (hgo) with crutches

A Pilot Study on Electromyographic Analysis of Single and Double Revolution Jumps in Figure Skating

Walking speemtmmkubjects and amputees: aspects of validity of gait analysis

The Effect of a Seven Week Exercise Program on Golf Swing Performance and Musculoskeletal Screening Scores

A QUALITATIVE ANALYSIS OF THE HIGH RACQUET POSITION BACKHAND DRIVE OF AN ELITE RACQUETBALL PLAYER

Proposed Paralympic Classification System for Va a Information for National federations and National Paralympic Committees

Analysis of Gait Characteristics Changes in Normal Walking and Fast Walking Of the Elderly People

REPRESENTATION OF HUMAN GAIT TRAJECTORY THROUGH TEMPOROSPATIAL IMAGE MODELLING

The Scientific Bulletin of VALAHIA University MATERIALS and MECHANICS Nr. 5 (year 8) 2010

EXSC 408L Fall '03 Problem Set #2 Linear Motion. Linear Motion

Inertial compensation for belt acceleration in an instrumented treadmill

REPLACING REDUNDANT STABILOMETRY PARAMETERS WITH RATIO AND MAXIMUM DEVIATION PARAMETERS

A New Approach to Modeling Vertical Stiffness in Heel-Toe Distance Runners

A Bio-inspired Behavior Based Bipedal Locomotion Control B4LC Method for Bipedal Upslope Walking

ITF Coaches Education Programme Biomechanics and stroke production: implications for the tennis coach

Motion Control of a Bipedal Walking Robot

SAMPLE MAT Proceedings of the 10th International Conference on Stability of Ships

Evaluation of gait symmetry in poliomyelitis subjects : Comparison of a

Development of a Simulation Model for Swimming with Diving Fins

video Outline Pre-requisites of Typical Gait Case Studies Case 1 L5 Myelomeningocele Case 1 L5 Myelomeningocele

Supplementary Figure 1 An insect model based on Drosophila melanogaster. (a)

REPORT. A comparative study of the mechanical and biomechanical behaviour of natural turf and hybrid turf for the practise of sports

Saturday, 15 July 2006 SAP-30: 10:45-11:15 CHANGE OF SPEED IN SIMULATED CROSS-COUNTRY SKI RACING: A KINEMATIC ANALYSIS

RELATIONSHIP OF SELECTED KINEMATIC VARIABLES WITH THE PERFORMANCE OF DOUBLE HANDEDBACKHAND IN TENNIS. Rajesh Kumar, M.P.Ed,

Running injuries - what are the most important factors

PURPOSE. METHODS Design

UvA-DARE (Digital Academic Repository) Hip and groin pain in athletes Tak, I.J.R. Link to publication

Can Asymmetric Running Patterns Be Predicted By Assessment of Asymmetric Standing Posture? A Case Study in Elite College Runners

Gait pattern and spinal movement in walking - A therapeutic approach in juvenile scoliosis

First Experimental investigations on Wheel- Walking for improving Triple-Bogie rover locomotion performances

Denny Wells, Jacqueline Alderson, Kane Middleton and Cyril Donnelly

VISUOMOTOR CONTROL OF STRAIGHT AND BREAKING GOLF PUTTS 1

ITF Coaches Education Programme Biomechanics of the forehand stroke

Transcription:

1 2 3 4 5 Original Article: Analysis of gait data in manifold space: The control of the centre of mass during gait. Enrica Papi a,1 ; Philip J. Rowe a and Valerie M. Pomeroy b 6 7 8 9 10 11 a Department of Bioengineering, University of Strathclyde, Wolfson Building, 106 Rottenrow, Glasgow, G4 0NW, UK b School of Allied Health Professions, Faculty of Medicine and Health Sciences, Queen s Building, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, UK 1 Department of Surgery and Cancer, Imperial College London, Room 7L16, Floor 7, Laboratory Block, Charing Cross Hospital, London, W6 8RF, UK 12 13 14 15 16 17 18 19 20 21 22 Corresponding author Enrica Papi Department of Surgery and Cancer, Imperial College London, Room 7L16, Floor 7, Laboratory Block, Charing Cross Hospital, London, W6 8RF, UK Phone: +44(0)20 3313 8833 Fax : +44(0)20 8383 8835 Email: e.papi@imperial.ac.uk 23 24 Keywords: Gait, Motor Control, Uncontrolled Manifold, Stroke, Centre Of Mass. 25 26 Word Count 3588 27 28 29 1

30 Abstract 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 This study investigated the feasibility of the uncontrolled manifold approach (UCM) to analyse gait data variability in relation to the control of the centre of mass (COM) in adults with and without neuropathology. The proposed method was applied to six able-bodied subjects to characterise mechanisms of normal postural control during stance phase. This approach was repeated on an early stroke patient, who attended the laboratory three times at three monthly intervals, to characterise the variability of COM movement during walking with and without an orthosis. Both able-bodied subjects and the stroke participant controlled COM movement during stance but utilized a different combination of lower limb joint kinematics to ensure that the COM trajectory was not compromised. Interestingly, the stroke subject, despite a higher variability in joint kinematics, was able to maintain a stable COM position throughout stance phase. The stabilisation of the COM decreased when the patient walked unaided without the prescribed orthosis but increased over the six months of study. The UCM analysis demonstrated how a stroke patient used a range of lower limb motion pattern to control the COM. It is suggested that this analysis can be used to track changes in these movement patterns in response to rehabilitation. As such we propose that this approach could have clinical utility to evaluate and prescribe rehabilitation in stroke patients. 46 47 48 49 50 51 52 53 54 55 56 57 58 59 2

60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 1. Introduction Many stroke survivors present with an altered ability to walk. Compensatory actions and gait strategies are often adopted to achieve a safe walking activity. Motor control is a key issue for those people who have suffered from a stroke and for whom limited joint coordination impairs their mobility. An understanding of how the central nervous system (CNS) compensates to control motion following a stroke may inform subsequent therapy. The theory of the uncontrolled manifold (UCM) has been recently introduced (Latash et al., 2007; Scholz and Sch ner, 1999) to investigate how the CNS acts with respect to selected motor tasks by choosing combinations of different musculoskeletal elements that are involved in the performance of the task. That is to say CNS may employ a variety of different approaches to achieve a task. Exploiting this approach it may be possible to predict which motor variables the CNS controls and what are the elements/degree of freedoms (DOFs) that it has to organise for that particular motor task to be performed. This theory can thus be seen as an analysis of the variability of a selected functional task in a multi-degree of freedom system. The variability can either be good, if the task goal remains unaltered, or bad, if deviations from it occur. The UCM itself is a subspace of all possible combinations of motor elements (elemental variables) that lead to a consistent value of a performance variable. For example all the different combinations of lower limb joint angles that together place the centre of mass (COM) in a certain position in 3D space define a UCM subspace. It is defined uncontrolled because the control of the variability within it is unnecessary as all the combinations (i.e., set of lower limb joint angles) within that subspace preserve the performance variable value (i.e., 3D position of the COM) (Scholz and Sch ner, 1999). Thereby, the UCM approach can also be seen as a method to quantify synergies. In this context, a synergy refers to an organization of elemental variables that stabilises a performance variable (Latash and Anson, 2006). A practical example could be to use the UCM approach to understand how the CNS organises joint angles (elemental variables) to allow a smooth COM movement (performance variable) and thus safe locomotion. It is important to mention that variability across trials is partitioned into two components: one that lies within the UCM and one that is perpendicular to the UCM. These two variabilities, expressed as indexes of variance across repetitions of the same task, are used to verify the hypothesis about the aspects of movement that are controlled. If the variance within the UCM is bigger than the one perpendicular to it, the hypothesis about the stabilisation of the selected motor task is accepted. This analysis can provide clinicians with a better understanding of motor coordination and its relationship with rehabilitation approaches providing an explanation on how different dynamic resources can lead to a successful motor performance. Having information on the behaviour of the system will allow a more specific and 3

90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 individualised treatment to accelerate recovery as the intervention target (musculoskeletal elements) can be identified. Which movement variations should be encouraged and which discouraged? An answer to this question will advance clinical practice and outcomes for stroke survivors. The UCM analysis method has recently been used to verify the control of motor task predominately related to the upper extremity, sit-to-stand, standing and hopping performances of able-bodied and impaired subjects (Auyang et al., 2009; Domkin et al., 2002; Freitas et al., 2006; Hsu et al., 2007; Reisman and Scholz, 2003; Scholz et al., 2003; Scholz and Sch ner, 1999; Yang et al., 2007; Yen and Chang, 2010). Less consideration has been given to gait and the relative motor redundancy it contains. Among the many studies that looked into gait and centre of mass stabilisation, only one used the UCM approach (Black et al., 2007) and one used covariation analysis, a method comparable to the UCM (Verrel et al., 2010). Both studies however limited their analysis to walking at times of heel strikes rather than the entire time history of the gait cycle. The ability to analyse the whole cycle rather than a single instant would have greater clinical impact in rehabilitation of impaired gait. Of the few studies (Krishnan et al., 2013; Robert et al., 2009; Rosenblatt et al., 2014) that analysed the temporal evolution of the UCM approach throughout the gait cycle, none have considered the COM trajectory as a performance variable. The exploratory study reported here, therefore, investigated the potential usefulness of UCM analysis of postural control during the stance phase of walking. It is known that stabilisation of the COM is key to walking ability so that is the component of postural control identified as the key focus. We hypothesised that different combinations of lower limb joint angles (kinematic synergy) can be used to control the COM movement while walking. The specific aims of this study were to (a) determine the feasibility of undertaking an UCM analysis of control of the COM during the stance phase of gait, (b) to find if such analysis could provide more knowledge of COM control than standard biomechanical analysis techniques, (c) explore whether the UCM analysis could identify differences between a stroke survivor with walking difficulty and adults without a brain lesion and the stroke survivor walking with and without a custom-made ankle-foot orthosis (AFO). 114 115 116 117 118 119 2. Methods 2.1 Participants Six adults (3 female, 3 male; height: 168.9 (± 10.5) cm, mass: 68.2 (± 9.9) kg, age: 29.8 (± 6.7) years ) with no known neurological pathology participated in the study. In addition, one 81 years old male (80 kg, 180 cm) was recruited 2 months after experiencing a stroke. He presented a left side hemiplegia of the upper and lower body. 4

120 121 122 He was prescribed a 5mm polypropylene AFO with carbon fibre reinforcement at the malleoli level. The AFO and shoes combination were tuned at 10 of forward inclination. All participants provided written consent for the study which was approved by the local ethics committee (West of Scotland REC3). 123 124 125 126 127 128 129 130 131 132 2.2 Equipment and Experimental Procedure A twelve-camera motion capture system (Vicon, Oxford Metrics Ltd., UK) was used to collect experimental data at 100 Hz while participants walked at comfortable speed on a flat surface of 6 m in length. The gait analysis protocol developed within the Bioengineering Department at University of Strathclyde was followed for data collection and processing (Papi et al., 2011). Ten walking trials were recorded with able-bodied subjects and the data derived from their left leg were used in the subsequent analysis. The stroke participant was assessed three times at three monthly time points. Six trials were collected during walking with and without AFO at each visit. Data from the hemiplegic leg were considered. 133 134 135 136 137 2.3 Data Processing Initial data processing was performed using Nexus software (Oxford Metrics Ltd., UK). Hip, knee and ankle sagittal angles, and the 3-D coordinates of anatomical landmarks were output. Data were time normalised to 100% of stance phase. 138 139 140 141 142 143 144 145 146 147 148 149 2.4 UCM Formulation The UCM method was applied to characterise the control of the COM during stance phase. The performance variable was COM movement and the elemental variables were the lower limb joint rotations. As a preliminary development of the UCM method for gait, it was decided to conduct the analysis for stance phase in the sagittal plane and to approximate the COM as a fixed point in the pelvis. This point was defined by the intersection of the diagonals connecting the anterior and posterior superior iliac spines. To estimate how the variability of joint angles influences the position of the COM in a global sagittal plane with x antero/posterior and y vertical axis, a geometric model (Fig. 1) that links hip, knee and ankle rotations to the COM position throughout the gait cycle was defined. Since the position of the COM depends also on the position of the foot on the ground and in particular on the angle between the sole of the foot and the ground, this angle was also considered an elemental variable: 5

150 151 (1) 152 153 Where: is the angle between the sole of the foot and the ground,, are the ankle, knee and hip 154 155 156 sagittal angles respectively. representing the ground: is the angle (Fig. 2) between the vectors, characterising the sole of the foot, and (2) 157 158 159 160 161 162 The first modelling requirement was to consider the position of the foot on the ground during walking to define the ankle joint centre. Three main cases were identified (Fig. 2): 1- the heel is in contact with the ground:, 2- foot flat:, 3- the fore part of the foot is on the ground:. The sagittal position of the ankle joint centre was defined as follow, for case 1 and 2: 163 164 165 (3) (4) 166 167 168 169 170 171 172 173 174 175 Where: α is the angle at the rear of the foot (Fig. 2) and centre (AJC) and the calcaneus. For case 3: Where: β is the angle at the fore part of the foot (Fig 2) and is the length of the segment joining the ankle joint (5) (6) is the length of the segment joining the 176 177 midpoint between the 1 st and 5 th metatarsal head and the AJC. Through a trigonometric analysis of the leg segment, the COM position in the sagittal plane can be expressed as: 178 179 (7) 6

180 181 (8) 182 183 184 185 186 Where: AK is the shank segment length, KH is the thigh segment length, HCM is the hip centre to COM segment length. The next step required for this approach was the linearization of the UCM (Latash et al., 2007). This was necessary because the concept of variance is a linear concept while, the UCM, and in particular the geometric 187 model defined, are not linear. The linearization implies the definition of the Jacobian matrix,, and the 188 189 190 computation of its null space, N(J). The Jacobian matrix is a matrix of all first-order partial derivatives of the COM coordinates with respect to the elemental variables. Changes in joint angles and changes of the COM trajectory are linked through this matrix. 191 The null space (Eq. 9) of the Jacobian matrix, spanned by the basis vectors, is the linear subspace of all 192 193 194 joint angles combinations that leave the COM coordinates unaffected. The dimension of this subspace is (n d) where n is the number of elemental variables and d is the number of dimensions of the performance variable. The null space in the current case had a dimensionality of 2. 195 196 (9) 197 198 199 The linearization was performed around a reference configuration, defined as the mean joint configuration across trials for each instant analysed (Latash et al., 2007). This is assumed to represent the set of angles that 200 leads to the desired COM position. The Jacobian matrix, was calculated with respect to this configuration. 201 202 The computation of the Jacobian matrix and then of its null space was performed for each time point of stance phase of each trial and hence they continuously varied. 203 204 In addition, the deviation, for each trial, from the mean joint configuration ( each instant: ) was calculated at 205 (10) 206 7

207 208 The obtained deviation vector (DV) was decomposed into a component that is within ( ) and perpendicular ( ) to the null space: 209 210 (11) 211 212 213 214 The scalar values obtained represent to what extent the trial joint configuration is consistent with the reference joint configuration. (12) 215 The variances of these projections (, ) were then calculated. Since they have different dimensions, the 216 217 variances were normalised per degree of freedom of each subspace. The variance across trials within the linearized UCM ( ) and perpendicular to it ( ) are: 218 219 (13) 220 221 (14) 222 223 224 225 226 227 228 Where: and are the squared length of the deviation vector component lying within and perpendicular the UCM respectively, N is the number of trials, n are the elemental variables and d is the dimensionality of the performance variable. The variances within and perpendicular the UCM were compared to verify the initial hypothesis about the control of the COM trajectory during walking. A ratio defined as in (Eq. 15) was used to summarise the results: 229 230 (15) 231 232 233 234 235 2.5 Data Analysis Matlab (The MathWorks Inc., Massachusetts, US) was used for data processing. Data analysis was conducted by first analysing the variability of joint angles and COM trajectories and then, by evaluating the structure of variances within the UCM framework. 8

236 237 238 239 240 241 For the analysis, stance phase was divided into three phases (i) contact phase (20%, initial contact and loading response), (ii) mid stance (30%), (iii) propulsive phase (50%, Terminal stance and pre swing) (Perry et al., ) to allow a more detailed analysis of the evolution of variances components during stance. To test the hypothesis that a kinematic synergy of the joint angles stabilise the COM trajectory and to evaluate if the strength of the synergy change within stance phase, a repeated measures ANOVA was performed with factors of variance components ( and ) and phases (contact phase, mid stance, propulsive phase). 242 One-tailed t-test was applied to verify if the ratio was significantly greater than 0, implying, and hence 243 244 245 246 accepting the hypothesis about the control of the COM movement during gait. Two-sample t-test was used for comparisons between stroke patient and able-bodied subjects parameters. Significance was set at 0.05. Prior to statistical analysis, variances components and ratios were averaged across each phase of the stance phase. 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 3. Results Overall sagittal joint kinematic mean pattern and variability (standard deviations bars) throughout stance phase are shown in Figure 3 for each able-bodied subject and in Figure 4 for the hemiplegic leg of the stroke participant during 3 assessments for walking with AFO and without AFO. Higher variability in joint kinematic is observed for the stroke participant and in particular during baseline test. COM displacements in the x and y directions are illustrated in Figure 5 and 6 for able-bodied participants and stroke patient respectively. The latter showed higher variability in COM trajectories. The variability in joint kinematics was related to the variability of the COM position through the UCM analysis. Time series of the variance within and perpendicular to the linearized UCM are presented in the graphs on the right side and ratios on the left side of Figure 7 and 8 for able-bodied participants and stroke patient respectively. Ratio values above or equal to 0 indicate that the hypothesis about the control of the COM can be accepted. Ratios for able-bodied participants and the stroke patient at the three assessments were all statistically different from 0 (p<0.000) and no statistical difference existed between stroke and able-bodied participant ratios (p=0.216). 262 The results of the ANOVA indicated a significant main effect of variance component ( ) supporting the 263 264 265 hypothesis that a kinematic synergy existed to stabilise the COM (F 1, 11 =18.5; p=0.001) but it showed that the synergy does not change within different periods of the stance phase (F 2, 22 =1.1; p=0.330). The stroke patient despite showing variability in joint angles higher without the AFO (Fig. 4) and in comparison to able-bodied 9

266 adults, maintained a stabilised COM position during walking, with and without the AFO. was statistically 267 different between stroke and able-bodied participants, although did not reach significance. 268 269 270 271 272 273 274 The control of COM when the patient used the AFO was reduced over time (from first to last visit); on the contrary, the control imposed without the AFO was always high. Confidence and stability acquired when walking with the AFO allowed the patient to be less vigilant with regards to COM displacement. For the patient a progression towards the ratio waveform seen in able-bodied subjects can be observed particularly at the 6-months follow-up, although differences were present due to altered gait events timing in hemiplegic walking (i.e., Prolonged terminal stance in which the hemiplegic leg waited for the contralateral foot to completely strike the ground before going into swing phase). 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 4. Discussion 4.1 Interpretation of UCM method findings The UCM hypothesis proposes that variability in movement patterns is essential for the performance of everyday functional tasks in a variety of environmental contexts. In other words, what was termed DOF problem might actually be a DOF advantage especially in the presence of a stroke lesion in the brain. Availability of different movement patterns to achieve the same functional goal enables achievement of that goal albeit with what is consider sub-optimal movement. Utilising the UCM approach, in this study, the variability obtained in the COM displacements in the sagittal plane during stance phase depending on sagittal kinematics was classified accordingly to the subject s ability to control the position of the COM. Lower limb sagittal kinematics showed variability through stance phase that was more evident in the stroke participant and particularly for walking without an AFO. All participants had, overall, a variance within the UCM greater than its complement in all phases of stance phase and hence maintained the COM position despite the stroke patient presenting with a higher variability of joint kinematics. This indicated that the patient adopted a way of walking which utilised more variation in joint configurations (Fig. 4) without altering the COM position. A tendency to increase the variance within the linearized UCM was observed when the patient walked without the AFO, indicating a greater focus was placed on the mainteinance of the COM position than when walking with an AFO. The AFO appeared to give the patient more confidence such that he was able to lower the control imposed on the COM position and approach a more normal gait. 10

294 295 296 297 298 299 300 301 302 303 304 The method highlights how adults with and without brain lesion used differently the available DOFs at the joint level but all achieving a stable COM movement during stance of the gait and how this can change by the use of AFO and with time. The aim of applying the UCM approach was to introduce a method, that has been applied to more stationary tasks, to gait in order to improve our understanding of how able-bodied and impaired populations control walking and thus providing additional explanation to standard kinematics and trajectories curves. This approach provides explanation of the variability in joint kinematics and how these are functionally employed with the purpose of achieving a successful task performance by controlling a particular parameter, in this case the trajectory of the COM. In rehabilitation practice it could guide clinicians as to how intervene and verify if selected interventions hinder or improve the achievement of locomotion by identifying which movement patterns should be encouraged to achieve a succesfull task performance. 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 4.2 Methodological Considerations The UCM approach was formulated to test the hypothesis of the control of COM position during walking. The advantages of the method proposed over studies already conducted and looking into COM trajectory are the ability to perform the analysis through a period of time rather than only instances of the motor task performed and its applicability to gait whereas most of the previous studies concentrated on upper body performances with only few in lower limbs tasks and gait (Auyang et al., 2009; Black et al., 2007; Domkin et al., 2002; Hsu et al., 2007; Krishnan et al., 2013; Reisman and Scholz, 2003; Robert et al., 2009; Rosenblatt et al., 2014; Scholz et al., 2003; Scholz and Sch ner, 1999; Yang et al., 2007; Yen and Chang, 2010). The method introduced showed to be feasible as reasonable and interpretable results were obtained although not without limitations. The geometric model used assumed the COM as a fixed point in the pelvis without taking into account the contributions of the arm and trunk movement, the analysis was confined to the sagittal plane and to the duration of stance phase. These choices were made to simplify the formulations of the UCM approach in the first instance. Calculation of the true COM as a sum of each body segment centre of mass weighted with respect to segments mass could be included as further development of this approach as well as the COM 3D coordinates. The analysis in swing phase is more complex as no fixed support is available for the swinging leg. A further development could be increasing the complexity of the model adopted by accounting for DOFs that were not considered in the current study. For example considering both legs will give a better description of the 11

323 324 3D COM movement and will overcome the lack of a fixed base in swing as the analysis will switch to the contralateral leg. 325 326 327 328 329 330 331 332 333 334 335 336 5. Conclusion The feasibility of the application of the UCM method to gait data time series was shown in both adults with and without neuropathology. Further research on a larger population group is required to strengthen the clinical meaning of the findings. However, this approach shows great promise with respect to understanding postural control mechanism during walking and the relation to rehab approaches. The UCM analysis allows a classification of the variability of a selected control variable with respect to the elemental variables and it can be seen as a tool to appraise the effect of rehabilitative techniques and rehabilitation over time. Its utility is to provide an explanation of how the CNS acts to cope with the different DOFs available for the walking task and how it compensates if impairments at the elemental variable level are present. It is another form of analysis to be added to the results of 3-D gait analysis in the process of evaluating stroke gait kinematics for which a retrospective analysis is conducted through the UCM. 337 338 339 Conflict of interest We declare that we have no conflict of interest 340 341 342 Acknowledgement Nonthing to declare 343 344 345 346 347 348 349 350 351 References Auyang, A.G., Yen, J.T., Chang, Y.H., 2009. Neuromechanical stabilization of leg length and orientation through interjoint compensation during human hopping. Experimental Brain Research 192, 253-264. Black, D.P., Smith, B.A., Wu, J., Ulrich, B.D., 2007. Uncontrolled manifold analysis of segmental angle variability during walking: preadolescents with and without Down syndrome. Experimental Brain Research 183, 511-521. Domkin, D., Laczko, J., Jaric, S., Johansson, H., Latash, M.L.2002. Structure of joint variability in bimanual pointing tasks. Experimental Brain Research 143, 11-23. 12

352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 Freitas, S.M.S.F., Duarte, M., Latash M.L., 2006. Two kinematic synergies in voluntary whole-body movements during standing. Journal of Neurophysiology 95,636-645. Hsu, W.L., Scholz, J.P., Schöner, G., Jeka, J.J., Kiemel, T., 2007. Control and estimation of posture during quiet stance depends on multijoint coordination. Journal of Neurophysiology 97, 3024-3035. Krishnan, V., Rosenblatt, N.J,. Latash, M.L., Grabiner, M.D., 2013. The effects of age on stabilization of the mediolateral trajectory of the swing foot. Gait and Posture 38, 923-928. Latash, M.L., Anson, J.G., 2006. Synergies in health and disease: relations to adaptive changes in motor coordination. Physical Therapy 86,1151-1160. Latash, M.L., Scholz, J.P., Schöner, G., 2007. Toward a new theory of motor synergies. Motor Control 11, 276-308. Perry, J., 1992. Gait Analysis: Normal and pathological Function. SLACK Incorporated, Thorofare, New Jersey. Papi, E., Ugbolue, U.C., Solomonidis, S.E., Rowe, P.J., 2011. Development of a gait analysis protocol for the evaluation of stroke patients. In Proceedings of 23 rd International Society of Biomechanics Congress, Brussels. Reisman, D., Scholz, J.P., 2003. Aspects of joint coordination are preserved during pointing in persons with post-stroke hemiparesis. Brain 126, 2510-2527. Robert, T., Bennett, B.C., Russell, S.D., Zirker, C.A,. Abel, M.F., 2009. Angular momentum synergies during walking. Experimental Brain Research 197, 185-197. Rosenblatt N.J., Hurt C.P., Latash M.L., Grabiner M.D., 2014).An apparent contradiction: Increasing variability to achieve greater precision? Experimental Brain Research 232, 403-413. Scholz, J.P., Kang, N., Patterson, D., Latash, M.L., 2003. Uncontrolled manifold analysis of single trials during multi-finger force production by persons with and without Down syndrome. Experimental Brain Research 153, 45-58. Scholz, J.P., Sch ner, G., 1999. The uncontrolled manifold concept: Identifying control variables for a functional task. Experimental Brain Research 126, 289-306. Verrel, J., Lovden, M., Linderberger, U., 2010. Motor-equivalent covariation stabilizes step parameters and center of mass position during treadmill walking. Experimental Brain Research 207, 13-26. Yang, J.F., Scholz, J.P., Latash, M.L. 2007. The role of kinematic redundancy in adaptation of reaching. Experimental Brain Research 176, 54-69. Yen, J.T., Chang, Y.H., 2010. Rate-dependent control strategies stabilize limb forces during human locomotion. Journal of the Royal Society Interface; 7, 801-810. 13

382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 Figures Captions: Fig. 1. Leg and foot stick model. Fig. 2. Positions of the foot on the ground at the three identified key points depending on θ G value. Fig.3. Mean and standard deviations shade areas of sagittal hip, knee and ankle kinematics of six able-bodied subject over 10 walking trials during stance phase. Fig. 4. Mean and standard deviations shade areas over 6 gait cycles of sagittal hip, knee and ankle joint kinematics with (dashed lines) and without AFO (solid lines) for the stroke patient during the three assessments. Fig. 5. Mean COM displacement during stance in antero/posterior (x) and vertical (y) directions for the six ablebodied subjects. Standard deviation shade areas are shown. Fig. 6. Mean COM displacement during stance in x and y directions for the stroke patient during walking with (dashed lines) and without AFO (solid lines). Standard deviation shade areas are shown. Fig. 7. Variance components within (Vucm, solid lines) and perpendicular (Vort, dashed lines) to the linearized UCM for able-bodied subjects (left). Mean ratio (± standard deviation) across six able-bodied subjects (ratio). Fig. 8. Variance components within (Vucm, solid lines) and perpendicular (Vort, dashed lines) to the linearized UCM and mean ratio for the stroke patient at three assessments with (solid lines) and without AFO (dashed lines). 14