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

ORIGINAL RESEARCH article

Front. Med., 27 October 2025

Sec. Nephrology

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

Ensemble machine learning for predicting renal function decline in chronic kidney disease: development and external validation

Hong Chen
Hong Chen1*Yuping HuangYuping Huang1Lizhen ChenLizhen Chen2
  • 1Department of Nephrology, the 95th Hospital of Putian in China RongTong Medical Health Corporation, Putian, China
  • 2Department of Rheumatology and Immunology, the 95th Hospital of Putian in China RongTong Medical Health Corporation, Putian, China

Introduction: Chronic kidney disease (CKD) poses a significant global health challenge, requiring timely interventions to manage renal function decline. Traditional predictive models often lack accuracy and generalizability. This study aimed to develop and validate a machine learning model to enhance risk prediction of renal function decline in CKD patients, enabling early and personalized interventions.

Methods: We developed an ensemble machine learning model using Random Forest, XGBoost, and LightGBM algorithms, incorporating advanced feature selection and hyperparameter tuning. The model was trained and validated on data from 1,200 CKD patients across multiple clinics, selected through stringent inclusion and exclusion criteria. Clinical, demographic, and laboratory data were processed with rigorous quality control. Model performance was assessed using area under the curve (AUC), calibration metrics, and five-fold cross-validation, with external validation across three medical centers.

Results: The ensemble model achieved an AUC of 0.89 (95% CI: 0.87-0.91), outperforming traditional Cox models (AUC: 0.82, 95% CI: 0.79-0.85) and standard machine learning approaches (AUC: 0.85, 95% CI: 0.83-0.87). Key predictors identified via SHAP analysis included estimated glomerular filtration rate (eGFR), age, and urinary protein-creatinine ratio. The model demonstrated excellent calibration (slope: 0.96, 95% CI: 0.94-0.98) and robust performance across diverse patient subgroups, with a 60.6% reduction in computational resource use compared to traditional methods.

Discussion: This machine learning model offers a significant advancement in predicting CKD progression, providing a reliable, generalizable tool for early risk stratification. Its superior accuracy and efficiency support integration into clinical workflows, potentially transforming CKD management by enabling proactive, data-driven interventions. Future research should focus on incorporating novel biomarkers and expanding multicenter validation to further enhance clinical applicability.

1 Introduction

Chronic kidney disease (CKD) has become a particularly urgent health challenge worldwide. In developed countries, the annual medical cost for chronic kidney disease exceeds 120 billion US dollars, and the prevalence of this disease is still on the rise, which has put pressure on the global healthcare system (1, 2). Traditional methods used to predict the progression of chronic kidney disease Relying heavily on some scattered clinical indicators and simple linear models, the accuracy of prediction is relatively poor, with an AUC less than 0.75 (35). Moreover, the universality of this prediction method is also relatively limited among different patient groups. Compared with traditional statistical methods (68), the accuracy of machine learning in healthcare applications has increased by 15 to 20%. However, the existing chronic kidney disease prediction models have some key problems, such as a lack of interpretability, inability to capture the dynamic changes of the disease over time, and insufficient validation in different clinical Settings (911).

In the early research on the risk prediction of chronic kidney disease, researchers used traditional statistical methods to predict CKD. In (12), the authors constructed a structural equation model for risk prediction to predict CKD. In (13), the author selected factors associated with renal failure from a large number of variables and then established a Cox proportional hazards regression model, using this model to predict and evaluate the risk of renal failure. Although traditional statistical methods can predict the risk of CKD, their accuracy is relatively low. With the continuous development of machine learning, researchers have begun to explore the application of machine learning prediction methods. In (14), the authors used machine learning techniques such as random forests and decision trees, effectively improving the performance of prediction. In (15), the authors combined five different machine learning methods, such as Naive Bayes and random Forest, with feature selection techniques and ensemble learning, and used SHAP and LIME to demonstrate the visualization of personalized CKD prediction models, thereby enhancing the interpretability of the models. It has provided a brand-new perspective for CKD medical research. In (16), the authors trained the medical records of 400 patients using different machine learning methods such as Cat Boost, AdaBoost, and Extra Trees. Finally, the accuracy rate reached 97.5%, which shows that the ensemble learning model has potential in the early diagnosis of CKD. In (17), the author proposes an interpretability strategy that uses five machine learning methods to predict CDK datasets and utilizes LIME features to enhance the interpretability of the model. Our code is publicly available at: https://gitee.com/forest-AI/CDK-Model.

This study addresses these fundamental challenges by leveraging four key innovative points, which enable CKD risk prediction to exceed the current capacity. We introduce a brand-new temporal feature engineering framework (18), which can systematically capture the short-term changes and long-term development trends of the disease. It has made great progress compared with the static snapshot methods used to describe existing models in the past. The previous static snapshot methods were rather limited. However, this new framework enables the model not only to simply assess the risk status at a certain moment but also to understand the progression pattern of diseases. We have optimized the integration architecture, which has significantly reduced the demand for computing resources by 60.6% while still maintaining a good prediction effect. Nowadays, many complex machine learning systems encounter some practical implementation obstacles in clinical applications. Our optimization directly addresses these issues. We have developed a comprehensive multi-center external validation strategy in three medical centers and conducted detailed analyses of resource utilization and scalability. In theoretical machine learning research, Most of the time, there is a lack of strong evidence regarding real-world deployment, and our study provides such evidence. We designed from a clinical perspective, combining domain knowledge with advanced feature selection methods to create an interpretable decision support framework. There is a critical gap between complex computational methods and actual clinical applications, and this framework fills this gap.

2 Methods

2.1 Study subjects

Choosing the right study population for the creation of a robust machine learning model to predict decline in renal function is quite the painstaking process. To achieve greater generalizability and external validity, our multi-center study devised a selection protocol to create a representative dataset with greater accuracy (9, 10).

Screening of the population’s initial sample involved 2,500 potential candidates across five tertiary health care centers referred to in Figure 1. Inclusion criteria were methodically developed based on clinical guidelines, specifically the Kidney Disease: Improving Global Outcomes (KDIGO) 2012 guidelines. Eligible participants were adults aged 18–75 years with documented chronic kidney disease (CKD), defined as an estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73m2 for at least 3 months, confirmed by at least two measurements, or the presence of persistent proteinuria (urine protein-to-creatinine ratio [UPCR] ≥ 0.2 g/g for at least 3 months) or other markers of kidney damage (e.g., abnormal renal imaging or biopsy findings) as recorded in electronic health records (EHRs) with standardized diagnostic codes (e.g., ICD-10 codes N18.1-N18.5). This operational definition ensures that CKD diagnosis is not solely reliant on eGFR but incorporates additional clinical and laboratory evidence consistent with KDIGO criteria, enhancing diagnostic specificity and reproducibility. Several longitudinal record requirements were established: a minimum of 2 years of electronic health records (EHRs) spanning from January 2021 to December 2024, and a minimum of four serum creatinine tests conducted within the 12 months prior to the study’s end date (December 2024). These requirements ensured robust longitudinal data to capture renal function trends while reflecting contemporary clinical practices and standardized assay technologies during the study period. After applying these criteria, 1,800 participants qualified for further evaluation.

Figure 1
Flowchart detailing participant selection for a study. Initial screening includes 2,500 individuals. Inclusion criteria: age 18-75, documented CKD, over two years of medical records, complete clinical data, informed consent, totaling 1,800. Exclusion criteria: acute kidney injury, active cancer treatment, severe comorbidities, incomplete or inconsistent records, reducing the sample by 600. Final analyzed sample is 1,200.

Figure 1. Study flow diagram.

To protect data quality and minimize potential confounding factors, a detailed exclusion criterion was employed (19). Some crucial exclusion criteria included active malignancy or chemotherapy within the previous 3 months of the study, kidney transplant, acute renal failure, some other severe comorbid condition that may affect renal function, and poor or incomplete medical file documentation. After applying these criteria, 600 participants were removed as a result of the process.

With the criteria applied, the total number of remaining participants is 1,200, determined through a power analysis (0.05, 0.10) based on expected model performance and complexity. Here, α = 0.05 represents the significance level (Type I error rate, false positive rate), i.e., the probability threshold for rejecting the null hypothesis. β = 0.10 represents the Type II error rate (false negative rate), corresponding to a statistical power of 1–β = 0.90 (90% power), i.e., the probability of correctly detecting a true effect (20). This comprehensive and systematic selection, along with stringent inclusion and exclusion criteria, enhances the quality of the constructed dataset, making it highly suitable for developing and validating complex predictive models for kidney function decline. Only participants clinically and demographically suitable for the model can be considered the “target population,” thereby increasing its utility in clinical practice.

2.2 Data collection and processing

A thorough strategy of data collection and processing was constructed in order to provide clean and usable data for the machine learning model. The quality control system comprised three sequential phases: data collection, preprocessing, feature engineering, with quality checks integrated at each stage (Figure 2).

Figure 2
Flowchart depicting the data processing workflow. It begins with data collection, including clinical features and laboratory tests, followed by data preprocessing involving missing value handling and normalization. Feature engineering comprises feature selection and extraction. Detailed processing steps include raw data cleaning, processing, feature selection, and culminate in the final dataset.

Figure 2. Data processing pipeline.

The data collection protocol was meticulously designed to ensure consistency and accuracy across all participating centers. The study employed a retrospective data collection approach, leveraging electronic health records (EHRs) from five tertiary healthcare centers. The data spanned a period from January 2021 to December 2024, capturing comprehensive clinical and laboratory parameters relevant to kidney function assessment and laboratory parameters which were defined according to existing protocols for assessing kidney function (21, 22). The clinical data consisted of demographic features, comorbid conditions, medication usage, and other clinically relevant observations, and were extracted through standardized electronic health record protocols. Laboratory measurements included comprehensive metabolic panels, complete blood counts, and specific renal function measurements such as serum creatinine, eGFR, and urine protein to creatinine ratio. Additional biochemical parameters such as hemoglobin, albumin, and electrolytes were collected to capture the multifaceted nature of kidney disease progression.

Data preprocessing was executed within a strict quality assurance framework, as described in (23). Participant flow was meticulously tracked, with 2,500 potential candidates initially screened across five tertiary healthcare centers, resulting in 1,800 eligible participants after applying inclusion criteria and 1,200 final participants after exclusion criteria were enforced (see Figure 1 for the study flow diagram). Missing data patterns were analyzed, revealing that missingness was primarily missing at random (MAR), with serum creatinine and urine protein-to-creatinine ratio (UPCR) missing in approximately 8 and 12% of cases, respectively, due to variations in clinical testing frequency. Advanced imputation methods were employed: multiple imputation by chained equations (MICE) for continuous data and mode imputation for categorical data, with validation tests confirming imputation precision (mean absolute error <5% for continuous variables). Continuous data was normalized with z-scores to ensure comparability, and categorical data was encoded using preservation-optimized schemes, such as one-hot encoding for nominal variables. Outlier detection and validation were performed through a combination of statistical techniques (e.g., interquartile range method) and clinical judgment to ensure clinical plausibility. Outcome assessment was conducted using a blinded approach, where evaluators determining renal function decline (defined as eGFR decline ≥30% or progression to dialysis) were unaware of the model’s predictions to minimize bias. Predictors, including eGFR, age, UPCR, comorbidities, and serum creatinine, were pre-specified based on clinical guidelines (KDIGO 2012) and prior literature (1, 13), ensuring alignment with established nephrology knowledge. This comprehensive preprocessing strategy, coupled with rigorous quality checks, ensured data integrity and supported robust model development.

Feature engineering integrates knowledge from related fields to enhance the performance of prediction. As stated in references (24, 25), we have constructed several temporal features, such as the time intervals between consecutive serum creatinine tests, represented by Δt, the changes in serum creatinine, represented by ΔSCr, and the rate of change of eGFR over time. Represented by ΔeGFR/Δt, these features are relied upon to capture the dynamic changes of renal function. In addition, we have also created interaction features, such as the interaction term between age and eGFR, represented as age × eGFR, and the interaction between urine protein-creatinine ratio, that is, UPCR and diabetes status, represented as UPCR × diabetes. These are relied upon to assess how different factors interact with each other. After introducing these features, the ability of the model to predict renal function deterioration has been significantly improved. We also noticed that if the proportion of missing data input is relatively high, it may cause bias. We calculated the False Negative Rate, that is, the False negative Rate, and used this to evaluate the performance of the model in high-risk patients. And optimize the model to reduce the possibility of false negatives. This integrated approach has enhanced the prediction accuracy of the model, made it more reliable, and made it more practical in clinical applications.

The quality control policy is in line with the modern standards of machine learning, as cited in (26). For instance, it will check the quality of the data, record every change made to the data, and also record all corresponding databases, as well as apply automatic verification programs. This approach can ensure the replicability of the research. It also provides a foundation for the sustainable improvement of the model and lays the groundwork for future machine learning analysis. To ensure the consistency and coherence of the research, all participating centers followed the standardized data collection protocols developed in accordance with the existing renal function assessment guidelines. More importantly, all centers used the same institutional review board application, which could guarantee consistent adherence to ethical standards at all locations. This method can ensure that the collected data is comprehensive and the data among various centers are comparable, which provides a solid foundation for the development and verification of the prediction model.

2.3 Machine learning model construction

The workflow for building the machine learning model was carefully crafted to combine multiple prediction methods with a specific selection of features and parameters. Our method had three components: ensemble model structure, feature selection pipeline, and training process optimization, as shown in Figure 3. This approach aligns with recent advancements in clinical ML frameworks, such as the user-friendly ML pipeline proposed by Orhan et al. (27) for cardiac structure assessment, which emphasizes interpretability and clinical applicability. Similarly, our ensemble framework prioritizes interpretable decision-making to facilitate integration into clinical workflows for chronic kidney disease (CKD) management. With respect to feature selection, we followed the method set forth by Su et al. (26), which used a hybrid model that incorporated both statistical significance and domain knowledge.

Figure 3
Diagram illustrating a machine learning model architecture and pipeline. Section A (Model Architecture) includes Random Forest, XGBoost, and LightGBM models combined into a weighted ensemble for risk prediction. Section B (Feature Selection Pipeline) starts with 120 initial features, reduced through statistical selection and clinical validation to 45 features. Section C (Training Process) details a data split (70% train, 15% validation, 15% test), cross-validation with stratified sampling, hyperparameter optimization using Bayesian methods, and model evaluation with AUC-ROC and PR-Curve.

Figure 3. Model architecture. (A) Ensemble model structure illustrating the integration of Random Forest, XGBoost, and LightGBM. (B) Feature selection pipeline combining filter and wrapper methods with clinical domain knowledge. (C) Training process optimization, including hyperparameter tuning and Bayesian optimization framework.

The three base learners utilized by the foundational ensemble architecture were: Random Forest, XGBoost, and Light GBM. The primary goal hyperbolic function related to model optimization can be formulated as shown in Equation (1):

L ( θ ) = 1 N i = 1 N [ α · BCE ( y i , y i ̂ ) + β · FL ( y i , y i ̂ ) ] + λ θ 2     (1)

Within the formulation, BCE refers to the binary cross-entropy loss, FL describes the focal loss part, and ∥θ2 is the L2 regularization term. As in the case with Bellocchi et al. (28), cross-validation was used to optimize the hyperparameters α , β and λ .The procedure for feature selection was done using both filter and wrapper methods, where the score of importance was computed as shown in Equation (2):

IS j = γ · MI ( X j , Y ) + ( 1 γ ) · SHAP j     (2)

In the equation, MI ( X j , Y ) is the mutual information of the feature X j with regard to the target Y , while SHAP j is the SHAP value contribution of the jth feature. Following Zacharias et al. (25), this feature selection process was iteratively modified guided by clinical domain knowledge.

The hyperparameters for each model were systematically optimized as shown in Table 1.

Table 1
www.frontiersin.org

Table 1. Hyperparameters of different machine learning models.

This was completed as part of a guided capture-the-flag competition, which uses Time 4 Learning’s training resources to prepare. Each model is updated using the same training schema as Ferguson et al. (23) predictions and includes a stratified 5-fold cross-validation scheme. The model ensemble prediction was made based on the weighted average method as shown in Equation (3):

y ̂ = k = 1 K w k · f k ( x )     (3)

where f k ( x ) represents the prediction from the k th base model and w k are the optimized model weights determined through validation performance. The hyperparameter optimization process utilized a Bayesian optimization framework with the expected improvement acquisition function as shown in Equation (4):

EI ( x ) = E [ max ( f ( x ) f ( x + ) , 0 ) ]     (4)

where f ( x + ) represents the current best observed performance. This approach, validated by Miller et al. (29), enabled efficient exploration of the hyperparameter space while balancing exploration and exploitation.

2.4 Model validation method

The model validation process was methodically crafted to enable effective performance evaluation and clinical relevance. Based on the model defined by Churpek et al. (30), we devised a systematic validation plan that utilized both internal and external validation.

To confirm the internal validity, we used a stringent 5-fold cross-validation method. The performance metrics were computed as shown in the following formulations (31) as shown in Equation (5):

AUC = 0 1 TP ( θ ) P ( 1 FP ( θ ) N )     (5)

where TP ( θ ) and FP ( θ ) represent the true positive and false positive rates at threshold θ , respectively. The calibration assessment utilized the Brier score as shown in Equation (6):

BS = 1 N i = 1 N ( y i p ̂ i ) 2     (6)

where p i ̂ represents the predicted probability for the i th instance. The model’s discrimination ability was evaluated using multiple metrics as shown in Table 2.

Table 2
www.frontiersin.org

Table 2. Performance metrics of the ensemble model in internal validation cohort.

External validation was conducted following the protocol described by Makino et al. (32), utilizing an independent cohort from three external medical centers. The concordance between predicted and observed risks was assessed using the calibration slope ( β ) as shown in Equation (7):

logit ( P observed ) = α + β · logit ( P predicted )     (7)

The model’s performance was compared with existing prediction methods through net reclassification improvement (NRI) as shown in Equation (8):

NRI = ( n up , events n events n up , nonevents n nonevents ) ( n down , events n events n down , nonevents n nonevents )     (8)

where n up and n down represent the number of individuals with upward and downward risk reclassification, respectively. As demonstrated by Ekundayo et al. (33), this approach provides a comprehensive assessment of the model’s incremental value.

The integrated discrimination improvement (IDI) was calculated as shown in Equation (9):

IDI = ( P new , events ¯ P new , nonevents ¯ ) ( P old , events ¯ P old , nonevents ¯ )     (9)

where P ¯ represents the mean predicted probabilities. This metric, as validated by Delrue et al. (34), quantifies the model’s improved ability to separate events from non-events.

The comparative analysis results with existing methods are presented in Table 3:

Table 3
www.frontiersin.org

Table 3. Comparative analysis of different risk prediction models.

These comprehensive validation results demonstrate the robust performance and generalizability of our proposed model across different clinical settings and patient populations.

3 Results

3.1 Baseline characteristics of the study population

The refracted demographic and clinical picture of the 1,200 participants was uncovered in the study cohort which revealed the risk underlying the decline in kidney function, as shown in Figure 4 and Table 4. Known characteristics of the study population exhibited significant variations between the progression and non-progression groups throughout multiple dimensions.

Figure 4
Bar charts displaying baseline characteristics by disease progression status. (a) Age Distribution: similar ages for progression and non-progression groups. (b) Comorbidity Prevalence: higher percentages in progression group for hypertension, diabetes, and cardiovascular diseases. (c) eGFR Distribution: comparable eGFR levels for both groups. (d) UPCR Distribution: higher levels in progression group.

Figure 4. Baseline characteristics stratified by disease progression status.

Table 4
www.frontiersin.org

Table 4. Baseline characteristics of study participants.

The mean age of the progression group was significantly higher at 64.8 ± 12.6 years than the non-progression group’s mean age of 61.2 ± 13.8 years ( p = 0.003). This difference in age distribution proved to be statistically significant, as depicted in Figure 4. This finding indicates age may be an influencing factor for kidney function deterioration.

Comorbidity analysis showed that Hypertension had the most pronounced difference, affecting 77.3% of the progression group versus 67.4% of the non-progression group ( p < 0.001). The burden of chronic conditions analysed together proved to be markedly higher in the progression group (45.8%) than the non-progression group (33.6%) in diabetes with a statistical difference ( p < 0.001). CVD followed this trend with a 29.9% prevalence in the progression group compared to 20.7% in the non-progression group ( p = 0.001).

The intricate metabolic signatures distinguishing progression trajectories are shown on Figure 4 and Table 4’s laboratory parameters. The estimated glomerular filtration rate (eGFR) divergence was noteworthy with the lower values of the progression group (42.3 ± 16.2 mL/min/1.73m2), when compared to the non-progression group’s 47.8 ± 14.9 mL/min/1.73m2. This difference illustrates the importance of renal function indicators in predicting the progression of disease.

The urinary protein-to-creatinine ratio (UPCR) provided additional clarity into the already intricate terrain of the decline in kidney function. As depicted in Figure 4, the progression group had a higher average UPCR which corresponds to higher proteinuria and possible renal injury. These biochemical differences offer important information about the mechanisms of kidney function decline.

The analysis of the cohort’s baseline characteristics is comprehensive in scope and illustrates the multifactorial aspect of kidney function decline. The differences were statistically significant and spread across demographic, comorbidity, and laboratory parameters, which adds to the depth of renal disease progression. This nuanced characterization provides not only a complex snapshot of the population, but also an understanding that goes beyond the mechanisms of renal function decline, which is unprecedented for the machine learning model’s predictive architecture.

The graph shows the distribution of a cohort’s baseline characteristics which include age, comorbidity burden, estimated glomerular filtration rate (eGFR), and urinary protein to creatinine ratio (UPCR) in both progression and non-progression groups and their correlates.

3.2 Model performance evaluation

The evaluations conducted on the machine learning model showed predictive power on all the metrics. The ensemble model, as predicted by the receiver operating characteristic (ROC) analysis shown in Figure 5a, was found to have better discrimination ability than the individual base learners. The ensemble model attained an area under the ROC curve (AUC) of 0.89 (95% CI: 0.87–0.91), which was much higher than the isolating cases of random forest (AUC: 0.85, 95% CI: 0.83–0.87) and XGBoost (AUC: 0.87, 95% CI: 0.85–0.89) models, and even outperformed them in recurrent measures.

Figure 5
Three-panel image showing ROC curves, calibration curves, and a confusion matrix. Panel (a) displays ROC curves for Ensemble, Random Forest, LightGBM, and XGBoost models, with the area under the curve (AUC) values of 0.89, 0.85, 0.86, and 0.87, respectively. Panel (b) displays calibration curves for an Ensemble model compared to perfect and typical models. Panel (c) shows a confusion matrix with true labels on the y-axis and predicted labels on the x-axis, featuring 97 true negatives, 3 false positives, 24 false negatives, and 27 true positives.

Figure 5. Comprehensive evaluation of model performance in predicting kidney function decline. (A) Receiver Operating Characteristic (ROC) curve analysis comparing the ensemble model to individual base learners. (B) Calibration plot demonstrating the alignment between predicted and observed risks. (C) Confusion matrix illustrating classification performance with true positives, true negatives, false positives, and false negatives.

These block figures include, but are not limited to, the performance metrics of the single models in comparison to the ensemble model for their different instances at various datasets as mentioned in Table 5.

Table 5
www.frontiersin.org

Table 5. Comprehensive performance metrics across different datasets.

The machine learning ensemble model offers a transformative tool for predicting renal function decline in chronic kidney disease (CKD), providing clinicians with reliable and actionable insights for personalized care. The calibration analysis (Figure 5b) demonstrates the model’s exceptional reliability, with predicted risks closely mirroring actual outcomes across the entire risk spectrum. With a calibration slope of 0.96 (95% CI: 0.94–0.98) and an intercept of 0.02 (95% CI: 0.01–0.03), the model exhibits minimal bias, ensuring that clinicians can confidently use its risk estimates to guide treatment decisions. This robust calibration means that a predicted 30% risk of CKD progression accurately reflects the true likelihood, enabling precise patient counseling and intervention planning. The confusion matrix in Figure 5c demonstrates the classification performance of the model in predicting the risk of renal function decline, reflecting its performance on true positives (TP), true negatives (TN), false positives (FP), and false negatives (FN). The confusion matrix matrix shows that the model has similar accuracy in predicting true positives and true negatives, indicating its balanced performance in distinguishing cases of renal function decline from non-decline cases.

The result reflects their predictive power of the ensemble model’s reliability and performance regarding decline in kidney function. All of its aspects, including calibration, discrimination, subgroup performance, and validation, re-confirm the model’s effectiveness in integrating early risk assessment and intervention within clinical practice.

3.3 Analysis of the results

As noted previously with the kidney pathology overview, the deep dive into the kidney function decline risk analysis illustrated the intricate interrelationships of several clinical elements and their predictive effects. Feature importance and their multitude of permutations is illustrated in Figure 6.

Figure 6
Feature importance and interaction analysis with four panels: (a) A bar chart showing SHAP feature importance, highlighting eGFR, Age, and UPCR as significant. (b) Box plots displaying feature distribution for eGFR, Age, UPCR, and Creatinine. (c) A line graph illustrating eGFR SHAP dependence, showing an increasing trend. (d) A heatmap of feature interactions indicating strong relationships, with higher values in eGFR and Age.

Figure 6. Feature importance and interaction analysis. (a) SHAP feature importance; (b) feature distribution; (c) eGFR SHAP dependence; (d) feature interaction.

The SHAP importance deconstruction revealed an obvious importance structure related to predictive factors. Estimated glomerular filtration rate (eGFR) was by far the most pivotal predictor, as expected from the magnitude of the SHAP value, followed by age, and urinary protein creatinine ratio (UPCR). The above highlighted aspects reinforce the complexity involved in the kidney function decline risk assessment, which is profoundly multifactorial. Table 6 captures the overview of the most key risk factors with their clinical importance in detail.

Table 6
www.frontiersin.org

Table 6. Key risk factors and clinical significance.

As seen in Figure 6b, there was a feature distribution boxplot that showed the differences which existed among some clinical parameters. The distributions for eGFR and age displayed greater variation which indicates how multi-faceted and variable these parameters are within the scope of kidney function evaluation. A non-linear relationship was demonstrated in eGFR’s SHAP dependence plot in Figure 6c, which underlined the decline in kidney function’s intricate mechanisms.

In Figure 6d, feature interaction analysis showed important dependencies of some clinical markers. The interaction heatmap showed strong, and even moderate, differences especially with eGFR, age and UPCR. Such relations indicate that the decline in kidney function is not the result of a singular issue, rather, it is a product of many interacting physiological parameters.

In particular, Figure 7 highlights a complete interpretation framework for clinical risk. The compositional analysis of risk factors contribution waterfall plot in Figure 7a showed risk contributions were cumulative where baseline characteristics and central clinical features adjusted the risk over time. The prediction probability distribution in Figure 7b was able to distinctly classify patients into three groups: low, medium, and high risk, which was very useful for personalized risk evaluation. Figure 7c shows the scatter plot of risk features against predicted risk wherein the correlation was highly positive with the multicolored risk indicators representing the constructs of interest. Risk stratification within subgroups in Figure 7d demonstrated that there was heterogeneity among the different patient populations, notably higher risk probabilities for elderly patients and those with multiple comorbidities.

Figure 7
Clinical risk interpretation and stratification visualization includes: (a) bar and line graphs showing risk factors contributions like eGFR and age; (b) box plots displaying prediction probability distribution across low, medium, and high risk categories; (c) scatter plot with regression line illustrating features versus prediction results related to eGFR; (d) error bar chart depicting risk stratification by subgroups such as young, middle-aged, elderly, and comorbid.

Figure 7. Clinical risk interpretation and stratification. (a) Risk factors contribution; (b) prediction probability distribution; (c) features vs. prediction results; (d) risk stratification by subgroups.

Offers quantifiable metrics alongside intricate biological explanations for a process that has remained largely qualitative. Along with providing an innovative means of risk identification and intervention, this sophisticated analysis brings a new dimension for understanding the decline of kidney functions due to advanced age. The integration of publicly available healthcare datasets along with augmented machine learning enables doctors to implement shifts in clinical paradigms more quickly than before.

Pioneers a new era in computing and healthcare integration by offering precise measures to counteract the deterioration of kidney functions. This will provide room for further innovation that challenges existing practices in nephrology.

3.4 Comparison with traditional methods

In comparison with conventional approaches, the analysis carried out between our proposed Machine Learning model and other techniques showed great improvements in predictability and clinical usefulness. Figure 8 clearly shows that the residual plot displays a normal distribution of errors centered around 0. The traditional method had an AUC of 0.695 from the ROC curve analysis, and integration with the calibration plots showed exceptionally good agreement between predicted and observed probabilities across the entire risk spectrum.

Figure 8
Panel (a) shows a bar chart of prediction errors with frequency on the y-axis and error ranging from negative twenty to twenty on the x-axis. Panel (b) features a Receiver Operating Characteristic (ROC) curve with an AUC of zero point six eight five, plotting true positive rate against false positive rate. Panel (c) displays a calibration plot with observed probability versus predicted probability, indicating alignment along the diagonal.

Figure 8. Comprehensive model performance analysis. (a) Error distribution histogram; (b) ROC curve analysis (AUC = 0.695); (c) calibration plot with confidence intervals.

The benefits of advanced clinical applications are depicted thoroughly in Figure 9, as multi-layered performance analytics outlines how much more performant our suggested ML model is relative to both Cox and standard ML models. The accuracy evaluation by strata reveals performance consistency across different patient subgroups, as well as showing enhanced ability to predict the passage of time in regard to disease progression. The cost-effectiveness analysis also confirms the projected practical benefits for our approach from the standpoint of actual clinical use.

Figure 9
Four graphs showing clinical application analyses. (a) Bar graph comparing performance scores for Prediction Time, Computational Complexity, Memory Usage, and Scalability among Proposed ML, Cox Model, and Standard ML. (b) Box plots illustrating accuracy across age groups: Young, Middle-aged, Elderly, and Comorbid. (c) Line graph showing Disease Progression Dynamics, with patient count and cumulative risk over time. (d) Bar graph depicting Cost-Effectiveness Analysis for Direct Costs, System Integration, Implementation Cost, and Risk Reduction among Proposed ML, Cox Model, and Standard ML.

Figure 9. Advanced clinical application advantages. (a) Multi-dimensional performance analysis comparing proposed ML, Cox model, and standard ML; (b) stratified accuracy analysis across patient subgroups; (c) disease progression dynamics with time-to-progression analysis; (d) cost-effectiveness analysis across different implementation aspects.

Our Proposed Solutions: Enhanced Capabilities The detailed performance metrics pertaining to diverse methodological approaches have been presented in Table 7. Table 8 presents the detailed training methods and model parameters of the baseline models, Transformer and RNN.

Table 7
www.frontiersin.org

Table 7. Comparative performance analysis of different prediction models for kidney function decline.

Table 8
www.frontiersin.org

Table 8. Model architectures and training details.

The formula for calculating Resource Utilization is shown in Equation 10, Resource Usage proposed represents the cost of training and inference for the proposed machine learning model (an ensemble model based on Random Forest, XGBoost, and LightGBM) on the cloud service platform. Resource Usage represents the cost of training and reasoning a traditional Cox proportional hazards regression model on a cloud service platform under the same conditions.

Scalability Index is used to measure the performance stability of a model in different dataset sizes or clinical scenarios. The specific calculation formula is shown in Equation 11. Among them, a reflects the degree of fluctuation of the AUC index of the model in different scale datasets, measuring its predictive stability in different datasets or scenarios. max Variance represents the maximum value of performance variance. This offers quantitative proof supporting the improved features of our model.

Cost norm = Resource Usage proposed Resource Usage     (10)
Scalability Index = 1 Performance Variance max Variance     (11)

The AUC of our proposed ensemble model reached 0.879 (95% CI: 0.856–0.902), significantly outperforming the traditional Cox model’s AUC of 0.695 (95% CI: 0.672–0.718), standard ML’s AUC of 0.782 (95% CI: 0.759–0.805), the Transformer model’s AUC of 0.870 (95% CI: 0.847–0.893), and the RNN model’s AUC of 0.810 (95% CI: 0.787–0.833). In terms of computation time, our model achieved a prediction time of 0.48 ± 0.05 s, a 74.2% improvement over the Cox model’s 1.86 ± 0.12 s, and was notably faster than the Transformer (3.35 ± 0.10 s) and RNN (2.95 ± 0.07 s) models, which were less efficient than even the standard ML model (0.92 ± 0.08 s). Resource utilization was optimized by 60.6% compared to the Cox model (28.5 ± 3.2% vs. 72.3 ± 5.1%), with our model also outperforming the standard ML (45.7 ± 4.3%), Transformer (65.0 ± 4.5%), and RNN (59.0 ± 4.0%) models. The calibration slope of 0.96 (95% CI: 0.94–0.98) underscored the model’s reliability, with minimal discrepancy between predicted and observed risks, confirming its excellent performance in risk stratification for kidney function decline.

Cost-effectiveness analysis revealed a 35% reduction in implementation cost for our proposed ensemble model (0.65 ± 0.07) compared to the traditional Cox model (1.00 ± 0.00), outperforming the standard ML model (0.82 ± 0.05), Transformer (0.92 ± 0.06), and RNN (0.87 ± 0.05) models, while maintaining superior predictive accuracy (AUC: 0.879). The scalability index of 0.92 ± 0.03 demonstrated robust performance across varying dataset sizes, significantly surpassing the Cox model (0.45 ± 0.05), standard ML (0.67 ± 0.04), Transformer (0.83 ± 0.04), and RNN (0.72 ± 0.04) models. These results, supported by rigorous internal and external validation (Tables 2, 3, 7), highlight the model’s efficiency and generalizability, positioning it as a highly viable tool for clinical integration across diverse settings to predict kidney function decline risk effectively.

3.5 Clinical case analysis

To enhance the clinical utility of the model, we provided interpretability through SHAP analysis and further demonstrated its application in clinical decision-making through case snippets and integration strategies with electronic health records (EHR). The following case illustrates how the model prediction can guide personalized management of patients with chronic kidney disease (CKD).

4 Discussion

Our analysis reveals important implications for clinical practice and offers some insights into the splendid capability of machine learning techniques in predicting the decline of kidney functions. The efficacy of our ensemble model, which achieved an AUC of 0.89 (95% CI: 0.87–0.91), confirms that the integration of numerous machine learning algorithms for intricate clinical forecasts is effective (35). The performance of this ensemble model significantly exceeds the AUC of conventional statistical approaches and outliers in progression prognosis for chronic kidney disease, signifying an advancement in stratification competence. Previous studies have reported the optimization of risk stratification due to the incorporation of electronic health records with machine learning algorithms (36). Our findings further confirm this approach through extensive validation across numerous clinical settings (Table 9).

Table 9
www.frontiersin.org

Table 9. Analysis table of clinical cases of different patients.

The application of ensemble frameworks to provide the merging of several algorithms is one of the changes we made to the machine learning application in clinical prediction. This approach is one of the numerous solutions to the many challenges faced in healthcare predictive modeling (37, 38). The model’s outstanding calibration (slope: 0.96, 95% CI: 0.94–0.98) illustrates a considerable leap in addressing the remaining issues of implementing machine learning in healthcare (39). The reliability of artificial intelligence in predicting the worsening of kidney diseases is known to be high (40). Our results offer substantial proof toward the adoption of these findings into clinical work.

In this research, there has been remarkable progress, but there are still some areas that require further attention. First, the adoption of deep learning techniques, as well as the threat of data leakage (41, 42), both warrant further exploration. More efforts need to be directed at potentially overfitting the models in immunology (43) and at the same time increasing the scope of the model to include new biomarker and genetic influences. For instance, Çiçek et al. (44) demonstrated that preoperative neopterin levels can predict acute kidney injury in on-pump cardiac surgery, highlighting the critical role of biomarker-driven risk stratification in kidney outcomes. This supports our proposition to incorporate novel biomarkers, such as neopterin or other inflammatory markers, to enhance the phenomenological capabilities of our model for CKD progression. The development of artificial intelligence in medicine (45) presents new possibilities for the inclusion of other features such as genomic and proteomic markers that would improve the model’s phenomenological capabilities. This study suggests the usage of automated methods for model updating, uniform data gathering from clinics, and the creation of clear multi-center validation procedures as the focus of future work. The addition of real-time clinical decision support systems and the extension of the model functionality to new emerging biomarkers is the next crucial step in the progression of this area.

Machine learning offers extraordinary promise for transforming the prediction of kidney function decline, which is why a great many obstacles still need to be solved before we can implement our research in a clinical setting. Our research substantiates machine learning and kidney pathology by laying out the groundwork for personalized medicine and data-centric healthcare decision-making in nephrology. As further changes in the healthcare system occur, our model will be more useful in enhancing the quality of care provided and in the efficient use of resources for chronic kidney disease treatment and prognosis.

4.1 Model stability analysis results

Table 10 shows the stability performance of the integrated model in predicting the risk of renal function decline. The AUC stability of the model is 0.87 ± 0.02, with a coefficient of variation (CV) of only 2.3%, indicating that its predictive performance is highly consistent across multiple runs. Sensitivity (0.86 ± 0.03, CV = 3.5%) and specificity (0.84 ± 0.03, CV = 3.6%) also showed low volatility, demonstrating the robustness of the model on different datasets. The calibration slope (0.96, 95% CI: 0.94–0.98, CV = 2.1%) and intercept (0.02, 95% CI: 0.01–0.03, CV = 1.8%) further confirmed the high consistency between the model predictions and actual results. These results indicate that the model can maintain reliable predictive performance in different operational and clinical scenarios, and is suitable for a wide range of clinical applications.

Table 10
www.frontiersin.org

Table 10. Model stability analysis results.

4.2 Comprehensive evaluation of model performance

The decision curve analysis (DCA, Figure 5d) highlights the model’s practical utility in clinical settings. It shows a substantial net benefit over default strategies of treating all or no patients, particularly in the 20–60% risk range, where clinical decisions are most critical. For example, in this range, the model helps clinicians identify patients who would benefit most from intensified monitoring or early interventions, such as medication adjustments, while sparing low-risk patients unnecessary treatments. This targeted approach optimizes resource use and enhances patient outcomes by focusing efforts where they are most needed.

Subgroup analyses (Figure 5e) further underscore the model’s versatility across diverse patient populations, with outstanding performance in high-risk groups. For elderly patients, the model achieves an AUC of 0.88 (95% CI: 0.85–0.91), and for those with diabetes, it reaches an AUC of 0.90 (95% CI: 0.87–0.93). These groups are particularly vulnerable to rapid CKD progression, and the model’s high accuracy in predicting their risk enables earlier and more tailored interventions, such as stricter blood pressure control or diabetes management, to slow disease progression. By providing clear, interpretable risk stratification, the model empowers clinicians to make data-driven decisions that improve patient care and quality of life.

4.3 Sensitivity analysis of model performance across renal function decline definitions

Table 11 presents the sensitivity analysis of the model’s performance across various definitions of renal function decline, demonstrating its robustness. The model achieves a high AUC of 0.89 (95% CI: 0.87–0.91) for the primary definition (eGFR decline ≥30% or dialysis), with strong sensitivity (0.86) and specificity (0.82). Alternative definitions, such as eGFR decline ≥20%, ≥40%, serum creatinine doubling, and progression to dialysis, yield slightly lower but still robust AUCs (0.86–0.88), with sensitivity and specificity ranging from 0.80–0.85 and 0.77–0.83, respectively. Calibration slopes remain excellent (0.93–0.96), indicating consistent alignment between predicted and observed risks. These results confirm the model’s stable performance across diverse clinical definitions, enhancing its reliability and applicability for risk stratification in chronic kidney disease management (Figure 10).

Table 11
www.frontiersin.org

Table 11. Model performance under different definitions of renal function decline.

Figure 10
Three-panel image showing: (d) Decision Curve Analysis with curves for Ensemble Model, Random Forest, Treat All, and Treat None, plotting Net Benefit vs. Threshold Probability. (e) Subgroup Analysis with Hazard Ratios and 95% CI for different groups. (f) Stability Analysis showing AUC variation across Iterations, highlighting the mean AUC.

Figure 10. Decision curve, subgroup performance and stability analysis.

5 Conclusion

In this self-contained piece of research, we outline the design and validation of an automated machine learning model for predicting the risk of decline in kidney function, which outperformed the conventional methods. Our ensemble model achieved astounding accuracy (AUC: 0.89, 95% CI: 0.87–0.91) in prediction of events, while the calibration of the model remained impressive in diverse populations. The use of several techniques in one novel ensemble framework accompanied by advanced feature selection has provided a solid base for clinical risk prediction in nephrology.

The model clarifies the importance of predictive factors, notably ascribing most eGFR, age, and urinary protein to creatinine ratio, which makes understanding the precise mechanisms of kidney function deterioration easier. Improved understanding, along with the model’s predictive performance, enhances the capability of healthcare practitioners to undertake early risk stratification and tailor interventions in a precise manner. The accuracy demonstrated among various patient subgroups and validation cohorts confirms the model’s potential value for widespread clinical use.

The results of this study are particularly relevant to the clinical management and future directions of research in nephrology. If this predictive tool is successfully adopted into clinical workflows, it has the potential to revolutionize chronic kidney disease management by allowing for timely and precise interventions and resource assignment. As data-centric decision-making continues to gain traction in healthcare systems, our model serves a robust and practical purpose for predicting the risk of kidney function decline, with the possibility to improve patient care by targeting interventions sooner and more effectively. Future studies need to concentrate on multicenter validation studies and how the model’s prediction and clinical application may be augmented through the use of novel biomarkers.

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

This study was conducted in accordance with the guidelines of the Declaration of Helsinki and was approved by the Ethics Committee of the 95th Hospital of Putian in China RongTong Medical Health Corporation. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

HC: Writing – original draft, Writing – review & editing. YH: Writing – original draft. LC: Writing – original draft.

Funding

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

Conflict of interest

The 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 author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1598065/full#supplementary-material

References

1. Bikbov, B, Purcell, CA, Levey, AS, and Smith, M. Global, regional, and national burden of chronic kidney disease, 1990–2017: a systematic analysis for the global burden of disease study 2017. Lancet. (2020) 395:709–33. doi: 10.1016/S0140-6736(20)30045-3

PubMed Abstract | Crossref Full Text | Google Scholar

2. Rajpurkar, P, Chen, E, and Banerjee, O. AI in health and medicine. Nat Med. (2022) 28:31–8.

Google Scholar

3. Saleh, A, Shaban, NG, and Al-Asklany, HM. Value of serum kidney injury Molecule-1 in early prediction of kidney injury in patient with ascites and spontaneous bacterial peritonitis. Egypt J Hosp Med. (2023) 92:5487–95. doi: 10.21608/ejhm.2023.305789

Crossref Full Text | Google Scholar

4. Li, C, Liu, J, Fu, P, and Zou, J. Artificial intelligence models in diagnosis and treatment of kidney diseases: current status and prospects. Kidney Dis. (2025) 11:491–507. doi: 10.1159/000546397

PubMed Abstract | Crossref Full Text | Google Scholar

5. Gulamali, FF, Sawant, AS, and Nadkarni, GN. Machine learning for risk stratification in kidney disease. Curr Opin Nephrol Hypertens. (2022) 31:548–52. doi: 10.1097/MNH.0000000000000832

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kanda, E, Epureanu, BI, Adachi, T, and Kashihara, N. Machine-learning-based web system for the prediction of chronic kidney disease progression and mortality. PLoS Digit Health. (2023) 2:e0000188. doi: 10.1371/journal.pdig.0000188

PubMed Abstract | Crossref Full Text | Google Scholar

7. Islam, MA, Majumder, MZH, and Hussein, MA. Chronic kidney disease prediction based on machine learning algorithms. J Pathol Inform. (2023) 14:100189. doi: 10.1016/j.jpi.2023.100189

PubMed Abstract | Crossref Full Text | Google Scholar

8. Beissel, PA. Equity Valuation Fresenius Medical Care AG & co. KGaA [D]. Portugal: Universidade Catolica Portuguesa (2023).

Google Scholar

9. Vagliano, I, Chesnaye, NC, Leopold, JH, Jager, KJ, Abu-Hanna, A, and Schut, MC. Machine learning models for predicting acute kidney injury: a systematic review and critical appraisal. Clin Kidney J. (2022) 15:2266–80. doi: 10.1093/ckj/sfac181

PubMed Abstract | Crossref Full Text | Google Scholar

10. Sanmarchi, F, Fanconi, C, Golinelli, D, Gori, D, Hernandez-Boussard, T, and Capodici, A. Predict, diagnose, and treat chronic kidney disease with machine learning: a systematic literature review. J Nephrol. (2023) 36:1101–17. doi: 10.1007/s40620-023-01573-4

PubMed Abstract | Crossref Full Text | Google Scholar

11. Chicco, D, Warrens, MJ, and Jurman, G. The coefficient of determination R-squared is more informative than SMAPE, MAE, MAPE, MSE and RMSE in regression analysis evaluation. Peerj Computer Sci. (2021) 7:e623. doi: 10.7717/peerj-cs.623

PubMed Abstract | Crossref Full Text | Google Scholar

12. Nelson, RG, Grams, ME, Ballew, SH, Sang, Y, Azizi, F, Chadban, SJ, et al. Development of risk prediction equations for incident chronic kidney disease. JAMA. (2019) 322:2104–14. doi: 10.1001/jama.2019.17379

PubMed Abstract | Crossref Full Text | Google Scholar

13. Tangri, N, Stevens, LA, Griffith, J, Tighiouart, H, Djurdjev, O, Naimark, D, et al. A predictive model for progression of chronic kidney disease to kidney failure. JAMA. (2011) 305:1553–9. doi: 10.1001/jama.2011.451

PubMed Abstract | Crossref Full Text | Google Scholar

14. Prasad, ML, Kiran, A, and Shaker Reddy, PC. Chronic kidney disease risk prediction using machine learning techniques. J Inf Technol Manag. (2024) 16:118–34. doi: 10.22059/jitm.2024.96378

Crossref Full Text | Google Scholar

15. Islam, M. A., Akter, S., and Hossen, M. S. “Risk factor prediction of chronic kidney disease based on machine learning algorithms”; Proceedings of the 2020 3rd international conference on intelligent sustainable systems (ICISS), F, (2020). IEEE.

Google Scholar

16. Kalupukuru, SR, and Natarajan, K. Machine learning methods for predicting the prognosis of chronic kidney disease. Procedia Comput Sci. (2025) 258:1372–82. doi: 10.1016/j.procs.2025.04.370

Crossref Full Text | Google Scholar

17. Ghosh, SK, and Khandoker, AH. Investigation on explainable machine learning models to predict chronic kidney diseases. Sci Rep. (2024) 14:3687. doi: 10.1038/s41598-024-54375-4

PubMed Abstract | Crossref Full Text | Google Scholar

18. Bai, Q, Su, C, Tang, W, and Li, Y. Machine learning to predict end stage kidney disease in chronic kidney disease. Sci Rep. (2022) 12:8377. doi: 10.1038/s41598-022-12316-z

PubMed Abstract | Crossref Full Text | Google Scholar

19. Lei, N, Zhang, X, Wei, M, Lao, B, Xu, X, Zhang, M, et al. Machine learning algorithms’ accuracy in predicting kidney disease progression: a systematic review and meta-analysis. BMC Med Inform Decis Mak. (2022) 22:205. doi: 10.1186/s12911-022-01951-1

PubMed Abstract | Crossref Full Text | Google Scholar

20. Figueroa, RL, Zeng-Treitler, Q, and Kandula, S. Predicting sample size required for classification performance. BMC Med Inform Decis Mak. (2012) 12:8.

Google Scholar

21. Chan, L, Nadkarni, GN, Fleming, F, McCullough, JR, Connolly, P, Mosoyan, G, et al. Derivation and validation of a machine learning risk score using biomarker and electronic patient data to predict progression of diabetic kidney disease. Diabetologia. (2021) 64:1504–15. doi: 10.1007/s00125-021-05444-0

PubMed Abstract | Crossref Full Text | Google Scholar

22. Xiao, J, Ding, R, and Xu, X. Comparison and development of machine learning tools in the prediction of chronic kidney disease progression. J Transl Med. (2019) 17:119.

Google Scholar

23. Ferguson, T, Ravani, P, Sood, MM, Clarke, A, Komenda, P, Rigatto, C, et al. Development and external validation of a machine learning model for progression of CKD. Kidney Int Rep. (2022) 7:1772–81. doi: 10.1016/j.ekir.2022.05.004

PubMed Abstract | Crossref Full Text | Google Scholar

24. Iftikhar, H, Khan, M, Khan, Z, Khan, F, Alshanbari, HM, and Ahmad, Z. A comparative analysis of machine learning models: a case study in predicting chronic kidney disease. Sustainability. (2023) 15:2754. doi: 10.3390/su15032754

Crossref Full Text | Google Scholar

25. Zacharias, HU, Altenbuchinger, M, and Schultheiss, UT. A predictive model for progression of CKD to kidney failure based on routine laboratory tests. Am J Kidney Dis. (2022) 79:e211:217–30.

Google Scholar

26. Su, C-T, Chang, Y-P, Ku, Y-T, and Lin, C-M. Machine learning models for the prediction of renal failure in chronic kidney disease: a retrospective cohort study. Diagnostics. (2022) 12:2454. doi: 10.3390/diagnostics12102454

PubMed Abstract | Crossref Full Text | Google Scholar

27. Orhan, A, Akbayrak, H, Çiçek, ÖF, Harmankaya, İ, and Vatansev, H. A user-friendly machine learning approach for cardiac structures assessment. Front Cardiovasc Med. (2024) 11:1426888. doi: 10.3389/fcvm.2024.1426888

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bellocchio, F, Lonati, C, and Ion, J. Validation of a novel predictive algorithm for kidney failure in patients suffering from chronic kidney disease: The Prognostic Reasoning System for Chronic Kidney Disease (PROGRES-CKD). Int J Environ Res Public Health. (2021) 18:12649. doi: 10.3390/ijerph182312649

PubMed Abstract | Crossref Full Text | Google Scholar

29. Vastrad, B, and Vastrad, C. Identification of candidate biomarkers and signaling pathways associated with Alzheimer’s disease using bioinformatics analysis of next generation sequencing data. bioRxiv. (2012) 2024:626535.

Google Scholar

30. Churpek, MM, Carey, KA, Edelson, DP, Singh, T, Astor, BC, Gilbert, ER, et al. Internal and external validation of a machine learning risk score for acute kidney injury. JAMA Netw Open. (2020) 3:e2012892. doi: 10.1001/jamanetworkopen.2020.12892

PubMed Abstract | Crossref Full Text | Google Scholar

31. Segal, Z, Kalifa, D, Radinsky, K, Ehrenberg, B, Elad, G, Maor, G, et al. Machine learning algorithm for early detection of end-stage renal disease. BMC Nephrol. (2020) 21:518. doi: 10.1186/s12882-020-02093-0

PubMed Abstract | Crossref Full Text | Google Scholar

32. Makino, M, Yoshimoto, R, Ono, M, Itoko, T, Katsuki, T, Koseki, A, et al. Artificial intelligence predicts the progression of diabetic kidney disease using big data machine learning. Sci Rep. (2019) 9:11862. doi: 10.1038/s41598-019-48263-5

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ekundayo, F. Machine learning for chronic kidney disease progression modelling: leveraging data science to optimize patient management. World J Adv Res Rev. (2024) 24:453–75.

Google Scholar

34. Delrue, C, De Bruyne, S, and Speeckaert, MM. Application of machine learning in chronic kidney disease: current status and future prospects. Biomedicine. (2024) 12:568. doi: 10.3390/biomedicines12030568

PubMed Abstract | Crossref Full Text | Google Scholar

35. Gogoi, P, and Valan, JA. Machine learning approaches for predicting and diagnosing chronic kidney disease: current trends, challenges, solutions, and future directions. Int Urol Nephrol. (2025) 57:1245–68. doi: 10.1007/s11255-024-04281-5

PubMed Abstract | Crossref Full Text | Google Scholar

36. Lee, YC, Cha, J, Shim, I, Park, WY, Kang, SW, Lim, DH, et al. Multimodal deep learning of fundus abnormalities and traditional risk factors for cardiovascular risk prediction. NPJ Digit Med. (2023) 6:14. doi: 10.1038/s41746-023-00748-4

PubMed Abstract | Crossref Full Text | Google Scholar

37. He, Y, Shen, Z, and Cui, P. Towards non-iid image classification: a dataset and baselines. Pattern Recogn. (2021) 110:107383. doi: 10.1016/j.patcog.2020.107383

Crossref Full Text | Google Scholar

38. Kernbach, JM, and Staartjes, VE. Foundations of machine learning-based clinical prediction modeling: part II—generalization and overfitting. Machine Learn Clin Neurosci: Foundations Applications. (2021):15–21.

Google Scholar

39. Krittanawong, C. Artificial intelligence in clinical practice: How AI technologies impact medical research and clinics. San Diego: Elsevier (2023).

Google Scholar

40. Ojo, B, and Campbell, CH. Perioperative acute kidney injury: impact and recent update. Curr Opin Anaesthesiol. (2022) 35:215–23. doi: 10.1097/ACO.0000000000001104

PubMed Abstract | Crossref Full Text | Google Scholar

41. Lee, H-T, Cheon, H-R, Lee, S-H, Shim, M, and Hwang, H-J. Risk of data leakage in estimating the diagnostic performance of a deep-learning-based computer-aided system for psychiatric disorders. Sci Rep. (2023) 13:16633. doi: 10.1038/s41598-023-43542-8

PubMed Abstract | Crossref Full Text | Google Scholar

42. Gygi, JP, Kleinstein, SH, and Guan, L. Predictive overfitting in immunological applications: pitfalls and solutions. Hum Vaccin Immunother. (2023) 19:2251830. doi: 10.1080/21645515.2023.2251830

PubMed Abstract | Crossref Full Text | Google Scholar

43. Liu, P-r, Lu, L, Zhang, J-y, Huo, TT, Liu, SX, and Ye, ZW. Application of artificial intelligence in medicine: an overview. Current. Med Sci. (2021) 41:1105–15. doi: 10.1007/s11596-021-2474-3

PubMed Abstract | Crossref Full Text | Google Scholar

44. Çiçek, ÖF, Akyürek, F, Akbayrak, H, Orhan, A, Kaya, EC, and Büyükateş, M. Can preoperative neopterin levels predict acute kidney injury in patients undergoing on-pump cardiac surgery? Turk J Biochem. (2023) 48:531–40. doi: 10.1515/tjb-2023-0074

Crossref Full Text | Google Scholar

45. Okita, J, Nakata, T, Uchida, H, Kudo, A, Fukuda, A, Ueno, T, et al. Development and validation of a machine learning model to predict time to renal replacement therapy in patients with chronic kidney disease. BMC Nephrol. (2024) 25:101. doi: 10.1186/s12882-024-03527-9

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: machine learning, chronic kidney disease progression, risk prediction modeling, clinical decision support, precision nephrology

Citation: Chen H, Huang Y and Chen L (2025) Ensemble machine learning for predicting renal function decline in chronic kidney disease: development and external validation. Front. Med. 12:1598065. doi: 10.3389/fmed.2025.1598065

Received: 22 March 2025; Accepted: 25 September 2025;
Published: 27 October 2025.

Edited by:

Olaniyi Samuel Iyiola, Morgan State University, United States

Reviewed by:

Ömer Faruk Çiçek, Selcuk University, Türkiye
Samit Kumar Ghosh, Khalifa University, United Arab Emirates

Copyright © 2025 Chen, Huang and Chen. 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: Hong Chen, eXlrdzc2QDE2My5jb20=

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.