images
Computer Modeling in Engineering & Sciences
images

DOI: 10.32604/cmes.2021.015472

ARTICLE

Modeling Dysentery Diarrhea Using Statistical Period Prevalence

Fouad A. Abolaban*

Department of Nuclear Engineering, College of Engineering, King Abdulaziz University, Jeddah, 21589, Saudi Arabia
*Corresponding Author: Fouad A. Abolaban. Email: fabolaban@kau.edu.sa
Received: 21 December 2020; Accepted: 07 February 2021

Abstract: Various epidemics have occurred throughout history, which has led to the investigation and understanding of their transmission dynamics. As a result, non-local operators are used for mathematical modeling in this study. Therefore, this research focuses on developing a dysentery diarrhea model with the use of a fractional operator using a one-parameter Mittag–Leffler kernel. The model consists of three classes of the human population, whereas the fourth one belongs to the pathogen population. The model carefully deals with the dimensional homogeneity among the parameters and the fractional operator. In addition, the model was validated by fitting the actual number of dysentery diarrhea infected cases covering 52 weeks in 2017, which occurred in Ethiopia. The biological parameters were fitted, and fractional order ν was optimized. The basic reproduction numbers R0 were 1.7031 and 1.9581, which correspond to the fractional and integer-order models, respectively. The fractional model showed smaller discrepancies compared to the integer-order model when the models were fitted and validated with the actual cases of infected humans. Qualitative theory for the existence and uniqueness of the solutions was extensively explored, and the stability of equilibria was investigated. The most influential parameters were found via forward sensitivity indices. This multidisciplinary research shows and explains the results of numerical simulations applied to a biological issue.

Keywords: Dysentery diarrhea model; Mittag–Leffler kernel; fractional operator

1  Introduction

The literature provides mathematical models for the transmission of infectious diseases. These models play a significant role in quantifying and evaluating the effective control and preventive measures of infectious diseases [1–3]. Furthermore, mathematical modeling has been used in several ways as a versatile and effective way of studying the dynamics of infectious disease transmission. This can include the classic susceptible, infected, and recovered (SIR) model or more advanced models [4]. Mathematical analysis and numerical simulations can be collectively used for the development and evaluation of persuasive control measures.

It is popularly known that mathematical models can predict the emergence of infectious diseases and epidemics, which are beneficial for public health planning and initiatives. By using compartmental models as a simple mathematical structure, the complex dynamics of epidemiological processes can be examined [5]. These compartmental models divide the population into two distinct health categories. The first category is depicted by S, and involves those susceptible to pathogen infection, while the second, denoted by I, involves the pathogen-infected humans. The manner in which these two populations interact is based on phenomenological assumptions used to develop the models. Ordinary differential equations (ODEs) are typically used to develop these models.

Additionally, other populations denoted by R, which is the image of the immune/removed/ recovered compartment, are considered to make these models more practical. A significant challenge here is to obtain sufficient parameters for a particular disease, which would determine the factors affecting potential control measures, such as medication or vaccination. The crucial question is about the execution of such measures from an optimal viewpoint. Several notable attempts have recently been made to introduce this research program for various diseases with integer compartmental models [6–11].

Over the past few decades, many scientists have shown that fractional models can effectively represent natural phenomena compared to integer-order differential equations. Therefore, fractional calculus has gained more importance and popularity for modeling realistic cases, especially memory effects [12–14]. Due to this particular function, various fractional operators have been developed to accurately model the memory effects of various types of diseases [15–22]. Nevertheless, further research is required to explain such complex dynamics. Classical fractional models with singular operators cannot effectively model the non-locality of real-world applications. To overcome this challenge, we investigated and examined a new fractional version of an epidemiological model for the dysentery diarrhea involving the ABC operator, known to have a non-singular kernel with memory effects.

2  Formulation of the Dysentery Diarrhea System

This section demonstrates the formulation of the dysentery diarrhea disease in terms of a deterministic model based on a nonlinear system of ODEs over a finite time interval [0, T], 0 < t < T. The model comprises four population classes. Three classes are reserved for samples from the human population, including those that can be infected with the disease called susceptible S(t). In addition, those that transmit the disease are denoted by I(t), and those that successfully recover from the disease are denoted by R(t). The other class is the dynamics for the pathogen population (concentration of Shigella dysenteriae), denoted by B(t).

The model is designed based on the following assumptions:

•   Transmission of the dysentery diarrhea disease occurs through multiple pathways.

•   There is a homogeneously mixed population

•   Π is the rate of recruitment of those susceptible either by birth or immigration

•   Standard incidence is assumed in the human to human interaction

•   Logistic phenomenon is taken into consideration in the human to environment interaction, which is represented below:

Πh=βhINΠB=βeBB+C,(1)

where C stands for the Shigella concentration that causes a 25% likelihood of getting the disease

•   βh and βe stand for rates of ingesting Shigella from human to human interaction and through a contaminated environment, respectively. βh=ap (a = contact rate and p = probability of disease transmission per contact)

•   After losing immunity, individuals return to S(t) at a rate of α

•   Infected individuals cause concentration of Shigella at a rate of ε

•   Shigella population dies at a rate of σ

•   Rate of recovery of infected humans is γ

•   Natural mortality rate of humans is μ where death due to dysentery diarrhea disease occurs at a rate of d

•   Non-negativity is assumed for all biological parameters introduced within the model

Hence, after incorporating all the above assumptions and considering the Atangana–Baleanu differential operator taken in the Caputo sense [15], we have obtained the following coupled nonlinear system of ordinary differential equations:

ABC(ν)1-ν∫0tS′(ζ)Eν[-ν1-ν(t-ζ)ν]dζ=Πν+ανR(t)-(Πh+ΠB+μν)S(t),ABC(ν)1-ν∫0tI′(ζ)Eν[-ν1-ν(t-ζ)ν]dζ=(Πh+ΠB)S(t)-(μν+γν+dν)I(t),ABC(ν)1-ν∫0tR′(ζ)Eν[-ν1-ν(t-ζ)ν]dζ=γνI(t)-(μν+αν)R(t),ABC(ν)1-ν∫0tB′(ζ)Eν[-ν1-ν(t-ζ)ν]dζ=ενI(t)-σνB(t),(2)

subject to the following initial conditions:

S(0)=S0≥0,I(0)=I0≥0,R(0)=R0≥0,B(0)=B0≥0,

where Πh=βhνI(t)N(t),ΠB=βeνB(t)B(t)+C, and N(t) = S(t)+I(t)+R(t).

3  Existence and Uniqueness

This section presents the existence and uniqueness of solutions of the proposed model using the techniques of fixed point theory. Here, we denote E=C([0,T],ℝ), the Banach space of all continuous real-valued function equipped with the norm defined by:

∥(S,I,R,B)∥=∥S(t)∥+∥I(t)∥+∥R(t)∥+∥B(t)∥,

where

|S∥= supt∈[0,T]|S(t)|,∥I∥= supt∈[0,T]|I(t)|,∥R∥= supt∈[0,T]|R(t)|,∥B∥= supt∈[0,T]|B(t)|.

Thus, the proposed fractional model takes the forms shown below:

S(t)-S(0)=ABCI0+ν{Πν+ανR(t)-(Πh+ΠB+μν)S(t))},I(t)-I(0)=ABCI0+ν{(Πh+ΠB)S(t)-(μν+γν+dν)I(t)},R(t)-R(0)=ABCI0+ν{γνI(t)-(μν+αν)R(t)},B(t)-B(0)=ABCI0+ν{ενI(t)-σνB(t)}.(3)

We then obtain:

S(t)=S(0)+1-νN(ν)F1(t,S(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F1(x,S(x))dx,I(t)=I(0)+1-νN(ν)F2(t,I(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F2(x,I(x))dx,R(t)=R(0)+1-νN(ν)F3(t,R(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F3(x,R(x))dx,B(t)=B(0)+1-νN(ν)F4(t,B(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F4(x,B(x))dx,(4)

where,

F1(t,S(t))=Πν+ανR(t)-(Πh+ΠB+μν)S(t),F2(t,I(t))=(Πh+ΠB)S(t)-(μν+γν+dν)I(t),F3(t,R(t))=γνI(t)-(μν+αν)R(t),F4(t,B(t))=ενI(t)-σνB(t).(5)

The kernels in Eq. (5) satisfy the Lipschitz condition for 0≤Mi<0, i=1,2,⋯4. When S(t) and S*(t) are two functions, we get:

∥F1(t,S(t))-F1(t,S*(t))∥=∥Πν+ανR(t)-(Πh+ΠB+μν)S(t)-(Πν+ανR(t)-(Πh+ΠB+μν)S*(t))∥=∥(Πh+ΠB+μν)(S*(t)-S(t))∥≤(Πh+ΠB+μν)∥S(t)-S*(t)∥=M1∥S(t)-S*(t)∥,(6)

where M1=(Πh+ΠB+μν).

Thus,

∥F1(t,S(t))-F1(t,S*(t))∥≤M1∥S(t)-S*(t)∥,(7)

Repeating the same procedure above, yields:

∥F2(t,E(t))-F2(t,E*(t))∥≤M2∥E(t)-E*(t)∥,∥F3(t,Q(t))-F3(t,Q*(t))∥≤M3∥Q(t)-Q*(t)∥,∥F4(t,IA(t))-F4(t,IA*(t))∥≤M4∥IA(t)-IA*(t)∥.(8)

Subsequently, Eq. (4) gives:

Sn(t)=S(0)+1-νN(ν)F1(t,Sn-1(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F1(x,Sn-1(x))dx,In(t)=E(0)+1-νN(ν)F2(t,In-1(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F2(x,In-1(x))dx,Rn(t)=Q(0)+1-νN(ν)F3(t,Rn-1(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F3(x,Rn-1(x))dx,Bn(t)=B(0)+1-νN(ν)F4(t,Bn-1(t))+νN(ν)1Γ(ν)∫0t(t-x)ν-1F4(x,Bn-1(x))dx,(9)

where S(t)≥0, I(t)≥0, R(t)≥0, B(t)≥0. The difference between successive components can be denoted by Φni,i=1,2,⋯6, respectively. Thus in view of Eqs. (7)–(9), we obtain:

∥Φn1(t)∥=1-νN(ν)M1∥S(t)-S*(t)∥+νN(ν)M1Γ(ν)∫0t(t-x)ν-1∥S(t)-S*(t)∥dx,∥Φn2(t)∥=1-νN(ν)M2∥I(t)-I*(t)∥+νN(ν)M2Γ(ν)∫0t(t-x)ν-1∥I(t)-I*(t)∥dx,∥Φn3(t)∥=1-νN(ν)M3∥R(t)-R*(t)∥+νN(ν)M3Γ(ν)∫0t(t-x)ν-1∥R(t)-R*(t)∥dx,∥Φn4(t)∥=1-νN(ν)M4∥B(t)-B*(t)∥+νN(ν)M4Γ(ν)∫0t(t-x)ν-1∥B(t)-B*(t)∥dx,(10)

Theorem 3.1 The fractional proposed model possesses a unique solution for t∈[0,T] if the condition is satisfied

(1-νN(ν)Mi+1N(ν)Miγ(ν)Tν)<1,i=1,2,…,6.(11)

Proof. Based on the assumptions that S(t), E(t), Q(t), IA(t), IS(t), R(t) are bounded functions, it is therefore clear that the kernels F1, F2, F3, F4, F5, F6 from Eqs. (7)–(8) satisfy the Lipschitz condition. Hence, Eq. (10) can be viewed as:

∥Φn1(t)∥≤(1-νN(ν)M1+1N(ν)M1Γ(ν)Tν)n,∥Φn2(t)∥≤(1-νN(ν)M2+1N(ν)M2Γ(ν)Tν)n,∥Φn3(t)∥≤(1-νN(ν)M3+1N(ν)M3Γ(ν)Tν)n,∥Φn4(t)∥≤(1-νN(ν)M4+1N(ν)M4Γ(ν)Tν)n.(12)

Hence, the sequences above exist as n→∞,∥Φni(t)∥→0,i=1,2,⋯6. Also, using the triangular inequality for any k value, Eq. (12) yields:

∥Sn+k(t)-Sn(t)∥≤∑i=n+1n+kP1i=P1n+1-P1n+k+11-P1,∥In+k(t)-In(t)∥≤∑i=n+1n+kP2i=P2n+1-P2n+k+11-P2,∥Rn+k(t)-Rn(t)∥≤∑i=n+1n+kP3i=P3n+1-P3n+k+11-P3,∥Bn+k(t)-Bn(t)∥≤∑i=n+1n+kP4i=P4n+1-P4n+k+11-P4,(13)

4  Stability and Basic Reproductive Number

In this section, the stability of the ailment-free equilibrium and its analytic conditions will be discussed. From the derivation of the existence of equilibria in [23], the basic reproduction number in the fractional model is given by:

R0=βhνμν+γν+dν+(Πβeε)νμνσν(μν+γν+dν)C.(14)

Theorem 4.1 The ailment-free equilibrium E0 is locally asymptotically stable if R0≤1 and unstable if R0>1.

Proof. Through the concept of the Jacobian matrix, local stability at E0 can be achieved by:

J[E0]=[-μν-βhνφν-βBΠKμ0βhν-(μν+dν+γν)0βBΠKμ0γν-(φν+μ)00ε0-σν.](15)

The associated eigenvalues are Π1=-μν,Π2=-(μν+φν) and the solutions of Π2+((μν+dν+γν+σν)-βh)Π+(μν+dν+γν)σν(1-R0)=0. Now, if R0<1, ((μν+dν+γν+σν)-βh)<0 and ((μν+dν+γν+σν)-βh)(μν+dν+γν)σν(1-R0)>0, in accordance with the Hurwitz principle, the expression Π2+((μν+dν+γν+σν)-βh)Π+(μν+dν+γν)σν(1-R0)=0, possesses negative real eigenvalues. Thus, E0 is locally asymptotically stable. On the other hand, when R0=1, one eigenvalue of the expression will be 0. By the concept of Hartman Grobman [24], the ailment-free equilibrium is nonhyperbolic. The absence of non hyperbolicity indicates that the linearized system cannot describe the local features of the system.

4.1 Global Stability of the Ailment-Free Equilibrium

By employing the approach used in [25], the asymptotic stability for the ailment-free equilibrium will be derived in a global sense. This takes the following into account:

ABCDνA′(T)=f(A,y),ABCDνy′(T)=g(A,y),g(A,0)=0,(16)

where A∈ℜm indicates the non-infectious ones and y∈ℜn indicates the infectious ones. Assume U0=(A*,0) represents the ailment-free equilibrium and suppose that:

•   ABCDνA′(T)=f(A,0), A* is globally asymptotically stable,

•   g(A,y)=Xy-g(A,y),g(A,y)≥0 for (A,y)∈X, and X=D1G(A*,0) denotes an M-matrix.

Theorem 4.2 If D1 and D2 are satisfied, and U0 = (A*, 0) is a fixed point, then the model is globally asymptotically stable only if Ro≤1.

Proof. The proof can be shown in a similar way as the process carried out in [24].

5  Parameters Estimation

For the validation of an epidemiological model, it is extremely important to compare the results of simulations with the actual data of infected individuals. This increases the reliability of the proposed disease model. Similar values from the simulations and actual data give better information on the disease being investigated. In addition, unknown values of the working parameters that contribute to the model can be determined.

There are different techniques including maximum likelihood estimation, Bayesian technique, nonlinear least-squares approach, and probability plotting, which can be used to obtain the best parameters. In this research, we utilized the nonlinear least-squares approach for computing the best-fitted parameters, including Π,βh,βe,d,ν,ε, and σ, along with the most important parameter of our proposed model, called the fractional order ν (one of the major components of the study). We also obtained the best-fitted parameters for the classical dysentery system and ABC fractional dysentery system, which are shown in Tab. 1. The actual data for the dysentery diarrhea infected individuals was from Ethiopia, which covered a period of 52 weeks in 2017 [26].

When the least-squares technique is utilized, we need to minimize the objective function. This is achieved by tuning the system’s parameters to fit the available data points accurately. Real data cases for the dysentery diarrhea disease in this study are denoted by m points (xn,yn),n=1,…,m, where xn stands for independent quantity and yn shows dependent quantity. The system function has the structure h(x, r), where s tuned parameters are shown in a vector of parameters, r. Hence, the objective is to identify the parameters which ensure that the system’s simulations for the infected cases fit well with the actual data points.

In the present study, we obtained the best fit by measuring the difference between the real data and the simulations. This is shown below:

En=yn-h(xn,r)(17)

Finally, the optimal set of parameters is obtained, as shown in Tab. 1, with the least-squares approach while minimizing the absolute relative error, on average, as shown below:

ARE=1N∑n=1m|yn-h(xn,r)|yn,(18)

where N stands for the total data value, which is 52 in this study. Moreover, real prevalent cases of the dysentery diarrhea disease, along with the classical and the ABC system’s simulations for the infected individuals, are listed in Tab. 2. In addition, Fig. 1 shows the best fit of the classical and ABC system with the real cases. The ABC system showed an average absolute relative error of 3.4284e −02, and the classical system was 3.4432e −02. Therefore, it shows that the ABC system had some advantages compared to the classical dysentery diarrhea system. In addition, the basic reproductive numbers R0 were 1.7031 and 1.9581 for the ABC and classical systems, respectively. This clearly shows that the disease can be well prevented if the ABC operator is taken into consideration while modeling the epidemic.

Figure 1: Plots and comparison of the real incidence cases, classical (ν=1) and ABC (ν=9.9410e-01) systems for the dysentery diarrhea

6  Sensitivity

In this section, the concept of sensitivity analysis is used to discover the robust significance of the generic parameters present in the base reproduction number R0. Furthermore, both the analytical and numerical values of the R0 parameters are derived from precise assumptions using parameter values.

If the dynamics follow the model (1), the analytical expressions can be used to explain the process of tracking the model’s onset at various locations. The threshold value, R0 can reduce and stop the ailment spread by reducing the number to less than unity. The sensitivity index technique is used to measure the most sensitive parameters in the model. The parameters with a positive sign are considered highly and proportionally sensitive to R0, while those with a negative sign are less sensitive to R0 decreasing. The other category is neutrally sensitive (with zero relative sensitivity). The cause of the transmission of the infringement is directly linked to the specific reproduction number R0. The R0 elasticity indices [27] are shown below:

ΥPiR0=∂R0∂Pi×PiR0,(19)

where R0 denotes the basic reproduction ratio and Pi is stated above. Following the described formula, we get:

Υβh=βhνν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν+γν+dν)-1,Υμ=μ(-βhνμνν(μν+γν+dν)2μ-(πβeε)ννμνσν(μν+γν+dν)Cμ-(πβeε)ννσν(μν+γν+dν)2Cμ)×(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1,Υγ=γ(-βhνγνν(μν+γν+dν)2γ-(πβeε)νγννμνσν(μν+γν+dν)2Cγ)×(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1,Υd=d(-βhνdνν(μν+γν+dν)2d-(πβeε)νdννμνσν(μν+γν+dν)2Cd)×(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1,Υσ=-(πβeε)νν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν)-1(σν)-1(μν+γν+dν)-1C-1,ΥC=-(πβeε)ν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν)-1(σν)-1×(μν+γν+dν)-1C-1,Υπ=(πβeε)νν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν)-1(σν)-1×(μν+γν+dν)-1C-1,Υε=(πβeε)νν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν)-1(σν)-1×(μν+γν+dν)-1C-1,Υβe=(πβeε)νν(βhνμν+γν+dν+(πβeε)νμνσν(μν+γν+dν)C)-1(μν)-1×(σν)-1(μν+γν+dν)-1C-1.(20)

The numerical values indicating the relative significance of R0 are given in Tab. 3. Some parameters are positive while some are negative. Parameters with positive values mean that an increase in the parameter’s values will have a major effect on the frequency of the ailment spread. On the other hand, parameters with negative values mean an increase in such parameters would decrease the effect of the disease. A representation of the values given in Tab. 3, is shown in Fig. 2.

Figure 2: Elasticity indices for various parameters of R0

7  Simulations for the ABC Model

In this section, an algorithm is first developed to obtain the approximate solution of the ABC dysentery diarrhea model, wherein the operator uses non-local and non-singular types of the kernel. The algorithm being developed is discussed in [28], which entails the combination of the fundamental theorem of fractional calculus and two-step Lagrange type polynomial. Therefore, the fundamental theorem of fractional calculus on the Cauchy type initial value problem is given below:

ABCD0νw(t)=G(t,w(t)),w(0)=w0,0<t<∞,(21)

This leads to:

w(t)-w(0)=(1-ν)ABC(ν)G(t,w(t))+νABC(ν)×Γ(ν)∫0tG(τ,w(τ))(t-τ)ν-1dτ.(22)

At t=tn+1,n=0,1,2,…, we have:

w(tn+1)-w(0)=(1-ν)ABC(ν)G(tn,w(tn))+νABC(ν)×Γ(ν)∫0tn+1G(τ,w(τ))(tn+1-τ)ν-1dτ. (23)

w(tn+1)-w(0)=(1-ν)ABC(ν)G(tn,w(tn))+νABC(ν)×Γ(ν)∑v=0n∫tvtv+1G(τ,w(τ))(tn+1-τ)ν-1dτ. (24)

With the help of interpolation polynomial, we approximate function G(τ,w(τ))over[tv,tv+1]:

G(τ,w(τ))≈Pk(τ)=G(tv,w(tv))h(τ-tv-1)-G(tv-1,w(tv-1))h(τ-tv).(25)

Eq. (24) becomes:

w(tn+1)=w(0)+(1-ν)ABC(ν)G(tn,w(tn))+νABC(ν)×Γ(ν)∑v=0n(G(tv,w(tv))h∫tvtv+1(τ-tv-1)(tn+1-τ)ν-1dτ-G(tv-1,w(tv-1))h∫tvtv+1(τ-tv)(tn+1-τ)ν-1dτ).(26)

By solving the above integrals, we obtain the approximate solution shown below:

wn+1=w(t0)+(1-ν)ABC(ν)G(tn,w(tn))+νABC(ν)∑v=0n (27)

[hνG(tv,w(tv))Γ(ν+2)((n+1-v)ν(n-v+2+ν)-(n-v)ν(n-v+2+2ν)) (28)

hνG(tv-1,w(tv-1))Γ(ν+2)((n+1-v)ν+1-(n-v)ν(n-v+1+ν))]. (29)

Hence, the proposed dysentery diarrhea model becomes:

Sn+1=S(t0)+(1-ν)ABC(ν)G1(tn,w(tn))+νABC(ν)∑v=0n (30)

[hνG1(tv,w(tv))Γ(ν+2)((n+1-v)ν(n-v+2+ν)-(n-v)ν(n-v+2+2ν)) (31)

hνG1(tv-1,w(tv-1))Γ(ν+2)((n+1-v)ν+1-(n-v)ν(n-v+1+ν))], (32)

In+1=I(t0)+(1-ν)ABC(ν)G2(tn,w(tn))+νABC(ν)∑v=0n (33)

[hνG2(tv,w(tv))Γ(ν+2)((n+1-v)ν(n-v+2+ν)-(n-v)ν(n-v+2+2ν)) (34)

hνG2(tv-1,w(tv-1))Γ(ν+2)((n+1-v)ν+1-(n-v)ν(n-v+1+ν))], (35)

Rn+1=R(t0)+(1-ν)ABC(ν)G3(tn,w(tn))+νABC(ν)∑v=0n (36)

[hνG3(tv,w(tv))Γ(ν+2)((n+1-v)ν(n-v+2+ν)-(n-v)ν(n-v+2+2ν)) (37)

hνG3(tv-1,w(tv-1))Γ(ν+2)((n+1-v)ν+1-(n-v)ν(n-v+1+ν))], (38)

Bn+1=B(t0)+(1-ν)ABC(ν)G4(tn,w(tn))+νABC(ν)∑v=0n (39)

[hνG4(tv,w(tv))Γ(ν+2)((n+1-v)ν(n-v+2+ν)-(n-v)ν(n-v+2+2ν)) (40)

hνG4(tv-1,w(tv-1))Γ(ν+2)((n+1-v)ν+1-(n-v)ν(n-v+1+ν))], (41)

where,

G1=Πν+ανR-(Πh+ΠB+μν)S,G2=(Πh+ΠB)S-(μν+γν+dν)I,G3=γνI-(μν+αν)R,G4=ενI-σνB.(42)

Using numerical simulations, the ABC dysentery diarrhea model (2) uses the developed algorithm shown above. Different values of the key parameters are chosen to investigate their effects on the dynamics of the disease. These parameter values are taken from Tab. 1. With simulations, we can identify important parameters that significantly affect the model’s dynamics. In order to investigate the transmission rate of dysentery due to human to human interaction (βh) and environment to human interaction (βe), we vary both parameters, as shown in Fig. 3.

Figure 3: Behavior of diarrhea infected people with increasing values of (a) effective transmission rate of dysentery due to human to human interaction (βh) and (b) effective transmission rate of dysentery due to the environment to human interaction (βe), while taking the remaining parameters from Tab. 1

This shows that diarrhea disease is principally due to the environment to human interaction. It means that humans must take care of their hygiene and surroundings in order to avoid the spread of dysentery. In order to investigate the effect of the concentration of Shigella pathogen (C) and the net death rate of Shigella pathogen (σ), we have varied both parameters, as shown in Fig. 4. This demonstrates that decreasing the values of the net death of Shigella pathogen leads to a reduction in the infection, while its concentration (C) follows the same pattern as the normal curve.

Figure 4: Behavior of diarrhea infected people with increasing values of (a) concentration of Shigella pathogen (C), and (b) net death rate of Shigella pathogen (σ), while taking the remaining parameters from Tab. 1

In order to investigate the effect of the pathogen shedding rate of infected humans (ε) and the recovery rate of dysentery (γ), we have varied both parameters, as shown in Fig. 5. This demonstrates that increasing the values of ε causes the infection rate to rise. This is because its shedding rate on humans also increases, whereas a slight improvement in the recovery rate decreases the infection rate substantially.

Figure 5: Behavior of diarrhea infected people with increasing values of (a) pathogen shedding rate of infected humans (ε) and (b) recovery rate of dysentery (γ), while taking the remaining parameters from Tab. 1

Figure 6: Behavior of basic reproductive number, R0 under the influence of different parameters of the proposed ABC system (2)

Finally, in order to investigate the effects of different parameters on the basic reproductive number, R0, we have presented various contour plots, as shown in Fig. 6. This shows that dysentery diarrhea can be prevented effectively if some strategies, including environment to human interaction, the concentration of Shigella pathogen, and pathogen shedding rate, are carefully considered.

8  Conclusion

In the present research, one of the robust non-local and non-singular fractional operator, called Atangana-Baleanu, was used to model dysentery diarrhea. The employed fractional operator was suitable for the investigation of transmission dynamics of a disease from the literature. The fractionalized order is ν, and the dimensional consistency between the rest of the parameters has been considered. Consequently, several important features of the proposed fractional variant of the model, such as the formation of the model, existence and uniqueness of the solution by means of fixed point theorem, stability analysis, sensitivity analysis and most importantly, the estimation of the parameters have been reported. It was observed that the fractional variant of the model represents the behavior of the ailment more accurately than the integer-order variant. The dimensional homogeneity among parameters and the ABC operator have been authenticated by fitting the dysentery diarrhea infected cases from Ethiopia covering a period of 52 weeks in 2017. In addition, the ABC operator’s fractional order, ν was optimized. This returns a basic reproduction number, R0 of 1.7031 and R0 of 1.9581, for the integer model. It is shown that the ABC model has smaller values than the integer-order, whereas the proposed models match the actual cases of infected humans.

Furthermore, in order to shed more light on the features of the model, various numerical simulations were carried out using an effective numerical scheme. In future studies, we plan to apply the techniques used in [29–34] to understand dysentery diarrhea dynamics in greater detail. In addition, optimal control theory will be utilized to devise effective control strategies to eliminate the epidemic.

Acknowledgement: This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under Grant No. RG-6-135-40. The authors, therefore, gratefully acknowledge DSR technical and financial support.

Funding Statement: This research is supported by the King Abdulaziz University RG-6-135-40.

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

References

 1.  Rachah, A., Torres, D. F. M. (2016). Dynamics and optimal control of Ebola transmission. Mathematics in Computer Science, 10(3), 331–342. DOI 10.1007/s11786-016-0268-y.

 2.  Ndaïrou, F., Area, I., Nieto, J. J., Silva, C. J., Torres, D. F. M. (2018). Mathematical modeling of Zika disease in pregnant women and newborns with microcephaly in Brazil. Mathematical Methods in the Applied Sciences, 41(18), 8929–8941. DOI 10.1002/mma.4702.

 3.  Djordjevic, J., Silva, C. J., Torres, D. F. M. (2018). A stochastic SICA epidemic model for HIV transmission. Applied Mathematics Letters, 84(1), 168–175. DOI 10.1016/j.aml.2018.05.005.

 4.  Brauer, F., Castillo-Chavez, C., Feng, Z. (2019). Mathematical models in epidemiology. USA: Springer.

 5.  Rao, F., Mandal, P. S., Kang, Y. (2019). Complicated endemics of an SIRS model with a generalized incidence under preventive vaccination and treatment controls. Applied Mathematical Modelling, 67, 38–61. DOI 10.1016/j.apm.2018.10.016.

 6.  Li, L., Sun, C., Jia, J. (2018). Optimal control of a delayed SIRC epidemic model with saturated incidence rate. Optimal Control Applications and Methods, 40(2), 367–374. DOI 10.1002/oca.2482.

 7.  Ullah, S., Khan, M. A., Gómez-Aguilar, J. F. (2019). Mathematical formulation of hepatitis B virus with optimal control analysis. Optimal Control Applications and Methods, 40(3), 529–544. DOI 10.1002/oca.2493.

 8.  Lin, J., Xu, R., Tian, X. (2018). Global dynamics of an age-structured cholera model with both human-to-human and environment-to-human transmissions and saturation incidence. Applied Mathematical Modelling, 63, 688–708. DOI 10.1016/j.apm.2018.07.013.

 9.  Bonyah, E., Khan, M. A., Okosun, K. O., Gómez-Aguilar, J. F. (2019). On the co-infection of dengue fever and Zika virus. Optimal Control Applications and Methods, 40, 394–421. DOI 10.1002/oca.2483.

10. Ghosh, I., Tiwari, P. K., Chattopadhyay, J. (2019). Effect of active case finding on dengue control: Implications from a mathematical model. Journal of Theoretical Biology, 464(7), 50–62. DOI 10.1016/j.jtbi.2018.12.027.

11. Eikenberry, S. E., Gumel, A. B. (2018). Mathematical modeling of climate change and malaria transmission dynamics: A historical review. Journal of Mathematical Biology, 77(4), 857–933. DOI 10.1007/s00285-018-1229-7.

12. Kumar, D., Singh, J., Tanwar, K., Baleanu, D. (2019). A new fractional exothermic reactions model having constant heat source in porous media with power, exponential and Mittag–Leffler laws. International Journal of Heat and Mass Transfer, 138(8), 1222–1227. DOI 10.1016/j.ijheatmasstransfer.2019.04.094.

13. Singh, J., Kumar, D., Baleanu, D. (2018). On the analysis of fractional diabetes model with exponential law. Advances in Difference Equations, 2018(1), 1–15. DOI 10.1186/s13662-018-1680-1.

14. Kumar, D., Singh, J., Baleanu, D. (2018). A new analysis of the Fornberg–Whitham equation pertaining to a fractional derivative with Mittag–Leffler-type kernel. The European Physical Journal Plus, 133(2), 2–10. DOI 10.1140/epjp/i2018-11833-3.

15. Atangana, A., Baleanu, D. (2016). New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Thermal Science, 20(2), 763–769. DOI 10.2298/TSCI160111018A.

16. Pinto, C. M. A., Carvalho, A. R. M. (2017). A latency fractional order model for HIV dynamics. Journal of Computational and Applied Mathematics, 312(1), 240–256. DOI 10.1016/j.cam.2016.05.019.

17. Jajarmi, A., Baleanu, D. (2018). A new fractional analysis on the interaction of HIV with CD4 + T-cells. Chaos, Solitons Fractals, 113(1), 221–229. DOI 10.1016/j.chaos.2018.06.009.

18. Rosa, S., Torres, D. F. M. (2018). Optimal control of a fractional order epidemic model with application to human respiratory syncytial virus infection. Chaos, Solitons Fractals, 117(7–8), 142–149. DOI 10.1016/j.chaos.2018.10.021.

19. Caputo, M., Fabrizio, M. (2015). A new definition of fractional derivative without singular kernel. Progress in Fractional Differentiation and Applications, 2(1), 73–85. DOI 10.12785/pfda/010201.

20. Yusuf, A., Qureshi, S., Inc, M., Aliyu, A. I., Baleanu, D. et al. (2018). Two-strain epidemic model involving fractional derivative with Mittag-Leffler kernel. Chaos, 28(12), 123121. DOI 10.1063/1.5074084.

21. Qureshi, S., Yusuf, A., Shaikh, A. A., Inc, M., Baleanu, D. (2019). Fractional modeling of blood ethanol concentration system with real data application. Chaos, 29, 13143. DOI 10.1063/1.5082907.

22. Qureshi, S., Yusuf, A., Ali Shaikh, A., Inc, M., Baleanu, D. (2020). Mathematical modeling for adsorption process of dye removal nonlinear equation using power law and exponentially decaying kernels. Chaos: An Interdisciplinary Journal of Nonlinear Science, 30(4), 43106. DOI 10.1063/1.5121845.

23. Berhe, H. W., Makinde, O. D., Theuri, D. M. (2019). Parameter estimation and sensitivity analysis of dysentery diarrhea epidemic model. Journal of Applied Mathematics, 2019(3), 1–13. DOI 10.1155/2019/8465747.

24. Castillo-Chavez, C., Song, B. (2004). Dynamical models of tuberculosis and their applications. Mathematical Biosciences and Engineering, 1(2), 361–404. DOI 10.3934/mbe.2004.1.361.

25. Berhe, H. W., Makinde, O. D., Theuri, D. M. (2019). Co-dynamics of measles and dysentery diarrhea diseases with optimal control and cost-effectiveness analysis. Applied Mathematics and Computation, 347(7), 903–921. DOI 10.1016/j.amc.2018.11.049.

26. Berhe, H. W., Qureshi, S., Shaikh, A. A. (2020). Deterministic modeling of dysentery diarrhea epidemic under fractional Caputo differential operator via real statistical analysis. Chaos, Solitons & Fractals, 131(1), 1–13. DOI 10.1016/j.chaos.2019.109536.

27. Asamoah, J. K. K., Jin, Z., Sun, G. Q., Li, M. Y. (2020). A Deterministic model for Q fever transmission dynamics within dairy cattle herds: Using sensitivity analysis and optimal controls. Computational and Mathematical Methods in Medicine, 2020(3), 1–18. DOI 10.1155/2020/6820608.

28. Toufik, M., Atangana, A. (2017). New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models. The European Physical Journal Plus, 132(10), 2–16. DOI 10.1140/epjp/i2017-11717-0.

29. Triet, N. A., Van Au, V., Long, L. D., Baleanu, D., Tuan, N. H. (2020). Regularization of a terminal value problem for time fractional diffusion equation. Mathematical Methods in the Applied Sciences, 43(6), 3850–3878. DOI 10.1002/mma.6159.

30. Kumar, P., Qureshi, S. (2020). Laplace–Carson integral transform for exact solutions of non-integer order initial value problems with Caputo operator. Journal of Applied Mathematics and Computational Mechanics, 19(1), 57–66. DOI 10.17512/jamcm.2020.1.05.

31. Tuan, N. H., Huynh, L. N., Baleanu, D., Can, N. H. (2020). On a terminal value problem for a generalization of the fractional diffusion equation with hyper-Bessel operator. Mathematical Methods in the Applied Sciences, 43(6), 2858–2882. DOI 10.1002/mma.6087.

32. Qureshi, S., Chang, M., Shaikh, A. A. (2020). Analysis of series RL and RC circuits with time-invariant source using truncated M, atangana beta and conformable derivatives. Journal of Ocean Engineering and Science, 18(5), 507. DOI 10.1016/j.joes.2020.11.006.

33. Bao, N. T., Baleanu, D., Minh, D., Huy, T. N. (2020). Regularity results for fractional diffusion equations involving fractional derivative with Mittag–Leffler kernel. Mathematical Methods in the Applied Sciences, 43(12), 7208–7226. DOI 10.1002/mma.6459.

34. Tuan, N. H., Tri, V. V., Baleanu, D. (2020). Analysis of the fractional corona virus pandemic via deterministic modeling. Mathematical Methods in the Applied Sciences, 44(1), 1086–1102. DOI 10.1002/mma.6814.

images
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.