Next Article in Journal
Smartphone-Based Ecological Momentary Assessment for Collecting Pain and Function Data for Those with Low Back Pain
Next Article in Special Issue
Gait Characteristics Associated with Fear of Falling in Hospitalized People with Parkinson’s Disease
Previous Article in Journal
Blind Fault Extraction of Rolling-Bearing Compound Fault Based on Improved Morphological Filtering and Sparse Component Analysis
Previous Article in Special Issue
MaD GUI: An Open-Source Python Package for Annotation and Analysis of Time-Series Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gait Detection from a Wrist-Worn Sensor Using Machine Learning Methods: A Daily Living Study in Older Adults and People with Parkinson’s Disease

by
Yonatan E. Brand
1,2,†,
Dafna Schwartz
1,3,†,
Eran Gazit
1,
Aron S. Buchman
4,
Ran Gilad-Bachrach
2,3,5 and
Jeffrey M. Hausdorff
1,2,6,7,*
1
Center for the Study of Movement, Cognition and Mobility, Neurological Institute, Tel Aviv Sourasky Medical Center, Tel Aviv 6492416, Israel
2
Sagol School of Neuroscience, Tel Aviv University, Tel Aviv 6997801, Israel
3
Department of Biomedical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv 6997801, Israel
4
Rush Alzheimer’s Disease Center, Department of Neurological Sciences, Rush University Medical Center, Chicago, IL 60612, USA
5
Edmond J. Safra Center for Bioinformatics, Tel-Aviv University, Tel Aviv 6997801, Israel
6
Rush Alzheimer’s Disease Center and Department of Orthopedic Surgery, Rush University, Chicago, IL 60612, USA
7
Department of Physical Therapy, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv 6997801, Israel
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 1 August 2022 / Revised: 31 August 2022 / Accepted: 12 September 2022 / Published: 19 September 2022
(This article belongs to the Special Issue Inertial Sensors for Clinically Relevant Mobility Outcome Measures)

Abstract

:
Remote assessment of the gait of older adults (OAs) during daily living using wrist-worn sensors has the potential to augment clinical care and mobility research. However, hand movements can degrade gait detection from wrist-sensor recordings. To address this challenge, we developed an anomaly detection algorithm and compared its performance to four previously published gait detection algorithms. Multiday accelerometer recordings from a wrist-worn and lower-back sensor (i.e., the “gold-standard” reference) were obtained in 30 OAs, 60% with Parkinson’s disease (PD). The area under the receiver operator curve (AUC) and the area under the precision–recall curve (AUPRC) were used to evaluate the performance of the algorithms. The anomaly detection algorithm obtained AUCs of 0.80 and 0.74 for OAs and PD, respectively, but AUPRCs of 0.23 and 0.31 for OAs and PD, respectively. The best performing detection algorithm, a deep convolutional neural network (DCNN), exhibited high AUCs (i.e., 0.94 for OAs and 0.89 for PD) but lower AUPRCs (i.e., 0.66 for OAs and 0.60 for PD), indicating trade-offs between precision and recall. When choosing a classification threshold of 0.9 (i.e., opting for high precision) for the DCNN algorithm, strong correlations (r > 0.8) were observed between daily living walking time estimates based on the lower-back (reference) sensor and the wrist sensor. Further, gait quality measures were significantly different in OAs and PD compared to healthy adults. These results demonstrate that daily living gait can be quantified using a wrist-worn sensor.

1. Introduction

Wrist-worn sensors are widely used in research studies and clinical trials to quantify movement. Varied behaviors of daily living, including sleep and circadian metrics, and the diverse facets of physical activity can be extracted and quantified from recordings obtained from an inertial measurement unit (IMU) such as an accelerometer and gyroscope. These metrics have been shown to predict morbidity, mortality, cognitive decline, and the development of dementia [1,2,3] and have also been employed for the early detection and monitoring of neurodegenerative disorders that commonly manifest impaired gait (e.g., Alzheimer’s and Parkinson’s disease (PD)) [4,5,6].
Many previous studies using wrist-worn sensors concentrated on quantifying the level of physical activity, predominantly a reflection of the total amount of daily walking [7,8,9]. As gait is a relatively complex behavior, recent efforts have focused on developing algorithms that can extract additional facets of gait quality (e.g., gait speed, cadence, and step regularity) from these same recordings to more fully characterize gait function [10,11,12]. Recent works suggest that quantifying more diverse facets of gait may improve predictions of specific adverse health outcomes [13,14]. Moreover, multiday recordings of gait during daily living may also enhance functional monitoring and the prediction of adverse outcomes such as falls, as compared to short assessments of walking obtained within or outside the lab setting [15,16,17,18,19,20].
To extract measures of gait quality and quantity from acceleration recordings of daily living, it is first necessary to detect the multiple bouts of gait that are contained in a multiday recording. However, reliable detection of walking from wrist-sensor data is challenging due to the many degrees of freedom of wrist movement and because wrist movements may occur in the absence of walking (e.g., when washing dishes or brushing teeth). In other words, the continuous acceleration signals obtained from a wrist sensor do not, a priori, contain reliable labels of gait/non-gait behaviors during daily living. To circumvent this limitation, many previous studies on gait have placed sensors on the lower back or on a lower extremity [4,21,22,23], where it is much easier to reliably identify the bouts of gait.
Existing gait detection algorithms distinguish gait from other behaviors using features from the time and frequency domain, either individually or combined, using conventional statistics [12,24,25]. Once these features are extracted, they can be used as the input for a machine learning model to automatically detect bouts of gait. These include both supervised [26,27] and unsupervised [28] methods. More recently, deep learning techniques have been applied to human activity recognition (HAR) tasks to classify human activity (e.g., walking, sitting, and driving) [29,30,31,32]. Deep learning models can deal with complex and large datasets and offer an approach for automatic gait detection from wrist recordings of daily living [31].
Only a few studies have developed and validated gait detection algorithms using wrist-sensor recordings of daily living [11,26,27,28,33,34]. One of the challenges of developing daily living gait detection algorithms is that walking during daily living occurs for only a small portion of the overall recording, as non-gait behaviors constitute the bulk of the overall recording. Thus, false-positive predicted bouts of gait occur frequently when analyzing daily living recordings. Hence, there is a significant trade-off between precision (i.e., a positive predictive value) and sensitivity (i.e., the percent of correctly detected gait bouts, also referred to as recall) of the algorithm. Thus, focusing on algorithm performance statistics, such as accuracy and specificity, without assessing precision, as has been done in the past, [28,30,31], will not provide an adequate assessment of the algorithm’s performance.
Many prior algorithms have been developed for healthy, younger adults. Gait detection is more difficult in OAs who commonly manifest diverse gait impairments even in the absence of specific diseases. Aging is often accompanied by reduced frequency of gait bouts, and the walking patterns in OAs are less periodic and consistent, increasing the likelihood of false-negative detections of gait. In addition, the presence of arm tremor or reduced arm swing during walking [35] may reduce the strength of the signal obtained from the wrist, making automatic detection of gait more difficult.
The use of unsupervised anomaly detection algorithms has become a useful approach in cases where labeling is not feasible and examples of one class are sparse [36]. Therefore, it seems reasonable to adopt this approach for the detection of gait from wrist-worn sensors in OAs with and without overt neurological diseases such as Parkinson’s disease (PD). It is possible to address anomaly detection as a one-class classification algorithm that first learns a model of one class during training (which can be referred to as the normal data in classical anomaly detection tasks) and then at the inference time computes an anomaly score of unseen data samples under the learned model [37,38]. Data samples are flagged as belonging to the learned model or not by comparing the anomaly score to a learned decision boundary. In recent years, the practice of normalizing flows, which are deep neural networks based on density estimators, has shown promising performance in different unsupervised anomaly detection tasks [38,39,40,41].
The present study adopted an anomaly detection algorithm to detect gait from wrist-worn 3D accelerometer recordings of daily living. We first trained and tested our model using accelerometer recordings from healthy young adults (HYAs) to examine the feasibility of our model. Then, we trained and tested the model using wrist recordings from OAs with and without PD. We compared the performance of the anomaly detection approach to four previously published gait detection algorithms. We extracted measures of gait quantity and gait quality from the best model and compared its performance to those obtained via a sensor worn on the lower back (a form of “gold-standard” reference).

2. Materials and Methods

2.1. Participants

Wrist-worn recordings of HYAs were obtained from the publicly available PhysioNet dataset [42,43], which included 32 healthy participants (39 ± 9 years, 19 women) who carried out various activities of daily living. The older adult participants were recruited from the Center for the Study of Movement, Cognition, and Mobility (CMCM) at the Neurological Institute in Tel Aviv Sourasky Medical Center, Israel, using a convenience sample of 12 community older adults (OAs) and 18 individuals with PD. The diagnosis of PD was obtained by a movement disorders specialist.
To characterize the OAs and the people with PD, several widely used measures were applied. In the lab, all subjects performed a 2 min walk test (2MWT) while walking at their comfortable speed along a 25 m corridor. A Zeno Walkway Gait Analysis System (Protokinetics LLC, Havertown, PA, USA) with a length of 7.92 m was placed in the middle of the walkway and used to measure the subject’s walking speed. The timed up and go test (TUG) measures the time required to stand up from a seated position, walk 3 m, turn around, and return to a seated position. TUG times above 13.5 s suggest the presence of gait impairments and are associated with an increased risk of falls [44]. The Montreal Cognitive Assessment (MOCA) was used as a general measure of cognitive function. Hoehn and Yahr staging (performed by a movement disorders clinician) was used to describe disease severity among the people with PD.

2.2. Wearable Sensors

In the HYA dataset, accelerometry data were collected simultaneously at four body locations—left wrist, left hip, left ankle, and right ankle—at a sampling frequency of 100 Hz (ActiGraph GT3X+, Pensacola, FL, USA). We used only the left wrist data, as we aimed to detect gait from a wrist-worn sensor. The data included labels of activity types performed for each time point in the data collection.
The OAs and people with PD were asked to wear a tri-axial accelerometer on their wrist for up to 10 days, at a sampling frequency of 25 Hz (Garmin, Olathe, KS, USA). The OA participants wore the sensor on their nondominant hand, and the PD participants wore it on their impaired hand (with the idea that the more impaired hand would roughly behave similarly to the nondominant hand and be less involved with non-gait activities). The raw data were uploaded online to cloud services. An additional sensor was placed on the lower back at the level of the fifth lumbar vertebrae (roughly near the center of mass) and the 3D accelerometer signals were sampled at 100 Hz (Axivity AX3, York, UK). Data from the lower back were stored locally and uploaded offline to the lab computers.
Labeling of the OA and PD datasets to gait and non-gait segments was conducted based on “gold-standard” annotations from the lower-back sensor [4,21,23,45]. The signal from the lower back was less noisy than the wrist and had validated algorithms for detecting gait [46,47]. We used nearest-neighbor interpolation to downsample the annotations from the lower back to match the wrist-sensor sampling rate. We considered only walking bouts with a duration of at least 6 s as “gait”, since shorter duration bouts are extremely difficult to reliably detect, even from a lower-back sensor, and since most gait quality measures cannot be applied to very short walking bouts [48].

2.3. Gait Detection Algorithms

In this study, we used both supervised and unsupervised algorithms for gait detection. We split the data into training and testing sets in the supervised models. For the HYA dataset, we used a 70/30 split validation. For the OAs and PD, we used a five-fold cross-validation to split the daily living data into the testing and training sets. Each fold consisted of data from different participants (i.e., the model did not train and test on data from the same participant in each fold) to minimize overfitting of the model. The models were trained and evaluated on each cohort (i.e., HYAs, OAs, and PD) separately, for customizing and optimizing the classifier to the characteristics of each cohort. The different folds included similar percentages of gait, such that they were relatively balanced (15% ± 4.8% for the OAs and 13% ± 2.7% for the PD group).

2.3.1. Anomaly Detection

We evaluated here, for the first time, an approach of unsupervised anomaly detection for the task of gait detection. We applied the algorithm to the HYA cohort as an initial feasibility test, which allowed us to examine the performance on a simpler and cleaner dataset, and then we applied it to the more complex cohorts of PD and OAs. We addressed the gait detection problem as a one-class classification problem in which the neural network was trained using samples from one class only (gait/non-gait), and this class was considered the “normal” class. The dataset was split into a training set, which consisted of 90% of the “normal” class data, and a validation set, which consisted of the remaining 10% of that “normal” class used in inference time, and a test set which consisted of all the data from the second class (i.e., the “abnormal” class). We examined the model both when gait data were considered normal and used for training, and when non-gait data were considered normal and used for training. At the inference time, to decide if a sample belonged to the “normal” class (i.e., the class that the model was trained on) or to the “abnormal” class (i.e., that the model was not exposed to during training), an anomaly score was calculated, and using a threshold, the samples were classified so that samples that obtained a higher anomaly score than the threshold were considered abnormal, and a sample that obtained a lower score than the threshold was considered as normal. We based our anomaly detection model on the masked autoregressive flow (MAF), which is a type of normalizing flow suitable for density estimation [49].
Normalizing flows are neural networks that are used to learn bijective transformations between data distributions and well-defined base densities. The bijectivity is ensured by stacking layers of affine transforms that are fixed or autoregressive. Under this invertibility assumption, the learned density p x can be calculated as:
    p x = π u f 1 x det f 1 x
where π u is the well-defined base density, which is usually chosen to be a standard Gaussian density; f is the invertible differential transformation. The MAF model uses a stack of MADE layers, which are a type of autoregressive flow [50]. Our architecture consisted of a stack of 5 MADE layers. The goal during training was to maximize the likelihood or, equivalently, the log-likelihood of the training data with respect to the network parameters. Low-activity segments were removed from the accelerometer signal; a moving window STD was used to remove windows with a low STD (i.e., low activity). We used moving windows with a length of 3 s and set the STD threshold to 0.1. After, principal component analysis (PCA) was conducted for dimensionality reduction. The signals were then split into sliding windows at a fixed width of 6 s with an overlap of 50%, and the DC offset was removed by subtracting the mean value in each window on the three axes. We calculated the mean and standard deviation (STD) of the training set and used it to normalize the training set, validation set, and testing set. Hyperparameters were found for each model of the three cohorts through a set of distinct experiments. Finally, we used an ADAM optimizer with an exponential decay learning rate scheduler with an initial value of 0.01 for the HYA, OA, and PD datasets and an l 2 regularization of 0.01 for the PD and OA datasets and 0.001 for the HYA dataset. At the inference time, the log-likelihoods of the data samples from the validation set and testing set were computed and used to classify the data according to a selected likelihood threshold.

2.3.2. Feature-Based Model

Another approach for gait detection is a feature-based (FB) approach. In this method of gait detection, features are extracted from time windows of fixed width for classification [12,24,25]. These features are usually based on the time and frequency domain, either individually or combined. Once features are extracted, we can define thresholds to predict if the window includes gait. In this study, we implemented an FB model published by Keren et al. [24] (Figure 1) for gait detection using a wrist-worn sensor. The original model was developed and validated in an HYA group; here, we applied it to the OAs and people with PD.

2.3.3. Autoregressive Infinite Hidden Markov Model (AR–iHMM)

Classical feature-based models have some challenges when using windows-based features. First, gait features vary between different windows, even in the same participant, depending on the hour of the day, environment, and medication status. Further, extracting frequency domain features relies on the assumption that the window is stationary, which is usually not the case in windows with a fixed width. These problems motivated Raykov et al. [28] to develop a data-driven approach, which did not require the use of windows. In their model, data were automatically segmented into states with variable time lengths. Each state represents a different gait and non-gait pattern. This unsupervised method combines an autoregressive (AR) process, used for parametric estimation of the frequency spectrum, with a hidden Markov model (HMM) used to model the sequence of the different states represented in the data. The input of this model were the raw 3D accelerometer signals. The Euclidean norms of the 3-axis accelerometer signals were calculated as mentioned in the model above. A moving window STD was used to remove windows with a low STD (i.e., low activity). We used moving windows with a length of 3 s and set the STD threshold to 0.1. The output of this algorithm was a vector of the same size as the accelerometry data in which each sample belonged to a specific state. Each state represented an AR model of the power spectral density (PSD). Then, we computed the total PSD of the gait-related frequencies for each state and selected a threshold that maximized the distinction between gait and non-gait states.

2.3.4. Deep Convolutional Neural Network (DCNN)

Another approach we used for gait detection was a DCNN, a type of artificial neural network, commonly applied in the computer vision and image processing fields. Semantic segmentation is a DCNN task used for detecting different clusters in an image, which assigns a separate label for each pixel according to its cluster. The tri-axial accelerometer data were a matrix with a size 3 × times of the samples; therefore, it can be represented as a two-dimensional image. Then, we used a semantic segmentation model to cluster the data into gait and non-gait segments. Zou et al. [31] used this method for gait detection on data that were recorded by smartphones equipped with an IMU system (Figure 2). They treat the tri-axial accelerometer and gyroscope data as a two-dimensional picture and used the DCNN model to apply semantic segmentation. Their model was inspired by the U-Net semantic segmentation algorithm [51]. In this study, we implemented their algorithm with some modifications. We used only the accelerometer data (without gyro), such that the input of the model was the number of mini-batches × 1 × 3 × 512, where each batch included 64 windows consisting of 512 samples for each of the 3 axes of the accelerometer. Further, we used the batch normalization method [52] for stabilizing the learning process and dropout regularization for preventing overfitting [53].

2.3.5. Long Short-Term Memory (LSTM)

The last model that we used for gait detection was based on LSTM. LSTM cells (Figure 3) are variants of recurrent neural networks (RNNs), which are used for problems of sequential data and can handle long-term dependencies using gate functions [54]. Hence, the final algorithm chosen for the gait detection task was an LSTM network. Each LSTM cell consisted of an input gate, a forget gate, and an output gate, which controlled the important information on the current state to be passed forward and which information was redundant. In bidirectional LSTM architecture, the network is trained simultaneously in both time directions.
The improvement of the bidirectional LSTM is that the current output is not only related to previous information in the time series but also to subsequent information in the time series. It was also shown by Wu, Zhang, and Zong (2016) [56] that adding skip connections to the cell outputs with a gated identity function could improve network performance. Yu Zhao et al. [30] proposed a deep network architecture using residual bidirectional long short-term memory (LSTM) cells for HAR. The dataset that was used in their work was a public domain UCI dataset that consisted of 30 young healthy participant readings of smartphone-embedded accelerometers and a gyroscope placed on the waist of the participants. As in the DCNN model, we applied the Deep-Res-Bidir-LSTM architecture to our datasets with several modifications. The input features were only the three axes of the accelerometer data, and the signals were split into sliding windows at a fixed width of 6 s with an overlap of 50%. The dimensions of the classification task were reduced to a binary classification task (i.e., gait/non-gait).

2.4. Algorithm Evaluations and Statistics

Gait detection algorithms aim to perform a binary classification at any point in time to determine whether a given window is gait or non-gait activity. The performance of these algorithms was evaluated using the following statistics.
True positive (TP) was the number of samples that included gaits that were correctly predicted as gaits. True negative (TN) was the number of samples that did not include gaits that were correctly predicted as non-gaits. False-positive (FP) was the number of samples that did not include gaits that were wrongly predicted as gaits. False-negative (FN) was the number of samples that included gaits that were wrongly predicted as non-gaits. From these statistics, we extracted measures to assess the model’s performance:
Accuracy (Acc) is the percent of correct predictions from the total number of samples:
T P + T N T P + F N + T N + F P × 100
Specificity (Sp) represents the ability of the algorithm to correctly predict samples as non-gait among all the non-gait samples:
T N T N + F P × 100
Sensitivity (Se) represents the ability of the algorithm to detect the gait samples:
T P T P + F N × 100
Precision (Pr) is the percent of correct predicted gait samples from the total number of samples predicted as gait (including non-gait samples that were predicted as gaits (i.e., FP)):
T P T P + F P × 100
In addition, we used the receiver operating characteristics (ROC) curve and the precision–recall curve (PR) to visualize the models’ performance. The ROC curve shows the relationship between the true positive rate (i.e., sensitivity) and the false positive rate (i.e., specificity). The PR curve shows the relationship between the precision and the recall (i.e., sensitivity) of the model. In an imbalanced dataset, such as daily living data, the ROC is overly optimistic; thus, the PR curve can be used as an alternative metric to evaluate a model’s performance [57]. We used the area under the curve (AUC) statistic to evaluate the performance in these curves.

2.5. Gait Quantity Measures

To examine the ability of the gait detection models to evaluate gait quantity measures, we computed the correlation between the amount of daily walking according to the lower-back gold standard and according to the wrist-based algorithm. We selected a classification threshold of 0.9 for classifying gait. Only days with at least 8 h of IMU recording were used for each participant. For each day, we calculated the number of minutes spent walking obtained via the lower-back reference value and the wrist-based models. Then, we computed the Pearson correlation between the two measures.

2.6. Gait Quality Measures

Four previously described measures of gait quality were extracted from the best wrist-based model and compared to those obtained via the lower-back sensor to begin to evaluate the ability of the wrist-worn measures to characterize the daily living gait of the study participants. The frequency measures included the dominant frequency (Hz) and amplitude of the dominant harmony of the PSD in the locomotor band (0.5–3 Hz) of the acceleration signal. The amplitude is related to the dominance (or strength) of the frequency in the signal [58]. The root mean square (RMS, m / s 2 ) quantifies the magnitude of the 3D acceleration signal as a measure of gait intensity [59]. Stride regularity is a measure of the consistency of the walking pattern (higher values reflect greater stride-to-stride consistency, and lower values reflect greater stride-to-stride variability) [60].
To increase the certainty that we were extracting gait quality measures from gait bouts and not from false-positive segments, we set a high classification threshold of 0.9. As a result, only samples with a predicted probability greater than 0.9 were classified as gaits, increasing the precision of the model.

3. Results

3.1. Subject Characteristics

Table 1 summarizes the characteristics of the daily living study participants. Both groups (PD and OAs) had relatively high timed up and go (TUG) scores and gait speeds that are typical of OAs. The Hoehn and Yahr scores indicated that the people with PD generally had moderate disease severity.

3.2. Comparison of the Gait Detection Algorithms’ Performances

Table 2 summarizes the performance of the new anomaly detection method compared to four previously published methods. We assessed the models’ performance on the test set. For the OAs and PD, we averaged the performance metrics of the different folds.
Figure 4 reports the ROC and PR curves of all four gait detection models. The curves were calculated from the test set predictions of the five-fold cross-validation. Each curve shows the performance of the model in the three study cohorts: PD (blue), OAs (orange), and HYAs (green). The feature-based model’s output was a binary decision (i.e., gait/non-gait) and not a continuous number that represented probability as in the other models. Therefore, the ROC and PR curves are not presented for this model.

3.3. Gait Quantity Measured in Daily Living

Figure 5 shows the correlation between the amount of daily walking according to the lower-back gold standard and the activity according to the best algorithm (i.e., DCNN). It shows the correlations in the study cohorts recorded under daily living conditions, i.e., the PD and OA groups. Correlations obtained using the other gait detection algorithms are presented in Table 3 and Table 4.

3.4. Gait Quality Measures in Daily Living

Gait quality measures were extracted from the gait bouts that were classified the DCNN model. Group differences in daily living gait quality are summarized in Table 5, Table 6 and Table 7. All the gait quality measures were significantly higher (i.e., better) in the HYAs dataset than in the PD and OA groups as expected, since the PD and OA groups had mobility deficits (recall their TUG times).
In addition, we examined the correlation between the daily quality measures obtained by the gold standard and according to the DCNN model. Table 8 and Table 9 report these correlations.

4. Discussion

In this work, we evaluated algorithms for automatic gait detection from wrist-worn sensor accelerometer data collected from OAs with and without PD in real-world, daily living settings. Data from 12 OAs and 18 PD who wore wrist sensors for up to 10 days were used to train and evaluate an unsupervised anomaly detection algorithm and four previously published gait detection methods. All algorithms were first applied to a dataset of HYAs that was collected in a controlled setting to first evaluate and compare the models in a less challenging cohort and conditions. We then applied these five gait detection algorithms to detect gait during daily living in OAs and people with PD.
As presented in Table 2, all models evaluated obtained high performance when applied to the HYA dataset. This indicates that, on some level, they are suitable for the gait detection task. These results also provide a baseline for the performance of the other models, which enabled a standardized comparison between the models. The TUG times of the OA group and the PD group were higher than 13.5 s (recall Table 1). This value is often used as a cut-off point to identify patients with an abnormal gait and an increased risk of falls [44]. This implies that the OAs and people with PD both had mobility deficits. Thus, it was not surprising that the gait detection results obtained for these two groups were not as good as those obtained for the HYAs. This highlights the challenges of gait detection in populations with impaired gait and movement disorders as well as the challenges that are encountered in detecting gait in real-world, daily living settings using a wrist sensor.
The comparisons between the present results to the previously published findings of the evaluated algorithms were limited and not necessarily direct since some of the performance metrics were lacking in those studies, especially those metrics that were more informative regarding imbalanced data (i.e., precision and recall). For example, Raykov et al. reported only the balanced accuracy (i.e., average of sensitivity and specificity) for gait detection from the wrist sensor, which was higher than the balanced accuracy obtained with our PD dataset. Specifically, we obtained 72.6% in our PD dataset, and they reported 83% with their PD dataset [28]. Their dataset included sensor recordings during uninterrupted and unscripted daily life activities in the participants’ homes and for approximately at least one hour. In our dataset, participants wore the wrist-worn sensor for up to 10 days including during various activities (e.g., driving, sleeping, and cycling) that likely complicated and challenged the gait detection algorithms. In the DCNN paper, only accuracy was reported, and it was 90.22%, which is lower than the accuracy obtained in our experiments (i.e., a 91.7–96.2% accuracy for the different cohorts) [31]. For the LSTM model, the recall and precision reported for walking windows detection were 95.77% and 93.33%, respectively, in healthy adults who wore the sensors on their waist [30]. Several factors can explain the difference in the lower results obtained with our dataset. First, the signal from the wrist sensor is much noisier than the signal from the waist due to the complexities of hand movement. Second, the sensor in their experiment consisted of both accelerometers and gyroscopes which provided more information about the movement and allowed better detection of the activities. Another explanation may be that the classification task was not binary, and the model learned better to distinguish walking from other classes.
The DCNN approach obtained the highest results among the five methods that we evaluated. This might be anticipated since this is a supervised method in which the labels are provided during the training phase. The relatively low precision obtained among all the models reveals the problem of imbalanced data that are embedded in general gait detection algorithms, and this is even more prominent among OAs and people with PD. Further research that addresses the imbalanced data challenge may help to improve the performance of the algorithms in these populations.
The anomaly detection-based algorithm achieved high accuracy (99.4%) and relatively high precision (79.1%) in the HYA dataset. This suggests that the anomaly detection approach is feasible for the gait detection task and has the potential to provide a flexible and robust method that does not require labels. When applied to the other datasets, the anomaly detection method also presented comparable results in terms of accuracy, specificity, and sensitivity to the other methods evaluated in the work; however, the precision was lower. The complexity of both the PD and OA gait data (e.g., due to the unrhythmic gait, reduced arm swing) and the non-gait data (i.e., due to the combination of the relatively high duration of stationary activities together with more active ones) challenges the task of density estimation on which this approach is based. More techniques to improve the ability to learn the complicated density, such as modification to the flow architecture to achieve a better semantic structure of the target data, are needed to further explore the possibility of improving the anomaly detection approach [61].
In conclusion, the results of the present study suggest that detecting gait from a wrist-worn accelerometer can be achieved, albeit with some limitations. The imbalance embedded in a daily living dataset was expressed in our results by a significant trade-off between precision and recall. Hence, for identifying all walking bouts during daily living, additional work is needed. Still, for some purposes, the current approaches may already be adequate. For example, if the goal is to evaluate gait quality, it might be sufficient to use a model with high precision (i.e., relatively few false positives) to correctly identify some of the gait segments, even though not all the segments will be identified (i.e., low recall) (recall Table 5, Table 6 and Table 7). In addition, the relatively high correlation between the predicted and gold-standard daily living walking amounts (Figure 5) suggests that we can use this algorithm to estimate gait quantity measures such as daily walking time. We can utilize the ability to extract daily living gait measures from the wrist sensor for more accurate clinical assessment and as a predictor for health outcomes, such as falls [13,14,15,16,17,18].
Nonetheless, it is important to note that, as mentioned in Section 2.2, we considered only walking bouts with a duration of at least 6 s as “gait”, since shorter duration bouts are difficult to reliably detect, in general, but especially using a wrist-worn device, and since most gait quality measures cannot be applied to very short walking bouts [48]. Moreover, several papers that aimed to detect gait from a wrist-worn device also used this approach [27]. However, many shorter walking bouts might be missed in this approach. Future studies should address this limitation by improving the resolution of the detector. Future work should also include more participants to increase the sample size and allow for a more informed comparison of model performance. Different cohorts should also be included to increase model generalizability and to examine the clinical and health value of these daily living gait measures. In addition, it would be interesting to compare the results based on the machine learning methods described here to other approaches that have been used to extract certain aspects of gait from daily living accelerometer data [26,33]. The relatively high compliance that can be obtained with wrist-worn sensors [62] and the present results set the stage for future studies.

Author Contributions

Conceptualization, Y.E.B., D.S., E.G., R.G.-B. and J.M.H.; Formal analysis, Y.E.B. and D.S.; Funding acquisition, J.M.H.; Methodology, Y.E.B. and D.S.; Software, Y.E.B. and D.S.; Supervision, J.M.H.; Validation, E.G.; Visualization, Y.E.B. and D.S.; Writing—original draft, Y.E.B. and D.S.; Writing—review and editing, Y.E.B., D.S., A.S.B., R.G.-B. and J.M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported, in part, by grants from the Israel Innovation Authority and from the NIH (R01AG017917; R01AG056352). E.G. and J.M.H. were supported, in part, by Mobilise-D. Mobilise-D has received funding from the Innovative Medicines Initiative 2 Joint Undertaking (under grant agreement no. 820820). This joint undertaking received support from the European Union’s Horizon 2020 Research and Innovation Programme and the EFPIA. This work was supported by a grant from the Tel Aviv University Center for AI and Data Science (TAD).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of the Tel Aviv Sourasky Medical Center (approval No. 0648-18-TLV, 1 July 2020).

Informed Consent Statement

Informed written consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author and upon consideration of human studies and Helsinki approvals.

Acknowledgments

We thank the study participants for their time. Thank you to Jordan Raykov for the invaluable discussions regarding unsupervised methods.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Buchman, A.S. Total Daily Physical Activity and Longevity in Old Age. Arch. Intern. Med. 2012, 172, 444. [Google Scholar] [CrossRef] [PubMed]
  2. Buchman, A.S.; Boyle, P.A.; Yu, L.; Shah, R.C.; Wilson, R.S.; Bennett, D.A. Total daily physical activity and the risk of AD and cognitive decline in older adults. Neurology 2012, 78, 1323–1329. [Google Scholar] [CrossRef] [PubMed]
  3. Huisingh-Scheetz, M.J.; Kocherginsky, M.; Magett, E.; Rush, P.; Dale, W.; Waite, L. Relating wrist accelerometry measures to disability in older adults. Arch. Gerontol. Geriatr. 2016, 62, 68–74. [Google Scholar] [CrossRef] [PubMed]
  4. Galperin, I.; Hillel, I.; Del Din, S.; Bekkers, E.M.J.; Nieuwboer, A.; Abbruzzese, G.; Avanzino, L.; Nieuwhof, F.; Bloem, B.R.; Rochester, L.; et al. Associations between daily-living physical activity and laboratory-based assessments of motor severity in patients with falls and Parkinson’s disease. Parkinsonism Relat. Disord. 2019, 62, 85–90. [Google Scholar] [CrossRef]
  5. Mahadevan, N.; Demanuele, C.; Zhang, H.; Volfson, D.; Ho, B.; Erb, M.K.; Patel, S. Development of digital biomarkers for resting tremor and bradykinesia using a wrist-worn wearable device. NPJ Digit. Med. 2020, 3, 1–12. [Google Scholar] [CrossRef]
  6. Powers, R.; Etezadi-Amoli, M.; Arnold, E.M.; Kianian, S.; Mance, I.; Gibiansky, M.; Trietsch, D.; Alvarado, A.S.; Kretlow, J.D.; Herrington, T.M.; et al. Smartwatch inertial sensors continuously monitor real-world motor fluctuations in Parkinson’s disease. Sci. Transl. Med. 2021, 13, eabd7865. [Google Scholar] [CrossRef]
  7. Ekblom, O.; Nyberg, G.; Bak, E.E.; Ekelund, U.; Marcus, C. Validity and Comparability of a Wrist-Worn Accelerometer in Children. J. Phys. Act. Health 2012, 9, 389–393. [Google Scholar] [CrossRef]
  8. Migueles, J.H.; Rowlands, A.V.; Huber, F.; Sabia, S.; van Hees, V.T. GGIR: A Research Community–Driven Open Source R Package for Generating Physical Activity and Sleep Outcomes From Multi-Day Raw Accelerometer Data. J. Meas. Phys. Behav. 2019, 2, 188–196. [Google Scholar] [CrossRef]
  9. Adamowicz, L.; Christakis, Y.; Czech, M.D.; Adamusiak, T. SciKit Digital Health: Python Package for Streamlined Wearable Inertial Sensor Data Processing. JMIR mHealth uHealth 2022, 10, e36762. [Google Scholar] [CrossRef]
  10. Soltani, A.; Dejnabadi, H.; Savary, M.; Aminian, K. Real-World Gait Speed Estimation Using Wrist Sensor: A Personalized Approach. IEEE J. Biomed. Health Inform. 2020, 24, 658–668. [Google Scholar] [CrossRef]
  11. Karas, M.; Urbanek, J.K.; Illiano, V.P.; Bogaarts, G.; Crainiceanu, C.M.; Dorn, J.F. Estimation of free-living walking cadence from wrist-worn sensor accelerometry data and its association with SF-36 quality of life scores. Physiol. Meas. 2021, 42, 065006. [Google Scholar] [CrossRef]
  12. Fasel, B.; Duc, C.; Dadashi, F.; Bardyn, F.; Savary, M.; Farine, P.A.; Aminian, K. A wrist sensor and algorithm to determine instantaneous walking cadence and speed in daily life walking. Med. Biol. Eng. Comput. 2017, 55, 1773–1785. [Google Scholar] [CrossRef]
  13. Studenski, S. Gait Speed and Survival in Older Adults. JAMA 2011, 305, 50. [Google Scholar] [CrossRef]
  14. Perera, S.; Patel, K.V.; Rosano, C.; Rubin, S.M.; Satterfield, S.; Harris, T.; Ensrud, K.; Orwoll, E.; Lee, C.G.; Chandler, J.M.; et al. Gait Speed Predicts Incident Disability: A Pooled Analysis. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2015, 71, 63–71. [Google Scholar] [CrossRef]
  15. Urbanek, J.K.; Roth, D.L.; Karas, M.; Wanigatunga, A.A.; Mitchell, C.M.; Juraschek, S.P.; Cai, Y.; Appel, L.J.; Schrack, J.A. Free-living gait cadence measured by wearable accelerometer: A promising alternative to traditional measures of mobility for assessing fall risk. J. Gerontol. Ser. A 2022, glac013. [Google Scholar] [CrossRef]
  16. Mc Ardle, R.; Del Din, S.; Donaghy, P.; Galna, B.; Thomas, A.J.; Rochester, L. The Impact of Environment on Gait Assessment: Considerations from Real-World Gait Analysis in Dementia Subtypes. Sensors 2021, 21, 813. [Google Scholar] [CrossRef]
  17. Weiss, A.; Herman, T.; Giladi, N.; Hausdorff, J.M. Objective Assessment of Fall Risk in Parkinson’s Disease Using a Body-Fixed Sensor Worn for 3 Days. PLoS ONE 2014, 9, e96675. [Google Scholar] [CrossRef]
  18. Ihlen, E.A.F.; Weiss, A.; Bourke, A.; Helbostad, J.L.; Hausdorff, J.M. The complexity of daily life walking in older adult community-dwelling fallers and non-fallers. J. Biomech. 2016, 49, 1420–1428. [Google Scholar] [CrossRef]
  19. van Schooten, K.S.; Pijnappels, M.; Rispens, S.M.; Elders, P.J.M.; Lips, P.; Daffertshofer, A.; Beek, P.J.; van Dieën, J.H. Daily-Life Gait Quality as Predictor of Falls in Older People: A 1-Year Prospective Cohort Study. PLoS ONE 2016, 11, e0158623. [Google Scholar] [CrossRef]
  20. Weiss, A.; Brozgol, M.; Dorfman, M.; Herman, T.; Shema, S.; Giladi, N.; Hausdorff, J.M. Does the Evaluation of Gait Quality During Daily Life Provide Insight Into Fall Risk? A Novel Approach Using 3-Day Accelerometer Recordings. Neurorehabil. Neural Repair 2013, 27, 742–752. [Google Scholar] [CrossRef]
  21. Hausdorff, J.M.; Hillel, I.; Shustak, S.; Del Din, S.; Bekkers, E.M.J.; Pelosin, E.; Nieuwhof, F.; Rochester, L.; Mirelman, A. Everyday Stepping Quantity and Quality Among Older Adult Fallers With and Without Mild Cognitive Impairment: Initial Evidence for New Motor Markers of Cognitive Deficits? J. Gerontol. A Biol. Sci. Med. Sci. 2018, 73, 1078–1082. [Google Scholar] [CrossRef]
  22. Shema-Shiratzky, S.; Hillel, I.; Mirelman, A.; Regev, K.; Hsieh, K.L.; Karni, A.; Devos, H.; Sosnoff, J.J.; Hausdorff, J.M. A wearable sensor identifies alterations in community ambulation in multiple sclerosis: Contributors to real-world gait quality and physical activity. J. Neurol. 2020, 267, 1912–1921. [Google Scholar] [CrossRef]
  23. Galperin, I.; Herman, T.; Assad, M.; Ganz, N.; Mirelman, A.; Giladi, N.; Hausdorff, J.M. Sensor-Based and Patient-Based Assessment of Daily-Living Physical Activity in People with Parkinson’s Disease: Do Motor Subtypes Play a Role? Sensors 2020, 20, 7015. [Google Scholar] [CrossRef]
  24. Keren, K.; Busse, M.; Fritz, N.E.; Muratori, L.M.; Gazit, E.; Hillel, I.; Scheinowitz, M.; Gurevich, T.; Inbar, N.; Omer, N.; et al. Quantification of Daily-Living Gait Quantity and Quality Using a Wrist-Worn Accelerometer in Huntington’s Disease. Front. Neurol. 2021, 12, 1754. [Google Scholar] [CrossRef]
  25. Femiano, R.; Werner, C.; Wilhelm, M.; Eser, P. Validation of open-source step-counting algorithms for wrist-worn tri-axial accelerometers in cardiovascular patients. Gait Posture 2022, 92, 206–211. [Google Scholar] [CrossRef]
  26. Willetts, M.; Hollowell, S.; Aslett, L.; Holmes, C.; Doherty, A. Statistical machine learning of sleep and physical activity phenotypes from sensor data in 96,220 UK Biobank participants. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef] [PubMed]
  27. Soltani, A.; Paraschiv-Ionescu, A.; Dejnabadi, H.; Marques-Vidal, P.; Aminian, K. Real-World Gait Bout Detection Using a Wrist Sensor: An Unsupervised Real-Life Validation. IEEE Access 2020, 8, 102883–102896. [Google Scholar] [CrossRef]
  28. Raykov, Y.P.; Evers, L.J.W.; Badawy, R.; Bloem, B.R.; Heskes, T.M.; Meinders, M.J.; Claes, K.; Little, M.A. Probabilistic Modelling of Gait for Robust Passive Monitoring in Daily Life. IEEE J. Biomed. Heal. Informatics 2021, 25, 2293–2304. [Google Scholar] [CrossRef]
  29. Vepakomma, P.; De, D.; Das, S.K.; Bhansali, S. A-Wristocracy: Deep learning on wrist-worn sensing for recognition of user complex activities. In Proceedings of the 2015 IEEE 12th International Conference on Wearable and Implantable Body Sensor Networks, BSN 2015, Cambridge, MA, USA, 9–12 June 2015; Institute of Electrical and Electronics Engineers, Inc.: Piscataway, NJ, USA, 2015. [Google Scholar]
  30. Zhao, Y.; Yang, R.; Chevalier, G.; Xu, X.; Zhang, Z. Deep Residual Bidir-LSTM for Human Activity Recognition Using Wearable Sensors. Math. Probl. Eng. 2018, 2018, 1–13. [Google Scholar] [CrossRef]
  31. Zou, Q.; Wang, Y.; Wang, Q.; Zhao, Y.; Li, Q. Deep Learning-Based Gait Recognition Using Smartphones in the Wild. IEEE Trans. Inf. Forensics Secur. 2020, 15, 3197–3212. [Google Scholar] [CrossRef] [Green Version]
  32. Shavit, Y.; Klein, I. Boosting Inertial-Based Human Activity Recognition with Transformers. IEEE Access 2021, 9, 53540–53547. [Google Scholar] [CrossRef]
  33. Karas, M.; Stra Czkiewicz, M.; Fadel, W.; Harezlak, J.; Crainiceanu, C.M.; Urbanek, J.K. Adaptive empirical pattern transformation (ADEPT) with application to walking stride segmentation. Biostatistics 2021, 22, 331–347. [Google Scholar] [CrossRef] [PubMed]
  34. CRAN—Package GENEAclassify. Available online: https://cran.r-project.org/web/packages/GENEAclassify/index.html (accessed on 15 July 2022).
  35. Huang, X.; Mahoney, J.M.; Lewis, M.M.; Guangwei, D.; Piazza, S.J.; Cusumano, J.P. Both coordination and symmetry of arm swing are reduced in Parkinson’s disease. Gait Posture 2012, 35, 373–377. [Google Scholar] [CrossRef]
  36. Seliya, N.; Abdollah Zadeh, A.; Khoshgoftaar, T.M. A literature review on one-class classification and its potential applications in big data. J. Big Data 2021, 8, 1–31. [Google Scholar] [CrossRef]
  37. Bergmann, P.; Fauser, M.; Sattlegger, D.; Steger, C. MVTec AD-A Comprehensive Real-World Dataset for Unsupervised Anomaly Detection. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR), Long Beach, CA, USA, 15–20 June 2019; pp. 9584–9592. [Google Scholar]
  38. Schmidt, M.; Simic, M. Normalizing flows for novelty detection in industrial time series data. In Proceedings of the International Conference on Machine Learning (ICML), Long Beach, CA, USA, 9–15 June 2019. [Google Scholar] [CrossRef]
  39. Yu, J.; Zheng, Y.; Wang, X.; Li, W.; Wu, Y.; Zhao, R.; Wu, L. FastFlow: Unsupervised Anomaly Detection and Localization via 2D Normalizing Flows. In Proceedings of the Conference on Computer Vision and Pattern Recognition (CVPR), Nashville, TN, USA, 19–25 June 2021. [Google Scholar] [CrossRef]
  40. Gudovskiy, D.; Ishizaka, S.; Kozuka, K. Cflow-ad: Real-time unsupervised anomaly detection with localization via conditional normalizing flows. In Proceedings of the IEEE/CVF Winter Conference on Applications Of Computer Vision, Waikoloa, HI, USA, 3–8 January 2021; pp. 98–107. [Google Scholar]
  41. Rudolph, M.; Wandt, B.; Rosenhahn, B. Same Same but DifferNet: Semi-Supervised Defect Detection With Normalizing Flows. In Proceedings of the IEEE/CVF Winter Conference on Applications of Computer Vision (WACV), Waikoloa, HI, USA, 5–9 January 2021; pp. 1907–1916. [Google Scholar]
  42. Karas, M.; Urbanek, J.; Crainiceanu, C.; Harezlak, J.; Fadel, W. Labeled Raw Accelerometry Data Captured during Walking, Stair Climbing and Driving (Version 1.0.0). June 2000. Available online: https://physionet.org/content/accelerometry-walk-climb-drive/1.0.0/ (accessed on 11 September 2022). [CrossRef]
  43. Goldberger, A.L.; Amaral, L.A.N.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.-K.; Stanley, H.E. PhysioBank, PhysioToolkit, and PhysioNet. Circulation 2000, 101, E215–E220. [Google Scholar] [CrossRef]
  44. Shumway-Cook, A.; Brauer, S.; Woollacott, M. Predicting the Probability for Falls in Community-Dwelling Older Adults Using the Timed Up & Go Test. Phys. Ther. 2000, 80, 896–903. [Google Scholar] [CrossRef]
  45. Hillel, I.; Gazit, E.; Nieuwboer, A.; Avanzino, L.; Rochester, L.; Cereatti, A.; Della Croce, U.; Rikkert, M.O.; Bloem, B.R.; Pelosin, E.; et al. Is every-day walking in older adults more analogous to dual-task walking or to usual walking? Elucidating the gaps between gait performance in the lab and during 24/7 monitoring. Eur. Rev. Aging Phys. Act. 2019, 16. [Google Scholar] [CrossRef]
  46. Iluz, T.; Gazit, E.; Herman, T.; Sprecher, E.; Brozgol, M.; Giladi, N.; Mirelman, A.; Hausdorff, J.M. Automated detection of missteps during community ambulation in patients with Parkinson’s disease: A new approach for quantifying fall risk in the community setting. J. Neuroeng. Rehabil. 2014, 11, 1–9. [Google Scholar] [CrossRef]
  47. Czech, M.D.; Patel, S. GaitPy: An Open-Source Python Package for Gait Analysis Using an Accelerometer on the Lower Back. J. Open Source Softw. 2019, 4, 1778. [Google Scholar] [CrossRef]
  48. Bao, L.; Intille, S.S. Activity recognition from user-annotated acceleration data. Lect. Notes Comput. Sci. 2004, 3001, 1–17. [Google Scholar] [CrossRef]
  49. Papamakarios, G.; Pavlakou, T.; Murray, I. Masked Autoregressive Flow for Density Estimation. NIPS. 2017, pp. 2338–2347. Available online: http://papers.nips.cc/paper/6828-masked-autoregressive-flow-for-density-estimation.pdf (accessed on 27 December 2019).
  50. Germain, M.; Gregor, K.; Murray, I.; Larochelle, H. MADE: Masked Autoencoder for Distribution Estimation. In Proceedings of the 32nd International Conference on Machine Learning, Lille, France, 6–11 July 2015; Bach, F., Blei, D., Eds.; PMLR: Lille, France, 2015; Volume 37, pp. 881–889. [Google Scholar]
  51. Ronneberger, O.; Fischer, P.; Brox, T. U-net: Convolutional networks for biomedical image segmentation. Lect. Notes Comput. Sci. 2015, 9351, 234–241. [Google Scholar] [CrossRef]
  52. Ioffe, S.; Szegedy, C. Batch normalization: Accelerating deep network training by reducing internal covariate shift. In Proceedings of the 32nd International Conference on Machine Learning, Lille, France, 11 February 2015. [Google Scholar]
  53. Srivastava, N.; Hinton, G.; Krizhevsky, A.; Salakhutdinov, R. Dropout: A Simple Way to Prevent Neural Networks from Overfitting. J. Mach. Learn. Res. 2014, 15, 1929–1958. [Google Scholar]
  54. Yu, Y.; Si, X.; Hu, C.; Zhang, J. A Review of Recurrent Neural Networks: LSTM Cells and Network Architectures. Neural Comput. 2019, 31, 1235–1270. [Google Scholar] [CrossRef]
  55. Zarzycki, K.; Ławryńczuk, M.; Ławryńczuk, Ł.; Lenci, S. LSTM and GRU Neural Networks as Models of Dynamical Processes Used in Predictive Control: A Comparison of Models Developed for Two Chemical Reactors. Sensors 2021, 21, 5625. [Google Scholar] [CrossRef]
  56. Wu, H.; Zhang, J.; Zong, C. An Empirical Exploration of Skip Connections for Sequential Tagging. arXiv Prepr. 2016, 203–212. [Google Scholar] [CrossRef]
  57. Saito, T.; Rehmsmeier, M. The Precision-Recall Plot Is More Informative than the ROC Plot When Evaluating Binary Classifiers on Imbalanced Datasets. PLoS ONE 2015, 10, e0118432. [Google Scholar] [CrossRef]
  58. Weiss, A.; Sharifi, S.; Plotnik, M.; Van Vugt, J.P.P.; Giladi, N.; Hausdorff, J.M. Toward automated, at-home assessment of mobility among patients with Parkinson disease, using a body-worn accelerometer. Neurorehabil. Neural Repair 2011, 25, 810–818. [Google Scholar] [CrossRef]
  59. Huijben, B.; van Schooten, K.S.; van Dieën, J.H.; Pijnappels, M. The effect of walking speed on quality of gait in older adults. Gait Posture 2018, 65, 112–116. [Google Scholar] [CrossRef] [Green Version]
  60. Moe-Nilssen, R.; Helbostad, J.L. Estimation of gait cycle characteristics by trunk accelerometry. J. Biomech. 2004, 37, 121–126. [Google Scholar] [CrossRef]
  61. Kirichenko, P.; Izmailov, P.; Wilson, A.G. Why Normalizing Flows Fail to Detect Out-of-Distribution Data. In Advances in Neural Information Processing Systems; Larochelle, H., Ranzato, M., Hadsell, R., Balcan, M.F., Lin, H., Eds.; Curran Associates, Inc.: Nice, France, 2020; Volume 33, pp. 20578–20589. [Google Scholar]
  62. Burq, M.; Rainaldi, E.; Ho, K.C.; Chen, C.; Bloem, B.R.; Evers, L.J.W.; Helmich, R.C.; Myers, L.; Marks, W.J.; Kapur, R. Virtual exam for Parkinson’s disease enables frequent and reliable remote measurements of motor function. npj Digit. Med. 2022, 5, 1–9. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flowchart based on Keren’s gait detection algorithm [24]. The raw data were divided into six-second windows. The Euclidean norm of the 3-axis accelerometer signal was calculated to minimize any effects of the sensor’s placement and angle or orientation dependency. An activity threshold was set to 1 ( m / s 2 ) to remove areas in the recording with low activity. Then, each window was classified as a walking bout based on the time and frequency domain features.
Figure 1. Flowchart based on Keren’s gait detection algorithm [24]. The raw data were divided into six-second windows. The Euclidean norm of the 3-axis accelerometer signal was calculated to minimize any effects of the sensor’s placement and angle or orientation dependency. An activity threshold was set to 1 ( m / s 2 ) to remove areas in the recording with low activity. Then, each window was classified as a walking bout based on the time and frequency domain features.
Sensors 22 07094 g001
Figure 2. The architecture of the DCNN model. Adapted from Zou et al. (2020) [31].
Figure 2. The architecture of the DCNN model. Adapted from Zou et al. (2020) [31].
Sensors 22 07094 g002
Figure 3. LSTM cell diagram. Adapted from Zarzycky and Ławryńczuk (2021) [55].
Figure 3. LSTM cell diagram. Adapted from Zarzycky and Ławryńczuk (2021) [55].
Sensors 22 07094 g003
Figure 4. PR (left) and ROC curves (right) illustrate the performance of the four gait detection models with the three study cohorts: PD, OAs, and HYAs: (a) PR curve (left) and ROC curve (right) of the AR-iHMM model; (b) PR curve (left) and ROC curve (right) of the DCNN model; (c) PR curve (left) and ROC curve (right) of the Bidir-LSTM model; (d) PR curve (left) and ROC curve (right) of the anomaly detection model.
Figure 4. PR (left) and ROC curves (right) illustrate the performance of the four gait detection models with the three study cohorts: PD, OAs, and HYAs: (a) PR curve (left) and ROC curve (right) of the AR-iHMM model; (b) PR curve (left) and ROC curve (right) of the DCNN model; (c) PR curve (left) and ROC curve (right) of the Bidir-LSTM model; (d) PR curve (left) and ROC curve (right) of the anomaly detection model.
Sensors 22 07094 g004aSensors 22 07094 g004b
Figure 5. Association between the gold standard and the DCNN model output. Each point represents the sum of the activity for each day. (a) Pearson’s correlation coefficient = 0.83; p-value < 0.001. (b) OAs: Pearson’s correlation coefficient = 0.95; p-value < 0.001. When removing the outlier point in the upper right corner, the Pearson’s correlation coefficient changed to 0.9.
Figure 5. Association between the gold standard and the DCNN model output. Each point represents the sum of the activity for each day. (a) Pearson’s correlation coefficient = 0.83; p-value < 0.001. (b) OAs: Pearson’s correlation coefficient = 0.95; p-value < 0.001. When removing the outlier point in the upper right corner, the Pearson’s correlation coefficient changed to 0.9.
Sensors 22 07094 g005
Table 1. Characteristics of the older adults analyzed in this study.
Table 1. Characteristics of the older adults analyzed in this study.
No Parkinson’s Disease
(No PD)
Parkinson’s Disease
(PD)
No. of Subjects (N)1218
Age (years)76.1 ± 7.372.7 ± 8.4
Gender (M:F)3:910:8
Education (years)15.8 ± 3.415.9 ± 3.1
Montreal Cognitive Assessment (0–30)24.3 ± 5.322.4 ± 5.3
Hoehn and Yahr Stage (1–4)-2.3 ± 0.8
Disease Duration in Years-5.56 ± 4.05
Walking aid (No. of Subjects)3 (25%)4 (22%)
Timed Up and Go (s)14.8 ± 3.515.1 ± 8.8
Dual-Task Timed Up and Go (s)17.5 ± 5.119.5 ± 10.7
Gait Speed (cm/s)99 ± 20 110 ± 27
Days Sensor Worn6.0 ± 2.48.1 ± 2.3
No. of Gait Bouts per Recording864 ±5481264 ± 841
Average Length of Gait Bouts (s)19.5 ± 4.120.4 ± 5.4
Gait Percent (%) per Recording7.2 ± 2.66.8 ± 4.7
Table 2. Comparing a new anomaly detection algorithm with prior methods. The algorithm’s performance in gait detection (values in parentheses indicate the STD between the different validation folds).
Table 2. Comparing a new anomaly detection algorithm with prior methods. The algorithm’s performance in gait detection (values in parentheses indicate the STD between the different validation folds).
Accuracy (%)Precision (%)Sensitivity (%)Specificity (%)
HYAs
FB71.036.377.869.5
AR-iHMM73.240.094.668.4
DCNN96.090.983.898.4
Bidir-LSTM90.970.684.392.3
New Anomaly Detection99.479.144.699.4
OAs
FB76.2 (2.4) 15.4 (5.8)59.8 (6.3)77.3 (2.5)
AR-iHMM85.2 (8.1)27.3 (8.6)78.7 (17.9)85 (9.43)
DCNN96.2 (4.2)61.4 (3.4)69.8 (9.8)97.6 (4.9)
Bidir-LSTM81.7 (6.6)38.8 (11.9)81.7 (5.7)81.8 (7.4)
New Anomaly Detection76.2 (0.5)29.9 (0.5)67.6 (0.1)78.1 (0.6)
PD
FB77.6 (5.7)15.2 (11.3)57.4 (15.6)57.4 (5.5)
AR-iHMM92.6 (4.7)38.1 (24.0)50 (23.3)95.1 (5.0)
DCNN91.7 (7.4)40.9 (6.0)72.9 (7.4)92.9 (9.2)
Bidir-LSTM76.6 (4.6)38.3 (7.1)88.6 (10.9)76.5 (6.1)
New Anomaly Detection83.3 (3.8)20.9 (2.3)74.7 (8.2)83.8 (5.4)
Table 3. Correlation between the gold standard and the gait detection models in the OAs cohort.
Table 3. Correlation between the gold standard and the gait detection models in the OAs cohort.
ModelPearson’s Correlation Coefficientp-Value
FB0.80 < 0.001   *
AR-iHMM0.78 < 0.001   *
DCNN0.95 < 0.001   *
Bidir-LSTM0.92 < 0.001   *
New Anomaly Detection0.84 < 0.001   *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Table 4. Correlation between the gold standard and the gait detection models in the PD cohort.
Table 4. Correlation between the gold standard and the gait detection models in the PD cohort.
ModelPearson’s Correlation Coefficientp-Value
FB0.37 < 0.001   *
AR-iHMM0.57 < 0.001   *
DCNN0.85 < 0.001   *
Bidir-LSTM0.88 < 0.001   *
New Anomaly Detection0.67 < 0.001   *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Table 5. Daily living gait quality measures in OAs and healthy adults (HYAs).
Table 5. Daily living gait quality measures in OAs and healthy adults (HYAs).
Daily Living Gait Quality MeasureOAsHYAsp-Value
RhythmDominant Frequency (Hz)0.86 ± 0.21.53 ± 0.49 < 0.001   *
MagnitudeRMS ( m / s 2 )1.36 ± 0.253.36 ± 0.73 < 0.001   *
Regularity/ConsistencyStride regularity (unitless)0.13 ± 0.020.47 ± 0.23 < 0.001   *
Amplitude dominant Frequency (unitless)0.43 ± 0.070.73 ± 0.22 < 0.001   *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Table 6. Daily living gait quality measures in people with PD and healthy adults (HYAs).
Table 6. Daily living gait quality measures in people with PD and healthy adults (HYAs).
Daily Living Gait Quality MeasurePDHYAsp-Value
RhythmDominant Frequency (Hz)1.21 ± 0.451.53 ± 0.49 < 0.001 *
MagnitudeRMS ( m / s 2 )1.16 ± 0.593.36 ± 0.73 < 0.001 *
Regularity/ConsistencyStride Regularity (unitless)0.14 ± 0.050.47 ± 0.23 < 0.001 *
Amplitude Dominant Frequency (unitless)0.33 ± 0.150.73 ± 0.22 < 0.001 *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Table 7. Daily living gait quality measures in OAs compared to people with PD.
Table 7. Daily living gait quality measures in OAs compared to people with PD.
Daily Living Gait Quality MeasureOAsPDp-Value
RhythmDominant Frequency (Hz)0.86 ± 0.21.21 ± 0.45 < 0.001 *
MagnitudeRMS ( m / s 2 )1.36 ± 0.251.16 ± 0.59 0.0293
Regularity/ConsistencyStride Regularity (unitless)0.13 ± 0.020.14 ± 0.05 0.1094
Amplitude Dominant Frequency (unitless)0.43 ± 0.070.33 ± 0.15 < 0.001 *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) . Entries are reported as the mean ± SD determined for each participant across all walking bouts larger than 10 s. In the HYA dataset, the values used for comparison were the measures from all the walking bouts. In the OA and the PD datasets, we calculated the gait quality measures for each day’s walking bouts. Then, we took the median of these measures for comparison with the HYA dataset.
Table 8. Correlations between the gold standard and the DCNN model in the OA cohort.
Table 8. Correlations between the gold standard and the DCNN model in the OA cohort.
Daily Living Gait Quality MeasurePearson’s Correlation Coefficientp-Value
Dominant Frequency (Hz)0.48 < 0.001   *
RMS ( m / s 2 )0.70 < 0.001   *
Stride Regularity (unitless)0.48 < 0.001   *
Amplitude Dominant Frequency (unitless)0.49 < 0.001   *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Table 9. Correlations between the gold standard and the DCNN model in the PD cohort.
Table 9. Correlations between the gold standard and the DCNN model in the PD cohort.
Daily Living Gait Quality MeasurePearson’s Correlation Coefficientp-Value
Dominant Frequency (Hz)0.46 < 0.001   *
RMS ( m / s 2 )0.48 < 0.001   *
Stride Regularity (unitless)0.34 0.016
Amplitude Dominant Frequency (unitless)0.44 0.011   *
* Significant after Bonferroni correction ( p = α c o m p o n e n t s ;   α = 0.05 ) .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brand, Y.E.; Schwartz, D.; Gazit, E.; Buchman, A.S.; Gilad-Bachrach, R.; Hausdorff, J.M. Gait Detection from a Wrist-Worn Sensor Using Machine Learning Methods: A Daily Living Study in Older Adults and People with Parkinson’s Disease. Sensors 2022, 22, 7094. https://0-doi-org.brum.beds.ac.uk/10.3390/s22187094

AMA Style

Brand YE, Schwartz D, Gazit E, Buchman AS, Gilad-Bachrach R, Hausdorff JM. Gait Detection from a Wrist-Worn Sensor Using Machine Learning Methods: A Daily Living Study in Older Adults and People with Parkinson’s Disease. Sensors. 2022; 22(18):7094. https://0-doi-org.brum.beds.ac.uk/10.3390/s22187094

Chicago/Turabian Style

Brand, Yonatan E., Dafna Schwartz, Eran Gazit, Aron S. Buchman, Ran Gilad-Bachrach, and Jeffrey M. Hausdorff. 2022. "Gait Detection from a Wrist-Worn Sensor Using Machine Learning Methods: A Daily Living Study in Older Adults and People with Parkinson’s Disease" Sensors 22, no. 18: 7094. https://0-doi-org.brum.beds.ac.uk/10.3390/s22187094

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop