Abstract Tractor is the most prominent off-road agricultural machinery that is significant to the global food security. The tractive modelling of tyre–soil interaction and agricultural implement dynamics is a complex phenomenon that require holistic approach. Terramechanics techniques such as empirical, semi-empirical, analytical, and numerical methods such as finite element models and discrete element models have gained traction in tractive performance studies. Some of these approaches are premised on large arrays of variables for modelling tractive performance based on the soil–tyre and tools interactions. In this study, soft computing in R software domain was used to model the tractor tractive performance during ploughing operations on a tropical Alfisol. The research farm at the National Centre for Agricultural Mechanization was used for the field experiment. The experimental design was a nested-factorial under a Randomized Complete Block Design having three replications. The input factors were tractor power size, T, (60, 65, and 70 hp); tyre inflation pressure, P, (83, 124, and 165 kPa); implement configuration, I, (2 and 3 bottoms disc plough); and operational speed, S, (6.31, 7.90, 9.47, 11.05, and 12.63 km/h). Standard procedures were followed to obtain the measured parameters in the field, which were statistically analysed. Correlation analysis and analysis of variance of the measured parameters at 5% significance level were established. Multiple linear regression was used to develop the model, validated using the 10-fold cross-validation method. The results revealed that the evaluated variables have a range of 1.56–7.79 kN, 5.15–27.20%, 9.10–32.00 cm, 4.50–13.94%, 1.31–1.67 g/cm3, 95.89–207.78 kPa, and 98.67–295.56 for draught, wheel slip, depth of cut, moisture content, bulk density, cone index (CI), and shear stress, respectively. A positive correlation exists between the towing force (TF) and the measured variables except for the shear stress and CI. The final developed model has seven variables for predicting TF with a 6.5% error and an average of 0.4735 cross validation root mean square error. The model quality of fit achieved an R Adj 2 = 0.8754 {R}_{\text{Adj}}^{2}=0.8754 which satisfactorily described the response variable. The study provides insights into tractive dynamic systems modelling of machine, tractive medium (soil), and agricultural tools anchored on soft computing approach. Its adoption will assist in quality ploughing operation integrating the variables established in the model.