Rectangular concrete-filled steel tubular (RCFST) columns are widely used in structural engineering due to their excellent load-carrying capacity and ductility. However, existing design equations often yield different design results for the same column properties, leading to uncertainty for engineering designers. Furthermore, basic regression analysis fails to precisely forecast the complicated relation between the column properties and its compressive strength. To overcome these challenges, this study suggests two machine learning (ML) models, including the Gaussian process (GPR) and the extreme gradient boosting model (XGBoost). These models employ a range of input variables, such as the geometric and material properties of RCFST columns, to estimate their strength. The models are trained and evaluated based on two datasets consisting of 958 axially loaded RCFST columns and 405 eccentrically loaded RCFST columns. In addition, a unitless output variable, termed the strength index, is introduced to enhance model performance. From evolution metrics, the GPR model emerged as the most accurate and reliable model, with nearly 99% of specimens with less than 20% error. In addition, the prediction results of ML models were compared with the predictions of two existing standard codes and different ML studies. The results indicated that the developed ML models achieved notable enhancement in prediction accuracy. In addition, the Shapley additive interpretation (SHAP) technique is employed for feature analysis. The feature analysis results reveal that the column length and load end-eccentricity parameters negatively impact compressive strength.