Next Article in Journal
Multiobjective Optimization of Fuzzy System for Cardiovascular Risk Classification
Next Article in Special Issue
Diffusion Kinetics Theory of Removal of Assemblies’ Surface Deposits with Flushing Oil
Previous Article in Journal
Simultaneous Integration of D-STATCOMs and PV Sources in Distribution Networks to Reduce Annual Investment and Operating Costs
Previous Article in Special Issue
Extended Online DMD and Weighted Modifications for Streaming Data Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Dynamics of Tuberculosis in Algeria Using a Compartmental VSEIT Model with Evaluation of the Vaccination and Treatment Effects

by
Bouchra Chennaf
1,
Mohammed Salah Abdelouahab
1,* and
René Lozi
2
1
Laboratory of Mathematics and Their Interactions, Abdelhafid Boussouf University Center, Mila 43000, Algeria
2
Laboratory J.A. Dieudonné, CNRS, Université Côte d’Azur, 06108 Nice, France
*
Author to whom correspondence should be addressed.
Computation 2023, 11(7), 146; https://doi.org/10.3390/computation11070146
Submission received: 9 June 2023 / Revised: 5 July 2023 / Accepted: 16 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Mathematical Modeling and Study of Nonlinear Dynamic Processes)

Abstract

:
Despite low tuberculosis (TB) mortality rates in China, Europe, and the United States, many countries are still struggling to control the epidemic, including India, South Africa, and Algeria. This study aims to contribute to the body of knowledge on this topic and provide a valuable tool and evidence-based guidance for the Algerian healthcare managers in understanding the spread of TB and implementing control strategies. For this purpose, a compartmental mathematical model is proposed to analyze TB dynamics in Algeria and investigate the vaccination and treatment effects on disease breaks. A qualitative study is conducted to discuss the stability property of both disease-free equilibrium and endemic equilibrium. In order to adopt the proposed model for the Algerian case, we estimate the model parameters using Algerian TB-reported data from 1990 to 2020. The obtained results using the proposed mathematical compartmental model show that the reproduction number ( R 0 ) of TB in Algeria is less than one, suggesting that the disease can be eradicated or effectively controlled through a combination of interventions, including vaccination, high-quality treatment, and isolation measures.

1. Introduction

The recent outbreak of the coronavirus disease, known as COVID-19 has indeed highlighted the critical role of epidemic research, particularly mathematical modeling, in understanding and combating infectious diseases, since it provides a powerful tool to analyze the dynamic transmission of diseases and assess the potential impact of various interventions and control measures.
Tuberculosis is a contagious infection caused by bacteria called Mycobacterium tuberculosis that primarily affects the lungs. It can also spread to other body parts, including the brain and spine. Importantly, it can be contracted not only through direct contact with an infected individual but also through the inhalation of airborne droplets containing the bacteria.
TB is one of the top 10 killers worldwide and causes 1.8 million deaths each year. Of all new TB cases recorded in 2020, 86 % occurred in the 30 countries with the highest disease burden. Two-thirds of the cases are concentrated in eight countries, with India leading, followed by China, Indonesia, Philippines, Pakistan, Nigeria, Bangladesh and South Africa (According to the World Health Organization (WHO)) [1]. This demonstrates that TB poses a threat to human health and has a detrimental impact on social and economic life.
Although Algeria may not be among the top eight countries with the highest concentration of TB cases globally, it is still a significant concern in Algeria. Thus, it is imperative that government agencies and scientists work together to manage and combat the spread of TB epidemics.
Mathematical modeling plays a critical role in the planning and implementation of TB control programs. Although Bernoulli used mathematical models for smallpox in 1760 [2], the research on infectious diseases using deterministic mathematical models actually started in the 20th century. Other major works in mathematical epidemiology are due to P.D En’ko between 1873 and 1894. However, it can be said that the foundations of mathematical epidemiology based on compartmental models are due to Sir Ronald Ross, who gave the first mathematical model of malaria transmission in 1911 [3]. The early research in this field are available in [4,5,6].
The main difference between compartmental models and other models of diseases is that compartmental models explicitly consider the different stages of disease progression and the transitions between them. This allows for a more detailed understanding of how diseases spread and how interventions can be implemented to control their spread [7]. Other models, such as statistical models or network models, may not explicitly include this level of detail.
Susceptible-Infected-Recovered (SIR) is a deterministic model that Kermack and McKendrick proposed in 1927 to characterize the behavior of epidemic spread [8]. Despite the fact that this model has been used successfully to represent the behavior of disease, it is unrealistic by ignoring other compartments and control techniques, such as vaccination, treatment, isolation, and the impact of age and sex.
Epidemiological compartmental models can be broadly classified into two categories: differential equation-based models that describe the dynamics of infectious diseases using continuous functions which can capture the continuous changes in the state variables over time, typically represented by systems of ordinary differential equations and difference equation-based models, often used when data are collected at discrete time points or when the population dynamics are better captured in a discrete manner.
The first mathematical model of TB was developed in 1962 by Waaler and Anderson, who divided the entire population into different groups [9]. Since then, numerous academics have created various mathematical models to investigate and control TB in countries heavily impacted by the disease; see for instance [10,11]. These models have been instrumental in guiding public health policies and interventions aimed at reducing the burden of TB.
Vaccination is one of the most vital factors in stopping and controlling the spread of TB. The tuberculosis vaccination against Bacillus Calmette–Guérin (BCG) was first given to a human in 1921. The World Health Organization (WHO) currently advises immunizing newborns with a single intradermal injection of BCG as soon as possible after the birth [12]. To address the prevailing epidemiological situation, the Algerian Health Care Administration implemented a dedicated vaccination schedule and mandatory vaccination campaigns for children. As a result, the BCG vaccination coverage reached a remarkable rate of over 98% across newborns. However, BCG vaccination is not typically recommended for adults, as its effectiveness in this age group is limited. Therefore, this paper neglects it.
The mathematical modeling of tuberculosis relies on vaccination for its importance in giving predictions to eradicate the disease. There are many previous studies concerned with this topic; for example, in [13] the goal was to determine the dynamics of tuberculosis in Turkey, and the impact of vaccine therapy on the disease. Yang et al. [7] formulated a mathematical model to investigate the effects of immunization and treatment on the dynamics of tuberculosis transmission. Egonmwan et al. [14] developed a mathematical model that includes immunization of newborn children and older susceptible people in the dynamics of TB transmission in a population, with the goal of providing protection to older susceptible people. Revelle et al. [15] formulated models for the economic allocation of activities to control tuberculosis in developing countries.
To the best of our knowledge, the proposed model is not considered elsewhere in its present form and there is no research on modeling the dynamic transmission of tuberculosis in Algeria using a compartmental model while simultaneously estimating the relevant biological parameters specific to the country. Therefore, conducting research in this domain would make a valuable contribution to the field of TB modeling and control in Algeria.
In this study, we propose a VSEIT epidemiological model to investigate the dynamics of TB disease in Algeria. To confirm its performance we estimated the biological model’s parameters using specific TB data, including disease incidence, prevalence, and other relevant epidemiological information from 1990 to 2020 from the WHO Global TB Report [1].
This paper is organized as follows: Section 2 presents the formulation of the VSEIT TB model and analyzes its dynamic properties. In Section 3, the estimation of model parameters is conducted, along with their sensitivity analysis. Section 4 is focused on the discussion of the obtained results. Finally, the paper is concluded in Section 5.

2. Mathematical Model and Dynamic Analysis

In this section, we present the proposed mathematical model for TB infection, and we examine its dynamic.

2.1. Model Formulation

The population is divided into five distinct subgroups: vaccinated individuals ( V ) , susceptible individuals ( S ) , exposed or exposed individuals ( E ) , infected individuals with active TB ( I ) , and individuals currently receiving treatment ( T ) . Hence, the total population is
N ( t ) = V ( t ) + S ( t ) + E ( t ) + I ( t ) + T ( t ) .
This model aims to provide a comprehensive understanding of the spread and progression of TB within a population, allowing for more effective prevention and treatment strategies to be developed.
The number of people that have received vaccination ( V ) is increased through a small proportion of immunized newborns, p Λ . The vaccinated population decreases as vaccinated individuals become susceptible at a rate k (the vaccine’s efficacy wanes over time), during the protection period, they will not become infected even if they contact infected individuals as long as the vaccination provides immunity to all of them. The natural death rate in the class V is μ . Hence, the population of vaccinated individuals is given by the first equation in system (1).
The population of susceptible individuals ( S ) is increased by the small proportion of newborns who are not immunized from TB, ( 1 p ) Λ , and also increases as vaccinated individuals become susceptible, at a rate k. This population decreases when there is contact with infected people, at a rate β . As older people die naturally at a rate μ , the population of susceptible individuals decreases. Hence, the population of susceptible individuals is given by the second equation in system (1).
We assume that the population of exposed individuals E increases when the susceptible population makes effective contact with infected individuals and decreases as latently infected peoples progress from exposed to active TB, at a rate ϵ , and die naturally, at a rate μ . The population of exposed individuals E is also increased by an inflow of a fraction, δ ( 1 α ) of individuals under treatment, where the parameter α represents the treatment failure rate. Particularly, α = 0 means that all treated individuals will move to a exposed state, whereas α = 1 , means that the treatment has failed, and all treated individuals will remain infectious. Hence, the population of exposed individuals is given by the third equation in system (1).
The population of infected individuals ( I ) increases as latently infected individuals progress from exposed to active TB, and the effectively treated patients return to active TB at a rate α δ , which significantly increases the population I. As infected people receive treatment, the population I decreases at a rate γ . Both natural death and TB disease kill people at rates σ and μ , respectively. Hence, the population of infected individuals is given by the fourth equation in system (1).
Finally, as infected people are treated, the population of treated individuals ( T ) grows at a rate γ . As individuals who have been successfully treated become reinfected, the population T decreases at a rate δ . The population continues to decline due to natural mortality ( μ ) and deaths from TB ( η ). Hence, the population of treated individuals is given by the fifth equation in system (1).
The dynamic of TB infection is described by the following system of differential equations:
d V ( t ) d t = p Λ ( k + μ ) V ( t ) , d S ( t ) d t = ( 1 p ) Λ + k V ( t ) β S ( t ) I ( t ) μ S ( t ) , d E ( t ) d t = β S ( t ) I ( t ) ( ϵ + μ ) E ( t ) + ( 1 α ) δ T ( t ) , d I ( t ) d t = ϵ E ( t ) + α δ T ( t ) ( γ + μ + σ ) I ( t ) , d T ( t ) d t = γ I ( t ) ( μ + δ + η ) T ( t ) .
With: V ( 0 ) 0 , S ( 0 ) 0 , E ( 0 ) 0 , I ( 0 ) 0 and T ( 0 ) 0 with N ( 0 ) > 0 .
The flowchart of the model is shown in Figure 1. The model variables are presented in Table 1 and the model parameters are presented in Table 2.
The values of parameters in Table 2 will be given in Section 3.1.

2.2. Model Analysis

In this subsection the proposed model (1), will be qualitatively analyzed.

2.2.1. Invariance of the Feasible Region

The TB model (1) will be studied in a biologically feasible region Ω R + 5 given by
Ω = ( V ( t ) , S ( t ) , E ( t ) , I ( t ) , T ( t ) ) R + 5 : N ( t ) Λ μ .
 Lemma 1.
For all t > 0 and non-negative initial conditions, the solution of TB model (1) is positive whenever it exists. Furthermore, if 0 N ( 0 ) Λ μ , then
0 N ( t ) Λ μ , for all t > 0 .
 Proof. 
For positive values of V ( t ) ,   S ( t ) ,   E ( t ) ,   I ( t ) and T ( t ) we have
V V = 0 = p Λ 0 , S S = 0 = ( 1 p ) Λ + k V 0 , E E = 0 = β S I + ( 1 α ) δ T 0 , I I = 0 = ϵ E + α δ T 0 , T T = 0 = γ I 0 .
Hence, for non-negative initial conditions, the solution remains positive ∀ t 0 .
It follows from the addition of the VSEIT model Equations (1) that
d N ( t ) d t = Λ μ V ( t ) + S ( t ) + E ( t ) + I ( t ) + T ( t ) σ I ( t ) + η T ( t ) ,
                                                  = Λ μ N ( t ) ( σ I ( t ) + η T ( t ) ) Λ μ N ( t ) .
For N ( t ) Λ μ , we have d N ( t ) d t 0 .
  • Thus, for 0 N ( 0 ) Λ μ , we obtain 0 N ( t ) Λ μ for all t 0 .
It follows that the feasible region Ω is positively invariant. □

2.2.2. Equilibrium Points and Their Stability

To find equilibrium points of the model (1), one solves the equations:
d V d t = d S d t = d E d t = d I d t = d T d t = 0 .
  • One gets two equilibrium points:
  • The disease-free equilibrium point “DFE”
    E 1 = ( V 1 * , S 1 * , E 1 * , I 1 * , T 1 * ) = p Λ k + μ , ( k + μ μ p ) Λ μ ( k + μ ) , 0 , 0 , 0 ,   at   which   N = V 1 * + S 1 * = Λ μ .
    and the endemic equilibrium point “EE”
E 2 = ( V 2 * , S 2 * , E 2 * , I 2 * , T 2 * ) = p Λ k + μ , ( k + μ μ p ) Λ ( k + μ ) ( β I 2 * + μ ) , ( γ + μ + σ ) ( μ + δ + η ) α δ γ ϵ ( μ + δ + η ) I 2 * , I 2 * , γ μ + δ + η I 2 * ,
where:
I 2 * = ( k + μ μ p ) ϵ Λ ( μ + δ + η ) ( k + μ ) ( ( ϵ + μ ) ( γ + μ + σ ) ( μ + δ + η ) ( ϵ + μ ) α δ γ ( 1 α ) γ δ ϵ ) μ β . = μ β ( R 0 1 ) .
Then, the endemic equilibrium point E 2 , exists for R 0 > 1 , and at this equilibrium, one has N = Λ ( σ I 2 * + η T 2 * ) μ < Λ μ .

2.2.3. The Basic Reproduction Number R 0

The basic reproduction number, denoted R 0 , is the expected number of secondary cases produced in a completely susceptible population, by a typical infectious individual during its infective period [16].
If R 0 < 1 , the disease will not be able to spread among the population, whereas, if R 0 > 1 , the disease has the potential to spread among the population and become endemic.
Using the next generation matrix one gets the basic reproduction number R 0 for the proposed model (1) as
R 0 = ϵ ( k + μ μ p ) Λ β k 3 μ ( k + μ ) ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) .
The details of the calculation of R 0 are presented in the Appendix A and Appendix B.

2.2.4. Local Stability Analysis of DFE

 Theorem 1.
The disease-free equilibrium point E 1 (DFE) of model (1) is locally asymptotically stable if R 0 < 1 and unstable if R 0 > 1 .
 Proof. 
The Jacobian matrix J E 1 of system (1) at the DFE E 1 is given by
J E 1 = ( k + μ ) 0 0 0 0 k μ 0 β S 1 * 0 0 0 k 1 β S 1 * ( 1 α ) δ 0 0 ϵ k 2 α δ 0 0 0 γ k 3
The characteristic equation of J E 1 is
( ( k + μ ) λ ) ( μ λ ) [ λ 3 + a 1 λ 2 + a 2 λ + a 3 ] = 0 ,
where
a 1 = [ k 1 + k 2 + k 3 ] , a 2 = [ k 1 k 2 + k 1 k 3 + k 2 k 3 + α δ γ + ϵ β S 1 * ] , a 3 = [ ϵ β S 1 * k 3 + k 1 k 2 k 3 k 1 α δ γ ϵ ( 1 α ) δ γ ] = [ ϵ β S 1 * k 3 + ϵ β S 1 * k 3 R 0 ] = ϵ β S 1 * k 3 ( 1 R 0 1 ) .
Then, all the eigenvalues of the characteristic Equation (6) have a negative real part if the coefficients a i , i = 1 , 2 , 3 fulfill the Routh–Hurwitz conditions, which are a 1 > 0 , a 3 > 0 and a 1 a 2 a 3 > 0 .
  • Hence, the disease-free equilibrium of model (1) is locally asymptotically stable, providing that R 0 < 1 . □

2.2.5. Global Stability Analysis of DFE

 Theorem 2.
The disease-free equilibrium point E 1 (DFE) of model (1) is globally asymptotically stable if R 0 < 1 and unstable if R 0 > 1 .
 Proof. 
To prove the theorem, we consider the following Lyapunov function:
W ( V , S , E , I , T ) = b 1 E + b 2 I + b 3 T ,
where b i , for i = 1 , 2 , 3 , are positive constants to be chosen later. Calculating the derivative of W with respect to time along the solutions of system (1), we obtain:
d W d t = b 1 d E d t + b 2 d I d t + b 3 d T d t = b 1 β S I k 1 E + ( 1 α ) δ T + b 2 ϵ E + α δ T k 2 I + b 3 γ I k 3 T b 1 Λ β μ I k 1 E + ( 1 α ) δ T + b 2 ϵ E + α δ T k 2 I + b 3 γ I k 3 T , because S Λ μ = b 1 Λ β μ + b 3 γ b 2 k 2 I + b 2 ϵ b 1 k 1 E + b 1 ( 1 α ) δ + b 2 α δ b 3 k 3 T ( b 2 k 2 b 3 γ ) ( k + μ μ p ) ( k + μ ) b 1 Λ β μ ( b 2 k 2 b 3 γ ) 1 I + b 2 ϵ b 1 k 1 E + b 1 ( 1 α ) δ + b 2 α δ b 3 k 3 T .
Choosing b 1 = ϵ k 3 ( k + μ μ p ) ( k + μ ) , b 2 = k 1 k 3 ( k + μ μ p ) ( k + μ ) , and b 3 = ( 1 α ) ϵ + α δ k 1 ( k + μ μ p ) ( k + μ ) , one obtains
d W d t b 1 Λ β μ R 0 ( R 0 1 ) I .
Hence, if R 0 < 1 , then d W d t is negative. By LaSalle’s invariant principle (A1), this implies that E 1 is globally asymptotically stable. □

2.2.6. Local Stability Analysis of EE

The Jacobian matrix J E 2 of system (1) at the EE E 2 is given by
J E 2 = ( k + μ ) 0 0 0 0 k β I 2 * μ 0 β S 2 * 0 0 β I 2 * k 1 β S 2 * ( 1 α ) δ 0 0 ϵ k 2 α δ 0 0 0 γ k 3 .
The characteristic equation of J E 2 is
( ( k + μ ) λ ) [ λ 4 + a 1 λ 3 + a 2 λ 2 + a 3 λ + a 4 ] = 0 ,
where
a 1 = ( μ + k 1 + k 2 + k 3 + β I * ) , a 2 = [ k 1 ( μ + β I * ) + k 3 ( μ + k 1 + k 2 + β I * ) + k 2 ( μ + k 1 + β I * ) γ α δ ϵ β S * ] , a 3 = [ ϵ ( β 2 I * S * + k 1 β S * ) + k 3 ( k 2 ( μ + k 1 + β I * ) + k 1 ( μ + β I * ) β ϵ S * ) + k 1 k 2 ( μ + β I * ) + γ ( α δ k 2 + ϵ δ ( α 1 ) ) γ α δ ( μ + k 1 + k 2 + β I * ) β ϵ S * ( μ + k 1 + β I * ) ] , a 4 = [ γ α δ ( k 2 ( μ + k 1 + β I * ) + k 1 ( μ + β I * ) β ϵ S * ) + k 3 ( k 1 k 2 ( μ + β I * ) + ϵ ( β 2 I * S * + k 1 β S * ) + β ϵ S * ( μ + k 1 + β I * ) ) γ ( α δ k 2 δ ϵ ( α 1 ) ) ( μ + k 1 + k 2 + β I * ) γ ( ϵ ( δ k 1 ( α 1 ) + β α δ S * ) ) + k 2 ( α δ k 2 + δ ϵ ( α 1 ) ) ] .
Then, all the eigenvalues of the characteristic Equation (7) have negative real parts if the coefficients a i , i = 1 , 2 , 3 , 4 fulfill the Routh–Hurwitz conditions, which are a i > 0 for i = 1 , 2 , 3 , 4 and a 1 a 2 a 3 > a 3 2 + a 1 2 a 4 . Hence, the endemic equilibrium of model (1) is locally asymptotically stable if R 0 > 1 .

2.2.7. Global Stability Analysis of EE

In this section, we prove the global asymptotic stability of the EE of model (1). Using the method described in [17], at the EE, from system (1) we obtain
p Λ = ( k + μ ) V 2 * , k V 2 * = ( 1 p ) Λ + β S 2 * I 2 * + μ S 2 * , k 1 E 2 * = β S 2 * I 2 * + ( 1 α ) δ T 2 * , k 2 I 2 * = ϵ E 2 * + α δ T 2 * , γ I 2 * = k 3 T 2 * .
 Theorem 3.
The endemic equilibrium point E 2 (EE) of model (1) is globally asymptotically stable if R 0 > 1 .
 Proof. 
To prove the theorem, we consider the following Lyapunov function:
W ( V , S , E , I , T ) = k V ( t ) V * V * ln V ( t ) V * + ϵ S ( t ) S * S * ln S ( t ) S * + ϵ E ( t ) E * E * ln E ( t ) E * + k 1 I ( t ) I * I * ln I ( t ) I * + δ T * k 1 α + ϵ ( 1 α ) γ I * T ( t ) T * T * ln T ( t ) T * .
Calculating the derivative of W with respect to time along the solutions of system (1), we obtain:
d W d t = k 1 V * V V + ϵ 1 S * S S + 1 E * E E + k 1 1 I * I I + δ T * k 1 α + ϵ ( 1 α ) γ I * 1 T * T T .
A straightforward simple calculation gives
k 1 V * V V = k 1 V * V p Λ ( k + μ ) V = k 1 V * V ( k + μ ) V * ( k + μ ) V = k ( k + μ ) V * 1 V * V 1 V V * = k ( k + μ ) V * 2 V * V V V * .
ϵ 1 S * S S = ϵ 1 S * S ( 1 p ) Λ + k V β S I μ S = ϵ μ S * 2 S * S S S * + ϵ β S * I * 1 S * S + I I * ϵ β S I .
ϵ 1 E * E E = ϵ 1 E * E β S I k 1 E + ( 1 α ) δ T = ϵ β S I ϵ β S I E * E k 1 ϵ E + k 1 ϵ E * + ( 1 α ) ϵ δ T ( 1 α ) δ ϵ T E * E = ϵ β S I ϵ β S I E * E k 1 ϵ E + ϵ β S * I * + ( 1 α ) δ ϵ T * + ( 1 α ) ϵ δ T ( 1 α ) δ ϵ T E * E .
k 1 1 I * I I = k 1 1 I * I ϵ E + α δ T k 2 I = k 1 ϵ E k 1 ϵ E I * I + k 1 α δ T k 1 α δ T I * I k 1 k 2 I + k 1 k 2 I * = k 1 ϵ E k 1 ϵ E * I * E I E * + k 1 α δ T k 1 α δ T I * I + ϵ β S * I * + k 1 α δ T * + ϵ ( 1 α ) δ T * ϵ β S * I * I I * ( 1 α ) ϵ δ T * I I * k 1 α δ T * I I * = k 1 ϵ E ϵ β S * I * E I * E * I ( 1 α ) ϵ δ T * E I * E * I + k 1 α δ T k 1 α δ T I * I + ϵ β S * I * + k 1 α δ T * + ϵ ( 1 α ) δ T * ϵ β S * I * I I * ( 1 α ) ϵ δ T * I I * k 1 α δ T * I I * .
δ T * k 1 α + ϵ ( 1 α ) γ I * 1 T * T γ I k 3 T = δ T * k 1 α + ϵ ( 1 α ) I * 1 T * T I I * T * T = δ α k 1 T * I I * + δ ( 1 α ) ϵ T * I I * δ α k 1 T * I T * I * T δ ( 1 α ) ϵ T * I T * I * T δ α k 1 T δ ( 1 α ) ϵ T + δ α k 1 T * + δ ( 1 α ) ϵ T * .
Using Equations (8)–(12), one gets
d W d t = k ( k + μ ) V * 2 V * V V V * + ϵ μ S * 2 S * S S S * + ϵ β S * I * 3 S * S I I * E I * E * I S I E * S * I * E 1 E S * E * S + ( 1 α ) ϵ δ T * 3 I * E I E * E * T E T * T * I T I * + δ α k 1 T * 2 I * T I T * I T * I * T .
Using the properties of the geometric and arithmetic means in Equation (13), one obtains
2 V V * V * V 0 , 2 S S * S * S 0 , 3 S * S I I * E I * E * I S I E * S * I * E 1 E S * E * S 0 , 3 I * E I E * E * T E T * T * I T I * 0 , 2 I * T I T * I T * I * T 0 .
Since none of the parameters are negative, it follows that d W d t 0 when R 0 > 1 . As a result, according to LaSalle’s Invariance Principle (A1), ( V , S , E , I , T ) V * , S * , E * , I * , T * as t . □

2.2.8. Transcritical Bifurcation Analysis

Here, we discuss the existence of transcritical bifurcation of system (1). At R 0 = 1 , if we take β as a bifurcation parameter, we obtain
β * = β = μ ( k + μ ) ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) ϵ ( k + μ μ p ) Λ k 3 .
The following modification are made in the variables ofsystem (1) so that V = x 1 , S = x 2 , L = x 3 , I = x 4 , and T = x 5 . Further using vector notation x = ( x 1 , x 2 , x 3 , x 4 , x 5 ) T , model (1) can then be written in the form d x d t = F , with F = ( f 1 , f 2 , f 3 , f 4 , f 5 ) T as shown below
x ˙ 1 = p Λ ( k + μ ) x 1 , x ˙ 2 = ( 1 p ) Λ + k x 1 β x 2 x 4 μ x 2 , x ˙ 3 = β x 2 x 4 ( ϵ + μ ) x 3 + ( 1 α ) δ x 5 , x ˙ 4 = ϵ x 3 + α δ x 5 ( γ + μ + σ ) x 4 , x ˙ 5 = γ x 4 ( μ + δ + η ) x 5 .
with N = i = 1 5 x i .
The Jacobian matrix evaluated at the disease-free equilibrium E 1 (DFE) for β = β * is
J E 1 = ( k + μ ) 0 0 0 0 k μ 0 β * S * 0 0 0 k 1 β * S * ( 1 α ) δ 0 0 ϵ k 2 α δ 0 0 0 γ k 3
The Jacobian matrix J E 1 has a simple zero eigenvalue calculated at β * .
The right and left eigenvectors denoted by: Y = ( y 1 , y 2 , y 3 , y 4 , y 5 ) and Z = ( z 1 , z 2 , z 3 , z 4 , z 5 ) , respectively, are obtained as follows
y 1 = 0 ,   y 2 = ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) μ ϵ k 3 y 4 ,   y 3 = k 1 k 2 k 3 α γ δ k 1 ϵ k 3 y 4 ,   y 4 > 0 ,   y 5 = γ k 3 y 4 .
and
z 1 = 0 ,   z 2 = 0 ,   z 3 > 0 ,   z 4 = k 1 ϵ z 3 ,   z 5 = ( 1 α ) δ ϵ + α δ k 1 ϵ k 3 z 3 .
We have
Y t D β F E 1 , β * = 0 , Y t D x D β F E 1 , β * Z = S * y 2 z 4 + S * y 3 z 4 , Y t D x 2 F E 1 , β * ( Z , Z ) = ( β * y 2 z 4 , β * y 2 z 4 ) .
The following conditions are satisfied:
Y t D β F E 1 , β * = 0 , Y t D x D β F E 1 , β * Z 0 , Y t D x 2 F E 1 , β * ( Z , Z ) 0 .
Hence, there is a transcritical bifurcation at β = β * as illustrated in Figure 2.

3. Parameters Estimation and Numerical Simulation

In this section, six model parameters will be estimated based on TB incidence data from the WHO Global TB Report [1] between 1990 and 2020 (see Table 3) and the other parameters will be inspired by the statistical data in the literature.

3.1. Parameters Estimation

Using the data of Algeria’s population from [18], one takes the death rate μ as the average death rate per year from 1990 to 2020, μ = 0.00498 , and the recruitment rate Λ , as the average birth per year from 1990 to 2020, Λ = 811 , 085 .
The child immunization rate, BCG, is the ratio of children aged 12–23 months who have received BCG vaccination. Figure 3 displays the percentage of one-year-old children who have received the BCG immunization in Algeria between 1990 and 2020, according to data from officially recognized sources compiled by the World Bank [19]. Hence, one gets the average vaccination rate p = 0.977 . The BCG has shown an overall efficacy of 70% to 80% against childhood TB, namely meningitis and miliary TB [20]. Hence, in this paper one takes the average rate of moving from V to S as the BCG immunization failure k = 1 0.75 = 0.25 . The treatment success from 2000 to 2020 [21] is used to calculate the treatment failure rate α = 1 0.8905 = 0.1095 .
The initial conditions were carefully selected as follows: The total initial population, N ( 0 ) , was set to 25,518,074, which corresponds to the population of Algeria in 1990, as reported in [18]. The initial infected population, I ( 0 ) = 11,607, was obtained from the WHO Global TB Report [1]. The initial exposed population, E ( 0 ) , was assumed to be 8852. Additionally, the initial treated population, T ( 0 ) , was set to 20,000, whereas the number of vaccinated individuals, V ( 0 ) , was determined to be 8,109,389. As a result of these values, the initial susceptible population can be calculated as S ( 0 ) = N ( 0 ) E ( 0 ) I ( 0 ) T ( 0 ) V ( 0 ) = 17,368,226. These initial conditions were carefully chosen to ensure accurate and reliable numerical simulations of the system under study.
One estimates the parameters β , γ , ϵ , σ , α , δ , k, and η by minimizing the error between actual TB incidence data and the solution of the proposed model (1). The objective function used in this parameter estimation is given by
ψ = i = 1 n ( I t i I t i * ) 2 ,
where I t i * denotes the actual TB-infected case, I t i is the corresponding model solution at time t i , and n is the number of available actual data. The MATLAB function ’fitnlm’, which solves nonlinear regression problems based on the Levenberg–Marquardt algorithm in MATLAB R2020b, was employed to minimize the function (17).
Table 3. Parameters and initial data of the model (1).
Table 3. Parameters and initial data of the model (1).
ParametersDescriptionAlgeria’s ParametersReferences
V ( 0 ) Initial number of vaccinated8,109,389Assumed
S ( 0 ) Initial number of susceptible17,368,226Calculated
E ( 0 ) Initial number of exposed8852Assumed
I ( 0 ) Initial number of infected11,607[1]
T ( 0 ) Initial number of treated20,000Assumed
Λ Recruitment rate811,085[18]
μ Natural death rate0.00498[18]
kRate of moving from V to S0.25[19]
β Transmission rate 6.6752 × 10 11 Fitted
γ Treatment rate0.0043Fitted
ϵ Progression rate0.0656Fitted
α Treatment failure rate0.1095[21]
δ Rate at which the treated0.1325Fitted
population leaves the class T
σ Disease death rate in I0.0136Fitted
η Disease death rate in T 4.2327 × 10 6 Fitted
pVaccination rate0.977[19]
In Figure 4 and Table 4, the incidence data are shown along with the model-fitted values, obtained using the values in Table 3.

3.2. Sensitivity Analysis

By using sensitivity, the spread and prevalence of diseases can be analyzed for each parameter. As a result of errors in data collection and assumed parameters, it is commonly used to measure the robustness of model predictions. To determine the relative significance of these parameters on disease transmission, we examined the impact of various model parameters. It has been determined how the model parameters β , γ , k, and ϵ affect the partial derivatives of the basic reproduction number R 0 . Since R 0 β > 0 , it follows that the transmission rate can be lowered to lessen the infection. However, given the partial derivatives R 0 γ < 0 , it is implied that TB infection can be controlled by raising the parameter γ . Increasing k results in an increase in R 0 because R 0 k > 0 . It demonstrates that the diseased population will grow more quickly. Parameter α is the failure rate of treatment and R 0 α > 0 . Therefore, by lowering the treatment failure rate α , the cumulative number of infected people can be minimized.
Sensitivity indices should be computed to estimate the relative change in a variable when parameters change. These indications were calculated using the following definition.
 Definition 1.
For a certain value σ, the normalized forward sensitivity index of R 0 is determined by
S σ R 0 = σ R 0 R 0 σ
For the baseline model parameters are determined by Formula (18), the derived sensitivity indices of the basic reproduction number R 0 are shown in Table 5.
We observe from Table 5 that the values of S β R 0 and S Λ R 0 are exactly + 1 . This indicates that a rise in β , and Λ will result in an increase in R 0 that is proportionate to both parameters. Additionally, we demonstrate that the parameters k, ϵ , α , and δ are exactly proportional to R 0 because S k R 0 > 0 , S ϵ R 0 > 0 , S α R 0 > 0 , and S δ R 0 > 0 . Moreover, the terms S μ R 0 < 0 , S γ R 0 < 0 , S σ R 0 < 0 , S p R 0 < 0 , and S η R 0 < 0 denote that the parameters μ , γ , σ , p, and η are inversely proportional to R 0 .

4. Results and Discussion

The results of parameters estimation are reported in Table 3, and Figure 4 illustrates the incidence data along with the model-fitted curve, obtained using the values in Table 3. The goodness fit of our model is supported by a height value coefficient of determination, namely R 2 = 0.9016 ; this indicates that the model fits the reported real data well.
Using the estimate parameters value, one obtains R 0 = 0.5228 , which is less than 1. This suggests that there is a possibility of decreasing or eliminating the disease by maintaining effective treatment and isolation measures in the future as illustrated by the model fitted curve for the period times 2020–2050 in Figure 4.
On the contrary, assuming that the government stops enforcing strict health measures against tuberculosis for children, such as vaccination, effective treatment strategies, and isolation of infected individuals after the year 2020, let us consider a hypothetical scenario. For instance, let us assume that p = 10 2 , γ = 10 3 , and β = 4 × 10 10 . The other parameters are taken from Table 3. With this parameter set, the basic reproduction number is found to be greater than one R 0 = 3.248 > 1 , indicating that the non-endemic equilibrium E 1 is unstable, whereas the endemic equilibrium E 2 is asymptotically stable. Obviously, the solutions of model (1) converge to E 2 as shown in Figure 5.
In order to gain a deeper understanding of how certain parameters affect the spread of the disease, one plots R 0 versus six parameters as shown in Figure 6. Obviously, there is a proportional relationship between the basic reproduction number R 0 and the three parameters β , ϵ , and α . This suggests that an increase in any of these parameters will result in an increase in the basic reproduction number, which in turn will lead to a greater spread of the disease.
On the other hand, this study found that there is an inverse relationship between the basic reproduction number R 0 and the other three parameters γ , p and μ . This implies that an increase in any of these parameters will result in a decrease in the basic reproduction number, which will lead to a slower spread of the disease. It is important to note that these results are in good agreement with reality.
The findings indicate that any plan aimed to prevent the spread of TB must consider four factors:
  • Enhancing the precision and quality of TB diagnosis to facilitate appropriate actions towards affected individuals.
  • Imposing isolation measures on infected individuals and monitoring their families medically to minimize contact with contagious patients.
  • Sustaining a high rate of vaccination for children to provide immunity.
  • Increasing the treatment rate by training specialized doctors, acquiring the most potent medicines, and establishing dedicated facilities to combat this disease.

5. Conclusions

In this study, we have developed a mathematical VSEIT model, which takes into account the biological factors of TB and certain realistic assumptions to analyze the transmission dynamics of this disease in Algeria. By using the reported infection data, we have estimated the model parameters using the least squares method. Our research has revealed that controlling the spread of TB is heavily dependent on some key factors, especially the contact parameter, β , the treatment parameter γ , and the vaccination parameter p. By identifying these crucial elements, we can better understand how to prevent and treat TB in Algeria. Using the estimated model parameters for Algeria we found that the reproduction number of the disease is less than one, meaning that TB can be eradicated by maintaining effective vaccination, high-quality treatment, and isolation measures. This research will have important implications for public health policymakers and healthcare professionals who are working to combat the spread of TB in Algeria and other regions where the disease is prevalent. The insights gained from this study can be used to develop more effective prevention and treatment strategies, which can ultimately help to reduce the burden of TB on affected communities.

Author Contributions

Conceptualization, M.S.A.; Methodology, M.S.A. and R.L.; Software, B.C.; Validation, R.L.; Investigation, B.C. and M.S.A.; Data curation, B.C.; Writing—original draft, B.C. and M.S.A.; Writing—review & editing, M.S.A. and R.L.; Supervision, M.S.A. and R.L.; Project administration, M.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Algerian Directorate General for Scientific Research and Technological Development OF FUNDER grant number C00L03CU430120230003.

Data Availability Statement

The data presented in this study are available in [1,18,19,21].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Basic Reproduction Number R 0

The basic reproduction number R 0 can be obtained using the next generation method [22], as R 0 = ρ ( F V 1 ) .
The associated matrices F for new infection in the infected compartments and V for the remaining transfer terms are given respectively by
F = β S I 0 0 , V = ( ϵ + μ ) E ( 1 α ) δ T ϵ E α δ T + ( γ + μ + σ ) I γ I + ( μ + δ + η ) T .
Let k 1 = ( ϵ + μ ) , k 2 = ( γ + μ + σ ) , k 3 = ( μ + δ + η ) , then
V = k 1 E ( 1 α ) δ T ϵ E α δ T + k 2 I γ I + k 3 T .
Next, we evaluate F and V which are the Jacobian of F and V respectively at E 1 such that F is non-negative and V is a non-singular matrix.
We denote the Jacobian of F and V by J ( F ) and J ( V ) respectively. Thus,
J ( F ) = 0 β S 0 0 0 0 0 0 0 , J ( V ) = k 1 0 ( 1 α ) δ ϵ k 2 α δ 0 γ k 3
F = J ( F ) and V = J ( V ) at E 1 , Thus,
F = 0 ( k + μ μ p ) Λ β μ ( k + μ ) 0 0 0 0 0 0 0 , and V = k 1 0 ( 1 α ) δ ϵ k 2 α δ 0 γ k 3
It follows that
V 1 = 1 ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) k 2 k 3 α δ γ γ ( 1 α ) k 2 ( 1 α ) ϵ k 3 k 1 k 2 k 1 α δ ( 1 α ) ϵ ϵ γ γ k 1 k 1 k 2 ,
then
F V 1 = 1 ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) ϵ ( k + μ μ p ) Λ β μ ( k + μ ) k 3 ϵ ( k + μ μ p ) Λ β μ ( k + μ ) k 1 k 2 ϵ ( k + μ μ p ) Λ β ( k 1 α δ ϵ ( 1 α ) ) μ ( k + μ ) 0 0 0 0 0 0
Thus,
R 0 = ρ ( F V 1 ) = ϵ ( k + μ μ p ) Λ β k 3 μ ( k + μ ) ( k 1 k 2 k 3 α γ δ k 1 ( 1 α ) δ γ ϵ ) ,

Appendix B. (LaSalle’s Invariance Principle)

Instead of solely focusing on the stability of a specific equilibrium point, this principle offers insights into the overall behavior and trajectories of the system’s solutions.
  • Consider the autonomous system
    x ˙ = f ( x ) , x R n ,
    where f is of class C 1 . The LaSalle’s invariance principle [23] is stated in the following theorem.
 Theorem A1.
Let Ω R n be a bounded closed (compact) set with the property that every solution of (A1) which begins in Ω remains for all future time in Ω.
Suppose there is a scalar function V ( x ) which has continuous first partials in Ω and is such that V ˙ ( x ) 0 in Ω. Let E be the set of all points in Ω where V ˙ ( x ) = 0 .
Let M be the largest invariant set in E. Then, every solution starting in Ω approaches M as t + .

References

  1. WHO. Global Tuberculosis Report; World Health Organization: Geneva, Switzerland, 2023; Available online: https://extranet.who.int/tme/generateCSV.asp?ds=notifications (accessed on 1 January 2023).
  2. Bernoulli, D. Essai d’une Nouvelle Analyse de la Mortalité Causée par la Petite Vérole, et des Avantages de l’inoculation Pour la Prévenir. Histoire de l’Acad. Roy. Sci. (Paris) avec Mem. 1760, pp. 1–45. Available online: https://gallica.bnf.fr/ark:/12148/bpt6k3558n/f223.item.r=daniel%20bernoulli (accessed on 1 January 2023).
  3. Ross, R. The Prevention of Malaria. Nature 1910, 85, 263–264. [Google Scholar] [CrossRef]
  4. Hamer, W.H. The Milroy Lectures on Epidemic Diseases in England: The Evidence of Variability and of Persistency of Type; The Bedford Press: London, UK, 1906; pp. 733–739. [Google Scholar]
  5. Martini, E. Berechnungen und Beobachtungen zur Epidemiologie und Bekampfung der Malaria; W. Gente: Hamburg, Germany, 1921. [Google Scholar]
  6. Lotka, A.J. Elements of Physical Biology; Williams and Wilkens: Philadelphia, PA, USA, 1925. [Google Scholar]
  7. Yang, Y.; Tang, S.; Ren, X.; Zhao, H.; Guo, C. Global stability and optimal control for a tuberculosis model with vaccination and treatment. Discret. Contin. Dyn. Syst. Ser. B 2016, 21, 1009–1022. [Google Scholar] [CrossRef] [Green Version]
  8. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character 1927, 115, 700–721. [Google Scholar]
  9. Waaler, H.; Geser, A.; Andersen, S. The use of mathematical models in the study of the epidemiology of tuberculosis. Am. J. Public Health Nation’s Health 1962, 52, 1002–1013. [Google Scholar] [CrossRef] [PubMed]
  10. Ullah, S.; Khan, M.A.; Farooq, M.; Gul, T. Modeling and analysis of Tuberculosis (TB) in Khyber Pakhtunkhwa, Pakistan. Math. Comput. Simul. 2019, 165, 182–184. [Google Scholar] [CrossRef]
  11. Abdelouahab, M.S.; Arama, A.; Lozi, R. Bifurcation analysis of a model of tuberculosis epidemic with treatment of wider population suggesting a possible role in the seasonality of this disease. Chaos Interdiscip. J. Nonlinear Sci. 2021, 31, 123125. [Google Scholar] [CrossRef] [PubMed]
  12. Andersen, P.; Doherty, T.M. The success and failure of BCG-implications for novel tuberculosis vaccine. Nature 2005, 3, 656–662. [Google Scholar] [CrossRef] [PubMed]
  13. Ucakan, Y.; Gulen, S.; Koklu, K. Analysing of Tuberculosis in Turkey through SIR, SEIR and BSEIR Mathematical Models. Math. Comput. Model. Dyn. Syst. 2021, 27, 179–202. [Google Scholar] [CrossRef]
  14. Egonmwan, A.O.; Okuonghae, D. Mathematical analysis of a tuberculosis model with imperfect vaccine. Int. J. Biomath. 2019, 12, 1950073. [Google Scholar] [CrossRef]
  15. Revelle, C.S.; Lynn, W.R.; Feldmann, F. Mathematical models for the economic allocation of tuberculosis control activities in developing nations. Am. Rev. Respir. Dis. 1967, 96, 893–909. [Google Scholar]
  16. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, Y.; Cao, J. Global stability of general cholera models with nonlinear incidence and removal rates. J. Frankl. Inst. 2015, 352, 2464–2485. [Google Scholar] [CrossRef]
  18. Population Growth in Algeria. Available online: https://www.donneesmondiales.com/afrique/algerie/croissance-population.php (accessed on 10 January 2022).
  19. Trading Economics. Immunization, BCG (% of One-Year-Old Children) Algeria. Available online: https://www.tradingeconomics.com/algeria/immunization-bcg-percent-of-one-year-old-children-wb-data.html/ (accessed on 20 January 2022).
  20. Katelaris, A.L.; Jackson, C.; Southern, J.; Gupta, R.K.; Drobniewski, F.; Lalvani, A.; Lipman, M.; Mangtani, P.; Abubakar, I. Effectiveness of BCG Vaccination Against Mycobacterium tuberculosis Infection in Adults: A Cross-sectional Analysis of a UK-Based Cohort. J. Infect. Dis. 2020, 221, 146–155. [Google Scholar] [CrossRef]
  21. The World Bank Group. Tuberculosis Treatment Success Rate (% of New Cases)—Algeria. Available online: https://data.worldbank.org/indicator/SH.TBS.CURE.ZS?locations=DZ (accessed on 5 January 2022).
  22. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
  23. LaSalle, J.P. Some extensions of Liapunov’s second method. IRE Trans. Circuit Theory 1960, 7, 520–527. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of VSEIT model.
Figure 1. Flowchart of VSEIT model.
Computation 11 00146 g001
Figure 2. Transcritical bifurcation of the model (1).
Figure 2. Transcritical bifurcation of the model (1).
Computation 11 00146 g002
Figure 3. Percentage of one-year-old children in Algeria receiving BCG vaccination during the time period 1990–2020.
Figure 3. Percentage of one-year-old children in Algeria receiving BCG vaccination during the time period 1990–2020.
Computation 11 00146 g003
Figure 4. Data fitting of the number of TB cases in Algeria.
Figure 4. Data fitting of the number of TB cases in Algeria.
Computation 11 00146 g004
Figure 5. The expected situation if the government decides to abandon strict health measuresafter 2020, where the used parameters are β = 4 × 10 10 , γ = 10 3 , σ = 0.136 , η = 4.2326 × 10 6 , δ = 0.1325 , α = 0.1095 , ϵ = 0.0656 , k = 0.25 , and p = 10 2 .
Figure 5. The expected situation if the government decides to abandon strict health measuresafter 2020, where the used parameters are β = 4 × 10 10 , γ = 10 3 , σ = 0.136 , η = 4.2326 × 10 6 , δ = 0.1325 , α = 0.1095 , ϵ = 0.0656 , k = 0.25 , and p = 10 2 .
Computation 11 00146 g005
Figure 6. Influence of the parameters β , γ , p , ϵ , α , and μ , respectively, on the basic reproduction number R 0 .
Figure 6. Influence of the parameters β , γ , p , ϵ , α , and μ , respectively, on the basic reproduction number R 0 .
Computation 11 00146 g006
Table 1. Description of variables of the model (1).
Table 1. Description of variables of the model (1).
VariableDescription
V ( t ) The vaccinated population at time t.
S ( t ) The susceptible population which is able to be infected at any time t.
E ( t ) The exposed population which is not yet infectious.
I ( t ) The infected population at time t.
T ( t ) The treated population at time t.
Table 2. Parameters of model (1).
Table 2. Parameters of model (1).
Model ParameterDescriptionUnit
Λ Recruitment rate y e a r 1
μ Natural death rate y e a r 1
kRate of moving from V to S y e a r 1
β Transmission rate y e a r 1
γ Treatment rate y e a r 1
ϵ Progression rate y e a r 1
α Treatment failure rate y e a r 1
δ Rate at which the treated population leave the class T y e a r 1
σ Disease death rate in I y e a r 1
η Disease death rate in T y e a r 1
pVaccination rate y e a r 1
Table 4. The reported data and the model fitted values of TB cases in Algeria.
Table 4. The reported data and the model fitted values of TB cases in Algeria.
YearReported DataNumerical ValueYearReported DataNumerical Value
199011,60711,607200621,14320,613
199111,33212,162200721,36920,884
199211,42812,793200820,58821,116
199313,34513,471200921,70121,313
199413,34514,137201022,33621,474
199513,50714,880201121,42921,604
199615,32915,578201221,88021,705
199716,52216,255201320,70121,778
199815,32416,255201422,51721,825
199916,64717,517201523,70521,850
200018,57218,090201622,80121,854
200118,25018,621201723,07721,838
200218,93419,108201823,46521,805
200319,73019,550201920,87921,757
200419,92919,947202017,21221,694
200521,33620,3012021-21,619
Table 5. Parameters and sensitivity index.
Table 5. Parameters and sensitivity index.
ParameterSensitivity Index
Λ + 1
μ 1.6502
k+0.0012
β + 1
γ 0.1671
ϵ + 0.0005
α + 2.1364 × 10 10
δ + 1.4003 × 10 09
σ 0.4043
η 2.5311 × 10 11
p 0.0194
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chennaf, B.; Abdelouahab, M.S.; Lozi, R. Analysis of the Dynamics of Tuberculosis in Algeria Using a Compartmental VSEIT Model with Evaluation of the Vaccination and Treatment Effects. Computation 2023, 11, 146. https://doi.org/10.3390/computation11070146

AMA Style

Chennaf B, Abdelouahab MS, Lozi R. Analysis of the Dynamics of Tuberculosis in Algeria Using a Compartmental VSEIT Model with Evaluation of the Vaccination and Treatment Effects. Computation. 2023; 11(7):146. https://doi.org/10.3390/computation11070146

Chicago/Turabian Style

Chennaf, Bouchra, Mohammed Salah Abdelouahab, and René Lozi. 2023. "Analysis of the Dynamics of Tuberculosis in Algeria Using a Compartmental VSEIT Model with Evaluation of the Vaccination and Treatment Effects" Computation 11, no. 7: 146. https://doi.org/10.3390/computation11070146

APA Style

Chennaf, B., Abdelouahab, M. S., & Lozi, R. (2023). Analysis of the Dynamics of Tuberculosis in Algeria Using a Compartmental VSEIT Model with Evaluation of the Vaccination and Treatment Effects. Computation, 11(7), 146. https://doi.org/10.3390/computation11070146

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop