Least squares support vector machines (LS-SVM) were developed for the analysis and prediction of the relationship between the cutting conditions and the corresponding fractal parameters of machined surfaces in face milling operation. These models can help manufacturers to determine the appropriate cutting conditions, in order to achieve specific surface roughness profile geometry, and hence achieve the desired tribological performance (e.g. friction and wear) between the contacting surfaces. The input parameters of the LS-SVM are the cutting parameters: rotational speed, feed, depth of milling. The output parameters of the LS-SVM are the corresponding calculated fractal parameters: fractal dimension D and vertical scaling parameter G. The LS-SVM were utilized successfully for training and predicting the fractal parameters D and G in face milling operations. Moreover, Weierstrass-Mandelbrot（W–M ）fractal function was integrated with the LS-SVM in order to generate an artificially fractal predicted profiles at different milling conditions. The predicted profiles were found statistically similar to the actual measured profiles of test specimens and there is a relationship between the scale-independent fractal coefficients(D and G).