This research focuses on developing an effective stress parameter equation for unsaturated soils using multi-gene genetic programming (MGGP). Utilizing a dataset of 101 test points from various experiments, the study demonstrates that the MGGP-derived equation accurately reflects laboratory results and performs well against existing models. Key findings indicate that air entry suction and net stress have significant impacts on effective stress variations, while the influence of other factors is less pronounced.