Surface roughness is an important indicator of the quality of the machined surface. One of the methods that can be applied to improve surface roughness is ball burnishing. Ball burnishing is a finishing process in which a ball is rolled over the workpiece surface. Defining adequate input variables of the ball burnishing process to ensure obtaining the required surface roughness is a typical problem in scientific research. This paper presents the results of experiments to investigate ball burnishing of AISI 4130 alloy steel with a high-stiffness tool and a ceramic ball. The experiments were conducted following a randomized full factorial design for different levels of input variables. The input variables included the initial arithmetic mean roughness (the initial surface roughness), the depth of ball penetration, the burnishing feed, and the burnishing ball diameter, while the output variable was the arithmetic mean roughness after ball burnishing (the final surface roughness). The surface roughness modeling was performed based on the experimental results, using regression analysis (RA), artificial neural network (ANN), and support vector regression (SVR). The regression model displayed large prediction errors at low surface roughness values (below 1 μm), but it proved to be reliable for higher roughness values. The ANN and SVR models have excellently predicted roughness across a range of input variables. Mean percentage error (MPE) during the experimental training research was 29.727%, 0.995%, and 1.592%, and MPE in the confirmation experiments was 34.534%, 1.559%, and 2.164%, for RA, ANN, and SVR, respectively. Based on the obtained MPEs, it can be concluded that the application of ANN and SVR was adequate for modeling the ball burnishing process and prediction of the roughness of the treated surface in terms of the possibility of practical application in real industrial conditions.