Next Article in Journal
Non-Destructive Investigation of Intrinsic Magnetic Field of Austenitic Biomaterials by Magnetic Field Sensors
Next Article in Special Issue
Channel Reduction for an EEG-Based Authentication System While Performing Motor Movements
Previous Article in Journal
Detection of Image Level Forgery with Various Constraints Using DFDC Full and Sample Datasets
Previous Article in Special Issue
Test-Retest Reliability of Acoustic Emission Sensing of the Knee during Physical Tasks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection of Parkinson’s Disease Using Wrist Accelerometer Data and Passive Monitoring

1
Department of Business Intelligence and Analytics, Business College, Creighton University, Omaha, NE 68178, USA
2
Department of Biomedical Informatics, College of Information Systems and Technology, University of Nebraska at Omaha, Omaha, NE 68182, USA
3
Department of Biomechanics, College of Education, Health and Human Sciences, University of Nebraska at Omaha, Omaha, NE 68182, USA
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(23), 9122; https://doi.org/10.3390/s22239122
Submission received: 1 October 2022 / Revised: 11 November 2022 / Accepted: 18 November 2022 / Published: 24 November 2022
(This article belongs to the Collection Sensors for Gait, Posture, and Health Monitoring)

Abstract

:
Parkinson’s disease is a neurodegenerative disorder impacting patients’ movement, causing a variety of movement abnormalities. It has been the focus of research studies for early detection based on wearable technologies. The benefit of wearable technologies in the domain rises by continuous monitoring of this population’s movement patterns over time. The ubiquity of wrist-worn accelerometry and the fact that the wrist is the most common and acceptable body location to wear the accelerometer for continuous monitoring suggests that wrist-worn accelerometers are the best choice for early detection of the disease and also tracking the severity of it over time. In this study, we use a dataset consisting of one-week wrist-worn accelerometry data collected from individuals with Parkinson’s disease and healthy elderlies for early detection of the disease. Two feature engineering methods, including epoch-based statistical feature engineering and the document-of-words method, were used. Using various machine learning classifiers, the impact of different windowing strategies, using the document-of-words method versus the statistical method, and the amount of data in terms of number of days were investigated. Based on our results, PD was detected with the highest average accuracy value (85% ± 15%) across 100 runs of SVM classifier using a set of features containing features from every and all windowing strategies. We also found that the document-of-words method significantly improves the classification performance compared to the statistical feature engineering model. Although the best performance of the classification task between PD and healthy elderlies was obtained using seven days of data collection, the results indicated that with three days of data collection, we can reach a classification performance that is not significantly different from a model built using seven days of data collection.

1. Introduction

Parkinson’s disease (PD) is a neurological disorder that occurs due to a decrease in the brain dopamine level, impacting patients’ movement. Symptoms of PD include movement abnormalities such as tremor, dysfunction of magnitude, and symmetry of arm swing, and gait characterized by shorter strides and slower movements compared to healthy individuals [1,2]. These symptoms have been shown to be detectable with a good amount of accuracy using body sensor networks and wearable technologies including wearable accelerometers [3,4]. The benefit of wearable technologies in the domain rises by continuous monitoring of this population’s movement patterns over time. Continuous monitoring seems to be a necessary step toward real-life applications of movement analysis in the context of the smart home and healthcare domain [5,6]. Monitoring patients in their own home over time provides opportunities for early detection of the disease and useful insights into the disease evolution, which in turn helps healthcare providers to assess the impact of various medications and interventions and come up with effective ways of slowing down the disease progression.
The wrist is the most common and acceptable body location for people to wear sensors on [7]. Advanced technologies such as smartwatches and activity trackers, which are being used by many people around the world, have provided us a platform to be able to monitor individuals’ wrist movement patterns continuously and remotely over time. Analyzing wrist movement patterns has been shown to be helpful in detection of freezing of gait events (FoG) and detecting tremor and dyskinesia in the PD population [7,8,9,10]. Early detection and tracking the evolution of the disease over time could be feasible using wrist accelerometers, as suggested by earlier studies [11,12,13]. In [12], authors reported a reduction in the amount and intensity of activity in patients with PD compared to healthy individuals using data from a sports watch. Wrist accelerometers have been shown to provide effective features for identifying PD patients from their healthy counterparts as well as the severity of the disease [13]. Early detection of PD has been shown to be feasible using wrist-worn accelerometry while passive monitoring the participants over a week [11].
The contributions of this study are listed as follows: (1) Early detection of PD using wrist accelerometer and passive monitoring, (2) using the document-of-word feature engineering approach to improve the performance over the traditional statistical feature engineering method, (3) identifying the best window size for segmentation of the wrist accelerometer data so that it increases the accuracy of the early detection of PD, and (4) identifying the optimal amount of monitoring time for the early detection of PD. Each of these goals that lead to our contributions in this paper are explained in the following paragraphs.
Passive monitoring of individuals over time using an easy-to-wear and comfortable sensor on the wrist could provide us with initial information on whether a person is likely to be in the early stages of a disease with movement disorder. This is the case for patients with PD, because PD is a progressive neurodegenerative disorder of the brain and central nervous system, affecting patients’ motor system, causing reduced movement, tremor, postural instability, and postural rigidity [14]. Fatigue, small shuffling steps, freezing of gait, and bradykinesia are some of the symptoms of PD [15]. If passive monitoring indicates any sign of the disease, then further investigation could be performed in the next step [11]. There are many research studies on detecting particular features of PD such as tremor, freezing of gait, and bradykinesia [3,16]. There is also a good amount of research studies on the early detection of PD in a control environment [17,18]. The only research study on the early detection of PD in an uncontrolled environment using wrist accelerometers is published by Williamson et al. [11], in which they used the Biobank dataset [19] collected from the participants’ dominant wrist. Our work will be among the very few studies using wrist accelerometers for the early detection of the disease. Plus, as suggested in [11], we should be able to rely on the early detection of PD using wearable technology regardless of which side of the body participants wear the wrist accelerometer (the most vs. less affected side and dominant vs. non-dominant side). While wearing an accelerometer on participants’ dominant wrist has been shown to provide good results [11], no study has reported on the early detection of PD by passive monitoring of participants’ non-dominant wrist. Moreover, using and combining symptom-based features with the other wrist movement features that are not based on the disease symptom could improve the sensitivity and specificity of the classification technique [11]. Therefore, our first contribution is the early detection of PD using wrist sensor data collected from the non-dominant wrist while combining symptom-based features with the other features.
The standard statistical feature engineering method has been the most common feature extraction method in this domain [20,21]. In this method, statistical features (e.g., vector magnitude) are extracted from signal subsequences of usually the same size. The document-of-words approach is a more recent approach in the domain of wearable signal processing. This approach has been shown to outperform the standard statistical approach in related domains such as real-life activity recognition, using data collected from accelerometers [22], and the early detection of PD, using accelerometers on participants’ left and right ankles [23]. Although the document-of-words feature engineering approach has been shown to provide good performance in PD detection using data collected from participants’ ankles, it was performed on a relatively small dataset with twenty participants. Moreover, it has never been used on accelerometers data collected from the participants’ wrist. Both approaches will be implemented and applied to our data. Then, the results will be compared to determine which approach has the highest performance.
Our third contribution is to test different window size values and identify the best one for data segmentation so that it leads to a classification model with the best performance for the early detection of PD. Choosing the right window size from which features are calculated has been shown to have a significant impact on detecting bradykinesia in patients with PD [3]. Choosing the right window size/sizes leads to a better identification of PD symptoms [3,4], which could in turn lead to a better model for the early detection of PD. Although the effect of feature window size has been studied for detecting PD symptoms, it has not been studied for the early detection of PD using wrist accelerometers, while it matters more in the context of passive monitoring in which there is no label for symptoms. In this study, we investigate the impact of different feature window sizes, when splitting the data, on our diagnostic models’ accuracy and efficiency.
The last contribution of this paper is to provide information on what could be the optimal amount of data needed for the early detection of PD when passive monitoring participants using a wrist accelerometer. The inertial sensors’ battery life depends on factors such as the number of sensors utilized, sampling frequency, and recording time, and even in the best-case scenario, when just accelerometer is utilized, there is a constraint on long-term data recording due to these devices’ limited battery lives [24]. Therefore, identifying the best possible amount of data needed for the early detection of PD without compromising the accuracy is a must-do task in this area. Compared to the time scale of hours, the time scale of days for monitoring PD patients has been found to better correlate with PD clinical metrics [25,26]. The total accuracy of a classification model for the classification of PD and healthy counterparts has been shown to improve gradually over a one full week time course [11]. The question would be when to stop collecting data without compromising the accuracy. In this study, we answer the question of how many days of passive data collection would provide us the best possible model for the early detection of patients with PD.

2. Related Works

Our main aim in this study is the early detection of PD in a free-living environment using wrist accelerometers. Many research studies have been conducted on PD clinical test automation using wearable technologies in the clinical setting [17,27,28,29]. A great deal of research focuses on automatic detection of PD symptoms, e.g., bradykinesia, tremor, freezing of gait, and rigidity in controlled environments such as clinics and biomechanics labs, using wearable devices [16,30,31]. To help PD patients recover from freezing of gait episodes, Gokul H. et al. proposed an embedded machine learning system from accelerometer data [32]. Work also has been performed for continuous monitoring of PD treatments such as levodopa [33]. Using sensors in a controlled environment, researchers have characterized wrist movement by variability measures in PD patients [34]. Using inertial sensors integrated in a smartphone combined with machine learning classification methods researchers could identify patients with PD [35]. Early detection of patients with PD is performed by researchers in controlled environments and mainly using gait features [36].
Our second aim is to investigate the impact of using two feature engineering methods on the early detection of PD: the statistical and document-of-words feature engineering methods. There are only a few research studies in which document-of-words feature engineering was performed. Kheirkhahan et al. [37] used the latter method in their work to identify activities of daily living and assess the energy expenditure per activity using wrist-worn accelerometers. The document-of-words feature engineering method was used by Rastegari et al. [23] to distinguish between PD and healthy elderlies using accelerometer data from participants’ ankles. Both studies showed a good outcome of the document-of-words feature engineering method. However, both used data collected in a controlled environment.
For related work, we consider studies in which a wrist sensor was used to collect data on individuals in free-living conditions, because the mean value of features calculated from data captured in a free-living environment has been shown to be different from that captured in a controlled environment [38,39].
Monitoring PD patients using wearable sensors and in a free-living environment has been attracting attention recently [40,41]. Monitoring in free-living conditions serves the following purposes for the PD population: early detection of the disease, identifying the disease severity level, and assessing the impact of a treatment or an intervention. In this study, our focus is passive monitoring of PD. While in active monitoring individuals are provided with a set of instructions on which tasks to perform and how to perform them, in passive monitoring individuals perform their activities of daily living without being instructed. Only three of the papers discussed in this section focus on the early detection of PD in free-living conditions using wearable technologies.
In a research study by San Segundo et al., wrist accelerometry was used to detect tremor in PD both using annotated laboratory data and free-living data [41]. Their result indicated that tremor in PD can be detected with high performance using laboratory data and not using the free-living data, mainly because of the weak self-reported labels. However, the percentage of time detected with tremor was in good agreement with the self-reported labels. The result of this study cannot be compared with ours because it is not about the early detection of PD, and the performance could not be reported for the free-living environment.
In a study conducted by Lipsmeier et al., individuals with PD and their healthy age-matched counterparts were monitored in free-living conditions using smartphones to investigate differences in features among participants [27]. Active (rest tremor, finger tapping, balance, and gait) and passive monitoring of individuals resulted in significant differences between all active and passive features from PD and controlled participants. Although they achieved 75% sensitivity and 81% specificity using the passive monitoring features, these results were not obtained using an independent test set.
Zhan et al. conducted a longitudinal study on PD subjects over a period of 6 months in which data were collected actively through five smartphone activities and used to develop a mobile PD score [17]. The score developed by them corelated well with several clinical tests including the Hoehn and Yahr scale and UPDRS (Unified PD Rating Scale) score. The result of this study cannot be compared to ours because they reported the coefficient of determination and not any accuracy, sensitivity, or specificity values.
The authors of the Parkinson@Home validation study [42] used accelerometers and gyrometers at five participants’ body locations, including the wrist, to distinguish between PD and healthy individuals. They reported accuracy and AUC values for the most and least affected wrist while the most affected wrist provided a better performance (AUC = 0.75 and accuracy = 57%). We use this study and its results to be compared with ours.
Habets J.G. conducted a study to monitor bradykinesia in patients with PD using a wrist accelerometer [3]. Their data were taken from the Parkinson@Home validation study, and they tried different windowing strategies to investigate the impact of those in the automatic identification of bradykinesia. They also compared group models (in which all subjects’ data were used) with individual models (in which each subject’s data were used separately and then the results averaged). Their results indicate that both group models and individual models perform at the same level (AUC = 0.70), which is the best performance obtained by increasing the window size to 300 s. This study can be used to compare our results with regards to different windowing strategies and with theirs.
The study of Williamson et al. on PD [11] is the most similar work to ours. The authors used the Biobank dataset [20], collected from the participants’ dominant wrist, and a Gaussian mixture model for the early detection of PD. In this study, analysis was performed based on two types of segments: gait segments and low-movement segments. They used gait dispersion features, low-movement correlation structure features, gait segments incidents, and the incidence of low-movement segments. They achieved good results with an area under the curve (AUC) of 0.85. The impact of different amounts of data in terms of number of days on model performance was also investigated in this study. Therefore, this study is a good comparison point with ours.

3. Materials and Methods

The following five subsections of this section explain this study data collection’s protocol and methodology. Figure 1 shows the steps taken to identify the best feature engineering method and the best windowing strategy/strategies. Figure 2 indicates the steps taken to identify the optimal amount of data that need to be collected in terms of days.

3.1. Data Collection

The data were collected from two groups of subjects: (1) healthy elderlies (64.2 ± 7 years and no movement deficiencies) and (2) individuals with various stages of PD (71 ± 6.2). PD patients reported no other illnesses impacting their gait and movement. The Hoehn and Yahr (H&Y) score used in our study is one of the commonly used measures of progression of PD and has been shown to correlate with motor decline happening in PD patients. The data were collected from participants while they had an Actigraph triaxial accelerometer on their non-dominant wrist for seven consecutive days. To be more specific, we used ActiGraph GT3X with ActiLife v6.13.3 Firmware v1.7.1 at 100 Hz. The sensor could be taken off while taking a shower or swimming.
The dataset includes three phases of data collection with 6–12 months between each two consecutive phases. The Actigraph device was recorded at 100 Hz with a range of ±8 g. In each phase, the H&R score and relevant clinical measures were recorded for each participant. The data from the first phase of the data collection is used for our analysis in this paper. Table 1 shows the participants’ characteristics. Figure 3 and Figure 4 show the first five minutes of raw data collected from the wrist of an individual with PD and a healthy individual, respectively.

3.2. Signal Segmentation

In this study, an epoch-based segmentation approach was used in which we converted each subject’s wrist movement signal into subsequences of a predefined length. Since a passive data collection was performed without labeling the data, each segment may represent both activities of daily living and symptoms of PD in the case of PD patients, which should be informative when it comes to distinguishing between the PD population and their healthy counterparts. Feature window size has been shown to significantly impact identifying symptoms of PD, including bradykinesia, resting tremor, and freezing of gait [3,7,16,43,44]. Therefore, different window sizes that are representing different signal granularity levels and could represent different symptoms of the disease were considered in our analysis. The classification task between PD and healthy elderlies was performed using the features obtained from each window size and also different combinations of different window size strategies.
The window lengths considered in this study were as follows: (S ∈ {3 s, 10 s, 60 s, 300 s, and 900 s}). Each individual time series was split into smaller subsequences using a sliding window of S—second with overlapping windows of size S 2 seconds. These window lengths were chosen because they are likely to capture symptoms of PD. Three-seconds window length has been shown to be effective in capturing bradykinesia, tremor, and FOG using wrist accelerometer data [7,16], while increasing the window size to 300 s improves the average AUC of bradykinesia detection in PD patients [3]. In addition, using 5 min (300 s) epochs, researchers found that UPDRS III score is highly corelated with wrist accelerometer metrics [45]. A study of wrist accelerometer data for early detection of PD from a free-living setting used a sliding window of size 10 s and achieved a good AUC value [11]. A window length of 60 s has been commonly used in PD studies with the purpose of symptom detection [3,46]. We also consider a window length of 900 s (15 min) in our analysis, since investigating a window size larger than 300 s has been mentioned as a limitation of other PD studies due to their small sample size [3], which is not the case in our study. Moreover, bigger segments could be representative of activities of daily living.

3.3. Feature Extraction

Using the epoch-based segmentation approach explained in Section 3.2, we used subsequences of data to extract features using two methods: standard statistical feature engineering and document-of-words method. To make the features independent from the accelerometer’s orientation, they were all calculated either over all three axes or the signal vector magnitude. All features were calculated over the subsequences of each windowing strategy.

3.3.1. Epoch-Based Statistical Features

To be able to compare the statistical feature engineering approach with the document-of-words method, a set of statistical features was used to create the epoch-based statistical features. Below is a list of features included in the analysis:
Signal vector magnitude: For the segments of each segmentation strategy, we calculated signal vector magnitude as √(X2 + Y2+ Z2) in which x, y, and z represent the accelerometers’ axes. This measure is dependent from the accelerometer’s orientation. This feature is used by researchers both for identifying PD symptoms [3] and early diagnosis of the disease [11] to identify a segment with hand movement versus a segment with low hand movement. Therefore, we assumed that this feature would help us in differentiating PD patients from their healthy counterparts.
RMS of acceleration in the anterior posterior (AP), mediolateral (ML), and vertical (Ver) directions: Root mean square (RMS) of acceleration is a time domain feature used in many studies focusing on PD symptom detection or early diagnosis of the disease [3,16,44]. Using the following equation (Equation (1)), we normalized the RMS values in each direction:
R M S R d i r = R M S d i r R M S A P 2 + R M S M L 2 + R M S V e r 2
The denominator is signal vector magnitude and dir represents each of the following directions: {AP, ML, Ver}.
Movement dispersion feature: this feature is introduced in [11] and only calculated on the segments identified as gait segments, not on the ones identified as low-movement segments. In this work, dispersion feature was calculated for all segments. Movement dispersion feature quantifies the average amount of acceleration variability within each segment. To calculate the dispersion feature vector for each segment, we used the method introduced by Williamson et al. [11]. First, we converted the acceleration values of each movement segment into their z-score values (zi(t) represents z-score of data point, t, in axis i within the nth segment). Then, we removed the outliers, defined as the points for which the standard deviation of either of the axes is greater than two. Then, using the valid points (non-outlier points) of each segment, dispersion in axis i can be calculated as follows in Equations (2) and (3):
D i s p e r s i o n i ( n ) =   | Zi ( t 1 ) Zi ( t 2 ) | S ( n ) :       t 1 ,   t 2 V ( n )
S ( n ) = t 1 t 2 v n 1  
where V(n) represents a set of valid points within the nth segment, Zi(t) represents the z-score of data point t in axis I, and S(n) is the number of the pairs of valid points within the nth frame. Therefore, the vector of dispersion features for each segment looks like {Dx(n), Dy(n), Dz(n)}.
Acceleration range: This feature is reported among the features used for identifying tremor in PD patients [16,23]. To calculate this set of features, we took the difference between the minimum and maximum signal amplitude from three axes (ARAP, ARML, and ARVer).
The average of these features across all subsequences was calculated in the statistical feature engineering method and considered for classification model development.

3.3.2. Epoch-Based Document-of-Words Features

The document-of-words feature representation approach mainly includes three steps: data segmentation, vocabulary generation, and feature calculation, as described in [23]. Every subject’s data (i.e., 100 HZ accelerometer data for duration of one week per individual) are a time series sequence. Throughout the three main steps of this approach, every subject’s data were split into overlapping segments using a sliding window. Then, a clustering method grouped the similar segments into an appropriate number of clusters (explained in the following paragraphs), and similar segments within the same cluster were assigned the same word. Using this approach, each subject’s wrist movement time series was converted to a document of words. More details on this method can be found in the following paragraphs.
Data segmentation was performed using the method explained in Section 3.2. Once the segments or subsequences of data were generated, they needed to be represented using a feature descriptor set. This step was mainly needed for the next step (vocabulary generation), in which a clustering method identified the optimum number of clusters (words) based on the segments’ feature descriptor set. Signal vector magnitude has been used as a subsequence descriptor for assessing activities of daily living [37], while various combinations of different features are used to study the performance of the document-of-words model on identifying PD using ankle accelerometer data [23]. In this study, we used the best-identified feature set for the statistical features as our subsequences’ descriptor set (SDS). The number of total segments calculated for all subjects (considering that we have n number of subjects) per windowing strategy was calculated as follows (Equation (4)):
N = i = 1 n 2 × ( t i s ) 1
where ti denotes the duration of data collection and s is the window size.
In the second step, we generated a vocabulary set in which each subsequence was considered as a word. Since it is unlikely to find two identical subsequences with the same feature descriptors, we put all the segments obtained from all subjects together and grouped them using K-medoid clustering [23]. Therefore, segments that share similar characteristics are grouped together. To identify the number of clusters (k), we used the elbow method based on within-cluster sum of square (WCSS). The optimum number of clusters determines the number of words needed to represent all individuals’ wrist movement patterns. More details on how the elbow method using WCSS works can be found in [23]. Figure 5 depicts how the second step was performed.
In the last step, every segment was assigned a word, which was assigned to its cluster centroid. Using these words, we could convert every subject’s wrist data into a document of words and construct a vector D = {d0, d1, …, dk} per subject’s data, in which di is a normalized value of word frequency based on the total number of words per individual (Equation (5)):
t f ( w i , D ) = d i = f r e q u e n c y   o f   w o r d   w i   i n   D t o t a l   n u m b e r   o f   w o r d s   i n   D

3.4. Feature Selection, Classification Algorithms, and Validation

Feature selection is mainly performed in the analytics domain so that relevant features are selected for model development while irrelevant and redundant features are removed to reduce the model complexity and avoid overfitting. Since genetic algorithm (GA) has been shown to be one of the best algorithms for optimizing the classification of movement data in the healthcare domain, it was used for feature selection in our work. GA is a method inspired by natural selection to solve optimization problems by taking a population of individuals and trying to produce better approximations. The function to optimize here is the performance of our classification models. Each individual was considered to be a binary representation of features. A total of 40 random individuals were generated as the initial population. Our fitness function was defined as a combination of the average accuracy of all classifiers in this study combined with the number of features. The best fitness score is the one with the highest accuracy and the lowest number of features. We ran GA with 50 generations. In each generation, 20 individuals were selected for crossover, half of whom were the fittest and the other half random. This way, we avoided the issue of falling into a local optimum. Crossover was performed four times using random pairs, generating four feature sets. Using crossover, each feature set was generated by the features randomly selected from the parent feature sets. Then, each individual was given a 14% chance of being mutated, creating the new feature sets. More details on how genetic algorithm was applied to our feature set can be found in [47].
To create a model that distinguishes between healthy subjects and patients with PD, we applied different classifiers, including support vector machine (SVM with an RBF kernel), decision tree (DT), random forest (RF), K nearest neighbors (KNN), AdaBoost (AB), and naive Bayes (NB). These classifiers have been used in the gait analysis and document classification literature [5,23,48].
To evaluate each classification, we performed 100 runs of that classifier with 10-fold cross-validation, such that in each run of cross validation, nine folds (56 individuals’ data) were used for training and the remaining one fold (six subjects’ data) was considered as the validation set. Reported results are based on the average accuracy values and standard deviation around the average. This validation strategy is performed to make sure that the results are robust enough.

3.5. Hypotheses Testing

3.5.1. Comparing Statistical and Document-of-Words Feature Engineering Methods

To reach our second goal, in addition to the observations and results obtained from the classification models, we developed a hypothesis to compare the average performance of the two feature engineering methods. The following hypothesis was developed:
H0: 
meanStatisticalFeatures ≥ meanWordFeatures.
H0a: 
meanStatisticalFeatures < meanWordFeatures.
in which meanStatisticalFeatures and meanWordFeatures are the average accuracy values across all classifiers in the case of using statistical features and word features, respectively. The feature set used for this hypothesis test was generated using the best windowing strategy/the best combination of windowing strategies identified prior to this step.
Considering one set of features, 100 runs per classifier, and six classifiers, we must have 600 accuracy values per method.

3.5.2. Examining the Impact of Different Windowing Strategies on Classification Performance

To investigate the impact of window size on classification performance, we applied different classification methods using features calculated from subsequences of each windowing strategy (3 s, 10 s, 60 s, 300 s, and 900 s), different combinations of all features together, and a reduced feature set selected from the combination of all. The performance of the classifiers based on each single window strategy was considered as baseline. Considering 100 runs of each classifier and 10-fold cross-validation, firstly, all authors investigated the results by observation and identified windowing strategy/strategies that seemed to provide the best performance values. Then, we took the best ones and statistically checked for any significant difference between them and baseline windowing strategies. Here is an example of one hypothesis test developed:
H0: 
PerformanceCombination_of_Aall ≤ Performance3 Seconds.
H0a: 
PerformanceCombination_of_All > Performance3 Seconds.
in which PerformanceCombination_of_All represents the average accuracy obtained using the features calculated from all data granularity levels, and Performance3 Seconds represents the average accuracy of obtained using features calculated out of subsequences of 3 s length. Similar hypotheses were developed for the comparison of each pair of best performing window size/sizes and for the comparison of them with the baseline windowing strategies.

3.5.3. Examining the Impact of the Data Size on Classification Performance

To identify the optimal number of days for data collection for early detection of PD, we selected the best windowing strategy/strategies and the best set of features (Section 3.5.2). Classification of PD against healthy elderlies was performed considering all possible number of days of data collection, consisting of one, two, three, four, five, six, and seven consecutive days. The best performing classification model and the best subset of features were used across all possible classification tasks based on different numbers of consecutive days. The classification task using each data size was performed 100 times with 10-fold cross-validation. Then, the accuracy values were taken for statistical analysis to compare the classification performance across different amounts of data. The developed hypothesis is as follows:
H0: 
performance_1 Day = performance_2 Days = performance_3 Days = performance_4 Days = performance_5 Days = performance_6 Day = performance_7 Day.
H0a: 
not all the performance values from different amounts of data are the same.
The null hypothesis states that the classification performance values are not significantly different across different numbers of days from which the data were collected.

4. Results

4.1. Classification Results Using Statistical Features

Using the epoch-based statistical feature representation, the average value of each feature across all subsequences was calculated. First, all possible combinations of the statistical features were investigated by running the classifiers. The best results were obtained using SDS = {RMSX, RMSY, RMSZ, Vector Magnitude, RangeAP, RangeML, RangeVER, DX, DY, DZ}. Then, using this SDS, each classifier was run 100 times using features from various subsequence lengths. The best results were obtained using a reduced set of features containing at least one feature from each windowing strategy.
Table 2 shows the average and standard error of accuracy values, and Figure 6 indicates the distribution of accuracy values obtained using features from different windowing strategies, a combination of features from all windowing strategies, and a reduced set of features. While we investigated all the combinations of different windowing strategies and their reduced subsets, we only listed the results for the baseline windowing strategies and two combinations with the highest performance.
As it appears, the best classification performance was obtained using a reduced set of features containing features from all windowing strategies. This set includes the following features: {DispersionAP_10 s, DispresionAP _60 s, RMSRAP-300, RMSRML-300, signal_vector_magnitude_300 s, signal_vector_magnitude_900 s, RMSRAP-3 s, ARML_3 s, ARAp_3 s, ARVer _3 s}.

4.2. Classification Results Using Document-of-Words Features

To create the features using the document-of-words method, we used the same SDS introduced in the statistical feature method. Various values of subsequence lengths and various numbers of clusters were investigated. The optimum number of clusters/values of k were selected based on the elbow method for each and all windowing strategies. This number ranged between 4 and 32 for each window size. The value of k determines the number of words needed to represent all individuals. After investigating different values of k (k ∈ {4,8,12,16,20,24,28,32}), the optimum number of clusters for subsequences of 3 s, 10 s, 60 s, 300 s, and 900 s were found to be 22, 22, 20, 20, and 16, respectively. Therefore, each subject’s wrist movement signal could be represented as follows:
Pi-3 s ={D1, D2, D3, …, D22}
Pi-10 s = {D1, D2, D3, …, D22}
Pi-1 min = {D1, D2, D3, …, D20}
Pi-5 min = {D1, D2, D3, …, D20}
Pi-15 min = {D1, D2, D3, …, D16}
where Pi-3 s, Pi-10 s, Pi-1 min, Pi-5 min, and Pi-15 min represent the ith person’s wrist data when the signal is split into subsequences of 3 s, 10 s, 1 min, 5 min, and 15 min, respectively. Combining all the words, each individual’s wrist data could be represented by a document of 100 words as follows:
Pi = {D1, D2, D3, …, D100}
Table 3 shows the average and standard error of accuracy values, and Figure 7 indicates the distribution of accuracy values obtained using the document-of-words features from different windowing strategies, a combination of features from all windowing strategies, and a reduced set of features.

4.3. Impact of Windowing Strategies

For each hypothesis test, we first checked the normality, and since our data were normally distributed, we used ANOVA to compare the average accuracy values based on the different windowing strategies listed in Table 3. We only used the data from the document-of-words approach. Our result from a one-way ANOVA indicated that the means of the seven conditions were unequal—F(6, 593) = 4.07, p = 0.022. Tukey’s HSD test for multiple comparisons found that the mean accuracy value from the reduced feature set containing features from all granularity levels was significantly higher than that of all other features sets obtained from the other windowing strategies (p = 0.024, 95% confidence level), except for the one containing the combination of all features.

4.4. Comparing the Models Performance between Document-of-Words and Statistical Feature Engineering

In general, if we look at the results in Table 2 and Table 3, we can see that the classification performance based on the document-of-words features is higher than that of the statistical features. To validate this observation statistically using a hypothesis test, we only selected the best classification results based on the windowing strategy/strategies out of each feature engineering method. Since the highest performance of each method was obtained using the reduced subset of features out of the combination of all window sizes, that feature set was used for our hypothesis test.
The normality test was demonstrated that the distribution of accuracy values per group was normal. Therefore, a paired-sample t-test was performed. The test outcome indicates that the average accuracy value of the classification using word features is significantly greater that the average accuracy value of the classification using statistical features (t = 6.85, df = 4797.8, p-value = 1.27 × 10−6).

4.5. Optimum Number of Days of Data Collection

In order to find the optimum number of days of data collection for the identification of PD patients from the healthy elderlies, we identified the best classification performance, which was obtained from the reduced subset of twenty-words features. This subset of word features calculated over 7 days of data were used as our baseline so that we could compare the classification performance values based on different amounts of data collected over various numbers of days. This way, we could identify the optimum number of days for data collection that provides us with the same differentiating power as seven days of the data would. We only used the results of our SVM classifier with the RBF kernel for this section of the study.
Upon the result of the normality (p = 0.001) and homogeneity (constant variance) check, we could not use parametric test. Thus, the Kruskal–Wallis non-parametric test was selected and performed, which resulted in a chi-squared = 40.951, with df = 6, and p-value = 2.96 × 10−7. There is enough evidence that the classification performance considering different amounts of data (day-wise) is significantly different.
To perform a pairwise comparison, Dunn’s test with Bonferroni correction was performed. The results of Dunn’s test can be seen in Table 4. As can be seen in the table, considering one day of data provides significantly different results (lower accuracy) compared to other amounts of data, except for two days of data. Therefore, at least three days of data collection are needed.

5. Discussion

The contributions of this study are as follows: (1) Early detection of patients with PD using accelerometers worn on the most common and accepted body location, the wrist. Our results indicate that PD can be detected using wrist-worn accelerometry and passive monitoring of individuals with good accuracy values. (2) Application of the document-of-words feature engineering method, which is novel for PD detection using data collected from wrist sensors. Although this method was used for PD detection using data collected from participants’ ankles [23], the dataset was relatively small, and the data were not collected from wrist. Moreover, this is the first time that the impact of the document-of-words feature engineering method on early detection of patients with PD is statistically tested, showing that this feature engineering method significantly improves the classification performance compared to the traditional statistical feature engineering method. (3) Identifying the best windowing strategy/strategies for the early detection of PD using wrist accelerometers and passive monitoring of individuals. Although different windowing strategies have been tested to find the best one/ones for identifying the symptoms of the disease [3,7,11,16,45], there is no study in the literature that has investigated the impact of window size on the early detection of PD. Furthermore, the window sizes in this study are not picked randomly but are selected based on the best window sizes reported for the identification of symptoms of the disease [3,7,11,16,45]. (4) Identifying the optimum amount of data collection, in terms of the number of days, for the early detection of PD using wrist accelerometer data and passive monitoring. Considering continuous monitoring of individuals over time, which is a necessary step toward real-life applications of movement analysis in the context of smart homes and the healthcare domain [5,6], including a huge amount of data in analysis leads to a more computationally intensive process. If we can show that the same task can be completed using less amount of data and without compromising accuracy, it would help with saving memory and computational resources and speed up the early detection of PD. Therefore, an additional important finding of our study is that although collecting data over seven days improves people’s movement profiles and provide us with the best classification results, three days of data collection is not significantly different.
In this study, we used two feature engineering methods: epoch-based traditional statistical feature engineering method and the document-of-words approach. Our statistical feature set, which was also used as the subsequence descriptor set in the document-of-words approach, included the following set of ten features: {signal vector magnitude, RMSRAP, RMSRML, RMSRVer, DispersionAP, DispresionML, DispersionVer, ARAp, ARML, ARVer}. This set was selected as the best-performing set of statistical features. Various classification models were applied to our dataset using 10-fold cross-validation, and the results were reported per 100 runs of each classifier.
Based on our results, the document-of-words method significantly improves the classification performance compared to the statistical method. Interestingly, our results indicate that the highest performance of both the epoch-based statistical approach and the document-of-words method was obtained when we used a combination of features from all windowing strategies (3 s, 10 s, 60 s, 300 s, and 900 s). This set of features was a reduced feature set out of all possible features generated from each windowing strategy. The statistical set of features that contributed to the best classification results is as follows: {DispersionAP_10 s, DispresionAP _60 s, RMSRAP-300, RMSRML-300, signal_vector_magnitude_300 s, signal_vector_magnitude_900 s, RMSRAP-3 s, ARML_3 s, ARAp_3 s, ARVer _3 s}. The reason why AR in all three directions using a window size of three is included in the selected feature set could be because of its importance in identifying tremor and bradykinesia [16]. The existence of RMSRAP from 3 s window sizes could be explained using the same rationale. Dispersion features using a window size of 10 has been shown to contribute positively to the diagnosis of PD [11], and that could be why DispersionAP_10 s appeared as one of the important features in our reduced statistical feature set. Vector magnitude features from windows of sizes 300 and 900, and RMSRAP-300 and RMSRML-300, can be also seen in the selected subset of features, which could be indicating the amount of hand movement over longer periods of time. Using the document-of-words features, the best performance was obtained using a reduced set of 20 words while the set had words generated from all windowing strategies. This indicates that considering data from each granularity level reveals some useful information about the populations in the study, and that contributes to a better performance of classification tasks.
This study is a step towards the early diagnosis of Parkinson’s disease using wrist-worn accelerometry in a real-life setting and would benefit patients and healthcare providers. The document-of-words feature engineering method eliminates the need for precise step/stride segmentation, which is a challenging task specifically for populations with a health condition.
It is not easy to compare our work and the results of this study with those from other studies because there has not been a great number of works that have attempted the early detection of PD using wrist-worn sensors and passive monitoring. Moreover, no study has reported on the impact of window size on PD diagnosis. Table 5 includes information on related works and specifies the similarities and differences between our work and theirs. As it may be noticed in the table, two studies that applied the document-of-words feature engineering method were both performed in controlled environments, one for the identifying of daily activities and the other one for PD diagnosis using accelerometry data from participants’ ankles [23,37]. By using participants’ ankle accelerometry, authors could achieve 80% accuracy [23] for the early diagnosis of PD, which is less than the average accuracy achieved in our study based on the document-of-words method and combining features from all windowing strategies. Lipsmeier analyzed data collected from PD patients and healthy age-matched individuals collected through passive monitoring and active scripted activates [27]. They reported 75% sensitivity and 81% specificity only using the turning speed feature but did not obtain these results using an independent test set. Moreover, they neither investigated the impact of different windowing strategies nor the impact of different amounts of data on models’ performance values. Evers et al.’s study focused on the diagnosis of PD using wrist-worn accelerometry and passive activity data collection [42]. By using data from the most affected wrist, they could achieve 0.75 AUC, which is less than the AUC value obtained in our analysis. The closest research study to our work is the one published by Williamson et al. [11], in which they used wrist-worn accelerometers and passive monitoring of the subjects. Since they reported AUC as the measure of performance in their study, we averaged AUC values obtained from 100 runs of the SVM classifier with 10-fold cross-validation and using the best subset of word features. We detected PD with average AUC = 86, which is very close to the results obtained by Williamson et al., who gained AUC = 0.85. This suggests that PD detection performance might be improved by considering low-movement and high-movement segments and their associated features, as reported in Williamson’s study [11]. Innovations of our work compared to Williamson et al.’s study [11] are as follows: we used the document-of-words feature engineering method, the impact of a different windowing strategy on classification models’ performance was investigated in our work, and we statistically identified the optimal amount of data that need to be collected for the early diagnosis of PD.
The methods and dataset used in this study have strengths and limitations. One of the strengths of the dataset is that it was passively collected from the most common and acceptable body location of individuals in a free-living condition. Moreover, our methodology does not require participants to wear the sensor on the most-affected wrist. Another strength of this work is that PD can be detected without explicit labeling of the data, e.g., symptoms detection and labeling. We did not have to label the data for activities of daily living and the symptoms of the disease. Furthermore, we had a balanced dataset in terms of the number of patients with PD and healthy elderlies. The other strength of our approach is that we had a set of features that was independent from sensor orientation.
Our dataset’s limitations are as follows: a limited number of participants and lack of additional information from each subject (e.g., if they presented any symptoms, and labeling the data or logs of their daily activities).
There are several opportunities for future work with the same dataset. Since we have three phases of data collection with 6–12 months in between, we should be able to test our results by combining participants from all three phases and see how that improves the results. Another possible task would be to select our test sets completely independently from the training set by selecting individuals for our test set from the other phase of the study. Another promising next step would be to conduct a longitudinal study and predict each individual’s health status in the third phase based on the first two phases of data. Another next step would be to develop models based on more sophisticated machine learning models to distinguish PD patients with different disease severities.
As mentioned in [11], wrist-worn accelerometry could be one modality of data collection that can be paired with other modalities for detecting PD and tracking disease severity. Combining wrist-worn accelerometry with other data sources, including sleep quality and speech analysis, could improve the power of PD detection algorithms and consequently improve this population’s quality of life.

6. Conclusions

This research study had four aims, including the early diagnosis of Parkinson’s disease in a free-living environment using wrist accelerometry, investigating the impact of two feature engineering methods on classification models’ performance, investigating the impact of different windowing strategies and their associated features on models’ performance, and investigating the impact of the amount of data on classification models’ performance. Our results indicate the feasibility of using data from wrist-worn accelerometry collected in free-living conditions for the early diagnosis of PD. The performance of our PD diagnosis models were improved using the document-of-words feature engineering method, and a subset of features including at least one feature from each windowing strategy was included. This indicates that including data from different granularity levels adds to our model performance. By comparing our models’ performance values using different amounts of data in terms of days, we noticed an increase in the model performance as the amount of data increased. However, our models’ performance values using three, four, five, and six days of data collection were not significantly different from that of seven days, suggesting the possibility of using three days of data collection for the diagnosis of PD without compromising the diagnosis model performance.

Author Contributions

Conceptualization, E.R.; methodology, E.R.; validation, E.R.; formal analysis, E.R.; resources, V.M., E.R. and H.A.; data curation, E.R. and V.M.; writing—original draft preparation, E.R.; writing—review and editing, E.R., V.M. and H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was approved by the Institutional Review Board of University of Nebraska Medical Center (UNMC) (protocol code: 0057-16-EP, approved on 5 November 2016).

Informed Consent Statement

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

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stamford, J.A.; Schmidt, P.N.; Friedl, K.E. What engineering technology could do for quality of life in Parkinson’s disease: A review of current needs and opportunities. IEEE J. Biomed. Health Inform. 2015, 19, 1862–1872. [Google Scholar] [CrossRef] [PubMed]
  2. Postuma, R.B.; Berg, D. Advances in markers of prodromal Parkinson disease. Nat. Rev. Neurol. 2016, 12, 622. [Google Scholar] [CrossRef]
  3. Habets, J.G.; Herff, C.; Kubben, P.L.; Kuijf, M.L.; Temel, Y.; Evers, L.J.; Bloem, B.R.; Starr, P.A.; Gilron, R.E.; Little, S. Rapid dynamic naturalistic monitoring of bradykinesia in Parkinson’s disease using a wrist-worn accelerometer. Sensors 2021, 21, 7876. [Google Scholar] [CrossRef] [PubMed]
  4. Patel, S.; Lorincz, K.; Hughes, R.; Huggins, N.; Growdon, J.; Standaert, D.; Akay, M.; Dy, J.; Welsh, M.; Bonato, P. Monitoring motor fluctuations in patients with Parkinson’s disease using wearable sensors. IEEE Trans. Inf. Technol. Biomed. 2009, 13, 864–873. [Google Scholar] [CrossRef] [Green Version]
  5. Gubbi, J.; Buyya, R.; Marusic, S.; Palaniswami, M. Internet of Things (IoT): A vision, architectural elements, and future di-rections. Future Gener. Comput. Syst. 2013, 29, 1645–1660. [Google Scholar] [CrossRef] [Green Version]
  6. Alemdar, H.; Ersoy, C. Wireless sensor networks for healthcare: A survey. Comput. Netw. 2010, 54, 2688–2710. [Google Scholar] [CrossRef]
  7. Mazilu, S.; Blanke, U.; Calatroni, A.; Gazit, E.; Hausdorff, J.M.; Tröster, G. The role of wrist-mounted inertial sensors in detecting gait freeze episodes in Parkinson’s disease. Pervasive Mob. Comput. 2016, 33, 1–16. [Google Scholar] [CrossRef]
  8. Mazilu, S.; Blanke, U.; Tröster, G. Gait, wrist, and sensors: Detecting freezing of gait in Parkinson’s disease from wrist movement. In Proceedings of the 2015 IEEE International Conference on Pervasive Computing and Communication Workshops (PerCom Workshops), St. Louis, MI, USA, 23–27 March 2015; pp. 579–584. [Google Scholar]
  9. Roy, S.H.; Cole, B.T.; Gilmore, L.D.; De Luca, C.J.; Thomas, C.A.; Saint-Hilaire, M.M.; Nawab, S.H. High-resolution tracking of motor disorders in Parkinson’s disease during unconstrained activity. Mov. Disord. 2013, 28, 1080–1087. [Google Scholar] [CrossRef]
  10. Iosa, M.; Picerno, P.; Paolucci, S.; Morone, G. Wearable inertial sensors for human movement analysis. Expert Rev. Med. Devices 2016, 13, 641–659. [Google Scholar] [CrossRef]
  11. Williamson, J.R.; Telfer, B.; Mullany, R.; Friedl, K.E. Detecting Parkinson’s disease from wrist-worn accelerometry in the UK Biobank. Sensors 2021, 21, 2047. [Google Scholar] [CrossRef]
  12. Pradhan, S.; Kelly, V.E. Quantifying physical activity in early Parkinson disease using a commercial activity monitor. Park. Relat. Disord. 2019, 66, 171–175. [Google Scholar] [CrossRef] [PubMed]
  13. Mirelman, A.; Ben Or Frank, M.; Melamed, M.; Granovsky, L.; Nieuwboer, A.; Rochester, L.; Del Din, S.; Avanzino, L.; Pelosin, E.; Bloem, B.R.; et al. Detecting Sensitive Mobility Features for Parkinson’s Disease Stages Via Machine Learning. Mov. Disord. 2021, 36, 2144–2155. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, S.D.; Allen, N.E.; Canning, C.G.; Fung, V.S. Postural instability in patients with Parkinson’s disease. CNS Drugs 2012, 27, 97–112. [Google Scholar] [CrossRef] [PubMed]
  15. Jankovic, J. Parkinson’s disease: Clinical features and diagnosis. J. Neurol. Neurosurg. Psychiatry 2008, 79, 368–376. [Google Scholar] [CrossRef] [Green Version]
  16. 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, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhan, A.; Mohan, S.; Tarolli, C.; Schneider, R.B.; Adams, J.L.; Sharma, S.; Elson, M.J.; Spear, K.L.; Glidden, A.M.; Little, M.A.; et al. Using smartphones and machine learning to quantify Parkinson disease severity: The mobile Parkinson disease score. JAMA Neurol. 2018, 75, 876–880. [Google Scholar] [CrossRef] [PubMed]
  18. Balaji, E.; Brindha, D.; Balakrishnan, R. Supervised machine learning based gait classification system for early detection and stage classification of Parkinson’s disease. Appl. Soft Comput. 2020, 94, 106494. [Google Scholar]
  19. Doherty, A.; Jackson, D.; Hammerla, N.; Plötz, T.; Olivier, P.; Granat, M.H.; White, T.; Van Hees, V.T.; Trenell, M.I.; Owen, C.G.; et al. Large scale population assessment of physical activity using wrist worn accelerometers: The UK biobank study. PLoS ONE 2017, 12, e0169649. [Google Scholar] [CrossRef] [Green Version]
  20. Rastegari, E.; Ali, H. A Correlation Network Model Utilizing Gait Parameters for Evaluating Health Levels. In Proceedings of the 8th ACM International Conference on Bioinformatics, Computational Biology, and Health Informatics, Boston, MA, USA, 20–23 August 2017; ACM: New York, NY, USA, 2017; pp. 568–574. Available online: http://dl.acm.org/citation.cfm?id=3107487 (accessed on 20 August 2017).
  21. Aich, S.; Youn, J.; Chakraborty, S.; Pradhan, P.M.; Park, J.H.; Park, S.; Park, J. A supervised machine learning approach to detect the on/off state in Parkinson’s disease using wearable based gait signals. Diagnostics 2020, 10, 421. [Google Scholar] [CrossRef] [PubMed]
  22. Kheirkhahan, M.; Mehta, S.; Nath, M.; Wanigatunga, A.A.; Corbett, D.B.; Manini, T.M.; Ranka, S. A bag-of-words approach for assessing activities of daily living using wrist accelerometer data. In Proceedings of the 2017 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Kansas City, MO, USA, 13–16 November 2017; pp. 678–685. [Google Scholar]
  23. Rastegari, E.; Ali, H. A bag-of-words feature engineering approach for assessing health conditions using accelerometer data. Smart Health 2020, 16, 100116. [Google Scholar] [CrossRef]
  24. Habib, M.A.; Mohktar, M.S.; Kamaruzzaman, S.B.; Lim, K.S.; Pin, T.M.; Ibrahim, F. Smartphone-based solutions for fall detection and prevention: Challenges and open issues. Sensors 2014, 14, 7181–7208. [Google Scholar] [CrossRef] [PubMed]
  25. Powers, R.; Etezadi-Amoli, M.; Arnold, E.M.; Kianian, S.; Mance, I.; Gibiansky, M.; Trietsch, D.; Alvarado, A.S.; Kretlow, J.D.; Her-rington, 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] [PubMed]
  26. Ossig, C.; Gandor, F.; Fauser, M.; Bosredon, C.; Churilov, L.; Reichmann, H.; Horne, M.K.; Ebersbach, G.; Storch, A. Correlation of quantitative motor state assessment using a kinetograph and patient diaries in advanced PD: Data from an observational study. PLoS ONE 2016, 11, e0161559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lipsmeier, F.; Taylor, K.I.; Kilchenmann, T.; Wolf, D.; Scotland, A.; Schjodt-Eriksen, J.; Cheng, W.Y.; Fernandez-Garcia, I.; Siebourg-Polster, J.; Jin, L.; et al. Evaluation of smartphone-based testing to generate exploratory outcome measures in a phase 1 Parkinson’s disease clinical trial. Mov. Disord. 2018, 33, 1287–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Rodríguez-Molinero, A.; Samà, A.; Pérez-López, C.; Rodríguez-Martín, D.; Alcaine, S.; Mestre, B.; Quispe, P.; Giuliani, B.; Vainstein, G.; Browne, P.; et al. Analysis of correlation between an accelerometer-based algorithm for detecting parkinsonian gait and UPDRS subscales. Front. Neurol. 2017, 8, 431. [Google Scholar] [CrossRef] [Green Version]
  29. Czech, M.; Demanuele, C.; Erb, M.K.; Ramos, V.; Zhang, H.; Ho, B.; Patel, S. The impact of reducing the number of wearable devices on measuring gait in parkinson disease: Noninterventional exploratory study. JMIR Rehabilitation Assist. Technol. 2020, 7, e17986. [Google Scholar] [CrossRef]
  30. Silva de Lima, A.L.; Evers, L.J.; Hahn, T.; Bataille, L.; Hamilton, J.L.; Little, M.A.; Okuma, Y.; Bloem, B.R.; Faber, M.J. Freezing of gait and fall detection in Parkinson’s disease using wearable sensors: A systematic review. J. Neurol. 2017, 264, 1642–1654. [Google Scholar] [CrossRef] [Green Version]
  31. Raiano, L.; di Pino, G.; di Biase, L.; Tombini, M.; Tagliamonte, N.L.; Formica, D. PDMeter: A wrist wearable device for an at-home assessment of the Parkinson’s disease ri-gidity. IEEE Trans. Neural Syst. Rehabil. Eng. 2020, 28, 1325–1333. [Google Scholar] [CrossRef] [PubMed]
  32. Gokul, H.; Suresh, P.; Vignesh, B.H.; Kumaar, R.P.; Vijayaraghavan, V. Gait recovery system for parkinson’s disease using machine learning on embedded platforms. In Proceedings of the 2020 IEEE International Systems Conference (SysCon), Montreal, QC, Canada, 24 August–20 September 2020; pp. 1–8. [Google Scholar]
  33. Goud, K.Y.; Moonla, C.; Mishra, R.K.; Yu, C.; Narayan, R.; Litvan, I.; Wang, J. Wearable electrochemical microneedle sensor for continuous monitoring of levodopa: Toward Parkinson management. ACS Sensors 2019, 4, 2196–2204. [Google Scholar] [CrossRef]
  34. Nóbrega, L.R.; Cabral, A.M.; Oliveira, F.H.; de Oliveira Andrade, A.; Krishnan, S.; Pereira, A.A. Wrist Movement Variability Assessment in Individuals with Parkinson’s Disease. Healthcare 2022, 10, 1656. [Google Scholar] [CrossRef] [PubMed]
  35. Juutinen, M.; Wang, C.; Zhu, J.; Haladjian, J.; Ruokolainen, J.; Puustinen, J.; Vehkaoja, A. Parkinson’s disease detection from 20-step walking tests using inertial sensors of a smartphone: Machine learning approach based on an observational case-control study. PLoS ONE 2020, 15, e0236258. [Google Scholar] [CrossRef] [PubMed]
  36. Rastegari, E.; Azizian, S.; Ali, H. Machine Learning and Similarity Network Approaches to Support Automatic Classification of Parkinson’s Diseases Using Accelerometer-Based Gait Analysis. 2019. Available online: http://scholarspace.manoa.hawaii.edu/handle/10125/59861 (accessed on 17 March 2019).
  37. Kheirkhahan, M.; Chakraborty, A.; Wanigatunga, A.A.; Corbett, D.B.; Manini, T.M.; Ranka, S. Wrist accelerometer shape feature derivation methods for assessing activities of daily living. BMC Med. Inform. Decis. Mak. 2018, 18, 124. [Google Scholar] [CrossRef]
  38. Del Din, S.; Godfrey, A.; Mazzà, C.; Lord, S.; Rochester, L. Free-living monitoring of Parkinson’s disease: Lessons from the field. Mov. Disord. 2016, 31, 1293–1313. [Google Scholar] [CrossRef]
  39. Warmerdam, E.; Hausdorff, J.M.; Atrsaei, A.; Zhou, Y.; Mirelman, A.; Aminian, K.; Espay, A.J.; Hansen, C.; Evers, L.J.; Keller, A.; et al. Long-term unsupervised mobility assessment in movement disorders. Lancet Neurol. 2020, 19, 462–470. [Google Scholar] [CrossRef] [PubMed]
  40. Pasluosta, C.F.; Gassner, H.; Winkler, J.; Klucken, J.; Eskofier, B.M. An Emerging Era in the Management of Parkinson’s Disease: Wearable Technologies and the Internet of Things. IEEE J. Biomed. Health Inform. 2015, 19, 1873–1881. [Google Scholar] [CrossRef]
  41. San-Segundo, R.; Zhang, A.; Cebulla, A.; Panev, S.; Tabor, G.; Stebbins, K.; Massa, R.E.; Whitford, A.; De la Torre, F.; Hodgins, J. Parkinson’s disease tremor detection in the wild using wearable accelerometers. Sensors 2020, 20, 5817. [Google Scholar] [CrossRef] [PubMed]
  42. Evers, L.J.; Raykov, Y.P.; Krijthe, J.H.; De Lima, A.L.; Badawy, R.; Claes, K.; Heskes, T.M.; Little, M.A.; Meinders, M.J.; Bloem, B.R. Real-life gait performance as a digital bi-omarker for motor fluctuations: The Parkinson@ Home validation study. J. Med. Internet Res. 2020, 22, e19068. [Google Scholar] [CrossRef]
  43. Kleanthous, N.; Hussain, A.J.; Khan, W.; Liatsis, P. A new machine learning based approach to predict Freezing of Gait. Pattern Recognit. Lett. 2020, 140, 119–126. [Google Scholar] [CrossRef]
  44. Shawen, N.; O’Brien, M.K.; Venkatesan, S.; Lonini, L.; Simuni, T.; Hamilton, J.L.; Ghaffari, R.; Rogers, J.A.; Jayaraman, A. Role of data measurement characteristics in the accurate detection of Parkinson’s disease symptoms using wearable sensors. J. Neuroeng. Rehabil. 2020, 17, 52. [Google Scholar] [CrossRef] [Green Version]
  45. Hssayeni, M.D.; Jimenez-Shahed, J.; Burack, M.A.; Ghoraani, B. Ensemble deep model for continuous estimation of Unified Parkinson’s Disease Rating Scale III. Biomed. Eng. Online 2021, 20, 32. [Google Scholar] [CrossRef]
  46. Kuijf, M.L.; Kubben, P.L.; Herff, C. Evaluation of Parkinson’s Disease at Home: Predicting Tremor from Wearable Sensors. In Prediction and Real-Life Monitoring of DBS Motor Response in Parkinson’s Disease; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar]
  47. Rastegari, E.; Orn, D.; Ali, H. Smart Computational Approaches with Advanced Feature Selection Algorithms for Optimiz-ing the Classification of Mobility Data in Health Informatics. In Proceedings of the 11th ACM International Conference on Bioinformatics, Computational Biology and Health Informatics, Virtual, 21–24 September 2020; pp. 1–9. [Google Scholar]
  48. Pla, A.; Mordvanyuk, N.; López, B.; Raaben, M.; Blokhuis, T.J.; Holstlag, H.R. Bag-of-steps: Predicting lower-limb fracture rehabilitation length by weight loading analysis. Neurocomputing 2017, 268, 109–115. [Google Scholar] [CrossRef]
Figure 1. Steps to identify the best feature engineering method and the best windowing strategy/strategies.
Figure 1. Steps to identify the best feature engineering method and the best windowing strategy/strategies.
Sensors 22 09122 g001
Figure 2. Steps to identify the optimal amount of data in terms of days.
Figure 2. Steps to identify the optimal amount of data in terms of days.
Sensors 22 09122 g002
Figure 3. First five minutes of raw data collected from the wrist of an individual with PD.
Figure 3. First five minutes of raw data collected from the wrist of an individual with PD.
Sensors 22 09122 g003
Figure 4. First five minutes of raw data collected from the wrist of a healthy individual.
Figure 4. First five minutes of raw data collected from the wrist of a healthy individual.
Sensors 22 09122 g004
Figure 5. Epoch-based document-of-words feature engineering method—vocabulary generation phase.
Figure 5. Epoch-based document-of-words feature engineering method—vocabulary generation phase.
Sensors 22 09122 g005
Figure 6. Distribution of accuracy values using various statistical feature sets.
Figure 6. Distribution of accuracy values using various statistical feature sets.
Sensors 22 09122 g006
Figure 7. Distribution of accuracy values using various Document-of-Words feature sets.
Figure 7. Distribution of accuracy values using various Document-of-Words feature sets.
Sensors 22 09122 g007
Table 1. Subjects’ characteristics—phase 1.
Table 1. Subjects’ characteristics—phase 1.
Healthy ElderliesPD
Subjects3228
Gender (M/F)10/2220/5
Age64.2 ± 771 ± 6.2
H&Y-1.73 ± 0.83
Table 2. Classification of PD vs. HE accuracy values—statistical feature engineering.
Table 2. Classification of PD vs. HE accuracy values—statistical feature engineering.
Window SizeFeaturesSVMNBDTRFKNNAB
3 s1067.5 ± 1868.4 ± 1561.3 ± 1867.1 ± 1670.2 ± 1862.2 ± 16
10 s1068.1 ± 1667.1 ± 1462.3 ± 1568.2 ± 1771.6 ± 1863.4 ± 15
60 s1068.4 ± 1769.4 ± 1662.2 ± 1969.2 ± 1771 ± 1862.3 ± 15
300 s1068.8 ± 1668.8 ± 1863.5 ± 1868.1 ± 1674.5 ± 1766.2 ± 19
900 s1069.2 ± 1661.9 ± 1463.6 ± 1669.4 ± 1573.7 ± 1660.8 ± 21
All 5074.9 ± 1874.2 ± 1366.5 ± 1674.8 ± 1577.4 ± 1461.8 ± 17
All1080.2 ± 1481.1 ± 1671.5 ± 1779.2 ± 1482.3 ± 1471.1 ± 16
Table 3. Classification of PD vs. HE accuracy values using document-of-words feature engineering.
Table 3. Classification of PD vs. HE accuracy values using document-of-words feature engineering.
Window Size FeaturesSVMNBDTRFKNNAB
3 s2276.2 ± 1775.2 ± 1561.1 ± 1976.6 ± 1374.4 ± 1569.8 ± 17
10 s2276.0 ± 1774.1 ± 1664.2 ± 1775.3 ± 1572.5 ± 1570.5 ± 18
60 s2078.1 ± 1677.2 ± 1663.5 ± 1677.7 ± 1476.5 ± 1869.0 ± 19
300 s2076.1 ± 1675.3 ± 1873.1 ± 1976.2 ± 1475.0 ± 1770.1 ± 17
900 s1679.2 ± 1579.1 ± 1663.2 ± 1778.0 ± 1576.0 ± 1674.3 ± 17
All 10077.8 ± 1775.6 ± 1869.3 ± 1678.2 ± 1573.0 ± 1772.0 ± 18
All-Reduced 2088.5 ± 1084.4 ± 1474.8 ± 1382.1 ± 1484.7 ± 1481.2 ± 15
Table 4. Dunn’s test pairwise comparison between the accuracy of classification using different amounts of data (day-wise). (Significant difference exist where “*” is observed).
Table 4. Dunn’s test pairwise comparison between the accuracy of classification using different amounts of data (day-wise). (Significant difference exist where “*” is observed).
Col-Mean
Row Mean
Day1Day2Day3Day4Day5Day6
Day2−2.74
0.064
-----
Day3−3.96
0.0008 *
−1.22
1.000
----
Day4−4.42
0.0001 *
−1.84
0.68
−0.69
1.000
---
Day5−5.07
0.000 *
−2.70
0.07
−1.64
1.000
−0.96
1.000
--
Day6−4.20
0.0003 *
−2.13
0.34
−1.20
1.000
−0.63
1.000
0.20
1.000
-
Day7−4.32
0.0002 *
−2.25
0.25
−1.33
1.000
−0.75
1.000
0.095
1.000
−0.10
1.000
Table 5. Related Works with a Focus on Similarities and Differences between them and our Work.
Table 5. Related Works with a Focus on Similarities and Differences between them and our Work.
StudyFree Living ConditionSensor LocationTaskFeaturesInvestigating Different Windowing StrategiesInvestigating the Impact of Amount of Data on Models’ PerformancePerfomrance MeasurePerformance Measure Value
This studyYes- passive monitoringWristPD diagnosisDocument-of-Words (20 features from all and each window size)YesYes (collected 7 days of data)Accuracy and AUC0.88 and 0.86, respectively
[11]Yes- passive monitoringWristPD diagnosis{Dispersion, Correlation Structure Features}NoYes (collecrted 7 days of data)AUC0.85
[27]Yes- passive and active monitoringSmartphonePD diagnosis{Turning speed, sit-to-stand transitions per hour, activity ratio}NoNo (collected data over 6 months for PD and 45 days for Healty individualsSpecificity and sensitivity81% and 75%, respectively, reported only for turning speed feature
[3]Yes-active monitoringWristBradykinesia detectionStatistical features, {maximum acceleration, coefficent of determination, root mean square, spectral power}Yes No (collected data for one hour pre medication and one hour post medication)AUC0.7
[37]NoWrist Identification of activities of daily livingDocument-of-words featuresNoNo (collected data from 4 visits performing activities of daily living as instructed)F1-score0.89
[23]NoAnkle PD diagnosisDocument-of-words featuresYes No (collected data from one session walking 10 meters for four timesAccuracy, precision and recall0.8, 0.7, and 0.9, respectively
[42]Yes- passive monitoringWristPD diagnosis{total power in 0.5–10 Hz, cadence, height of dominant peak, width of dominant peak}NoNo (collected data for one hour pre medication and one hour post medication)AUC0.75 for the most affected wrist and 0.49 for the least affected one
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rastegari, E.; Ali, H.; Marmelat, V. Detection of Parkinson’s Disease Using Wrist Accelerometer Data and Passive Monitoring. Sensors 2022, 22, 9122. https://doi.org/10.3390/s22239122

AMA Style

Rastegari E, Ali H, Marmelat V. Detection of Parkinson’s Disease Using Wrist Accelerometer Data and Passive Monitoring. Sensors. 2022; 22(23):9122. https://doi.org/10.3390/s22239122

Chicago/Turabian Style

Rastegari, Elham, Hesham Ali, and Vivien Marmelat. 2022. "Detection of Parkinson’s Disease Using Wrist Accelerometer Data and Passive Monitoring" Sensors 22, no. 23: 9122. https://doi.org/10.3390/s22239122

APA Style

Rastegari, E., Ali, H., & Marmelat, V. (2022). Detection of Parkinson’s Disease Using Wrist Accelerometer Data and Passive Monitoring. Sensors, 22(23), 9122. https://doi.org/10.3390/s22239122

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