The typical textile dyeing process calls for a wide range of operational parameters, and it has always been difficult to pinpoint which of these qualities is the most important in dyeing performance. Consequently, this research used a combined design of experiments and machine learning prediction models’ method to offer a sustainable and beneficial reactive cotton fabric dyeing process. To be more precise, we built a least square support vector regression (LSSVR) model based on Taguchi's statistical orthogonal design (L27) to predict exhaustion percentage (E%), fixation rate (F%), and total fixation efficiency (T%) and color strength (K/S) in the reactive cotton dyeing process. The model's prediction accuracy was assessed using many measures, including root mean square error (RMSE), mean absolute error (MAE), and the coefficient of determination (R2). Principal component regression (PCR), partial least square regression (PLSR), and fuzzy modelling were some of the other types of regression models used to compare results. Our findings reveal that the LSSVR model greatly outperformed competing models in predicting the E%, F%, T%, and K/S. This is shown by the LSSVR model's much smaller RMSE and MAE values. Overall, it provided the highest possible R2 values, which reached 0.9819.