Next Article in Journal
Improving the Cybersecurity Awareness of Young Adults through a Game-Based Informal Learning Strategy
Previous Article in Journal
Short-Term Water Demand Forecasting from Univariate Time Series of Water Reservoir Stations
Previous Article in Special Issue
Identification of Optimal Data Augmentation Techniques for Multimodal Time-Series Sensory Data: A Framework
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fall Detection Based on Data-Adaptive Gaussian Average Filtering Decomposition and Machine Learning

1
Department of Automatic Control Engineering, Feng Chia University, Taichung 40724, Taiwan
2
Master’s Program of Biomedical Informatics and Biomedical Engineering, Feng Chia University, Taichung 40724, Taiwan
3
Department of Electrical Engineering, National Cheng Kung University, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
Information 2024, 15(10), 606; https://doi.org/10.3390/info15100606
Submission received: 27 August 2024 / Revised: 25 September 2024 / Accepted: 30 September 2024 / Published: 3 October 2024
(This article belongs to the Special Issue Human Activity Recognition and Biomedical Signal Processing)

Abstract

:
Falls are a significant health concern leading to increased morbidity and healthcare costs, especially for the elderly. Early and accurate detection of fall events is critical for timely intervention and preventing severe complications. This study presents a novel approach to triaxial accelerometer signals by employing data-adaptive Gaussian average filtering (DAGAF) decomposition in conjunction with machine learning techniques for fall detection. The triaxial accelerometer signals from the FallAllD dataset were decomposed into intrinsic mode functions (IMFs) and a residual component, from which feature vectors were extracted to train support vector machine (SVM) and k-nearest neighbor (kNN) classifiers. Experimental results demonstrate that the combination of the first and the third IMFs with the residual component yields the highest classification accuracy of 96.34%, with SVM outperforming kNN across all performance metrics. This approach significantly improves fall detection accuracy compared to using raw accelerometer signals, highlighting its potential in enhancing wearable fall detection systems. The proposed DAGAF decomposition method not only enhances feature extraction but also provides a promising advancement in the field, suggesting its potential to increase the reliability and accuracy of fall detection in practical applications.

Graphical Abstract

1. Introduction

According to the World Health Organization (WHO), a fall is defined as an unexpected event that results in a person inadvertently coming to rest on the ground, the floor, or another lower level [1]. Falls and their related injuries pose significant health risks, particularly for the elderly. These incidents adversely affect functional independence and quality of life, and are closely associated with increased morbidity, mortality, and healthcare costs [2]. With at least 28% of individuals aged 65 or older experiencing one or more falls annually [3,4], falls have become a critical public health concern. Fatal falls rank as the second leading cause of unintentional injury deaths worldwide, surpassed only by traffic injuries [1]. Early and accurate fall detection is essential for timely intervention and preventing severe complications. Numerous approaches have been explored for fall detection, as extensively reviewed in several articles [5,6,7]. In addition to wearable devices, ambient devices, and cameras [5], other technologies such as acoustics [8], WiFi [9,10], light detection and ranging (LiDAR) technology [11,12], and millimeter-wave (mmWave) radar [13] have also been utilized for fall detection. Multimodal data fusion is another promising technique that has been investigated for this purpose [14]. Given that over 80% of fall-related fatalities occur in low- and middle-income countries [1], there is a critical demand for cost-effective fall detection solutions. Advancements in microelectromechanical system (MEMS) technology have made accelerometers widely available, making wearable devices equipped with these sensors a particularly promising approach for addressing this issue [15]. Feature extraction plays a crucial role in identifying fall events from accelerometer signals, which is essential for effective fall detection. Palmerini et al. proposed a wavelet-based approach for fall detection, in which the mother wavelet was adaptively optimized based on the input accelerometer signal [16]. The main limitations of this study include a relatively high number of false alarms and the use of a small test sample size. Time-frequency analysis using the continuous wavelet transform (CWT) has also been applied to the analysis of accelerometer signal [17]. The CWT spectrogram enables the analysis of instantaneous frequency variations over time while preserving good frequency resolution. Another method that can simultaneously present both time and frequency information is the Hilbert–Huang transform (HHT) [18], which was applied by Erfianto et al. to accelerometer signal for gait analysis [19]. However, time–frequency analysis generally requires more computing resources in terms of data storage and processing speed, which becomes especially critical in scenarios involving wearable fall alarm devices. This demand is particularly pronounced with the HHT due to its reliance on empirical mode decomposition (EMD) [18] or ensemble empirical mode decomposition (EEMD) [20]. The repeated interpolation processes required to extract the intrinsic mode functions (IMFs) significantly increase the computational burden, making the implementation of such an algorithm on a microcontroller unit (MCU) for the wearable device a challenging task. In addition to feature extraction from triaxial accelerometer signals [21], incorporating machine learning [7] and deep learning [6] to improve recognition accuracy has also become a popular research direction. Moreover, Šeketa et al. employed event-centered data segmentation to evaluate detection performance in different window configurations [22]. Santoyo-Ramón et al. utilized one-class classifiers (OCCs) by training the model solely on the majority data from activities of daily living (ADLs), achieving detection performance comparable to that obtained with supervised binary classification [23]. Liu et al. proposed a deep learning-based approach for enhancing accelerometer signals which significantly improved the fall detection accuracy at low sampling rates [24]. Despite these advancements, there remains a demand for innovative techniques that can further improve the accuracy and reliability of fall detection.
This paper explores this issue from a new perspective that has not been investigated in previous studies by utilizing a signal decomposition algorithm for triaxial accelerometer signals before feature extraction. The novel data-adaptive Gaussian average filtering (DAGAF) algorithm [25] is employed to decompose triaxial accelerometer signals into IMFs and the residual component. The performance (including accuracy, sensitivity, specificity, precision, and F 1 -score) is compared across various combinations of the decomposed components and the raw triaxial accelerometer signals using machine learning-based classifiers. Performance is also evaluated using different event-based time segments. The results demonstrate that the DAGAF decomposition improves fall detection performance, offering a promising advancement in the field and providing a new direction for future research on signal processing in wearable monitoring devices.
The rest of this paper is organized as follows. Section 2 primarily contains the dataset used, the proposed DAGAF algorithm, feature extraction, classifiers, and the performance metrics. The experimental results and the related discussion are presented in Section 3. Finally, the conclusions are given in Section 4.

2. Materials and Methods

This section is divided into two parts. The first part (Section 2.1) provides details on the dataset used, along with the computing resources (including the software and computing environment) employed in this research. The second part (Section 2.2) introduces the proposed DAGAF decomposition and the associated methods used for verification.

2.1. Materials

2.1.1. Dataset

This study utilized the FallAllD dataset, which was developed by Saleh et al. [15] and can be accessed on IEEE DataPort™ [26]. To capture movement during the experiment, inertial measurement units (IMUs) were placed on the subjects’ neck, chest, and waist. Each IMU contains a triaxial accelerometer (sampling frequency: 238 Hz, measurement range: ± 8 g), a triaxial gyroscope (sampling frequency: 238 Hz, angular rate: ± 2000 dps), a triaxial magnetometer (sampling frequency: 80 Hz, full-scale magnetic field: ± 4 Gauss), and a barometer (sampling frequency: 10 Hz). Özdemir’s study [27] found that fall detection systems using waist-worn accelerometers perform better than those with accelerometers placed at other locations. This study thus focused exclusively on triaxial accelerometer signals collected from the waist.
This dataset was originally collected from 15 healthy young adults (8 males, 7 females) with an average age of 32 years, an average height of 171 cm, and an average weight of 67 kg. To ensure data consistency and relevance, only the data from 10 subjects who performed both fall and ADL scenarios were included in the analysis. The dataset used in this research consists of 1053 ADL instances and 423 fall trials.
Each signal in the original dataset has a duration of 20 s. A previous investigation identified that an impact-defined window consisting of 2 s before the impact (backward) and 1.23 s after the impact (forward) yielded the best performance on the FallAllD dataset [28]. Accordingly, the same impact-defined data segmentation was adopted in this study. Specifically, the peak summarized energy from triaxial accelerometer signals occurred at the 2-s mark, resulting in a length of 3.23 s for each data instance in the dataset used.

2.1.2. Computing Resources

The codes for DAGAF decomposition have been developed in both Python and MATLAB (MathWorks Inc., Natick, Massachusetts, USA), and are available on GitHub [29]. The Python implementation was created using Python 3.9.12, along with the NumPy (version 1.21.5), SciPy (version 1.7.3), and Matplotlib (version 3.5.3) libraries. The MATLAB code was developed using MATLAB 2020b. Additionally, the Statistics and Machine Learning Toolbox is required for classification analyses involving support vector machine (SVM) and k -nearest neighbor ( k NN) classifiers. In this study, all of the statistical tests and computational experiments were conducted in the MATLAB environment on a MacBook Pro (Apple M1 Max, 64 GB of memory, Apple Inc., Cupertino, CA, USA).

2.2. Methods

2.2.1. Data-Adaptive Gaussian Average Filtering (DAGAF) Decomposition

The process of DAGAF decomposition shares similarities with the widely used empirical mode decomposition (EMD) [18] and its modified variant, ensemble empirical mode decomposition (EEMD) [20]. The primary difference is that DAGAF employs an iterative lowpass filtering process with the Gaussian window [25], instead of the sifting procedure that extracts the instantaneous mean from the interpolated local maxima and the minima [18,20]. Beyond addressing the commonly mentioned limitations associated with EMD and EEMD (such as the boundary effect, mode mixing, unnecessary redundant decomposition, and the lack of a rigorous mathematical formulation [30]), DAGAF has demonstrated its effectiveness in analyzing nonstationary biomedical signal across a wide range of scenarios [25]. This study is the first attempt to use DAGAF to decompose triaxial accelerometer signals into IMFs and a residual component to extract the subtle fluctuations that are difficult to observe in the original signals.
Let the signal to be analyzed be denoted as s ( n ) for 0 n N 1 . Note that s ( ) may represent the accelerometer signal from the x-, y-, and z-axis, or any other signal required for further decomposition. The discrete Gaussian window with a length of ( 2 M + 1 ) points is given by
w g ( m ) = e ( α m / M ) 2 / 2 ,   for M m M ,
where α is a parameter inversely proportional to the standard deviation of the Gaussian distribution. In this study, α is set to 4.0728, resulting in the end values being as small as 0.025% of the maximum window value. Under this condition, the characteristics of the continuous-time Gaussian function are approximately preserved in its discrete-time counterpart. A notable feature is that the corresponding spectrum retains the bell shape of the Gaussian distribution centered at 0 [31]. To maintain energy consistency after filtering, a normalized Gaussian window is used in practical applications. The normalized Gaussian window is obtained as follows:
w n ( m ) = w g ( m ) / ( l = M M w g ( l ) ) ,   for   M m M .
As previously mentioned, the signal s ( n ) is defined over the interval 0 n N 1 , while the normalized Gaussian window, as shown in Equation (2), is defined over the interval M m M . To perform the filtering operation using w n ( ) , the signal must be extended at both boundaries. The developed codes provide four types of signal extension methods [29]. Following the type used in [25], this paper also employs the “double-symmetrical reflection” extension method, which has demonstrated good performance in various scenarios. This extension is applied to both boundaries in all computational experiments. The signal extended by this method is given as follows:
s e ( n ) = 2 E { s ( ) } s ( n ) , for M n 1 s ( n ) , for   0 n N 1 2 E { s ( ) } s ( 2 N n 2 ) , for N n N + M 1 ,
where E { s ( ) } denotes the mean value of s ( n )   for   0 n N 1 . The instantaneous mean is then obtained using the moving-average operation as follows:
m i ( n ) = m = M M w n ( m ) s e ( m + n ) ,   for   0 n N 1 .
Due to the symmetric nature of the Gaussian window w n ( ) , the computation in Equation (4) is essentially a convolution sum between w n ( ) and s e ( ) . In the frequency domain, this operation corresponds to the direct multiplication between the spectra of w n ( ) and s e ( ) . Given that the spectrum of the Gaussian window forms a bell-shaped curve symmetrically centered at 0, this indicates that the moving-average operation in Equation (4) inherently functions as a lowpass filtering on s e ( ) by the normalized Gaussian window w n ( ) . The key aspect in the filtering operation is to determine the appropriate window length for w n ( ) . For this issue, the criterion proposed by Cicone et al. [32] is adopted in this study, and the value of M is determined as follows:
M = 2 ν N / κ ,
where N is the length of the signal to be analyzed, κ is the number of local extrema (including the local maxima and the local minima) present in the signal to be analyzed, ν is a parameter chosen around 1.6 (the value 1.6 is used in all computational experiments), and the symbol denotes the floor operation, which rounds a positive number down to the nearest integer. In the developed codes [29], the second derivative test is used to identify the locations of local extrema in the analyzed signal.
Once the instantaneous mean m i ( ) has been derived using Equation (4), the intrinsic mode function (IMF) can be obtained using the following expression:
r ( n ) = s ( n ) m i ( n ) ,   for   0 n N 1 .
Such a decomposition process is applied iteratively to r ( ) using the same procedure until any predefined stop criterion is met. There are three stop criteria built in the developed codes [29]. The first one is based on the energy ratio (in dB) between the original signal s ( ) and the i-th IMF, which must fall within a predefined threshold [33]. If this energy ratio exceeds the threshold, it indicates that the energy of the new IMF is significantly smaller than that of the original signal, suggesting that further decomposition is unnecessary. The second criterion assesses the energy ratio between the residual component and the separated IMF after each decomposition. If this ratio drops below a specified threshold, indicating that the residual component is nearly flat, the decomposition process will also be terminated. In the developed codes [29], the default threshold values are 20 dB for the original signal-to-IMF ratio and 0.001 for residue-to-IMF ratio, and these default values were used in all computational experiments. The third stop criterion is associated with the length of Gaussian window. Based on Equation (5), the window length is determined from the number of local extrema κ present in the analyzed signal. If the value of κ is too small, the estimated window length may exceed or be equal to the signal length. In this case, the signal to be analyzed appears to be a simple mode that is nearly flat in pattern, and it is not necessary to continue the decomposition process. Accordingly, the third stop criterion that will terminate the decomposition process is mathematically expressed as:
M ( N 1 ) / 2 .
Given the aforementioned background, the computational procedure of DAGAF is summarized as follows:
(i) Based on the number of local extrema in the signal to be analyzed, calculate the value of M using Equation (5) and verify the value of M according to Equation (7). If the stop criterion is met, terminate the process. Otherwise, construct the Gaussian window w n ( ) of length ( 2 M + 1 ) according to Equations (1) and (2).
(ii) Extend the signal at both ends according to Equation (3) and compute the instantaneous mean m i ( ) using the moving-average operation as specified in Equation (4).
(iii) Derive the IMF r ( ) according to Equation (6) and check if any energy ratio-based stop criterion has been fulfilled. If not, repeat steps (i) to (iii). Otherwise, terminate the entire procedure.
In this study, the DAGAF algorithm was employed to perform a three-level decomposition on all triaxial accelerometer signals. Specifically, each signal was processed through three iterations of decomposition, resulting in three IMFs labeled as imf 1 (the first decomposed IMF), imf 2 (the second decomposed IMF), and imf 3 (the third decomposed IMF), respectively. The residual component, which is denoted as Res in this paper, represents the remaining part of the signal after extracting the three IMFs. Figure 1 presents an example of this three-level decomposition for the z-axis accelerometer signal, corresponding to a jogging to stumbling/tripping event, recorded at the waist position of subject 1 from the FallAllD dataset.
By examining the subfigures of Figure 1 from top (Step 1) to bottom (Step 3), it can be observed that the patterns of the solid black line (representing the original raw signal in Step 1 and imf 1 , imf 2 in Steps 2 and 3) along with the dash-dotted line (representing the corresponding instantaneous mean at each step) become simpler and flatter progressively. This observation supports the decision to limit the decomposition process to three iterations in this study. The consistent effectiveness of this three-level decomposition across all triaxial accelerometer signals further reinforces the validity of this choice.

2.2.2. Feature Extraction

In this study, each sample consists of the accelerometer signals collected from three orthogonal axes. The raw data, without conversion to acceleration unit, were used directly in figure demonstration and computational analysis. Let s x ( n ) , s y ( n ) , and s z ( n ) represent either the raw signals or the combination of components obtained through the DAGAF algorithm at time index n along the x-, y-, and z-axis, respectively, where 0 n N 1 . Before extracting the features, the following energy-related vectors are derived:
Norm ( n ) = s x 2 ( n ) + s y 2 ( n ) + s z 2 ( n ) ,   Hori ( n ) = s y 2 ( n ) + s z 2 ( n ) , Vert ( n ) = s x 2 ( n ) + s y 2 ( n ) , for   0 n N 1 .
Having obtained s x ( ) , s y ( ) , s z ( ) , Norm ( ) , Hori ( ) , and Vert ( ) , the following features are used to construct the feature vector:
(1) Statistical measures: The mean, standard deviation, and variance are computed for each of the six data, resulting in 18 features.
(2) Extreme values: The maximum, minimum, and range (the difference between the maximum and minimum) are determined for each of the six data, yielding an additional 18 features.
(3) Higher-order statistics: The skewness (the third central moment) and kurtosis (the fourth central moment) are calculated for each of the six data, adding 12 more features to the vector.
(4) Pearson’s correlation coefficients: The correlation coefficients are derived for the following pairs: s x ( ) versus s y ( ) , s x ( ) versus s z ( ) , s y ( ) versus s z ( ) , Norm ( ) versus Hori ( ) , Norm ( ) versus Vert ( ) , and Hori ( ) versus Vert ( ) , adding 6 more features.
In summary, a total of 54 features are included in the feature vector for each signal to be analyzed. These features have been verified to be effective in fall detection [24,34], and therefore are also adopted in this study.

2.2.3. Classifiers

To assess the effectiveness of DAGAF decomposition in fall detection, two classical machine learning classifiers, SVM and k NN, were employed in this study. These classifiers have been demonstrated to be more reliable than other methods, such as naïve Bayes and decision trees, in fall detection tasks [24,34]. The classifiers used in this study are briefly described as follows:
(1) SVM is a powerful supervised learning algorithm designed to find the optimal hyperplane that separates data points from different classes with the maximum margin. The core idea of SVM is to maximize the margin between the closest data points (called support vectors) of different classes. This optimization problem can be formulated by minimizing
1 2 w 2
subject to the following constraint:
y k ( w X k + b ) 1 ,   k 1 ,   2 ,   ,   K ,
where w is the weight vector that defines the hyperplane, the symbol w denotes the Euclidean norm of the weight vector w , b is the bias term, X k is the feature vector of the k-th sample, y k 1 ,   1 is the corresponding class label for the k-th sample, and K is the total number of samples. The goal of SVM is to find the weight vector w and the bias b that maximize the margin between classes. When the data are not linearly separable, a kernel function can be applied to map the input features into a higher-dimensional space, allowing for the separation of more complex patterns. In this study, we use the radial basis function (RBF) kernel, which is defined as:
K RBF ( X i ,   X j ) = e γ X i X j 2 ,
where γ is a positive-valued parameter that controls the spread of the kernel, and X i X j represents the Euclidean distance between feature vectors X i and X j . The RBF kernel is particularly effective for handling nonlinear relationships in high-dimensional feature spaces. To enhance the classification performance, we used the Bayesian optimization to fine-tune the parameters of the RBF kernel, specifically the kernel parameter γ and the regularization parameter C . The regularization term controls the trade-off between maximizing the margin and minimizing the classification error. The optimization objective function can be expressed as:
C k = 1 K max ( 0 ,   1 y k ( w X k + b ) ) + 1 2 w 2 .
In this study, the SVM classifier was used to classify ADLs and fall activities based on the feature set consisting of 54 features (refer to Section 2.2.2 for details on feature extraction). These features capture diverse characteristics of the data, allowing the SVM classifier to detect the subtle differences between the ADL and fall classes. Both parameters γ and C were initially set to 1.
(2) k NN is an instance-based learning method used for classification. It classifies a new data point by identifying the k closest neighbors in the training set and assigning the most frequent class label among these neighbors. The core concept of k NN is that data points close to each other in feature space are likely to belong to the same class. Typically, proximity between data points is measured using Euclidean distance. Once the distances between the query point and all points in the training set have been obtained, the algorithm selects the k nearest data points. The class label of the query point is then determined by a majority vote among these k neighbors, with the most frequent class being assigned to the query point. k NN is a nonparametric method, making it effective for datasets where little prior knowledge is available. Previous research [28] has shown that fall detection achieves the highest accuracy with k = 3 when using Euclidean distance. Consequently, the 3NN classifier with Euclidean distance was also employed in this study. To enhance classification performance, a comprehensive set of 54 features was extracted from the proposed signal combination or the raw accelerometer signal. These features include statistical measures (mean, standard deviation, and variance), extreme values (maximum, minimum, and the range between them), higher-order statistics (skewness and kurtosis), and Pearson’s correlation coefficients. These features capture diverse characteristics of the data, enabling the k NN classifier to more effectively distinguish between ADLs and fall activities with acceptable accuracy.

2.2.4. Performance Metrics

In this study, the leave-one-subject-out (LOSO) cross-validation method was employed to evaluate the classification performance. In the LOSO approach, data from a single subject are reserved as the testing dataset, while data from the remaining subjects form the training dataset. Given that the dataset used in this study consisted of 10 subjects, this process was repeated 10 times, with each subject serving as the test dataset once. The overall performance was then calculated as the average of the results across all 10 iterations.
Fall detection outcomes are classified into two categories: ADL and fall, making this a binary classification problem. In this context, a fall is considered a positive outcome, while recognizing an ADL is considered a negative outcome. The performance of a fall detection system is quantitatively assessed by categorizing the outcomes into four types of counts as follows:
(1) True positives (TP): The number of instances where the system correctly identifies a fall when the fall has actually occurred.
(2) False positives (FP): The number of instances where the system incorrectly identifies an ADL as a fall.
(3) True negatives (TN): The number of instances where the system correctly identifies an ADL when no fall has occurred.
(4) False negatives (FN): The number of instances where the system fails to identify a fall, instead incorrectly classifying it as an ADL.
The objective of a fall detection system is to minimize the occurrence of both false positives (FP) and false negatives (FN) to enhance overall performance. Common metrics for evaluating binary classification performance include accuracy, sensitivity, specificity, precision, and the F 1 -score [35]. These metrics are widely used in the literature to assess the effectiveness of fall detection systems [15,21,22,23,24,28,34]. Higher values of these metrics generally indicate better system performance.
Accuracy is defined as the ratio of correctly classified samples (both positive and negative) to the total number of samples, as given by the following equation:
Accuracy = TP + TN TP + FP + TN + FN .
Sensitivity (or recall) represents the proportion of correctly classified positive samples (i.e., falls) to the total number of actual positive samples, and is defined as follows:
Sensitivity = TP TP + FN .
Specificity (or inverse recall) is defined as the ratio of correctly classified negative samples (i.e., ADLs) to the total number of actual negative samples, and is given by
S pecificity = TN FP + TN .
Precision represents the proportion of correctly classified positive samples to the total number of predicted positive samples, as indicated by the following equation:
Precision = TP FP + TP .
F 1 -score is the harmonic mean of sensitivity and precision. It balances the trade-off between sensitivity and precision, particularly in datasets with an imbalanced class distribution. Given that the dataset used in this study comprises 1053 ADL instances and 423 fall trials, the F 1 -score offers a balanced evaluation on classification performance by accounting for both FP and FN. The F 1 -score is defined as follows:
F 1 - score = 2 TP 2 TP + FP + FN = 2 Sensitivity Precision Sensitivity + Precision .

2.2.5. Statistical Analysis

The objective of the statistical analysis was twofold. First, its aim was to examine whether there is statistical significance between the proposed signal combination and the raw accelerometer signal for each feature under the same activity (ADL or fall). Second, its objective was to determine whether the same feature shows statistical significance between ADL and fall activities within the same signal combination. Several statistical tests were employed to achieve the desired objectives. The first test was the Shapiro–Wilk normal test, which was used to assess whether the derived features of each signal combination, under the same activity, follow a normal distribution. For the first objective, a paired-samples test was conducted. If the values of the same feature are normally distributed, a paired-samples t-test was used. Otherwise, the nonparametric Wilcoxon signed-rank test was applied. For the second objective, an independent sample test was used. If the values of the same feature followed a normal distribution, an independent sample t-test was performed. Otherwise, the nonparametric Wilcoxon rank-sum test was applied.

3. Results and Discussion

In this study, the DAGAF algorithm was employed to perform a three-level decomposition on triaxial accelerometer signals. Each signal underwent this decomposition process three times, resulting in components labeled as imf 1 (the first derived IMF), imf 2 (the second derived IMF), imf 3 (the third derived IMF), and Res (the residual component). In addition, according to the study by Liu et al. [28], the same impact-defined window, which spans 2 s before the impact (backward) and 1.23 s after the impact (forward), was adopted for the datasets in this study.
The first experiment aimed to evaluate fall detection performance by comparing the classification results obtained from various combinations of the decomposed components with those directly from the raw signals. In this experiment, fifty-four features were extracted from each signal to be analyzed (see Section 2.2.2 for details on feature extraction). Two classifiers, SVM and k NN ( k = 3 ), were employed to assess and compare their effectiveness in the fall detection task. The test results are summarized in Table 1, with the highest values highlighted in bold for each metric. The results demonstrate that combining the residual component with the raw signal or certain IMFs can enhance classification performance across both classifiers. Notably, the combinations “ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ” consistently outperformed other combinations as well as the raw signals across all tests. Specifically, the highest accuracy of 96.34% was achieved using the combination “ imf 1 + imf 3 + Res ” with the SVM classifier, which also exhibited the highest sensitivity (92.91%), specificity (98.29%), precision (95.56%), and F 1 -score (94.22%) compared to other combinations and the raw signals. Additionally, the SVM classifier consistently outperformed the k NN ( k = 3 ) classifier across all metrics. For this reason, the SVM classifier was used in subsequent computational experiments. Another notable result is that the performance of the combination “ raw + Res ” was comparable to “ imf 1 + imf 3 + Res ”. It even surpassed the latter in terms of sensitivity, precision, and F 1 -score when using the k NN ( k = 3 ) classifier. Although the computation load of GADAF decomposition is relatively lighter than that of wavelet-based spectrogram and HHT, the proposed DAGAF decomposition requires an iterative filtering procedure to extract IMFs, making its implementation on an MCU for real-time applications a challenging task. However, the good performance of “ raw + Res ” offers an alternative solution for this issue. Since Res is essentially a very low-frequency component contained in the raw signal, developing a lowpass filter that can adaptively extract a residual component similar to that in DAGAF decomposition seems to be a promising direction for future research.
These findings suggest that incorporating the residual component along with the raw signal or the specific IMFs (the first IMF and either the second or third IMF) enhances the performance in classifying ADL and fall activities. The consistent improvement across all performance metrics highlights the effectiveness of these combinations in capturing critical information embedded in the raw signals, thereby enabling more reliable fall detection.
As described previously, the combinations “ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ” consistently outperform other combinations and the raw signals across all tests. By examining the patterns of the separated components after DAGAF decomposition, valuable insights can be acquired for the information captured by these decomposed components. Figure 2 presents the three-level decomposition of triaxial accelerometer signals recorded at the waist position from subject 1 of the FallAllD dataset during (a) a daily activity (sitting down) and (b) a forward fall (walking to stumbling/tripping). This decomposition reveals the IMFs and residual component for each axis, highlighting the signal characteristics associated with both ADLs and fall activities. In this figure, the vertical axis ranges were intentionally adjusted to ensure consistency at the same decomposition level for signals from both ADLs and fall activities. Notably, a marked difference was observed in imf 1 , where the fall signals exhibited more densely oscillated patterns compared to the ADL signals, likely reflecting the abrupt and rapid movements appeared in a fall event. Similar distinct characteristics were also observed in imf 2 , imf 3 , and Res across all axes. Even for the low-frequency residual component (labeled as Res in the figure), it can be observed that the dynamic ranges during fall were greater than that during ADL activity across all axes. These differing patterns likely contribute to a more varied distribution of the extracted features in the feature space, thereby enhancing the ability of classifiers to distinguish between falls and ADLs. Moreover, as observed in Table 1, the combinations “ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ” demonstrated superior performance compared to other combinations or single decomposed components. This indicates that the richness in the oscillatory characteristics within various IMFs provides a diverse set of discriminative features, which in turn facilitates a more effective separation in the feature space when using a combination of decomposed components for the fall detection task.
Even though a clear distinction can be observed with the help of DAGAF decomposition, analyzing the patterns of accelerometer signals in fall activities and ADLs remain a complex task. Figure 3 presents the examples of TP, TN, FP, and FN by both SVM and 3NN classifiers. The triaxial accelerometer signals shown in this figure are from subject 1 of the FallAllD dataset during a forward fall (walking to stumbling/tripping) for TP, a daily activity (sitting down) for TN, a daily activity (standing up) for FP, and from subject 2 of the FallAllD dataset during a forward fall (jogging to stumbling/tripping) for FN, respectively. From this figure, it is evident that feature engineering is necessary because the activity cannot be easily identified solely based on the patterns in the selected examples.
To examine the characteristics of the selected features across different signal combinations in ADLs and fall activities, we also performed two types of statistical analysis. The first analysis was to verify whether there exists a statistical significance between the proposed signal combination and the raw accelerometer signal for each feature under the same activity. The second analysis was to determine whether the same feature shows a statistical significance between ADLs and fall activities within the same signal combination. The results of these tests are summarized in Table 2. The results of the Shapiro–Wilk normal test showed that none of the features in the proposed signal combinations or the raw accelerometer signals followed a normal distribution. Therefore, nonparametric methods were used for both tests. For comparing features between the proposed signal combination and the raw accelerometer signal within the same activity, the Wilcoxon signed-rank test was employed for statistical testing. The asterisk (*) in the table indicates a significant difference between the proposed signal combinations and the “raw” signal for the same feature in the same activity, evaluated at a significance level of 0.05 ( p < 0.05 ) using the Wilcoxon signed-rank test. For comparing features between ADL and fall activities within each signal, the Wilcoxon rank-sum test was applied. The dagger (†) in the table denotes a significant difference between ADLs and fall activities for the same feature in the same combination, evaluated at a significance level of 0.05 ( p < 0.05 ) using the Wilcoxon rank-sum test. The results in Table 2 indicate that the statistical significance is observed across different features, with each feature demonstrating statistical significance in at least one signal combination. This finding supports that the proposed signal combinations indeed introduce measurable statistical deviations in the feature characteristics.
Based on the study by Liu et al. [28], the same impact-defined window (2 s before the impact and 1.23 s after the impact) was adopted for the datasets in this study. In this context, the second experiment was designed to compare the performance of raw signals and three out-performing combinations (“ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ”) across five different time segments using the SVM classifier. The experimental results are summarized in Table 3, where the highest values for each metric across five different time segments are highlighted in bold. It can be observed that using signals of the whole length consistently achieves the highest performance across all metrics when compared to the signals of shorter segment. Furthermore, the combinations “ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ” generally enhanced classification performance across different time segments when compared to the raw signals. Among them, the combination “ raw + Res ” or “ imf 1 + imf 3 + Res ” consistently achieved the best performance across all metrics and time segments. For instance, when analyzing signals of the whole segment (length of 3.23 s), the combination “ imf 1 + imf 3 + Res ” achieved an accuracy of 96.34%, sensitivity of 92.91%, specificity of 98.29%, precision of 95.56%, and an F 1 -score of 94.22%. Interestingly, the shorter segments containing the impact (which occurs at 2 s) demonstrated slightly lower yet comparable performance to the whole segment across all metrics. However, excluding the impact from the analysis segment (i.e., using the 0~1.8 s window in the experiments) results in a significant performance drop across various metrics. In this scenario, sensitivity is the most significantly impacted metric, dropping to a low value of 52.01% in the raw triaxial accelerometer signals. The low sensitivity is nearly equivalent to a random guess in classifying ADLs and fall activities. However, sensitivity increases to 68.79% while using the combination “ raw + Res ”. This significant improvement in sensitivity underscores the advantage of using DAGAF decomposition for fall detection, even when the signal does not fully encompass the impact event. These results also suggest that including the impact event within the analysis segment is crucial for enhancing the effectiveness of fall detection.
The proposed DAGAF decomposition for triaxial accelerometer signals introduces a novel approach to feature extraction in fall detection. To evaluate its effectiveness, it is essential to compare the performance of this method with that of other studies employing machine learning classifiers for fall detection. In this regard, Table 4 presents a summary of the performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score), allowing for a direct comparison between the proposed method and other studies that have used the same database (the FallAllD dataset) for fall detection. The highest values for each metric are highlighted in bold. The results in Table 4 clearly demonstrate that the proposed DAGAF method, especially the combination “ imf 1 + imf 3 + Res ”, outperforms previous studies across all metrics. Specifically, this combination achieved an accuracy of 96.34%, a sensitivity of 92.91%, a specificity of 98.29%, a precision of 95.56%, and an F 1 -score of 94.22%. These results surpass those reported in previous studies, such as Liu et al. [24], which achieved an accuracy of 95.60% and an F 1 -score of 92.16% using an enhanced preprocessing technique on triaxial accelerometer signals. The studies referenced in [15,21,22,23] did not report all performance metrics. Among these, the most competitive result is the sensitivity of 92.90% obtained by Ramón et al. [23] using the one-class k NN (OC- k NN) classifier trained with the majority data from ADLs, which is only slightly lower than the sensitivity achieved by the proposed approach. Please note that only the result covering the impact-window (denoted as W1) is presented in Table 4 for Šeketa et al. [22].
In summary, the results clearly demonstrate that the proposed DAGAF decomposition consistently outperforms conventional raw signal-based approaches as well as other related studies across various performance metrics. This indicates that DAGAF decomposition significantly enhances the reliability of fall detection systems. Moreover, the superior performance of the proposed method highlights its potential for practical integration into wearable devices, where accurate and reliable fall detection is crucial.

4. Conclusions

This study presents a novel approach for fall detection by applying DAGAF decomposition to triaxial accelerometer signals in combination with machine learning classifiers. Although conceptually similar to EMD, the proposed DAGAF algorithm stands out for its unique use of Gaussian average filtering to compute the instantaneous mean. The length of Gaussian filter is dynamically adjusted based on the number of extrema in the analyzed signal, making the proposed method a data-adaptive decomposition. While the DAGAF algorithm has been successfully applied in various biomedical scenarios [25], this study marks its first application in the domain of fall detection. To facilitate further research, the detailed mathematical formulation of the DAGAF algorithm is provided, along with the codes implemented in both MATLAB and Python, which are accessible on GitHub [29].
The three-level DAGAF decomposition was consistently applied throughout all computational experiments in this study. This choice was driven by its applicability to all signals in the FallAllD dataset, thereby ensuring consistency across all analyses. After applying the three-level DAGAF decomposition, each signal was separated into components labeled as imf 1 (the first derived IMF), imf 2 (the second derived IMF), imf 3 (the third derived IMF), and Res (the residual component). The same impact-defined window used by Liu et al. [28], which covers 2 s before the impact (backward) and 1.23 s after the impact (forward), was adopted for the dataset in this study. To evaluate the effectiveness of the proposed algorithm in fall detection, computational experiments were conducted comparing the classification performance across various component combinations against the raw signals (refer to Table 1). The results indicate that specific component combinations significantly improve classification performance in fall detection. Specifically, the combination of imf 1 , imf 3 , and Res with an SVM classifier achieves the highest performance metrics, with an accuracy of 96.34%, a sensitivity of 92.91%, a specificity of 98.29%, a precision of 95.56%, and an F 1 -score of 94.22%. This combination outperforms all other tested combinations, including the use of raw signals.
This study also examined the characteristics of the selected features across different signal combinations in ADLs and fall activities using two types of statistical analysis. The first analysis assessed whether there exists a statistical significance between the proposed signal combination and the raw accelerometer signal for each feature within the same activity. The second analysis evaluated whether the same feature exhibits a statistical significance between ADLs and fall activities within the same signal combination. Nonparametric methods were used in both statistical tests since none of the features in the proposed signal combinations or the raw accelerometer signals followed a normal distribution. The Wilcoxon signed-rank test was used to compare features between the proposed signal combination and the raw accelerometer signal within the same activity, while the Wilcoxon rank-sum test was employed to compare features between ADLs and fall activities within each signal. The results of these tests are summarized in Table 2. The table shows that statistical significance was observed across different features, with each feature demonstrating statistical significance in at least one signal combination. This finding indicates that the proposed signal combinations introduce statistical deviations in the feature characteristics.
To assess the impact of selected time segments on fall detection, computational experiments were conducted across five different time segments using the SVM classifier. The performance of the two top-performing combinations (“ imf 1 + imf 2 + Res ”, “ raw + Res ”, and “ imf 1 + imf 3 + Res ”) were compared against the raw signals in classifying ADLs and fall activities (see Table 3). It was found that the combination “ raw + Res ” or “ imf 1 + imf 3 + Res ” consistently outperformed in all performance metrics across all time segments. The signals of the whole segment (length of 3.23 s) consistently achieved higher performance metrics than those of the shorter segment. Moreover, the performance dropped significantly if the impact event was not included in the signal segment to be analyzed (0~1.8 s window in the experiments); the sensitivity dropped to 52.01% using the raw signals. However, the sensitivity increased to 68.79% for the combination “ raw + Res ” in this scenario.
As compared with other studies using the same database, the combination “ imf 1 + imf 3 + Res ” with the SVM classifier outperforms previous studies across all performance metrics (refer to Table 4). These findings suggest that the DAGAF decomposition not only improves the feature extraction process but also provides a promising advancement in the field of fall detection. This research holds significant potential for advancing the development of effective and reliable wearable monitoring systems, which are crucial for improving the safety monitoring on elderly individuals.
From Table 1, it is noteworthy that the performance of the combination “ raw + Res ” is comparable to “ imf 1 + imf 3 + Res ”. It even surpasses the latter in terms of sensitivity, precision, and F 1 -score when using the k NN ( k = 3 ) classifier. Although the computation load of DAGAF decomposition is lower than that of wavelet-based spectrogram and HHT, the DAGAF algorithm requires an iterative filtering process to extract IMFs, which makes its implementation on an MCU challenging for real-time applications. However, the good performance of “ raw + Res ” offers a viable solution to this challenge. Since Res is essentially a very low-frequency component within the raw signal, developing a lowpass filter that can adaptively extract a residual component similar to that in DAGAF is a promising direction for future research. In addition, the proposed DAGAF decomposition has potential to be incorporated in the analysis of multimodal data sources to further improve the accuracy and reliability of fall detection tasks. Using multiple time segment windows, similar to how [22] was conducted, in combination with DAGAF decomposition is also another promising direction for further development.

Author Contributions

Conceptualization, Y.-D.L.; methodology, Y.-D.L., C.-J.L., M.-H.S. and J.-H.H.; software, Y.-D.L., C.-J.L., M.-H.S. and J.-H.H.; validation, C.-J.L., M.-H.S. and J.-H.H.; formal analysis, C.-J.L., M.-H.S. and J.-H.H.; investigation, C.-J.L., M.-H.S. and J.-H.H.; resources, Y.-D.L.; data curation, C.-J.L., M.-H.S. and J.-H.H.; writing—original draft preparation, Y.-D.L.; writing—review and editing, Y.-D.L.; visualization, Y.-D.L., C.-J.L., M.-H.S. and J.-H.H.; supervision, Y.-D.L.; project administration, Y.-D.L.; funding acquisition, Y.-D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science and Technology Council, Taiwan (contract number: NSTC 110-2221-E-035-006-MY3 and NSTC 113-2221-E-035-011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The FallAllD dataset is available on IEEE DataPort™ [26], and the Python and MATLAB codes for DAGAF decomposition are accessible on GitHub [29].

Acknowledgments

The authors would like to thank Kai-Chun Liu for his kind sharing of his related knowledge and his experience in fall detection. This article could not have been finished without his help.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. World Health Organization. Falls. Available online: https://www.who.int/news-room/fact-sheets/detail/falls (accessed on 24 September 2024).
  2. Montero-Odasso, M.; Van Der Velde, N.; Martin, F.C.; Petrovic, M.; Tan, M.P.; Ryg, J.; Aguilar-Navarro, S.; Alexander, N.B.; Becker, C.; Blain, H. World guidelines for falls prevention and management for older adults: A global initiative. Age Ageing 2022, 51, afac205. [Google Scholar] [CrossRef] [PubMed]
  3. World Health Organization. WHO Global Report on Falls Prevention in Older Age. Available online: https://www.who.int/publications/i/item/9789241563536 (accessed on 24 September 2024).
  4. Ganz, D.A.; Latham, N.K. Prevention of falls in community-dwelling older adults. N. Engl. J. Med. 2020, 382, 734–743. [Google Scholar] [CrossRef] [PubMed]
  5. Mubashir, M.; Shao, L.; Seed, L. A survey on fall detection: Principles and approaches. Neurocomputing 2013, 100, 144–152. [Google Scholar] [CrossRef]
  6. Islam, M.M.; Tayan, O.; Islam, M.R.; Islam, M.S.; Nooruddin, S.; Kabir, M.N.; Islam, M.R. Deep learning based systems developed for fall detection: A review. IEEE Access 2020, 8, 166117–166137. [Google Scholar] [CrossRef]
  7. Usmani, S.; Saboor, A.; Haris, M.; Khan, M.A.; Park, H. Latest research trends in fall detection and prevention using machine learning: A systematic review. Sensors 2021, 21, 5134. [Google Scholar] [CrossRef]
  8. Principi, E.; Droghini, D.; Squartini, S.; Olivetti, P.; Piazza, F. Acoustic cues from the floor: A new approach for fall classification. Expert Syst. Appl. 2016, 60, 51–61. [Google Scholar] [CrossRef]
  9. Wang, H.; Zhang, D.; Wang, Y.; Ma, J.; Wang, Y.; Li, S. RT-Fall: A real-time and contactless fall detection system with commodity WiFi devices. IEEE Trans. Mob. Comput. 2017, 16, 511–526. [Google Scholar] [CrossRef]
  10. Hu, Y.; Zhang, F.; Wu, C.; Wang, B.; Liu, K.R. DeFall: Environment-independent passive fall detection using WiFi. IEEE Internet Things J. 2021, 9, 8515–8530. [Google Scholar] [CrossRef]
  11. Frøvik, N.; Malekzai, B.A.; Øvsthus, K. Utilising LiDAR for fall detection. Healthc. Technol. Lett. 2021, 8, 11–17. [Google Scholar] [CrossRef]
  12. Piñeiro, M.; Araya, D.; Ruete, D.; Taramasco, C. Low-cost LIDAR-based monitoring system for fall detection. IEEE Access 2024, 12, 72051–72061. [Google Scholar] [CrossRef]
  13. Rezaei, A.; Mascheroni, A.; Stevens, M.C.; Argha, R.; Papandrea, M.; Puiatti, A.; Lovell, N.H. Unobtrusive human fall detection system using mmwave radar and data driven methods. IEEE Sens. J. 2023, 23, 7968–7976. [Google Scholar] [CrossRef]
  14. Qi, P.; Chiaro, D.; Piccialli, F. FL-FD: Federated learning-based fall detection with multimodal data fusion. Inf. Fusion 2023, 99, 101890. [Google Scholar] [CrossRef]
  15. Saleh, M.; Abbas, M.; Le Jeannès, R.B. FallAllD: An open dataset of human falls and activities of daily living for classical and deep learning applications. IEEE Sens. J. 2020, 21, 1849–1858. [Google Scholar] [CrossRef]
  16. Palmerini, L.; Bagalà, F.; Zanetti, A.; Klenk, J.; Becker, C.; Cappello, A. A wavelet-based approach to fall detection. Sensors 2015, 15, 11575–11586. [Google Scholar] [CrossRef] [PubMed]
  17. Tocco, F.; Solinas, R.; Velluzzi, F.; Massidda, M.; Mattana, D.V.; Fois, A.; Melis, L.; Bertetto, A.M.; Bonisoli, E.; Venturini, S. A mechatronic tool for revealing inverse relationships among heart’s stroke volume and head’s linear acceleration induced by moored boats rolling in elderly sailors with unchamged body sizes: A non-drug anti-hypertensive advantage? Int. J. Mech. Control 2024, 25, 133–142. [Google Scholar] [CrossRef]
  18. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  19. Erfianto, B.; Rizal, A.; Hadiyoso, S. Empirical mode decomposition and Hilbert spectrum for abnormality detection in normal and abnormal walking transitions. Int. J. Environ. Res. Public Health 2023, 20, 3879. [Google Scholar] [CrossRef]
  20. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  21. Silva, C.A.; García− Bermúdez, R.; Casilari, E. Features selection for fall detection systems based on machine learning and accelerometer signals. In Proceedings of the Advances in Computational Intelligence: 16th International Work-Conference on Artificial Neural Networks, IWANN 2021, Virtual Event, 16–18 June 2021; pp. 380–391. [Google Scholar]
  22. Šeketa, G.; Pavlaković, L.; Džaja, D.; Lacković, I.; Magjarević, R. Event-centered data segmentation in accelerometer-based fall detection algorithms. Sensors 2021, 21, 4335. [Google Scholar] [CrossRef]
  23. Santoyo-Ramón, J.A.; Casilari, E.; Cano-García, J.M. A study of one-class classification algorithms for wearable fall sensors. Biosensors 2021, 11, 284. [Google Scholar] [CrossRef]
  24. Liu, K.-C.; Hung, K.-H.; Hsieh, C.-Y.; Huang, H.-Y.; Chan, C.-T.; Tsao, Y. Deep-learning-based signal enhancement of low-resolution accelerometer for fall detection systems. IEEE Trans. Cogn. Dev. Syst. 2021, 14, 1270–1281. [Google Scholar] [CrossRef]
  25. Lin, Y.-D.; Tan, Y.K.; Tian, B. A novel approach for decomposition of biomedical signals in different applications based on data-adaptive Gaussian average filtering. Biomed. Signal Process. Control 2022, 71, 103104. [Google Scholar] [CrossRef]
  26. Saleh, M.; Abbas, M.; Le Jeannès, R.B. FallAllD: An Open Dataset of Human Falls and Activities of Daily Living for Classical and Deep Learning Applications. Available online: https://ieee-dataport.org/open-access/fallalld-comprehensive-dataset-human-falls-and-activities-daily-living (accessed on 24 September 2024).
  27. Özdemir, A.T. An analysis on sensor locations of the human body for wearable fall detection devices: Principles and practice. Sensors 2016, 16, 1161. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, K.-C.; Hsieh, C.-Y.; Huang, H.-Y.; Hsu, S.J.-P.; Chan, C.-T. An analysis of segmentation approaches and window sizes in wearable-based critical fall detection systems with machine learning models. IEEE Sens. J. 2020, 20, 3303–3313. [Google Scholar] [CrossRef]
  29. Liu, K.-C.; Lin, Y.-D. DAGAF. Available online: https://github.com/t22302856/DAGAF (accessed on 3 August 2024).
  30. Feng, Z.; Zhang, D.; Zuo, M.J. Adaptive mode decomposition methods and their applications in signal analysis for machinery fault diagnosis: A review with examples. IEEE Access 2017, 5, 24301–24331. [Google Scholar] [CrossRef]
  31. Harris, F.J. On the use of windows for harmonic analysis with the discrete Fourier transform. Proc. IEEE 1978, 66, 51–83. [Google Scholar] [CrossRef]
  32. Cicone, A.; Liu, J.; Zhou, H. Adaptive local iterative filtering for signal decomposition and instantaneous frequency analysis. Appl. Comput. Harmon. Anal. 2016, 41, 384–411. [Google Scholar] [CrossRef]
  33. Rato, R.; Ortigueira, M.D.; Batista, A. On the HHT, its problems, and some solutions. Mech. Syst. Signal Process. 2008, 22, 1374–1394. [Google Scholar] [CrossRef]
  34. Hsieh, C.-Y.; Liu, K.-C.; Huang, C.-N.; Chu, W.-C.; Chan, C.-T. Novel hierarchical fall detection algorithm using a multiphase fall model. Sensors 2017, 17, 307. [Google Scholar] [CrossRef]
  35. Tharwat, A. Classification assessment methods. Appl. Comput. Inform. 2021, 17, 168–192. [Google Scholar] [CrossRef]
Figure 1. Iterative procedure of three-level decomposition using the proposed DAGAF decomposition algorithm, demonstrated with the z-axis accelerometer signal (jogging to stumbling/tripping) collected at the waist position from subject 1 of the FallAllD dataset. In each subfigure, the solid line (in black color) shows the raw accelerometer signal (in Step 1) or the IMF (in Steps 2 and 3), while the dash-dotted line (in three different colors) indicates the instantaneous mean sifted by a Gaussian window.
Figure 1. Iterative procedure of three-level decomposition using the proposed DAGAF decomposition algorithm, demonstrated with the z-axis accelerometer signal (jogging to stumbling/tripping) collected at the waist position from subject 1 of the FallAllD dataset. In each subfigure, the solid line (in black color) shows the raw accelerometer signal (in Step 1) or the IMF (in Steps 2 and 3), while the dash-dotted line (in three different colors) indicates the instantaneous mean sifted by a Gaussian window.
Information 15 00606 g001
Figure 2. Three-level decomposition of triaxial accelerometer signals collected at the waist position from subject 1 of the FallAllD dataset during (a) a daily activity (sitting down) and (b) a forward fall (walking to stumbling/tripping). In both (a,b), each subfigure (from top to bottom) presents the raw accelerometer signal (solid black line), three IMFs, and the residual component (labeled as Res) after three-level decomposition for the x-, y-, and z-axis signals (from left to right). The top subfigures in both (a,b) also display the patterns for the combinations of “ imf 1 + imf 2 + Res ” (blue dashed line), “ raw + Res ” (green dotted line), and “ imf 1 + imf 3 + Res ” (magenta dash-dotted line).
Figure 2. Three-level decomposition of triaxial accelerometer signals collected at the waist position from subject 1 of the FallAllD dataset during (a) a daily activity (sitting down) and (b) a forward fall (walking to stumbling/tripping). In both (a,b), each subfigure (from top to bottom) presents the raw accelerometer signal (solid black line), three IMFs, and the residual component (labeled as Res) after three-level decomposition for the x-, y-, and z-axis signals (from left to right). The top subfigures in both (a,b) also display the patterns for the combinations of “ imf 1 + imf 2 + Res ” (blue dashed line), “ raw + Res ” (green dotted line), and “ imf 1 + imf 3 + Res ” (magenta dash-dotted line).
Information 15 00606 g002aInformation 15 00606 g002b
Figure 3. Examples for TP, TN, FP, and FN (from top to bottom) by both SVM and 3NN classifiers. The triaxial accelerometer signals of the examples are from subject 1 of the FallAllD dataset during a forward fall (walking to stumbling/tripping) for TP, a daily activity (sitting down) for TN, a daily activity (standing up) for FP, and from Subject 2 of the FallAllD dataset during a forward fall (jogging to stumbling/tripping) for FN, respectively. The signals collected at the waist position are used in the examples.
Figure 3. Examples for TP, TN, FP, and FN (from top to bottom) by both SVM and 3NN classifiers. The triaxial accelerometer signals of the examples are from subject 1 of the FallAllD dataset during a forward fall (walking to stumbling/tripping) for TP, a daily activity (sitting down) for TN, a daily activity (standing up) for FP, and from Subject 2 of the FallAllD dataset during a forward fall (jogging to stumbling/tripping) for FN, respectively. The signals collected at the waist position are used in the examples.
Information 15 00606 g003
Table 1. Performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) for SVM and k NN classifiers across various signal combinations.
Table 1. Performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) for SVM and k NN classifiers across various signal combinations.
Signals SVM k NN   ( k = 3 )
acc sen spe pre F 1 acc sen spe pre F 1
raw 95.80 90.78 97.82 94.35 92.53 94.51 91.02 95.92 89.95 90.48
imf 1 88.69 76.83 93.45 82.49 79.56 87.47 78.96 90.88 77.67 78.31
imf 2 92.68 85.58 95.54 88.51 87.02 91.46 83.92 94.49 85.96 84.93
imf 3 90.65 82.98 93.73 84.17 83.57 90.72 84.16 93.35 83.57 83.86
Res 93.16 86.29 95.92 89.46 87.85 92.62 86.29 95.16 87.74 87.01
imf 1 + imf 2 91.46 82.74 94.97 86.85 84.75 91.26 84.16 94.11 85.17 84.66
imf 2 + imf 3 92.28 86.76 94.49 86.35 86.55 92.01 86.05 94.40 86.05 86.05
imf 1 + imf 3 90.38 79.43 94.78 89.53 82.55 90.79 82.74 94.02 84.75 83.73
imf 1 + Res 95.46 90.07 97.63 93.84 91.92 94.51 88.89 96.58 90.95 89.91
imf 2 + Res 95.26 90.07 97.34 93.15 91.58 94.51 89.13 96.68 91.50 90.30
imf 3 + Res 94.99 88.65 97.53 93.52 91.02 94.24 88.65 96.39 91.02 89.82
imf 1 + imf 2 + imf 3 92.07 83.92 95.35 87.87 85.85 91.73 84.87 94.49 89.09 86.93
imf 2 + imf 3 + Res 95.46 90.07 97.63 93.84 91.92 94.11 89.60 95.92 89.81 89.70
imf 1 + imf 2 + Res 95.80 91.49 97.72 94.72 93.0894.65 90.31 96.77 91.61 90.96
raw + Res 96.21 91.02 98.29 95.53 93.22 95.39 91.73 96.87 92.16 91.94
imf 1 + imf 3 + Res 96.34 92.91 98.29 95.56 94.2295.53 91.02 96.99 91.71 91.36
Legends: (1) “raw” represents the raw triaxial accelerometer signals collected at the waist position. (2) “ imf i ” (i is 1, 2, or 3) denotes the i-th decomposed IMF, respectively, from the triaxial accelerometer signals using the proposed DAGAF decomposition algorithm, whereas “Res” is the residual component after three IMFs have been decomposed. (3) “acc”, “sen”, “spe”, “pre”, and “ F 1 ” are the abbreviations for accuracy, sensitivity, specificity, precision, and F 1 -score, respectively. (4) All of the values are shown in percentage (%).
Table 2. Statistical analysis for 54 features in ADLs and fall activities across the proposed signal combinations.
Table 2. Statistical analysis for 54 features in ADLs and fall activities across the proposed signal combinations.
Feature ADL Fall
raw comb 1 raw + Res comb 2 raw comb 1 raw + Res comb 2
mean ( s x ) 3504.29 ± 1042.75 3503.85 ± 1043.21 700.78 ± 2085.44 *3504.12 ± 1042.511948.93 ± 959.09 1950.05 ± 957.52 3899.96 ± 1918.00 *1949.59 ± 960.26
mean ( s y ) −219.14 ± 1138.26 −218.40 ± 1137.37 −437.49 ± 2275.70 *−219.14 ± 1138.3659.93 ± 1347.92 60.82 ± 1349.06 120.77 ± 2697.32 59.91 ± 1348.27
mean ( s z ) 323.37 ± 893.38 323.27 ± 891.83 646.44 ± 1785.05 *323.22 ± 893.32−193.63 ± 1329.81 −193.98 ± 1329.18 −387.45 ± 2658.92 *−193.50 ± 1329.75
mean ( Norm ) 4269.45 ± 341.78 4200.36 ± 312.91 *8205.06 ± 205.63 *4151.02 ± 151.58 *4416.80 ± 202.35 4213.22 ± 202.15 *8233.10 ± 230.63 *4222.65 ± 99.88 *
mean ( Hori ) 1709.89 ± 869.67 1607.53 ± 861.58 *2977.65 ± 1783.87 *1565.76 ± 883.33 *2805.04 ± 651.61 2643.61 ± 629.21 *5097.20 ± 1314.20 *2632.97 ± 665.63 *
mean ( Vert ) 4036.98 ± 461.99 3989.90 ± 464.76 *7816.40 ± 720.46 *3940.04 ± 364.36 *3679.13 ± 560.79 3516.78 ± 563.90 *6869.80 ± 1088.86 *3515.36 ± 544.05 *
std ( s x ) 1353.28 ± 1145.44 1190.55 ± 1092.13 *1533.13 ± 1204.68 *910.87 ± 729.31 *2491.16 ± 788.09 2229.85 ± 746.97 *3957.79 ± 1256.71 *2111.31 ± 658.33 *
std ( s y ) 1005.24 ± 628.32 882.19 ± 587.25 *1318.10 ± 937.82 *770.77 ± 531.67 *2077.83 ± 627.27 1848.58 ± 614.00 *3167.59 ± 1208.14 *1703.18 ± 618.32 *
std ( s z ) 810.28 ± 514.33 660.27 ± 481.05 *1118.76 ± 825.78 *664.44 ± 461.63 *2045.12 ± 678.39 1802.99 ± 632.38 *3184.20 ± 1225.53 *1742.81 ± 633.91 *
std ( Norm ) 1259.15 ± 1070.41 1095.81 ± 1029.02 *1385.21 ± 1216.81 *835.04 ± 740.38 *2103.03 ± 725.14 1735.41 ± 702.96 *2286.82 ± 782.81 *1370.20 ± 468.99 *
std ( Hori ) 804.21 ± 492.36 696.62 ± 465.68 *1033.86 ± 630.42 *621.29 ± 396.20 *2094.33 ± 534.76 1850.64 ± 499.75 *3042.34 ± 792.11 *1665.32 ± 432.32 *
std ( Vert ) 1286.82 ± 1043.17 1134.79 ± 1006.16 *1478.81 ± 1196.52 *878.19 ± 718.55 *1936.37 ± 678.10 1651.49 ± 663.08 *2407.55 ± 828.55 *1359.60 ± 448.53 *
var ( s x ) (3142.16 ± 5580.57) × 103(2609.02 ± 5171.44) × 103 *(3800.37 ± 5999.48) × 103 *(1361.08 ± 2415.73) × 103 *(6825.47 ± 4240.33) × 103 (5528.90 ± 3717.09) × 103 *(17,239.72 ± 9937.40) × 103 *(4890.00 ± 2750.42) × 103 *
var ( s y ) (1404.93 ± 2008.09) × 103(1122.79 ± 1697.03) × 103 *(2616.06 ± 4355.61) × 103 *(876.49 ± 1354.46) × 103 *(4709.92 ± 2800.99) × 103 (3793.34 ± 2461.39) × 103 *(11,489.80 ± 8806.44) × 103 *(3282.25 ± 2377.06) × 103 *
var ( s z ) (920.84 ± 1371.21) × 103(667.14 ± 1168.61) × 103 *(1932.88 ± 3832.36) × 103 *(654.38 ± 1113.51) × 103 *(4641.64 ± 3008.49) × 103 (3649.73 ± 2513.84) × 103 *(11,637.50 ± 9060.17) × 103 *(3438.28 ± 2531.63) × 103 *
var ( Norm ) (2730.16 ± 4692.12) × 103(2258.67 ± 4368.69) × 103 *(3398.01 ± 6160.73) × 103 *(1244.92 ± 2416.51) × 103 *(4947.31 ± 3576.56) × 103 (3504.61 ± 3139.28) × 103 *(5840.89 ± 4349.90) × 103 *(2096.87 ± 1486.17) × 103 *
var ( Hori ) (888.95 ± 1291.45) × 103(701.93 ± 1090.62) × 103 *(1465.92 ± 1925.45) × 103 *(542.82 ± 771.32) × 103 *(4671.50 ± 2307.95) × 103 (3674.04 ± 1900.01) × 103 *(9881.77 ± 4651.09) × 103 *(2959.73 ± 1428.19) × 103 *
var ( Vert ) (2743.08 ± 4598.80) × 103(2299.15 ± 4313.69) × 103 *(3617.18 ± 6058.64) × 103 *(1287.04 ± 2349.72) × 103 *(4208.28 ± 3305.33) × 103 (3166.06 ± 3004.31) × 103 *(6476.22 ± 4666.46) × 103 *(2049.22 ± 1411.53) × 103 *
max ( s x ) 9047.06 ± 5639.30 8557.67 ± 5374.96 *12,744.19 ± 6011.98 *7536.47 ± 4406.26 *12,221.09 ± 5692.71 11,033.52 ± 5290.21 *15,177.50 ± 6026.64 *9340.90 ± 4308.40 *
max ( s y ) 2831.59 ± 2968.11 2583.76 ± 2829.12 *3026.86 ± 3792.53 *2115.59 ± 2550.40 *9890.79 ± 6781.39 9066.59 ± 6321.60 *11,098.97 ± 7896.69 *7228.08 ± 5337.65 *
max ( s z ) 2793.95 ± 2248.62 2521.34 ± 2203.06 *3523.89 ± 2833.92 *2275.51 ± 1928.96 *8717.96 ± 7369.19 8136.75 ± 6675.19 *9452.84 ± 8361.46 *6534.16 ± 5463.11 *
max ( Norm ) 10,340.01 ± 5842.69 9703.58 ± 5552.35 *14,352.34 ± 5802.07 *8547.32 ± 4562.59 *22,998.50 ± 7458.42 20,639.10 ± 6975.41 *27,101.05 ± 7485.76 *17,059.38 ± 5690.31 *
max ( Hori ) 5484.36 ± 3858.98 5017.62 ± 3618.10 *7033.17 ± 4164.83 *4405.15 ± 3109.22 *20,024.90 ± 7858.64 17,908.54 ± 7230.11 *23,636.80 ± 8400.44 *14,771.54 ± 5950.01 *
max ( Vert ) 9965.43 ± 5647.65 9378.83 ± 5367.82 *13,910.36 ± 5659.37 *8247.72 ± 4377.37 *18,456.92 ± 6439.53 16,670.26 ± 6179.58 *21,733.75 ± 6572.93 *13,679.53 ± 5030.26 *
min ( s x ) 381.84 ± 2415.20 569.49 ± 2295.26 3573.95 ± 3219.18 *1201.07 ± 2034.64 *−6454.05 ± 4912.55 −5651.08 ± 4550.60 *−6421.52 ± 5833.33 −4546.28 ± 3824.39 *
min ( s y ) −3517.72 ± 3039.80 −3211.51 ± 2893.10 *−4162.69 ± 3696.03 *−2692.58 ± 2454.99 *−9086.35 ± 7007.68 −8393.14 ± 6355.89 *−9988.98 ± 8463.90 *−6465.27 ± 5418.55 *
min ( s z ) −2412.30 ± 3072.41 −2104.12 ± 2895.14 *−2510.63 ± 3661.93−1816.97 ± 2613.41 *−10,726.19 ± 8518.68 −9598.53 ± 7598.53 *−12,439.80 ± 10185.63 *−8098.34 ± 6551.50 *
min ( Norm ) 1870.10 ± 1100.73 1917.83 ± 1115.63 5175.05 ± 1982.90 *2328.19 ± 1113.32 *1117.59 ± 668.78 1176.69 ± 647.77 *3535.62 ± 1463.73 *1618.60 ± 726.18 *
min ( Hori ) 388.88 ± 616.35 433.36 ± 645.47 *1093.97 ± 1459.72 *504.67 ± 709.95 *214.17 ± 287.13255.92 ± 319.72 *758.07 ± 939.51 *345.71 ± 465.32 *
min ( Vert ) 1558.98 ± 1131.97 1609.05 ± 1156.71 4623.38 ± 2159.88 *2004.20 ± 1183.65 *477.28 ± 484.51 487.23 ± 461.89 1789.39 ± 1491.02 *790.09 ± 690.67 *
range ( s x ) 8665.22 ± 7357.40 7988.18 ± 6977.40 *9170.25 ± 7444.89 *6335.40 ± 5684.58 *18,675.13 ± 7402.31 16,684.60 ± 7019.46 *21,599.02 ± 7574.02 *13,887.18 ± 5690.56 *
range ( s y ) 6349.31 ± 4774.01 5795.27 ± 4613.56 *7189.56 ± 4973.51 *4808.17 ± 3916.64 *18,977.14 ± 8026.56 17,459.73 ± 7952.71 *21,087.95 ± 7951.05 *13,693.35 ± 6403.44 *
range ( s z ) 5206.25 ± 4397.04 4625.46 ± 4315.33 *6034.52 ± 4597.18 *4092.47 ± 3771.19 *19,444.15 ± 9267.60 17,735.29 ± 8949.08 *21,892.64 ± 9417.59 *14,632.50 ± 7350.06 *
range ( Norm ) 8469.92 ± 6688.76 7785.75 ± 6423.03 *9177.29 ± 7523.98 *6219.13 ± 5460.71 *21,880.91 ± 7606.79 19,462.42 ± 7194.53 *23,565.43 ± 8009.41 *15,440.77 ± 5971.52 *
range ( Hori ) 5095.48 ± 3936.14 4584.26 ± 3714.41 *5939.21 ± 4208.27 *3900.48 ± 3209.34 *19,810.73 ± 7867.74 17,652.62 ± 7245.98 *22,878.74 ± 8400.39 *14,425.83 ± 5965.70 *
range ( Vert ) 8406.46 ± 6456.81 7769.78 ± 6194.49 *9286.98 ± 7331.84 *6243.52 ± 5257.74 *17,979.64 ± 6520.75 16,183.03 ± 6301.68 *19,944.36 ± 6891.24 *12,889.45 ± 5216.05 *
kurt ( s x ) 6.62 ± 6.10 7.54 ± 7.50 *5.63 ± 4.81 *6.87 ± 5.95 *7.89 ± 10.207.61 ± 10.50 *4.07 ± 6.28 *5.72 ± 8.43 *
kurt ( s y ) 5.38 ± 6.56 5.95 ± 8.29 *4.56 ± 5.06 *5.66 ± 7.9913.87 ± 14.08 14.61 ± 17.11 7.13 ± 8.65 *10.60 ± 14.87 *
kurt ( s z ) 5.37 ± 7.09 6.51 ± 10.31 *4.22 ± 4.02 *5.02 ± 6.93 *16.23 ± 16.28 16.32 ± 17.65 *7.41 ± 8.19 *10.12 ± 11.02 *
kurt ( Norm ) 9.87 ± 11.68 11.15 ± 13.42 *8.09 ± 7.82 *10.17 ± 13.12 *37.50 ± 25.70 45.43 ± 31.05 *30.93 ± 22.85 *38.13 ± 28.53
kurt ( Hori ) 7.40 ± 9.43 8.31 ± 10.74 *5.57 ± 5.50 *7.25 ± 9.96 *24.71 ± 19.25 24.10 ± 21.57 *10.96 ± 11.68 *16.45 ± 17.01 *
kurt ( Vert ) 8.29 ± 8.70 9.28 ± 10.69 *6.62 ± 6.11 *8.18 ± 8.58 *24.50 ± 21.29 27.12 ± 23.86 *15.24 ± 15.64 *21.19 ± 20.35 *
skew ( s x ) 0.75 ± 1.02 0.75 ± 1.06 0.63 ± 0.92 *0.79 ± 1.02 *0.27 ± 1.36 0.21 ± 1.28 *−0.03 ± 0.94 *0.09 ± 1.09 *
skew ( s y ) −0.18 ± 0.93 −0.16 ± 0.94 *−0.16 ± 0.85 *−0.16 ± 0.97 *0.23 ± 2.14 0.24 ± 2.07 0.19 ± 1.38 0.22 ± 1.67
skew ( s z ) −0.07 ± 0.99 −0.10 ± 1.05 *−0.05 ± 0.82 *−0.05 ± 0.89 *−0.23 ± 2.54 −0.16 ± 2.46 *−0.15 ± 1.49 *−0.19 ± 1.72
skew ( Norm ) 1.42 ± 1.37 1.38 ± 1.47 1.01 ± 1.17 *1.35 ± 1.354.55 ± 2.01 5.01 ± 2.22 *3.70 ± 2.00 *4.27 ± 2.10 *
skew ( Hori ) 1.24 ± 1.17 1.29 ± 1.28 *0.83 ± 1.00 *1.07 ± 1.27 *3.25 ± 1.68 3.02 ± 1.78 *1.62 ± 1.23 *2.18 ± 1.52 *
skew ( Vert ) 1.19 ± 1.21 1.17 ± 1.29 0.95 ± 1.00 *1.11 ± 1.143.14 ± 1.94 3.15 ± 2.13 1.77 ± 1.78 *2.42 ± 1.97 *
corr ( s x ,   s y ) −0.03 ± 0.38 −0.02 ± 0.39 *−0.02 ± 0.41 *−0.02 ± 0.41 *0.05 ± 0.44 0.04 ± 0.49 *0.04 ± 0.600.04 ± 0.55
corr ( s x ,   s z ) −0.01 ± 0.37 −0.02 ± 0.38 *−0.04 ± 0.41 *−0.03 ± 0.38 *0.00 ± 0.430.02 ± 0.49 *−0.01 ± 0.590.01 ± 0.54
corr ( s y ,   s z ) 0.07 ± 0.45 0.07 ± 0.46 0.05 ± 0.48 *0.05 ± 0.47 *0.05 ± 0.430.06 ± 0.460.03 ± 0.53 *0.04 ± 0.49 *
corr ( Norm ,   Hori ) 0.50 ± 0.26 0.46 ± 0.27 *0.41 ± 0.30 *0.43 ± 0.29 *0.76 ± 0.14 0.72 ± 0.15 *0.63 ± 0.18 *0.67 ± 0.16 *
corr ( Norm ,   Vert ) 0.93 ± 0.18 0.92 ± 0.18 *0.91 ± 0.21 *0.90 ± 0.21 *0.84 ± 0.15 0.81 ± 0.19 *0.73 ± 0.24 *0.76 ± 0.22 *
corr ( Hori ,   Vert ) 0.35 ± 0.35 0.30 ± 0.37 *0.23 ± 0.40 *0.24 ± 0.39 *0.46 ± 0.29 0.36 ± 0.34 *0.16 ± 0.39 *0.24 ± 0.35 *
Legends: “raw” represents the raw triaxial accelerometer signals collected at the waist position. “Res” is the residual component after three IMFs have been decomposed. (1) “comb 1” is a brief notation for the combination “ imf 1 + imf 2 + Res ”, and “comb 2” is used for the combination “ imf 1 + imf 3 + Res ”. (2) “std” refers to the computation of standard deviation, whereas “var” refers to the computation of variance. (3) “max” represents the computation of the maximum value, whereas “min” represents the computation of the minimum value. (4) “range” denotes the difference between the maximum and minimum values. (5) “kurt” refers to the computation of kurtosis, and “skew” refers to the computation of skewness. (6) “corr” represents the computation of Pearson’s correlation coefficient. (7) All of the values are presented as mean ± S . D .   ( standard   deviation ) . The accelerometer signals are analyzed in their original digital forms without conversion to acceleration unit (m/s2). Some big values are presented in scientific notation. (8) The asterisk (*) indicates a significant difference between the proposed combinations and the “raw” signal for the same feature in the same activity, evaluated at a significance level of 0.05 ( p < 0.05 ) using the Wilcoxon signed-rank test. (9) The dagger () indicates a significant difference between fall activities and ADLs for the same feature in the same combination, evaluated at a significance level of 0.05 ( p < 0.05 ) using the Wilcoxon rank-sum test.
Table 3. Performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) for SVM classifier applied to three different signal compositions across different time segments.
Table 3. Performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) for SVM classifier applied to three different signal compositions across different time segments.
Signals Metrics Segment for Analysis
Whole 1.5~3 s 0.5~2.5 s 0.5~3 s 0~1.8 s
rawacc 95.80 94.31 94.72 95.05 82.52
imf 1 + imf 2 + Res 95.80 94.92 94.85 95.46 86.59
raw + Res 96.21 95.33 95.12 95.60 87.53
imf 1 + imf 3 + Res 96.34 94.99 95.05 95.66 87.47
raw sen 90.78 88.42 89.60 89.83 52.01
imf 1 + imf 2 + Res 91.49 89.60 89.83 90.31 65.01
raw + Res 91.02 90.07 89.83 90.31 68.79
imf 1 + imf 3 + Res 92.91 90.07 89.95 90.38 65.48
raw spe 97.82 96.68 96.77 96.96 94.78
imf 1 + imf 2 + Res 97.72 96.87 97.53 97.53 95.06
raw + Res 98.29 97.15 97.72 98.10 95.06
imf 1 + imf 3 + Res 98.29 97.15 97.91 98.20 96.49
raw pre 94.35 91.44 91.77 92.27 80.00
imf 1 + imf 2 + Res 94.72 92.03 93.53 94.06 84.19
raw + Res 95.53 92.54 93.98 94.97 84.84
imf 1 + imf 3 + Res 95.56 92.67 94.37 95.21 88.14
raw F 1 92.53 89.90 90.67 91.03 63.04
imf 1 + imf 2 + Res 93.08 90.80 91.64 92.15 73.37
raw + Res 93.32 91.28 91.86 92.58 75.98
imf 1 + imf 3 + Res 94.22 91.35 92.11 92.73 75.14
Legends: (1) “raw” represents the raw triaxial accelerometer signals collected at the waist position. (2) “ imf i ” (i is 1, 2, or 3) denotes the i-th decomposed IMF, respectively, from the triaxial accelerometer signals using the proposed DAGAF decomposition algorithm, whereas “Res” is the residual component after three IMFs have been decomposed. (3) “acc”, “sen”, “spe”, “pre” and “ F 1 ” are the abbreviations for accuracy, sensitivity, specificity, precision, and F 1 -score, respectively. (4) “Whole” refers to 768 data points captured from each triaxial accelerometer signal (sampled at 238 Hz), with an impact-defined window consisting of 476 points (2 s) before the impact (backward) and 292 points (1.23 s) after the impact (forward). The impact event occurs precisely at point 476. The indicated segments can be converted to time indices by multiplying by the sampling frequency and applying the ceiling operation. (5) All of the values are shown in percentage (%).
Table 4. Comparison of performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) with other studies.
Table 4. Comparison of performance metrics (accuracy, sensitivity, specificity, precision, and F 1 -score) with other studies.
Method Signals Classifier Metrics
acc sen spe pre F 1
Saleh et al. [15]raw SAE 93.46 N/A N/A N/A N/A
Šeketa et al. [22]raw SVM N/A N/A N/A N/A 89.50
Ramón et al. [23]raw OC- k NN N/A 92.90 90.62 N/A N/A
Silva et al. [21]raw SVM N/A 88.57 95.54 N/A N/A
Liu et al. [24]enhanced triaxial acceleration SVM 95.60 90.31 97.72 94.09 92.16
Proposed imf 1 + imf 2 + Res SVM 95.80 91.49 97.72 94.72 93.08
raw + Res SVM 96.21 91.02 98.29 95.53 93.22
imf 1 + imf 3 + Res SVM 96.34 92.91 98.29 95.56 94.22
Legends: (1) “raw” represents the triaxial accelerometer signals collected at the waist position. (2) “enhanced triaxial acceleration” means the triaxial accelerometer signals are preprocessed by a dedicated signal enhancement model based on a deep learning approach [24]. (3) “SAE”, “SVM”, and “OC- k NN” are the abbreviations for stacked autoencoder, support vector machine, and one-class k -nearest neighbors, respectively. (4) “ imf i ” (i is 1, 2, or 3) denotes the i-th decomposed IMF, respectively, from the triaxial accelerometer signals using the proposed DAGAF decomposition algorithm, whereas “Res” is the residual component after three IMFs have been decomposed. (5) N/A means not applicable. (6) “acc”, “sen”, “spe”, “pre”, and “ F 1 ” are the abbreviations for accuracy, sensitivity, specificity, precision, and F 1 -score, respectively. (7) All of the values are shown in percentage (%).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lin, Y.-D.; Lu, C.-J.; Sun, M.-H.; Hung, J.-H. Fall Detection Based on Data-Adaptive Gaussian Average Filtering Decomposition and Machine Learning. Information 2024, 15, 606. https://doi.org/10.3390/info15100606

AMA Style

Lin Y-D, Lu C-J, Sun M-H, Hung J-H. Fall Detection Based on Data-Adaptive Gaussian Average Filtering Decomposition and Machine Learning. Information. 2024; 15(10):606. https://doi.org/10.3390/info15100606

Chicago/Turabian Style

Lin, Yue-Der, Chi-Jen Lu, Ming-Hsuan Sun, and Ju-Hsuan Hung. 2024. "Fall Detection Based on Data-Adaptive Gaussian Average Filtering Decomposition and Machine Learning" Information 15, no. 10: 606. https://doi.org/10.3390/info15100606

APA Style

Lin, Y. -D., Lu, C. -J., Sun, M. -H., & Hung, J. -H. (2024). Fall Detection Based on Data-Adaptive Gaussian Average Filtering Decomposition and Machine Learning. Information, 15(10), 606. https://doi.org/10.3390/info15100606

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