Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Analysis of Stochastic COVID-19 and Hepatitis B Co-infection Model with Brownian and Lévy Noise

View ORCID ProfileMichael A. Pobbi, View ORCID ProfileStephen E. Moore
doi: https://doi.org/10.1101/2024.08.12.24311861
Michael A. Pobbi
1Department of Mathematics, University of Cape Coast, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael A. Pobbi
  • For correspondence: michael.pobbi{at}stu.ucc.edu.gh
Stephen E. Moore
1Department of Mathematics, University of Cape Coast, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stephen E. Moore
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

In this article, we formulate and analyze a mathematical model for the coinfection of HBV and COVID-19 that incorporates the effects of Brownian and Lévi noise. We studied the dynamics and effects of these diseases in a given population. First, we establish the basic reproduction number of the disease-free equilibrium point of the stochastic model by means of a suitable Lyapunov function. Additionally, we provided sufficient conditions for the stability of the model around the disease-free equilibrium points. Finally, using a few simulation studies, we demonstrate our theoretical results. In particularly, we derived threshold values for HBV only Embedded Image, COVID-19 only, Embedded Image, and coinfection Embedded Imagefor the stochastic model around disease-free equilibrium point. Next, the conditions for stability in the stochastic sense for HBV only, COVID-19 only submodels, and the full model are established. Furthermore, we devote our concentrated attention to sufficient conditions for extinction and persistence using each of these reproductive numbers. Finally, by using the Euler–Murayama scheme, we demonstrate the dynamics of the coinfection by means of numerical simulations.

1 Introduction

Both Coronavirus Disease 2019 (COVID-19) and Hepatitis B virus (HBV) infections remain major public health threats [4]. COVID-19 is caused by Severe Acute Respiratory Syndrome Corona virus 2 (SARS-CoV-2) and affect the functionality of the respiratory system. It is spread from one person to another by coming into direct contact with respiratory droplets from an infected individual. It can also spread by touching surfaces contaminated with the virus and then touching the face. On the other hand, HBV disease is caused by the hepatitis B virus and is associated with serious liver infection, leading to liver malfunction, including cirrhosis, liver fibrosis, and hepatocellular carcinoma. Both COVID-19 and HBV are life-threatening diseases that have been declared global health problems [18]. For example, [32] reported that more than 3 billion people are infected with HBV globally. An estimated 7 million COVID-19 deaths and 700 million infections have been reported since December 2019 [18]. These statistics demonstrate how devastating HBV and COVID-19 diseases are worldwide. Presently, global efforts have been made to develop and distribute COVID-19 vaccines, aiming to reduce severe illness, hospitalizations, and deaths.

Over the years, mathematical models have provided frameworks for explaining transmission dynamics and for making necessary interventions for several infectious diseases; see, e.g. [2, 11, 28, 5, 7, 16]. Several models of COVID-19, HBV and other infectious disease transmission dynamics have been developed in the literature.

The transmission dynamics of COVID-19, HBV, and other diseases have been extensively studied, mostly in isolation. Both deterministic and stochastic models have been developed to understand the spread of these diseases [13, 11, 9]. In [11], the authors deemed the susceptible-exposed-infectious-recovered (SEIR) model introduced by [34] to describe the transmission dynamics and forecast the national and global spread of COVID-19 disease to be simplistic and largely unrealistic because it excluded the relapse from recovered class to susceptible class. Like [34], authors including [22, 30] also made an assumption of no relapse from the recovered to the susceptible class, an assumption that is critiqued in [28]. [26] examine the extinction and persistence for a HBV only model and also conducted simulations on the model.

Amidst inconclusive findings from clinical studies on whether one disease increases the severity of another disease in the event of coinfection (see e.g. in [4, 29]), there has been a spike in interest for both clinical and theoretical scientific research on disease coinfections in recent times. The theoretical study of coinfection through mathematical modelling will be valuable to shaping our understanding of transmission dynamics of these diseases. A number of studies that have been done on co-dynamics of COVID-19 including Dengue [23], cholera [10], malaria, [20, 1, 21, 3, 19] and co-dynamics s TB, see, e.g. [24, 14]. [29] employed a deterministic model and examined the effects of optimal control strategies on the Covid-19, HBV coinfection transmission. On the other hand, [8] examined COVID-19 and HBV coinfection by employing Gaussian noise to explain the perturbation in observed data for the COVID-19 and Malaria coinfection model. Stochastic coinfection models are very scarce in literature, and more so for COVID-19 and HBV Stochastic coinfection models. While mathematical models, see, e.g., [17, 28], have been proposed to describe the dynamics of diseases including COVID-19 and HBV, there is a dearth in studies imprint the effect stochastic perturbations on COVID-19-HBV coinfection. In this paper, we employ both Gaussian noise and Non-Gaussian noise to explain the noise in the Coinfection of COVID -19 and HBV in the presence of COVID-19 Vaccination, and determine conditions for extinction and persistence for the coinfection model.

Stochastic coinfection SIRS models unlike classic or deterministic SIRS models, incorporate the effects of random fluctuations in the transmission dynamics of disease and recovery rates, as well as the possibility of coinfection by different diseases. In [35, 15], the authors formulated techniques for examining existence, uniqueness and persistence of solution. Understanding and analysing the dynamics of infectious illness transmission in both humans and animals can benefit from mathematical modelling. The persistence or extinction of the diseases depends on several factors, such as the basic reproduction numbers, the coinfection rates, and the stochastic perturbations; hence, determining the conditions for the persistence or extinction of the diseases is a cumbersome task. One of the methods proposed to analyze the persistence or extinction of stochastic coinfection SIRS models is the use of Lyapunov functions, which are scalar functions that measure the distance from an equilibrium point [33]. For example, [33] employed Lyapunov functions to show that if the basic reproduction numbers of both diseases are less than one, then the disease-free equilibrium is globally asymptotically stable for a stochastic coinfection SIRS model with nonlinear incidence rates and cross-immunity effects. In this paper we develop condition for extinction and persistence for coinfection model.

In this article, we will formulate a mathematical model which builds on (2) of [8] by firstly assuming relapse to susceptible state. The new model includes COVID-19 vaccination state, and also assume the presence of both standard random fluctuation and massive disturbances which arise due to sudden environmental shocks in the transmission dynamics. The novel coinfection model will open up some new avenues for further research in stochastic modelling dynamics of coinfection of diseases. Specifically, our study seeks to investigate the following objectives:

  1. To develop a compartmental model for explaining the transmission dynamics of HBV-COVID-19 coinfection.

  2. To establish the Invariance and Positivity of the stochastic system for coinfection Model.

  3. To determine the basic Reproductive number for Stochastic HBV Only, COVID-19 Only and Coinfection models respectively.

  4. To determine conditions for stability at Disease Free Equilibrium for Stochastic HBV Only, COVID-19 Only and coinfection models respectively using Brownian Noise and Poisson Lévi Noise.

  5. To present numerical simulations and discussions of the model.

The rest of the paper is organised as follows; Section 2 discusses the Mathematical Model and the underlying assumptions. In Section 3, we analyze the COVID-19, and HBV sub-models as well as the Coinfection model. Section 4 is devoted to discussing the numerical method and numerical simulation of the model. Finally, we present the conclusion in Section 5.

2 Mathematical Model Formulation

In this section, we consider a stochastic compartmental mathematical model that uses differential equations, see Figure 1.

Figure 1:
  • Download figure
  • Open in new tab
Figure 1:

Compartmental diagram for the COVID-19 Hepatitis B coinfection cohort.

The sub-populations of the human population at time t, denoted by N (t), include susceptible individuals S(t), COVID-19 vaccinated individuals V (t), COVID-19 infected individuals IC(t), Hepatitis-B infected individuals IH (t), individuals infected with both COVID-19 and Hepatitis-B IHC(t), and recovered individuals R(t).

We will consider the following assumptions;

(A1) The total human population (N (t)) is a sum of all the different compartments, i.e. Embedded Image

(A2) Covid-19 Vaccination wears off after some period of time. (A3) There is relapse to susceptible state,

(A4) There is no vertical transmission. Embedded Image with initial conditions S(0) > 0, IC(0) ≥0, IH (0) ≥ 0, IHC(0) ≥0, R(0) ≥0, V (0) ≥ 0. βC, βH are the effective contact rate for COVID-19, and HBV respectively.

Two popular ways to introduce stochastic factors into epidemic models are: (i) to assume some small and standard random fluctuation, (ii) to assume massive disturbances caused by sudden environmental shocks. The first, is described by Gaussian white noise while the second is described by the Lévy noise. It is assumed that fluctuations in the environment will manifest mainly as fluctuations in the parameter β, i.e. β → β + ϵidBi(t), (i = 1, 2, …, 6), where dBi(t) is a one dimensional standard Brownian motion with Bi(0) = 0, and ϵi is the intensity of the white noise. Also a multiplicative noise is considered in our model, as random term depends on the state space, e.g. [8]. The Lévy jump process, thus, is divided into linear drift term, Brownian motion and compensated Poisson Process. The stochastic model of the system (2) takes the following form: Embedded Image where dBi(t) is independent of Embedded Image. 𝒩 is a Poisson counting measure with compensating martingale Embedded Image and characteristic measure v on a measurable subset 𝒰 ∈ (0, ∞) where v(𝒰) ≤ ∞ and Ψi : 𝒰 → ℜ (i = 1, 2, 3, 4, 5, 6) denotes the effects of random jumps and is bounded and continuous. It is assumed that v is a Levy measure so that Embedded Image. We note that S(t−), IH (t−), IC(t−), IHC(t−), R(t−), V (t−) denote the left limits of S(t), IH (t), IC(t), IHC(t), R(t), V (t) respectively.

The following definitions are necessary for the analysis of the stochastic model.

Definition 1

(Itô-Lévi [26]). First, set Embedded Image. Suppose a complete probability space (Ω, F, Ft ≥ 0, P) with filtration Ft ≥ 0, which satisfies the usual conditions. Bi(t) is defined on this probability space.

We assume an Itô -Lévi process, X(t) ∈ ℜ+, of the form Embedded Image where, F : ℜn × ℜ+ × S → ℜn, G : ℜn × ℜ+ × S → ℜn, and Ψ : ℜn × ℜ+ × S × Y → ℜn are measurable functions and X(t−) denotes the left limit of X(t). F (X(t), t) represents the linear drift term, G(X(t), t) the Brownian noise, Embedded Image (dt, du) represents the compensated Poisson random measure, and Ψ(X, t, u) represents the intensity of jumps. It is assumed that the conditionEmbedded Image

Definition 2

(Itô-Lévi Formula [26]). We consider the process X expressed by (4) and let 𝒱 ∈ C2,1(ℜn × ℜ+ × S; ℜ+) such that 𝒴 ≡ 𝒱(t, X(t)). Then, 𝒴(t) is again an Itô-Lévi process and Embedded Image

Which is the representation of generalized Itô’s formula with jumps. Embedded Image

To gain insights into the underlying mechanisms that drive the spread of infectious diseases to identify the long-term behaviour of disease transmission we conduct qualitative analysis of the stochastic model (3).

3 Qualitative Analysis of Model

In this section, we present the Positivity of the Solution to (3), determine the Reproductive Number for Single Disease Models and Coinfection model, determine the conditions for Local and Global equilibrium for COVID-19 only, HBV only, and for Coinfection model. Further we determine the conditions for Extinction and Persistence for all three models.

3.1 Positivity of Model

A first step is to determine the positivity of solution to the stochastic system (3) as a negative state will not be biologically meaningful. We will show the positivity of the first equation in(3).

Theorem 1

Let ((S(t), IH (t), …, V (t)) be a solution of system (3), given initial values ((S(0), IH (0), …, V (0)) ∈ Ω, where Ω ∈ ℜ+, then limt→∞ S(t) ≥ 0, limt→∞ IH (t) ≥ 0, …, limt→∞ V (t) ≥ 0.

Proof. Let us consider (3). By eliminating terms and rearrangement, we obtain Embedded Image where r = βH IH (t)+βCIC(t)+βHCIHC(t) − (µ+θ). Next, we take a function g(t, S(t)) = ln(t, S(t)) and apply Itô’s formula to get, Embedded Image

Combining (6) and (7), solving for S(t) and taking the limit we obtain, Embedded Image

Hence S(t) is positive definite. Similarly, limt→∞ IH (t) ≥ 0, limt→∞ IC(t) ≥ 0, limt→∞ IHC(t) ≥ 0, limt→∞ R(t) ≥ 0, limt→∞ V (t) ≥ 0.

3.2 Invariant Region of Stochastic COVID-19 and Hepatitis B Coinfection Model

An investigation of the invariant region of the model will help us determine the set of states for which the disease will eventually either persist or die out infinitely.

Theorem 2

Let N(t) be the population of system (3), given initial values N (0) > 0, and Embedded Image, then Embedded Image, and (S(t), IH (t), …, V (t)) ∈ Ω.

Proof. By adding and simplifying terms in (3) and choosing ϵ = max(ϵi) we obtain; Embedded Image

Applying (5) we have, Embedded Image

We let Embedded Image. The complementary solution to (3)) was obtained as Embedded Image

Next, we determine the particular solution of (8) and obtain the general solution of (8) as Embedded Image

We then find the limt→∞ of the general equation of (8) and apply comparison to obtain Embedded Image

3.3 Disease Free Equilibrium Point

We consider the first equation of (3). For disease free equilibrium, we have IC(t) = 0, R(t) = 0, and let Embedded Image, solve for S(t), and apply comparison to obtain, Embedded Image

Since 0 < ϵ1(dB1) + ∫ 𝒰 Ψ1(u) < m <∞, where m is a positive constant, the final inequality is obtained by comparison. Hence, the disease free equilibrium is given by Embedded Image

3.4 Reproductive Number of COVID-19 only Model

Theorem 3

Given D.F.E. point Embedded Image, the Reproductive Number of COVID-19 only model Embedded Image.

Proof. Applying Itô-Lévi formula to the COVID-19 infected class of (3), we obtain Embedded Image

Considering initial infections F and secondary infections V, we determine the basic reproduction number by means of the next generation matrix as follows Embedded Image

Taking the spectral radius of F V −1, i.e., ρ(F V −1), we obtain, Embedded Image

3.4.1 Local Stability of COVID-19 Only Model

Theorem 4

If Embedded Image, then for any initial values of Embedded Image, IC(t) satisfies Embedded Image

Proof. By integrating both sides of (11) and evaluating at disease free-equilibrium point point, we have, Embedded Image where Embedded Image are a martingales given Embedded Image. Dividing through by t and taking lim supt→∞ of both sides, we get, Embedded Image

The final inequality holds if Embedded Image, which implies,Embedded Image.

3.4.2 Global Stability of COVID-19 Only Model

Theorem 5

If Embedded Image, then E0C is globally asymptotically stable in Ω.

Proof. We consider the Lyapunov function, Embedded Image For Embedded Image and Q =(γ + μ + ρC) Embedded Image since Embedded Image. Hence L1 is globally stable in the domain Ω.

3.5 Basic Reproduction Number for HBV Model

We consider only HBV transmission in human. The model of the system is given as (3).

Theorem 6

Given D.F.E. point Embedded Image then the Reproductive Number of HBV only Model Embedded Image.

Proof. Applying Itô-Lévi formula to the HBV infected class of (3), we obtain, Embedded Image

Considering initial infections F and secondary infections V, we determine the basic reproduction number by means of the next generation matrix as follows; Embedded Image

Taking the spectral radius of F V −1, i.e., ρ(F V −1), we obtain, Embedded Image

3.5.1 Local Stability of HBV Only Model

Theorem 7

If Embedded Image, then for any initial values of Embedded Image satisfies Embedded Image.

Proof. We integrate both sides of (12) and evaluate at D.F.E. point to obtain, Embedded Image

Next, we divide through by t, find lim sup → ∞ and apply the strong law of Martingales as follows, Embedded Image which completes the proof. Embedded Image as (α + µ + ρH) > 0, which implies, Embedded Image where Embedded Image, and Embedded Image are a martingales given Embedded Image.

3.5.2 Global Stability of HBV Only Model

Theorem 8

If Embedded Image, then E0H is globally asymptotically stable in Ω

Proof. We consider the Lyapunov function L2, Embedded Image

For Embedded Image,and P = (α+ μ+ ρH). Embedded Image hence L2 is globally stable in the domain Ω.

3.6 Reproductive Number of Coinfection Model

Theorem 9

Given D.F.E. point Embedded Image, the Reproductive Number of Coinfection Embedded Image.

Proof. Applying Itô-Lévi formula to the Coinfected equation of (3), we obtain, Embedded Image

Considering initial infections F and secondary infections V, and since Embedded Image Embedded Image

Taking the spectral radius of FV −1, i.e., Embedded Image we obtain, Embedded Image

3.6.1 Local Stability of Disease Free Equilibrium point in the Coinfection Model

Theorem 10

If Embedded Image, then for any initial values of (S(0), IH (0), IC(0), IHC(0), R(0), V (0)) ∈ ℜ6, IHC(t) satisfies Embedded Image.

Proof. Integrating (14), eliminating terms, choosing constants, and setting G = (η + µ + ρHC) yields, Embedded Image

Next, we divide through by t, find lim sup and apply the strong law of Martingales as follows, Embedded Image

The final inequality holds if Embedded Image as (η + µ + ρC) > 0, which implies, Embedded Image where Embedded Image, and Embedded Image are a martingales given Embedded Image.

3.6.2 Global Stability of Coinfection Model

Theorem 11

If Embedded Image, then E0HC is globally asymptotically stable in Ω

Proof. We consider the Lyapunov function, Net we differential and employ comparision to determine conditions for Global Stability of Coinfection. Embedded Image hence L3 is globally stable in the domain Ω.

3.7 Extinction of HBV-COVID-19 Coinfection disease

Next, we determine conditions under which the disease will eventually die out in the population with a probability of one. In this section we study the conditions of extinction for the coinfection model.

Definition 3

[8] We define Embedded Image, and Embedded Image.

Lemma 1

[25] Let (S(t), IC(t), IH (t), IHC(t), R(t), V (t)) be the positive solution of system(4) with given initial condition Embedded Image, Let also X(t) be the positive solution of equation (4) with given initial condition given condition X(0) = N (0) = S(0) + IC(0) + IH (0) + IHC(0) + R(0) + V (0) ∈ ℜ+. Then

  1. Embedded Image

  2. Embedded Image

  3. Embedded Image,and

  4. Embedded Image

Theorem 12

Let (S(t), IH (t), …, V (t)) be the solution of (3) with initial values (S(0), IH (0), …, V (0) ∈ Ω), the coinfection disease of model (3) goes extinct almost surely (limt→∞ IHC(t) = 0) a.s if one of the following assumptions holds:

  1. Embedded Image and Embedded Image,

  2. Embedded Image and Embedded Image,

  3. Embedded Image.

Proof. By integrating (3), dividing by t, adding terms, followed by some algebraic manipulations we obtain the following equations: Embedded Image Embedded Image where Embedded Image

Next, we integrate both sides (14). Since Embedded Image we get, Embedded Image

We divide through (17) by t. Next we substituting (15) and (16) into (17) we obtain, Embedded Image Embedded Image

Next, we find lim supt→∞ and apply a Lemma 2 and Lemma 1 to obtain following results, Embedded Image Embedded Image hence Embedded Image

Similarly, we can prove Embedded Image when Embedded Image, and Embedded Image when Embedded Image.

3.8 Persistence in mean

Persistence of a disease means that the disease will remain endemic in the population with positive probability. We study the disease persistence for the system reported in (3) and derive that the disease persists under certain some conditions.

Definition 4

[8] Persistent for Model (3) will hold if, Embedded Image

Lemma 2

[6] Set g ∈ ℂ [[0, ∞] × Ω, (0, ∞)], assume that there exist ξ0, ξ > 0 such that Embedded Image such that G ∈ (ℂ[[0, ∞] × Ω, (0, ∞)]) satisfying Embedded Image. Then Embedded Image.

Theorem 13

Let (S(t), IH (t), IC(t), IHC(t), R(t), V (t)) be a solution of system (3) with initial values (S(0), IH (0), IC(0), IHC(0), R(0), V (0)) ∈ Ω.

  1. If Embedded Image and Embedded Image, the disease IH persists in mean. In addition, IH holds Embedded Image

  2. If Embedded Image, max (Embedded Imageand Embedded Image) < 1, then the disease IC is persistent in mean. In addition, IC satisfies Embedded Image

Proof. We begin with the first statement of the theorem. By rearranging terms in (15), we get: Embedded Image

Next, we integrate both sides of (12), and eliminate terms to obtain, Embedded Image

We divide (24) by t, substitute (23) into (24), and find lim inft→+∞ to get, Embedded Image

Then, by Lemma 2 and since Embedded Image, we conclude, Embedded Image which concludes the proof. Next we prove the persistence in mean for COVID-19. Let us integrate both sides of (11), and then eliminate terms to obtain, Embedded Image

We divide (25) by t, substitute (23) into (25), and find lim inft→+∞ to get, Embedded Image

Then, by Lemma 2 and since Embedded Image, we conclude, Embedded Image

Finally in a similar way, we can prove persistence for coinfection of diseases; Embedded Image

4 Numerical Results

In this section, we employ a standard numerical procedure In this section, we employ a standard numerical procedure to obtain numerical simulation results for the HBV-only, COVID-19-only, and HBV-COVID-19 coinfection models. The objective here is to verify the theoretical results obtained earlier in this paper. The parameter values used along with their sources are displayed in Table 2. A few values were assumed. As we sought to verify theoretical results, we selected initial values of S(0) = 3000, IH (0) = 10, IC(0) = 10, IHC(0) = 10, R(0) = 0, V (0) = 40 The Euler-Maruyana scheme [12] was employed in conducting the numerical simulations for all models. The results from numerical simulation for HBV only, COVID-19 only and Confection model are presented in Figures 2 - 8.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1: Parameters and description
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2: Parameters Values of Source
Figure 2:
  • Download figure
  • Open in new tab
Figure 2:

Solution for HBV Only Model.

4.1 Hepatitis B Virus Only Model

By using fixed values of all parameters shown in Table 2 we plot the curves representing the dynamics of the susceptible, infectious, and recovered classes in the HBV-only model (3). The graphs of the compartments are presented in Figure 2 reveal. From the parameter values in Table 2, we obtained our basic reproduction numbers Embedded Image.

We observe initially, from Figure 2 that as number of infectious population increase in time, the susceptible population decrease. However, in the long run, the solution curves for both the susceptible and infective compartments behave like decreasing functions and approach a fixed point. A decrease in the number of susceptible individuals occurred within a relatively short period of the epidemic, after which there was no marked change over time. On the other hand, the decline in the number of infective compartments tends to increase but eventually declines to a fixed point. Similarly, the recovered class also approach a fixed point in the long run. In other words, all compartments will reach to the respective endemic fixed point of model (3) within a finite time and in long run.

The numerical outcomes verifies the theorem statement that when Embedded Image, then the disease free equilibrium is locally and globally asymptotically stable. The impact of various parameter on the transmission dynamics of HBV is explored. Numerical results are presented in Figure 3 and Figure 4. We proceeded by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (βH = 0.001, 0.002, 0.003), while all remaining parameter values of the model were held constant. From figure 3, we observe that as values of contact rate increases, there is an increased possibility for the population to be infected by HBV, and vice versa.

Figure 3:
  • Download figure
  • Open in new tab
Figure 3:

HBV Infection with varied effective contact rate.

Figure 4:
  • Download figure
  • Open in new tab
Figure 4:

Recoveries with varied effective recovery rate.

Next, we investigated the effects of recovery rate on population of Recovered individuals (R(t)). Similarly, we held all other parameters constant, and selected varied recovery rate (ρH = 0.05, 0.06, 0.07) for this investigation. We observe Figure 4, that more individuals recover with increasing value in ρH, thus, the reducing the number of infectious individuals. From this we can conclude that by treating the infectious individuals the infected population goes to recovered state and HBV is eliminated from the community.

4.2 COVID-19 Only Model

In Figure 5 simulation results representing the dynamics of Susceptible, Infectious, Vaccinated and Recovered classes in the COVID-19 only model of(3) are displayed. Figure 5 is plotted using combination of all compartments and by taking fixed values of parameters in Table 2. We observe initially, from Figure 5 that as number of infectious population increase in time, susceptible population sharply decrease. However, in the long run, the solution curves for both Susceptible and infective compartments behave like decreasing functions and approach reach to a fixed point. The decline in the number of infective compartment occurs during the after the twentieth time step of the epidemic, after which there is no remarkable change over time, but for fluctuations around the deterministic solution.

Figure 5:
  • Download figure
  • Open in new tab
Figure 5:

Numerical Solution for COVID-19 Only Model.

Similarly, the remaining classes, the vaccinated, and recovered populations also approach fix points in the long run. Thus, all compartments will reach to the respective endemic fixed point of the COVID-19 only model in (3) within a finite time and in long run.

The numerical outcomes verifies the theorem’ statement that when Embedded Image, then the disease free equilibrium is locally and globally asymptotically stable.

The impact of various parameter on the transmission dynamics of COVID-19 was explored. Numerical results are presented in Figure 6 and Figure 7. We proceeded by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (βC = 0.005, 0.006, 0.007), while all remaining parameter values of the model were held constant. Figure 6, shows that as the contact rate increases, the probability that the population will be infected by diseases increases, and vice versa.

Figure 6:
  • Download figure
  • Open in new tab
Figure 6:

COVID-19 Infection for different effective contact rates.

Figure 7:
  • Download figure
  • Open in new tab
Figure 7:

COVID-19 recoveries for different treatment rates.

Next, we examine the effects of recovery rate on population of COVID-19 recoveries RC(t). Again, while we held all other parameters constant, we used various recovery rate (ρC = 0.3, 0.4, 0.5). We observe in Figure 7, that more individuals recover with increasing value in ρHC, thus, the reducing the number of COVID-19 infections in the population.

4.3 Coinfection Model

Figure 8 reveals that the curves representing the dynamics of Susceptible, Infectious, Vaccinated and Recovered classes in the coinfection model (3). Figure 8 is plotted using fixed values of all parameters from Table 2 we obtained our basic reproduction numbers Embedded Image.

Figure 8:
  • Download figure
  • Open in new tab
Figure 8:

Numerical Solution for coInfection Model.

We observe initially, from Figure 8 that as number of infectious population increase in time, susceptible population will decrease, however, in the long run, the solution curves for Susceptible and infective compartments IC, IH, and IHC behave like decreasing functions and approach reach to a fixed point over time, but for oscillation around their deterministic path due to the presence of the noise.

Similarly, the all remaining classes either decrease to zero or also approach to the fixed points in the long run. In other words, all compartments will reach to the respective endemic fixed point of model (3) within a finite time and in long run. The numerical outcomes verifies Theory 10 which states that when Embedded Image, then the disease free equilibrium is locally and globally asymptotically stable.

The impact of various parameter on the transmission dynamics of COVID-19 -HBV coinfection is explored. Numerical results are presented in Figure 9 and Figure 10. We proceed by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (βHC = 0.001, 0.002, 0.003), while all remaining parameter values of the model were held constant. From figure 9, we observe that as values of contact rate increases, there is an increased chance for the population to be infected by both diseases, and vice versa.

Figure 9:
  • Download figure
  • Open in new tab
Figure 9:

Coinfection for different effective contact rates.

Figure 10:
  • Download figure
  • Open in new tab
Figure 10:

Coinfection recoveries for different treatment rates.

Next, we examine the effects of recovery rate on population of Coinfected recoveries RHC(t). Again, while we held all other parameters constant, we made use of varied recovery rate (ρHC = 0.1, 0.15, 0.2). We observe in Figure 10, that more individuals recover with increasing value in ρHC, thus, the reducing the number of coinfectious in the population.

5 Conclusion

The majority of issues in the real world are not deterministic. Because of its proximity to ambient sounds, the stochastic models allow us to model epidemic diseases in a more realistic manner. In this paper, two distinct diseases were used to study a stochastically perturbed SIRS-type model. Emphasis was given to explaining the effect of small fluctuations and perturbations due to sudden environmental shocks in the coinfection transmission dynamics, hence, the inclusion of the white noise and the Lévy noise. To determine the conditions for the local asymptotic stability of the equilibria, we first computed the reproductive number and the equilibria for the underlying stochastic model (3). Furthermore, the conditions of extinction and persistence of HBV only, COVID-19 only, and coinfection model we investigated, and it was shown that persistence for all models were dependent on the intensity of the white noise and the values of epidemic parameters involved in disease transmission.

Additionally, we demonstrated the influence of different parameters on the infectious compartments and demonstrated a broad variety of simulation results for all three stochastic models.

Data Availability

All data produced in the present work are contained in the manuscript

Footnotes

  • This version of the manuscript has been revised to update the following: i) Number of Authors ii) definition 1 and 2 iii) R_{0C}, R_{0H}, R_{0HC} notations in Theorem 12 iv) inclusion of two new citations.

References

  1. [1].↵
    A. I. Abioye, O. J. Peter, H. A. Ogunseye, F. A. Oguntolu, T. A. Ayoola, and A. O. Oladapo. A fractional-order mathematical model for malaria and COVID-19 co-infection dynamics. Healthcare Analytics, 4:100210, 2023.
    OpenUrl
  2. [2].↵
    L. J. Allen. An introduction to stochastic processes with applications to biology. CRC press, 2010.
  3. [3].↵
    P. Boonyarangka, K. Phontham, S. Sriwichai, K. Poramathikul, K. Harncharoenkul, W. Kuntawunginn, N. Maneesrikhum, S. Srisawath, C. Seenuan, C. Thanyakait, K. Inkabajan, S. Pludpiem, K. Pidtana, S. Demons, B. Vesely, M. Wojnarski, J. S. Griesenbeck, and M. Spring. Co-infection with plasmodium vivax and covid-19 in thailand. Tropical Medicine and Infectious Disease, 7(8), 2022.
  4. [4].↵
    X. Chen, Q. Jiang, Z. Ma, J. Ling, W. Hu, Q. Cao, P. Mo, L. Yao, R. Yang, Gao S., X. Gui, W. Hou, Y. Xiong, J. Li, and Y. Zhang. Clinical characteristics of hospitalized patients with sars-cov-2 and hepatitis b virus co-infection. Virol Sin., 35:842–845, 2020.
    OpenUrl
  5. [5].↵
    G. Contreras, M. R. Mora, and P. J. Gómez. Estimation of the epidemiological evolution through a modelling analysis of the covid-19 outbreak. Microbiology, 3:152–158, 2020.
    OpenUrl
  6. [6].↵
    P. Das, R. K. Upadhyay, A. K. Misra, F. A. Rihan, P. Das, and D. Ghosh. Mathematical model of covid-19 with comorbidity and controlling using non-pharmaceutical interventions and vaccination. Nonlinear Dyn., 106:1213–1227, 2021.
    OpenUrl
  7. [7].↵
    M. L. Diagne, H. Rwezaura, S. Y. Tchoumi, and J. M. Tchuenche. A Mathematical Model of COVID-19 with Vaccination and Treatment. Computational and Mathematical Methods in Medicine, 2021:1250129, 2021.
    OpenUrl
  8. [8].↵
    A. Din, S. Amine, and A. Allali. A stochastically perturbed co-infection epidemic model for covid-19 and hepatitis b virus. Nonlinear Dynamics, 111(2), 2023.
  9. [9].↵
    M. F. Hama, R. R. Q. Rasul, Z. Hammouch, K. A. H. Rasul, and J. Danane. Analysis of a stochastic seis epidemic model with the standard brownian motion and lévy jump. Results in Physics, 37:105477, 2022.
    OpenUrl
  10. [10].↵
    I. M Hezam, A. Foul, and A. Alrasheedi. A dynamic optimal control model for COVID-19 and cholera co-infection in yemen. Adv. Differ. Equ., 2021(1), December 2021.
  11. [11].↵
    E. A. Iboi, C. N. Ngonghala, and A. B. Gumel. Will an imperfect vaccine curtail the covid-19 pandemic in the us? Infectious Disease Modelling, 5:510–524, 2020.
    OpenUrl
  12. [12].↵
    T. Jahnke, T. Udrescu, and M. Mikl. Numerical methods in mathematical finance, 2020.
  13. [13].↵
    C. Ji, D. Jiang, and N. Shi. The behavior of an sir epidemic model with stochastic perturbation. Stochastic Analysis and Applications, 30(5):755–773, 2012.
    OpenUrl
  14. [14].↵
    K. G. Mekonen and L. L. Obsu. Mathematical modeling and analysis for the co-infection of COVID-19 and tuberculosis. Heliyon, 8(10):e11195, October 2022.
    OpenUrl
  15. [15].↵
    E. F. Mohamed, E. K. Mohamed, G. Richard, L. Aziz, and T. Regragui. Stationary distribution and threshold dynamics of a stochastic sirs model with a general incidence. Physica A-statistical Mechanics and Its Applications, 534:120696, 2019.
    OpenUrl
  16. [16].↵
    S. E. Moore, H. L. Nyandjo-Bamen, O. Menoukeu-Pamen, J. K. K. Asamoah, and Z. Jin. Global stability dynamics and sensitivity assessment of covid-19 with timely-delayed diagnosis in ghana. Computational and Mathematical Biophysics, 10(1):87–104, 2022.
    OpenUrl
  17. [17].↵
    M. M. Ojo and D. F. E. Goufo. The impact of covid-19 on a malaria dominated region: A mathematical analysis and simulations. Alexandria Engineering Journal, 65:23–39, 2023.
    OpenUrl
  18. [18].↵
    World Health Organization. COVID-19 vaccination in humanitarian settings. https://www.who.int/publications/i/item/9789240079434, 2023. [Online; accessed 25 October 2023].
  19. [19].↵
    M. A. Pobbi, S. M. Naandam, and S. E. Moore. Mathematical modelling and analysis of stochastic malaria and covid-19 co-infection model. doi:10.21203/rs.3.rs-3591875/v1, 2023.
    OpenUrlCrossRef
  20. [20].↵
    A. Pusparani, J. Henrina, and A. Cahyadi. Co-infection of covid-19 and recurrent malaria. The Journal of Infection in Developing Countries, 15(05):625–629, 2021.
    OpenUrl
  21. [21].↵
    M. Ray, A. Vazifdar, and S. Shivaprakash. Co-infection with malaria and coronavirus disease-2019. Journal of Global Infectious Diseases, 12(3), 2020.
  22. [22].↵
    J. M. Read, J. R. E. Bridgen, D. A. T. Cummings, A. Ho, and C. P Jewell. Novel coronavirus 2019-ncov (covid-19): early estimation of epidemiological parameters and epidemic size estimates. Philosophical Transactions of the Royal Society B, 376(1829):20200265, 2021.
    OpenUrl
  23. [23].↵
    A. Rehman, R. Singh, and P. Agarwal. Modeling, analysis and prediction of new variants of covid-19 and dengue co-infection on complex network. Chaos Solitons Fractals, 150(111008):111008, September 2021.
    OpenUrl
  24. [24].↵
    H. Rwezaura, M. L. Diagne, A. Omame, A. L. de Espindola, and J. M. Tchuenche. Mathematical modeling and optimal control of sars-cov-2 and tuberculosis co-infection: a case study of indonesia. Modeling Earth Systems and Environment, 8:5493–5520, 2022.
    OpenUrl
  25. [25].↵
    Y. Sabbar. Mathematical Analysis of Some Stochastic Infectious Disease Models with White Noises and Lévy Jumps. PhD thesis, Université Sidi Mohamed Ben Abdellah de Fès (Maroc), 2021.
  26. [26].↵
    S. M. A. Shah, Y. Nie, A. Din, and A. Alkhazzan. Dynamics of hepatitis b virus transmission with a lévy process and vaccination effects. Electron. J. Differ. Equ., 12:16–45, 2024.
    OpenUrl
  27. [27].
    L. Tang, Y. Zhou, L. Wang, S. Purkayastha, L. Zhang, J. … He, and P. X. K. Song. A review of multi-compartment infectious disease models. International Statistical Review, 88, 2020.
  28. [28].↵
    S. Y. Tchoumi, M. L. Diagne, H. Rwezaura, and J. M. Tchuenche. Malaria and covid-19 co-dynamics: A mathematical model and optimal control. Applied mathematical modelling, 99:294–327, 2021.
    OpenUrl
  29. [29].↵
    S. W. Teklu. Analysis of hbv and covid-19 coinfection model with intervention strategies. Comput Math Methods Med., 2023.
  30. [30].↵
    A. W. Tesfaye and T. S. Satana. Stochastic model of the transmission dynamics of covid-19 pandemic. Advances in Difference Equations, 2021:1–21, 2021.
    OpenUrl
  31. [31].
    UNICEF. How COVID-19 caused malaria to soar – and what countries and the global community should continue to do in response. https://data.unicef.org/data-for-action/how-covid-19-caused-malaria-to-soar, 2023. [Online; accessed 25 October 2023].
  32. [32].↵
    J. Wang and X. Tian. Global stability of a delay differential equation of hepatitis b virus infection with immune. Electron. J. Differ. Equ., 94:1–11, 2013.
    OpenUrl
  33. [33].↵
    X. Wang, C. Wang, and K. Wang. Extinction and persistence of a stochastic sica epidemic model with standard incidence rate. Advances in Difference Equations, 2021.
  34. [34].↵
    J. T. Wu, K. Leung, and G. M. Leung. Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: a modelling study. The Lancet, 395(10225):689–697, 2020.
    OpenUrlCrossRef
  35. [35].↵
    Z. Zhang, A. Zeb, S. Hussain., and E. Alzahrani. Dynamics of COVID-19 mathematical model with stochastic perturbation. Advances in Difference Equations, 2020(1):451, 2020.
    OpenUrl
Back to top
PreviousNext
Posted August 14, 2024.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Analysis of Stochastic COVID-19 and Hepatitis B Co-infection Model with Brownian and Lévy Noise
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Analysis of Stochastic COVID-19 and Hepatitis B Co-infection Model with Brownian and Lévy Noise
Michael A. Pobbi, Stephen E. Moore
medRxiv 2024.08.12.24311861; doi: https://doi.org/10.1101/2024.08.12.24311861
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Analysis of Stochastic COVID-19 and Hepatitis B Co-infection Model with Brownian and Lévy Noise
Michael A. Pobbi, Stephen E. Moore
medRxiv 2024.08.12.24311861; doi: https://doi.org/10.1101/2024.08.12.24311861

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)