Studies comparing numerous sorption curve models and different error functions are lacking completely for soil-metal adsorption systems. We aimed to fill this gap by studying several isotherm models and error functions on soil-metal systems with different sorption curve types. The combination of fifteen sorption curve models and seven error functions were studied for Cd, Cu, Pb, and Zn in competitive systems in four soils with different geochemical properties. Statistical calculations were carried out to compare the results of the minimizing procedures and the fit of the sorption curve models. Although different sorption models and error functions may provide some variation in fitting the models to the experimental data, these differences are mostly not significant statistically. Several sorption models showed very good performances (Brouers-Sotolongo, Sips, Hill, Langmuir-Freundlich) for varying sorption curve types in the studied soil-metal systems, and further models can be suggested for certain sorption curve types. The ERRSQ error function exhibited the lowest error distribution between the experimental data and predicted sorption curves for almost each studied cases. Consequently, their combined use could be suggested for the study of metal sorption in the studied soils. Besides testing more than one sorption isotherm model and error function combination, evaluating the shape of the sorption curve and excluding non-adsorption processes could be advised for reliable data evaluation in soil-metal sorption system.