Background <p>Total joint arthroplasty (TJA) complications necessitate the development of accurate risk prediction models; however, interpretability in machine learning remains a challenge. While Shapley Additive Explanations (SHAP) offers insights at the individual level, partial dependence plots (PDPs) may provide a better understanding at the population level for developing clinical guidelines. This study compared PDPs and SHAP in explaining machine learning-based 30-day complication risk prediction following TJA.</p> Methods <p>We conducted a retrospective cohort study using the American College of Surgeons National Surgical Quality Improvement Program (NSQIP) database (2019–2023), including 517,826 primary TJA cases. Binary classification models (Random Forest, Gradient Boosting) predicted composite 30-day complications based on 20 clinical predictors. A comprehensive interpretability analysis employed directional concordance validation between PDP and SHAP, permutation importance thresholding (5% relative influence), followed by one- and two-dimensional partial dependence analyses with explicit interaction modeling.</p> Results <p>The cohort comprised 517,826 primary TJA procedures with a complication rate of 6.67%. The baseline Random Forest model achieved test AUC = 0.678. Directional concordance analysis demonstrated 97.8% weighted agreement between PDP trends and SHAP attributions, validating methodological comparison. Threshold analysis identified seven significant features, with interaction effects accounting for 49.9% of total model influence (71.9% among top features). PDPs showed actionable dose–response relationships, including critical thresholds for preoperative hematocrit (&lt; 38%), operative time (&gt; 120&#xa0;min), and complementary interactions, such as age × ASA classification (19.1% importance), operative time × ASA classification (10.1%), and hematocrit × diabetes (6.4%). Comparative patient analysis demonstrated that while SHAP quantified individual contributions, only PDPs provided population thresholds directly translatable to institutional protocols.</p> Conclusion <p>PDPs appear more methodologically appropriate than SHAP for population-level clinical guideline development, offering actionable dose–response relationships and population risk thresholds that SHAP’s individualized attribution framework cannot provide. The dominance of interaction effects among the most influential predictors validates that PDPs accurately capture complementary relationships while presenting them in a format directly applicable to evidence-based perioperative protocols and institutional quality improvement initiatives.</p> <p><MediaObject ID="MOESM3"> <VideoObject FileRef="MediaObjects/42836_2025_360_MOESM3_ESM.mp4" VideoID="Dcq8D6aHjyr-cxqU1UHHwW"> <Caption Language="En" xml:lang="en"> <CaptionContent> <p>Video Abstract</p> </CaptionContent> </Caption> </VideoObject> </MediaObject></p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Evaluating the methodological suitability of partial dependence plots and Shapley additive explanations for population-level interpretation of machine learning models in total joint arthroplasty

  • Kole Joachim,
  • Othneil Sparks,
  • Amanda Perrotta,
  • Adrian Lin,
  • Brandon Gettleman,
  • Christopher Hamad,
  • Sumin Jeong,
  • Ezekiel Dingle,
  • Alexandra Stavrakis,
  • Alexander B. Christ

摘要

Background

Total joint arthroplasty (TJA) complications necessitate the development of accurate risk prediction models; however, interpretability in machine learning remains a challenge. While Shapley Additive Explanations (SHAP) offers insights at the individual level, partial dependence plots (PDPs) may provide a better understanding at the population level for developing clinical guidelines. This study compared PDPs and SHAP in explaining machine learning-based 30-day complication risk prediction following TJA.

Methods

We conducted a retrospective cohort study using the American College of Surgeons National Surgical Quality Improvement Program (NSQIP) database (2019–2023), including 517,826 primary TJA cases. Binary classification models (Random Forest, Gradient Boosting) predicted composite 30-day complications based on 20 clinical predictors. A comprehensive interpretability analysis employed directional concordance validation between PDP and SHAP, permutation importance thresholding (5% relative influence), followed by one- and two-dimensional partial dependence analyses with explicit interaction modeling.

Results

The cohort comprised 517,826 primary TJA procedures with a complication rate of 6.67%. The baseline Random Forest model achieved test AUC = 0.678. Directional concordance analysis demonstrated 97.8% weighted agreement between PDP trends and SHAP attributions, validating methodological comparison. Threshold analysis identified seven significant features, with interaction effects accounting for 49.9% of total model influence (71.9% among top features). PDPs showed actionable dose–response relationships, including critical thresholds for preoperative hematocrit (< 38%), operative time (> 120 min), and complementary interactions, such as age × ASA classification (19.1% importance), operative time × ASA classification (10.1%), and hematocrit × diabetes (6.4%). Comparative patient analysis demonstrated that while SHAP quantified individual contributions, only PDPs provided population thresholds directly translatable to institutional protocols.

Conclusion

PDPs appear more methodologically appropriate than SHAP for population-level clinical guideline development, offering actionable dose–response relationships and population risk thresholds that SHAP’s individualized attribution framework cannot provide. The dominance of interaction effects among the most influential predictors validates that PDPs accurately capture complementary relationships while presenting them in a format directly applicable to evidence-based perioperative protocols and institutional quality improvement initiatives.

Video Abstract