Modelling Soil Compaction Parameters Using an Enhanced Hybrid Intelligence Paradigm of ANFIS and Improved Grey Wolf Optimiser

: The criteria for measuring soil compaction parameters, such as optimum moisture content and maximum dry density, play an important role in construction projects. On construction sites, base/sub-base soils are compacted at the optimal moisture content to achieve the desirable level of compaction, generally between 95% and 98% of the maximum dry density. The present technique of determining compaction parameters in the laboratory is a time-consuming task. This study proposes an improved hybrid intelligence paradigm as an alternative tool to the laboratory method for estimating the optimum moisture content and maximum dry density of soils. For this purpose, an advanced version of the grey wolf optimiser (GWO) called improved GWO (IGWO) was integrated with an adaptive neuro-fuzzy inference system (ANFIS), which resulted in a high-performance hybrid model named ANFIS-IGWO. Overall, the results indicate that the proposed ANFIS-IGWO model achieved the most precise prediction of the optimum moisture content (degree of correlation = 0.9203 and root mean square error = 0.0635) and maximum dry density (degree of correlation = 0.9050 and root mean square error = 0.0709) of soils. The outcomes of the suggested model are noticeably superior to those attained by other hybrid ANFIS models, which are built with standard GWO, Moth-ﬂame optimisation, slime mould algorithm, and marine predators algorithm. The results indicate that geotechnical engineers can beneﬁt from the newly developed ANFIS-IGWO model during the design stage of civil engineering projects. The developed MATLAB models are also included for determining soil compaction parameters.


Introduction
In the parlance of geotechnical engineering, soil compaction is a method of compressing soil particles by reducing air voids while maintaining steady water content [1,2]. Compaction can be used to enhance the mechanical qualities of soils in a number of different ways. Proctor [3] recommended compacting soil with different water contents at an appropriate compaction energy. As a result, the compaction curve can be used to determine the optimum moisture content (OMC) and maximum dry density (MDD) of soils. To sustain the long-term performance of various engineering structures, such as embankments of railways, highways, and airport runways, these two compaction parameters are commonly deals with knowledge that can be obtained and comprehended, neural networks deal with knowledge that can be obtained only via optimised learning [28,30]. However, like many other MLTs, ANFIS has some limitations, such as overfitting issues. Additionally, because it is hard to define the exact global optimum, it may produce undesirable outcomes during the validation phase [14,31].
To solve these issues, researchers have employed a number of meta-heuristic algorithms (MHAs), such as GA, PSO, GWO, etc. [28,29,32], and a number of hybrid models of traditional MLTs and MHAs were built for the estimation of desired output(s). It is important to note that construction of an effective ANFIS model requires optimum selection of its consequent and antecedent (C&A) and fuzzy inference system (FIS) parameters. These two parameters significantly affect how the learning phase turns out, which in turn affects how well a hybrid ANFIS model can predict the desired variables. Due to the robust global search capabilities of MHAs, the C&A parameters of ANFIS are iteratively adjusted, resulting in improved performance. Over the past decade, several hybrid ANFIS models have seen widespread use in addressing a wide range of engineering problems, including compressive strength estimation [28,33], flood assessment [34], prophecy of groundwater level [35], and so on.
Nevertheless, a detailed review of the literature reveals that no previous study has employed hybrid ANFIS models constructed with a specific group of MHAs to predict soil compaction parameters. On the other hand, it is important to highlight that no algorithm provides perfect solutions for all optimisation problems due to improper exploration and exploitation (E&E) processes [31,36]. Therefore, implementing a standard version of MHA in hybrid modelling does not ensure optimum hybrid model generation. It may also be noted that researchers reported modified versions of MHAs and demonstrated that the performance of standard MHA could be improved by implementing different strategies [37][38][39]. Considering these points as a reference, and to fill the gap in the literature, an enhanced hybrid technique of ANFIS and an improved grey wolf optimiser (IGWO), i.e., ANFIS-IGWO, has been constructed and presented in this study for the estimation of OMC and MDD of soils. The performance of the ANFIS-IGWO model was compared to that of three hybrid ANFIS models built using moth-flame optimisation (MFO), slime mould algorithm (SMA), and marine predator algorithm (MPA). The performance of the ANFIS-IGWO model was also compared with the standard hybrid model of ANFIS and GWO, i.e., ANGIS-GWO. Thus, as a part of ongoing research and to extend the work of Bardhan and Asteris [14], a suitable database of various soils was compiled from the studies of Günaydın [15], Wang and Yin [1], and Bardhan and Asteris [14] and a modified database was prepared. Specifically, a total of 251 datasets from 15 different soils were acquired and utilised in the current study for the estimation of the OMC and MDD of soils.
The remainder of this work is organised as follows. The significance of the present study is presented in Section 2. Section 3 details the methodological development of ANFIS-based hybrid models. Section 4 discusses data collection, descriptive details, and the computer modelling procedure. Section 5 provides and discusses the realisations of the developed models, followed by Section 6 with the limitations and future scope of the study. At the end, summary and conclusions are presented.

Research Significance
In the last two decades, a multitude of modern computational methods, techniques, and algorithms have been proposed and published with the aim of predicting the response of complex phenomena whose strongly non-linear nature and behaviour make impossible the widely accepted use of deterministic techniques [25,26]. In these methods, artificial intelligence, machine learning, and MHAs have a dominant position. In fact, despite the fact that these techniques started with the first applications in medicine [40], they were particularly applied in the fields of sciences [30,[41][42][43] and engineering [44][45][46][47][48]. The use of contemporary intelligence techniques in geotechnical and geological engineering domains, such as landslide susceptibility mapping [49], reliability analysis [50], and es-timation of various geotechnical parameters [24,27] can also be found in the literature. However, the existing literature in the geotechnical engineering area does not demonstrate sufficient implementation of enhanced/improved versions of MHAs in estimating various geotechnical parameters. Taking the above discussion as a reference, this study proposes a high-performance intelligence paradigm built using an upgraded version of MHA for estimating the OMC and MDD of soils.

Methodology
This section presented the theoretical details of GWO and IGWO, followed by a short discussion on MFO, SMA, and MPA. Subsequently, the methodological development of hybrid ANFIS models is presented and discussed. However, before presenting the above details, the working principles of the ANFIS are briefly presented.

Adaptive Neuro-Fuzzy Inference System
ANFIS, proposed by Jang [51], is an ANN-FIS integration, which was intended to eliminate the drawbacks of the individual ANN and FIS approaches. ANFIS is grounded in fuzzy logic and rules produced in the particular training procedure of the model. These inference systems contain five layers (see Figure 1). The nodes of layer 0 are the inputs, while the nodes of layer 5 represent the output in the connection-based structure. The fixed adaptable nodes of the hidden layers stand for the membership functions (MFs).

Grey Wolf Optimiser
GWO [52] is comes under the category of evolutionary algorithm developed for timisation based on the imitation of the grey wolves' social behaviour. Specifically, algorithm mimics the process that grey wolves utilise to capture their prey, along with structure of their leadership. For the recreation of the hierarchical structure in GWO, g wolves of four different types are assumed for every wolf pack. The leader and the m significant wolf in the pack are called , and wolves. wolves with minimum sponsibility are placed at the bottom of the food. In GWO, the entire hunting process be classified as searching, encircling, hunting, and attacking. The mathematical express for encircling prey is given by: For a summarised description of the ANFIS approach, let x 1 and x 2 be the inputs.
Additionally, let f be the output. The ANFIS represents the relationship between the inputs and output by fuzzy if-then rules. The Takagi-Sugeno fuzzy rules in the model are shown as: where A 1 , A 2 , B 1 , and B 2 are linguistic symbols, while p 1 , q 1 , r 1 , p 2 , q 2 , and r 2 are the consequent variables. The layers include: Layer 1: Fuzzification layer-it is assumed that node i has an adaptive function as: where O 1,i is the output of node i, while µ A i denotes the MF. Layer 2: Ruler layer-Node i within this layer is assumed to be fixed (II). In addition, the node output is generated by incoming signals, such as where O 2,i denotes the output of the second layer, and w i represents the firing strength of rule i. Layer 3: Normalisation layer-Node i undergoes normalisation in the third layer (firing strengths). The ratio of the firing force of rule i to the total firing force can be obtained as: i is the output of the third layer, and w i stands for the normalised firing strength. Layer 4: Defuzzification layer-In this layer, each of the nodes is adaptive and has a function representing the contribution of rule i to the total output. Layer 5: Output layer-Eventually, this layer yields the final output.
ANFIS can be equipped for MF parameter identification through hybrid learning approaches. Such approaches define the parameters of the defuzzification layer using the forward least squares technique. Errors undergo backpropagation to modify a i , b i, and c i as the premise parameters through gradient descent. For more details, the works of Paryani et al. [49], Piro et al. [30], Golafshani et al. [28], can be referred to.

Grey Wolf Optimiser
GWO [52] is comes under the category of evolutionary algorithm developed for optimisation based on the imitation of the grey wolves' social behaviour. Specifically, this algorithm mimics the process that grey wolves utilise to capture their prey, along with the structure of their leadership. For the recreation of the hierarchical structure in GWO, grey wolves of four different types are assumed for every wolf pack. The leader and the most significant wolf in the pack are called α, β and δ wolves. ω wolves with minimum responsibility are placed at the bottom of the food. In GWO, the entire hunting process can be classified as searching, encircling, hunting, and attacking. The mathematical expression for encircling prey is given by: where X and X p are the position vectors of the grey wolf and the prey, respectively; t and t + 1 represent current and subsequent epochs, respectively. A and C are two vectors given by: where r 1 and r 2 are the two random vectors that are uniformly distributed [0 1], and the components of a are linearly decreased from 2 to 0. When |A| > 1, the exploration of prey location is possible by diverting the search agents. Conversely, with |A| < 1, convergence of search agents can be used to achieve exploitation. The hunting process in GWO can be mathematically modelled as follows: In GWO, E&E is handled using parameters a and C, in which the parameter a is decreased from 2 to 0. Additionally, it is seen that the final position would be in a random place within a circle, which is defined by the positions of α, β, and δ in the search space. More mathematical details can be found in the original work of Mirjalili et al. [52].

Improved Grey Wolf Optimiser
In GWO, α, β, and δ guide ω wolves toward regions of the search space where the optimal solution is likely to be located. This approach may result in entanglement in a locally optimal solution. Another drawback is the decline in population diversity, which causes GWO to approach the local optimum. Nadimi-Shahraki et al. [53] proposed IGWO to address these problems. According to the study of Nadimi-Shahraki et al. [53], the enhancements involve a new search strategy involving a step of selecting and upgrading. Therefore, IGWO consists of three phases, as discussed below.
Initialising phase: During the initialisation phase, N wolves are randomly distributed in [l j , u j ], as: The position of the i-th wolf in the t-th iteration, represented by X i (t) = {x i1 , x i2 , . . . , x iD }, where D is the dimension number. The fitness value of X i (t) is calculated using f (X i (t)).
Movement phase: The IGWO, proposed by Nadimi-Shahraki et al. [53], includes a different mobility tactic known as the dimension learning-based hunting (DLH) method, in which each wolf is learned by its neighbours to be a different contender for the new position, X i (t).
Selecting and updating phase: During this stage, the best candidate is first chosen by contrasting the fitness ratings between two candidates X i-GWO (t + 1) and X i-DLH (t + 1), given by: Then, to update the position of X i (t + 1), if the fitness of the selected candidate is less than X i (t), X i (t) is updated by the selected candidate. Otherwise, X i (t) remains unchanged. After this procedure, the iteration count is increased by 1, and the search operation is repeated until the predetermined number of epochs has been reached.

Brief Overview of MFO, SMA, and MPA
The other employed MHAs, viz., MFO, SMA, and MPA, are briefly discussed in this sub-section. All of these MHAs are swarm-based and they have been widely used in different engineering disciplines [54][55][56][57].
MFO, proposed by Mirjalili [58], is an innovative MHA that draws inspiration from the intriguing behaviour of moths attracted to flames. MFO incorporates the unique phenomenon of moths spiralling around a flame into its search strategy. This behaviour, while seemingly irrational and perilous for the moths, serves as a metaphor for E&E in optimisation problems. MFO leverages a chaotic search mechanism that emulates the unpredictable flight patterns of moths around a flame. This mechanism enables the MFO to efficiently explore diverse solution spaces, avoiding stagnation in the local optima. By introducing chaos, the MFO promotes global exploration while maintaining its ability to exploit promising regions of the search space. The core idea behind MFO is to strike a balance between E&E, mimicking the trade-off faced by moths as they navigate the dangerous allure of flames. By dynamically adjusting the balance between E&E strategies, the MFO adaptively evolves its search behaviour, allowing it to effectively handle complex optimisation problems with varying landscapes. The effectiveness of MFO has been demonstrated across a wide range of applications, including engineering design, data mining, and image processing. Its ability to handle both continuous and discrete optimisation problems make it a versatile tool in the field of MHAs.
SMA [59] simulates the nutritive phase of a slime mould as a unique approach that is grounded in nature (a single-celled eukaryote). The foraging behaviour of slime moulds is simulated by this programme. By smelling potential food sources, slime moulds locate them, wrap them, and then digest them by secreting enzymes. In SMA, the phase of iterations to produce the highest smell concentration is the theoretical description of how to approach the optimal solution. The slime mould's flexible weight ensures rapid convergence and prevents it from becoming stranded in regional extremes. This approach enables the slime mould to advance along any viable path in the direction of the ideal outcome, which mimics the slime mould's eating-related architecture. The next stage is wrapping the meal using contractions of the intravenous framework inside the upper and lower limitations. The vein with the maximum contraction of food generates more bio-oscillator waves, which cause the cytoplasm to flow more quickly through the vein, increasing its thickness. The search patterns in SMA are altered in response to the opposing signals from veins regarding the concentration of food.
MPA [60], a MHA inspired by the natural principles governing optimal foraging strategies and encounter rates between predator and prey in marine ecosystems. Marine predators adopt a Lévy strategy when navigating environments with scarce prey, while employing Brownian movement in areas abundant with prey. Throughout their lifetimes, these predators exhibit a consistent balance of Lévy and Brownian movement as they traverse diverse habitats. Environmental factors, such as eddy formation, influence their behaviour, prompting adaptive changes to explore regions with varying prey distributions. Leveraging their remarkable memory capabilities, they capitalise on the recollection of successful foraging locations and associations with other individuals. MPA harnesses these concepts to guide its search process, mimicking the adaptive foraging behaviour of marine predators. By integrating these nature-inspired mechanisms, MPA demonstrates a powerful optimisation approach capable of addressing complex problems in diverse domains.
Note that the detailed working principles of these OAs are not presented in this study because they are well established, and the original studies of MFO [58], SMA [59], and MPA [60] can be referred to for more details.

Hybrid Modelling of ANFIS and MHAs
In this work, the C&A parameters of ANFIS were optimised using MHAs. It is important to note that proper setting of the FIS and C&A parameters is necessary for creating an optimum ANFIS model because learning parameters have a significant impact on the model's performance. Notably, the selection of all of the ANFIS hyperparameters at once is a challenging operation because they must be searched in continuous domains, leading to an infinite number of parameters sets. As a result, it is possible to define the problem of ANFIS parameter tweaking as an optimisation problem. Thus, the values of the FIS and C&A parameters were optimised using IGWO, GWO, MFO, SMA, and MPA, and five hybrid ANFIS models, i.e., ANFIS-IGWO, ANFIS-GWO, ANFIS-MFO, ANFIS-SMA, and ANFIS-MPA, were created. A flow chart of the construction procedure of the hybrid ANFIS models is presented in Figure 2.

Data Description and Modelling
A broad variety of experimental results of soil compaction parameters were acquired from the studies of Günaydın [15], Wang and Yin [1], and Bardhan and Asteris [14]. Specifically, a total of 372 results were obtained, the details of which are presented in Table 1. The work of Günaydın [15] consists of 126 compaction results of nine distinct soil types (CH, CI, CL, GC, GM, MH, MI, ML, and SC) with six influencing parameters, viz., fines content (F), sand content (S), gravel content (G), specific gravity, liquid limit (LL), and plastic limit (PL). Wang and Yin [1] gathered a total of 226 records from the literature. The database includes G, S, F, LL, PL, and compaction energy of various soil types, such as CL, CL-ML, CH, MH, ML, SC, SP-SC, SW-SC, SM, GC, GP-GC, GW-GC, and GM. Bardhan and Asteris [14] presented 20 experimental records of soil compaction parameters, including four distinct soil types (CH, CI, CL, and SC) and six influencing parameters, identical to Günaydn [15]. According to the study of Wang and Yin [1], the majority of 226 soil compaction experiments were conducted using either the conventional Proctor or the reduced compaction energy. Additionally, thirty modified Proctor compaction tests were incorporated into the database. In this study, the database presented by Wang and Yin [1] has been revised, and a total of 105 records were chosen. Additionally, all datasets of Günaydın [15] and Bardhan

Data Description and Modelling
A broad variety of experimental results of soil compaction parameters were acquired from the studies of Günaydın [15], Wang and Yin [1], and Bardhan and Asteris [14]. Specifically, a total of 372 results were obtained, the details of which are presented in Table 1. The work of Günaydın [15] consists of 126 compaction results of nine distinct soil types (CH, CI, CL, GC, GM, MH, MI, ML, and SC) with six influencing parameters, viz., fines content (F), sand content (S), gravel content (G), specific gravity, liquid limit (LL), and plastic limit (PL). Wang and Yin [1] gathered a total of 226 records from the literature. The database includes G, S, F, LL, PL, and compaction energy of various soil types, such as CL, CL-ML, CH, MH, ML, SC, SP-SC, SW-SC, SM, GC, GP-GC, GW-GC, and GM. Bardhan and Asteris [14] presented 20 experimental records of soil compaction parameters, including four distinct soil types (CH, CI, CL, and SC) and six influencing parameters, identical to Günaydn [15]. According to the study of Wang and Yin [1], the majority of 226 soil compaction experiments were conducted using either the conventional Proctor or the reduced compaction energy. Additionally, thirty modified Proctor compaction tests were incorporated into the database. In this study, the database presented by Wang and Yin [1] has been revised, and a total of 105 records were chosen. Additionally, all datasets of Günaydın [15] and Bardhan and Asteris [14] were used. Therefore, 126, 105, and 20 experimental records were acquired from the studies of Günaydın [15], Wang and Yin [1], and Bardhan and Asteris [14], respectively. The details of data dimension are also presented in Table 2. Therefore, the final database includes 251 records and five influential parameters viz., F in %, S in %, G in %, LL in %, and PL in %, of 15 different soil types. These five influential parameters were used to estimate the OMC and MDD of soils. Descriptive details of the final dataset are given in Table 2. In addition, the minimum and maximum values of influential (soil-type wise) and compaction parameters are presented in Table 3. Note that the abbreviations of soil types are presented as per the Indian Standard Soil Classification System (ISSCS) and ASTM [61].   Figure 3 shows the comparative histograms for each influential variable. To better illustrate, the correlation matrices between influential variables and compaction parameters are presented in Figure 4. From the information given in Table 2, Figures 3 and 4, it can be seen that the OMC has a negative correlation with the contents of S and G, whereas F, LL, and PL show a positive correlation. In contrast, F, LL, and PL exhibit negative correlations, while S and G contents have positive correlations with MDD. Notably, these figures are particularly useful, as they indicate the range of values of the parameters for which the reliability is limited, and further experimental investigation is required for values of the parameters included in these regions and not with the aim of updating the database in the future.
After finalising the database, it was divided into two subsets: a training (TR) subset that contained 80% of the overall dataset and a testing (TS) subset that contained the remaining 20% of the data. The following steps can be used to describe the computational modelling process for estimating soil compaction parameters: (a) choosing the main dataset; (b) data normalisation; (c) data partitioning and selection of TR and TS subsets; (d) model  After finalising the database, it was divided into two subsets: a training (TR) subset that contained 80% of the overall dataset and a testing (TS) subset that contained the remaining 20% of the data. The following steps can be used to describe the computational modelling process for estimating soil compaction parameters: (a) choosing the main dataset; (b) data normalisation; (c) data partitioning and selection of TR and TS subsets; (d) model construction using a training subset; (e) check model performance; (f) check terminating criteria; (g) model validation if terminating criteria are satisfied; and (h) performance assessment. The steps of computational modelling are illustrated in Figure 5.  After finalising the database, it was divided into two subsets: a training (TR) subset that contained 80% of the overall dataset and a testing (TS) subset that contained the remaining 20% of the data. The following steps can be used to describe the computational modelling process for estimating soil compaction parameters:

Results and Discussion
The outcomes of the hybrid ANFIS models used to estimate soil compaction parameters are described in this section. As previously mentioned, the primary dataset was divided into training (201 samples) and testing (50 samples) subsets before the models were built. Note that all models were constructed and validated using identical training and testing subsets. The output of the developed models was then assessed using a number of indices, namely performance index (PFI), correlation coefficient (R), variance account factor (VAF), Willmott's Index of agreement (WI), mean absolute error (MAE), root mean square error (RMSE), RMSE to observation's standard deviation ratio (RSR), and weighted mean absolute percentage error (WMAPE). Notably, these indices are frequently used to evaluate the generalisability of any prediction model from a variety of perspectives, including correlation accuracy, related error, amount of variation, and so on.
In contrast, the deterministic parameters of MHAs, such as swarm size, maximum iteration number, and upper and lower bounds, play a vital part in hybrid modelling; thus, they were calibrated throughout the optimisation process. The details of the deterministic and hyper-parameters of hybrid ANFIS models in estimating soil compaction parameters are described in the following sub-section, followed by a comparative assessment of the results.

Results and Discussion
The outcomes of the hybrid ANFIS models used to estimate soil compaction parameters are described in this section. As previously mentioned, the primary dataset was divided into training (201 samples) and testing (50 samples) subsets before the models were built. Note that all models were constructed and validated using identical training and testing subsets. The output of the developed models was then assessed using a number of indices, namely performance index (PFI), correlation coefficient (R), variance account factor (VAF), Willmott's Index of agreement (WI), mean absolute error (MAE), root mean square error (RMSE), RMSE to observation's standard deviation ratio (RSR), and weighted mean absolute percentage error (WMAPE). Notably, these indices are frequently used to evaluate the generalisability of any prediction model from a variety of perspectives, including correlation accuracy, related error, amount of variation, and so on.
In contrast, the deterministic parameters of MHAs, such as swarm size, maximum iteration number, and upper and lower bounds, play a vital part in hybrid modelling; thus, they were calibrated throughout the optimisation process. The details of the deterministic and hyper-parameters of hybrid ANFIS models in estimating soil compaction parameters are described in the following sub-section, followed by a comparative assessment of the results.
where y i = actual ith value;ŷ i = estimated ith value; n = is the number of samples; and y mean = mean of the actual value. Note that for a perfect predictive model, the values of the aforementioned indices should be identical to their identical values, the details of which can be obtained from the literature [14,17].

Model Performance
The results of the hybrid ANFIS models that were built to estimate soil OMC and MDD are presented in this sub-section. As stated above, the optimum selection of hyperparameters is a challenging operation, and hence, proper tuning of FIS and C&A parameters was performed during the course of hybrid modelling. The number of FIS parameters (N FIS ) were investigated between 2 and 15. Using Gaussian MF and RMSE as fitness functions, the most appropriate value of N FIS was determined to be 5. Notably, Gaussian and linear MFs were used in the input and output layers, respectively. A total of 60 C&A membership functions of ANFIS were optimised for the nine-dimensional input space. Note that the optimised values of N FIS and C&A were chosen following a trial-and-error approach and according to the performance during the testing phase. The convergence behaviour and computational time of the developed hybrid ANFIS models are presented in Figures 6 and 7, respectively. It should be noted that the computational time of the developed ANFIS-IGWO model was found to be longer due to the use of an upgraded version of GWO that required changed mathematical calculations to handle E&E operations. Moreover, it is seen that all the developed hybrid models converge within 500 epochs; hence, they are considered to be sufficient as the maximum iteration count.
The performance of the developed ANFIS models is presented in Tables 4 and 5, respectively, for the OMC and MDD estimations. The abilities of the constructed models for training, testing, and total outputs are shown here. It should be underlined that the training subset performance was used to define the goodness of fit of the developed models, while the testing dataset was used to evaluate their generalisation potential. According to Table 4, it is seen that the developed ANFIS-MPA achieved the highest R and lowest RMSE values of 0.9335 and 0.0590, respectively, during the training phase of OMC prediction. However, during the testing phase, the constructed ANFIS-IGWO achieved the most precise precision, with R = 0.8645 and RMSE = 0.0754. According to the overall results of the OMC estimation, the ANFIS-IGWO was determined to be the best-fitted model with R = 0. On the contrary, the results of Table 5  and linear MFs were used in the input and output layers, respectively. A total of 60 C&A membership functions of ANFIS were optimised for the nine-dimensional input space. Note that the optimised values of NFIS and C&A were chosen following a trial-and-error approach and according to the performance during the testing phase. The convergence behaviour and computational time of the developed hybrid ANFIS models are presented in Figures 6 and 7, respectively. It should be noted that the computational time of the developed ANFIS-IGWO model was found to be longer due to the use of an upgraded version of GWO that required changed mathematical calculations to handle E&E operations. Moreover, it is seen that all the developed hybrid models converge within 500 epochs; hence, they are considered to be sufficient as the maximum iteration count.  The performance of the developed ANFIS models is presented in Tables 4 and 5, respectively, for the OMC and MDD estimations. The abilities of the constructed models for training, testing, and total outputs are shown here. It should be underlined that the training subset performance was used to define the goodness of fit of the developed models,  The performance of the developed ANFIS models is presented in Tables 4 and 5, respectively, for the OMC and MDD estimations. The abilities of the constructed models for training, testing, and total outputs are shown here. It should be underlined that the training subset performance was used to define the goodness of fit of the developed models, while the testing dataset was used to evaluate their generalisation potential. According to   To better demonstrate the performance of the developed ANFIS models, scatterplots are presented in Figures 8 and 9 for the OMC and MDD estimations, respectively. Herein, the illustrations of actual and estimated values for the best three models (based on RMSE value) are shown. The amount of variance in these diagrams can be visualised by viewing red-coloured dotted lines put at 10% levels. The performance of the generated hybrid models is compared in the following sub-section, and a comparative assessment is presented using a variety of graphical illustrations.
Mathematics 2023, 11, x FOR PEER REVIEW 15 of 23 models is compared in the following sub-section, and a comparative assessment is presented using a variety of graphical illustrations.  models is compared in the following sub-section, and a comparative assessment is presented using a variety of graphical illustrations.

Discussion of Results
It is critical to highlight that a data-driven model is incomplete without a visual representation of the results. Visualisations enable the identification of the degree of accuracy and associated errors in a model that is easier to comprehend. Therefore, the results are displayed in the form of an accuracy matrix, Taylor diagrams, and radar plots. Notably, the Taylor diagrams are shown for the testing dataset only, since the performance of a data-driven model during the testing dataset should be accepted with more certainty. For thoroughly assessing a model's overall correctness, these diagrams are quite beneficial. An accuracy matrix is a heat map matrix that is used to measure the level of accuracy that a model achieves in terms of certain performance criteria. This matrix makes it simple to evaluate a model's correctness without having to look up each index's value. As previously stated, a number of indices must be created to assess a model's accuracy from multiple angles; however, interpreting results by looking at each index's values takes time and necessitates in-depth observation. Figure 10 shows the accuracy matrix for the models created for the OMC and MDD predictions. The accuracy matrix demonstrates that the constructed ANFIS-IGWO attained higher predictive precision against each index during the testing phase.
On the other hand, the Taylor diagram [62] is used to provide a quick assessment of a model's accuracy in terms of the coefficient of correlation, ratio of standard deviations, and RMSE index. Generally, a point inside a Taylor diagram indicates a model. The position of the point should line up with the reference point for an ideal model. The Taylor diagrams for the hybrid ANFIS models developed for OMC and MDD predictions are shown in Figure 11. In addition to the accuracy matrix and Taylor diagrams, radar plots representing the R value are also presented in Figure 12 for the training, testing, and total cases of OMC and MDD estimations. A ridgeline chart and distribution with Kernel smooth of error between actual and estimated values are presented in Figure 13. From these diagrams, the predictive capability of hybrid ANFIS models can be assessed from different perspectives.
However, according to the aforementioned results, the ANFIS-IGWO model was found to be the best-obtained model in both instances of prediction. As indicated previously, eight indices were used to evaluate the performance of the developed ANFIS models. Based on the overall results against OMC prediction, the developed ANFIS-IGWO achieved the highest level of accuracy, with R = 0.9203 and RMSE = 0.0635, whereas against MDD prediction, R = 0.9050 and RMSE = 0.0709 achieved the highest level of accuracy. Therefore, the suggested ANFIS-IGWO model can be used to approximate the OMC and MDD ranges for various soil types. This will aid engineers and practitioners in reducing the operational time required for laboratory compaction experiments. The developed MATLAB models, as well as the employed dataset, are included as Supplementary Materials for future use. The details of MATLAB implementations of the developed models are also presented in Appendix A. For better demonstration, the steps of OMC and MDD estimations using basic soil parameters and the developed MATLAB models are illustrated in Figure 14.
evaluate a model's correctness without having to look up each index's value. As previously stated, a number of indices must be created to assess a model's accuracy from multiple angles; however, interpreting results by looking at each index's values takes time and necessitates in-depth observation. Figure 10 shows the accuracy matrix for the models created for the OMC and MDD predictions. The accuracy matrix demonstrates that the constructed ANFIS-IGWO attained higher predictive precision against each index during the testing phase.   diagrams for the hybrid ANFIS models developed for OMC and MDD predictions are shown in Figure 11. In addition to the accuracy matrix and Taylor diagrams, radar plots representing the R value are also presented in Figure 12 for the training, testing, and total cases of OMC and MDD estimations. A ridgeline chart and distribution with Kernel smooth of error between actual and estimated values are presented in Figure 13. From these diagrams, the predictive capability of hybrid ANFIS models can be assessed from different perspectives.   However, according to the aforementioned results, the ANFIS-IGWO model was found to be the best-obtained model in both instances of prediction. As indicated previously, eight indices were used to evaluate the performance of the developed ANFIS models. Based on the overall results against OMC prediction, the developed ANFIS-IGWO achieved the highest level of accuracy, with R = 0.9203 and RMSE = 0.0635, whereas against MDD prediction, R = 0.9050 and RMSE = 0.0709 achieved the highest level of accuracy. Therefore, the suggested ANFIS-IGWO model can be used to approximate the OMC and MDD ranges for various soil types. This will aid engineers and practitioners in reducing the operational time required for laboratory compaction experiments. The developed MATLAB models, as well as the employed dataset, are included as supplementary materials for future use. The details of MATLAB implementations of the developed models are also presented in Appendix A. For better demonstration, the steps of OMC and MDD estimations using basic soil parameters and the developed MATLAB models are illustrated in Figure 14.

Limitations and Future Research
In this section, the limitations of the proposed ANFIS-IGWO model are presented, as well as the main points that need to be further investigated. Regarding the estimation of the OMC and MDD of soils, it is worth noting that despite the excellent prediction it achieves, and which are presented in the previous section, it only applies to values of the input parameters between the minimum and the maximum values (refer to Table 2) that have the corresponding parameters of the database used to train and develop the model. Additionally, the reliability of the proposed model is limited for parameter value ranges, where the number of experimental data and soil types are not sufficient. That is, in cases in which the number of data is very small, and thus unable to satisfactorily describe the soil compaction. For example, based on the histograms of Figure 3, it was found that there is not enough data in some regions. Such value ranges, where there is a shortage of each of the input parameters, should be studied experimentally in the near future and updated the database with the aim of future development of more efficient forecasting soft computing models.

Summary and Conclusions
Soil compaction parameters play a vital role in construction projects. They are crucial for comparing the level of compaction achieved in the field. However, the traditional laboratory method for determining OMC and MDD is time consuming. Therefore, the main objective of this study is to sidestep the need for multiple laboratory tests by leveraging the predictive capabilities of high-performance hybrid intelligence paradigms. Taking these points into consideration, the current study proposes a high-performance hybrid model to sidestep the operation of typical laboratory testing of soil compaction parameters. To achieve this goal, an ANFIS-IGWO model was constructed, and the performance of this model was compared with four hybrid ANFIS models, namely ANFIS-GWO, AN-FIS-MFO, ANFIS-SMA, and ANFIS-MPA. The experimental results clearly demonstrate that the proposed ANFIS-IGWO model effectively predicts soil compaction parameters. With an accuracy range of 92.5% to 94%, according to the RMSE index, the developed ANFIS-IGWO model exhibits superior generalisation abilities for estimating soil compaction parameters.

Limitations and Future Research
In this section, the limitations of the proposed ANFIS-IGWO model are presented, as well as the main points that need to be further investigated. Regarding the estimation of the OMC and MDD of soils, it is worth noting that despite the excellent prediction it achieves, and which are presented in the previous section, it only applies to values of the input parameters between the minimum and the maximum values (refer to Table 2) that have the corresponding parameters of the database used to train and develop the model. Additionally, the reliability of the proposed model is limited for parameter value ranges, where the number of experimental data and soil types are not sufficient. That is, in cases in which the number of data is very small, and thus unable to satisfactorily describe the soil compaction. For example, based on the histograms of Figure 3, it was found that there is not enough data in some regions. Such value ranges, where there is a shortage of each of the input parameters, should be studied experimentally in the near future and updated the database with the aim of future development of more efficient forecasting soft computing models.

Summary and Conclusions
Soil compaction parameters play a vital role in construction projects. They are crucial for comparing the level of compaction achieved in the field. However, the traditional laboratory method for determining OMC and MDD is time consuming. Therefore, the main objective of this study is to sidestep the need for multiple laboratory tests by leveraging the predictive capabilities of high-performance hybrid intelligence paradigms. Taking these points into consideration, the current study proposes a high-performance hybrid model to sidestep the operation of typical laboratory testing of soil compaction parameters. To achieve this goal, an ANFIS-IGWO model was constructed, and the performance of this model was compared with four hybrid ANFIS models, namely ANFIS-GWO, ANFIS-MFO, ANFIS-SMA, and ANFIS-MPA. The experimental results clearly demonstrate that the proposed ANFIS-IGWO model effectively predicts soil compaction parameters. With an accuracy range of 92.5% to 94%, according to the RMSE index, the developed ANFIS-IGWO model exhibits superior generalisation abilities for estimating soil compaction parameters.
According to the overall outcomes, the proposed ANFIS-IGWO model offers a significant advantage by transforming the C&A parameters of the model into the coordinates of individual wolves within the community. Each wolf's position represents a result of the ANFIS-IGWO model, with N FIS = 5 and 500 epochs utilised. However, a drawback of the ANFIS-IGWO model is its high computational cost due to the implementation of a modified approach. Multiple runs were performed to identify the most suitable search space for accurate output estimation, further increasing the time required. Other limitations encompass no external validation performed and the exclusion of factors such as compaction energy and the parental significance of soils during modelling. Therefore, additional research is needed to expand the application of the suggested ANFIS-IGWO model in estimating soil compaction parameters.

Conflicts of Interest:
The authors declare no conflict of interest.