Application of KRR, K-NN and GPR Algorithms for Predicting the Soaked CBR of Fine-Grained Plastic Soils

California bearing ratio (CBR) test is one of the comprehensive tests used for the last few decades to design the pavement thickness of roadways, railways and airport runways. Laboratory-performed CBR test is considerably rigorous and time-taking. In a quest for an alternative solution, this study utilizes novel computational approaches, including the kernel ridges regression, K-nearest neighbor and Gaussian process regression (GPR), to predict the soaked CBR value of soils. A vast quantity of 1011 in situ soil samples were collected from an ongoing highway project work site. Two data divisional approaches, i.e., K-Fold and fuzzy c-means (FCM) clustering, were used to separate the dataset into training and testing subsets. Apart from the numerous statistical performance measurement indices, ranking and overfitting analysis were used to identify the best-fitted CBR prediction model. Additionally, the literature models were also tried to validate through present study datasets. From the results of Pearson’s correlation analysis, Sand, Fine Content, Plastic Limit, Plasticity Index, Maximum Dry Density and Optimum Moisture Content were found to be most influencing input parameters in developing the soaked CBR of fine-grained plastic soils. Experimental results also establish the proficiency of the GPR model developed through FCM and K-Fold data division approaches. The K-Fold data division approach was found to be helpful in removing the overfitting of the models. Furthermore, the predictive ability of any model is considerably influenced by the geological location of the soils/materials used for the model development.

GPR K GPR model at K-Fold data division approach IOA Index of agreement IOS Index of scatter I P

Introduction
Road transportation network facilitates transferring goods from one place to another and door-to-door services for passengers throughout the world.As of this, the road transport infrastructures majorly govern the economy of the country.In this aspect, many new expressways and green highways are being constructed in India by the Ministry of Road Transport and Highways (MoRTH) department through various infrastructure development plans.The pavement thickness design construction of these roads is based on the strength of the material used in the subgrade and subbase layer.Therefore, highway engineers always desire that the material used in the subgrade layer should fulfill some of the engineering and technical properties such as swell criteria, plasticity properties, soil settlement conditions, subgrade reaction, bearing capacity etc.A method for recognizing the strength of such layers is of utmost requisite in highway engineering.
In general, the California bearing ratio (CBR) test is espoused to measure the stiffness modulus and the shear strength of subgrade material [1,2] which may be performed on either re-compacted samples in the laboratory or undisturbed samples cut from the field or in situ surface of subgrade formation [3].The test is an indirect measure that compares the strength of subgrade material (at known density and moisture content) to standard crushed rock material [4,5].Both laboratory and in situ tests are based on the principle of penetrating a standard dimension plunger into a soil specimen at a deformation rate of 1.25 mm/min.The laboratory and field engineers always encounter several difficulties in obtaining the CBR value in the laboratory.Laboratory-soaked CBR test requires a large amount of materials (almost 6 kg), more effort to prepare the test specimen, and lastly, 96 h of the soaking period to simulate the field conditions.Consequently, all those activities make the CBR test more tedious, laborious, and time-consuming.Additionally, if the properties of soil change for each small stretch of highway, then preserving such a huge quantity of soil and conducting the CBR test in the laboratory is laborious and time-consuming.Laboratories are also often packed due to the long queue of materials testing, which causes a delay in testing as well as the testing reports, ultimately the design of construction projects.Furthermore, the test method includes the material transportation cost (from construction site to testing laboratory), testing charge, and finally, the dumping of tested materials, which became more exhausted and increased the final cost of the projects.
Owing to the aforementioned problems, many researchers considered that CBR needs to be replaced either partially or entirely.Although not a fundamental material property, it has a long history in pavement design, and it is reasonably correlated with the index and engineering properties of soil by several investigators in the past.To the author's knowledge, the first fame in predicting the CBR value was earned by Kleyn [6].Earlier, he attempted to address the discrepancy in the CBR test and later prepared a chart based on a nest of straight lines that relate CBR to PI and grading module for over 1000 soaked CBR tests obtained from road and airport work throughout central and southern Africa.Black [7] suggested that the relationship between CBR and ultimate bearing capacity depends on the type of soil and compaction method, i.e., static or dynamic.Agarwal and Ghanekar [8] tried to generate the correlation equation through statistical analysis between CBR and Atterberg limits for 48 soil samples collected from different parts of India.However, they could not find any significant correlation between these parameters.But when LL and OMC were added, they observed an improved correlation with adequate accuracy for the preliminary identification of materials.National Cooperative Highway Research Program [9] attempted to develop the correlation equation for CBR from the index properties for clean and coarse-grained soil.Kin [10] tried to develop the correlation equation for the CBR value of fine-grained and coarse-grained soil through gradational properties.Taskiran [11] attempted to establish the correlation for 151 CBR test data of fine-grained soils, taken from 354 test samples, by ANN and gene expression programming (GEP) methods.Both techniques were found to exhibit promising results.Using 124 datasets, Yildirim and Gunaydin [12] studied the estimation of CBR by regression and ANN approach.They observed that the ANN technique is better than the regression analysis.Erzin and Turkoz [13] tried to predict the CBR value of Aegean sand from the results of mineralogical properties through ANN and regression approach.They also found that ANN is superior to the regression technique.Farias, Araujo [14] used the local polynomial regression (LPR) and radial basis network (RBN) techniques for developing the predictive equations for the CBR of soil samples.Using 207 CBR test results of granular soil, Taha, Gabr [15] observed that the correlation obtained through ANN is of excellent accuracy and lower bias than the regression analysis.A comparative study conducted by Tenpe and Patel [16] for 389 datasets collected from City and Industrial Development Corporation, Maharashtra state in India, reveals that ANN and GEP are efficient in predicting the CBR value.Later, in another study, Tenpe and Patel [17] found that SVM can better predict the CBR value than GEP.Recently, Bardhan, Samui [18] attempted to predict the soaked CBR value of 312 soil datasets through a particle swarm optimization (PSO) algorithm with adaptive and time-varying acceleration coefficients.The comparative analysis of various extreme learning machine (ELM) based adaptive neuro swarm intelligence (ANSI) such as ELM coupled-modified PSO (ELM-MPSO), ELM coupledtime-varying acceleration coefficients PSO (ELM-TPSO) and ELM coupled-improved PSO (ELM-IPSO) reveals that the modified and improved version of PSO has high accuracy at early iterations than the standard PSO.In another investigation, Bardhan, Gokceoglu [19] observed that multivariate adaptive regression splines with piecewise linear (MARS-L) demonstrate a higher accuracy in predicting the soaked CBR as compared to MARS with piecewise-cubic (MARS-C), Gaussian process regression and genetic programming.Hassan, Alshameri [20] attempted to predict the CBR value of plastic fine-grained soil from their index properties and compaction parameters through multi-linear regression analysis (MLR).The study was conducted for the standard proctor compaction energy level, whereas the engineers always prefer the modified proctor compactive energy level to construct the highways and expressways.It is observed from the above literature investigations (also shortened in Table 1) that several artificial intelligence (AI)-based models were used to predict the soaked CBR value which demonstrates the precision from 80 to 100% (R 2 values 0.8 to 1.0).However, there are still some advanced computation approaches which have proven their competency in solving many problems of civil engineering.The literature studies also omitted the investigation of statistical analysis over the obtained results of the model.The deep insight view of literature studies reveals that the range of geotechnical parameters and quantity of dataset are limited.Using a large amount of dataset is always considered to be much worthwhile from generalization point of view [18,19,[21][22][23][24].

Research Significance and Contributions
The main contribution of this study is to develop an efficient model for predicting one of the challenging real-world problems of highway engineers, i.e., estimation of soil California bearing ratio value.For this purpose, Kernel Ridge Regression (KRR), K-Nearest Neighbor (K-NN) and Gaussian Process Regression (GPR) algorithms were adopted.This study utilizes 1011 in situ samples of fine-grained plastic soils with an extensive range of index and engineering properties.Numerous geotechnical parameters were extracted from the laboratory experiments conducted through Bureau of Indian Standards (BIS) specifications.Two data divisional approaches, viz.K-fold and FCM, were adopted to investigate the influence of the training data features on the predictive ability of the developed model.Furthermore, the influence of employed machine learning algorithms as well as the data division approaches was investigated on the predictive ability of the model.Lastly, the literature models were attempted to validate through the present study datasets.

Machine Learning Algorithms and Statistical Assessment Indices
The term "Machine Learning (ML)" is a subfield/type of Artificial Intelligence (AI) and is referred to as predictive analytics or predictive modeling.ML is the development of computer systems that can learn and adapt without following explicit instructions by using algorithms and statistical models to analyze and draw inferences from patterns in data.In the recent past, numerous ML algorithms have been adopted by several researchers for solving many significant engineering problems .This section briefly introduces the most prominent ML algorithms used to develop the model for predicting the soaked CBR value of fine-grained plastic soils and several indices to measure their performances.

Kernel Ridge Regression (KRR)
KRR is the nonlinear regression approach that is based on the "kernel trick" in which datasets are nonlinearly transformed into some high-dimensional (or even infinite-dimensional) feature space determined by the kernel functions satisfying Mercer's theorem [50][51][52].Consider a TR set of (x  [51,53,54].The OLS minimizes the squared loss function: where .indicates the L 2 norm.In order to control the tradeoff between bias and variance of the estimate, a shrinkage or ridge parameter λ is added to the above expression which is represented below: Using the "kernel trick," the KRR extends the linear regression into nonlinear and high-dimensional space.The data x i in X is replaced with the feature vectors: x i → F F (x i )induced by the kernel where K i j k x i , x j F (x i ) F x j .Therefore, the predicted class label of a new example x is now represented as: where k (k 1 , k 2 , . . . . . ., k N ) T , k N x N .xand n 1, 2, . . . . . ., N .
In KRR, the kernel function is used to increase the computational power by mapping the data into a high-dimensional feature space which makes the data linear separable, and consequently, increases the stability, accuracy and generalization for both regression and classification problems.In general, the kernel functions are categorized into local and global functions.The local function influences the kernel function values within the vicinity like the Gaussian kernel function, whereas the global function allows the data points which are far from each other like the polynomial kernel function.The other kernel functions can also be employed, such as RBF, sigmoid and Laplacian.

K-Nearest Neighbor (K-NN)
K-NN is a nonparametric supervised ML algorithm that uses the k-number of most similar outputs from the TR dataset [55][56][57][58].The K-NN algorithm is also known as a lazy learner algorithm, because, instead of immediate prediction, it stores the dataset and categorizes it based on its similarities and then approaches an action on the dataset.K-NN predicts the new records/data for the regression and classification-type problems based on their Euclidean distances, estimated mean, median or model output variable [59][60][61][62].
A classic K-NN algorithm for a regression problem can follow the below-given steps [63]: Step 1: Consider a vector X with m independent variables, in this study X {S, FC, P L, P I , M D Dand O MC}, as predictors and Y dependent variable, i.e., Y {C B R}.
Step 2: Consider a TR set consisting of a dataset with t vectors of X n {S 1n , FC 1n , P L 1n , P I 1n , M D D 1n and O MC 1n } and a dependent variable of Y n {C B R 1n } associated with each vector in X n .
Step 3: The distance of the predictor vector with each of the n training vectors is calculated.The Euclidean distance (E d ) is most commonly used to estimate (see Eq. 4) which of the k output in the TR dataset is similar to the new input.However, other distances such as the Manhattan distance, which computes the distance between real vectors using the sum of absolute distance, Hamming distance, estimating the distance between binary vectors, and Minkowski distance, a generalization of the Euclidean and Manhattan distances [58,64] could also be adopted according to the availability of the datasets.
Step 4: Select the k training vector with the least distance to the predictor vector.
Step 5: Calculate the kernel function (using an equation) for each of the k-selected training vectors Step 6: The predicted dependent variable is calculated as where Y nk is the dependent variable of each selected neighbor.

Gaussian Process Regression (GPR)
GPR is a probabilistic, nonparametric Bayesian approach for generalizing the nonlinear and complex problems related to regression and classification-type data sets [65,66].Many supervised ML algorithms learn exact values from the dataset, whereas GPR infers a probability distribution of all admissible functions that could reasonably fit the data space regarding the problems [67].GPR is very efficient to handle nonlinear data due to the use of kernel functions.
A GPR model can make predictions by incorporating the prior knowledge through covariance (kernel) functions and provide uncertainty measures over predictions [66,[68][69][70].
The algorithm has recently received huge attention from researchers as of having ability to solve many complex engineering problems related to various disciplines [66,67,69,[71][72][73][74].
For a given TR data set of S {((x i , y i )|i 1, 2, . . . . . . . . ., n)}, the input data X ∈ R S×n is called the designed matrix and Y ∈ R n is the vector of the desired output.The main assumption of GPR is that the output can be estimated as follows: where ε signifies the noise term.GPR assumes that ε follows the Gaussian distribution with a mean of 0 and the variance of σ 2 n : In GPR methodology, the n observations in the data set of interest y {y 1 , y 2 , . . . . . . . . ., y n } are considered as a single point sampled from a multivariate Gaussian distribution.Moreover, it can be assumed that this Gaussian distribution has a mean of zeros.The covariance function K (x, x ) dictates the relation of one observation to another observation.
For a given TR data set, the ultimate goal of the learning process is to predict the output value y * of a new queried input pattern.To achieve such a goal, it is necessary to establish three covariance matrices as follows: ( 1 0 ) Due to the assumption that the data is sampled from multivariate Gaussian distribution, we have the following expression: Since, (y * |y) is proved to be generated from multivariate Gaussian distribution the mean and the variance of the estimated mean and variance of the predicted output are given as:

Hyperparameters Tuning Using a Grid Search
Several ML algorithms contain a set of parameters that control many aspects of the algorithm.These parameters are termed as hyperparameters that are fixed before the learning process of the algorithm.These hyperparameters are adjusted to enhance the quality of the model as well as the ability to correctly predict the unseen dataset.In this study, hyperparameters of the above-used algorithms are optimized through a cross-validation (CV)-based grid search method.In the grid search method, the dataset is divided into training sets and test sets using K-Fold cross-validation (CV) [75].For CV calculation, the grid points are assigned.Out of a total five number of folds, onefold is served as an independent testing set and the remaining k-1 (four) folds are designated as training sets.The grid search method is significantly superior to the random search method where only a few combinations are searched.The most significant advantage of this approach is that it provides better generalization performance for the respective model.Initially, some manifestly unreasonable values of the hyperparameters were bypassed.After several trial-and-error approaches, the expected values for each of the hyperparameters were passed through the grid search method.Based on the results of best fit, the values obtained for the selected hyperparameters of KRR, K-NN and GPR algorithms are given in Tables 2, 3 and 4.

Statistical Performance Measurement Indices
The precision of all the models was assessed through several statistical performance indicators.The widely used performance measurement indicators are coefficient of determination (R 2 ), adjusted R 2 (adj.R 2 ), coefficient of correlation (R), mean absolute error (MAE), mean absolute percentage error (MAPE), root-mean-square error (RMSE), variance accounted for (VAF), performance index (I P ), Willmott's index of agreement (IOA), index of scattering (IOS), a20index and performance strength (S P ) [17,[76][77][78][79][80][81].The mathematical expression for these indicators is given in Eqs.(15) to (26) along with their ranges.In order to select an efficient prediction model espousing many more statistical indicators might be useful for assessing the performance of predictive models in terms of error and trend point of view [19].Using the least amount of parameters might be challenging, especially when comparing the performance of two or more models in an instant.Often, the models show not much difference in their statistical indices value.In that particular situation considering many more statistical performance indices might be helpful in selecting the best-fitted model.

Parameters
Ideal value M AP E(%) V ar(y i (a)) × 100 100 ( 21) where y i (a) actual value (laboratory-obtained value); y i ( p) predicted value (value obtained through the developed model); y i (a) mean of actual value; y i ( p) mean of predicted value; N number of observations; P number of input parameters used to develop the model; n20 number of observations lies within error range of ± 20%; a20 percentage of observations having error ≤ 20% (Fig. 1).

Data Collection and Geographical Location
The Ministry of Road Transport and Highways (MoRTH) has decided to implement the Engineering Procurement and Construction (EPC) mode for constructing the National Highways (NH) in India.The Twelfth 5-Year Plan envisions the construction of 20,000 km of 4-lane National Highways projects through EPC mode in the various state of India.In this study, the Varanasi Gorakhpur section of NH-29 belongs to the state of Uttar Pradesh (UP), India, was selected as the study area [42].The geographical location of the study area is shown in Fig. 2.

Laboratory Experiments
From the aforementioned project worksite, a total of 1011 soil samples were collected and brought to the laboratory for experimental investigations.The laboratory tests were conducted as per BIS specifications.The test method includes IS 2720 (Part 4) [82], IS 2720 (Part 5) [83], IS 2720 (Part 8) [84] and IS 2720 (Part 16) [85] for the grain size distribution, Atterberg limits, modified Proctor compaction parameters and soaked CBR of fine-grained soil, respectively.Through these laboratory tests, numerous geotechnical parameters were collected such as gravel content (G), sand content (S), silt and clay content termed as fine content (FC), liquid limit (LL), plastic limit (PL), plasticity index (PI), maximum dry density (MDD), optimum moisture content (OMC) and CBR value of fine-grained soil.The laboratory-obtained database of fine-grained soil was further classified into various soil groups using the USCS soil classification system.Figure 3 presents the histogram plot for different soil groups of finegrained soil.The vertical column represents the amount of particular soil groups present in the fine-grained soil database.

Statistical Visualization and Correlation Analysis
Table 5 presents the descriptive statistic values for all the finegrained soil parameters.As seen from Table 5, the obtained database covers an extensive range of CBR values from 1.0 to 13.20.The gravel content varies from 0 to 28%, sand ranges between 2 and 49%, and fine content exists in the range of 50 to 96%.Similarly, from a plasticity point of view, the selected database shows the liquid limit from 24 to 85%, plastic limit from 11 to 50% and plasticity index from 1 to 39%.MDD and OMC of fine-grained soil range from 1.455 g/cc to 1.959 g/cc and 9% to 30%, respectively.
Pearson's correlation (R) is one of the commonly used measures of association between the parameters.The range of R varies from − 1 to 1, where ± 1 indicates the strong association between the parameters, and a value of 0 (zero) illustrates no relationship between the parameters.The positive or negative sign specifies the respective increase or decrease in the associated parameters simultaneously.The correlation matrix obtained for all the geotechnical parameters is presented in Fig. 4. As observed from Fig. 4, CBR is positively correlated with S and MDD, whereas negatively with FC, LL, PL, PI and OMC.The final selected input parameters for developing the CBR prediction model are S, FC, PL, PI, MDD and OMC.

Data Divisional Approaches
Data division is the process of separating the complete dataset into the TR and TS subsets.In this study, about 80% of Fig. 4 Correlation matrix for all the geotechnical parameters the entire dataset was considered to train the model and the remaining 20% was kept for testing the model.The basic problem with machine learning modeling is that we are unknown to the fact that how well a model performs or will perform until it is tested on the unseen dataset.One can build a perfect model on the TR data with 100% accuracy or 0 error, but it may fail to generalize for unseen data.It is not a good model as it over-fits the TR data.Machine learning is all about generalization, meaning that the model's performance can only be measured with data points that have never been used during the training process.To overcome this problem, two data divisional approaches were adopted, discussed below:

K-Fold Division Approach
Initially, shuffle the dataset randomly and split it into Knumber of folds (as seen in Fig. 5).Once the dataset is separated, the first fold is used as the TS dataset and the remaining k-1 folds are used for training purposes.The model with specific hyperparameters is trained with TR data (K-1 folds) and TS data as one fold.The performance of the model is recorded.The above steps are repeated until each k-fold got used for testing purposes.Using fivefold, the complete datasets were divided into TR and TS sets.The descriptive statistics value obtained for the TR and TS datasets is shown in Tables 6 and 7, respectively.

Fuzzy C-Means (FCM) Division Approach
Clustering is the process of separating the dataset into a number of groups.Each group represents the observations that are homogeneous to each other and the objects that are dissimilar to each other are clustered into different groups.In fuzzy clustering, each observation can belong to more than one cluster based on its membership value.The total membership value for any observation distributed over the entire cluster is 1.0.Brief information about this approach is given in Shi [86], Shahin, Maier [87] and Das [88].Initially, two clusters were taken and the silhouette value was estimated.The number of a cluster was increased gradually and silhouette values for each of the corresponding clusters were calculated.The silhouette score obtained corresponding to the number of clusters is depicted in Fig. 6.It can be clearly observed from Fig. 6 that the maximum silhouette score was obtained when two clusters were used  for the analysis.The dataset obtained in the first and second clusters is represented as C 1 and C 2 , respectively.The C 1 dataset was separated into TR and TS sets through the Kfold approach (discussed in Sect.8 and  9.

Statistical Performance of Developed Models
Table 10 illustrates the performance of the trained KRR, K-NN and GPR models in terms of several performance measurement indices.In Table 10, KRR K and KRR F indicate the KRR model developed through K-fold and FCM approaches, respectively.Similarly, models developed for K-NN and GPR algorithms are represented as K-NN K , K-NN F and GPR K , GPR F , respectively.It is observed from Table 10 that the R 2 values obtained for the KRR algorithm in The performance of developed models on the TS dataset is presented in Table 11.The R 2 value obtained in the TS dataset of the K-Fold approach is 0.680, 0.706 and 0.758 for KRR, K-NN and GPR models, respectively, whereas for KRR F , K-NN F and GPR F models are 0.407, 0.645 and 0.700, respectively.It is clearly observed from the deep closest view of all the statistical parameters in training (refer to Table 10) and testing (refer to Table 11) dataset that the GPR algorithm achieved the highest prediction in both K-Fold and FCM approaches which is followed by K-NN and KRR algorithms.

Visual Interpretation of Developed Models
Visual interpretation facilitates the viewer to find the insight features from the model which is represented in a graphical form such as a scatter plot, error plot and regression error characteristics curve, etc.

Trend and Error Plot for the Developed Models
Figures 7 and 8 present the actual versus predicted CBR value of the TR and TS datasets, respectively, for the KRR, K-NN and GPR models at K-Fold and FCM approaches.It is observed from the scatter plot results that the maximum number of observations follows a specific trend along with the line of equality.However, the closeness of data points toward the line of equality in the TR dataset (see Fig. 7) is maximum for GPR models, and an acceptable inclination is obtained for KRR and K-NN models.Moreover, the utmost precision in the GPR model is obtained when the dataset is trained through the FCM data division approach followed by the K-Fold approach.Similarly, the results obtained for the TS dataset (see Fig. 8) also reveal that GPR models are more efficient in predicting the unseen dataset, followed by K-NN and KRR algorithms.
The error distribution and bar plots for the TR and TS dataset of KRR, K-NN and GPR models at K-Fold and FCM data division approaches are shown in Figs. 9 and 11, respectively.The center horizontal line in the error distribution plot represents the zero error line, the datasets existing on that line have zero error, i.e., the difference in the actual and predicted CBR value is zero.At the same time, the upper and lower line specifies the + 20% and -20%, respectively, error or variation band.It is observed from Figs. 9 and 11 that some datasets are below the zero error line (displays negative error) and some are above the zero error line (shows positive error).This random pattern of error indicates that the developed models have a decent fit for the dataset.Furthermore, the existence of a dataset within ± 20% variation seems to be the maximum for the GPR F model (refer to Fig. 11).This can also be confirmed from the comparative analysis of Figs. 10 and 12, showing the error frequency plot for the TR and TS datasets of KRR, K-NN and GPR models at K-Fold and FCM approaches, respectively.As seen from Fig. 12, the GPR F model can predict almost 96% and 90% observations of the TR and TS dataset, respectively, within ± 10%, whereas 100% and 99% within ± 20% variations which are substantially higher than other models (Figs.11, 12).

Regression Error Characteristics (REC) curve
In regression problems, REC curves are equivalents to the receiver operating characteristics (ROC) curves in classification problems.The X-axis of the REC curve plot demonstrates the error tolerance, whereas the Y-axis represents the accuracy in terms of the percentage of points predicted within the tolerance [79,80,89].An ideal model's curve should pass through the upper left corner and therefore should have an area under the curve (AUC) value is 1.This means that the model can perfectly discriminate between all the positive and the negative class points.In general, an AUC of 0.5 suggests no discrimination, 0.7 to 0.8 is considered acceptable, 0.8 to 0.9 is deemed to be excellent, and more than 0.9 is considered outstanding.Figures 13, 14, 15 and 16 depict the REC curve obtained corresponding to the TR and TS dataset of KRR, K-NN and GPR models at K-Fold and FCM data division approaches.As seen from these figures, the AUC value obtained for all the models in TR and TS datasets is higher than 0.9, which means the developed models outperform very well and are stated to be reliable in predicting the soaked CBR value of fine-grained plastic soils.Furthermore, it is observed from the comparative analysis of both models in the training and testing set that the REC curve for the GPR F model exists more closely to the upper left corner compared to the other models as well as the AUC value achieved for the GPR F model is also higher than those models.

Accuracy Analysis
The accuracy analysis is a novel assessment used to evaluate the efficiency of the models.The analysis demonstrates the accuracy (%) of a model, which is obtained through the comparative analysis of the values obtained for different performance measurement parameters to their ideal values (as shown in Sect.2.2) using Eqs.( 27) and (28).
Fig. 9 Error distribution plot for the TR and TS dataset of KRR, K-NN and GPR models at K-Fold approach where A e and A t denote the error and trend-measuring performance parameters.m e and m t indicates the measured values of the error and trend-measuring performance parameters.The performance measurement parameters MAE, MAPE, RMSE and IOS belong to the class of error, whereas R 2 , Adj.R 2 , R, VAF, I P , IOA and a20-index belong to the trend.i t represents the ideal value of the respective trend parameters.Tables 12 and 13 present the accuracy of all the developed models in the TR and TS dataset, respectively.As seen from Tables 12 and 13, the accuracy of the GPR F model throughout numerous statistical performance measurement indices is significantly higher than other models.

Selection of Best-Fitted CBR Prediction Model
Numerous statistical performance indices were adopted to evaluate the model performances and the conclusions can easily be drawn by comparing their values.However, the situation becomes more complicated when the model demonstrates adequate accuracy in the TR dataset and failed to reach a good amount of accuracy in the TS dataset as well.Moreover, it also becomes complicated when the value of different statistical performance indices describes their own best models.In that particular situation, the overfitting analysis and ranking analysis might be useful for identifying the bestfitted model as they envisage the overall analysis of all the parameters of a particular model.

Ranking Analysis (RA)
According to RA, used by many researchers in the past [79,80,89,90], a maximum score of s (equal to the total number of corresponding models) is assigned to the model having the highest value in particular performance indices, minimum to the model with the lowest value and the score to the other intermediate models are assigned either in the ascending or descending order.Table 14 presents the RA results obtained for the KRR, K-NN and GPR models at K-Fold and FCM approaches.

Overfitting Ratio (OR)
OR is the ratio of the RMSE of the TS dataset to the TR dataset as shown in Eq. ( 29).According to this formula, a model with a lesser OR value or close to one is considered Fig. 12 Error frequency plot the TR and TS dataset of KRR, K-NN and GPR models at FCM approach less prone to overfitting, which implicitly enhances the generalization capacity of the models and provides the maximum realistic approach for the real-world applications [16].The obtained OR value for all developed models is shown in Table 14.

OR
RMSEofTSdataset RMSEofTRdataset (29) Initially, a score was assigned to each developed model corresponding to their statistical performance measurement parameters, as shown in Table 14.The total score for a particular model was achieved by combining the score obtained in its TR and TS dataset.The overall accuracy of the developed models was determined through the performance strength parameters (using Eq. ( 26)).Based on the S P value obtained for each model, a score was assigned to them as per the ranking analysis method.The final score of the model was accomplished by concatenating its total score value and the score obtained for the S P value.Accordingly, the rank was assigned to the models with respect to their final scoring value.As seen from Table 14, the maximum score was achieved by the GPR F model, therefore, designated as the first rank, followed by GPR K , K-NN K , K-NN F , KRR K , and KRR F models.The ranked models were passed through overfitting analysis.It is clearly observed from the results of the overfitting analysis, shown in Table 14, that the GPR F model exhibits a little amount of overfitting as the OR value is considerably higher than 1.Again, the analysis was performed on a second-ranked model, i.e., the GPR K model.The OR value obtained for the GPR K model is almost close to one, which means that the model is more generalized to the unseen dataset.However, the GPR F model could also be used as a predictive model as the R-value obtained in the TS dataset is 0.83.According to Smith [91], Taskiran [11], Yildirim and Gunaydin [12], Verma [92] and Tenpe and Patel [17], if the R-value of a model is ≥ 0.8, then it is established that the predicted and actual values are good in agreement and strongly correlated.Conclusively, GPR F and GPR K models were finally selected as the best-fitted model for predicting the soaked CBR value of fine-grained plastic soils.

Influence of ML Algorithms and Data Division Approaches on the Model Performance
The overall influence of adopted ML algorithms and data divisional approaches were identified using the final score value (refer to Table 14) obtained for each model.Figures 17  and 18 present the overall influence of ML algorithms and data division approaches, respectively, on the predictive ability of the developed models.It is perceived from Fig. 17 that the GPR algorithm proves its proficiency more than K-NN and KRR algorithms as the combined final score value is significantly higher than those.The combined final score value obtained for the K-Fold approach is slightly more than the FCM approach.However, the K-Fold approach has been found substantially useful in removing the overfitting of the models.

Validation of Literature Study Models
For this purpose, only those models were selected which were having input parameters similar to the present study's geotechnical parameters.Kin [10], Taskiran [11], Yildirim and Gunaydin [12] and Bardhan, Gokceoglu [19] models, belonging to various countries (as shown in Table 15), were attempted to validate through the present study datasets.Datasets from the present study were selected as per their minimum and maximum input and output parameter values.Table 16 depicts the comparative performance of literature study models on the present study datasets.It is observed from Table 16 that the R 2 value obtained for all the models is negative, which means that the selected models don't follow the specific trend of the dataset, leading to a worse fit than the horizontal line.The Bardhan, Gokceoglu [19] model can predict more than 65% observations within ± 20% variations, which is relatively higher than Kin [10], Taskiran [11] and Yildirim and Gunaydin [12] models.This is because the datasets used for the model development and validation are geological and almost nearby, i.e., India.Nagaraj and Suresh [93] also state that soils are likely to be quite variable depending on their geological locations.

Discussion of Results
In the above subsections, the prediction of the soaked CBR of fine-grained plastic soils was assessed through some machine learning algorithms.The collected dataset is comprised of various groups of soils such as ML, CL-ML, CL, MH and CH.Pearson's correlation analysis reveals that S, FC, PL, PI, MDD and OMC were the substantial input parameters in predicting the soaked CBR value of fine-grained plastic soils.Using these input parameters, six models for GPR, K-NN and KRR algorithms were developed through K-Fold and FCM data divisional approaches.The quantitative performance of these models was identified using twelve statistical measurement parameters.Based on the ranking analysis and overfitting ratio, GPR K and GPR F models were considered the best-fitted model, followed by K-NN K , K-NN F , KRR K , and KRR F models.The OR value obtained for GPR K , K-NN K and KRR K models was less than the OR value obtained for GPR F , K-NN F and KRR F models (as shown in Table 14).This means that the RMSE value in the TR and TS datasets of GPR K , K-NN K and KRR K models are almost close to each other.Therefore, the K-Fold approach was considered to be more significant in removing the overfitting of the models.Furthermore, the comparative results of ML algorithms exhibit the maximum proficiency of the GPR algorithm, followed by K-NN and KRR algorithms.The validation results of the literature models on the present study dataset ascertain that the prediction ability of any model in the field of geotechnical/highway engineering is significantly influenced by the geological location of the dataset.Indian 610

Conclusions
The current study offers novel applications of KRR, K-NN and GPR algorithms in predicting the soaked CBR value of fine-grained plastic soils.The analysis was performed on in situ soil samples collected from an ongoing NHAI project work site.Large datasets of 1011 soil samples were collected and laboratory tests were performed as per BIS specifications.The prepared datasets were divided into TR and TS sets using the K-Fold and FCM data divisional approaches.The competency of models was compared through various statistical performance measurement indices, accuracy analysis and REC curve analysis.Apart from them, the final selection of the best-fitted model was accomplished through ranking analysis and overfitting analysis.The obtained results of statistical performance indices reveal that the developed models can explain a maximum variability of 88.7% and a minimum of 64.7% in the CBR value of the TR dataset through S, FC, PL, PI, MDD and OMC as input parameters.The performance of models in the TS dataset was considerably less than that obtained in the TR dataset.Based on the overall analysis of the results, the final selected best-fitted models for predicting the soaked CBR value of fine-grained plastic soils were GPR F and GPR K models.The proposed GPR F and GPR K models can predict 99% and 98% observations, respectively, within ± 20% variations.Additionally, the comparative analysis of the adopted algorithm demonstrates that the GPR algorithm gives the highest proficiency compared to the K-NN and KRR algorithms.Results also indicate that K-Fold and FCM approaches are suitable methods for data division.Eventually, the validation results establish that the predictive ability of any model is substantially influenced by the geological location of the soils/materials used for model development.

Fig. 5
Fig. 5 Data splitting process in K-fold cross-validation approach

Fig. 6
Fig. 6 Silhouette value obtained for each of the clusters 3.4.1)and labeled as Train 1 and Test 1 , respectively.Similarly, C 2 TR and TS dataset was designated as Train 2 and Test 2 , respectively.The final TR dataset was obtained by concatenating the Train 1 and Train 2 datasets.Similarly, the TS dataset was achieved by concatenating the Test 1 and Test 2 datasets.The descriptive statistics value of the final TR and TS dataset is given in Tables

Fig. 7 Fig. 8
Fig. 7 Scatter plot for the TR dataset of KRR, K-NN and GPR models at K-Fold and FCM approaches

Fig. 10
Fig.10 Error frequency plot for the TR and TS dataset of KRR, K-NN and GPR models at K-Fold approach

Fig. 11
Fig. 11 Error distribution plot for the TR and TS dataset of KRR, K-NN and GPR models at FCM approach

Fig. 13
Fig. 13 REC curve for the TR dataset of KRR, K-NN and GPR models at K-Fold approach

Fig. 17 Fig. 18
Fig. 17 Influence of ML algorithms on the model performance

Table 1
Brief overview of the literature study attempted to predict the CBR value of various soil types

Table 2
Combination of KRR algorithm hyperparameters

Table 3
Combination of K-NN algorithm hyperparameters

Table 4
Combination of GPR algorithm hyperparameters

Table 5
Descriptive statistic details for all fine-grained soil parameters

Table 6
Descriptive statistics value for the selected input and output parameters of the TR dataset obtained through the K-fold approach

Table 8
Descriptive statistics value for the selected input and output parameters of the TR dataset obtained through the FCM approach

Table 10
This means that KRR K and KRR F models can explain 64.7%, and 68.1% variability in the soaked CBR value of fine-grained plastic soils.Similarly, the variability explained by K-NN K and K-NN F models are 71.5% and 72.7%, respectively, and by GPR K and GPR F models are 74.6% and 88.7%, respectively.The MAE values obtained for KRR K and KRR F models are 0.515 and 0.512, for K-NN K and K-NN F models are 0.432 and 0.437, and for GPR K and GPR F models are 0.438 and 0.308, respectively.Results of the a20-index demonstrate that all the models can predict almost 98% of observations within ± 20% variations.

Table 12
Accuracy of all the developed models for the TR dataset

Table 13
Accuracy of all the developed models for the TS dataset

Table 14
Rank analysis and overfitting ratio for all the developed models

Table 15
Literature models validation on the present study dataset

Table 16
Comparative performance of literature models on the present study datasets