Open Access

ARTICLE

# A Novel Fractional Dengue Transmission Model in the Presence of Wolbachia Using Stochastic Based Artificial Neural Network

1 Department of Mathematics, COMSATS University Islamabad, Islamabad, 45550, Pakistan

2 Department of Mathematics, Cankaya University, Ankara, 06790, Turkey

3 Institute of Space Sciences, Magurele-Bucharest, 077125, Romania

4 Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, 40402, Taiwan

5 Department of Computer Science and Mathematics, Lebanese American University, Beirut, 135053, Lebanon

6 Department of Mathematics, Mathematics Research Center, Near East University, Nicosia/Mersin, 99138, Turkey

* Corresponding Author: Shumaila Javeed. Email:

*Computer Modeling in Engineering & Sciences* **2024**, *139*(2), 1217-1238. https://doi.org/10.32604/cmes.2023.029879

**Received** 12 March 2023; **Accepted** 02 August 2023; **Issue published** 29 January 2024

## Abstract

The purpose of this research work is to investigate the numerical solutions of the fractional dengue transmission model (FDTM) in the presence of Wolbachia using the stochastic-based Levenberg-Marquardt neural network (LM-NN) technique. The fractional dengue transmission model (FDTM) consists of 12 compartments. The human population is divided into four compartments; susceptible humans (*S*), exposed humans (

_{h}*E*), infectious humans (

_{h}*I*), and recovered humans (

_{h}*R*). Wolbachia-infected and Wolbachia-uninfected mosquito population is also divided into four compartments: aquatic (eggs, larvae, pupae), susceptible, exposed, and infectious. We investigated three different cases of vertical transmission probability (

_{h}*η*), namely when Wolbachia-free mosquitoes persist only (

*η*= 0.6), when both types of mosquitoes persist (

*η*= 0.8), and when Wolbachia-carrying mosquitoes persist only (

*η*= 1). The objective of this study is to investigate the effectiveness ofWolbachia in reducing dengue and presenting the numerical results by using the stochastic structure LM-NN approach with 10 hidden layers of neurons for three different cases of the fractional order derivatives (

*α*= 0.4, 0.6, 0.8). LM-NN approach includes a training, validation, and testing procedure to minimize the mean square error (MSE) values using the reference dataset (obtained by solving the model using the Adams-Bashforth-Moulton method (ABM). The distribution of data is 80% data for training, 10% for validation, and, 10% for testing purpose) results. A comprehensive investigation is accessible to observe the competence, precision, capacity, and efficiency of the suggested LM-NN approach by executing the MSE, state transitions findings, and regression analysis. The effectiveness of the LM-NN approach for solving the FDTM is demonstrated by the overlap of the findings with trustworthy measures, which achieves a precision of up to 10

^{−4}.

## Keywords

Nomenclature

LM | Levenberg-Marquardt |

FSTM | Fractional dengue transmission model |

ANN | Artificial neural network |

MSE | Mean square error |

ABM | Adams-Bashforth-Moulton |

Dengue fever is widespread, with an approximated 3.9 billion people at threat and 390 million new dengue infections estimated each year [1]. The female Aedes aegypti and Aedes albopictus mosquitoes are the main vectors spreading dengue infection [2]. The virus that causes dengue has four serotypes: DEN1, DEN2, DEN3, and DEN4, with a probable fifth strain recently discovered [3]. Recent progress has been made in the development of dengue vaccines that provide protection against four strains of the virus. For example, the first vaccine called Dengvaxia has been approved in many countries. Other dengue vaccine candidates like TAK-003 and TV003/TV005 are also in development [4–6].

Dengue is transmitted to people through the bite of a dengue-infected female mosquito. The extrinsic incubation period (EIP) is the time it takes for the dengue virus to multiply in a mosquito’s body before it can reach the salivary glands and be transmitted to a human through a bite. Factors such as temperature, humidity, and rainfall affect the length of the EIP. A comprehensive approach that includes these factors can help to reduce the EIP of dengue virus in mosquitoes and ultimately reduce the transmission of the disease. Only mosquitoes that survive the EIP are capable of transmitting the virus to humans [7].

Wolbachia, a bacterium, is present in certain mosquito species and about 60% of all groups of insects. The potential use of Wolbachia-infected mosquitoes as an alternative approach to control dengue transmission. While vaccines and insecticides are effective, they have limitations such as incomplete protection and adverse effects on the environment. In contrast, Wolbachia-infected mosquitoes have shown promise in reducing dengue transmission in field trials by preventing the virus from being transmitted to humans. This method appears to have fewer environmental impacts and a lower likelihood of insecticide-resistant mosquito populations developing [8]. Wolbachia is transferred through vertical transmission and changes the reproductive phenotype of potentially infected insects, providing the bacterium a reproductive advantage over uninfected insects. Depending on the insect species and Wolbachia strain, the reproductive phenotype displayed might lead to feminization, male death, parthenogenesis, or, most frequently, cytoplasmic incompatibility. Cytoplasmic incompatibility (CI) is the characteristic of Wolbachia that is critical for the invasion of population and perseverance because it provides a reproductive advantage to female mosquitoes that carry Wolbachia [9]. Fractional models can capture the behavior of a system when classical differential equations fail, especially when memory effects or long-range interactions are present. To obtain the results in non-integer order we use the fractional model. Compared to ordinary differential equations, fractional differential equations can describe more complex dynamics and can be formulated in various ways, such as in terms of Caputo or Riemann-Liouville fractional derivatives. Fractional mathematical models can help to better understand the dynamics of Wolbachia transfer and how it affects the invasion of dengue.

The theory of fractional derivatives has a wide range of applications in medicine, finance, engineering, physics, and many other fields [10–12]. Many researchers analyzed a wide range of complicated phenomena based on the fractional derivatives, Kilicman et al. [13] worked on the SIR fractional model for dengue transmission, Jajarmi et al. [14] investigated a fractional mathematical model for dengue fever, Qureshi et al. [15] analyzed the derivatives of fractional-order. Zafar et al. [16] analyzed the fractional model for the Human papillomavirus (HPV). The Caputo derivatives are easier to utilize than the other fractional derivative concepts. The definition of Caputo derivative as defined in [11] is:

where

An artificial neural network (ANN) is a method for simulating biological neural activity using mathematics. The essential mathematical principles of an artificial neural network employing an artificial neuron are multiplication, summation, and activation. At the start of the neuron, the input values are multiplied by the specified weights. The sum function, which adds all weighted inputs and biases, is located in the middle layer of the artificial neuron. The summation of previously weighted inputs and bias is sent via an activation function, also known as a transfer function, at the end of an artificial neuron [17].

In ANN, an activation function is utilized to provide nonlinear features. In a neural network, the yn are input variables, wn are the weights, and fn shows the output in Fig. 1. Using a multilayer perceptron (MLP), the number of hidden units can be optimized. The number mentioned in the hidden layer shows how many neurons were utilized in the network [18].

In literature, many fractional-order mathematical models along with artificial neural network approach models for real-world problems have been proposed utilizing the ANNs approach. Pornsawad et al. [19] analyze fractional-order mathematical differential equations for the backward heat problems by developing neural networks procedure. Qu et al. [20] analyzed the stability of the Caputo fractional-order derivative by utilizing the cellular neural network approach. Sabir et al. [21] analyzed the fractional order mathematical model of computer viruses by employing the Levenberg approach. Ghazizadeh et al. [22] studied the fractional order inverse problem with the estimation of relaxation time by utilizing the neural network approach. Souayeh et al. [23] utilized the Bayesian regularization approach to understand the Allee effect on the food chain supply. Umar et al. [24] investigated the neural network approach for analyzing the behavior of the nonlinear coronavirus model. Noinang et al. [25] worked on a system of magnetohydrodynamic (MHD) past over an inclined plate by utilizing the artificial neural network approach. Zhuang et al. [26] analyzed parameter estimation of inverse problem heat conduction problem in a composite medium utilizing Levenberg-Marquardt. Sabir et al. [27] studied nonlinear singular second-order delay differential equation utilizing the ANN approach.

The aim of this work was to develop and numerically solve a fractional dengue transmission model (FDTM) with Wolbachia using neural network. In the present model, we incorporated Wolbachia-infected mosquitoes along with Wolbachia-free mosquitoes and humans. The impact of vertical transmission probability on the persistence of dengue transmission in humans and mosquitoes has been analyzed using the Levenberg-Marquardt neural network (LM-NN) approach. The convergence and accuracy of the proposed LM-NN technique have been checked through mean square error, regression analysis plots, and state transition results. The results obtained from the proposed LM-NN approach for FDTM have been compared with the reference solutions obtained using the Adams-Bashforth-Moulton (ABM) method and they were found to be in good agreement.

A fractional dengue transmission model in the presence of Wolbachia is developed to observe the dynamics of dengue between humans and mosquitoes. The model consists of 12 compartments, namely, susceptible human

The reproduction rate of Wolbachia-carrying and Wolbachia-free mosquitoes are denoted by

The model flowchart is given in Fig. 2.

All the parameters that appear in model (1) are defined in Table 1. The FDTM model in the presence of Wolbachia is given below:

According to the model, a susceptible human becomes exposed to a bite from an infected mosquito. The biting rates of Wolbachia-infected and Wolbachia-free mosquitoes are

The exposed Wolbachia-free and Wolbachia-carrying mosquitoes move to the infected class after biting the dengue-infected human with rate

In this section, we explained the proposed Levenberg-Marquardt neural network (LM-NN) method that is used to solve the FDTM model (1). The Pseudocode of the proposed LM-NN technique is also given in Table 2. The construction of neuron is given in Fig. 3. The details of steps involved in utilizing the LM-NN technique to solve the FDTM model in the presence of Wolbachia given in Fig. 4.

The novel properties of the proposed research work is as follows:

• Levenberg-Marquardt neural network (LM-NN) is based on a backpropagation process with a novel methodology or design is created for solving the fractional-order model (1) dealing with Wolbachia carrying and Wolbachia-free mosquitoes with humans.

• For the fractional dengue transmission model (FDTM), a stochastic-based back propagated numerical Levenberg-Marquardt neural network (LM-NN) procedure is utilized to obtain the best-approximated solution.

• Adams-Bashforth-Moulton (ABM) numerical method has been employed to generate the reference dataset.

• Numerical results obtained from the proposed LM-NN technique show good agreement with the reference solutions achieved using the ABM method.

• The reliability and consistency of the proposed LM-NN approach are investigated through statistical analysis.

• The best performance of the proposed LM-NN technique is analyzed using regression analysis (RA) plots, mean square error (MSE) graphs, fitness curves, and state transitions (STs) for different classes.

• Error analysis performed for different classes of FDTM is used to check the accuracy of the proposed LM-NN technique.

This section describes the numerical simulations based on the outcomes of three scenarios of the nonlinear fractional dengue transmission model (FDTM) performed with the suggested LM-NN approach. The simulations were conducted on scenarios with only Wolbachia-free mosquitoes present, both Wolbachia-free and Wolbachia-carrying mosquitoes present, and only Wolbachia-carrying mosquitoes present. The mathematical structure for each case is given as:

4.1 Case-1: Only Wolbachia-Free Mosquitoes Persist

For this case, we consider the following values of the parameters

4.2 Case-2: Both Type of Mosquitoes Wolbachia-Free and Wolbachia-Carrying Persist

For this case, the parameter values being examined are

4.3 Case-3: Only Wolbachia-Carrying Mosquitoes Persist

Similarly, after putting the values of parameters in the model (1)

The numerical results of the FDTM model is obtained through LM-NN technique which is implemented in MATLAB using the built-in package ‘nftool’ together with 10 neurons. To train the model using the Levenberg method, we used 80% of the data for training purposes, 10% of the data is used for validation and the remaining 10% of the data is used for testing. Fig. 3 shows the structure of a single neuron used for the best solution of the FDTM. The MSE performance for training, validation, and testing along with the number of epochs and gradient are presented in Table 3.

The mean square error results obtained from the proposed neural network approach for each scenario of the FDTM are presented in Fig. 5. At epochs 176, 426, and 146 the mean square error is very small and the best performance values are

Figs. 6–8 present the graphs of gradient, Mu, and validation performances of the proposed LM-NN to solve the fractional dengue transmission model (FDTM). The values of gradient, and Mu, for Case I are

Regression is a graphical representation of the precision of the LM-NN solution to the reference values for training, validation, and testing points separately. A straight line represents an available reference solution in these plots, while dots or small circles represent LM-NN results. The numerical value of regression can also be used to assess the computation’s accuracy and precision. R = 1 indicates that the LM-NN solution is very near to the reference values, whereas R = 0 indicates that the relationship between the reference and the LM-NN solution is very weak. The regression plots are shown in Figs. 12–14. From the regression plots, we observe that R value is 1 for all the three different cases. The testing, validation, and training plots demonstrate that the developed LM-NN technique is efficient for solving the FDTM model.

Fig. 15 illustrates the comparison of the numerical results of model (1) using the ABM approach to the solution of the FDTM using the LM-NN technique. Fig. 15a shows the comparison of the

Error plots were generated for three different scenarios involving the classes

The purpose of this research study was to develop the fractional dengue transmission model (FDTM) in the presence of Wolbachia-carrying mosquitoes. The model incorporated the vertical transmission of Wolbachia in mosquitoes. In this study, we proposed the numerical solution of the fractional dengue transmission model (FDTM) using the stochastic-based Levenberg-Marquardt neural networks (LM-NN) technique. The FDTM model was solved using the Adams-Bashforth-Moulton (ABM) technique to create a reference dataset. The created dataset was used for training (80%), testing (10%), and validation (10%) of the proposed LM-NN technique with 10 hidden layers of neurons. The graphs of mean square error were presented for three different cases, when Wolbachia-free mosquitoes persist only

Acknowledgement: The authors are thankful to the anonymous reviewers for improving this article.

Funding Statement: The authors received no specific funding for this study.

Author Contributions: The authors confirm their contribution to the paper as follows: study conception and design: SJ, IA, and DB; methodology: ZF; analysis and interpretation of results: SJ, IA, ZF, and DB; draft manuscript preparation: SJ, IA, ZF, and DB; validation: SJ and DB; review and editing: IF and ZF. All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials: Not applicable.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

1. Bhatt, S., Gething, P. W., Brady, O. J., Messina, J. P., Farlow, A. W. et al. (2013). The global distribution and burden of dengue. Nature, 496(7446), 504–507. [Google Scholar] [PubMed]

2. Adekunle, A. I., Meehan, M. T., McBryde, E. S. (2019). Mathematical analysis of a wolbachia invasive model with imperfect maternal transmission and loss of wolbachia infection. Infectious Disease Modelling, 4, 265–285. [Google Scholar] [PubMed]

3. Normile, D. (2013). First new dengue virus type in 50 years. Health. https://doi.org/10.1126/article.24052 [Google Scholar] [CrossRef]

4. Vaccine, D. (2018). WHO position paper–September 2018. The Weekly Epidemiological Record, 93, 457–476. [Google Scholar]

5. de Silva, A. (2023). Safety of dengue vaccine? Clinical Infectious Diseases, 76(2), 371–372. [Google Scholar] [PubMed]

6. Tricou, V., Low, J. G., Oh, H. M., Leo, Y. S., Kalimuddin, S. et al. (2020). Safety and immunogenicity of a single dose of a tetravalent dengue vaccine with two different serotype-2 potencies in adults in singapore: A phase 2, double-blind, randomised, controlled trial. Vaccine, 38(6), 1513–1519. [Google Scholar] [PubMed]

7. Tjaden, N. B., Thomas, S. M., Fischer, D., Beierkuhnlein, C. (2013). Extrinsic incubation period of dengue: Knowledge, backlog, and applications of temperature dependence. PLoS Neglected Tropical Diseases, 7(6), e2207. [Google Scholar] [PubMed]

8. Dorigatti, I., McCormack, C., Nedjati-Gilani, G., Ferguson, N. M. (2018). Using wolbachia for dengue control: Insights from modelling. Trends in Parasitology, 34(2), 102–113. [Google Scholar] [PubMed]

9. Turelli, M. (2010). Cytoplasmic incompatibility in populations with overlapping generations. Evolution, 64(1), 232–241. [Google Scholar] [PubMed]

10. Fatmawati, Jan, R., Khan, M. A., Khan, Y., Ullah, S. et al. (2020). A new model of dengue fever in terms of fractional derivative. Mathematical Biosciences and Engineering, 17(5), 5267–5288. [Google Scholar] [PubMed]

11. Ahmad, S., Javeed, S., Ahmad, H., Khushi, J., Elagan, S. et al. (2021). Analysis and numerical solution of novel fractional model for dengue. Results in Physics, 28, 104669. [Google Scholar]

12. Srivastava, H., Saad, K. M., Khader, M. (2020). An efficient spectral collocation method for the dynamic simulation of the fractional epidemiological model of the Ebola virus. Chaos, Solitons & Fractals, 140, 110174. [Google Scholar]

13. Hamdan, I. N., Kilicman, A. (2018). A fractional order sir epidemic model for dengue transmission. Chaos, Solitons & Fractals, 114, 55–62. [Google Scholar]

14. Jajarmi, A., Arshad, S., Baleanu, D. (2019). A new fractional modelling and control strategy for the outbreak of dengue fever. Physica A: Statistical Mechanics and its Applications, 535, 122524. [Google Scholar]

15. Qureshi, S., Atangana, A. (2019). Mathematical analysis of dengue fever outbreak by novel fractional operators with field data. Physica A: Statistical Mechanics and its Applications, 526, 121127. [Google Scholar]

16. Zafar, Z. U. A., Hussain, M., Inc, M., Baleanu, D., Almohsen, D. et al. (2022). Fractional-order dynamics of human papillomavirus. Results in Physics, 34, 105281. [Google Scholar]

17. Suzuki, K. (2011). Artificial neural networks: Methodological advances and biomedical applications,pp. 1–376. UK: InTech. [Google Scholar]

18. Alhakami, H., Umar, M., Sulaiman, M., Alhakami, W., Baz, A. (2022). A numerical study of the dynamics of vector-born viral plant disorders using a hybrid artificial neural network approach. Entropy, 24(11), 1511. [Google Scholar] [PubMed]

19. Pornsawad, P., Böckmann, C., Panitsupakamon, W. (2022). The levenberg–marquardt regularization for the backward heat equation with fractional derivative. Electronic Transactions on Numerical Analysis, 57, 67–79. [Google Scholar]

20. Qu, H., Zhang, T., Zhou, J. (2020). Global stability analysis of S-asymptotically ω-periodic oscillation in fractional-order cellular neural networks with time variable delays. Neurocomputing, 399, 390–398. [Google Scholar]

21. Sabir, Z., Raja, M. A. Z., Mumtaz, N., Fathurrochman, I., Sadat, R. et al. (2022). An investigation through stochastic procedures for solving the fractional order computer virus propagation mathematical model with kill signals. Neural Processing Letters, 55(2), 1783–1797. [Google Scholar]

22. Ghazizadeh, H. R., Azimi, A., Maerefat, M. (2012). An inverse problem to estimate relaxation parameter and order of fractionality in fractional single-phase-lag heat equation. International Journal of Heat and Mass Transfer, 55(7–8), 2095–2101. [Google Scholar]

23. Souayeh, B., Sabir, Z., Hdhiri, N., Al-Kouz, W., Alam, M. W. et al. (2022). A stochastic bayesian regularization approach for the fractional food chain supply system with allee effects. Fractal and Fractional, 6(10), 553. [Google Scholar]

24. Umar, M., Sabir, Z., Raja, M. A. Z., Shoaib, M., Gupta, M. et al. (2020). A stochastic intelligent computing with neuro-evolution heuristics for nonlinear sitr system of novel COVID-19 dynamics. Symmetry, 12(10), 1628. [Google Scholar]

25. Noinang, S., Sabir, Z., Javeed, S., Raja, M. A. Z., Ali, D. et al. (2022). A novel stochastic framework for the mhd generator in ocean. Computers, Materials & Continua, 73(2), 3383–3402. [Google Scholar]

26. Zhuang, Q., Yu, B., Jiang, X. (2015). An inverse problem of parameter estimation for time-fractional heat conduction in a composite medium using carbon–carbon experimental data. Physica B: Condensed Matter, 456, 9–15. [Google Scholar]

27. Sabir, Z., Guirao, J. L., Saeed, T. (2021). Solving a novel designed second order nonlinear Lane–Emden delay differential model using the heuristic techniques. Applied Soft Computing, 102, 107105. [Google Scholar]

28. Ndii, M. Z., Hickson, R. I., Allingham, D., Mercer, G. (2015). Modelling the transmission dynamics of dengue in the presence of wolbachia. Mathematical Biosciences, 262, 157–166. [Google Scholar] [PubMed]

29. Junsawang, P., Zuhra, S., Sabir, Z., Raja, M. A. Z., Shoaib, M. et al. (2022). Numerical simulations of vaccination and wolbachia on dengue transmission dynamics in the nonlinear model. IEEE Access, 10, 31116–31144. [Google Scholar]

## Cite This Article

**APA Style**

*Computer Modeling in Engineering & Sciences*,

*139*

*(2)*, 1217-1238. https://doi.org/10.32604/cmes.2023.029879

**Vancouver Style**

**IEEE Style**

*Comput. Model. Eng. Sci.*, vol. 139, no. 2, pp. 1217-1238. 2024. https://doi.org/10.32604/cmes.2023.029879

This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.