The success rate of matrix acidizing in hydraulic fractured sandstone formation is less than 55%, much lower compared to the more than 91% success rate in carbonate formation. The need for alternative approaches to help the success ratio in matrix acidizing is crucial. This paper demonstrates a modeling technique to improve the success ratio of matrix acidizing in a hydraulic fractured sandstone formation. Supervised machine learning with 4 models of a neural network, logistic regression, tree, and random forest was selected to predict the successfulness of matrix acidizing in hydraulic fracturing. In parallel, multivariate analysis of principal component regression and partial least square regression approach were utilized to predict the oil gain of the job. For qualitative prediction, the results showed that the random forest was the best model to predict the successfulness of the job with the area under the curve (AUC) of 0.68 and precision of 0.73 in the training model with 70% of the data. Subsequently, the validation test with the rest of the data (30% data) gave 0.51 AUC and 61% precision. For quantitative prediction, the net oil gain was evaluated by using principal component regression (PCR) and partial least square regression (PLS-R). The PCR and PLS-R model gave a coefficient of determination (Rsquare) of 0.22 and 0.35, respectively. The p-value of PLS-R was 0.047 (95% confidence interval) which indicates that the model is significant. The results of this work demonstrate the potential application of supervised machine learning, principal component regression, and partial least square regression to improve candidate selection of oil wells for matrix acidizing especially in hydraulic fractured wells with limited design data.