Global SHAP analysis showed that the most important input features in predicting bone metastasis were AJCC staging, the use of EGFR inhibitor, age, T‐staging, and lymphovascular invasion. This evidence concerns the gene EGFR and bone metastasis.