Skip to main content

Using a cohort study of diabetes and peripheral artery disease to compare logistic regression and machine learning via random forest modeling



This study illustrates the use of logistic regression and machine learning methods, specifically random forest models, in health services research by analyzing outcomes for a cohort of patients with concomitant peripheral artery disease and diabetes mellitus.


Cohort study using fee-for-service Medicare beneficiaries in 2015 who were newly diagnosed with peripheral artery disease and diabetes mellitus. Exposure variables include whether patients received preventive measures in the 6 months following their index date: HbA1c test, foot exam, or vascular imaging study. Outcomes include any reintervention, lower extremity amputation, and death. We fit both logistic regression models as well as random forest models.


There were 88,898 fee-for-service Medicare beneficiaries diagnosed with peripheral artery disease and diabetes mellitus in our cohort. The rate of preventative treatments in the first six months following diagnosis were 52% (n = 45,971) with foot exams, 43% (n = 38,393) had vascular imaging, and 50% (n = 44,181) had an HbA1c test. The directionality of the influence for all covariates considered matched those results found with the random forest and logistic regression models. The most predictive covariate in each approach differs as determined by the t-statistics from logistic regression and variable importance (VI) in the random forest model. For amputation we see age 85 + (t = 53.17) urban-residing (VI = 83.42), and for death (t = 65.84, VI = 88.76) and reintervention (t = 34.40, VI = 81.22) both models indicate age is most predictive.


The use of random forest models to analyze data and provide predictions for patients holds great potential in identifying modifiable patient-level and health-system factors and cohorts for increased surveillance and intervention to improve outcomes for patients. Random forests are incredibly high performing models with difficult interpretation most ideally suited for times when accurate prediction is most desirable and can be used in tandem with more common approaches to provide a more thorough analysis of observational data.

Peer Review reports


Machine learning algorithms, such as random forest (RF) models, are an ensemble regression tree method commonly used for prediction and measuring variable importance in predicting an outcome of interest [1]. Ensemble learning methods use multiple learning algorithms, in this case the RF consists of multiple regression trees, to obtain better predictive performance than an individual algorithm. Machine learning algorithms have long been applied in computational research including data mining, artificial intelligence, and genomics [2,3,4,5,6]. In recent years, machine learning has been championed for use in health care research [7,8,9]. However, due to the lack of detailed understanding of machine learning algorithms among many healthcare workers, this method has been warily utilized among physician investigators where practical application and interpretability for clinical decision-making are essential.

Logistic regression versus machine learning

The application of traditional regression approaches such as logistic regression as compared to machine learning approaches has been widely studied with no unqualified recommendation of which approach is best suited to analyzing data and making predictions [10, 11]. While the interpretability of logistic regression is considered an important feature, the predictive power of RF algorithms sway some in favor of its application. Multivariable regression approaches are frequently employed to produce prediction models. There is a growing consensus that machine learning and the advent of big data in health care offer new opportunities for understanding predictors of disease and treatment outcomes in medicine. Recent literature offers several new studies highlighting the application of RF in health services research, such as identifying multi-drug resistant tuberculosis (TB) versus drug-sensitive TB, predicting SARS COVID-19 infection, and predicting patients most likely to require dental implants. However, no studies to date have applied machine learning methods to examine predictors associated with outcomes for patients with peripheral arterial disease (PAD).

To fully understand the data at hand and how the end results of the analysis will be utilized, both approaches have their role in the comprehensive analysis of the data, particularly for prediction and understanding previously latent factors influencing outcomes. In this study we illustrate how the application of traditional generalized linear regression models and RF models can elucidate the role of sociodemographic and clinical covariates in predicting outcomes for patients with peripheral arterial disease and diabetes mellitus. This study is significant from past studies in that it leverages novel real world data sources, machine learning methods and multivariate logistic regression to provide complementary information for physicians caring for patients with PAD and diabetes.

Using research on peripheral artery disease to illustrate identification of risk factors

Peripheral artery disease (PAD) results from a partial or complete buildup of atherosclerotic plaque in arteries carrying blood to the brain, limbs, and/or organs [12]. PAD occurs at the highest rates in patients with risk factors such as diabetes mellitus (DM) [12, 13]. where patients with both are at higher risk of amputation than those with either condition alone [14]. Furthermore, disparities in amputation rates have been observed among patients with concomitant lower extremity PAD and diabetes across race and geographic region [15,16,17]. Integrated management strategies for diabetes and PAD consist of preventative treatments such as hemoglobin A1c (HbA1c) testing, diabetic foot care, and vascular assessment with imaging. Several studies have shown that these interventions reduce amputation rates in patients with PAD and DM [18,19,20,21,22,23,24,25]. Effective implementation of the recommended integrated management strategies requires identifying the observable and modifiable risk factors which are most predictive of outcomes such as amputation and death in this population for targeted delivery.

Goal of this manuscript

In this methodologic review, we use both traditional generalized linear regression models and practical approach to RF modeling to analyze the data and provide information on which sociodemographic and clinical covariates are most predictive of negative outcomes for patients with PAD and DM. By utilizing this approach we hope to provide insight into what observable factors influence the outcomes in this cohort of patients with PAD and DM. More importantly, we hope to provide a framework which outlines a more comprehensive analytic approach to applying traditional regression-based methods in tandem with interpretable machine learning techniques in order to provide the clearest insights in to risk adjustment for observational cardiovascular datasets.


Data sources and study population

We conducted an observational cohort study using a complete national sample of fee-for-service (FFS) Medicare beneficiaries in 2015 and 2016. We included all FFS patients who were newly diagnosed with concomitant peripheral artery disease and diabetes mellitus in 2015 and were United States residents between the ages of 65 and 95 at the time of diagnosis.

Their first claim in Medicare (the Physician Services/Carrier [Part B], Outpatient, or the Medicare Provider Analysis and Review [MedPAR] files) containing their first annual diabetes related International Classification of Diseases Ninth and Tenth Revisions (ICD-9 or ICD-10) code was used to define their index date. We also required patients to be FFS for 1 year following their index date to ensure we could observe their outcomes during the entire follow-up period. Additionally, patients had to have a diagnosis of an ulcer during the first 6 months of the index year so that we could ensure the ulcer occurred post diagnosis but before the outcomes of interest. Finally, patients had to be outcome-free (alive, no reinterventions or amputations) for at least 6 months following their index date so that the windows for exposure and outcome observation were distinct.


We gathered baseline health characteristics from the Medicare claims data, including patient-level comorbidities as determined by the individual disease indicators comprising the Charlson comorbidity index [26, 27]. We used the Medicare Master Beneficiary Summary File (MBSF) to determine each patient’s sex, race, age at diagnosis, Medicare-Medicaid dual-eligibility status, and used their ZIP code to determine whether they were rural or urban residing [28]. Our main clinical exposure variables include whether patients received at least one or more of the following integrated measures in the 6 months following their index date: HbA1c test, foot exam, or vascular imaging study.

Outcome measures

We identified outcomes for the cohort in the 6 months following the exposure window and excluded any patients that had outcome events during the exposure window. Outcomes studied include any lower extremity amputation (both minor and major), any reintervention, and death. Death was determined from the Medicare MBSF while reintervention and amputation were determined using Current Procedure Terminology codes in the Part B and Outpatient files.

Statistical analysis

We applied traditional statistical methods to describe the demographic and health characteristics of the cohort, including mean and standard deviations for continuous variables including age and comorbidity count, and counts and percentages for categorical variables (sex, race, exposures, Medicare-Medicaid dual-eligibility status, urban/rural indicator, and outcomes).

Logistic regression analysis

We fit multivariable logistic regression models to assess statistical association between the outcomes of interest (amputation, death, and reintervention) and patient characteristics including demographics and comorbidities. Model-estimated odds ratios (ORs) and p-values with a significance level of 0.05 are presented. We compared the relative importance of each predictor in the logistic regression models using the absolute value of the t-statistic for each model parameter. In traditional logistic regression, the t-statistic is the parameter estimate divided by the standard error; as the significance of the parameter estimate is based on the t distribution, often quantified by the p-value, the larger the t-statistic the more significant the predictor. The goodness-of-fit for each model was determined using the McFadden pseudo-R[29]. A McFadden pseudo-R2value between 0.2 and 0.4 represents excellent fit, where values close to 0 represents weak fit [30].

Random forest (RF) modeling

Random forest (RF) models are collections of prediction trees, wherein many trees form a “forest” which can be used to provide a large number of trees to divide data elements [1, 10, 31,32,33,34]. A prediction tree is a non-linear approach to modeling complex data, which partitions the covariates into optimal splits until it achieves partitions allowing for the most homogenous subnodes.

How random forest models work

Random forest algorithms grow many prediction trees to create the forest. First, the algorithm selects a training set. In this analysis, we selected 2/3rds of the data, to fit each tree. The remaining 1/3rd of the data is used as a test set to calculate the out-of-bag (OOB) error, an unbiased estimate of the classification error as trees are added to the forest with each iteration. That is, we used the estimated model to predict the outcome for those observations in the test set (those not used to estimate the model) and compared it to the observed truth. The OOB error is then the proportion of the test set that the model predicts incorrectly. We used the randomForest library [35] in R to fit the RF algorithm to our data. See Additional file 1: Appendix Figs. 1 and 2 for a depiction of the algorithm, it’s associated parameters, and application in prediction.

Fig. 1
figure 1

Variable importance by outcome from (A) logistic regression models and (B) averaged over the 500 iterations of the RF algorithm

Fig. 2
figure 2

Accumulated Local Effects (ALE) plots for the top 5 most predictive covariates in the RF model fit to amputation, hospitalization, and death. A Categorical ALE Plots. B Continuous ALE plot for Comorbidity Count for Death

Rare outcomes in RF modeling

A consideration in using a RF model is the use of an imbalanced classification dataset. That is, the observed outcomes in our dataset are rare events, thus the non-events are dominant in the dataset. If we fit the algorithm to the dataset as is, typically it would predict all observations in the test sample to have non-events, making the classification error appear small (equivalent to the event rate) and the resulting algorithm would appear to fit the data well while in reality, it would not identify any useful information in predicting events. With imbalanced datasets, the RF algorithm does not ascertain the necessary information about the rare events to make an accurate prediction. Hence, it is desirable to use RF algorithms with balanced data sets [36]. To fix the imbalance in the events, we performed a combination of oversampling (resampling patients with the event to balance the data) and under-sampling the non-events to obtain the same sample size as the original cohort. Both techniques use a random sampling algorithm to select cases to comprise the final analytic cohorts with balanced numbers of outcome events and non-events [37, 38].

Variable importance in RF modeling

One key feature of RF models is their ability to produce a measure of variable importance. Variable importance in the RF models is estimated by examining how the prediction error increases when the data for each individual covariate is permuted while all other covariates are left fixed. This approach was repeated for each tree grown in the forest and then averaged over all trees. The difference in the number of correctly predicted outcomes in the variable-permuted test data from the number of correctly predicted cases in the original test data gives the variable importance for the individual tree. The average of this difference over all the trees in the forest is the raw importance score for each variable. This prediction error is calculated for each iteration of the algorithm, for each tree generated, and normalized by the standard deviation then averaged over all trees [1, 35].

Variable effects in RF modeling

To describe how covariates influence the prediction ability of the random forest model, we use accumulated local effects (ALE) plots [39]. ALE plots describe how the covariates, commonly called features, influence the prediction in machine learning models, including RF models. For ALE plots, the cumulative effects of a given predictor value are calculated over a conditional distribution to quantify how the effect of a predictor on a target variable/outcome varies with the predictor’s value. The basic interpretation of an ALE plot is, conditional on the given value of a predictor, the relative effect of changing the feature on the prediction is the value given by the ALE plot, as all ALE plots are centered at zero (the sample means for the given predictor). Thus, the value for a given predictor on the ALE plot is interpreted as the difference from that covariate value to the sample mean prediction.

Comparison of logistic regression and RF approaches

To compare the results of the RF analysis to more conventional regression-based methods, we iteratively estimated multivariable logistic regression analysis for each of the outcomes of interest on the training and test data in each iteration of the RF algorithm. As with the RF models, we controlled for age, sex, race, comorbidities, urban/rural indicator, Medicare-Medicaid dual eligibility status, and each of the exposures of interest (HbA1c, foot exam, and vascular imaging). We fit the model on the same over/under sampled data used for the RF algorithm. We produced the prediction error rate of the estimated logistic regression model by predicting each outcome from the estimated probability for each model on all patients in the test set using a probability cutoff of 0.5.

We compared this to the prediction error rate for the same data using the RF methodology. We also compared the variable importance rankings of the covariates as determined by the two methods. This comparison emphasizes the importance of the multipronged approach to analyzing the data for both associations (regression) and predictive importance (RF) in a cohort such as this where areas for intervention are unclear.

All analyses were performed using SAS version 9.4 (Cary, NC) and R version 3.6.1 [40]. This study was approved by our Institutional Review Board, the Dartmouth-Hitchcock Health Human Research Protection Program (STUDY00030829) with waiver of informed consent.


Characteristics of the Cohort

From 2015, there were 88,898 FFS Medicare beneficiaries diagnosed with concurrent PAD and diabetes with an ulcer in our cohort (Table 1). In the original cohort, before over/under sampling, 25% (n = 22,235) of patients were Medicare-Medicaid dual-eligible in the year of their diabetes diagnosis. The majority of patients in the cohort were white (82%), female (53%), and urban residing (81%). The average age of the cohort was 76.6 years, with an average comorbidity count of 1.73 using the Charlson comorbidity index. The rate of preventative treatments for these patients in the first six months following diagnosis were 52% (n = 45,971) with foot exams, 43% (n = 38,393) had vascular imaging, and 50% (n = 44,181) had an HbA1c test. Finally, the rate of outcomes among those in the original cohort included a mortality rate of 4.5% (n = 3,977), a reintervention rate of 3.6% (n = 3,175), and an amputation rate of 1.6% (n = 1,407). The over/under sampling balanced the cases and non-cases for each outcome with the final analytic cohorts for each outcome given in Table 1.

Table 1 Cohort characteristics for patients with a concomitant diagnosis of PAD and diabetes with an ulcer in 2015, overall and after over/under sampling each outcome (N = 88,898)

Logistics regression results

Variable effects and variation explained

We examined the association of the covariates with each of the outcomes using multivariable-adjusted logistic regression models (Table 2). We found for amputation all covariates were significant predictors. A patient being urban-residing had a protective effect against amputation (OR = 0.59; p-val < 0.001). The covariate with the largest odds of amputation is receipt of HbA1c testing (OR = 1.36; p-val < 0.001) followed closely by foot exam (OR = 1.16; p-val < 0.001). Vascular imaging decreases the odds of amputation (OR = 0.90; p < 0.001). Age is protective against amputation (for 85 + year olds vs. 65–69 OR = 0.29; p-val < 0.001) with the likelihood of amputation decreasing as age decreases. Results were similar for reintervention and death.

Table 2 Logistic regression results for models fit on the over/under-sampled cohorts for each outcome

The McFadden R2 = 0.071 indicating approximately 7% of the total variation in amputation risk is explained by the model, falling outside the 0.2–0.4 range representing excellent model fit. For reintervention, the McFadden R2 = 0.046, which signifies less than 5% of the observed variation is explained by the model. Finally, the McFadden R2 = 0.124 for death, which means the model for death explains almost double the variation as explained for amputation and over double than explained for reintervention. Thus, the McFadden R2 values for each of the logistic regression models indicated that the chosen predictors did not explain a significant portion of the observed variance using these models.

Variable importance

Variable importance, as assessed by the absolute value of the t-statistic (Fig. 1a), tells us that age (85 + , t = -53.17), sex (t = -40.52), and urban/rural status (t = -31.33) are the most important predictors of amputation in the logistic regression model. While these factors were important, factors including dual eligibility status, receipt of a vascular image, and race are the least predictive of amputation. Similar findings for reintervention and death outcomes are shown as well (Fig. 1a).

Random forest results

Variable importance

In Fig. 1b, we see the results of the variable importance calculations from the RF estimation for each of the 3 outcomes. For the prediction of amputation in the cohort, the most important variable is whether they are rural- or urban-residing. We see the level of importance for the other covariates including age and sex follow closely behind. The first clinical variable, HbA1c testing, is the fourth most important variable in the prediction of amputation followed by dual-eligibility. As with our logistic regression models, similar results are shown for reintervention and death.

Variable effects

In Fig. 2a, the ALE plot for amputation indicates that being rural-residing increases the likelihood of amputation compared to the mean population, while being urban-residing decreases it. Similarly, as age (category) increases, the less likely amputation is to occur compared to the mean population. Females are less likely than the mean population to receive an amputation, as are those who are not dual-eligible. Those not receiving an HbA1c test are also less likely than the mean population to receive an amputation.

For reintervention, certain differences were evident. The younger the patient, the more likely they are to have a reintervention compared to the mean population. Patients who are non-white are more likely to have a reintervention than the mean population while white patients are less likely. Finally, for death, we see the older the patient is, the more likely they are to die. Receiving an HbA1c test, being dual-eligible, and male all increase the likelihood of death compared to the mean population. In Fig. 2b, we see the ALE plot for death as comorbidity count varies. We see the ALE plot crosses zero at approximately 2, meaning those with 2 or fewer comorbidities are less likely to die than the general population and with those with greater than 2 comorbidities having a greater likelihood of death than the mean population.

Comparison of logistic regression and random forest approaches

Variable importance: RF models can help highlight hidden themes

Assessing variable importance of each model produced similar results from each approach. The top 3 variables for each outcome in each model are the same, except for death where the logistic regression model indicated comorbidity count as a highly significant predictor, and where this predictor is ranked least important in the RF model, and sex was included in its place in the RF model, following a close 4th in the logistic regression model. However, the variation in the variable importance measure in the RF models is noticeably less than in the logistic regression models. This indicates that the predictive power of each variable according to the RF model is much closer, highlighting the necessity for comprehensive analyses and importance of non-dominant subgroup analysis. For example, age (or at least 1 level of age in the logistic regression models) is the most significant predictor (or the second most) in both models for amputation, death, and reintervention. Notably, age is not the most significant predictor for the RF model for amputation. The sensitivity of the approach was able to tease out a key difference in amputation rates for rural- and urban-residing residence; the logistic regression model did identify this difference as well, following age and sex. This illustrates how the RF model can potentially tease out differences highlighting significant and non-dominant subpopulations for intervention or further study.

Variable effects: RF models and logistic regression often have similar effect sizes

For amputation, for the top 5 most predictive covariates determined in both techniques, we see exact concordance in the directionality of the conclusions. For example, in the RF model, rurality increases the risk of amputation compared to the general population, while the logistic regression results show a ratio less than 1 for urban vs. rural indicating rurality decreases the odds of amputation compared to urban-residing population.For death and reintervention, we see similar concordance among the directionality of conclusions according to the logistic regression results and the ALE plots for the RF models.

Predictive capabilities: RF models often illustrate lower error rates

Finally, in Fig. 3, we see the average results of prediction on the test sample created in each iteration of the RF algorithm. The OOB error rate, on average, when predicting amputation is 31%, while for the same test datasets is 63% using the estimated logistic regression models. Similarly, for reintervention the RF OOB error rate is 36% while for logistic regression it is 60%. For death, we see an OOB error rate of 30% while logistic regression has a prediction error rate of 68% when predicting for test data, averaged over 500 iterations.

Fig. 3
figure 3

Average prediction error rates for the random forest (RF) and logistic regression models. RF prediction error rate is the out-of-bag (OOB) error rate. Logistic regression test sample prediction error rate calculated using average error rate across 500 iterations


In this study, we analyzed factors predicting amputation, death, and reintervention in a cohort of patients with peripheral artery disease and diabetes mellitus, and an ulcer using both a logistic regression approach and a random forest approach. In general, both approaches provided similar results. However, the RF approach illustrated differences in certain subgroups, highlighting significant and non-dominant subpopulations for intervention or further study. Specifically, in the RF approach, we found patient age 65–69 was most predictive of amputation, as shown on the categorical ALE plots in Fig. 2, but whether the patient is rural- or urban-residing is the most predictive observed covariate (Fig. 1b). In other words, though urban/rural status is the most important predictor of amputation, and being rural leads to an increased risk of amputation compared to being urban, being in the youngest age group confers the largest single increased risk of amputation. In the RF model, HbA1c testing and comorbidity count follow age in the most predictive covariates of death in this cohort, where the HbA1c test is given more often in those who die, a consequence of healthy patients not receiving monitoring as closely as those symptomatic patients, and the likelihood of death increases with the patient’s comorbidities.

Using advanced statistical methods to complement common approaches

RF algorithms have a number of statistical and computational strengths. RFs are versatile in terms of the structure and types of data analyzed including both regression and classification approaches to estimation [35]. The algorithms use of only a subset of features at a time allow it to process significantly faster than other machine learning algorithms and ultimately allows for rapid training and prediction across many trees. Additionally, unlike typical regression approaches including logistic regression, RFs are robust to outliers with their predictive power little influenced [41]. Finally, each tree in the RF has high variance and low bias. Averaging across many such trees to create the final RF model results in a final model with low bias and only moderate variance [42].

While the relationship between the outcomes and the covariates given by the logistic regression results and the random forest models are similar in directionality, the RF results offer two key advantages to the regression approach. First, the RF approach offers a substantive and sensitive approach to identifying which covariates are most predictive and have the greatest importance in predicting the outcomes of interest. Second, the RF models offered increased predictive power over the classic logistic regression while using the same covariates. In fact, the prediction error rate was twice as high in logistic regression models than in the RF models. These two properties, when paired with the common interpretability of the logistic regression approach, provide a framework for analyzing data in a population with clear magnitude of associations between the outcomes and covariates (regression) and the ability to deeply analyze which covariates are most predictive of the outcome and best suited for clinical consideration and/or intervention (RF). Potential extensions to the current analysis would be to include more clinical factors in the RF model to increase the predictability and more modifiable patient or care patterns to potentially prevent or delay negative outcomes for this cohort.

On the other hand, there are some limitations to using RF models. While RFs excel at classification, they do not predict beyond the range of the covariates in the training data. Additionally, RFs may overfit datasets that are particularly noisy, that is data with an abundance of unexplained variation. This is potentially problematic in healthcare research where we are limited by the information in the Medicare claims data. The use of over/under sampling adds a layer of complexity and data manipulation that many audiences may not be familiar or comfortable with. It overemphasizes cases (in the current study), stressing the influence of covariates that may ultimately impact a small portion of the true sample population. Finally, the lack of understanding of random forest plots among clinical audiences, may make interpretation and practical application in health care settings less likely. Thus, utilizing RF models in tandem with traditional regression approaches may best serve to establish the use of RF models in healthcare research while also improving the ability to predict potential outcomes and pathways for patients.


The use of RF models to analyze data and provide predictions for patients holds great potential in identifying modifiable patient-level and health-system factors and cohorts for increased surveillance and intervention to improve outcomes for patients. RFs are incredibly high performing models with difficult interpretation most ideally suited for times when accurate prediction is most desirable and can be used in tandem with more typical methods and tools (ALE plots, logistic regression) to provide a more thorough analysis of observational data.

Availability of data and materials

In accordance with our Data Use Agreement with the Centers for Medicare and Medicaid Services, analytic datasets cannot be shared. However, aggregate data elements and coding algorithms will be made freely available upon request. Please contact the corresponding author to request aggregate data elements and coding algorithms.



Critical limb ischemia


Chronic obstructive pulmonary disease


Diabetes Mellitus




Hemoglobin A1c


International Classification of Diseases


Machine Learning


Medicare Provider Analysis and Review


Medicare Master Beneficiary Summary File


Peripheral Arterial Disease


Random Forest


Variable Importance


Vascular Quality Initiative


  1. Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.

    Article  Google Scholar 

  2. Comulada WS, Desmond KA, Gildner JL, Leibowitz AA. Transitioning from medicaid disability coverage to long-term medicare coverage: the case of people living with HIV/AIDS in California. AIDS Educ Prev. 2017;29(1):49–61.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Eraslan G, Avsec Ž, Gagneur J, Theis FJ. Deep learning: new computational modelling techniques for genomics. Nat Rev Genet. 2019;20(7):389–403.

    Article  CAS  PubMed  Google Scholar 

  4. Huang S, Cai N, Pacheco PP, Narrandes S, Wang Y, Xu W. Applications of Support Vector Machine (SVM) learning in cancer genomics. Cancer Genomics Proteomics. 2018;15(1):41–51.

    CAS  PubMed  Google Scholar 

  5. Bellinger C, Mohomed Jabbar MS, Zaïane O, Osornio-Vargas A. A systematic review of data mining and machine learning for air pollution epidemiology. BMC Public Health. 2017;17(1):907.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Tack C. Artificial intelligence and machine learning | applications in musculoskeletal physiotherapy. Musculoskelet Sci Pract. 2019;39:164–9.

    Article  PubMed  Google Scholar 

  7. Deo RC. Machine learning in medicine. Circulation. 2015;132(20):1920–30.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Ngiam KY, Khor IW. Big data and machine learning algorithms for health-care delivery. Lancet Oncol. 2019;20(5):e262–73.

    Article  PubMed  Google Scholar 

  9. Beam AL, Kohane IS. Big data and machine learning in health care. JAMA. 2018;319(13):1317–8.

    Article  PubMed  Google Scholar 

  10. Kirasich K, Smith T, Sader B. Random forest vs logistic regression: binary classification for heterogeneous datasets. SMU Data Sci Rev. 2018;1(3):9.

    Google Scholar 

  11. Cerrito P. The difference between predictive modeling and regression. MWSUG conference. Indianapolis: 2008. Available online at: Accessed 11 Oct 2020.

  12. Peripheral Artery Disease. Health Topics Accessed 11 Dec 2019.

  13. Fowkes FG, Aboyans V, Fowkes FJ, McDermott MM, Sampson UK, Criqui MH. Peripheral artery disease: epidemiology and global perspectives. Nat Rev Cardiol. 2017;14(3):156–70.

    Article  PubMed  Google Scholar 

  14. Buso G, Aboyans V, Mazzolai L. Lower extremity artery disease in patients with type 2 diabetes. Euro J Prev Cardiol. 2019;26(2):114–24.

    Article  Google Scholar 

  15. Goldberg JB, Goodney PP, Cronenwett JL, Baker F. The effect of risk and race on lower extremity amputations among Medicare diabetic patients. J Vasc Surg. 2012;56(6):1663–8.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Goodney PP, Beck AW, Nagle J, Welch HG, Zwolak RM. National trends in lower extremity bypass surgery, endovascular interventions, and major amputations. J Vasc Surg. 2009;50(1):54–60.

    Article  PubMed  Google Scholar 

  17. Taylor SM, Kalbaugh CA, Blackhurst DW, et al. Postoperative outcomes according to preoperative medical and functional status after infrainguinal revascularization for critical limb ischemia in patients 80 years and older. Am Surg. 2005;71(8):640–5.

    Article  PubMed  Google Scholar 

  18. Mayfield JA, Reiber GE, Maynard C, Czerniecki JM, Caps MT, Sangeorzan BJ. Trends in lower limb amputation in the Veterans Health Administration, 1989-1998. J Rehabil Res Dev. 2000;37(1):23-30. PMID: 10847569.

    CAS  PubMed  Google Scholar 

  19. Mayfield JA, Reiber GE, Nelson RG, Greene T. Do foot examinations reduce the risk of diabetic amputation? J Fam Pract. 2000;49(6):499–504.

    CAS  PubMed  Google Scholar 

  20. Mayfield JA, Reiber GE, Sanders LJ, Janisse D, Pogach LM. Preventive foot care in diabetes. Diabetes Care. 2004;27(Suppl 1):S63-64.

    PubMed  Google Scholar 

  21. Czupryniak L. Guidelines for the management of type 2 diabetes: is ADA and EASD consensus more clinically relevant than the IDF recommendations? Diabetes Res Clin Pract. 2009;86(Suppl 1):S22-25.

    Article  PubMed  Google Scholar 

  22. Rooke TW, Hirsch AT, Misra S, et al. 2011 ACCF/AHA focused update of the guideline for the management of patients with peripheral artery disease (updating the 2005 guideline): a report of the American college of cardiology foundation/American heart association task force on practice guidelines. J Am Coll Cardiol. 2011;58(19):2020–45.

    Article  PubMed  Google Scholar 

  23. Bruckner M, Mangan M, Godin S, Pogach L. Project LEAP of New Jersey: lower extremity amputation prevention in persons with type 2 diabetes. Am J Manag Care. 1999;5(5):609–16.

    CAS  PubMed  Google Scholar 

  24. LEAP Program (Lower Extremity Amputation Prevention). Med Health R I. 1998;81(11):359-60. PMID: 15580793.

  25. American Diabetes Association. Standards of medical care in diabetes--2014. Diabetes Care. 2014;37 Suppl 1:S14-80.

  26. Quan H, Li B, Couris CM, et al. Updating and validating the charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173(6):676–82.

    Article  PubMed  Google Scholar 

  27. Quan H, Sundararajan V, Halfon P, et al. Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data. Med Care. 2005;43(11):1130–9.

    Article  PubMed  Google Scholar 

  28. Hadden, Louise S. and Michael S. Zdeb. “ZIP Code 411: A Well-Kept SAS® Secret.” 2005.

  29. McFadden D, University of California BIoU, Development R. Conditional Logit Analysis of Qualitative Choice Behavior. Institute of Urban and Regional Development, University of California; 1973. Available at:

  30. Mcfadden D. Quantitative Methods for Analyzing Travel Behaviour of Individuals: Some Recent Developments. 1977.

  31. Djekic D, Fagman E, Angerås O, Lappas G, Torén K, Bergström G, Rosengren A. Social Support and Subclinical Coronary Artery Disease in Middle-Aged Men and Women: Findings from the Pilot of Swedish CArdioPulmonary bioImage Study. Int J Environ Res Public Health. 2020;17(3):778.

  32. Biau G. Analysis of a random forest model. J Mach Learn Res. 2012;13:1063–95.

    Google Scholar 

  33. Denil M, Matheson D, Freitas N. Narrowing the gap: random forests in theory and in practice. 31st international conference on machine learning, PMLR. 2014;32(1):665-73. Available online at:

  34. Sarica A, Cerasa A, Quattrone A. Random Forest Algorithm for the Classification of Neuroimaging Data in Alzheimer's Disease: A Systematic Review. Front Aging Neurosci. 2017;9:329.

  35. Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2(3):18–22.

    Google Scholar 

  36. Team AV. Practical guide to deal with imbalanced classification problems in R. In: Analytics Vidhya. 2016.

    Google Scholar 

  37. Menardi G, Torelli N. Training and assessing classification rules with imbalanced data. Data Min Knowl Disc. 2014;28:92–122.

    Article  Google Scholar 

  38. Lunardon N, Menardi G, Torelli N. ROSE: a package for binary imbalanced learning. R Journal. 2014;6:82–92.

    Article  Google Scholar 

  39. Molnar C. Interpretable Machine Learning: A Guide for Making Black Box Models Explainable. Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License2021.

  40. Team RC. R: A language and environment for statistical computing. . 2019;

  41. Gunduz, Necla and Ernest P. Fokoue. “Robust Classification of High Dimension Low Sample Size Data.” arXiv: Applications (2015): n. pag. Available online at:

  42. Hastie T, Tibshirani R, Friedman J. Random forests. In: The elements of statistical learning. Springer New York; NY 2009:587–604.

Download references




The authors wish to acknowledge the American Heart Association and its support via its Strategically Focused Research Networks Grant (Bonaca PI, AHA 18SFRN3390008).

Author information

Authors and Affiliations



AA: concept and study design, analyses, interpretation and writing. NR: interpretation and critical revision. BG: analyses and interpretation; and critical revision of the manuscript. JAB: interpretation and writing. ME: interpretation and writing. KM: critical revision of the manuscript and administrative and technical support. MF: interpretation and critical revision. MC: interpretation and critical revision. MB: concept and study design, data acquisition, interpretation and critical revision. PG: concept and study design, data acquisition, interpretation, and writing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Philip P. Goodney.

Ethics declarations

Ethics approval and consent to participate

We accessed the data using Centers for Medicare and Medicaid Services (CMS) data use agreement RSCH-2018–52223. In accordance with this, all data was subject to standard CMS data suppression regulations. This study was approved by our Institutional Review Board, the Dartmouth-Hitchcock Health Human Research Protection Program (STUDY00030829) with waiver of informed consent.

Consent for publication

Not Applicable.

Competing interests


Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Appendix: Figure 1.

Random forest algorithm for prediction. Figure 2. Decision tree construction and prediction from randomforest algorithm. 

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Austin, A.M., Ramkumar, N., Gladders, B. et al. Using a cohort study of diabetes and peripheral artery disease to compare logistic regression and machine learning via random forest modeling. BMC Med Res Methodol 22, 300 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: