Abstract The objective of this study was to explore the potential of Machine Learning (ML) algorithms to increase the accuracy of predicting individual days in the herd (an indicator of stayability) using reproductive records and genomic information. A total of 6943 cows from 3 herds with reproductive performance were included in the study from which 696 cows had genomic information (genotyped using Illumina Bovine 50k SNP BeadChip). Different libraries based on R and Python were used to test various ML models including Lazy Predict, Scikit-learn, PyCaret, and H2O Flow. Genomic information was subjected to quality control by removing SNPs with an allele frequency less than 0.05 or with a call rate lower than 0.95. A total of 42,689 SNP remained for further analysis and accounted for 11% of phenotypic variation (heritability of 0.11±0.02) in DIH. Different numbers of SNPs (500 SNPs, 1K, 5K, 10K, and 15K) were selected based on their contribution to phenotypic variation from GWAS and were included in the models. Model performance measures, such as mean absolute error (MAE) and mean square of error (MSE), worsened with increased SNPs in the model. Bayesian Ridge algorithm using 500 top SNPs contributed to the phenotypic variance, had the best performance to predict DIH with MAE of 612.6 and R2 of 0.52 in the training population using PyCaret program. When BWT and WWT were added to the model, in addition to SNPs, little change was observed in the model’s performance. Overall, we concluded that ML models had better performance compared to the conventional modeling approach and genomic analysis; CatBoost model had 55% lower mean square of error compared to the simple linear regression (734650 vs 1637410). The results suggest that ML tools have the potential to improve the accuracy of predicting DIH compared to simple linear regression and conventional genomic analysis.