# A Mathematical Model of COVID-19 with Vaccination and Treatment

Authors: M. L. Diagne; H. Rwezaura; S. Y. Tchoumi; J. M. Tchuenche
Journal: Computational and Mathematical Methods in Medicine
 (2021)

---

## Abstract

We formulate and theoretically analyze a mathematical model of COVID-19 transmission mechanism incorporating vital dynamics of the disease and two key therapeutic measures—vaccination of susceptible individuals and recovery/treatment of infected individuals. Both the disease-free and endemic equilibrium are globally asymptotically stable when the effective reproduction numberR0v is, respectively, less or greater than unity. The derived critical vaccination threshold is dependent on the vaccine efficacy for disease eradication whenever R0v>1, even if vaccine coverage is high. Pontryagin’s maximum principle is applied to establish the existence of the optimal control problem and to derive the necessary conditions to optimally mitigate the spread of the disease. The model is fitted with cumulative daily Senegal data, with a basic reproduction number R0=1.31 at the onset of the epidemic. Simulation results suggest that despite the effectiveness of COVID-19 vaccination and treatment to mitigate the spread of COVID-19, when R0v>1, additional efforts such as nonpharmaceutical public health interventions should continue to be implemented. Using partial rank correlation coefficients and Latin hypercube sampling, sensitivity analysis is carried out to determine the relative importance of model parameters to disease transmission. Results shown graphically could help to inform the process of prioritizing public health intervention measures to be implemented and which model parameter to focus on in order to mitigate the spread of the disease. The effective contact rate b, the vaccine efficacy ε, the vaccination rate v, the fraction of exposed individuals who develop symptoms, and, respectively, the exit rates from the exposed and the asymptomatic classes σ and ϕ are the most impactful parameters.


---

## Body

## 1. Introduction

The December 2019 outbreak of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing COVID-19, was first reported in Wuhan, Hubei Province of China [1–4]. Coronaviruses can be extremely contagious and spread easily from person to person [5]. The disease, now a global pandemic, has spread rapidly worldwide, causing major public health concerns and economic crisis [3, 4, 6], having a massive impact on populations and economies and thereby placing an extra burden on health systems around the planet [7–9]. In fact, all social levels of the society have suffered major disruptions due to the COVID-19 pandemic [10].Starting with the work of Daniel Bernoulli in 1760 [11], the development of mathematical models has been critical in our understanding of the dynamics of infectious diseases [12]. With the work of Kermack and McKendrick [13], mathematical models have since been used to provide framework for understanding the dynamics of infectious diseases. COVID-19 transmission dynamics models are flourishing and abound in the literature [14–19], to cite a few and the references therein. With the availability of COVID-19 vaccine and its known high efficacy, there is an urgent need to assess the impact of such vaccines with imperfect transmission-blocking effects [6] and potentially refine previous mathematical models of COVID-19 that incorporated the potential impact of an imperfect vaccine [2, 8, 20]. A recent study by Pearson et al. [21] found that COVID-19 vaccination in low- and middle-income settings is highly cost-effective and even cost saving, when the vaccine is reasonably priced and efficacy is high. Prior to pharmaceutical measures such as treatment and vaccination being available, nonpharmaceutical intervention measures such as self-quarantine of confirmed cases, isolation, face masks, hand washing, social/physical distancing, and the most restrictive lockdowns, closure, or limited openings of shops and schools have been relied upon and continue to be widely implemented [22–24].As COVID-19 vaccines are being deployed worldwide, we formulate and qualitatively analyze a COVID-19 mathematical model, taking into consideration available therapeutic measures, vaccination of susceptible and treatment of hospitalized/infected individuals. Our proposed model incorporates some key epidemiological and biological features of COVID-19, including demographic parameters (recruitment/birth and death). Optimal control is carried out using Pontryagin’s maximum principle as described in [25] and applied in epidemiological models [26–35]. To identify the model parameters with greater influence on the initial disease transmission R0v when vaccination and treatment are implemented [36], a sensitivity analysis is carried out using partial rank correlation coefficients (PRCCs) and the results are shown graphically. This identification is crucial to inform policy decision on which parameters to focus either for data collection or to mitigate the spread of the disease. To the best of our knowledge, this study provides the first in-depth mathematical analysis of the qualitative dynamics of COVID-19 with an imperfect vaccine and treatment.The rest of the paper is organized as follows. The proposed COVID-19 model is formulated in Section2 and theoretically analyzed in Section 3. By applying Pontryagin’s maximum principle, optimal control of the model to mitigate the spread of COVID-19 is presented in Section 4. Numerical simulations performed to support theoretical results are presented in Section 5. The conclusion is provided in Section 6.


## 2. Model Formulation and Analysis

Consider a homogeneous mixing within the population, i.e., individuals in the population have equal probability of contact with each other. Using a deterministic compartmental modeling approach to describe the disease transmission dynamics, at any timet, the total population Nt is subdivided into several epidemiological states depending on individuals’ health status: susceptible St, vaccinated Vt, exposed Et, symptomatic infected individuals It, infected asymptomatic At, hospitalized Ht, and recovered Rt. The total human population Nt is given by
(1)Nt=St+Vt+Et+It+At+Ht+Rt.Figure1 depicts the schematic model flow. The description of the model variables and parameters is presented in Table 1.Figure 1
Compartment diagram of the human component of the model.Table 1
Model parameter values and source.
ParameterDescriptionValueReferenceΠRecruitment rate of individuals into the population1000059×365[26, 38]pProportion of recruited individuals who are vaccinated0.0001AssumedvVaccination rate0.4AssumedωAReduction in the transmission from asymptomatic0.3[39]ωIIncrease in the transmission from symptomatic1.8AssumedωHReduction in the transmission from hospitalized0.3AssumedμNatural death rate159×365[26, 38]δDisease-induced death rate0.018[40, 41]εInfection reduction of vaccinated individuals0.8AssumedσExit rate from the exposed class0.13[9]γExit rate from the infectious class0.0833[41]κProportion of infectious who recover naturally0.05[42]ψFraction of exposed who become infected0.7[39]bEffective contact rate1.12[41]τRecovery rate of hospitalized individuals0.0701[40]ϕProportion of asymptomatic who recover naturally0.14[9]λExit rate from the asymptomatic class0.13978[39, 41]ηRate at which individuals lose immunity0.011[19]Since COVID-19 vaccination is available, it is realistic to consider a specific vaccinated classV. The transition rates from susceptible and vaccinated to exposed is given by, respectively, by the following force of infection βS=bωAA+ωII+ωHH/N and βV=1−εβS, where b is the effective contact rate, ωX with X∈A,I,H represent the transmission probability after a contact with an individual in status X, and ε represents the infection reduction of vaccinated individuals. The nonlinearity of the force of infection is one of the key features of dynamic infectious disease models, because to model a population of individuals, the status of each individual is required [37]. Individuals are recruited into the population at a rate Π, with a fraction p vaccinated and the remaining 1−p susceptible. The latter are vaccinated at a rate v. Because vaccine could be imperfect, that is, vaccination providing only partial protection, we assume that vaccinated individuals can become exposed to the disease b1−ε. It has been shown that even a partially efficacious vaccine can offer a fundamental solution to the SARS-CoV-2 pandemic [20]. The parameter ε represents the vaccine efficacy (or infection reduction of vaccinated individuals).After close contacts with symptomatic, asymptomatic, and hospitalized individuals, susceptible become exposed with the disease. We assume that the rate of disease transmission from asymptomatic to susceptible individuals is less than that from symptomatic and hospitalized individuals. While outbreaks usually persist for a shorter period of time, the COVID-19 pandemic which started in December 2019 is still ongoing, and for this reason, we incorporate vital dynamics (recruitment and death). Letσ be the exit rate from exposed class where a fraction ψ develops infection while the remaining 1−ψ becomes asymptomatic. The exit rate from the asymptomatic class is λ. Asymptomatic individuals A are diminished by natural death at a rate μ (it is assumed that death due to the disease in this group of individuals is negligible), by those developing symptoms and moving to the asymptomatic class at a rate 1−ϕ, while a fraction ϕ may recover naturally from the asymptomatic infection and move to the recovered class R. Exit from the infected class is γ, where a fraction 1−κ are hospitalized, and a fraction κ recovers naturally. Finally, hospitalized individuals are treated and recovered at a rate τ. It is assumed that symptomatic and hospitalized individuals experience and additional disease-induced death rate δ, respectively. We also consider that the recovered individuals die at a rate μ, while a fraction η becomes susceptible again.From the aforementioned and the model flow diagram of the disease transmission mechanisms Figure1, we derive the following nonlinear system of ordinary differential equations that captures the transmission dynamics of COVID-19.From the model flow diagram in Figure1, we derive the following system of nonlinear ordinary differential equations:
(2)Ṡ=1−pΠ+ηR−βS+μ+vS,V̇=pΠ+vS−βV+μV,Ė=βSS+βVV−σ+μE,İ=σψE+λ1−ϕA−γ+μ+δI,Ȧ=σ1−ψE−λ+μA,Ḣ=γ1−κI−τ+μ+δH,Ṙ=γκI+λϕA+τH−η+μR,with initial conditions
(3)S0≥0,V0≥0,E0≥0,I0≥0,A0≥0,H0≥0,R0≥0,where
(4)βS=bωAA+ωII+ωHHN,βV=b1−εωAA+ωII+ωHHN.For simplicity, letg1=μ+v, g2=μ, g3=σ+μ, g4=γ+μ+δ, g5=λ+μ, g6=τ+μ+δ, and g7=η+μ. Then, model system 1 now reads
(5)Ṡ=1−pΠ+ηR−βS+g1S,V̇=pΠ+vS−βV+g2V,Ė=βSS+βVV−g3E,İ=σψE+λ1−ϕA−g4I,Ȧ=σ1−ψE−g5A,Ḣ=γ1−κI−g6H,Ṙ=γκI+λϕA+τH−g7R.All the model parameters and their description, values, and sources are presented in Table1.


## 3. Model Analysis

Well-posedness, nonnegativity, and boundedness of solutions of the proposed model can be shown using basic theory of dynamical systems as described in [43, 44]; also see [45, 46]. By adding all the equations of the system, we have Ṅ=Π−μN−δI+H≤Π−μN. Therefore, it follows that the biologically feasible region for model 1 is
(6)D=S,V,E,I,A,H,R∈ℝ+7:N≤Πμ.


### 3.1. Disease-Free Equilibrium and Basic Reproduction Number

Model system 1 admits a disease-free equilibrium (DFE) given byE0=S0,V0,0,0,0,0,0, where
(7)S0=Π1−pg1,V0=Πpg1+v1−pg1g2.The linear stability ofE0 is established using the next-generation method [47, 48]. The rate of appearance of new infections and the rate of transfer of individuals by all other means are given by the following at least twice continuously differentiable functions
(8)F=bωAA+ωII+ωHHNS+b1−εωAA+ωII+ωHHNV000,V=−g3EσψE+λ1−ϕA−g4Iσ1−ψE−g5Aγ1−κI−g6H.From [48], the nonnegative matrix F and the nonsingular M-matrix V for the new infection terms and the remaining transfer terms are given by
(9)F=0bωIS0+1−εV0N0bωAS0+1−εV0N0bωHS0+1−εV0N0000000000000,V=−g3000σψ−g4λ1−ϕ0σ1−ψ0−g500γ1−κ0−g6.Thus, the effective reproduction number is given by(10)R0v=λ1−ϕ1−ψ+g5ψg6ωI+γ1−κωH+g4g61−ψωAg21−p+1−εpg1+1−pvbμσg1g2g3g4g5g6=bμσG1G2+G3G4+1−εpμ+v+1−pvμ+vg2g3g4g5g6,where G1=λ1−ϕ1−ψ+g5ψ, G2=g6ωI+γ1−κωH, G3=g4g61−ψωA, and G4=g21−p.The effective reproduction numberR0v is defined as the average number of secondary infections generated by a single infectious individual during his entire duration infectiousness in a totally susceptible population when vaccination is implemented.


### 3.2. Stability of the Disease-Free Equilibrium

We now study the global stability of the DFE using the approach described in [49]. Consider a system of ordinary differential equations of the form
(11)dxdt=Fx,I,dIdt=Gx,I,Gx,0=0,where x∈ℝm denotes (its components) the number of uninfected individuals and I∈ℝn denotes (its components) the number of infected individuals including latent and infectious. Let U0=x∗,0 be the disease-free equilibrium of this system, where 0 is a zero vector. Global stability of the DFE is guaranteed when the following conditions (H1) and (H2) are satisfied.(H1) Fordx/dt=Fx,0, 0 is globally asymptotically stable (g.a.s.).(H2)Gx,I=AI−G^x,I,G^x,I⩾0 for x,I∈Ω, where A=DIGx∗,0 is an M-matrix (the off diagonal elements of A are nonnegative) and Ω is the region where the model makes biological sense.Corollary 1 (see [49]).
The fixed pointU0=x∗,0 is a globally asymptotic stable (g.a.s.) equilibrium of (11) provided that R0v<1 (l.a.s.) and that assumptions (H1) and (H2) are satisfied.Theorem 2 (global asymptotic stability of the DFE).
The DFEE0 of model 1 is globally asymptotically stable if R0v<1.Proof.
First, we rewrite model 1 in the form 6 by settingx=S,V and I=E,I,A,H,R.
Then, the DFE is given byU0=x∗,0=Π1−p/g1,Πpg1+v1−p/g1g2,0 and the system dx/dt=Fx,0 becomes
(12)Ṡ=1−pΠ−g1S,V̇=pΠ+vS−g2V.
This equation has a unique equilibrium point(13)x∗=Π1−pg1,Πpg1+v1−pg1g2,which is globally asymptotically stable. Therefore, the condition (H1) is satisfied.
We now verify the second condition (H2). For model 1, we have(14)Gx,I=βSS+βVV−g3EσψE+λ1−ϕA−g4Iσ1−ψE−g5Aγ1−κI−g6HγκI+λϕA+τH−g7R,(15)DIGx∗,0=−g3bωIN0S0+V01−εbωAN0S0+V01−εbωHN0S0+V01−ε0σψ−g4λ1−Φ00σ1−ψ0−g5000γ1−κ0−g600γκλΦ0−g7.
Clearly,A is an M-matrix. On the other hand,
(16)G^x,I=AI−Gx,I=βSS0N0−S+βvV0N0−V0000,which implies that G^x,I⩾0 for all x,I∈Ω. Therefore, the conditions (H1) and (H2) are satisfied. By Corollary 1, the global stability of the DFE is obtained. This completes the proof.☐☐Global stability of the DFE precludes the model to exhibit bistability also known as backward bifurcation [50, 51], a situation where both the disease-free and endemic equilibria coexist when R0v<1.


### 3.3. The Critical Vaccination Coverage

We investigate the critical vaccination coverage rate that could help eradicate the disease.R0v≔RE. When there is no vaccination in the community, that is, v=0, then the effective reproduction number reduces to
(17)R0=R00=bσG1G2+G3G4+1−εpμg2g3g4g5g6.In fact,R0 is the so-called basic reproduction number which is the average number of secondary cases arising from one infectious individual in a totally susceptible population [47, 52]. After some rearrangement, R0v can be written as
(18)R0v=μμ+vG4+1−εpμ+vG4+1−εpμR0=μμ+vμ1−p+1−εpμ+vμ1−p+1−εpμR0=1μ+vμ1−εp+1−εv1−εpR0=R0μ+vμ+1−εv1−εp=R0μμ+v+1−ε1−εpvμ+v.Thus,(19)R0∞≔limv⟶+∞R0v=1−ε1−εpR0.Taking the partial derivative ofR0v with respect to v yields
(20)∂R0v∂v=μμ+v2εp−11−εpR0=−ε1−p1−εpμR0μ+v2<0.Therefore,1−ε/1−εpR0⩽R0v⩽R0, and hence, R0<1 implies R0v<1, but the reverse is not true. For R0>1, it is important to note that
(21)R0∞<1⇔1−ε1−εpR0<1⇔ε>ε∗≔R0−1R0−p.This can be interpreted to mean that when the vaccine efficacyε is low and R0 is far greater than unity, the disease may not be eradicated even if the vaccine coverage is high. Figure 2 can be interpreted as additional efforts will be needed to reduce R0v below unity even when there is a high vaccine coverage v. In fact, the role of human/social behavior among those vaccinated (behavior compensation) such as careless increase in social/physical contacts can undermine vaccine impact [18, 20]. Thus, cautious phased relaxation of nonpharmaceutical interventions could substantially reduce population-level morbidity and mortality [23].Figure 2
Graphical representation ofε∗R0.The next result provides the critical vaccination thresholdv∗ for disease eradication.Lemma 3.
Assume that the basic reproduction numberR0>1. Then, there exists
(22)v∗=μ1−εpR0−1ε−ε∗R0−p>0,such that R0v∗=1. Furthermore, R0v><1 when v<>v∗.Figure3 depicts the vaccine coverage v as a function of R0, with a 91% vaccine efficacy ε>ε∗. The green surface represents the case when the vaccine coverage v exceeds the critical vaccination threshold v∗.Figure 3
Graphical representation ofv∗R0.


### 3.4. Stability of the Endemic Equilibrium

After some algebraic manipulations, the endemic equilibrium of model system 1 is obtained as(23)S∗=1−pΠ+ηR∗βS∗+g1,V∗=pΠ+vS∗1−εβS∗+g2,E∗=βS∗S∗+1−εV∗g3,A∗=σ1−ψE∗g5,I∗=σψE∗+λ1−ϕA∗g4,H∗=γ1−κI∗g6.From the last equation of model system 1 and using the definition ofβS∗, we obtain the following quadratic equation
(24)βS∗P2βS∗2+P1βS∗+P0=0,where P1 and P2 are positive constants given in the appendix and
(25)P0=g1g2g3g4g5g6R0v−1.Thus,βS∗=0 implies I=0, which represents the disease-free equilibrium, and hence, because P1 and P2 are positive, equation (24) has a unique positive solution if and only if P0>0, which is fulfilled when R0v>1. Therefore, we have the following result.Lemma 4.
IfR0v>1, model system 1 has a unique endemic equilibrium E∗.It can also be shown using the theory or permanence as described in [53] that when R0v>1, the disease class It is uniformly persistent, that is, there exists a positive constant k, such that
(26)liminft⟶∞It≥k.Consequently, the model system is uniformly persistent whenR0v>1.Figure4 depicts the impact of the reduction in transmission from asymptomatic ωA, and the increase in transmission from symptomatic ωI on the effective reproduction number R0v>1 is shown in Figure 4. Similar graphical representation for different model parameters space can be found in [54]. The white panel in Figure 4 shows the endemic equilibrium regions in the (ωI,ωA) space when R0v>1. The solid line corresponds to R0v=1. To mitigate the spread of the disease, it is important that the reduction in the transmission from asymptomatic ωA<0.36 and the increase in transmission from symptomatic ωI<1.16.Figure 4
Endemic equilibrium regions whenR0v>1 in the (ωI,ωA) space.


## 3.1. Disease-Free Equilibrium and Basic Reproduction Number

Model system 1 admits a disease-free equilibrium (DFE) given byE0=S0,V0,0,0,0,0,0, where
(7)S0=Π1−pg1,V0=Πpg1+v1−pg1g2.The linear stability ofE0 is established using the next-generation method [47, 48]. The rate of appearance of new infections and the rate of transfer of individuals by all other means are given by the following at least twice continuously differentiable functions
(8)F=bωAA+ωII+ωHHNS+b1−εωAA+ωII+ωHHNV000,V=−g3EσψE+λ1−ϕA−g4Iσ1−ψE−g5Aγ1−κI−g6H.From [48], the nonnegative matrix F and the nonsingular M-matrix V for the new infection terms and the remaining transfer terms are given by
(9)F=0bωIS0+1−εV0N0bωAS0+1−εV0N0bωHS0+1−εV0N0000000000000,V=−g3000σψ−g4λ1−ϕ0σ1−ψ0−g500γ1−κ0−g6.Thus, the effective reproduction number is given by(10)R0v=λ1−ϕ1−ψ+g5ψg6ωI+γ1−κωH+g4g61−ψωAg21−p+1−εpg1+1−pvbμσg1g2g3g4g5g6=bμσG1G2+G3G4+1−εpμ+v+1−pvμ+vg2g3g4g5g6,where G1=λ1−ϕ1−ψ+g5ψ, G2=g6ωI+γ1−κωH, G3=g4g61−ψωA, and G4=g21−p.The effective reproduction numberR0v is defined as the average number of secondary infections generated by a single infectious individual during his entire duration infectiousness in a totally susceptible population when vaccination is implemented.


## 3.2. Stability of the Disease-Free Equilibrium

We now study the global stability of the DFE using the approach described in [49]. Consider a system of ordinary differential equations of the form
(11)dxdt=Fx,I,dIdt=Gx,I,Gx,0=0,where x∈ℝm denotes (its components) the number of uninfected individuals and I∈ℝn denotes (its components) the number of infected individuals including latent and infectious. Let U0=x∗,0 be the disease-free equilibrium of this system, where 0 is a zero vector. Global stability of the DFE is guaranteed when the following conditions (H1) and (H2) are satisfied.(H1) Fordx/dt=Fx,0, 0 is globally asymptotically stable (g.a.s.).(H2)Gx,I=AI−G^x,I,G^x,I⩾0 for x,I∈Ω, where A=DIGx∗,0 is an M-matrix (the off diagonal elements of A are nonnegative) and Ω is the region where the model makes biological sense.Corollary 1 (see [49]).
The fixed pointU0=x∗,0 is a globally asymptotic stable (g.a.s.) equilibrium of (11) provided that R0v<1 (l.a.s.) and that assumptions (H1) and (H2) are satisfied.Theorem 2 (global asymptotic stability of the DFE).
The DFEE0 of model 1 is globally asymptotically stable if R0v<1.Proof.
First, we rewrite model 1 in the form 6 by settingx=S,V and I=E,I,A,H,R.
Then, the DFE is given byU0=x∗,0=Π1−p/g1,Πpg1+v1−p/g1g2,0 and the system dx/dt=Fx,0 becomes
(12)Ṡ=1−pΠ−g1S,V̇=pΠ+vS−g2V.
This equation has a unique equilibrium point(13)x∗=Π1−pg1,Πpg1+v1−pg1g2,which is globally asymptotically stable. Therefore, the condition (H1) is satisfied.
We now verify the second condition (H2). For model 1, we have(14)Gx,I=βSS+βVV−g3EσψE+λ1−ϕA−g4Iσ1−ψE−g5Aγ1−κI−g6HγκI+λϕA+τH−g7R,(15)DIGx∗,0=−g3bωIN0S0+V01−εbωAN0S0+V01−εbωHN0S0+V01−ε0σψ−g4λ1−Φ00σ1−ψ0−g5000γ1−κ0−g600γκλΦ0−g7.
Clearly,A is an M-matrix. On the other hand,
(16)G^x,I=AI−Gx,I=βSS0N0−S+βvV0N0−V0000,which implies that G^x,I⩾0 for all x,I∈Ω. Therefore, the conditions (H1) and (H2) are satisfied. By Corollary 1, the global stability of the DFE is obtained. This completes the proof.☐☐Global stability of the DFE precludes the model to exhibit bistability also known as backward bifurcation [50, 51], a situation where both the disease-free and endemic equilibria coexist when R0v<1.


## 3.3. The Critical Vaccination Coverage

We investigate the critical vaccination coverage rate that could help eradicate the disease.R0v≔RE. When there is no vaccination in the community, that is, v=0, then the effective reproduction number reduces to
(17)R0=R00=bσG1G2+G3G4+1−εpμg2g3g4g5g6.In fact,R0 is the so-called basic reproduction number which is the average number of secondary cases arising from one infectious individual in a totally susceptible population [47, 52]. After some rearrangement, R0v can be written as
(18)R0v=μμ+vG4+1−εpμ+vG4+1−εpμR0=μμ+vμ1−p+1−εpμ+vμ1−p+1−εpμR0=1μ+vμ1−εp+1−εv1−εpR0=R0μ+vμ+1−εv1−εp=R0μμ+v+1−ε1−εpvμ+v.Thus,(19)R0∞≔limv⟶+∞R0v=1−ε1−εpR0.Taking the partial derivative ofR0v with respect to v yields
(20)∂R0v∂v=μμ+v2εp−11−εpR0=−ε1−p1−εpμR0μ+v2<0.Therefore,1−ε/1−εpR0⩽R0v⩽R0, and hence, R0<1 implies R0v<1, but the reverse is not true. For R0>1, it is important to note that
(21)R0∞<1⇔1−ε1−εpR0<1⇔ε>ε∗≔R0−1R0−p.This can be interpreted to mean that when the vaccine efficacyε is low and R0 is far greater than unity, the disease may not be eradicated even if the vaccine coverage is high. Figure 2 can be interpreted as additional efforts will be needed to reduce R0v below unity even when there is a high vaccine coverage v. In fact, the role of human/social behavior among those vaccinated (behavior compensation) such as careless increase in social/physical contacts can undermine vaccine impact [18, 20]. Thus, cautious phased relaxation of nonpharmaceutical interventions could substantially reduce population-level morbidity and mortality [23].Figure 2
Graphical representation ofε∗R0.The next result provides the critical vaccination thresholdv∗ for disease eradication.Lemma 3.
Assume that the basic reproduction numberR0>1. Then, there exists
(22)v∗=μ1−εpR0−1ε−ε∗R0−p>0,such that R0v∗=1. Furthermore, R0v><1 when v<>v∗.Figure3 depicts the vaccine coverage v as a function of R0, with a 91% vaccine efficacy ε>ε∗. The green surface represents the case when the vaccine coverage v exceeds the critical vaccination threshold v∗.Figure 3
Graphical representation ofv∗R0.


## 3.4. Stability of the Endemic Equilibrium

After some algebraic manipulations, the endemic equilibrium of model system 1 is obtained as(23)S∗=1−pΠ+ηR∗βS∗+g1,V∗=pΠ+vS∗1−εβS∗+g2,E∗=βS∗S∗+1−εV∗g3,A∗=σ1−ψE∗g5,I∗=σψE∗+λ1−ϕA∗g4,H∗=γ1−κI∗g6.From the last equation of model system 1 and using the definition ofβS∗, we obtain the following quadratic equation
(24)βS∗P2βS∗2+P1βS∗+P0=0,where P1 and P2 are positive constants given in the appendix and
(25)P0=g1g2g3g4g5g6R0v−1.Thus,βS∗=0 implies I=0, which represents the disease-free equilibrium, and hence, because P1 and P2 are positive, equation (24) has a unique positive solution if and only if P0>0, which is fulfilled when R0v>1. Therefore, we have the following result.Lemma 4.
IfR0v>1, model system 1 has a unique endemic equilibrium E∗.It can also be shown using the theory or permanence as described in [53] that when R0v>1, the disease class It is uniformly persistent, that is, there exists a positive constant k, such that
(26)liminft⟶∞It≥k.Consequently, the model system is uniformly persistent whenR0v>1.Figure4 depicts the impact of the reduction in transmission from asymptomatic ωA, and the increase in transmission from symptomatic ωI on the effective reproduction number R0v>1 is shown in Figure 4. Similar graphical representation for different model parameters space can be found in [54]. The white panel in Figure 4 shows the endemic equilibrium regions in the (ωI,ωA) space when R0v>1. The solid line corresponds to R0v=1. To mitigate the spread of the disease, it is important that the reduction in the transmission from asymptomatic ωA<0.36 and the increase in transmission from symptomatic ωI<1.16.Figure 4
Endemic equilibrium regions whenR0v>1 in the (ωI,ωA) space.


## 4. Optimal Control Problem

We investigate the impact of implementing pharmaceutical interventions to mitigate the spread of COVID-19. To accomplish this, we introduce a set of time-dependent control variablesu1t,u2t where
(a)
u1t represents the implementation of continuous vaccination(b)
u2t represents treatment of infected (often hospitalized) individualsThe proposed COVID-19 model with optimal controlu1t,u2t consists of the following nonautonomous system of nonlinear ordinary differential equations.
(27)Ṡ=1−pΠ+ηR−βS+μ+u1S,V̇=pΠ+u1S−βV+μ+uV,Ė=βSS+βVV−σ+μE,İ=σψE+λ1−ϕA−γ+μ+δI,Ȧ=σ1−ψE−λ+μA,Ḣ=γ1−κI−u2+μ+δH,Ṙ=γκI+λϕA+u2H−η+μR.We wish to find the controls that minimize the total infected individuals, that is, to find an optimal control for the two control strategies while reducing their relative costs. In other words, we want to find the optimal values ofu1t,u2t that minimize the objective functional Ju1,u2 where
(28)Ju1,u2=∫0TA1E+A2I+A3A+A4H+B1u12t+B2u22tdt,subject to the differential equation (27), where T is the final time. This objective functional involves the total exposed, asymptomatic, infected, and hospitalized individuals, along with the cost of applying the controls u1t and u2t. We consider a quadratic objective functional for measuring the control cost as frequently used in the literature [29–32, 55]. The positive coefficients A1,A2,A3,A4,B1 and B2 are balancing weight parameters, while the controls u1tandu2t are bounded, Lebesgue integrable functions [31, 32]. We then seek to find optimal controls u1∗ and u2∗, such that
(29)Ju1∗,u2∗=minΩJu1,u2.To derive the necessary conditions that the two optimal controls and corresponding states must satisfy, we apply Theorem 5.1 (Pontryagin’s maximum principle [25]) in Fleming and Rishel [56] to develop the optimal system for which the necessary conditions that must be satisfied by an optimal control and its corresponding states are derived. In fact, Theorem 1 in Agusto [26] which is based on the boundedness of solution of model system 1 without control variables ensures the existence of the optimal control while the existence of an optimal control with a given control pair follows from Fleming and Rishel [56] and Caratheodory’s existence theorem [57]. For discussions on various forms of the objective functional (linear, quadratic), see [30, 58].
(30)ℍ=A1E+A2I+A3A+A4H+B1u12+B2u22+ξ11−pΠ+ηR−βS+μ+u1S+ξ2pΠ+u1S−βV+μ+uV+ξ3βSS+βVV−σ+μE+ξ4σψE+λ1−ϕA−γ+μ+δI+ξ5σ1−ψE−λ+μA+ξ6γ1−κI−u2+μ+δH+ξ7γκI+λϕA+u2H−η+μR,where ξi,i=1,⋯,7 are the adjoint variables or costate variables. The following result presents the adjoint system and control characterization.Theorem 5.
Given an optimal controlu1∗,u2∗ and corresponding state solutions S,V,E,I,A,H,R of the corresponding state system 1, there exists adjoint variables, ξi,i=1,⋯,7, satisfying
(31)ξ′1=1−SNβsξ1−ξ3+VNβvξ3−ξ2+u1ξ1−ξ2+μξ1,ξ′2=βsNSξ3−ξ1+1−VNβvξ3−ξ2+uξ2−ξ1+μξ2,ξ′3=βsNSξ3−ξ1+βvNVξ3−ξ2+σϕξ5−ξ4+σξ3−ξ5+μξ3−A1,ξ′4=SNβs−bωIξ3−ξ1+VNβv−b1−εωIξ3−ξ2+γκξ6−ξ7+γξ4−ξ6+μ+δξ4−A2,ξ′5=SNβs−bωAξ3−ξ1VNβv−b1−εωAξ3−ξ2+Φλξ4−ξ7+λξ5−ξ4+μξ5−A3,ξ′6=SNbωH−βsξ1−ξ3+VNb1−ε−βvξ2−ξ3+u2ξ6−ξ7+μ+δξ6−A4,ξ′7=SNξ3−ξ1+VNβvξ3−ξ2+ηξ7−ξ1+μξ7,ξ′iT=0fori=1,⋯7.
The controlsu1∗ and u2∗ satisfy the following optimality condition:
(32)u1∗=max0,min1,ξ2−ξ1S2B1,u2∗=max0,min1,ξ7−ξ6H2B2,Proof.
The differential equations governing the adjoint variables are obtained by differentiation of the Hamiltonian function, evaluated at the optimal control. Then, the adjoint system can be written as(33)ξ′1=−∂ℍ∂S,ξ′2=−∂ℍ∂V,ξ′3=−∂ℍ∂E,ξ′4=−∂ℍ∂I,ξ5′=−∂ℍ∂A,ξ′6=−∂ℍ∂H,ξ′7=−∂ℍ∂R,with zero final time conditions (transversality) ξiT=0. Replacing the derivatives of ℍ with respect to S,V,E,I,A,H,R in the above equations, we obtain the optimality condition (32). The optimal conditions for the Hamiltonian are given by ∂ℍ/∂u1∗=∂ℍ/∂u2∗=0 or equivalently ∂ℍ/∂u1∗=2B1u1∗+ξ2−ξ1S=0,∂ℍ/∂u2∗=2B2u2∗+ξ7−ξ6H=0.From the above equations, we obtain
(34)u1∗=ξ1−ξ2S2B1,u2∗=ξ6−ξ7H2B2.
Thus,u1∗ and u2∗ satisfy (32).☐☐


## 5. Numerical Simulations

To illustrate the theoretical results, numerical simulations are carried out. Model parameter values for the numerical simulations with their description and source are listed in Table1. Whenever parameter values were not available in the literature, we assumed realistic values for the purpose of illustration.


### 5.1. Model Fitting

Following the approach described in [14], model system 1 is fitted with daily new COVID-19 cases in Senegal from 29 March to 29 April 2020, which correspond to the first wave of infections in the country. It is important to note that while starting with an expanded model helps to account for the various disease classes explicitly, for the model fitting, model system 1 is reduced to a simpler version which retains the key compartments and characteristics of the complex one by eliminating similar classes [28]. This reduction involves removing the asymptomatic and hospitalized classes, which implies that the reduction in the transmission from hospitalized individuals is equated to zero, that is, ωH=0. The main key requirement for the simple model to approximate the complex one is that the susceptible class S should contain the same number of individuals in both models (when one compares both models compartment wise). The observed data were fitted in MatLab using the optimization function createOptimProblem. The estimated value of the basic reproduction number of the disease in Senegal during the period under consideration is R0=1.31, a value close to unity, which might explain why the first wave of COVID-19 infections did not pick up in Senegal. Figure 5 depicts the daily cumulative number of COVID-19 cases in Senegal. Red dots represent actual confirmed cases while the blue curve is the best fitting curve of the model.Figure 5
Model fit with cumulative daily COVID-19 cases in Senegal, 29 March–29 April 2021.


### 5.2. Long-Term Dynamics of the Disease

We investigated the impact of vaccination and treatment on mitigating the spread of COVID-19. An iterative fourth-order Runge-Kutta method (both forward and backward algorithms) is employed to compute the optimal controls and state values used. For more details on this approach, see [29, 32].The baseline weight parametersA1=A2=A3=A4=1,B1=12,B2=19 are chosen to illustrate the optimal control strategies, as well as the following nonnegative initial conditions S=2500,V=10,E=20,I=70,A=3,H=3,R=1. These weights do not necessarily have a significant meaning attached, but are only of theoretical sense to illustrate our proposed control strategies [31–33, 59]. Using parameter values in Table 1, the reproduction number R0v>1, indicating that the disease is endemic in the population. The positive constants A1, A2, A3, and A4 represent, respectively, the weight that balance off the COVID-19 exposed, infected, asymptomatic, and hospitalized individuals; B1 and B2 are, respectively, the weight constant for vaccination and treatment. Because low cost could potentially be associated with COVID-19 vaccination compared to its treatment (as the cost associated to u2t includes the cost of medical examination, hospitalization, and medicines), the weight factor B1 has been made lower than B2, while the two control strategies u1,u2 are all constrained between zero and one, that is, 0≤uit≤1,i=1,2. For instance, if u1=0, it implies no COVID-19 vaccination, and u1≠0 implies vaccination campaign measures are being implemented in the community. Table 1 provides all the model parameter values used for the simulations.Figure6 depicts the time series of model system 1 for the susceptible, vaccinated, exposed, and infected classes. Because of the potential limitations of vaccines in some settings, we consider in Figure 7 a saturated Holling type II vaccination rate
(35)v≔v1+ξS,where ξ represents the limitation of the availability of vaccine [60]. We observe that when vaccines are widely available, the number of susceptible individuals decreases pretty fast as more and more people are getting the vaccine (the vaccinated class V increases at the onset of implementation (Figure 7)).Figure 6
Profile ofS, V, E, and I without saturation.Figure 7
Profile ofS, V, E, and I with saturation.


### 5.3. Impact of Control Interventions: Vaccination and Treatment

In the next set of figures generated from model system 14, optimal control strategies are implemented. Using the model parameter values in Table1, the basic reproduction number R0=2.21>1. For community with no vaccination program v=0, the basic reproduction number R0=4.67 is almost double compared to the case when vaccination program is introduced. These values are in agreement with those estimated in recent COVID-19 modeling studies [61–63]. Figures 8–12 depict the graphical representations of the simulations of the COVID-19 model as a function of time without control and with optimal control. As can be observed on Figures 8–12, implementing control measures at the optimum level could help to mitigate the spread of the disease as shown by the significant decrease in the total number of individuals in the diseased classes and an increased in the vaccinated class. Figure 13 depicts the profile for the effect of the control functions u1t and u2t on the dynamics of the model system. The two control measures need to be optimally implemented for the first 75 days and maintained for at least 3 months (close to 100 days). It is important however to note that while therapeutic measures, vaccination, and treatment are very effective in curtailing the spread of the epidemic, more control efforts are required to eradicate the disease when R0v>1. Thus, continuously and concurrently applying both pharmaceutical and nonpharmaceutical public health interventions such as face mask, hand washing, and social distancing should be encouraged [64].Figure 8
Dynamics of the vaccinated classV.Figure 9
Dynamics of the exposed classE.Figure 10
Dynamics of the asymptomatic classA.Figure 11
Dynamics of the infected classI.Figure 12
Dynamics of the hospitalized classH.Figure 13
Controlsu1t and u2t.


### 5.4. Sensitivity of the Reproduction NumberR0v

Because mathematical models are symbolic/mechanistic representations of complex biological systems, some model parameter values are not often known with certainty due to natural and seasonal variations, potential measurement errors [18]. To show how changes in model parameters values affect R0v, we determine the relative importance of model parameters to disease transmission and graphically depict how sensitive the effective reproduction number is to the model parameters. Early ranking the intervention measures and other model parameters based on their impact could ideally partially inform the process of prioritizing public health intervention measures to be implemented, thereby helping policy and decision-makers to focus on those key impactful interventions. Figure 14 indicates that when the vaccine efficacy ε is low and R0v is greater than unity, the disease may not be eradicated even if the vaccine coverage v is high. The epidemiological implication is that high vaccine efficacy and vaccination coverage will drastically reduce the number of secondary infections in the community. This result also agrees with the conclusion from Figure 2. The graphical representation (contour plots) of R0v in the parameter space (v,τ) is similar to Figure 14, and this tends to suggest that the impact of the vaccine efficacy is similar to the effect of treatment on the initial disease prevalence.Figure 14
Contour plot ofR0v with respect to vaccine coverage v and efficacy ε.Partial rank correlation describes the relationship between two variables while at the same time removing the effects of several other variables from the relationship [65, 66]. We perform sensitivity analysis by employing the partial rank correlation coefficients (PRCCs) and the Latin hypercube scheme to identify the impact of each model parameters on the initial disease transmission R0v. PRCCs showing the effect of varying the input parameters on the effective reproduction number R0v are shown in Figure 15. All parameters with positive PRCCs will result in an increase on the number of initial disease transmission, while an increase in parameters with negative PRCCs will result in a reduction of R0v. By exploring Figure 15, the most influential parameters can be identified. Model parameters that should be targeted to reduce the spread of the disease are the effective contact rate b, the infection reduction (vaccine efficacy) of vaccinated individuals ε, the fraction of exposed individuals who develop symptoms, and, respectively, the exit rates from the exposed and the asymptomatic classes σ and ϕ.Figure 15
PRCCs showing the effect of varying the input parameters onR0v.Figures16 depicts the impact of the effective contact rate b and the vaccine efficacy on the effective reproduction number R0v. As expected, to reduce the value of R0v below unity, the effective contact rate must be very low, almost irrespective of the vaccine efficacy. That is, despite the availability of vaccines, care-free mixing should continue to be monitored to avoid excessive number of contacts. Similarly, Figure 17 indicates that irrespective of the treatment rate, effective contact rate should be minimized to mitigate the spread of COVID-19 in the population.Figure 16
Contour plot ofR0v with respect to effective contact rate b and vaccine efficacy ε.Figure 17
Contour plot ofR0v with respect to effective contact rate b and treatment rate τ.


## 5.1. Model Fitting

Following the approach described in [14], model system 1 is fitted with daily new COVID-19 cases in Senegal from 29 March to 29 April 2020, which correspond to the first wave of infections in the country. It is important to note that while starting with an expanded model helps to account for the various disease classes explicitly, for the model fitting, model system 1 is reduced to a simpler version which retains the key compartments and characteristics of the complex one by eliminating similar classes [28]. This reduction involves removing the asymptomatic and hospitalized classes, which implies that the reduction in the transmission from hospitalized individuals is equated to zero, that is, ωH=0. The main key requirement for the simple model to approximate the complex one is that the susceptible class S should contain the same number of individuals in both models (when one compares both models compartment wise). The observed data were fitted in MatLab using the optimization function createOptimProblem. The estimated value of the basic reproduction number of the disease in Senegal during the period under consideration is R0=1.31, a value close to unity, which might explain why the first wave of COVID-19 infections did not pick up in Senegal. Figure 5 depicts the daily cumulative number of COVID-19 cases in Senegal. Red dots represent actual confirmed cases while the blue curve is the best fitting curve of the model.Figure 5
Model fit with cumulative daily COVID-19 cases in Senegal, 29 March–29 April 2021.


## 5.2. Long-Term Dynamics of the Disease

We investigated the impact of vaccination and treatment on mitigating the spread of COVID-19. An iterative fourth-order Runge-Kutta method (both forward and backward algorithms) is employed to compute the optimal controls and state values used. For more details on this approach, see [29, 32].The baseline weight parametersA1=A2=A3=A4=1,B1=12,B2=19 are chosen to illustrate the optimal control strategies, as well as the following nonnegative initial conditions S=2500,V=10,E=20,I=70,A=3,H=3,R=1. These weights do not necessarily have a significant meaning attached, but are only of theoretical sense to illustrate our proposed control strategies [31–33, 59]. Using parameter values in Table 1, the reproduction number R0v>1, indicating that the disease is endemic in the population. The positive constants A1, A2, A3, and A4 represent, respectively, the weight that balance off the COVID-19 exposed, infected, asymptomatic, and hospitalized individuals; B1 and B2 are, respectively, the weight constant for vaccination and treatment. Because low cost could potentially be associated with COVID-19 vaccination compared to its treatment (as the cost associated to u2t includes the cost of medical examination, hospitalization, and medicines), the weight factor B1 has been made lower than B2, while the two control strategies u1,u2 are all constrained between zero and one, that is, 0≤uit≤1,i=1,2. For instance, if u1=0, it implies no COVID-19 vaccination, and u1≠0 implies vaccination campaign measures are being implemented in the community. Table 1 provides all the model parameter values used for the simulations.Figure6 depicts the time series of model system 1 for the susceptible, vaccinated, exposed, and infected classes. Because of the potential limitations of vaccines in some settings, we consider in Figure 7 a saturated Holling type II vaccination rate
(35)v≔v1+ξS,where ξ represents the limitation of the availability of vaccine [60]. We observe that when vaccines are widely available, the number of susceptible individuals decreases pretty fast as more and more people are getting the vaccine (the vaccinated class V increases at the onset of implementation (Figure 7)).Figure 6
Profile ofS, V, E, and I without saturation.Figure 7
Profile ofS, V, E, and I with saturation.


## 5.3. Impact of Control Interventions: Vaccination and Treatment

In the next set of figures generated from model system 14, optimal control strategies are implemented. Using the model parameter values in Table1, the basic reproduction number R0=2.21>1. For community with no vaccination program v=0, the basic reproduction number R0=4.67 is almost double compared to the case when vaccination program is introduced. These values are in agreement with those estimated in recent COVID-19 modeling studies [61–63]. Figures 8–12 depict the graphical representations of the simulations of the COVID-19 model as a function of time without control and with optimal control. As can be observed on Figures 8–12, implementing control measures at the optimum level could help to mitigate the spread of the disease as shown by the significant decrease in the total number of individuals in the diseased classes and an increased in the vaccinated class. Figure 13 depicts the profile for the effect of the control functions u1t and u2t on the dynamics of the model system. The two control measures need to be optimally implemented for the first 75 days and maintained for at least 3 months (close to 100 days). It is important however to note that while therapeutic measures, vaccination, and treatment are very effective in curtailing the spread of the epidemic, more control efforts are required to eradicate the disease when R0v>1. Thus, continuously and concurrently applying both pharmaceutical and nonpharmaceutical public health interventions such as face mask, hand washing, and social distancing should be encouraged [64].Figure 8
Dynamics of the vaccinated classV.Figure 9
Dynamics of the exposed classE.Figure 10
Dynamics of the asymptomatic classA.Figure 11
Dynamics of the infected classI.Figure 12
Dynamics of the hospitalized classH.Figure 13
Controlsu1t and u2t.


## 5.4. Sensitivity of the Reproduction NumberR0v

Because mathematical models are symbolic/mechanistic representations of complex biological systems, some model parameter values are not often known with certainty due to natural and seasonal variations, potential measurement errors [18]. To show how changes in model parameters values affect R0v, we determine the relative importance of model parameters to disease transmission and graphically depict how sensitive the effective reproduction number is to the model parameters. Early ranking the intervention measures and other model parameters based on their impact could ideally partially inform the process of prioritizing public health intervention measures to be implemented, thereby helping policy and decision-makers to focus on those key impactful interventions. Figure 14 indicates that when the vaccine efficacy ε is low and R0v is greater than unity, the disease may not be eradicated even if the vaccine coverage v is high. The epidemiological implication is that high vaccine efficacy and vaccination coverage will drastically reduce the number of secondary infections in the community. This result also agrees with the conclusion from Figure 2. The graphical representation (contour plots) of R0v in the parameter space (v,τ) is similar to Figure 14, and this tends to suggest that the impact of the vaccine efficacy is similar to the effect of treatment on the initial disease prevalence.Figure 14
Contour plot ofR0v with respect to vaccine coverage v and efficacy ε.Partial rank correlation describes the relationship between two variables while at the same time removing the effects of several other variables from the relationship [65, 66]. We perform sensitivity analysis by employing the partial rank correlation coefficients (PRCCs) and the Latin hypercube scheme to identify the impact of each model parameters on the initial disease transmission R0v. PRCCs showing the effect of varying the input parameters on the effective reproduction number R0v are shown in Figure 15. All parameters with positive PRCCs will result in an increase on the number of initial disease transmission, while an increase in parameters with negative PRCCs will result in a reduction of R0v. By exploring Figure 15, the most influential parameters can be identified. Model parameters that should be targeted to reduce the spread of the disease are the effective contact rate b, the infection reduction (vaccine efficacy) of vaccinated individuals ε, the fraction of exposed individuals who develop symptoms, and, respectively, the exit rates from the exposed and the asymptomatic classes σ and ϕ.Figure 15
PRCCs showing the effect of varying the input parameters onR0v.Figures16 depicts the impact of the effective contact rate b and the vaccine efficacy on the effective reproduction number R0v. As expected, to reduce the value of R0v below unity, the effective contact rate must be very low, almost irrespective of the vaccine efficacy. That is, despite the availability of vaccines, care-free mixing should continue to be monitored to avoid excessive number of contacts. Similarly, Figure 17 indicates that irrespective of the treatment rate, effective contact rate should be minimized to mitigate the spread of COVID-19 in the population.Figure 16
Contour plot ofR0v with respect to effective contact rate b and vaccine efficacy ε.Figure 17
Contour plot ofR0v with respect to effective contact rate b and treatment rate τ.


## 6. Conclusion

We formulated a deterministic model of the transmission dynamics of COVID-19 with an imperfect vaccine. The model is theoretically analyzed; its effective and basic reproduction numbers are derived. The disease-free equilibrium is globally asymptotically stable, and the disease could be eradicated when the reproduction number is below unity. The critical vaccination threshold is derived, and it is noted that if the vaccine efficacy is low and the disease reproduction number is high, the disease may not be eradicated even if a large proportion of the population is vaccinated. That is, additional efforts will be needed to reduceR0v below unity even if vaccine coverage is high.We then introduce into model system 1 time-dependent control variablesu1t representing vaccination and u2t representing treatment of hospitalized individuals and applied the Pontryagin maximum principle to determine the optimal control strategy for mitigating the spread of the disease. We analytically derived the optimality conditions for disease eradication. The model fits quite well the observed daily data from Senegal early COVID-19 epidemic. Numerical simulations of the optimal control of the full model are carried out using a set of model parameter values. Numerical simulations indicate that COVID-19 can be controlled in the community with the implementation of vaccination and treatment. While our results suggest that vaccination and treatment are very effective in mitigating the spread of COVID-19, more efforts are needed to eradicate the disease. Thus, a combination of or concurrently applying personal protection/preventive measures (nonpharmaceutical public health interventions) such as face masks, hand washing, and social distancing should continue to be encouraged.Finally, we performed a sensitivity analysis using the partial rank correlation coefficient in conjunction with the Latin hypercube sampling technique, to identify the model parameters that significantly influence the initial disease transmissionR0v. Early identification of model parameters with greater influence on disease transmission is important to inform policy decision on which parameters to focus either for data collection or to mitigate the spread of the disease.This study is not exhaustive, and future studies could investigate the impact of both therapeutic and (adherence to) nontherapeutic measures on the dynamics of COVID-19.



---

*Source: 1250129-2021-09-06.xml*