Abstract
In this paper, the deep eutectic solvent-functionalized carbon nanotube was used for arsenic removal from water solution. The adsorbent used was characterized using Raman spectroscopy, Fourier transform infrared (FTIR) and zeta potential. The effect of the parameters (adsorbent dosage, pH, initial concentration and contact time) was studied to find the optimum conditions for maximum adsorption capacity of the functionalized carbon nanotube. The pseudo-second-order, the pseudo-first-order and intraparticle diffusion kinetic models were applied to identify the adsorption rate and mechanism, the pseudo-second-order model best described the adsorption kinetics of the system. The non-linear autoregressive network with exogenous inputs (NARX) neural network strategy was used for the modelling and predicting of the adsorption capacity of functionalized carbon nanotube. Different indicators were used to determine the efficiency and accuracy of the NARX neural network model which were mean square error (MSE), root mean square error (RMSE), relative root mean square error (RRMSE) and mean absolute percentage error (MAPE). The sensitivity study of the used parameters in the experimental work was completed. Comparison of the NARX model results with the experimental data confirmed that the NARX model was able to predict the arsenic removal from water.
LIST OF ABBREVIATIONS
- ANN
Artificial neural network
- CNT
Carbon nanotube
- DAC
N,N-diethylethanolammonium chloride
- DES
Deep eutectic solvent
- FTIR
Fourier transform infrared
- Gly
Glycerol
- HBA
Hydrogen bond acceptor
- HBD
Hydrogen bond donor
- MAPE
Mean absolute percentage error
- MSE
Mean square error
- MWCNT
Multi-wall carbon nanotube
- NARX
Non-linear autoregressive network with exogenous inputs
- P-CNT
Pristine carbon nanotube
- Qc
Adsorption capacity
- RE
Relative error
- RMSE
Root mean square error
- RRMSE
Relative root mean square error
INTRODUCTION
The presence of heavy metal ions in water is considered as a major problem due to their non-biodegradability, toxicity and human health complications. Arsenic is one of the heavy metals known as a carcinogenic material to humans. It can be found in polluted groundwater as a result of industrial waste discharge, rock weathering, pesticides and arsenical herbicides use for agricultural purposes (Xia et al. 2014). The exposure and consumption of arsenic polluted drinking water causes numerous health problems in several countries such as Bangladesh (Wasserman et al. 2004), Bengal (Mazumder et al. 1997) and China (Wang et al. 2007). Due to the serious health problems related to arsenic in drinking water, the World Health Organization determined the maximum allowable arsenic in drinking water as 10 μg/L (Smith et al. 2002). Several hundreds of mg/L of arsenic concentration are available in groundwater. Consequently, millions of people are affected due to the high contamination of drinking water with arsenic (Argos et al. 2010). Many chemical, biological and physical processes have been utilized for heavy metal treatment such as precipitation, ion-exchange, reverse osmosis, biosorption, filtration and adsorption (Ding et al. 2014; Ahmadi et al. 2015; Cho et al. 2015; Luo et al. 2015). Adsorption is the most suitable technique due to its cost effectiveness, feasible operation and high removal efficiency (Kamble et al. 2007; Kocabaş-Ataklı & Yürüm 2013; Hu et al. 2015; Ramos et al. 2016). In addition, this method has the ability to remove small concentrations of heavy metals from a large amount of water solutions. The effectiveness of adsorption mainly depends on the selection of appropriate process conditions such as the mass of sorbent, pH, system temperature and the process duration (Lourie & Gjengedal 2011). Various studies have been conducted using different materials like clay minerals, biomaterials and activated carbon for heavy metal ions removal from water (Chen & Wu 2004; Gupta et al. 2006; Oubagaranadin & Murthy 2010). However, the use of traditional adsorbents has common drawbacks such as small adsorption capacity and low adsorption efficiency (Rao et al. 2007). Consequently, highly efficient adsorbents are necessary to remove the arsenic ions from water solution. Therefore, the development of new adsorbents has become the major interest of researchers in water treatment technology.
Nanoparticles are one of the most popular adsorbents for several pollutants, owing to their features, such as catalytic potential, small size, large surface area and high reactivity (Ali 2012). Carbon nanotubes (CNTs) possess different properties from other materials used, which make them suitable for many applications in electronics, water treatment, optics, nanotechnology and other material science fields (Atieh et al. 2011). However, carbon nanotubes have some limitations ascribed to various flaws in solubility, difficulty in manipulation and aggregation. Hence, surface modification of carbon nanotubes has magnificent affinity by interfacing with other compounds (Thostenson et al. 2001; Sun et al. 2002). The CNTs' surface charge can be enhanced by the oxidative functionalization. However, this method requires the use of strong acids, which is not environmentally friendly. Therefore, the need for an environmentally friendly modification technique is crucial for the widespread application of CNTs (Martinez et al. 2003; Hayyan et al. 2015).
Deep eutectic solvents (DESs) are identified as the liquid combination formed by the complexation of hydrogen bond acceptors (HBAs) and hydrogen bond donors (HBDs) (Abbott et al. 2004; Gorke et al. 2008; Zhang et al. 2012). DESs are new green solvents with many advantages as compared to the ionic liquids (ILs) (Xu et al. 2015). The prime advantages of DESs over conventional ILs are the diversity of physical properties and different molar ratios, ease of synthesis and cheaper price of raw materials. The DESs are synthesized from two or more inexpensive materials, consisting of non-flammable and usually non-toxic components which are able to connect together via hydrogen bonding (Paiva et al. 2014). The components mixture has a lower melting point than the individual compounds (Abbott et al. 2004).
Recently, DESs have been reported in many applications; examples include the uses of ChCl-based DES as a functional additive for starch-based plastics (Leroy et al. 2012), the synthesis of zeolite analogues (Cooper et al. 2004), mediums for the deposition of specific metals in electro and electroless plating of metals (Abbott et al. 2007, 2008), and most recently, in nanotechnology applications (Abo-Hamad et al. 2015).
In general, the adsorption process is considered complicated for heavy metal removal due to the influence of many variables such as contact time, adsorbent dosage, pH and initial heavy metal concentration. The conventional linear method for modelling this kind of process is hectic. On the other hand, the artificial neural networks (ANNs) modelling technique, which is known as a robust black-box modelling tool, is capable of transforming a given dataset into its target outputs. The ANN capability to generalize and learn the behaviour of any non-linear and complex process makes it a powerful tool. ANNs consist of a massive parallel numerical architecture which can solve complicated problems by the assistance of highly connected neurons organized in layers. Recently, the ANNs technique has been used for various engineering applications (Fayaed et al. 2013; Fiyadh et al. 2017). Some studies suggest the NARX neural network is suitable for non-linear systems modelling (McAvoy & Werbos 1992; Çoruh et al. 2014).
The objective of this study was to synthesize a DES by mixing glycerol (Gly) with N,N-diethylethanolammonium chloride, for CNTs functionalization. Subsequently, the DES-functionalized CNTs were utilized as an adsorbent for arsenic removal from water. Furthermore, the NARX neural network was used for modelling and establishing the relationship that exists between the operational variables.
MATERIALS AND METHODOLOGY
This section describes the materials used in DES synthesis, multi-wall carbon nanotubes (MWCNTs) functionalization by DES and the functionalized CNTs characterization. The materials used in the experimental work were hydrochloric acid (36.5–38%), Glycerol (Gly), potassium permanganate (KMnO4), MWCNTs with specification of L 5 μm × D 6–9 nm, >95% carbon, and sodium hydroxide pellets and were all provided by Sigma Aldrich. Arsenic standard solution of 1,000 mg/L and N,N-diethylethanolammonium chloride (DAC) > 99% purity were provided by Merck, Germany.
The DES synthesis was performed by mixing DAC and Gly at a molar ratio of 2:1 HBD:salt (AlOmar et al. 2016b), at 80 °C temperature for 3 hours until the DES became homogenous. The produced DES is referred to as D in this study. The prepared DES was kept in a tightly controlled environment to avoid the effect of humidity.
The first step of the functionalization was to dry the pristine MWCNTs (P-CNTs) at 100 °C overnight. Then, 200 mg of the dried P-CNTs was mixed with 7 mL of KMnO4 and sonicated for 2 h at 65 °C to produce K-CNTs. The functionalization by D-DES was achieved by mixing 7 mL of the prepared D-DES with 200 mg of K-CNTs and sonicating them at 65 °C for 3 h to produce DK-CNTs. After that, the filtration process was performed by washing the functionalized CNTs using distilled water and filtering with a PTFE 0.45 μm membrane until the pH of the filtered water reached neutral. Later, the washed functionalized CNTs were dried overnight at 100 °C before being used for the removal.
Raman spectroscopy, Fourier transform infrared (FTIR) and zeta potential were used for the characterization of P-CNTs, K-CNTs and DK-CNTs adsorbent. The Raman shift was obtained to identify the degree of functionalization using a Renishaw System 2000 Raman Spectrometer. The FTIR was used to identify the surface chemical modification of the adsorbent using FTIR spectroscopy via a PerkinElmer® FTIR spectrometer, USA, with a range of 400–4,000 wave number and four times repetition. Zetasizer (Malvern, UK) was used to recognize the adsorbent partial surface charge by measuring the zeta potential. The surface area of the adsorbent was measured using a fully Automated Gas Sorption System (micromeritics ASAP2020, TRISTAR II 3020 Kr®, USA) (AlOmar et al. 2016a).
Adsorption experiments
The functionalized DK-CNTs adsorbent was used for arsenic removal from water solution. The experiments were conducted with various dosages of DK-CNTs adsorbent (20 to 40 mg), arsenic concentration (1 to 5 mg/L), and with different pH values (3 to 8). The pH of the solution was controlled using NaOH and HCl. A volume of 50 ml of contaminated water was poured into 250-mL flasks, the flasks were shaken at 180 rpm using a mechanical system at room temperature. The number of samples prepared in this study was 156. The arsenic concentrations were tested at different time intervals to study the adsorption equilibrium time using inductively coupled plasma (ICP) with an OES OPTIMA 7000DV PerkinElmer®, USA.
NARX neural network modelling and evaluation indicators
An artificial neural network is a tool that can generate and capture the linear and non-linear relationship between dependent and independent variables (Agami et al. 2009). The NARX neural network is a dynamic network which contains various layers with a back-propagation connection (Chen et al. 1990). The NARX neural network is well known for its high speed of convergence as well as its high degree of generalization.
f is the non-linear function.
u(t) is the network inputs at time t.
y(t) is the network outputs at time t.
nu and ny are the order of inputs and outputs.
During the training, the outputs of the network regress on the target of the actual values as long as they are accessible. Within the training process, the values of the actual target are fed back to the network. This basically results in a better learning and training, and the network acts as feed-forward network which is always steady. The resulting system is known as a NARX network, when the f approach is used with multi-layer perception (Chen et al. 1990). In this paper, two-layers of NARX were used (presented in Figure 1) for the prediction of the adsorption capacity (Qc) of the DK-CNTs. The input layer of the network consists of four inputs (time, adsorbent dosage, pH and initial concentration) and there is one output layer (Qc). In Figure 1, bh is the network bias, wij is the network weight and z is the delay element.
Df(t) = the predicted value.
Da(t) = the actual value.
The RRMSE, MSE, RMSE, MAPE and RE are the indicators that were used for evaluating the model performance. The aim of using different indicators is to confirm the accuracy of the model. All the indicators are based on the obtained results by comparing the error between the actual and predicted results.
RESULTS AND DISCUSSION
In this study, a new adsorbent was prepared (DK-CNTs) and used for arsenic removal from water. The NARX neural network was used for the modelling of the adsorption capacity and different indicators were utilized to evaluate the proposed neural network model. A sensitivity study of the involved parameters in the experimental work, i.e. pH, adsorbent dosage and initial concentration, was implemented. In addition, the adsorption rate order was investigated using three different kinetic models.
Characterization of hybrid material
Studying the electric charge of any adsorbent is crucial due to its influence on the adsorption efficiency. The zeta potential is considered as the electrical potential between the bulk fluid and the surface across the dielectrical layer attached to the suspended particles in a solution. This potential is a source of balancing electrostatic forces that keep the micro- or nano-particles stable in suspension or emulsion. Herein, the absolute zeta potential has increased from 5.5 to −37.6 mV for P-CNTs and DK-CNTs respectively. In addition, the Raman spectra show that the ID/IG (Intensity defect/Intensity graphite) ratio also increased from 1.11 for the P-CNTs to 1.2 for DK-CNTs indicating the presence of new functional groups in sp3 direction resulting from the functionalization effect of D-DES. These functional groups play a significant role in increasing the adsorption capacity of DK-CNTs. FTIR results were in accordance with Raman results. The sp3 direction functionalization was observed by OH stretching appearing in the peaks around 3,400 cm−1. The N–H stretches were in the range of 3,207 cm–1. In addition, the presence of C–Cl bonds may overlap with other CO groups between 600 and 700 cm–1.
It is well known that the surface area of the adsorbent has a huge effect on the adsorption system. The introduction of D-DES as a functionalization agent of CNTs increased the surface area significantly from 123.5 to 200.5 m2/g. This significant increment is reflected in the maximum adsorption capacity of DK-CNTs (AlOmar et al. 2017). The functionalized CNTs showed a better result in the arsenic removal compared to the pristine CNT (AlOmar et al. 2016a).
NARX modelling and performance
Selection of the right NARX network structure with good productivity and accuracy is a complicated task, which includes many points such as the selection of the proper number of hidden layers and the neurons number at the hidden layer. In general, the NARX network structure contains input layer(s), hidden layer(s) and output layer(s). The network selection has been carried out based on the network performance and productivity, using the MSE value during the training phase.
The parameters used in this work were arsenic concentration (1 to 5 mg/L), adsorbent dosage (20 to 40 mg), pH (3 to 8), and contact time until the equilibrium of reaction. One hundred and fifty-six (156) combinations were prepared in lab scale and divided into two sets (training set and testing set), 136 data were used for the training and 20 data were used for the testing. The MATLAB R2014a computational platform was used in the current study to code and optimize the network structure. The optimum hidden layers used for the model creation are two hidden layers with 10 neurons in each hidden layer with one input layer consisting of four nodes and one output layer with one node. The back-propagation training algorithm (trainlm) was selected to update the bias and weight vector values corresponding to the momentum, and the tangent sigmoid transfer function (tansig) was selected as the neuron transfer function for the network. The node numbers at the hidden layer were selected by training and testing the network using different neuron numbers and checking the value of the MSE of the testing set. The network performance depends on the net input, weight of trainlm and tangent sigmoid transfer function, tansig. The minimum value of MSE achieved was 4.75 × 10−4 at the testing set, with correlation coefficient (R2) of 0.9922, which shows a good agreement between the actual and the predicted data, the correlation coefficient plot for the testing set is presented in Figure 2. Different indicators were used to evaluate the trained model such as the relative root mean square error (RRMSE), the root mean square error (RMSE), the mean square error (MSE) and mean absolute percentage error (MAPE). The results of all these indicators are presented in Table 1. The relative error is one of the indicators used in the modelling prediction and it compares the actual values to predicted values. Figure 3 shows the percentage relative error of the model. The maximum relative error value for the NARX model was 5.79%. The best prediction performance depends on the accuracy of the neural network training. This study aimed to get the mathematical approach benefit during the real-time experiment. The NARX model development is becoming a challenge for the real-time experiment. The prepared NARX model used for the sensitivity study involved parameters in the experimental work (initial concentration, adsorbent dosage and pH). Moreover, the kinetic models were applied to the NARX outputs in order to check the model accuracy.
Correlation coefficient of actual and predicted normalized arsenic removal (testing dataset).
Correlation coefficient of actual and predicted normalized arsenic removal (testing dataset).
Evaluation indicators
. | NARX . |
---|---|
MSE | 4.75 × 10−4 |
RMSE | 4.87 × 10−3 |
RRMSE | 2.78 × 10−3 |
MAPE | 2.05 |
. | NARX . |
---|---|
MSE | 4.75 × 10−4 |
RMSE | 4.87 × 10−3 |
RRMSE | 2.78 × 10−3 |
MAPE | 2.05 |
Illustration of the accuracy of the hybrid model based on the testing dataset.
Sensitivity study
Initial concentration
The effect of initial concentration on the adsorption was studied by varying the As3+ concentration from 1 to 5 mg/L, while all the other parameters such as pH (3), adsorbent dosage (20 mg) and contact time (5 min) were kept constant at their nominal levels indicated. The As3+ adsorption percentage is inversely proportional to the As3+ initial concentration. The adsorption capacity increases with increasing the initial metals concentration at a fixed adsorbent dosage. When the initial As3+ concentration was increased from 1 to 3 mg/L, the adsorption capacity also increased from 2.18 to 3.56 mg/g. Whereas, in increasing the initial As3+ concentration from 3 to 5 mg/L the adsorption capacity increased from 3.56 to 5.64 mg/g. This might be attributed to the increase in the driving force of the mass transfer which led to an increase in the uptake capacity of As3+ ions from the water solution. At low concentration, the As3+ ions interact with the adsorbent active sites. On the other hand, at higher As3+ concentration, the adsorbent active sites are saturated and the removal percentage decreases (Banerjee et al. 2016). The obtained experimental data are used for training in the NARX modelling techniques. The NARX model prediction results were found to be in good agreement with the experimental data observation. The experimental and the NARX outputs are presented in Figure 4.
Experimental and NARX outputs as a function of initial concentration.
pH effect
The aqueous solution pH is one of the most important parameters in controlling the adsorption and ion exchange. It is known that the pH can affect the functional groups (i.e. amino groups, phosphate and carboxyl) protonation in the biomass, and also the metal chemistry (i.e. its solubility) (Kazemipour et al. 2008; Witek-Krowiak et al. 2011). The effect of pH was examined by varying its value from 3 to 8, and fixing all the other involved parameters such as, adsorbent dosage (20 mg), initial concentration (1 mg/L) and contact time (50 min) at their nominal levels indicated.
It can be observed in Figure 5 that the adsorption capacity of DK-CNTs was increased with the increase of pH value until pH 6. Thereafter, the adsorption capacity became almost steady with increasing pH from 6 to 8. This increase may be due to the presence of the negative charge of the oxygen-containing functional groups, such as the carboxylic group, and the concentration of the negative electron charges enhanced by the presence of OH– in the solution. This pattern of negative charge intensive distribution on the surface of the adsorbent may be responsible for metal binding. It is well known that at pH greater than 7.0, the dominant species of As3+ are As (OH)+. This complexation may occur due to the extensive presence of OH− at this pH level resulting in a precipitation form, this phenomenon was also stated by (Gupta et al. 2011). In addition, the decreasing of H+ plays a significant role in the mechanism of As3+ adsorption due to the decrease of competition for the active sites of the adsorbent. The NARX technique was used for the modelling and prediction of the obtained data from the experimental work. The NARX-based ANN model prediction results showed good agreement with the experimental result trend. The NARX outputs and the experimental results as a function of pH versus the uptake capacity are presented in Figure 5.
Adsorbent dosage study
The effect of adsorbent dosage was studied by varying the adsorbent dosage from 20 to 40 mg under a fixed time of 90 min, initial concentration of 3 mg/L and pH of 6. It can be seen from Figure 6 that the adsorption capacity of As3+ ions decreased with increasing the DK-CNTs adsorbent dosage value. The adsorption capacity for 20 mg dosage was 7.56 mg/g then, as the DK-CNTs adsorbent was increased to 30 and 40 mg the adsorption capacity decreased to 6.23 and 5.89 mg/g, respectively. The decrease in the arsenic uptake capacity as the adsorbent dosage increases might be attributed to the increase of more active sites due to the addition of adsorbent surface area which was also reported by Das et al. (2014). The obtained experimental data were used in training and prediction by using the NARX modelling techniques. The NARX model prediction was found to be satisfactory for the experimental data observation. The experimental and predicted outputs of the NARX are presented in Figure 6.
Adsorption kinetics study
Three adsorption kinetic models were implemented in this study to investigate the As3+ adsorption rate and mechanism as well as the solute removal rate (Ayoob et al. 2007). The intraparticle diffusion, pseudo-first-order and pseudo-second-order models were used in this work. The kinetic studies were achieved at different pH values (3, 5.5 and 8) and initial concentration (1 and 3 mg/L) with different adsorbent dosage (20 and 30 mg) and agitation speed of 180 rpm. The correlation coefficient (R2) was used as the conformity indicator of the kinetic models between the experiment and the predicted adsorption values for each kinetic model. The pseudo-second-order model best described the adsorption kinetics of the system and the results are presented in Figure 7, as time over adsorption capacity (T/Q). The results of the intraparticle diffusion and pseudo-first-order models are presented in Table 2. The kinetics study results reveal that the amount of DK-CNTs adsorbent and its concentration are associated with the rate determination step, which indicates that the rate limiting step involves chemisorption. The same behaviour was reported elsewhere (Veličković et al. 2013).
Adsorption kinetics and correlation coefficient
Dose mg . | PH . | C0mg/L . | Pseudo-first-order ln (qe–qt) vs time (t) . | Pseudo-second-order (t/qt vs t) . | Intraparticle (qt vs t0.5) . | |||
---|---|---|---|---|---|---|---|---|
Experimental R2 . | NARX output R2 . | Experimental R2 . | NARX outputs R2 . | Experimental R2 . | NARX outputs R2 . | |||
20 | 3 | 1 | 0.8925 | 0.9052 | 0.9973 | 0.9979 | 0.7471 | 0.7172 |
20 | 8 | 1 | 0.5483 | 0.5997 | 0.9983 | 0.9988 | 0.6496 | 0.6244 |
20 | 3 | 3 | 0.656 | 0.6819 | 0.9711 | 0.963 | 0.7931 | 0.7311 |
20 | 5.5 | 3 | 0.7858 | 0.7017 | 0.9562 | 0.9897 | 0.6775 | 0.6884 |
30 | 5.5 | 3 | 0.8939 | 0.8682 | 0.9978 | 0.9909 | 0.5338 | 0.5848 |
30 | 8 | 3 | 0.6758 | 0.6069 | 0.9905 | 0.9945 | 0.7542 | 0.795 |
30 | 3 | 1 | 0.756 | 0.7888 | 0.9994 | 0.9974 | 0.5921 | 0.5229 |
30 | 5.5 | 1 | 0.8102 | 0.8608 | 0.9999 | 0.9994 | 0.674 | 0.6838 |
Dose mg . | PH . | C0mg/L . | Pseudo-first-order ln (qe–qt) vs time (t) . | Pseudo-second-order (t/qt vs t) . | Intraparticle (qt vs t0.5) . | |||
---|---|---|---|---|---|---|---|---|
Experimental R2 . | NARX output R2 . | Experimental R2 . | NARX outputs R2 . | Experimental R2 . | NARX outputs R2 . | |||
20 | 3 | 1 | 0.8925 | 0.9052 | 0.9973 | 0.9979 | 0.7471 | 0.7172 |
20 | 8 | 1 | 0.5483 | 0.5997 | 0.9983 | 0.9988 | 0.6496 | 0.6244 |
20 | 3 | 3 | 0.656 | 0.6819 | 0.9711 | 0.963 | 0.7931 | 0.7311 |
20 | 5.5 | 3 | 0.7858 | 0.7017 | 0.9562 | 0.9897 | 0.6775 | 0.6884 |
30 | 5.5 | 3 | 0.8939 | 0.8682 | 0.9978 | 0.9909 | 0.5338 | 0.5848 |
30 | 8 | 3 | 0.6758 | 0.6069 | 0.9905 | 0.9945 | 0.7542 | 0.795 |
30 | 3 | 1 | 0.756 | 0.7888 | 0.9994 | 0.9974 | 0.5921 | 0.5229 |
30 | 5.5 | 1 | 0.8102 | 0.8608 | 0.9999 | 0.9994 | 0.674 | 0.6838 |
The NARX neural network technique was used for modelling and prediction of the obtained data from the experimental work. The three kinetics models used for modelling the experimental data were also applied on the NARX outputs. The pseudo-second-order model best described the adsorption kinetics of this study as compared to the intraparticle diffusion and pseudo-first-order models. The results of the kinetics study are presented in Table 2. The NARX model shows good agreement with the experimental work, which shows that the NARX model has high accuracy.
CONCLUSION
A new adsorbent was developed using a DES system as the functionalization agent of carbon nanotubes. The amount of arsenic removal increased with the increasing of contact time, pH and initial concentration, whereas the arsenic removal decreased with increasing adsorbent dosage. Comparing the experiment results with the NARX model outputs, it can be concluded that the NARX model is able to predict the amount of arsenic removal from water sufficiently with acceptable error. The minimum value of MSE achieved was 4.75 × 10−4 at the testing set, with correlation coefficient (R2) of (0.9922), which shows a good agreement between the actual and the predicted data. Other indicators were also used such as the RMSE (4.87 × 10−3), RRMSE (2.78 × 10−3) and MAPE (2.05). All these indicators confirmed the high accuracy of the NARX model.
ACKNOWLEDGEMENTS
The authors express their thanks to the University of Malaya for funding this research. UMRG (RP044D-17AET) and (RP025A-18SUS).