Background: Detection of axillary metastases in breast cancer is critical for treatment options and prognosis. The aim of this study is to investigate the value of radiomic features obtained from short tau inversion recovery (STIR) sequences in magnetic resonance imaging (MRI) of primary tumor in breast cancer in predicting axillary lymph node metastasis (ALNM). Methods: Lesions of 165 patients with a mean age of 51.12 ±11 (range 28-82) with newly diagnosed invasive breast cancer who underwent breast MRI before treatment were manually segmented from STIR sequences in the 3D Slicer program in all sections. Machine learning (ML) analysis was performed using the extracted 851 features Python 3.11, Pycaret library program. Datasets were randomly divided into train (123, 80%) and independent test (63, 20%) datasets. The performances of ML algorithms were compared with area under curve (AUC), accuracy, recall, presicion and F1 scores. Results: Accuracy and AUC in the training set were in the range of 57 %-86 % and 0.50-0.95, respectively. The best model in the training set was the catBoost classifier with an AUC of 0.95 and 84% accuracy. The AUC, accuracy, recall, precision values and F1 score of the CatBoost classifier on the test set were 0.92, 84 %, 89 %, 85 %, 86 %, respectively. Conclusion: Radiomic features obtained from primary tumors on STIR sequences have the potential to predict ALNM in invasive breast cancer.