The resilient modulus (MR) of subgrade, which shows relationship between stress and unit deformation of a pavement systems under traffic loads, is a design parameter of the pavement structure. Although a cyclic triaxial test apparatus can be used to directly determine the MR of the subgrade in the laboratory, utilizing prediction models based on easily obtainable soil parameters, is a more efficient method when taking time and cost considerations into account. A comprehensive laboratory testing program is designed to create MR prediction models using machine learning (ML) algorithms. 70 undisturbed soil samples are subjected to MR tests, as well as physical and engineering soil properties tests (water content, field density, specific gravity, gradation, consistency limits, unconfined compressive strength, swell pressure, swell percentage). Soil samples are drilled from a highway that has been in operation for over five years.First, a linear model like MLR is used in the study. Next, nonlinear regression models like RF, GBM, LightGBM, CatBoost, and XGBoost algorithms are used. Research findings showed that nonlinear regression models outperformed linear regression models in predicting the MR (R2 > 0.85), with the XGBoost algorithm yielding the best accuracy (R2 = 0.90). Apart from the primary effects such as confining pressure (σ3) and deviatoric stress (σd), it was found that unconfined compressive strength (qu), natural water content (wn), and swelling percentage (SR) are significant parameters in the prediction of MR among all parameters.