BioMed Research International
Volume 2015 (2015), Article ID 402596, 11 pages
http://dx.doi.org/10.1155/2015/402596

Research Article

The Discriminant Value of Phase-Dependent Local Dynamic Stability of Daily Life Walking in Older Adult Community-Dwelling Fallers and Nonfallers

Espen A. F. Ihlen,¹ Aner Weiss,² Jorunn L. Helbostad,1,3 and Jeffrey M. Hausdorff2,4

¹Department of Neuroscience, Norwegian University of Science and Technology, 7491 Trondheim, Norway
²Center for the Study of Movement, Cognition, and Mobility, Department of Neurology, Tel Aviv Sourasky Medical Center, 64239 Tel Aviv, Israel
³Clinic for Clinical Services, St. Olavs Hospital, Trondheim University Hospital, 7006 Trondheim, Norway
⁴Department of Physical Therapy, Sackler School of Medicine and Sagol School of Neuroscience, Tel Aviv University, Tel Aviv, Israel

Received 2 June 2015; Accepted 9 September 2015

Academic Editor: Catherine Said

Copyright © 2015 Espen A. F. Ihlen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The present study compares phase-dependent measures of local dynamic stability of daily life walking with 35 conventional gait features in their ability to discriminate between community-dwelling older fallers and nonfallers. The study reanalyzes 3D-acceleration data of 3-day daily life activity from 39 older people who reported less than 2 falls during one year and 31 who reported two or more falls. Phase-dependent local dynamic stability was defined for initial perturbation at 0%, 20%, 40%, 60%, and 80% of the step cycle. A partial least square discriminant analysis (PLS-DA) was used to compare the discriminant abilities of phase-dependent local dynamic stability with the discriminant abilities of 35 conventional gait features. The phase-dependent local dynamic stability λ at 0% and 60% of the step cycle discriminated well between fallers and nonfallers (AUC = 0.83) and was significantly larger () for the nonfallers. Furthermore, phase-dependent λ discriminated as well between fallers and nonfallers as all other gait features combined. The present result suggests that phase-dependent measures of local dynamic stability of daily life walking might be of importance for further development in early fall risk screening tools.

1. Introduction

Falls among older people are an important reason for dependence in daily life, reduced quality of life, and admission to hospitals or nursing homes. At a European level, the annual costs of falls among older persons are estimated to be 30 Billion Euro [1]. Early prediction of falls amongst community-dwelling older persons could provide opportunities for early fall prevention. Thus, considerable efforts have been made for early fall risk assessment and fall prediction in older persons.

More than 400 risk factors for falls have been reported (e.g., [2]). Most risk factors have been assessed in laboratory settings or in clinical test situations, and fall risk assessment tools have been developed based on these assessments [3–6]. However, most of these screening tools reflect the performance of the older person at a specific moment in time or they are based on self-report. Furthermore, falls in older persons are often experienced during activities of daily living [7, 8]. Thus, monitoring of behaviour in daily life, rather than assessment by performance tests and self-report, may be important for furthering evidence-based recommendations for fall risk assessment and screening for fall prevention interventions.

Daily life activities, like lying, sitting, standing, and walking, can be identified by body-fixed sensors containing inertial sensors like accelerometers and gyroscopes [9]. Features of the acceleration signal within these activities and in the transition between activities might be important in fall risk assessment [10]. Several measures of gait stability and variability are significantly different between elderly fallers and nonfallers [11–13]. Amongst these is the measure of local dynamic stability that has been suggested to be one of the most sensitive measures of gait instability in older persons [14]. Local dynamic stability is defined as the rate of exponential increase in infinitely close distances between trajectories in the reconstructed state space of the gait dynamics [15]. These distances are considered as infinitesimal perturbations and, thus, local dynamic stability defines the reaction of the gait dynamics to these perturbations. The gait dynamics are local dynamic stable when , indicating an exponential decrease in the distance between neighbouring trajectories. In contrast, the gait dynamics are local dynamic unstable when , indicating an exponential increase in the distance between neighbouring trajectories. Recent extensions of computational methods for local dynamic stability λ indicate that λ is phase dependent and changes within the gait cycle [16, 17]. Despite promising results, the discriminating ability of phase-dependent has not been compared to the discriminating abilities of other features of daily life walking.

The main aim of the present study is to compare phase-dependent λ with conventional gait features in their ability to discriminate between the daily life walking of community-dwelling elderly fallers and nonfallers.

2. Methods

2.1. Participants and Data Collection

Inertial sensor data previously studied by Weiss et al. [13] were reanalysed in the present study. The data can be downloaded at http://www.physionet.org/. The data consist of 3 days of 3D-acceleration data from 71 community-dwelling older persons (mean age:  yrs; range: 65–87 yrs; gender: 64.79% women; mean height:  m; mean weight:  kg). None of the participants included had been diagnosed with gait or balance disorders or had cognitive impairments (i.e., Mini Mental State Examination score > 24). The participants were classified as fallers or nonfallers based on retrospective self-report. Participants reporting 2 or more falls in the year prior to testing were considered as fallers; this definition was used to ensure a clear distinction between the two groups and to focus on (multiple) fallers and nonfallers, excluding older adults who may be in an intermediate, less well-defined, and more ambiguous state with respect to their fall history. There was no difference between fallers and nonfallers in age, gender, years of education, height, weight, or body mass index, but a difference in in-lab preferred gait speed (nonfallers:  m/s; fallers:  m/s). The acceleration along the anterior-posterior (AP), mediolateral (ML), and vertical (V) axes was sampled at 100 Hz by a small inertial sensor (DynaPort Hybrid, McRoberts, The Hague, Netherlands; 87 × 45 × 14 mm, 74 g). The sensor had a range and resolution of ±6 g and ±1 mg, respectively. The acceleration signals were recorded on a Secure Digital (SD) card at a sample frequency of 100 Hz and later transferred to a personal computer for further analysis using Matlab (MathWorks, Natick, MA). The sensor was fitted without any difficulties on a belt on the center of lower back, at the L5 level. The sensor had to be removed during the shower and swimming and occasionally during sleep. The participants received a diary for tracking when and why they took off and put on the device. No specific problems were evident during data collection and retrieval.

2.2. Preprocessing of the Data

The classification procedure was restricted to the walking bouts of duration ≥ 60 seconds, identical to those originally analysed by Weiss et al. [13]. This size was chosen to ensure that these were indeed walking segments and that the acceleration derived measures would be robust. The walking bouts were identified by using two filters: one filter was based on the acceleration-magnitude, and the other filter was based on the energy in the frequency domain [13, 18]. The activity bouts were visually observed to ensure that these were indeed valid walking segments. A mean of 28.3 walking bouts (range: 5 to 90) with duration ≥ 60 seconds was identified for each of the participants. There was no significant difference in number of walking bouts between fallers and nonfallers. The reader is referred to Weiss et al. [13] for further details about the participants, protocols, and preprocessing of the 3D-acceleration data.

Intrastep 3D-velocity was estimated from the 3D-acceleration signal. The 3D-acceleration was detrended using an orthogonal wavelet procedure that preserved intrastep variation in the 3D-velocity but removed interstride nonlinear trends [19]. This detrending procedure provides stationary 3D-velocity signal necessary for computation of local dynamic stability [20]. The local maxima of the vertical velocity were defined as the beginning of a step. This step identification method provided similar results to the autocorrelation method used in previous studies based on comparison of the mean step time [12, 21].

2.3. State Space Construction Methods

Two 6D state spaces were constructed for each walking bout by the two following methods [20].

Method 1. Differential coordinate embedding was defined as , combining both acceleration signal and velocity signal in AP, ML, and V directions. The local dynamic stability computed from state space construction Method 1 has the subscript in the result section.

Method 2. Delayed coordinate embedding was defined as , where is the velocity signal,  s is the sampling interval, and is the time lag. This delayed coordinate embedding combines the velocity signal in AP, ML, and V directions for the velocity signal and uses a short time lag, , to prevent the blending of phases within the gait cycle. The local dynamic stability computed from state space construction Method 2 has the subscript in the result section.

2.4. Computation of Phase-Dependent Local Dynamic Stability λ

Phase-dependent local dynamic stability was defined according to a method developed by Ihlen et al. [16] and based on two equations:

where is the reaction curve of the initial perturbation and the outer brackets are the mean across all steps in the walking bout. The initial perturbation was defined as the distance between the reference point and the th neighbourhood trajectory within a small neighbourhood of predefined size (see Figure 1(a)). The initial perturbation was considered at 0%, 20%, 40%, 60%, and 80% of the step cycle (see Figures 1(b) and 1(c), e.g., for 0% and 60%). The reaction distance of the initial distance was traced to the next starting point of a step. The average reaction distance was computed across all th neighbourhood trajectories. The reaction distance, , for less than 10 neighbourhood trajectories or with instantaneous stride time outside the 5% and 95% percentiles was excluded for further analysis. The portion of excluded was less than 10% of the total number of strides for all participants and these strides were typically short periods of deviating patterns of the acceleration and velocity signals due to large deviations from normal patterns. In (1), the remaining was normalized to step time before λ was assessed as the linear regression slope for the first 10% of the step cycle (see Figure 1(b)). In (2), the remaining was not normalized and was the reaction distance at time equal to 10% of the step cycle. The first 10% of the step cycle was considered to prevent influence of curvatures in state space trajectories [22]. The median of λ was computed across all walking bouts for each participant. In addition, the conventional phase-independent was computed by the method of Wolf et al. [23]. Wolf’s method was applied to the acceleration signal in the AP, ML, and V directions, separately. A 6D delayed coordinate embedding was used with time lag, , which was the mean lag for the first minima of the average mutual information function [24]. The Matlab code for the phase-dependent measures, and , is available at http://www.physionet.org/.

(a)

(b)

(c)

Figure 1: (a) A schematic representation of the reaction distance based on trajectories within a small neighborhood (gray circle) of the state space. The initial average perturbation distance was computed from multiple distances (t) within the neighborhood. Note that the left panel illustrates a 3D state space reconstruction where the computations of are based on a 6D state space reconstruction. (b) A representative example of a log-reaction curve normalized to the step cycle for a faller (red) and a nonfaller (blue) for initial perturbation at 0% of the step cycle. (c) The same example of a log-reaction curve for initial perturbation at 60% of the step cycle. The slopes of the regression lines for the initial 10% (black lines in shaded areas of (a) and (b)) were defined as the local dynamic stability according to (1).

2.5. Test-Retest Reliability of Phase-Dependent Local Dynamic Stability λ

Different circumstances of walking, like turns, walking surfaces, obstacles, variations in walking speed, and dual tasking, may introduce random fluctuations in phase-dependent local dynamic stability between walking bouts. Even though the median of each selected feature across several walking bouts will reduce these fluctuations, it is uncertain if the median has sufficient reliability across walking bouts within a 3-day period. Thus, the test-retest reliability of the median of local dynamic stability was assessed by interclass correlation (ICC) absolute agreement for the first and the last 1/3 of the walking bouts within the 3-day recording period.

2.6. Partial Linear Square Discriminant Analysis (PLS-DA)

PLS-DA relates the predictor matrix of gait features with the response of fall status (see Figure 2). PLS-DA is able to identify a low-dimensional latent structure from a large number of gait features which discriminates between fallers and nonfallers. In contrast to other regression approaches, PLS-DA is designed to perform discriminate analyses based on a large set of noisy and collinear predictors and is therefore suitable for the large number of gait features investigated in the present study. The present study used a nonlinear iterative partial least square (NIPALS) algorithm extended by a target projection (TP) as summarized in Figure 2 [25–27]. The TP-loadings define the contribution of each gait feature in the PLS model. A TP-loading closer to −1 or 1 indicates that the gait feature has a strong influence in discriminating between fallers and nonfallers, whereas a TP-loading close to 0 indicates that the feature has little or no influence in discriminate analyses. Thus, the TP-loading provides a ranking list of the most influential gait features for the classification of fallers and nonfallers. Three different predictor matrices were defined for the PLS-DA to compare the discrimination ability of phase-dependent local dynamic stability with other gait features (see Table 1). All gait features in the predictor matrices were converted to -scores before the application of PLS-DA. A PLS-DA cross-validation procedure was used to estimate how well the model would generalize to new samples from the same population [28]. Four latent variables provided the minimum error of the cross-validation for all predictor matrices and were used in the PLS-DA. Sensitivity and specificity and area under the ROC curve (AUC) were defined based on the real and predicted outcome variables from PLS-DA. All analyses were performed in Matlab R2014a.

Table 1: The gait features contained in the three predictor matrices used in the partial least square discriminatory analysis (PLS-DA) of elderly fallers and nonfallers. The gait features written in italic style are the same features used in Weiss et al. (2013) [13].

Predictor matrix (46 gait features)
Predictor matrix (38 gait features)
Predictor matrix (8 gait features)

(phase: 0%, (1))
—
(phase: 0%, (1))

(phase: 0%, (2))
—
(phase: 0%, (2))

(phase: 60%, (1))
—
(phase: 60%, (1))

(phase: 60%, (2))
—
(phase: 60%, (2))

(phase: 0%, (1))
—
(phase: 0%, (1))

(phase: 0%, (2))
—
(phase: 0%, (2))

(phase: 60%, (1))
—
(phase: 60%, (1))

(phase: 60%, (2))
—
(phase: 60%, (2))

—

Acceleration range
Acceleration range
—

Acceleration root-mean-square
Acceleration root-mean-square
—

Amplitude of dominant frequency
Amplitude of dominant frequency
—

Average stride duration
Average step duration
—

Average step duration
Average step duration
—

Cadence
Cadence
—

Harmonic ratio
Harmonic ratio
—

Median walking bout duration
Median walking bout duration
—

Median number of steps for bout
Median number of steps for bout
—

Slope of dominant frequency
Slope of dominant frequency
—

Step symmetry
Step symmetry
—

Step regularity
Step regularity
—

Stride regularity
Stride regularity
—

Total number of steps
Total number of steps
—

Total number of walking bouts
Total number of walking bouts
—

Total percent of walking duration
Total percent of walking duration
—

Width of dominant frequency
Width of dominant frequency
—

Gait feature defined for AP, ML, and V direction, separately.

Figure 2: A schematic illustration of partial least square discriminant analysis (PLS-DA) and target projection (TP) used in the present study. The prediction matrix has a number of columns equal to the number of gait features and number of rows equal to the number of participants. The high number of gait features is projected to a small number of principle axes in the feature space and this projection is defined as the -scores (). The projection is provided by a weight matrix for the gait feature matrix and by a weight matrix for the categorical (faller and nonfaller) response variable . The target projection method combines the weight matrices and in order to define the influence of each gait feature in the discrimination between fallers and nonfallers. The cross-product of weights and is used to calculate the target projection scores , which is a variable containing one score for each older person which maximizes the discrimination between fallers and nonfallers. The target projection score is used to define the target projection loadings TP containing one loading for each gait feature that denotes its influence on the target projection score.

3. Results

Figure 3 shows the mean reaction curve of the fallers and nonfallers for initial perturbation at 0%, 20%, 40%, 60%, and 80% of the step cycle. The figure indicates that the reaction distance curve for both elderly fallers and nonfallers has a phase-dependent shape. The nonfallers had a significantly larger local dynamic stability, median λ, compared to the fallers at 0%, 20%, and 60% of the step cycle, irrespective of the state space reconstruction method and definition of (see upper and lower panels in Figure 3). The TP-loadings of all 46 included gait features indicate that phase-dependent at 0 and 60% of step cycle was most influential in discriminating between elderly fallers and nonfallers (see green bars in Figure 4). The conventional measures of local dynamic stability had less influence in the discrimination analysis compared to the phase-dependent (compare red bars of with green bars of and in Figure 4). Increased error, decreased specificity, and decreased AUC were found when the phase-dependent measures, and , were removed (see Table 2 and Figure 5). Furthermore, the eight phase-dependent measures, and , performed as good as all the 38 gait features together in classification of fallers and nonfallers (see third column in Table 2 and the red ROC curve in Figure 5). In addition, all the phase-dependent measures, median and , included in the discriminate analysis had high test-retest reliability (ICC coefficients > 0.80; see Table 3). Thus, phase-dependent local dynamic stability seems to be an important feature for the classification of fallers and nonfallers in community-dwelling older persons.

Table 2: Classification performance for predictor matrices , , and (see Table 1 for their definitions).

Predictors
Predictors
Predictors

Sensitivity
0.72
0.72
0.69

Specificity
0.90
0.79
0.87

AUC
0.93
0.84
0.83

Error (1 – accuracy)
0.18
0.24
0.21

Table 3: Interclass correlation (ICC) coefficient and its 95% confidence interval (CI) for the phase-dependent local dynamic stability measures, and .

Features
ICC
ICC (95% CI)

(phase: 0%, (1))
0.90
[0.84, 0.94]

(phase: 0%, (2))
0.89
[0.82, 0.93]

(phase: 60%, (1))
0.92
[0.88, 0.95]

(phase: 60%, (2))
0.90
[0.85, 0.94]

(phase: 0%, (1))
0.86
[0.77, 0.91]

(phase: 0%, (2))
0.85
[0.76, 0.91]

(phase: 60%, (1))
0.93
[0.88, 0.95]

(phase: 60%, (2))
0.92
[0.87, 0.95]

(a) State space construction Method 1 (differential coordinate embedding)

(b) State space construction Method 2 (delayed coordinate embedding)

Figure 3: (a) of the log-median for fallers (red) and nonfallers (blue) defined by the state space reconstruction Method 1 (differential coordinate embedding) for initial perturbation at 0%, 20%, 40%, 60%, and 80% of the step cycle. The smaller upper subplots show of local dynamic stability for fallers (red) and nonfallers (blue) together with values. (b) of the log-median for fallers (red) and nonfallers (blue) defined by the state space reconstruction Method 2 (delay coordinate embedding) for initial perturbation at 0%, 20%, 40%, 60%, and 80% of the step cycle. The smaller upper subplots show of local dynamic stability for fallers (red) and nonfallers (blue) together with values.

Figure 4: The TP-loading scores and corresponding values for 46 gait features (predictor matrix in Table 1). The feature numbers in the middle are linked to the list of gait features as follows: (1) total number of walking bouts, (2) ML harmonic ratio, (3) total percent of walking duration, (4) V step symmetry, (5) AP amplitude of dominant frequency, (6) ML step regularity, (7) V acceleration root-mean-square, (8) ML acceleration root-mean-square, (9) total number of steps, (10) V harmonic ratio, (11) ML step symmetry, (12) AP acceleration root-mean-square, (13) AP slope of dominant frequency, (14) AP stride regularity, (15) ML stride regularity, (16) V step regularity, (17) ML acceleration range, (18) V stride regularity, (19) AP step symmetry, (20) median walking bout duration, (21) V acceleration range, (22) AP harmonic ratio, (23) ML width of dominant frequency, (24) V , (25) AP step regularity, (26) V slope of dominant frequency, (27) V width of dominant frequency, (28) ML , (29) AP width of dominant frequency, (30) AP , (31) V amplitude of dominant frequency, (32) ML amplitude of dominant frequency, (33) ML slope of dominant frequency, (34) median number of steps for bout, (35) AP acceleration range, (36) (phase: 0%, (2)), (37) (phase: 0%, (2)), (38) Cadence, (39) average stride duration, (40) average step duration, (41) (phase: 60%, (2)), (42) (phase: 0%, (1)), (43) (phase: 0%, (1)), (44) (phase: 60%, (1)), (45) (phase: 60%, (1)), and (46) (phase: 60%, (2)). The phase-dependent local dynamic stability measures, and , are represented as green bars whereas conventional local dynamic stability measures, , are represented as red bars. The yellow bars represent gait features used in Weiss et al. (2013).

Figure 5: ROC curves summarizing the performance of the PLS-DA for the predictor matrix (blue), predictor matrix (green), and predictor matrix (red) of gait features (see Table 1 for their definitions).

4. Discussion

The main purpose of the present study was to compare phase-dependent local dynamic stability measures with more conventional gait features in their ability to discriminate between community-dwelling elderly fallers and nonfallers. The phase-dependent and at 0% and 60% of the step cycle had the best classification performance and considerably improved the PLS-DA, compared to the 38 conventional features of daily life walking.

In a treadmill study, Ihlen et al. [17] found that healthy older persons had larger phase-dependent and compared to young adults. In contrast, in the present study, the phase-dependent and were larger in older nonfallers compared to fallers. These contrasting findings may be due to several reasons. First, a recent study indicates that gait characteristics obtained in in-lab studies are different from those recorded for daily life walking [29]. Thus, the more unstable gait dynamics of elderly nonfallers might be due to a more heterogeneous and challenging walking environment for this group including more frequent turns and multitasking while walking. These factors might contribute to the larger and compared to the fallers but might also indirectly reflect subtle decline in balance and mobility in the fallers group. Second, the increased local dynamic stability found for the fallers might also reflect adaption in the gait dynamics towards more cautious gait including less frequent turns and multitasking while walking and less challenging walking environment. However, the present study cannot conclude whether the larger and in the nonfallers are due to differences in external environmental factors or internal neuromuscular factors or a combination of these two. Further studies are needed to include contextual information of daily life walking in community-dwelling older persons.

Rispens et al. [12] found that local dynamic stability, , computed by the method of Wolf et al. [23] for the acceleration signal in the V direction was able to discriminate between fallers and nonfallers. However, in the present study, the phase-dependent and were found to be more sensitive to falls status of community-dwelling older person compared to . In contrast to the findings of Rispens et al. [12], did not influence the classification of fallers and nonfallers in the present study.

The present study also shows the potential of PLS-DA for the comparison of the influence of different gait features to discriminate between fallers and nonfallers. Numerous features of gait stability and variability have been introduced in the last decades, but their abilities to discriminate between fallers and nonfallers are seldom compared [11]. The TP-loadings in Table 2 are able to rank the influence of the different gait features in the classification of elderly fallers and nonfallers. Table 2 indicates that many of these features have significantly different mean values for fallers and nonfallers, while the discriminatory power is low (i.e., TP-loading < 0.5). A consensus on a procedure to compare the abilities of different gait features in the classification of elderly fallers and nonfallers, like PLS-DA, might have important value for evaluation of new features of gait stability and variability. Furthermore, procedures like PLS-DA might also be helpful for the identification of fall risk profiles for different groups of older people, and the procedures might be extended to include clinical test scores and demographic variables.

The present study has several limitations. First, the present study did only distinguish between fallers and nonfallers based on retrospective fall reports from a relatively small sample of community-dwelling older adults. The specificity, sensitivity, and AUC reported in the present study are in the upper end of values that could be expected from a perfect fall prediction model [30]. Thus, it is likely that specificity, sensitivity, and AUC will decrease for and in a prediction model of prospective falls. Consequently, further studies on larger samples with prospective fall data are necessary before concluding that phase-dependent and will improve fall prediction models or early fall risk assessment in the population of community-dwelling older adults.

Second, demographic variables and variables of clinical tests used for fall risk assessment, like tests of balance and mobility performance, were not included in the classification procedure. Inertial sensor based tools for unsupervised in-home test of physical function, including mobility and balance, could also contribute to the improvement of early fall risk assessment in community-dwelling older adults [31]. However, former studies have shown that features of daily life walking improve the risk assessment when combined with instrumented tests of mobility performance [10, 13]. Nevertheless, falls have multifactorial causes including medication, urinary control, vision, footwear, environmental hazards, cognitive function, mental health, and fear of falling, to mention but a few, and it is therefore likely that a combination of outcomes of clinical tests and features of daily life activities will optimize fall risk assessment and fall prediction models. Even though the inclusion of and might further improve the fall risk assessment when combined with clinical tests, issues like cost (and maintenance cost) of accelerometers, unsupervised device handling in an in-home setting, provision and retrieval from patient in a clinical setting, and the potential for an easy-to-use online estimation of and will decide the feasibility of the use of and in fall risk assessment tools. Thus, further studies and cost-benefit analyses have to be conducted to determine the usability and feasibility of these analyses when implemented in smartphone and desktop application and the gain in accuracy of fall risk assessment needed to compensate potential decline in clinical feasibility.

Third, the relationship between phase-dependent stability and and variables related to the health status of the older adults was not investigated. Investigation of these relationships would be important to improve the clinical interpretation of and and thus should be included in further studies. In addition, assessment of and could be combined with experimental in-lab research of stability, like experimental perturbation studies of in-lab gait, as well as studies on neurophysiological mechanisms in animal models to improve the understanding of the underlying mechanisms of and [32, 33].

Fourth, the time consumption of the computational steps (i.e., gait bout identification, preprocessing, and estimation procedure) to assess and was not recorded. The time consumption of these steps would be important to decide the possibility for online computation of and which is necessary for clinical feasible smartphone and desktop application.

Fifth, even if the present study did investigate the test-retest reliability of and based on the first and the last 1/3 of the walking bouts during the 3-day recording, this is a considerable shorter test-retest interval compared to the 1-week test-retest interval considered in Rispens et al. [12]. ICC for 1-week test-retest might be weaker compared to the ICC found in the present study and further studies should investigate test-retest reliability of the local divergence features for longer test-retest intervals.

Sixth, the accuracy of the phase-dependent and is dependent on the reliability of the step identification. The inertial sensor was placed on the lower back which makes heel strike and toe-off events more difficult to identify within the gait cycle. Thus, the phase-dependent and in Figure 3 were not defined according to single and double support phases within the gait cycle, but according to the local peaks of the vertical velocity. The employment of advanced step identification algorithms might define the phase-dependent and according to heel strike and toe-offs, but further validation of these algorithms is necessary [34]. Furthermore, as inertial sensors become smaller and more wearable, further studies should include an additional sensor on the lower extremities and/or insole data to identify heel strikes and toe-offs and thereby single and double support phases.

Seventh, the sample size used in the present study is small. Rispens et al. [12] were not able to replicate the results of Weiss et al. [13] for some of the spectral features for another study with a larger sample size. A similar contrast in results might be present for phase-dependent local dynamic stability when replicated for different groups of community-dwelling older persons. Thus, further studies should replicate these initial findings on cohorts of community-dwelling older persons with different health status. Finally, the present study suggests that phase-dependent and are related to falls status and might be important to include in fall risk assessments and fall prediction models. Several studies indicate that measures of gait stability and variability improve fall risk assessment and fall prediction models when compared to assessments and models based on clinical tests and fall history [12, 13]. Thus, further application of phase-dependent and , to track changes in falls status and to prospectively identify fallers, is needed to determine their influence in fall risk assessments and fall prediction models.

5. Conclusions

The present study compared phase-dependent measures of local dynamic stability of elderly fallers and nonfallers in daily life walking with existing features of gait stability and variability. These phase-dependent measures had the best classification performance of all included gait features and improved the discrimination between elderly fallers and nonfallers compared to all other features of daily life walking. Thus, phase-dependent measures of local dynamic stability might be of importance for further development in early fall risk assessment, fall prediction, and fall prevention amongst community-dwelling older persons. The present results set the stage for follow-up prospective studies in larger cohorts and clinical feasibility studies to further assess the potential of these metrics.

Conflict of Interests

The authors declare that no conflict of interests is associated with the present study.

Acknowledgments

Both the long-term recordings on which the present analyses were made and the Matlab functions for the computation of phase-dependent and will be available at http://www.physionet.org/, the National Institutes of Health-Sponsored Research Resource for Complex Physiologic Signals. Local dynamic stability measure, , was computed by Matlab function available at http://www.mathworks.com/matlabcentral/fileexchange/48084-lyapunov-exponent-estimation-from-a-time-series. This work was funded by the Norwegian Research Council (FRIMEDBIO, Contract no. 230435) and in part by the European Commission (WIISEL, FP7-ICT-2011-7-ICT-2011.5.4-Contract no. 288878).

References

S. Heinrich, K. Rapp, U. Rissmann, C. Becker, and H.-H. König, “Cost of falls in old age: a systematic review,” Osteoporosis International, vol. 21, no. 6, pp. 891–902, 2010.

S. Deandrea, E. Lucenteforte, F. Bravi, R. Foschi, C. La Vecchia, and E. Negri, “Risk factors for falls in community-dwelling older people: a systematic review and meta-analysis,” Epidemiology, vol. 21, no. 5, pp. 658–668, 2010.

D. Oliver, M. Britton, P. Seed, F. C. Martin, and A. H. Hopper, “Development and evaluation of evidence based risk assessment tool (STRATIFY) to predict which elderly inpatients will fall: case-control and cohort studies,” British Medical Journal, vol. 315, no. 7115, pp. 1049–1053, 1997.

M. Raîche, R. Hébert, F. Prince, and H. Corriveau, “Screening older adults at risk of falling with the Tinetti balance scale,” The Lancet, vol. 356, no. 9234, pp. 1001–1002, 2000.

S. Nandy, S. Persons, C. Cryer et al., “Development and preliminary examination of the predictive validity of the Falls Risk Assessment Tool (FRAT) for use in primary care,” Journal of Public Health, vol. 26, no. 2, pp. 138–143, 2004.

A. M. Tromp, S. M. F. Pluijm, J. H. Smit, D. J. H. Deeg, L. M. Bouter, and P. Lips, “Fall-risk screening test: a prospective study on predictors for falls in community-dwelling elderly,” Journal of Clinical Epidemiology, vol. 54, no. 8, pp. 837–844, 2001.

C. Becker, L. Schwickert, S. Mellone et al., “Proposal for a multiphase fall model based on real-world fall recordings with body-fixed sensors,” Zeitschrift für Gerontologie und Geriatrie, vol. 45, no. 8, pp. 707–715, 2012.

S. N. Robinovitch, F. Feldman, Y. Yang et al., “Video capture of the circumstances of falls in elderly people residing in long-term care: an observational study,” The Lancet, vol. 381, no. 9860, pp. 47–54, 2013.

F. Horak, L. King, and M. Mancini, “Role of body-worn movement monitor technology for balance and gait rehabilitation,” Physical Therapy, vol. 95, no. 3, pp. 461–470, 2015.

T. Iluz, A. Weiss, E. Gazit et al., “Can a body-fixed sensor reduce Heisenberg’s uncertainty when it comes to the evaluation of mobility? Effects of aging and fall risk on transitions in daily living,” The Journals of Gerontology Series A: Biological Sciences and Medical Sciences, 2015.

D. Hamacher, N. B. Singh, J. H. van Dieën, M. O. Heller, and W. R. Taylor, “Kinematic measures for assessing gait stability in elderly individuals: a systematic review,” Journal of the Royal Society Interface, vol. 8, no. 65, pp. 1682–1698, 2011.

S. M. Rispens, K. S. van Schooten, M. Pijnappels, A. Daffertshofer, P. J. Beek, and J. H. Van Dieën, “Identification of fall risk predictors in daily life measurements: gait characteristics' reliability and association with self-reported fall history,” Neurorehabilitation and Neural Repair, vol. 29, no. 1, pp. 54–61, 2015.

A. Weiss, M. Brozgol, M. Dorfman et al., “Does the evaluation of gait quality during daily life provide insight into fall risk? A novel approach using 3-Day accelerometer recordings,” Neurorehabilitation and Neural Repair, vol. 27, no. 8, pp. 742–752, 2013.

S. M. Bruijn, O. G. Meijer, P. J. Beek, and J. H. van Dieën, “Assessing the stability of human locomotion: a review of current measures,” Journal of the Royal Society Interface, vol. 10, no. 83, Article ID 20120999, 2013.

J. B. Dingwell and J. P. Cusumano, “Nonlinear time series analysis of normal and pathological human walking,” Chaos, vol. 10, no. 4, pp. 848–863, 2000.

E. A. F. Ihlen, T. Goihl, P. B. Wik, O. Sletvold, J. Helbostad, and B. Vereijken, “Phase-dependent changes in local dynamic stability of human gait,” Journal of Biomechanics, vol. 45, no. 13, pp. 2208–2214, 2012.

E. A. F. Ihlen, O. Sletvold, T. Goihl, P. B. Wik, B. Vereijken, and J. Helbostad, “Older adults have unstable gait kinematics during weight transfer,” Journal of Biomechanics, vol. 45, no. 9, pp. 1559–1565, 2012.

A. Weiss, S. Sharifi, M. Plotnik, J. P. P. van Vugt, N. Giladi, and J. M. Hausdorff, “Toward automated, at-home assessment of mobility among patients with Parkinson disease, using a body-worn accelerometer,” Neurorehabilitation and Neural Repair, vol. 25, no. 9, pp. 810–818, 2011.

Y. Lin, W. Chong-Yang, and W. Zhi-Cheng, “Research on the selection of wavelet bases for wavelet-based signal trend elimination,” in Proceedings of the International Conference on Wavelet Active Media Technology and Information Processing (ICWAMTIP '12), pp. 20–24, Chengdu, China, December 2012.

H. Kantz and T. Schreiber, Nonlinear Time Series Analysis, Cambridge University Press, Cambridge, UK, 2nd edition, 2004.

R. Moe-Nilssen and J. L. Helbostad, “Estimation of gait cycle characteristics by trunk accelerometry,” Journal of Biomechanics, vol. 37, no. 1, pp. 121–126, 2004.

M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,” Physica D, vol. 65, no. 1-2, pp. 117–134, 1993.

A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena, vol. 16, no. 3, pp. 285–317, 1985.

H. D. Abarbanel, Analysis of Observed Chaotic Data, Springer, New York, NY, USA, 1996.

O. M. Kvalheim and T. V. Karstang, “Interpretation of latent-variable regression models,” Chemometrics and Intelligent Laboratory Systems, vol. 7, no. 1-2, pp. 39–51, 1989.

S. Wold, K. Esbensen, and P. Geladi, “Principal component analysis,” Chemometrics and Intelligent Laboratory Systems, vol. 2, no. 1–3, pp. 37–52, 1987.

S. Wold, M. Sjöström, and L. Eriksson, “PLS-regression: a basic tool of chemometrics,” Chemometrics and Intelligent Laboratory Systems, vol. 58, no. 2, pp. 109–130, 2001.

J. A. Westerhuis, H. C. J. Hoefsloot, S. Smit et al., “Assessment of PLSDA cross validation,” Metabolomics, vol. 4, no. 1, pp. 81–89, 2008.

S. M. Rispens, Walking on the wild side: fall risk assessment from daily-life gait [Ph.D. thesis], VU University, Amsterdam, The Netherlands, 2014.

P. Palumbo, L. Palmerini, and L. Chiari, “A probabilistic model to investigate the properties of prognostic tools for falls,” Methods of Information in Medicine, vol. 54, no. 2, pp. 189–197, 2015.

T. Shany, S. J. Redmond, M. Marschollek, and N. H. Lovell, “Assessing fall risk using wearable sensors: a practical discussion,” Zeitschrift für Gerontologie und Geriatrie, vol. 45, no. 8, pp. 694–706, 2012.

M. D. Chang, E. Sejdić, V. Wright, and T. Chau, “Measures of dynamic stability: detecting differences between walking overground and on a compliant surface,” Human Movement Science, vol. 29, no. 6, pp. 977–986, 2010.

L. Hak, H. Houdijk, F. Steenbrink et al., “Speeding up or slowing down? Gait adaptations to preserve gait stability in response to balance perturbations,” Gait & Posture, vol. 36, no. 2, pp. 260–264, 2012.

A. Godfrey, S. Del Din, G. Barry, J. C. Mathers, and L. Rochester, “Instrumenting gait with an accelerometer: a system and algorithm examination,” Medical Engineering & Physics, vol. 37, no. 4, pp. 400–407, 2015.