Development of a Machine Learning Algorithm for Prediction of Complications and Unplanned Readmission Following Primary Anatomic Total Shoulder Replacements.

Background: The demand and incidence of anatomic total shoulder arthroplasty (aTSA) procedures is projected to increase substantially over the next decade. There is a paucity of accurate risk prediction models which would be of great utility in minimizing morbidity and costs associated with major post-operative complications. Machine learning is a powerful predictive modeling tool and has become increasingly popular, especially in orthopedics. We aimed to build a ML model for prediction of major complications and readmission following primary aTSA.
Methods: A large California administrative database was retrospectively reviewed for all adults undergoing primary aTSA between 2015 to 2017. The primary outcome was any major complication or readmission following aTSA. A wide scope of standard ML benchmarks, including Logistic regression (LR), XGBoost, Gradient boosting, AdaBoost and Random Forest were employed to determine their power to predict outcomes. Additionally, important patient features to the prediction models were indentified.
Results: There were a total of 10,302 aTSAs with 598 (5.8%) having at least one major post-operative complication or readmission. XGBoost had the highest discriminative power (area under receiver operating curve AUROC of 0.689) of the 5 ML benchmarks with an area under precision recall curve AURPC of 0.207. History of implant complication, severe chronic kidney disease, teaching hospital status, coronary artery disease and male sex were the most important features for the performance of XGBoost. In addition, XGBoost identified teaching hospital status and male sex as markedly more important predictors of outcomes compared to LR models.
Conclusion: We report a well calibrated XGBoost ML algorithm for predicting major complications and 30-day readmission following aTSA. History of prior implant complication was the most important patient feature for XGBoost performance, a novel patient feature that surgeons should consider when counseling patients.

Introduction

Anatomic total shoulder arthroplasty (aTSA) has been shown to reliably improve pain and range of motion in patients with severe glenohumeral osteoarthritis (OA) and a functioning rotator cuff. Compared to hemiarthroplasty (HA) and reverse total shoulder arthroplasty (rTSA), aTSA has been shown to have lower surgical complication and readmission rates.[2-5] Over 40,000 primary aTSA procedures are performed in the United states each year. With an aging population and advancements in implants, the annual incidence of aTSA is projected to increase by up to 50% by 2025.[6,7] However, given the significant cost and morbidity inherently associated with complications and unplanned readmissions, accurate risk stratification of patients who undergo aTSA would be of great utility.

Use of machine learning (ML) methods to generate prediction models has become increasingly popular within orthopaedics due to their ability to detect complex non-linear relationships and identify novel predictive factors.[8-11] Patient records, especially billing based data sets, contain large amounts of quantitative and qualitative data which introduces too many variables for traditional regression models to perform optimally. Recent studies have shown that ML algorithms outperform linear models and commonly used indices in predicting outcome such as Patient Reported Outcome Measures (PROMs), readmissions, and extended length of stay.[13-17] An accurate ML based prediction model for aTSA has the potential for improving preoperative decision making, informed consent, post-operative outcomes and help guide outcome-based performance measures and reimbursement programs.

Our primary aim was to build a ML model for prediction of major perioperative complications and readmission following primary aTSA for OA. Secondarily, we aim to compare the performance of our ML models to traditional logistic regression (LR). Lastly, we aim to compare the relative importance of clinical patient features that predict outcomes in our best performing model to the most important predictive patient features from logistic regression models. We hypothesized that we will develop a ML model that outperforms LR and identify novel patient features that are important for prediction of major complications and readmissions following primary aTSA.

Methods

Data Source

Data were obtained from California's Office of Statewide Health Planning and Development (OSHPD) database, a mandatory statewide database containing codes for up to 24 diagnoses and 20 inpatient procedures per hospitalization from all licensed nonfederal hospitals in California. The OSHPD database includes patient and hospital characteristics including age, gender, race/ethnicity, insurance type, multiple comorbidities, and hospital volume. Patients in this database are assigned unique record linkage numbers that allows patients to be tracked longitudinally for complications regardless of whether future admission are at a different hospital in the database from where the index procedure was performed.

Inclusion and Exclusion Criteria

The OSHPD database was retrospectively reviewed to select patients older than 18 from October first 2015 to December 13th 2017 who underwent elective primary aTSA using International Classification of Diseases, Tenth Revision, procedural codes (ICD-10-CPS) codes. The exclusion criteria included patients with fracture of the upper extremity/shoulder girdle coded in the principal or secondary discharge diagnosis fields of index admission; concurrent revision, resurfacing, or implanted device/prosthesis removal procedure; mechanical complications coded in the principal discharge diagnosis field; malignant neoplasm of the upper extremities/shoulder girdle, bone/bone marrow or a disseminated malignant neoplasm coded in the principal discharge diagnosis fields. All principal inclusion codes were: 0RRJ0JZ, 0RRJ0KZ, 0RRK0JZ and ORRKOKZ. The extensive exclusion ICD-10 codes can be made available upon request to the authors.

Outcome and Explanatory Variables

The primary outcome of interest was any major complication or unplanned readmission after index primary aTSA (Table 2). Major complications were identified using ICD-10 codes adapted from performance measures developed by the Centers for Medicare and Medicaid (CMS) for total joint replacement. These include acute myocardial infarction, pneumonia, sepsis, pulmonary embolism, surgical site bleeding, wound infection and mechanical complication. Myocardial infarction, pneumonia, and sepsis were included if the complication occurred during the index admission or within seven days of start of index admission. Pulmonary embolism was included if it occurred during the index admission or within 30 days of admission. Surgical site bleeding, wound infection and mechanical complication were included during the index admission or within 90 days. Readmission for any cause within 30 days of index aTSA was also included as an outcome.

Table 2.

Major Complications and Readmission.

Complications	All Patients(n = 10,302)
	Number (%)
At least one complication or readmission	598 (5.8)
Readmission within 30 days	400 (3.88)
Wound infection	157 (1.52)
Sepsis	38 (0.37)
Mechanical complication	4 (0.04)
Pneumonia	83 (0.81)
Pulmonary embolism	27 (0.26)
Surgical site bleeding	34 (0.33)
Acute myocardial infarction	23 (0.22)

The patient features (explanatory variables) included in or derived from the OSHPD database include patient demographic characteristics (ie age, sex, race, ethnicity, body mass index, insurance type), hospital type (community vs. teaching hospital) and patient medical comorbidities using the CMS condition categories as defined by the CMS Hierarchical Condition Category (HCC) risk adjustment model (ie diabetes, coronary artery disease, chronic obstructive pulmonary disease, malignancy, renal failure).

Model Development

We utilized 5 standard ML benchmarks that cover different classes of ML modeling as follows: LR (linear classifier), random forest (a tree-based ensemble classifier), AdaBoost, gradient boosting machines (Gradient Boosting), and XGBoost (boosting ensemble classifiers). We implemented LR, Random Forest, AdaBoost, and Gradient Boosting machines using the scikit-learn Python library and XGBoost using the xgboost Python library. The hyperparameters (which define the mathematical limits of an ML algorithm) of each model were selected via grid search: for LR, the coefficient for L2 regularization was chosen from a set of values in a logarithmic scale between 1e-3 to 1e3; for Random Forest, Adaboost, Gradient Boosting, and XGBoost the number of trees and the maximum depth of each tree were selected from {50, 100, 200, 300} and {2, 3, 4, 5}, respectively.

Model Evaluation

In statistical modeling, discrimination refers to how well a model distinguishes patients who developed post-operative complications and those who did not., while calibration refers to the level of agreement between prediction and the observed outcomes. We evaluated the discriminative and calibration performances of the prognostic models via 5-fold stratified cross-validation to avoid overfitting. In every cross-validation fold, the training cohort (80% of the study population) was used to derive our 5 ML benchmark models, and then a held-out testing cohort (20% of the study population) was used for performance evaluation.

Discrimination was assessed using area under the receiver operating characteristic curve (AUROC). AUROC represents the probability that a randomly selected patient who experienced an outcome was assigned a higher risk by the model than a patient who did not experience the outcome. An AUROC of 0.5 indicates that a prognostic model has no discriminative power while an AUROC of 1 indicates that a prognostic model provides perfect discrimination. Calibration was assessed using Brier scores: a measure of the agreement between the observed binary outcome and the predicted probability of that outcome, which is equivalent to the mean squared error. Lower Brier scores indicates better calibration of the prognostic model.

In addition to AUROC, we also determined area under the precision-recall curve (AUPRC) values which is a useful performance metric when analyzing an imbalanced dataset; that is, a dataset where negative cases far outnumber positive cases. The precision-recall (PR) curve is constructed by plotting positive predictive value (precision) versus the sensitivity (recall). The PR curve focuses on identifying the ability of the model to correctly identify positive cases; it ignores true negatives, which is the dominant group in an imbalanced dataset.[25,26] Unlike the AUROC, the baseline AUPRC is the proportion of true positive cases. An ideal classifier predicts every positive case (perfect recall) without marking any negative case as positive (perfect precision) and will return an AUPRC of 1. Random prediction will result in the baseline AUPRC. The further the AUPRC is from the random prediction value, the better the model handles positive cases.

AUROC, AUPRC and Brier scores were reported as mean values with standard deviations (SD).

Feature Importance

We utilized the partial dependence function introduced by Friedman et al. 2001 to measure the importance of an individual feature by assessing the average effect in predicted risks when its value is perturbed (Appendix I). The continuous variables were standardized to zero mean and unit variance, and the categorical variables were one-hot encoded.

Results

Demographic Characteristics

Between 10/01/2015 to 12/13/2017, there was a total of 10, 302 primary aTSAs, the majority of which were females (54%). Patient age ranged from 45 to 98 years old with a median age of 73. Overall demographics and some of the most common medical comorbidities are summarized in Table 1. A total of 598 (5.8%) patients had at least one complication or readmission. There were 400 (3.9%) patients who required readmission within 30 days. The most common complications were wound infection, pneumonia and sepsis (Table 2).

Table 1.

Baseline Cohort Demographics.

Variable	All Patients(n = 10,302)
	Median (IQR)
Age (years)	71 (12)
Hospital volume†	103 (141)
	Number (%)
Male	4727 (45.88)
Race
White	8835 (85.76)
Black	333 (3.23)
Asian / Pacific Islander	215 (2.09)
Native American	51 (0.49)
Other	777 (7.54)
Unknown	91 (0.88)
Ethnicity
Non-Hispanic	8900 (86.39)
Hispanic	1309 (12.71)
Unknown	93 (0.90)
Insurance
Medicare	7433 (72.15)
Private	1831 (17.77)
Medi-Cal	393 (3.81)
Workers’ compensation	518 (5.02)
Other	127 (1.23)
Medical comorbidities
Diabetes mellitus without complications	739 (7.17)
Diabetes mellitus with complications	662 (6.43)
Coronary atherosclerosis	719 (6.98)
Morbid obesity	664 (6.45)
COPD	700 (6.79)
Chronic kidney disease, mild	682 (6.62)
Chronic kidney disease, moderate	621 (6.03)
Chronic kidney disease, severe	553 (5.37)
Chronic kidney disease requiring dialysis	549 (5.33)
Vascular disease	662 (6.43)
Other circulatory disease	623 (6.05)
Acute renal failure	650 (6.31)
Cardio-respiratory failure	603 (5.85)
Major depressive or bipolar disorder	636 (6.17)
Major fracture (except skull)	574 (5.57)
Hip fracture or dislocation	554 (5.38)
Protein-calorie malnutrition	573 (5.56)
Metastatic cancer or leukemia	544 (5.28)
Complications of implants	708 (6.87)
History of prior complications	595 (5.78)
Osteoarthritis of hip or knee	737 (7.15)
Osteoporosis	667 (6.47)
History of bone/joint/muscle infection	589 (5.72)
	Mean (SD)
Number of comorbidities	0.34 (1.19)

IQR = interquartile range; COPD = chronic obstructive pulmonary disease; SD = standard deviation

† Cases of primary aTSA performed between 10/1/2015 and 12/13/2017

Model Performance and Calibration

XGBoost demonstrated higher discrimination compared to LR (AUROC 0.689 vs. 0.662) as well as outperforming the other three standard benchmark models (Table 3). XGBoost is well-calibrated with Brier score of 0.051. The LR and standard ML models are similarly well-calibrated with the exception of AdaBoost. XGBoost and Gradient Boosting models had the highest AUPRC values of 0.207 and 0.214. These values are compared against a random classifier for this cohort of 0.058. The receiver operating characteristic curves and precision recall curves of the XGBoost and logistic regression models are depicted in Figures 1 and 2 respectively.

Table 3.

Discrimination and Calibration.

Model	AUROC	AUPRC	Brier score
XGBoost	0.689 ±  0.026	0.207 ±  0.0.44	0.051 + 0.002
Logistic Regression	0.662 ± 0.026	0.137 ± 0.024	0.055 ± 0
Gradient Boosting	0.687 ± 0.027	0.214 ± 0.049	0.051 ± 0.002
AdaBoost	0.677 ± 0.013	0.199 ± 0.049	0.245 ± 0.002
Random Forest	0.624 ± 0.022	0.121 ± 0.016	0.061 ± 0.001

Figure 1.

Area Under Receiver Operating Curve. Receiver operating characteristic curves for XGBoost and logistic regression

Figure 2.

Area Under Precision Recall Curve. Precision Recall Curves of XGBoost and logistic regression.

Relative Feature Importance

Given XGBoost had the highest discriminative performance based on AUROC, the relative importance of the categorical, continuous and top ten binary variables to the performance of the XGBoost model are listed in Table 4. A history of implant complications was the most important predictive feature for both XGBoost and LR. Specific to XGBoost, history of implant complication was more than twice as important for the model as the second most influential feature (chronic kidney disease, severe). Overall, the majority of variables that are important for XGBoost differed from those that were important for LR.

Table 4.

Relative Feature Importance for Complications or Readmission Following Primary aTSA.

Feature	Rank in XGBoost(Rank in Logistic Regression)	Change to Risk Prediction
Binary features
Complication of implants	1 (1)	0.0297
Chronic kidney disease, severe	2 (23)	0.0117
Teaching Hospital	3 (63)	0.0111
Coronary atherosclerosis	4 (2)	0.0076
Male sex	5 (64)	0.0070
Dimentia without complications	6 (21)	0.0060
Other circulatory diseases	7 (14)	0.0059
Osteoarthritis of hip or knee	8 (31)	−0.0058
Osteoporosis	9 (3)	0.0051
Morbid Obesity	10 (18)	−0.0029
Continuous features
Number of medical comorbidities	1 (1)	0.0650
Age	2 (3)	0.0278
Hospital volume	3 (2)	−0.0101
Insurance status
Medicare	Reference	0
Private	1 (1)	−0098
Medical	2 (2)	0.0051
Workers Comp	3 (3)	−0.0008
Other	3 (4)	−0.0008
Race
White	Reference	0
Asian/Pacific Islander	1 (2)	−0.0312
Black	2 (1)	0.0124
Other	3 (4)	0.0002
Native American	4 (3)	5.02E-05
Unknown	4 (5)	5.02E-05

Discussion

Due to an aging population and improved outcomes, the demand for shoulder arthroplasty (including aTSA) is projected to increase drastically over the next 5 to 10 years.[6,27] With this increase in prevalence, mitigating the morbidity and costs associated with post-operative complications and readmissions is increasingly important. The purpose of this study was to create an ML algorithm to predict perioperative complications following aTSA using a statewide retrospective database. We found that XGBoost produced the most accurate predictive model and, by analyzing relative feature importance, that a history of prior implant complications was the most important predictive patient feature.

Until recently, multivariate LR has been the most prevalent modeling method used for outcome prediction.[28-30] However, ML, a subset of artificial intelligence, has grown in popularity due to advanced detection of complex non-linear relationships and factor-factor interactions within a given dataset.[31,32] Accordingly, ML models have been shown to outperform LR in many cases across different medical and surgical specialties.[33-35] Substantial effort has gone into utilizing ML to develop various prediction models of outcomes in orthopaedic surgery.[36-38] Clinically, developing prediction models using ML can be very useful for physicians and patients when it comes to setting expectations, preparing for particular outcomes and managing complications. Accurate ML models can also be used for allocating physician reimbursement by differentiating the high risk, complicated patients for whom increased care and costs are expected for a particular procedure.

ML has been sparingly used to predict outcomes following aTSA. Using the American College of Surgeons-National Surgical Quality Improvement Program (ACS-NASQIP) database, Gowd et. al developed an ML tool to aid in patient selection for outpatient TSA based on medical comorbidities and demographic factors. With a cohort of 4500 patients, their random forest ML model was used to predict which patients had a length of stay of 1 day or less with an AUROC of 0.77. Most recently Arvind and colleagues analyzed 9043 patients who underwent primary total shoulder arthroplasty to predict 3-day unplanned readmission rates. They reported C-Statistic scores (analogous to AUROC) of 0.74 using Random Forest and 0.54 using support vector machine.

To our knowledge, this is the first study to utilize XGBoost, a gradient boosting machine learning method, in predicting complications/readmission following TSA. The model had an AUROC of 0.689, which is comparable to the aforementioned similar studies, is well-calibrated, and demonstrates superior discrimination than LR. As previously stated, AUPRC can be a useful performance metric to evaluate models built on imbalanced data set.[25,26] Our XGBoost model demonstrated good performance with an AUPRC of 0.207 compared to a random classifier of 0.058 for this cohort.

We also evaluated the relative importance of 64 different binary patient variables along with continuous variables (patient age and hospital volume) and categorical variables (patient insurance and race/ethnicity). Some ML methods may allow for detection of indirect nonlinear relationships and multivariate effects that others are not able to identify. Therefore, it is important to recognize that prediction models should not be interpreted as explanatory models. Specifically, the magnitude of feature importance should not be taken to imply causal relationships or lack thereof. Rather, inclusion of these features increases the predictive performance of the model. The most important binary feature for both the XGBoost and LR models was a history of implant-related complication. This suggests that surgeons should weigh perioperative complications following previous orthopedic implant surgeries more significantly when counseling prior to aTSA. Severe chronic kidney disease (CKD) was also found to be the second most important predictor in our model. Dacombe et al. reported CKD as a predictor for increased length of stay in their multiple linear regression analysis of 640 shoulder arthroplasty cases. Teaching hospital was the third most important feature, which may be secondary to more complex patients/cases being performed at teaching hospitals. Interestingly, teaching hospital as a feature was markedly less important for the LR model. Similarly, male sex and osteoarthritis of the hip or knee also were found to be important for XGBoost performance but much less so for LR. In regard to continuous variables, a patient's total number of CMS condition categories was the most important predictor. This is consistent with overall comorbidity burden contributing to the risk of complication or readmission following aTSA.

That there were significant differences in feature importance between XGBoost and LR underlines the fact that advanced machine learning methods treat the same features very differently than traditional LR-based methods. For example, male sex was the fifth most influential variable for our XGBoost algorithm whereas it was the 64th most influential variable for LR. This is possibly due to the ability of ML methods to capture occult relationships between variables that LR is unable to detect.

The retrospective nature of this study inherently lends itself to limitations. Though a de-identified state-wide code-based database has a large patient sample, the patient features and outcomes collected are limited. The reliance on ICD diagnoses and procedure codes is less reliable than thorough chart review. Code based searches of databases are dependent on accurate coding and can lead to exclusion of a patient of interest or underestimation of outcomes. With this database we were unable to assess mortality, patient-reported functional outcomes and patient satisfaction. Due to the low complication rate found in this cohort, our data may be imbalanced. However, we believe that predictive models trained with an artificially balanced dataset cannot be directly used in a clinical setting as they will be inherently poorly calibrated. To better address the concern of imbalanced data, we evaluated the five prognostic models in terms of area under precision-recall curve. Along the same lines, due to the low overall complication rate, secondary analysis of individual complications/outcomes was beyond the scope of our model (though we certainly recognize the clinical importance of such results). Lastly, we must acknowledge the black box nature of ML algorithms that can lead to non-physiologic predictors or predictors that effect a very small portion of the cohort having high significance which highlights that ML provides predictive modeling at the expense of statistical inference.

Conclusion

Here we show that the use of ML modeling, specifically using XGBoost allows for prediction of major complications and readmission following aTSA from state-wide claims-based retrospective data. Our model is well calibrated and superior in performance to a traditional LR model. Based on our relative feature importance, shoulder arthroplasty surgeons should inquire and consider a history of prior implant related complications during pre-operative counseling of patients. While further studies are needed to externally validate this model, we hope that this tool can be a building block for physicians to identify modifiable risk factors and help with preoperative counseling, managing patient expectations, informed consent, shared decision making and potentially be useful for risk-adjustment in reimbursement programs.

Click here for additional data file.

Supplemental material, sj-xlsx-2-sea-10.1177_24715492221075444 for Development of a Machine Learning Algorithm for Prediction of Complications and Unplanned Readmission Following Primary Anatomic Total Shoulder Replacements by Sai K Devana, Akash A Shah, Changhee Lee, Andrew R Jensen, Edward Cheung, Mihaela van der Schaar and Nelson F SooHoo in Journal of Shoulder and Elbow Arthroplasty