T2DFitTailor: A tool for type 2 diabetes patients to tailor exercise plan

Data Resources

Our research utilized an extensive data collection process, focusing on the effects of exercise on Type 2 Diabetes (T2D). The search for relevant data was conducted through the PubMed database, employing a comprehensive set of keywords to ensure the inclusion of pertinent studies. The keywords used for this search were (“Type 2 Diabetes” OR “Diabetes Mellitus Type 2” OR “T2DM”) AND (“Exercise” OR “Physical Activity” OR “Aerobic Training” OR “Resistance Training” OR “Exercise Therapy” OR “Physical Fitness” OR “Physical Exercise”). The raw data are available at supplementary materials.^{(11),(12),(13)}

Definition of Responder and Non-Responder

In our analysis, individuals were categorized based on their HbA1c response to the exercise interventions. A sample was defined as a responder if there was a positive reduction in HbA1c levels following the exercise intervention. Conversely, a sample was categorized as a non-responder if no change or an increase in HbA1c levels was observed post-intervention.⁽¹⁰⁾

LASSO Regression Analysis

Least Absolute Shrinkage and Selection Operator (LASSO) regression, a variation of linear regression, employs shrinkage to promote model simplicity and sparsity by driving some regression coefficients towards zero. This approach effectively eliminates variables with null coefficients, simplifying the model selection process and mitigating overfitting. In our research, LASSO regression was instrumental in isolating pivotal features relevant to HbA1c fluctuation post-exercise intervention. We implemented LASSO regression in Python using the “sklearn.linear_model” module, with the function of “Lasso”.

Leave-One-Out (LOO)

Leave-One-Out (LOO) cross-validation, a technique where each observation in the dataset is sequentially used as the validation data while the rest serve as the training set, offers a comprehensive assessment of a regression model’s predictive accuracy, particularly beneficial in contexts of limited data size. In our study on Type 2 Diabetes (T2D), LOO was crucial for validating the precision of our regression model in predicting HbA1c reduction following exercise interventions. This method enabled a close examination of how well the predicted HbA1c reductions mirrored actual outcomes, ensuring robustness and mitigating overfitting risks. Our implementation of LOO, aligned with the project’s aim to measure the effectiveness of personalized exercise plans in managing T2D, employed the function of “LeaveOneOut” in python’s “sklearn.model_selection” module.

K-Nearest Neighbors (KNN) Algorithm

The k-Nearest Neighbors (kNN) algorithm, a fundamental machine learning method for both classification and regression tasks, operates on the principle of utilizing the nearest data points in the feature space to predict the outcome for a new sample. In our research focused on Type 2 Diabetes (T2D), we employed the kNN algorithm to predict the HbA1c reduction in patients following various exercise interventions, identifying them as either responders or non-responders based on their individual profile variation. To fine-tune our model, we determined the optimal ‘k’ value through Leave-One-Out (LOO) cross-validation on the training set, ensuring a balance between bias and variance that is critical for the predictive accuracy of the model. Our analysis found that different exercise plan required distinct ‘k’ values for optimal results, with Taiji and Cycling performing best at k=2, Qigong, Stretching, and Walking at k=5, and Rugby at k=4. This LOO process not only facilitated the identification of the best ‘k’ values but also verified the model’s effectiveness on the training set. We further validated our model against a preserved test set based on the whole training set with above k values. The finalized kNN models on T2DFitTailor were built on the combination of training and test set with optimal k values. The implementation of the kNN regression was achieved utilizing the “sklearn.neighbors module” module with the function of “KNeighborsRegressor”.

Tailored exercise plan from T2DFitTailor

In this section, we elucidate the application of the T2DFitTailor package, which is designed to develop personalized exercise plans for individuals with Type 2 Diabetes (T2D) and forecast the HbA1c outcomes of these plans. The process requires inputting seven essential parameters for T2D patients: “Age”, “Sex”, “BMI”, “Duration_T2D (years)”, “HDL (mmol/L)”, “PCS”, and “WHtR”, with the “PCS” value calculable via the PCS_Calculation function in the package (details on these inputs are available in Table 1). To demonstrate, we create eight random samples based on the range of these parameters in raw data.

Upon utilizing the TailorExercisePlan function, an analysis is conducted to assess the suitability between the T2D patient’s profile and six predefined exercise types: “Taiji”, “Qigong”, “Stretching”, “Rugby”, “Cycling”, and “Walking” (See all the output parameters in Table 1). We identify and recommend exercise plans with a predicted reduction in HbA1c levels post-exercise intervention. Upon the execution of this recommendation criterion, subsequent visualization analysis is conducted with the aid of the VisualizeTailoredExercisePlan function, showcasing the outcomes for all above samples in Supplementary Figure 1.

Focusing on the sample “Xiaoli” (as depicted in Figure 1), we find that five exercise plans including Taiji, Stretching, Rugby, Cycling, and Walking, are recommended. Notably, Taiji and Stretching are forecasted to have the most substantial effect in lowering HbA1c levels, approximately by 1.2%. However, the difference in their recommendation scores, which reflect the proportion of responders, individuals who exhibit a reduction in HbA1c post-exercise intervention, in the nearest neighbor samples, is significant, being 1 and 0.6 respectively. In our analysis, the predicted fluctuation in HbA1c after each exercise plan, such as Stretching, is based on the average fluctuation observed in nearest neighbor samples. Therefore, despite having a lower recommendation score, the predicted reduction in HbA1c from Stretching exercises is similar to that from Taiji. This phenomenon might be attributable to the selection of optimal reference samples for Stretching, where responders demonstrate significantly higher positive reductions in HbA1c.

• Download figure
• Open in new tab

Figure 1:

Visualization of recommended exercise plans for “Xiaoli”.

• Download figure
• Open in new tab

• Download figure
• Open in new tab

Identification of feature panel for each exercise type

To select features meticulously enhances model generalizability and user-friendliness for T2D patients. Our initial step involved conducting an extensive literature review to identify variables associated with different forms of exercise.^{(11),(12),(13)} Subsequently, employing a 70-30 split strategy, we divided the entire dataset into training and test sets, ensuring a robust framework for model development and validation. Within this structured approach, we applied the LASSO technique to the training set to refine our variable selection, focusing on those with non-zero coefficients (Supplementary Table 1.1).

Figure 2 illustrates the finalized comprehensive feature panel, which includes seven critical variables: Age, Sex, BMI, Duration_T2D (years), HDL (mmol/L), PCS, and WHtR. Among these, Age, previously reported as a significant variable for adjusting exercise dosing, ⁽¹⁴⁾ emerges as the most frequently occurring variable across various exercise types, highlighting its critical role in tailoring personalized exercise plans for T2D patients.

• Download figure
• Open in new tab

Figure 2:

The finalized feature panel.

Nearest neighbor-based predictions for HbA1c fluctuation post-exercise intervention

Following the construction of the feature panel, we aimed to evaluate the effects of various exercise plans on individuals with Type 2 Diabetes (T2D) using the k-Nearest Neighbors (kNN) algorithm. Initially, we determined the optimal k-value for each exercise plan based on the training set, employing a Leave-One-Out strategy. This critical step facilitates the TailorExercisePlan function to identify the k most similar reference samples to a given patient, using the average fluctuation in HbA1c values post-exercise intervention as a reference for potential outcomes. Additionally, based on exercise plans detailed in the literature from which our raw data were derived, specific exercise plans and predicted reductions in HbA1c are provided.

To confirm the effectiveness of our comprehensive exercise recommendation framework, we utilized calibration curves. Analysis depicted in Figure 3 reveals a significant linear correlation between predicted and actual HbA1c reduction values across both the training and preserved test sets (R² = 0.28, MAD = 0.79, p = 5.29e-8 in the training set; R² = 0.29, MAD = 0.89, p = 2.09e-04 in the test set), solidifying the framework’s predictive accuracy.

• Download figure
• Open in new tab

Figure 3:

Calibration curves on (A) training set and (B) test set.

These findings collectively affirm the robustness and efficacy of the T2DFitTailor framework in tailoring exercise plans for T2D management. Through a strategic blend of feature panel selection, nearest neighbor-based predictions, and rigorous validation, we demonstrate the potential of a well-structured model to not only predict HbA1c reductions effectively but also to tailor exercise recommendations that resonate with the dynamic needs of T2D patients, thereby enhancing both the practicality and accessibility of diabetes care.