Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 26 June 2025

Sec. Precision Medicine

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1574546

A study on the risk prediction model for venous thromboembolism in orthopedic inpatients based on machine learning

Bo Zhang&#x;Bo Zhang1Yumei Qin&#x;Yumei Qin2Liandi JiuLiandi Jiu1Chunming QinChunming Qin2Jiangbo WangJiangbo Wang2Haiqing Zhao
Haiqing Zhao2*
  • 1Digital Health China Technologies Co., Ltd., Beijing, China
  • 2Nanxishan Hospital of Guangxi Zhuang Autonomous Region, The Second People’s Hospital of Guangxi Zhuang Autonomous Region, Guilin, China

Objective: To construct a venous thromboembolism (VTE) risk prediction model for orthopedic inpatients using machine learning modeling techniques, identify high-risk patients, and optimize clinical interventions.

Methods: This study involved a retrospective analysis of 286 orthopedic inpatients from Nanxishan Hospital of Guangxi Zhuang Autonomous Region (The Second People’s Hospital of Guangxi Zhuang Autonomous Region) from January 1, 2022 to December 31, 2022. To ensure patient information security, all data were fully anonymized before access. The collected data included basic information such as gender, age, ethnicity, and body mass index (BMI), lifestyle factors and medical history (including smoking, alcohol use, diabetes, hypertension, and personal and family history of VTE), clinical test results (such as thrombin time, plasma D-dimer, total bilirubin, and urinary protein via dry chemistry), as well as genetic test results related to VTE risk. Feature analysis and data mining were conducted, and eight different machine learning algorithms were used to build the prediction model. The SHapley Additive exPlanation (SHAP) method was used to rank the feature importance and explain the final model.

Results: Through a comprehensive evaluation and comparison of eight different machine learning models, the results clearly indicate that the XGBoost model outperforms the others across all performance metrics, achieving the highest accuracy of 0.828 and AUROC of 0.931, significantly surpassing the other models, particularly in prediction accuracy and discriminative ability. Compared to the traditional Caprini scoring model, XGBoost not only shows improvements in accuracy and specificity but also demonstrates a significant increase in Area Under the Curve (AUC), further validating its superior performance in VTE risk prediction.

Conclusion: This model can be effectively used for early risk prediction of VTE, helping to reduce the incidence of venous thromboembolism in orthopedic patients. Given its promising results, further validation and wider application of the model in clinical settings are warranted to enhance patient outcomes and improve preventive strategies.

Introduction

Venous thromboembolism (VTE), which includes deep vein thrombosis and pulmonary embolism, is the third most common cardiovascular disease worldwide, following myocardial infarction and stroke (1, 2). VTE is especially prevalent among hospitalized patients. In China, studies show that as many as 45.2% of hospitalized patients are at high risk for VTE, with 53.4% of surgical patients facing elevated risk (3). Orthopedic patients, in particular, are at significantly higher risk due to factors such as surgery, prolonged immobility, and common comorbidities (46). This puts them at a higher incidence and mortality rate for VTE, placing a considerable physical and economic burden on both patients and their families. Despite the availability of effective preventive measures, including anticoagulation therapy and mechanical prevention, the incidence of VTE remains high (7, 8). Given its high mortality rate and severe complications, early identification and accurate assessment of VTE risk, followed by personalized prevention strategies, is a critical clinical challenge. Therefore, research on VTE risk prediction for orthopedic inpatients is of urgent importance.

Currently, VTE risk assessment primarily relies on clinical experience and standardized scoring tools such as the Caprini, Padua, and Khorana scores (911). These scales evaluate risk based on a range of known factors, including age, gender, body mass index (BMI), and comorbidities. The Caprini score, in particular, is widely used for assessing VTE risk in surgical patients, especially in orthopedic inpatients (1214). However, despite their widespread use, these scoring systems have limited value in guiding specific preventive measures and often fail to fully account for individual patient differences and the complexity of clinical situations.

In recent years, with the rapid development of artificial intelligence, machine learning has made significant progress in various fields (15) such as disease risk prediction (16), drug dosage individualization (17), and treatment outcome evaluation (18). Machine learning can handle complex nonlinear relationships and extract potential key factors from vast amounts of data, thereby improving prediction accuracy. At the same time, the widespread use of Electronic Medical Record (EMR) systems in hospitals has made the collection of clinical data more precise and convenient, providing reliable data support for machine learning modeling. As a result, machine learning methods based on EMRs have gradually gained the attention and recognition of clinicians (19, 20).

In this study, we aimed to develop and validate explainable machine learning models for early and accurate prediction of VTE in orthopedic inpatients by analyzing their clinical characteristics, medical history, laboratory results, and genetic testing data. Based on this risk assessment, appropriate interventions will be implemented according to different risk stratifications to reduce the incidence of VTE-related complications. Through predictive analysis using machine learning, we aim to improve clinical outcomes, optimize healthcare resource utilization, enhance patient safety, and improve the quality of care.

In conclusion, personalized VTE risk assessment tools represent a significant advancement in the management of surgical patients. By integrating modern machine learning technologies, we aim to bridge the gap between traditional risk assessment methods and the needs of high-risk patients, supporting precision medicine and individualized care. The findings of this study have the potential to transform current VTE management practices, making them more aligned with patients’ specific needs, and driving the medical field toward a more precise and efficient future.

Methods

Study population

This is a single-center retrospective cohort study, with subjects consisting of 286 orthopedic inpatients at Nanxishan Hospital of Guangxi Zhuang Autonomous Region (The Second People’s Hospital of Guangxi Zhuang Autonomous Region) from January 1, 2022 to December 31, 2022. The inclusion criteria included the following: (1) aged 18 years or older; (2) orthopedic inpatients with a hospital stay > 3 days; (3) completed VTE risk gene polymorphism assessment; (4) no contraindications to anticoagulation; (5) voluntarily agreed to participate in the study and signed an informed consent form. The exclusion criteria were as following: (1) patients who were bedridden or had restricted mobility (e.g., hemiplegia) prior to admission; (2) patients with renal or hepatic dysfunction; (3) patients with hematologic disorders or coagulation dysfunction; (4) pregnant or breastfeeding women; (5) patients with severely missing gene or clinical phenotype data. After applying these criteria, the final study cohort was selected, ensuring the representativeness and scientific rigor of the research findings.

Data collection and processing

We collected demographic information (such as gender, age, race, BMI), lifestyle factors and medical history (including smoking, alcohol consumption, diabetes, hypertension, VTE history, family history of VTE), laboratory test results (such as thrombin time, plasma D-dimer, total bilirubin, urinary protein by dry chemistry, etc.), as well as genetic polymorphism data to identify and select features associated with VTE and construct a risk prediction model. All data were obtained from the EMR system.

Firstly, features with more than 40% missing values were excluded from subsequent analyses to minimize potential bias. A total of 34 features, including age, body mass index (BMI), sex, and others, were ultimately retained and missing data were addressed using median imputation. Outliers for each feature were identified using the IQR (Interquartile Range) method and replaced with the corresponding feature’s median value. Finally, Min-Max scaling was applied to normalize the data, rescaling it to a range of 0 to 1.

Feature selection

Selecting the most relevant and impactful features from the original dataset not only improves model performance and interpretability but also reduces storage and computational resource requirements.

Firstly, this study identified the optimal feature subset using Recursive Feature Elimination (RFE), based on Random Forest model and XGBoost model, to reduce dimensionality and improve model performance. Subsequently, by comparing the differences in each feature between the VTE group and the non-VTE group, features without significant differences (p < 0.1) were excluded.

Additionally, due to the potential impact of multicollinearity among features on prediction accuracy, when two features were highly correlated (correlation coefficient > 0.9) in Spearman’s correlation analysis, the feature less correlated with the outcome was eliminated from the subset, as shown in Figure 1. Finally, combining insights from VTE-related literature, the final features used to construct the model were determined.

Figure 1
A heatmap showing the correlation matrix of various medical variables. The variables are listed along both axes, such as BMI, age, and blood glucose. Correlation values range from -1 to 1, with color gradients from blue (negative correlation) to red (positive correlation). The diagonal displays perfect correlation values of 1. A color scale on the right indicates correlation strength.

Figure 1. Feature correlation heatmap. VTE, venous thromboembolism; BMI, body mass index; GOT, glutamic-oxaloacetic transaminase; APTT, activated partial thromboplastin time.

Model construction and selection

To develop a VTE risk prediction model, this study randomly selected 80% of the dataset as the training set for model training, while the remaining 20% was used as the test set for model performance evaluation (internal validation). To ensure class balance between positive and negative samples in the training set, Synthetic Minority Oversampling Technique (SMOTE) was applied to the training data.

A total of eight binary classification machine learning models were constructed to predict the VTE risk in orthopedic inpatients, including Naive Bayes (NB), K-Nearest Neighbors (KNN), Support Vector Machines (SVM), Logistic Regression (LR), Decision Tree (DT), Adaptive Boosting (AdaBoost), eXtreme Gradient Boosting (XGBoost), and Random Forest (RF).

To enhance model robustness and reduce the risk of overfitting, five-fold cross-validation was employed for training. Model performance was comprehensively evaluated using various metrics, including Accuracy, Sensitivity, Specificity, Precision, F1 Score, and AUC. The XGBoost model and Random Forest model, recognized as the two best-performing predictive model, is utilized for further model optimization.

Hyperparameter tuning and internal validation

Hyperparameter tuning can optimize the performance and generalization ability of the model. A combination of random search and manual fine-tuning was employed to optimize hyperparameters, ensuring the model achieved its best performance.

An internal data set consisting of 58 samples, including 41 negative cases and 17 positive cases was employed for the internal validation.

Model explanation

To enhance the transparency and reliability of the risk assessment model, this study uses SHAP (SHapley Additive exPlanations) to interpret the model’s predictions. SHAP is effective because it provides both a broad overview of feature importance and a detailed explanation of individual predictions. At the global level, SHAP evaluates the contribution of each feature across all samples, highlighting the most influential factors driving the model’s decisions. This insight aids in model optimization and identifying key features. At the local level, SHAP shows how each feature influences a specific prediction. This helps clarify how particular factors impact an individual patient’s risk score, providing valuable guidance for targeted interventions.

Statistical analysis

Continuous variables are presented as mean (standard deviation), while categorical variables are shown as counts and percentages. For the comparative analysis between negative and positive subgroups, continuous variables were assessed using the Mann–Whitney U test or T-test, and categorical variables were analyzed using the Chi-square test. Spearman correlation analysis was used to evaluate the relationships between continuous variables. The predictive power was assessed using the Area Under the Curve (AUC). All data analyses were conducted using Python 3.8.3 and the scikit-learn library (version 1.3.2). The corresponding source code is publicly available on GitHub at https://github.com/LDjiu/VTE_predict.

Results

Patient characteristics

This retrospective study included 286 patients, comprising 86 patients who developed VTE during hospitalization and 200 patients who did not. Among 973 patients admitted to the orthopedic department of Nanxishan Hospital from January 1, 2022 to December 31, 2022, 687 patients were excluded, including 15 patients who were under 18 years old, 43 patients whose hospital stay of less than 3 days, and 629 patients who did not undergo genetic testing for VTE risk. The 286 patients were divided into independent training and testing sets. The detailed design of the study was shown in Figure 2.

Figure 2
Flowchart illustrating a process in three stages: data preparation, model development, and model explanation. Data preparation involves selecting 286 eligible patients from 973, excluding 687 based on criteria, and preprocessing data. Model development includes dividing data into training (80%) and testing (20%) subsets, feature selection, and comparing machine learning models. The best-performing model is selected and validated. Model explanation involves global and local interpretability, using SHAP to explain the final model.

Figure 2. Flow diagram for patients screening, data processing, model development and model explanation.

Demographic characteristics, lifestyle habits, medical history, and laboratory test results were collected. Features with more than 40% missing values were removed, leaving 34 features. A comparison of these 34 features between VTE and non-VTE patients is detailed in Table 1. Continuous variables were described using mean (standard deviation), and categorical variables were described using frequencies. Significant differences (p < 0.05) were observed between VTE and non-VTE patients in terms of VTE history, comorbidities, total protein, albumin, D-dimer, and erythrocyte sedimentation rate.

Table 1
www.frontiersin.org

Table 1. Comparison of demographic and clinical characteristics and outcomes between non-VTE and VTE patients.

Feature selection

After removing features with more than 40% missing values, 34 features remained, including: polygenic risk, gender, race, age, BMI, smoking history, alcohol consumption history, VTE history, comorbidities, urine protein, urine leukocytes, urine specific gravity, glutamic-oxaloacetic transaminase (GOT), glutamic pyruvic transaminase, gamma-glutamyl transferase, alkaline phosphatase, total bilirubin, total protein, albumin, serum creatinine, blood urea nitrogen, serum uric acid, globulin, blood glucose, activated partial thromboplastin time (APTT), D-dimer, fibrinogen, thrombin time, prothrombin activity, prothrombin time, international normalized ratio, erythrocyte sedimentation rate, surgery, bed rest duration.

Based on Random Forest model and XGBoost model, Recursive Feature Elimination (RFE) was applied to select the optimal feature subset, as shown in Figure 3. Among the two models, a feature subset with 11 features achieved the optimal performance on the XGBoost model. The 11-feature subset consisted of D-dimer, erythrocyte sedimentation rate, blood glucose, urine specific gravity, serum creatinine, urine leukocytes, VTE history, activated partial thromboplastin time (APTT), bed rest duration, gender and age. Due to without differences (p < 0.1) between VTE and non-VTE patients, five features, including bed rest duration, blood glucose, urine specific gravity, serum creatinine and urine leukocytes, were eliminated, as shown in Table 1. We also substituted gender with comorbidities that showed stronger association with the outcome, as shown in Figure 1. Given that previous literature has indicated that the genetic polymorphisms MTHFR (C677T) and PAI-1(4G/5G) are associated with an increased risk of VTE (2123), the polygenic risk based on these two polymorphic sites were included as features in the model construction. Although statistical analysis did not show a significant difference between the VTE and non-VTE groups, these genetic characteristics still hold important value in VTE risk assessment. Additionally, according to the literature, bed rest duration is also a critical factor influencing the occurrence of VTE.

Figure 3
Line chart showing feature reduction of models with AUROC scores on the y-axis and number of features on the x-axis. Random forest and XGBoost are compared. Both models peak around 11 features, then decline. Random forest is in red, XGBoost in blue. A vertical dashed line marks 11 features.

Figure 3. The performance of the models during recursive feature elimination. XGboost, eXtreme gradient boosting.

Finally, 8 features were selected to construct the model: age, VTE history, comorbidities, APTT, D-dimer, erythrocyte sedimentation rate, bed rest duration, and polygenic risk.

Model construction and performance comparison

Based on the training dataset, eight machine learning prediction models were constructed, including SVM, LR, KNN, XGBoost, AdaBoost, RF, DT, and NB. Table 2 shows the average performance of these models under five-fold cross-validation, and the ROC curves are presented in Figure 4.

Table 2
www.frontiersin.org

Table 2. Performance of eight machine learning models and caprini score on the training set.

Figure 4
ROC curve comparing multiple models, showing true positive rate (TPR) versus false positive rate (FPR). Models include SVM (AUC=0.782), LR (AUC=0.765), and RF (AUC=0.873), among others. XGBoost (AUC=0.869) performs strongly, while DT (AUC=0.696) shows weaker performance.

Figure 4. Receiver operating characteristic curves of eight machine learning models and caprini score on the training set. NB, naive bayes; KNN, k-nearest neighbors; SVM, support vector machines; LR, logistic regression; DT, decision tree; AdaBoost, adaptive boosting; XGboost, eXtreme gradient boosting; RF, Random Forest; AUC, area under the curve.

Among the eight models, Random Forest (RF) and XGBoost outperformed the others, achieving the best predictive performance with AUROCs of 0.873 and 0.869, respectively, making them the two best-performing predictive models selected for further hyperparameter tuning.

Additionally, when comparing the eight machine learning models with the traditional Caprini score, seven of the models (except Decision Tree) achieved higher AUROCs than the Caprini score. Moreover, the accuracy, specificity, and precision of all eight machine learning models exceeded those of the Caprini score.

Hyperparameter tuning and validation of the models

The combination of different parameters can directly impact the predictive ability, generalization performance, and practical applicability of a model. In this study, the optimal hyperparameter combination was obtained through random search and manual fine-tuning, resulting in the following settings of the XGBoost model: subsample of 0.72, n_estimators of 50, min_child_weight of 1, max_depth of 6, learning_rate of 0.074, gamma of 0.25, and colsample_bytree of 0.5, and the following settings of the Random Forest model: n_estimators of 300, min_samples_split of 2, max_features of 3, max_depth of 5, and bootstrap of True. The discriminative ability of the two model on the test set was shown in Figure 5, revealing that the XGBoost model had the best predictive performance. So, the XGBoost model was selectd as the final model.

Figure 5
Receiver Operating Characteristic (ROC) curve comparing three models: XGBoost (red, area = 0.931), RF (blue, area = 0.872), and Caprini (purple, area = 0.737). The plot shows true positive rate against false positive rate, with a baseline diagonal line indicating random performance.

Figure 5. Receiver operating characteristic curves of the two best-performing machine learning models and caprini score on internal validation set.

The test set consisted of 58 samples, including 41 negative cases and 17 positive cases. The final model performed well on the test set, as shown in Figures 5, 6A. The final model achieved an accuracy of 0.828, sensitivity of 0.824, specificity of 0.829, precision of 0.667, F1 score of 0.737, and AUROC of 0.931. Among the 41 negative samples, 34 were correctly predicted, while 14 of the 17 positive samples were correctly predicted.

Figure 6
Confusion matrices labeled A and B comparing true versus predicted values for

Figure 6. The performance of the final model and caprini score on internal validation set. (A) Confusion matrix of the final model on internal validation set. (B) Confusion matrix of caprini score on internal validation set.

Additionally, the final model outperformed the Caprini score on the test set, as detailed in Figures 5, 6. The Caprini score achieved an AUROC of 0.737, correctly predicting 9 out of 41 negative samples and 17 out of 17 positive samples. The Caprini score tended to overestimate low-risk cases as high risk for VTE.

Model explanation

Since clinicians often find it difficult to accept predictive models that are not directly interpretable or understandable, the SHAP method was employed to explain the output of the final model by quantifying each feature’s contribution to the prediction. The main advantage of SHAP lies in its ability to provide both global and local interpretability. Global interpretation highlights the most influential features in the model’s decision-making process. Figures 7A,B present SHAP summary plots, where the SHAP mean value represents each feature’s contribution to the model predictions, ranked in descending order of importance. The order of importance is as follows: D-dimer, VTE history, erythrocyte sedimentation rate, APTT, bed rest duration, polygenic risk, age, and comorbidities.

Figure 7
Panel (A) displays a scatter plot of SHAP values for various features related to VTE risk, with color indicating feature value from low (blue) to high (pink). Panel (B) is a bar chart showing mean SHAP values for features like D-dimer and VTE history, highlighting their importance in the model.

Figure 7. Global interpretability for the final model prediction by SHAP method. (A) SHAP summary plot of 8 features included in the final prediction model for VTE. The horizontal axis represents the SHAP values, and the vertical axis represents the features. Each point corresponds to a sample, with the color of the points indicating the magnitude of the feature values: red represents higher feature values, while blue represents lower feature values. A positive SHAP value indicates an increased risk of VTE. Additionally, a positive correlation is shown when higher feature values result in larger SHAP values. (B) SHAP summary plot of 8 features ranked by the mean absolute SHAP values across all samples, representing the average impact of each feature on the prediction of VTE. VTE, venous thromboembolism; APTT, activated partial thromboplastin time.

The model predicts outcomes for specific individuals by assigning a SHAP value to each feature. As shown in Figures 8A,B, predictions for negative and positive individuals are visualized. The length of the bars represents the magnitude of the feature’s impact on the final prediction, with red bars indicating a positive contribution and blue bars indicating a negative contribution. For positive individuals, APTT, bed rest duration, and D-dimer drive the model to predict VTE; for negative individuals, erythrocyte sedimentation rate, VTE history, and D-dimer drive the model to predict non-VTE.

Figure 8
Graphical representation of two probabilistic profiles, labeled A and B. In A, the probability of a higher value is represented in red, including APTT at 20.8, bed rest duration at 10.0, and D-dimer at 1200.0, while the probability of a lower value is shown in blue, with erythrocyte sedimentation rate at 32.0. The overall function value is 1.62. In B, red indicates higher values for bed rest duration at 12.0, and APTT at 26.8. Blue indicates lower values for D-dimer at 610.0, erythrocyte sedimentation rate at 45.0, and VTE history at 0.0, with a function value of -2.09.

Figure 8. Local interpretability for the final prediction model by SHAP method. (A) Force plot of a patient with VTE, APTT, bed rest duration and D-dimer are major features contributing to a higher predicted risk of VTE. (B) Force plot of a patient without VTE, D-dimer, erythrocyte sedimentation rate and VTE history are major features contributing to a lower predicted risk of VTE. APTT, activated partial thromboplastin time; VTE, venous thromboembolis.

Discussion

In this study, we extracted data from the electronic medical record system of orthopedic inpatients and applied eight machine learning algorithms to build predictive models for the risk of VTE. The goal was to identify the optimal model for predicting the occurrence of VTE. Additionally, to validate the performance advantages of machine learning methods, we also used the traditional Caprini score to predict VTE risk in the same group of patients. A comparative analysis of the results revealed that the performance of the Caprini score model was significantly lower than the optimal XGBoost model. Specifically, the Caprini score model has high sensitivity but very low specificity, meaning a significant number of patients are misclassified as having a VTE risk, resulting in a high false positive rate. Even though the final AUROC is 0.737, the Caprini score performs poorly in predicting negative results. In contrast, the evaluation metrics of the XGBoost model were as follows: sensitivity of 0.824, specificity of 0.829, and an AUROC of 0.931, indicating that the XGBoost model demonstrated superior performance in identifying VTE risk, providing more accurate and effective support for clinical decision-making.

In addition, the innovation of this study lies in incorporating genetic factors related to VTE in the Chinese population into the predictive model, combined with clinical phenotype data to construct a comprehensive risk prediction model. Previous studies have shown that Factor V Leiden (rs6025) and prothrombin G20210A mutations are typical genetic risk factors for VTE (24, 25); however, these mutations are relatively rare in the Chinese population (26). Therefore, this study selected the MTHFR (C677T) and PAI-1(4G/5G) variants, which are more common in the Chinese population, for investigation (2123). Mutations in the MTHFR gene can lead to hyperhomocysteinemia, which in turn causes endothelial injury, increasing the risk of VTE (27, 28). Variations in the PAI-1 gene inhibit the fibrinolytic system, resulting in fibrinolysis dysfunction, which becomes an important trigger for thrombosis (29, 30). By incorporating the genetic variations of MTHFR and PAI-1 into the model, this study aims to further improve the accuracy and clinical applicability of the VTE risk prediction model.

However, this study has several limitations. First, the relatively small sample size may not fully capture the diverse risk characteristics of all orthopedic inpatients, which limits the broader applicability of the findings. Second, the study was conducted at a single center and lacks multi-center validation, which could affect the model’s generalizability across different clinical settings. Third, due to substantial missing data for some variables, these features could not be included in the analysis, which may restrict the model’s ability to fully elucidate the complex biological mechanisms underlying venous thromboembolism (VTE). As such, while the findings of this study show promising clinical potential, caution is necessary when interpreting the results. Future research should focus on larger, multi-center cohort studies and integrate additional biomarkers to address these limitations and enhance the robustness of the model. The ultimate goal is to integrate the model into clinical decision support systems, enabling real-time prediction and intelligent warning of VTE risk based on dynamic clinical data.

In conclusion, this study successfully developed a machine learning-based VTE risk prediction model for orthopedic inpatients. The model demonstrated strong predictive performance on both the training and testing datasets, highlighting its potential for early identification of high-risk patients in clinical practice. These findings underscore the importance of integrating advanced analytical methods into clinical risk assessment, laying the foundation for personalized preventive strategies in VTE management. Future research should focus on validating the model in larger, more heterogeneous populations and exploring the integration of additional clinical and molecular data to further enhance prediction accuracy and utility, improve model generalizability, and ultimately benefit patients, advancing the development of personalized medicine.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by the The Ethics Committee of Guangxi Zhuang Autonomous Region Nanxishan Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because this is a single-center retrospective cohort study.

Author contributions

BZ: Methodology, Writing – original draft, Data curation, Writing – review & editing. YQ: Methodology, Writing – original draft, Data curation, Writing – review & editing. LJ: Data curation, Investigation, Supervision, Writing – original draft. CQ: Investigation, Validation, Supervision, Writing – original draft. JW: Software, Supervision, Validation, Writing – original draft. HZ: Data curation, Project administration, Writing – original draft, Writing – review & editing, Investigation, Methodology, Supervision.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

BZ and LJ were employed by the Digital Health China Technologies Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Gregson, J, Kaptoge, S, Bolton, T, Pennells, L, Willeit, P, Burgess, S, et al. Cardiovascular risk factors associated with venous thromboembolism. JAMA Cardiol. (2019) 4:163–73. doi: 10.1001/jamacardio.2018.4537

PubMed Abstract | Crossref Full Text | Google Scholar

2. Goldhaber, SZ. Venous thromboembolism: epidemiology and magnitude of the problem. Best Pract Res Clin Haematol. (2012) 25:235–42. doi: 10.1016/j.beha.2012.06.007

PubMed Abstract | Crossref Full Text | Google Scholar

3. Zhai, Z, Kan, Q, Li, W, Qin, X, Qu, J, Shi, Y, et al. VTE risk profiles and prophylaxis in medical and surgical inpatients: the identification of Chinese hospitalized patients' risk profile for venous thromboembolism (DissolVE-2)-a cross-sectional study. Chest. (2019) 155:114–22. doi: 10.1016/j.chest.2018.09.020

PubMed Abstract | Crossref Full Text | Google Scholar

4. Whiting, PS, White-Dzuro, GA, Greenberg, SE, VanHouten, J, Avilucea, FR, Obremskey, WT, et al. Risk factors for deep venous thrombosis following Orthopaedic trauma surgery: an analysis of 56,000 patients. Arch Trauma Res. (2016) 5:e32915. doi: 10.5812/atr.32915

PubMed Abstract | Crossref Full Text | Google Scholar

5. Kahn, SR, and Shivakumar, S. What's new in VTE risk and prevention in orthopedic surgery. Res Pract Thromb Haemost. (2020) 4:366–76. doi: 10.1002/rth2.12323

PubMed Abstract | Crossref Full Text | Google Scholar

6. Wang, Y, Xu, X, and Zhu, W. Anticoagulant therapy in orthopedic surgery - a review on anticoagulant agents, risk factors, monitoring, and current challenges. J Orthop Surg (Hong Kong). (2024) 32:10225536241233473. doi: 10.1177/10225536241233473

PubMed Abstract | Crossref Full Text | Google Scholar

7. Lutsey, PL, and Zakai, NA. Epidemiology and prevention of venous thromboembolism. Nat Rev Cardiol. (2023) 20:248–62. doi: 10.1038/s41569-022-00787-6

PubMed Abstract | Crossref Full Text | Google Scholar

8. Onwuzo, C, Olukorode, J, Sange, W, Tanna, SJ, Osaghae, OW, Hassan, A, et al. A review of the preventive strategies for venous thromboembolism in hospitalized patients. Cureus. (2023) 15:e48421. doi: 10.7759/cureus.48421

PubMed Abstract | Crossref Full Text | Google Scholar

9. Khorana, AA, Kuderer, NM, Culakova, E, Lyman, GH, and Francis, CW. Development and validation of a predictive model for chemotherapy-associated thrombosis. Blood. (2008) 111:4902–7. doi: 10.1182/blood-2007-10-116327

PubMed Abstract | Crossref Full Text | Google Scholar

10. Barbar, S, Noventa, F, Rossetto, V, Ferrari, A, Brandolin, B, Perlati, M, et al. A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua prediction score. J Thromb Haemost. (2010) 8:2450–7. doi: 10.1111/j.1538-7836.2010.04044.x

PubMed Abstract | Crossref Full Text | Google Scholar

11. Sterbling, HM, Rosen, AK, Hachey, KJ, Vellanki, NS, Hewes, PD, Rao, SR, et al. Caprini risk model decreases venous thromboembolism rates in thoracic surgery Cancer patients. Ann Thorac Surg. (2018) 105:879–85. doi: 10.1016/j.athoracsur.2017.10.013

PubMed Abstract | Crossref Full Text | Google Scholar

12. Wilson, S, Chen, X, Cronin, MA, Dengler, N, Enker, P, Krauss, ES, et al. Thrombosis prophylaxis in surgical patients using the Caprini risk score. Curr Probl Surg. (2022) 59:101221. doi: 10.1016/j.cpsurg.2022.101221

PubMed Abstract | Crossref Full Text | Google Scholar

13. Zhang, X, Hao, A, Lu, Y, and Huang, W. Deep vein thrombosis and validation of the Caprini risk assessment model in Chinese orthopaedic trauma patients: a multi-center retrospective cohort study enrolling 34,893 patients. Eur J Trauma Emerg Surg. (2023) 49:1863–71. doi: 10.1007/s00068-023-02265-1

PubMed Abstract | Crossref Full Text | Google Scholar

14. Lin, Z, Sun, H, Chen, M, Li, D, Cai, Z, Wang, Y, et al. Utilization of the Caprini risk assessment model(RAM) to predict venous thromboembolism after primary hip and knee arthroplasty: an analysis of the healthcare cost and utilization project(HCUP). Thromb J. (2024) 22:68. doi: 10.1186/s12959-024-00633-4

PubMed Abstract | Crossref Full Text | Google Scholar

15. Alowais, SA, Alghamdi, SS, Alsuhebany, N, Alqahtani, T, Alshaya, AI, Almohareb, SN, et al. Revolutionizing healthcare: the role of artificial intelligence in clinical practice. BMC Med Educ. (2023) 23:689. doi: 10.1186/s12909-023-04698-z

PubMed Abstract | Crossref Full Text | Google Scholar

16. Matheson, MB, Kato, Y, Baba, S, Cox, C, Lima, JAC, and Ambale-Venkatesh, B. Cardiovascular risk prediction using machine learning in a large Japanese cohort. Circ Rep. (2022) 4:595–603. doi: 10.1253/circrep.CR-22-0101

PubMed Abstract | Crossref Full Text | Google Scholar

17. Li, QY, Tang, BH, Wu, YE, Yao, BF, Zhang, W, Zheng, Y, et al. Machine learning: a new approach for dose individualization. Clin Pharmacol Ther. (2024) 115:727–44. doi: 10.1002/cpt.3049

PubMed Abstract | Crossref Full Text | Google Scholar

18. Feuerriegel, S, Frauen, D, Melnychuk, V, Schweisthal, J, Hess, K, Curth, A, et al. Causal machine learning for predicting treatment outcomes. Nat Med. (2024) 30:958–68. doi: 10.1038/s41591-024-02902-1

PubMed Abstract | Crossref Full Text | Google Scholar

19. Schwartz, JT, Gao, M, Geng, EA, Mody, KS, Mikhail, CM, and Cho, SK. Applications of machine learning using electronic medical Records in Spine Surgery. Neurospine. (2019) 16:643–53. doi: 10.14245/ns.1938386.193

PubMed Abstract | Crossref Full Text | Google Scholar

20. Hu, J, Xu, J, Li, M, Jiang, Z, Mao, J, Feng, L, et al. Identification and validation of an explainable prediction model of acute kidney injury with prognostic implications in critically ill children: a prospective multicenter cohort study. EClinicalMedicine. (2024) 68:102409. doi: 10.1016/j.eclinm.2023.102409

PubMed Abstract | Crossref Full Text | Google Scholar

21. Wang, J, Wang, C, Chen, N, Shu, C, Guo, X, He, Y, et al. Association between the plasminogen activator inhibitor-1 4G/5G polymorphism and risk of venous thromboembolism: a meta-analysis. Thromb Res. (2014) 134:1241–8. doi: 10.1016/j.thromres.2014.09.035

PubMed Abstract | Crossref Full Text | Google Scholar

22. Zhang, P, Gao, X, Zhang, Y, Hu, Y, Ma, H, Wang, W, et al. Association between MTHFR C677T polymorphism and venous thromboembolism risk in the Chinese population: a meta-analysis of 24 case-controlled studies. Angiology. (2015) 66:422–32. doi: 10.1177/0003319714546368

PubMed Abstract | Crossref Full Text | Google Scholar

23. Wang, B, Xu, P, Shu, Q, Yan, S, and Xu, H. Combined effect of MTHFR C677T and PAI-1 4G/5G polymorphisms on the risk of venous thromboembolism in Chinese lung Cancer patients. Clin Appl Thromb Hemost. (2021) 27:10760296211031291. doi: 10.1177/10760296211031291

PubMed Abstract | Crossref Full Text | Google Scholar

24. Emmerich, J, Rosendaal, FR, Cattaneo, M, Margaglione, M, De Stefano, V, Cumming, T, et al. Combined effect of factor V Leiden and prothrombin 20210A on the risk of venous thromboembolism--pooled analysis of 8 case-control studies including 2310 cases and 3204 controls. Study Group for Pooled-Analysis in Venous Thromboembolism. Thromb Haemost. (2001) 86:809–16. doi: 10.1055/s-0037-1616136

Crossref Full Text | Google Scholar

25. Zoller, B, Svensson, PJ, Dahlbäck, B, Lind-Hallden, C, Hallden, C, Elf, J, et al. Genetic risk factors for venous thromboembolism. Expert Rev Hematol. (2020) 13:971–81. doi: 10.1080/17474086.2020.1804354

PubMed Abstract | Crossref Full Text | Google Scholar

26. Jun, ZJ, Ping, T, Lei, Y, Li, L, Ming, SY, Jing, W, et al. Prevalence of factor V Leiden and prothrombin G20210A mutations in Chinese patients with deep venous thrombosis and pulmonary embolism. Clin Lab Haematol. (2006) 28:111–6. doi: 10.1111/j.1365-2257.2006.00757.x

PubMed Abstract | Crossref Full Text | Google Scholar

27. Park, WC, and Chang, JH. Clinical implications of methylenetetrahydrofolate reductase mutations and plasma homocysteine levels in patients with thromboembolic occlusion. Vasc Specialist Int. (2014) 30:113–9. doi: 10.5758/vsi.2014.30.4.113

PubMed Abstract | Crossref Full Text | Google Scholar

28. Raghubeer, S, and Matsha, TE. Methylenetetrahydrofolate (MTHFR), the one-carbon cycle, and cardiovascular risks. Nutrients. (2021) 13:4562. doi: 10.3390/nu13124562

PubMed Abstract | Crossref Full Text | Google Scholar

29. Mukhopadhyay, S, Johnson, TA, Duru, N, Buzza, MS, Pawar, NR, Sarkar, R, et al. Fibrinolysis and inflammation in venous Thrombus resolution. Front Immunol. (2019) 10:1348. doi: 10.3389/fimmu.2019.01348

PubMed Abstract | Crossref Full Text | Google Scholar

30. Sillen, M, and Declerck, PJ. A narrative review on plasminogen activator Inhibitor-1 and its (Patho)physiological role: to target or not to target? Int J Mol Sci. (2021) 22:2721. doi: 10.3390/ijms22052721

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: venous thromboembolism, machine learning, risk assessment, orthopedic inpatients, clinical decision support

Citation: Zhang B, Qin Y, Jiu L, Qin C, Wang J and Zhao H (2025) A study on the risk prediction model for venous thromboembolism in orthopedic inpatients based on machine learning. Front. Med. 12:1574546. doi: 10.3389/fmed.2025.1574546

Received: 14 February 2025; Accepted: 16 June 2025;
Published: 26 June 2025.

Edited by:

Jiuping Ji, National Cancer Institute at Frederick (NIH), United States

Reviewed by:

Dimitrios Liakopoulos, General Hospital Nice Piraeus Saint Panteleimon, Greece
Lifan Zhang, West China Hospital, Sichuan University, China

Copyright © 2025 Zhang, Qin, Jiu, Qin, Wang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Haiqing Zhao, NjU3OTI1NDUxQHFxLmNvbQ==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.