Next Article in Journal
Multi-Lane Differential Variable Speed Limit Control via Deep Neural Networks Optimized by an Adaptive Evolutionary Strategy
Next Article in Special Issue
Prediction of Preterm Labor from the Electrohysterogram Signals Based on Different Gestational Weeks
Previous Article in Journal
Torque Measurement and Control for Electric-Assisted Bike Considering Different External Load Conditions
Previous Article in Special Issue
Estimation of Respiratory Rate during Biking with a Single Sensor Functional Near-Infrared Spectroscopy (fNIRS) System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Diabetes Mellitus Progression Using Supervised Machine Learning

1
Department of Mechanical Engineering, University of Nevada, Las Vegas, NV 89154, USA
2
Department of Mechanical Engineering, University of Washington, Seattle, WA 98195, USA
3
Department of Internal Medicine, University of Nevada, Las Vegas, NV 89154, USA
4
Department of Kinesiology and Nutrition Sciences, University of Nevada, Las Vegas, NV 89154, USA
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(10), 4658; https://doi.org/10.3390/s23104658
Submission received: 1 April 2023 / Revised: 8 May 2023 / Accepted: 9 May 2023 / Published: 11 May 2023
(This article belongs to the Special Issue Machine Learning Methods for Biomedical Data Analysis)

Abstract

:
Diabetic peripheral neuropathy (DN) is a serious complication of diabetes mellitus (DM) that can lead to foot ulceration and eventual amputation if not treated properly. Therefore, detecting DN early is important. This study presents an approach for diagnosing various stages of the progression of DM in lower extremities using machine learning to classify individuals with prediabetes (PD; n = 19), diabetes without (D; n = 62), and diabetes with peripheral neuropathy (DN; n = 29) based on dynamic pressure distribution collected using pressure-measuring insoles. Dynamic plantar pressure measurements were recorded bilaterally (60 Hz) for several steps during the support phase of walking while participants walked at self-selected speeds over a straight path. Pressure data were grouped and divided into three plantar regions: rearfoot, midfoot, and forefoot. For each region, peak plantar pressure, peak pressure gradient, and pressure–time integral were calculated. A variety of supervised machine learning algorithms were used to assess the performance of models trained using different combinations of pressure and non-pressure features to predict diagnoses. The effects of choosing various subsets of these features on the model’s accuracy were also considered. The best performing models produced accuracies between 94–100%, showing the proposed approach can be used to augment current diagnostic methods.

1. Introduction

Diabetes mellitus (DM) is a condition in which an individual is unable to either produce insulin, use naturally produced insulin correctly, or produce enough insulin to break down blood glucose [1]. The resulting high blood sugar levels can cause long-term organ damage. The World Health Organization estimates that diabetes is the seventh leading cause of death globally [2]. According to the American Diabetes Association, 11.3% of the U.S. population, or 37.3 million people, had DM in 2019 [3].
Health complications associated with diabetes include vision loss, heart disease, stroke, kidney failure, and amputation of lower extremities [4]. Those affected may also experience a higher mortality rate due to these health challenges [2]. Some of these complications can be related to peripheral neuropathy, which is the degradation of nerves in the extremities, and peripheral vascular disease, which causes reduced blood flow to the limbs [5]. The combined internal blood vessel damage coupled with the loss of sensation results in an increase in mechanical loads on the affected plantar tissue, leading to the development of plantar ulceration. These ulcerations, which are difficult to diagnose early, can lead to amputations if left untreated [5]. It was shown that DM was associated with more than half of the lower-limb amputations in the United States in 1991 [6]. It was also estimated that approximately 70% of all amputations were for people with DM [7] and about 60% of below-knee amputations were associated with DM [8]. While great strides were achieved in helping amputees cope with their losses, addressing the causes of amputation is highly desirable since approximately two million people in the United States live with limb loss. This number is expected to double by 2050 [9].
It has been reported that only 2% of the patients who had DM-related lower extremity amputations had regular rigorous foot examinations, while 29% had no examinations at all [10]. Therefore, early detection through prediction can provide information that will guide interventions to reduce the occurrence of foot ulcers and associated morbidity, mortality, and healthcare-related costs. The most widely used method for diagnosing neuropathy in a clinical setting involves testing the sensation of plantar foot tissue with a 10 g monofilament and tuning fork [11]. However, these tests may not be effective in identifying the early stages of neuropathy. Additionally, they are influenced by the subjectivity of the patient or clinician administering them [11,12,13]. More reliable techniques for diagnosing neuropathy, such as nerve conduction studies, are time-intensive and costly [13,14]. Recently, several diagnostic tests, including application of adhesive patches to measure skin hydration [15] and the use of infrared cameras to measure the foot’s thermal response to cold stimulus, have been proposed [16]. However, these techniques may also be both time-intensive and cost-prohibitive. As the number of patients with DM is still on the rise, it is imperative to predict the progression of DM in lower extremities.
A reliable, readily available, and cost-effective method to diagnose neuropathy can improve the accuracy and timeliness of detection [13]. Such a method may drive the development of effective interventions that can be applied to prevent ulcers and amputations. Early diagnosis of diabetic peripheral neuropathy will reduce the burden of diseases associated with it to the patient and the healthcare system. It was shown that measuring plantar pressures is an essential step in determining the response of tissues to stress during ambulation because it is considered to be a “surrogate measure of trauma, and an important contributing factor to skin breakdown” in patients with DM [17]. Several researchers have conducted clinical investigations to determine the variations in plantar pressures and stresses in patients with DM, since experimental data can provide insight on where high plantar pressures are applied and their potential effect on gait and development of foot ulcerations [18,19,20,21,22]. For example, the plantar tissue stress levels in the feet of patients with DM were compared for those with and without a history of foot ulcers. Results identified significantly decreased activity and higher stress levels in the DM patients with neuropathy compared to both DM patients without neuropathy and a control group [19]. It was suggested that DM patients with foot ulcers are at higher risk of plantar tissue injury even at relatively low levels of cumulative tissue stress [20]. Investigators have used wavelet transforms and neural-network classifiers with static plantar pressure measurements to classify individuals with diabetic neuropathy (DN), DM, or as non-diabetics [21,22]. Additionally, supervised machine learning algorithms and dynamic plantar pressure parameters have been used to successfully differentiate individuals with DN and DM from healthy individuals [23]. The researchers implemented binary support vector machine (SVM) classifiers to classify participants as individuals with diabetes with and without DN or as healthy controls. However, the literature shows that individuals with pre-diabetes (PD) have altered plantar pressure measures similar to those of individuals with DM when compared to healthy individuals, which can lead to individuals with PD being misclassified as having DM, or vice versa [24].
The purpose of this study was to determine the feasibility of using supervised machine learning algorithms to accurately classify persons who are diagnosed with pre-diabetes, diabetes, or diabetes with peripheral neuropathy. It was hypothesized that by combining several plantar pressure features along with patient-specific features, classification efficacy suitable for clinical use can be achieved.

2. Materials and Methods

2.1. Instrumentation

Medilogic® pressure-measuring insoles (Schönefeld, Germany) were used in this study. Six European insole sizes were used: 35–36, 37–38, 39–40, 41–42, 43–44, and 45–46. These insoles have a grid of sensors, ranging between 93 and 162 sensors depending on the size of the insole (Figure 1). The 0.75 × 1.5 cm rectangular sensors measure the change in electrical resistance, which is proportional to the applied normal pressure. The output of the sensors is a 0–255 digital scale, which is converted to 0–64 N/cm2 based on linear interpolation [25]. The sampling frequency of the insoles was set at 60 Hz.

2.2. Methods

Data were collected at either the UNLV Diabetes and Endocrinology Clinic or at the UNLV Sports and Injury Research Center. The research protocol was approved by the Institutional Review Board at the University of Nevada, Las Vegas (IRB# 777036). Based on the size of related studies [21,23], a convenience sample of 110 individuals was recruited into three groups:
  • Pre-DM (PD),
  • DM without peripheral neuropathy (D), and
  • DM with peripheral neuropathy (DN).
Participants were grouped into these three groups based on the diagnosis by their treating physicians.
Participants were briefed about the testing protocol and asked to provide informed consent. Informed consent was obtained from all subjects involved in the study. Individuals who had pre-DM and DM with or without neuropathy were invited to participate. Exclusion criteria for the study were inability to walk without assistance or support, unhealed foot ulcerations, or pregnancy beyond second trimester. The demographic data of the participants are listed in Table 1.
Participants were fitted with the Medilogic® insoles that best matched their foot size held by a pair of thin socks. The participants performed the following sequence of tasks for the walking trials: (1) sitting in a chair with both feet raised off the ground for 5 s; (2) standing up and staying still for 10 s; (3) sitting in the chair again with feet raised off the ground for 5 s; and (4) standing up and walking along a walkway (6–10 m) at their self-selected pace. The data were recorded when the participant reached a steady walking speed [26], as shown in Figure 2. Each participant was directed to perform five successful walking trials during the experiment.

2.3. Preprocessing

The sensor output data recorded during the walking trials were saved and processed with custom Matlab® (Mathworks, Natick, MA, USA) algorithms. Support phase pressure data were extracted and processed to obtain plantar pressure–time histories for each step in each trial. The data were filtered using a 5 Hz low-pass filter [17,27]. To address the anatomical structure and function of the foot relative to the main phases of the support phase of walking, the foot was divided into three regions by dividing the insole into three sections:
  • Rearfoot (RF),
  • Midfoot (MF), and
  • Forefoot (FF).
Sensors within each of these three regions of an insole size were grouped as shown in Figure 3, and the number of sensors in each region are shown in Table 2.
For each foot (left/right) and each region i, the peak plantar pressure (PPP) was recorded at each time instant. PPP was defined as:
P P P l , i = m a x p j , t , k
where l is limb (left (1)/right (2)), i is the region number, t is a time instant, and p j , t , k is pressure recorded by sensor j at time instant t of step k. Rearfoot, midfoot, and forefoot were labeled as 1, 2, and 3, respectively.
Similarly, the second set of extracted features was the peak pressure gradient (PPG) [28], which was computed at the location of the peak pressure sensor using the following equation:
P P G l , i = m a x p j , t , k p n , t , k Δ n , t , k
where, p n , t , k is pressure recorded at each of the eight sensors neighboring the peak pressure sensor p j , t , k of region i at time t of step k of the trial, and Δ n , t , k is the length of the vector connecting the center of the peak pressure sensor to each of its eight neighbors. These distances are functions of the size of the insole sensors: 0.75 × 1.5 cm. The magnitude of PPG was determined to be the maximum spatial change at the location of the peak pressure sensor.
The third set of extracted features was the pressure–time integral (PTI) for each region, which was computed by summing the pressures of all of the sensors in each region over each time instant for each participant using the following equation:
P T I l , i = t = 0 T j 1 s i p j , t , k
where t is the time instant recorded within a step, j is sensor number, and s i is the number of active sensors in region i.
From these equations, features corresponding only to the left foot were used as features or,
F 1 , i F 2 , i F 3 , i = P P P 1 , 1 : 3 P P G 1 , 1 : 3 P T I 1 , 1 : 3 i = 1 ,   2 ,   3
Additional features were introduced to assess the left/right asymmetry of the pressure features, which may be indicative of diabetic peripheral neuropathy. These features were defined as:
F 4 , i F 5 , i F 6 , i = P P P 1 , i P P P 2 , i P P P 1 , i , P P G 1 , i P P G 2 , i P P G 1 , i , P T I 1 , i P T I 2 , i P T I 1 , i i = 1 ,   2 ,   3
Therefore, a total of 18 pressure features were extracted: the PPP, PPG, and PTI for each of the three regions for the left foot, along with the corresponding asymmetry indices computed using Equation 5. These features were combined with the four non-pressure features: age, mass, height, and HbA1c (Table 1) to train, validate, and test the models using various supervised machine learning algorithms. These features were labeled as:
F 7 , 1 F 7 , 2 F 7 , 3 F 7 , 4 = A g e   M a s s   H e i g h t   H b A 1 c
Some participants completed fewer than five valid trials. For these participants, the total number of completed trials was used in the analysis. Additionally, applying the trial exclusion criteria discussed above resulted in a dataset consisting of 1621 observations or steps.

2.4. Feature Selection and Creation of Feature Subsets

A total of 22 features were included in the research: 18 pressure-related and four non-pressure features. The 18 pressure features included the values of PPP, PPG, and PTI for the three regions of the left foot, as well as their corresponding asymmetry indices. The four non-pressure features included age, mass, height, and HbA1c.
Multiple feature subsets were created using a combination of all or some of the 18 pressure and the four non-pressure features. The features that have the most impact on classifying the observations, and thereby are most useful for practical applications, can be determined by analyzing the performance of the models trained using these datasets of feature subsets.
Instead of studying the dataset with all features, it was decided to study various cases with various subsets to understand the importance of the different features in classification. Therefore, a total of 14 aggregated datasets were explored with various combinations of pressure and non-pressure features. In all of these datasets, PD, D, and DN were labeled as 0, 1, and 2 classes, respectively. The datasets were created by determining all possible combinations of the different pressure features, with the non-pressure features being included in each subset. HbA1c was of special interest, as some theorize it may not be a useful feature since some participants at an advanced stage of DM may have moderate HbA1c values if they have their blood sugar level under control. Consequently, HbA1c was included in the odd-numbered datasets only. Table 3 lists the features included in each dataset for training the models.

2.5. Classification

Each of the 14 datasets were used to train the 31 algorithms that are available within the Classification Learner App in MATLAB version R2022a (Mathworks, Natick, MA, USA), (Appendix A). The first 20% of the data were set aside for testing. The remaining 80% of the data were trained with 10 cross-validation folds, split at the individual level such that all observations related to a single individual’s trial were kept together. Models were also trained and tested using only the pressure features in the 14 datasets, as well as two datasets, 15 and 16, that tested only the non-pressure features (one with and one without HbA1c). The Classification Learner App allowed the analysis of 31 different machine learning algorithms. These included various ensemble algorithms that had been documented as being more effective than traditional machine learning models [29].
Precision, recall, the F1 score, and the false negative rate for DN were reported. Precision is the percentage of positively classified predictions that are true positive; recall is the percentage of correctly classified positive predictions out of the total true positive predictions; and the F1 score is the harmonic mean of these two metrics. A high F1 score indicates that the model is successful at predicting the positive cases, as well as a high certainty that the positive predictions are true. All performance metrics were computed based on formulas provided in the literature [30].
Since the primary purpose of testing the different datasets was to determine which features had a greater impact on classification, a principal component analysis (PCA) was carried out as a method of dimensionality reduction. The analysis was set to determine the principal components that explained at least 99% of the variability in the datasets. The analysis was performed on each dataset after the data were first standardized, and then a covariance matrix of the features was computed to calculate the eigenvectors and eigenvalues that were used to determine the principal components [31]. PCA is a reliable method of data reduction that can remove insignificant data, which can also prevent the model from overfitting the data. If PCA used with the feature subsets is able to produce high accuracy, it may provide a more generalizable and efficient model.

3. Results

The models trained using all 22 features (Dataset 1) showed that several algorithms yielded similar performance results for the test data (Table 4). The best performing algorithm along with the corresponding precision, recall, F1 score, and percentage of false negatives for DN for each dataset are listed in Table 5. The best performing algorithm produced the highest test classification accuracy with the test data. The best performing algorithms for the models trained using only the pressure features in each dataset are listed in Table 6. When only considering the pressure features, the dataset pairs became identical. For clarity, only the odd datasets are listed. Table 7 lists the best performing algorithm, precision, recall, F1 score, and false negative rate for DN recorded for the datasets after PCA was applied. Both datasets 15 and 16 had precision, recall, and F1 scores of 100%.

4. Discussion

The purpose of this study was to determine the feasibility of using supervised machine learning algorithms to accurately classify persons who are diagnosed with pre-diabetes, diabetes, or diabetes with peripheral neuropathy. The classification accuracies, as well as the false negative rates for DN in the test dataset, were used to determine which algorithms and features subsets performed best. The false negative rates for DN participants are especially important because of the difficulty in diagnosing DN using current techniques.
Table 5 indicates that all datasets yielded reasonably accurate classifications. However, since false negative rates are important for this type of analysis, it is reasonable to consider dataset 1 and datasets 11 through 14; all had zero percent false negative rates. Of these five datasets, datasets 13 and 14 had the highest F1 score, of 100%. The results may indicate that a subset of the tested features can be more successful in classification than using all available features. A recent study that used biomechanical data to train and test machine learning algorithms for DN diagnosis also found that using a subset of features, rather than the entire dataset, yielded greater accuracies [32]. Thus, the identification and inclusion of significant features, rather than all available features, could produce more effective and less computationally costly algorithms.
The only pressure features in datasets 13 and 14 were for PTI. This indicates that PTI may play a significant role in DN classification. A previous study indicated that rearfoot values of peak pressure are significant in classification [22]. PPP has also been reported to be a significant indicator for DN [17,22]. This was corroborated in this study by the high performance of datasets 9 and 10, which relied solely on PPP pressure features, and datasets 5 and 6, which used PPP and PTI features. For these specific datasets, ensemble bagged trees and ensemble subspace KNN classifiers produced the best classifications of the data. It is also of note that subspace KNN was the second-best performing algorithm for five of the seven datasets where it was not the best performing, excepting only datasets 3 and 7. Table 5 also indicates that the exclusion of HbA1c did not yield any significant difference in results.
Testing only the non-pressure features, with and without HbA1c, resulted in 100% in all three measures. This may be because during data collection, while the pressure values would change for each trial, the non-pressure features for each participant would remain the same, since the measures were not variable during a single data collection session. As such, the current dataset did not produce any significant trends or results when looking only at the non-pressure features. No trends in the non-pressure features when plotted across the three classes were found, either. Since the non-pressure features did not have any discernable pattern in delineating the three classes, it could not be concluded that non-pressure features are all that is needed to classify individuals with PD, D, or DN, and further exploration into the impact of these features on classification was unfruitful. The non-pressure features alone were thus concluded to not have the ability to accurately classify the data. Analyzing solely the pressure features allowed the features that played a more significant role in accurate classification to become clearer.
The results of analyzing only the pressure features (Table 6) indicate that datasets 1, 5, 7, and 13 had test accuracies of greater than 80% and false negative rates of less than 3%. It should be noted that dataset 1 included all pressure features; dataset 5 had PPP and PTI features; dataset 7 combined PPG and PTI features; and dataset 13 only contained PTI features. The recurrence of PTI indicated its significance in classifying the progression of diabetes. When considering only pressure features, the results indicated that including more features does produce a higher test accuracy, but only when PTI features are included. For example, dataset 3, which included PPP and PPG features, achieved a test accuracy of 65% and false negative rate of 39%. Similarly, dataset 9 included only PPP features, and dataset 11 only included PPG features. All of these datasets produced test accuracies of less than 70% and false negative rates of greater than 50%. Overall, the inclusion of PTI, even as the sole pressure feature, greatly enhanced all three measures. The results from Table 6 indicate that a variety of pressure features used alone can still provide adequately accurate classifications, as datasets 1 and 5 had precision, recall, and F1 scores of higher than 90%. By comparing the results from Table 5 and Table 6, however, it can be seen that if using only a limited number of pressure features, having non-pressure features will greatly enhance the feasibility of this method, exemplified best by the drastic performance difference in dataset 11.
A study that used center of pressure data to train and test a deep clustering model found that individuals with DN tend to apply less force while walking and have longer stance times during gait [33]. As such, individuals who have DN might have values for PPG or PPP similar to those of the other classes, but their PTI values are more distinct, since the feature combines pressure, foot size (due to the increased number of sensors for larger insoles), and time. Additionally, since PTI represents the combined effect of these three factors, it can avoid the variability due to gait variation better than PPP or PPG. Given the high presence of PTI in the best performing algorithms, further study into the relationship between foot size, walking speed, and PTI is desirable.
Table 7 shows that PCA achieved comparable results for the best datasets of Table 5: datasets 1 and 11 through 14. The results that were achieved using ensemble classifier with subspace KNN model produced the highest results in conjunction with PCA. Apart from datasets 8, 9, and 12–14, PCA also led to higher false negative rates. While dimensionality reduction was explored to see if the model performance could be enhanced, the results indicate that PCA was unsuccessful in this.
While ensemble classifier with bagged trees was a high performing algorithm throughout the results, various forms of KNN also performed well for half of the datasets (Table 5). KNN performing well has been documented before in a previous study that also analyzed multiple algorithms in classifying DN using EMG and gait data [34].

5. Conclusions

This paper aimed to determine the feasibility of using machine learning algorithms to accurately classify individuals with pre-diabetes, diabetes without neuropathy, and diabetes with neuropathy. Using plantar pressure features (PPP, PPG, PTI), along with the non-pressure features (age, mass, height, and HbA1c), 14 datasets with different feature combinations were trained to determine the combination of features that provided the highest overall accuracy with the lowest false negative rates for diabetic peripheral neuropathy.
It was determined that using multiple pressure features yielded high classification accuracies. Though non-pressure features provided no visible trends that could be used to differentiate between the three classes, the addition of non-pressure features significantly contributed toward improving the model’s classification ability. The results show that PTI features played a stronger role in classification of the three subject groups, with the cases that included PTI yielding test accuracies of greater than 90%. On the other hand, the best case without PTI had a test accuracy of about 70%. Additionally, cases including PTI features had false negative rates of less than 3%. These observations are reasonable, since PTI combines pressure, foot size, and time. Since it represents the combined effect of these three factors over an area, it avoids the potential anomalies in peak pressure and pressure gradients that may be due to gait variations. It was also observed that the lowest false negative rate from a dataset that did not include PTI was 39.5%. The results also show that including HbA1c in the features consistently improved the results slightly.
The cubic SVM model had a test classification accuracy of 98.5% with no false negative classifications when using all of the available features. However, when classifying with a subset of features, ensemble classifiers with bagged trees or subspace KNN models performed the best. Various forms of KNN also performed well in all the analyses, making KNN an attractive algorithm for further analysis.
The effect of using PCA to reduce the dimensionality of the data was minimal. Overall, PCA was successful, though it led to higher false negative rates in nine of the 14 datasets.
This study was limited by the relatively small sample size, which led to an imbalanced dataset. A larger dataset with a more even split between classes would create a more generalizable model. Additionally, participants were asked to wear socks in this study to eliminate the uncertainty that could be caused by the type of shoe, but this does not represent daily life. A follow-up study that considers the effect of shoes on plantar pressure warrants investigation. Another limitation is that only HbA1c was included as a measure of DM severity. As mentioned before, HbA1c is a marker that can be controlled with a medication regime and lifestyle choices and is not always indicative of an individual’s diabetic progression. A future study including different means to assess DM severity, such as the length of diagnosis, would be beneficial. Similarly, arch type, foot size, and walking speed should be considered in future studies.
Overall, the research indicates that supervised machine learning algorithms are well-suited for classifying the progression of DM. It would be of interest to repeat the study using a larger number of participants.

Author Contributions

Conceptualization, A.S.C., M.B.T. and J.S.D.; data curation: A.S.C., M.S.V. and M.B.T.; methodology, M.S.V., M.B.T. and J.S.D.; software, A.S.C., M.B.T. and M.S.V.; validation, A.S.C., M.B.T. and J.S.D.; formal analysis, A.S.C.; investigation, M.S.V.; resources, M.B.T., J.S.D. and K.I.; writing—original draft preparation, A.S.C.; writing—review and editing, A.S.C., M.B.T., J.S.D. and M.S.V.; visualization, A.S.C., M.B.T. and J.S.D.; supervision, M.B.T. and J.S.D.; project administration, A.S.C., M.S.V., M.B.T. and J.S.D. 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 conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board (or Ethics Committee) of the University of Nevada, Las Vegas (IRB#777036, originally approved 2 August 2017, recently amended 5 February 2023).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent was obtained from the patient(s) to publish this paper.

Data Availability Statement

The data presented in this study will be available in the future upon request to the authors.

Acknowledgments

Thank you to all of the volunteers who participated in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Algorithms included in MATLAB Classification Learner App and their hyperparameters.
Table A1. Algorithms included in MATLAB Classification Learner App and their hyperparameters.
AlgorithmHyperparameters
Fine TreeMaximum number of splits: 100
Split Criterion: Gini’s diversity index
Surrogate decision splits: Off
Medium TreeMaximum number of splits: 20
Split Criterion: Gini’s diversity index
Surrogate decision splits: Off
Coarse TreeMaximum number of splits: 4
Split Criterion: Gini’s diversity index
Surrogate decision splits: Off
Linear DiscriminantCovariance structure: Full
Quadratic DiscriminantCovariance structure: Full
Gaussian Naïve BayesDistribution name for numeric predictors: Gaussian
Distribution name for categorical predictors: Not Applicable
Kernel Naïve BayesDistribution name for numeric predictors: Kernel
Distribution name for categorical predictors: Not Applicable
Kernel type: Gaussian
Support: Unbounded
Linear SVMKernel function: Linear
Kernel scale: Automatic
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Quadratic SVMKernel function: Quadratic
Kernel scale: Automatic
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Cubic SVMKernel function: Cubic
Kernel scale: Automatic
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Fine Gaussian SVMKernel function: Gaussian
Kernel scale: 1.2
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Medium Gaussian SVMKernel function: Gaussian
Kernel scale: 4.8
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Coarse Gaussian SVMKernel function: Gaussian
Kernel scale: 19
Box constraint level: 1
Multiclass method: One-vs-One
Standardize data: true
Fine KNNNumber of neighbors: 1
Distance metric: Euclidean
Distance weight: Equal
Standardize data: true
Medium KNNNumber of neighbors: 10
Distance metric: Euclidean
Distance weight: Equal
Standardize data: true
Coarse KNNNumber of neighbors: 100
Distance metric: Euclidean
Distance weight: Equal
Standardize data: true
Cosine KNNNumber of neighbors: 10
Distance metric: Cosine
Distance weight: Equal
Standardize data: true
Cubic KNNNumber of neighbors: 10
Distance metric: Minkowski (cubic)
Distance weight: Equal
Standardize data: true
Weighted KNNNumber of neighbors: 10
Distance metric: Euclidean
Distance weight: Squared inverse
Standardize data: true
Boosted TreesEnsemble method: AdaBoost
Learner type: Decision tree
Maximum number of splits: 20
Number of learners: 30
Learning rate: 0.1
Number of predictors to sample: Select All
Bagged TreesEnsemble method: Bag
Learner type: Decision tree
Maximum number of splits: 1296
Number of learners: 30
Number of predictors to sample: Select All
Subspace DiscriminantEnsemble method: Subspace
Learner type: Discriminant
Number of learners: 30
Subspace dimension: 12
Subspace KNNEnsemble method: Subspace
Learner type: Nearest neighbors
Number of learners: 30
Subspace dimension: 12
RUSBoosted TreesEnsemble method: RUSBoost
Learner type: Decision tree
Maximum number of splits: 20
Number of learners: 30
Learning rate: 0.1
Number of predictors to sample: Select All
Narrow Neural NetworkNumber of fully connected layers: 1
First layer size: 10
Activation: ReLU
Iteration limit: 1000
Regularization strength (Lambda): 0
Standardize data: Yes
Medium Neural NetworkNumber of fully connected layers: 1
First layer size: 25
Activation: ReLU
Iteration limit: 1000
Regularization strength (Lambda): 0
Standardize data: Yes
Wide Neural NetworkNumber of fully connected layers: 1
First layer size: 100
Activation: ReLU
Iteration limit: 1000
Regularization strength (Lambda): 0
Standardize data: Yes
Bilayered Neural NetworkNumber of fully connected layers: 2
First layer size: 10
Activation: ReLU
Iteration limit: 1000
Regularization strength (Lambda): 0
Standardize data: Yes
Trilayered Neural NetworkNumber of fully connected layers: 3
First layer size: 10
Activation: ReLU
Iteration limit: 1000
Regularization strength (Lambda): 0
Standardize data: Yes
SVM KernelLearner: SVM
Number of expansion dimensions: Auto
Regularization strength (Lambda): Auto
Kernel scale: Auto
Multiclass method: One-vs-One
Iteration limit: 1000
Logistic Regression KernelLearner: Logistic Regression
Number of expansion dimensions: Auto
Regularization strength (Lambda): Auto
Kernel scale: Auto
Multiclass method: One-vs-One
Iteration limit: 1000

References

  1. American Diabetes Association. Common Terms. Available online: https://www.diabetes.org/resources/students/common-terms (accessed on 23 August 2022).
  2. Word Health Organization. Diabetes. Available online: https://www.who.int/news-room/fact-sheets/detail/diabetes (accessed on 20 April 2023).
  3. Statistics About Diabetes: American Diabetes Association. Available online: http://www.diabetes.org/diabetes-basics/statistics/ (accessed on 20 April 2023).
  4. Centers for Disease Control and Prevention. Prevent Diabetes Complications. Available online: https://www.cdc.gov/diabetes/managing/problems.html#:~:text=Common%20diabetes%20health%20complications%20include,how%20to%20improve%20overall%20health (accessed on 23 August 2022).
  5. Cavanagh, P.R.; Simoneau, G.G.; Ulbrecht, J.S. Ulceration, unsteadiness, and uncertainty: The biomechanical consequences of diabetes mellitus. J. Biomech. 1993, 26 (Suppl. 1), 23–40. Available online: http://www.ncbi.nlm.nih.gov/pubmed/8505350 (accessed on 23 August 2022).
  6. Pecoraro, R.; Reiber, G.; Burgess, E. Pathways to Diabetic Basis for Prevention. Diabetes Care. 1990, 13, 513–521. [Google Scholar] [CrossRef]
  7. Boulton, A.; Vileikyte, L.; Ragnarson-Tennvall, G.; Apelqvist, J. The global burden of diabetic foot disease. Lancet 2005, 366, 1719–1724. [Google Scholar] [CrossRef] [PubMed]
  8. Zayed, M.; Bech, F.; Hernandez-Boussard, T. National review of factors influencing disparities and types of major lower extremity amputations. Ann. Vasc. Surg. 2014, 28, 1157–1165. [Google Scholar] [CrossRef]
  9. Sheehan, T.; Gondo, G. Impact of Limb Loss in the United States. Phys. Med. Rehabil. Clin. N. Am. 2014, 25, 9–28. Available online: http://www.sciencedirect.com/science/article/pii/S1047965113000740 (accessed on 22 October 2014). [CrossRef] [PubMed]
  10. Cook, E.; Cook, J.; Labre, M.; Givens, H.; Di Resta, J. The Amputation Prevention Initiative. J. Am. Podiatr. Med. Assoc. 2014, 104, 1–10. Available online: http://www.japmaonline.org/doi/abs/10.7547/0003-0538-104.1.1 (accessed on 22 October 2014). [CrossRef] [PubMed]
  11. Dros, J.; Wewerinke, A.; Bindels, P.J.; Van Weert, H.C. Accuracy of monofilament testing to diagnose peripheral neuropathy: A systematic review. Ann. Fam. Med. 2009, 7, 555–558. [Google Scholar] [CrossRef]
  12. Wang, F.; Zhang, J.; Yu, J.; Liu, S.; Zhang, R.; Ma, X.; Yang, Y.; Wang, P. Diagnostic Accuracy of Monofilament Tests for Detecting Diabetic Peripheral Neuropathy: A Systematic Review and Meta-Analysis. J. Diabetes Res. 2017, 2017, 8787261. [Google Scholar] [CrossRef]
  13. Selvarajah, D.; Kar, D.; Khunti, K.; Davies, M.J.; Scott, A.R.; Walker, J.; Tesfaye, S. Diabetic peripheral neuropathy: Advances in diagnosis and strategies for screening and early intervention. Lancet Diabetes Endocrinol. 2019, 7, 938–948. [Google Scholar] [CrossRef]
  14. Fuller, G. How to get the most out of nerve conduction studies and electromyography. Neurol. Pract. 2005, 76, ii41–ii46. [Google Scholar] [CrossRef]
  15. Papanas, N.; Boulton, A.J.M.; Malik, R.; Manes, C.; Schnell, O.; Spallone, V.; Tentolouris, N.; Tesfaye, S.; Valensi, P.; Ziegler, D.; et al. A simple new non-invasive sweat indicator test for the diagnosis of diabetic neuropathy. Diabet. Med. 2013, 30, 525–534. [Google Scholar] [CrossRef]
  16. Liu, C.; van der Heijden, F.; Klein, M.E.; van Baal, J.G.; Bus, S.A.; van Netten, J.J. Infrared dermal thermography on diabetic feet soles to predict ulcerations: A case study. SPIE BiOS. Int. Soc. Opt. Photonics 2013, 8572, 102–110. [Google Scholar] [CrossRef]
  17. Zou, D.; Mueller, M.J.; Lott, D.J. Effect of peak pressure and pressure gradient on subsurface shear stresses in the neuropathic foot. J. Biomech. 2007, 40, 883–890. [Google Scholar] [CrossRef] [PubMed]
  18. Cavanagh, P.R.; Hewitt, F.; Perry, J. In-shoe plantar pressure measurement: A review. Foot 1992, 2, 185–194. [Google Scholar] [CrossRef]
  19. Perry, J.; Ulbrecht, J.; Derr, J.; Cavanagh, P. The use of running shoes to reduce plantar pressures in patients who have diabetes. J. Bone Jt. Surg. 1995, 77 Pt A, 1819–1828. Available online: http://jbjs.org/article.aspx?articleid=22857 (accessed on 3 September 2014).
  20. Maluf, K.; Mueller, M. Comparison of physical activity and cumulative plantar tissue stress among subjects with and without diabetes mellitus and a history of recurrent plantar ulcers. Clin. Biomech. 2003, 18, 567–575. [Google Scholar] [CrossRef]
  21. Acharya, U.R.; Tan, P.H.; Subramaniam, T.; Tamura, T.; Chua, K.C.; Goh, S.C.E.; Lim, C.M.; Goh, S.Y.D.; Chung, K.R.C.; Law, C.; et al. Automated Identification of Diabetic Type 2 Subjects with and without Neuropathy Using Wavelet Transform on Pedobarograph. J. Med. Syst. 2008, 32, 21–29. [Google Scholar] [CrossRef]
  22. Rahman, M.A.; Aziz, Z.; Acharya, U.R.; Ha, T.P.; Kannathal, N.; Ng, E.; Law, C.; Subramaniam, T.; Shuen, W.Y.; Fang, S.C. Analysis of plantar pressure in diabetic type 2 subjects with and without neuropathy. ITBM-RBM 2006, 27, 46–55. [Google Scholar] [CrossRef]
  23. Botros, F.S.; Taher, M.F.; Elsayed, N.M.; Fahmy, A.S. Prediction of diabetic foot ulceration using spatial and temporal dynamic plantar pressure. In Proceedings of the 2016 8th Cairo International Biomedical Engineering Conference (CIBEC), Cairo, Egypt, 15–17 December 2016; Volume 2016, pp. 43–47. [Google Scholar] [CrossRef]
  24. Robinson, C.C.; Balbinot, L.F.; Silva, M.F.; Achaval, M.; Zaro, M.A. Plantar Pressure Distribution Patterns of Individuals with Prediabetes. J. Diabetes Sci. Technol. 2013, 7, 1113–1121. [Google Scholar] [CrossRef]
  25. T&T Medilogic Medizintechnik GmbH. Manual: Medilogic Pressure Measurement; T&T Medilogic Medizintechnik GmbH: Schöenfeld, Germany, 2012; pp. 48–63. [Google Scholar]
  26. De Berardinis, J. The Development of a Viscoelastic Ellipsoidal Model for Use in Measuring Plantar Tissue Material Properties during Walking. Ph.D. Thesis, Department of Mechanical Engineering, UNLV, Las Vegas, NV, USA, 2019. Available online: https://www.proquest.com/docview/2312584828?pq-origsite=gscholar&fromopenview=true (accessed on 3 September 2022).
  27. De Berardinis, J.; Dufek, J.S.; Trabia, M.B.; Lidstone, D.E. Assessing the validity of pressure-measuring insoles in quantifying gait variables. J. Rehabil. Assist. Technol. Eng. 2018, 5, 2055668317752088. [Google Scholar] [CrossRef]
  28. Mueller, M.J.; Zou, D.; Lott, D.J. “Pressure gradient” as an indicator of plantar skin injury. Diabetes Care 2005, 28, 2908–2912. [Google Scholar] [CrossRef]
  29. Khan, Y.F.; Kaushik, B.; Chowdhary, C.L.; Srivastava, G. Ensemble Model for Diagnostic Classification of Alzheimer’s Disease Based on Brain Anatomical Magnetic Resonance Imaging. Diagnostics 2022, 12, 3193. [Google Scholar] [CrossRef]
  30. Hossin, M.; Sulaiman, M.N. A Review on Evaluation Metrics for Data Classification Evaluations. Int. J. Data Min. Knowl. Manag. Process 2015, 5, 1–12. [Google Scholar] [CrossRef]
  31. BuiltIn. A Step-by-Step Explanation of Principal Component Analysis (PCA). Available online: https://builtin.com/data-science/step-step-explanation-principal-component-analysis (accessed on 23 July 2021).
  32. Guiotto, A.; Bortolami, G.; Ciniglio, A.; Spolaor, F.; Guarneri, G.; Avogaro, A.; Cibin, F.; Silvestri, F.; Sawacha, Z. Machine learning approach to diabetic foot risk classification with biomechanics data. Gait Posture 2022, 97, 30–31. [Google Scholar] [CrossRef]
  33. Ji, X.; Akiyarna, Y.; Yamada, Y.; Okamoto, S.; Hayashi, H. Development of deep clustering model to stratify occurrence risk of diabetic foot ulcers based on foot pressure patterns and clinical indices. In Proceedings of the 2020 IEEE In-ternational Joint Conference on Biometrics (IJCB), Houston, TX, USA, 28 September–1 October 2020; pp. 1–8. [Google Scholar] [CrossRef]
  34. Haque, F.; Reaz, M.B.I.; Chowdhury, M.E.H.; Ezeddin, M.; Kiranyaz, S.; Alhatou, M.; Ali, S.H.M.; A Bakar, A.A.; Srivastava, G. Machine Learning-Based Diabetic Neuropathy and Previous Foot Ulceration Patients Detection Using Electromyography and Ground Reaction Forces during Gait. Sensors 2022, 22, 3507. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) A typical Medilogic insole. (b) Corresponding sensor map.
Figure 1. (a) A typical Medilogic insole. (b) Corresponding sensor map.
Sensors 23 04658 g001
Figure 2. A participant fitted with socks conducting trial.
Figure 2. A participant fitted with socks conducting trial.
Sensors 23 04658 g002
Figure 3. Representation of the typical division of the plantar foot overlaid on a size 35 insole map.
Figure 3. Representation of the typical division of the plantar foot overlaid on a size 35 insole map.
Sensors 23 04658 g003
Table 1. Participant demographic data for the pre-DM (PD), DM without peripheral neuropathy (D), and DM with peripheral neuropathy (DN) groups (mean ± standard deviation).
Table 1. Participant demographic data for the pre-DM (PD), DM without peripheral neuropathy (D), and DM with peripheral neuropathy (DN) groups (mean ± standard deviation).
PDDDN
No. of participants196229
Sex (M/F)9/1029/3314/15
Age (Years)59.6 ± 11.458 ± 15.863.8 ± 10.2
Body Mass (kg)93.7 ± 27.690.8 ± 26.296.9 ± 25.9
Height (m)1.68 ± 0.091.66 ± 0.091.70 ± 0.10
HbA1c (%) 5.9 ± 0.47.6 ± 1.67.2 ± 1.3
Table 2. Number of sensors in each foot region, for each insole size.
Table 2. Number of sensors in each foot region, for each insole size.
European Insole SizeRegionNumber of Sensors
35–36RF29
MF29
FF35
37–38RF38
MF32
FF37
39–40RF34
MF32
FF50
41–42RF45
MF34
FF50
43–44RF47
MF38
FF66
45–46RF48
MF36
FF78
Table 3. List of datasets with the features ( F i , j ) included: Pressure features (i = 1:6) at the regions (j = 1:3): rearfoot (1), midfoot (2), and forefoot (3) regions, and non-pressure features (i = 7; j = 1:4).
Table 3. List of datasets with the features ( F i , j ) included: Pressure features (i = 1:6) at the regions (j = 1:3): rearfoot (1), midfoot (2), and forefoot (3) regions, and non-pressure features (i = 7; j = 1:4).
Training DatasetDescriptionNo. of FeaturesFeatures
1Pressure Features: PPP, PPG, PTI, Corresponding Asymmetry18 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 2 , 1 , F 2 , 2 , F 2 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3
F 4 , 1 , F 4 , 2 , F 4 , 3 , F 5 , 1 , F 5 , 2 , F 5 , 3 , F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
2Pressure Features: PPP, PPG, PTI, Corresponding Asymmetry18 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 2 , 1 , F 2 , 2 , F 2 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3
F 4 , 1 , F 4 , 2 , F 4 , 3 , F 5 , 1 , F 5 , 2 , F 5 , 3 , F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
3Pressure Features: PPP, PPG, Corresponding Asymmetry12 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 2 , 1 , F 2 , 2 , F 2 , 3 , F 4 , 1 , F 4 , 2 , F 4 , 3 F 5 , 1 , F 5 , 2 , F 5 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
4Pressure Features: PPP, PPG, Corresponding Asymmetry12 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 2 , 1 , F 2 , 2 , F 2 , 3 ,   F 4 , 1 , F 4 , 2 , F 4 , 3 F 5 , 1 , F 5 , 2 , F 5 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
5Pressure Features: PPP, PTI, Corresponding Asymmetry12 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3 , F 4 , 1 , F 4 , 2 , F 4 , 3 F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
6Pressure Features: PPP, PTI, Corresponding Asymmetry12 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3 , F 4 , 1 , F 4 , 2 , F 4 , 3 F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
7Pressure Features: PPG, PTI, Corresponding Asymmetry12 F 2 , 1 , F 2 , 2 , F 2 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3 , F 5 , 1 , F 5 , 2 , F 5 , 3 F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
8Pressure Features: PPG, PTI, Corresponding Asymmetry12 F 2 , 1 , F 2 , 2 , F 2 , 3 , F 3 , 1 , F 3 , 2 , F 3 , 3 , F 5 , 1 , F 5 , 2 , F 5 , 3 F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
9Pressure Features: PPP and Corresponding Asymmetry6 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 4 , 1 , F 4 , 2 , F 4 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
10Pressure Features: PPP and Corresponding Asymmetry6 F 1 , 1 , F 1 , 2 , F 1 , 3 , F 4 , 1 , F 4 , 2 , F 4 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
11Pressure Features: PPG and Corresponding Asymmetry6 F 2 , 1 , F 2 , 2 , F 2 , 3 F 5 , 1 , F 5 , 2 , F 5 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
12Pressure Features: PPG and Corresponding Asymmetry6 F 2 , 1 , F 2 , 2 , F 2 , 3 , F 5 , 2 , F 5 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
13Pressure Features: PTI and Corresponding Asymmetry6 F 3 , 1 , F 3 , 2 , F 3 , 3 , F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
14Pressure Features: PTI and Corresponding Asymmetry6 F 3 , 1 , F 3 , 2 , F 3 , 3 , F 6 , 1 , F 6 , 2 , F 6 , 3
Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
15Non-Pressure Features4 F 7 , 1 , F 7 , 2 , F 7 , 3 , F 7 , 4
16Non-Pressure Features3 F 7 , 1 , F 7 , 2 , F 7 , 3
Table 4. Precision, recall, F1 score, and false negative rates for the best performing algorithms, ranked by test accuracy for Dataset 1.
Table 4. Precision, recall, F1 score, and false negative rates for the best performing algorithms, ranked by test accuracy for Dataset 1.
AlgorithmPrecision (%)Recall (%)F1 ScoreFalse Negative Rate (%)
Cubic Support Vector Machine (SVM)97.998.598.20
Subspace K-Nearest Neighbors (KNN) *96.998.597.71.9
Bagged Trees *98.096.897.41.9
Wide Neural Network96.498.497.41.9
Boosted Trees97.696.497.01.9
Fine KNN96.696.796.62.8
Weighted KNN96.097.296.60.9
Fine Tree96.395.495.91.9
Medium Neural Network94.297.395.70
Trilayered Neural Network94.995.695.22.8
Quadratic SVM94.095.494.70
RUSBoosted Trees *96.693.194.81.9
* Ensemble.
Table 5. Best performing algorithm for each dataset with precision, recall, F1 score, and false negative rate.
Table 5. Best performing algorithm for each dataset with precision, recall, F1 score, and false negative rate.
DatasetBest Performing AlgorithmPrecision (%)Recall (%)F1 ScoreFalse Negative Rate (%)
1Cubic SVM97.998.598.20
2Subspace KNN *98.398.798.42.4
3Bagged Trees *98.499.398.74.7
4Bagged Trees *98.299.598.80
5Subspace KNN *99.099.299.12.3
6Subspace KNN *99.299.299.22.3
7Bagged Trees *98.598.598.53.5
8Subspace KNN *98.098.998.53.5
9Bagged Trees *1001001000
10Subspace KNN *98.699.599.02.4
11Bagged Trees *99.299.499.30
12Bagged Trees *99.899.699.70
13Subspace KNN *1001001000
14Subspace KNN *1001001000
* Ensemble.
Table 6. Best performing algorithm, precision, recall, F1 score, and false negative rate for each dataset with only pressure features.
Table 6. Best performing algorithm, precision, recall, F1 score, and false negative rate for each dataset with only pressure features.
DatasetBest Performing AlgorithmPrecision (%)Recall (%)F1 Score (%)False Negative Rate (%)
1Fine KNN91.292.391.70
3Fine KNN64.665.164.839.5
5Quadratic SVM92.793.493.10
7Medium Neural Network88.689.288.92.4
9Weighted KNN *60.964.162.552.3
11Bagged Trees *49.957.553.467.1
13Fine Tree82.988.685.72.4
* Ensemble.
Table 7. Best performing algorithm, number of components, precision, recall, F1 score, and false negative rate for each dataset with features selected through PCA.
Table 7. Best performing algorithm, number of components, precision, recall, F1 score, and false negative rate for each dataset with features selected through PCA.
DatasetAlgorithmNumber of ComponentsPrecision (%)Recall (%)F1 Score (%)False Negative Rate (%)
1Subspace KNN *993.195.594.24.7
2Subspace KNN *994.494.794.64.7
3Subspace KNN *893.093.893.45.8
4Subspace KNN *895.596.796.13.5
5Subspace KNN *696.597.196.88.1
6Subspace KNN *696.997.597.24.7
7Subspace KNN *794.895.695.25.9
8Bilayered Neural Network697.396.296.81.2
9Subspace KNN *697.998.098.01.2
10Bagged Trees *695.597.496.43.5
11Bilayered Neural Network696.197.196.68.2
12Subspace KNN *694.895.595.21.2
13Subspace KNN *51001001000
14Subspace KNN *51001001000
* Ensemble.
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

Chauhan, A.S.; Varre, M.S.; Izuora, K.; Trabia, M.B.; Dufek, J.S. Prediction of Diabetes Mellitus Progression Using Supervised Machine Learning. Sensors 2023, 23, 4658. https://doi.org/10.3390/s23104658

AMA Style

Chauhan AS, Varre MS, Izuora K, Trabia MB, Dufek JS. Prediction of Diabetes Mellitus Progression Using Supervised Machine Learning. Sensors. 2023; 23(10):4658. https://doi.org/10.3390/s23104658

Chicago/Turabian Style

Chauhan, Apoorva S., Mathew S. Varre, Kenneth Izuora, Mohamed B. Trabia, and Janet S. Dufek. 2023. "Prediction of Diabetes Mellitus Progression Using Supervised Machine Learning" Sensors 23, no. 10: 4658. https://doi.org/10.3390/s23104658

APA Style

Chauhan, A. S., Varre, M. S., Izuora, K., Trabia, M. B., & Dufek, J. S. (2023). Prediction of Diabetes Mellitus Progression Using Supervised Machine Learning. Sensors, 23(10), 4658. https://doi.org/10.3390/s23104658

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