Abstract
The spread and control of coronavirus disease 2019 (COVID-19) present a worldwide economic and medical burden to public health. It is imperative to probe the effect of vaccination and infectivity reductions in minimizing the impact of COVID-19. Therefore, we analyze a mathematical model incorporating different infectivity reductions. This work provides the most economical and effective control methods for reducing the impact of COVID-19. Using data from Ghana as a sample size, we study the sensitivity of the parameters to estimate the contributions of the transmission routes to the effective reproduction number We also devise optimal interventions with cost–benefit analysis that aim to maximize outcomes while minimizing COVID-19 incidences by deploying cost-effectiveness and optimization techniques. The outcomes of this work contribute to a better understanding of COVID-19 epidemiology and provide insights into implementing interventions needed to minimize the COVID-19 burden in similar settings worldwide.
1 Introduction
Coronavirus disease 2019 (COVID-19) is one of the infectious diseases that has caused a global medical and financial burden. Due to its high transmission rate, the World Health Organization (WHO) officially announced the prevalence of it as a global pandemic on 8th March 2020. The coronavirus spreads directly or indirectly from one infectious living host cell to another by replicating and causing infectious disease in the host (humans and/or animals) [1]. Therefore, it is imperative to investigate how the virus spreads in the community and devise strategies to halt the disease’s transmission.
Globally, COVID-19 has caused over 6 million deaths, and by March 2022, approximately 480 million incidences had been reported, according to [2]. However, the COVID-19 reinfection scenario using the SIR model shows that transmission dynamics could arise as a result of immunity waning, even in cases where the force of reinfection is relatively weak [3]. As mentioned in [4], the COVID-19 model, which entails infection through objects contaminated with SARS-CoV-2, is suggested to be made public knowledge.
COVID-19 spreads quickly, threatening global health and igniting a pandemic. This pandemic has had a major effect on other sectors, especially the socioeconomic sector [5]. As of 30 April 2022, there were 512,466,045 reported incidences and 6,257,512 COVID-19 fatalities worldwide, according to the Worldometer. With 6,046,467 COVID-19 cases reported overall, Indonesia ranks seventh in Asia and 18th globally [6].
These findings highlight the need to investigate the spread of COVID-19 in order to minimize its transmission. Using a deterministic model is one way that mathematics is crucial to simulating the epidemic phenomenon of the disease’s spread. The analysis of the COVID-19 model with declining immunity has advanced significantly. For the latter, the natural immunity period is defined by the vaccine efficacy level, which dictates when to start the mass vaccination strategy based on models involving symptomatic and asymptomatic infected populations, as explained in [6].
Mathematical models have been an important tool in epidemiology since the 18th century and can be used to determine the level of spread of infectious disease, aside from medical and biological research studies and strategies for controlling diseases. A modified susceptible–exposed–infectious–removed (SEIR) model with vaccination, quarantine, and isolated SVEQIMR is used to analyze the spread of COVID-19 disease in [7]. In this article, compartmental models were used to demonstrate the pattern of vaccination and control strategies for infectious diseases.
Furthermore, the compartmental model is used to study the behavior of COVID-19 after the introduction of vaccines in [8]. In this article, the population is divided into distinct compartments to demonstrate the control efficacy of infectious diseases in India. Again, a parameterized nonlinear SEIHR model to analyze the transmission of coronavirus disease in Indonesia using the compartmental model is analyzed in [9].
The transmission dynamics, basic reproduction number, and control measures of COVID-19 are analyzed using the compartmental model in [10]. In this article, effective public health interventions were proposed to control the spread of the disease in Wuhan, China. The cost-effective analysis, global stability, and control strategy for the spread of COVID-19 in Ghana are analyzed using the compartmental model in [11]. In this article, control measures were outlined by formulating the SEAIRV model to investigate human–environment–human transmission. However, these studies do not capture the effects of reducing the infectivity rate on the disease’s transmission.
The aforementioned issues have driven this study to explore the effects of vaccination and infectivity reductions on the spread of COVID-19, devise optimal control interventions, and analyze the cost–benefits of implementing the interventions. The subsequent sections are as follows: in Section 2, we present the design of the epidemiological model together with the definition of the parameters and variables in the model. In Section 3, we analyze the positivity of the model’s solutions, computation of the models’ reproduction numbers, and stability of the model. In Section 4, we present the model parameter estimation and sensitivity analysis. In Section 5, we present the numerical simulations and model analysis. In Section 6, we present the optimal control strategies and analyze the cost–benefit analysis. We finally present the concluding remarks of this study in Section 7.
2 model formulation
The flowchart below illustrates the mechanism of the model. We denote as the total population divided into seven different compartments: susceptible individuals at a given time , vaccinated individuals at a given time , exposed individuals at a given time , individuals under quarantine at a given time , infected individuals at a given time isolated individuals at a given time , and individuals recovered at a given time . We assume that the vaccinated individuals become susceptible again due to the vaccine’s inefficacy. Again, a portion of susceptible individuals undergo self-quarantine, while others enter the exposed class. Some of the exposed individuals become infected and then recover naturally without any special treatment. The remaining proportion of individuals either go through quarantine or become infected. It is also assumed that the individuals in the susceptible and vaccinated compartments would come into direct or indirect contact with the individuals in the exposed, quarantined, and isolated compartments. Therefore, we introduce infectivity reductions , and attributed to the exposed, quarantined, and isolated compartments, respectively. The introduction of these different infectivity reductions is to minimize transmission from the potential carriers of the virus. Again, we assume that the quarantined and isolated individuals are not part of the active population at the time of incidence, so the total active population for the standard incidence at any given time is Therefore, the forces for infection of susceptible and vaccinated classes are given as and , respectively. Recovered individuals become susceptible again once they come into contact with the possible carriers of the COVID-19 virus. The description of the model’s parameters is given in Table 1. The optimal control interventions , and are explained in the subsequent section.
TABLE 1
| Parameter | Definition |
|---|---|
| Transmission rate from to | |
| Transmission rate from to | |
| Infectivity reduction of the exposed class | |
| Infectivity reduction of the quarantined class | |
| Infectivity reduction of the isolated class | |
| Λ | Recruitment rate |
| Natural death rate | |
| Vaccine inefficacy | |
| Vaccination rate | |
| Movement rate from to | |
| Movement rate from to | |
| Movement rate from to | |
| Movement rate from to | |
| Movement rate from to | |
| Incubation (latent) period | |
| Movement rate from to | |
| Movement rate from to | |
| Movement rate from to | |
| Movement rate from to |
Definition of model parameters.
The total population in Eq. (1) is defined based on Figure 1 as follows:
FIGURE 1
The following nonlinear ordinary differential equations illustrate the model.where , , and
With initial conditions all initial parameters of model (2) are nonnegative.
3 model mathematical analysis
The detailed analysis of the model is conducted mathematically to show the positivity of solutions and their boundedness within a specific domain.
3.1 Positivity and boundedness of solutions
The variables and parameters of model (2) are nonnegative since it is based on a population with COVID-19.
Define ; if then and its solutions and initial values are nonnegative for and bounded in the region
Proof: Let us consider the following instance where there exists an initial time such that
and for all
Here, Without the loss of generalization, .
Therefore, , and for all However,
This contradicts the claim Therefore, for all This shows that all the solutions are positive for in all other cases.
To prove the boundedness of model (2), define a positive invariant set as and attract positive solutions.
Proof: Considering Eq. 1.
The rate of change in the total population is given as
From Eq. (3), it follows that,
Then, Thus, is bounded, and all solutions in the approach enter or remain in If shows that is a set of positive invariant and is in the region This theorem proves the existence of COVID-19 at a given time in an area that was not infected with COVID-19 disease, and all the initial state variables are positive. This completes the proof.
3.2 Disease-free equilibrium and effective reproduction number
The disease equilibrium point, , of the model is achieved by equating system (2) to 0. The result is defined as follows:
At disease-free equilibrium (DFE), we set , and the following results are obtained from Eq. (4):
At disease-free equilibrium,
The corresponding Jacobian matrix of the system (2) evaluated at (see, Eq. (5)) to obtain the disease-free equilibrium Jacobian matrix is given as follows:where using ; hence,
Using , Eq. (6) the following matrices can be deduced to evaluate the effective reproduction number . We define the next-generation matrix as the square matrix, which consists of matrix representing new infections and matrix V.where
The effective reproduction number, , of the model is evaluated as the spectral radius of matrix ; that is, given aswhere and .
Substituting , and into Eq. (7), gives
Therefore, there is a unique equilibrium in the model, which implies that
in Eq. (
8) is unique. The effective reproduction number,
indicates the number of secondary infections that one infected person can produce if they come into contact with people living in a safe zone.
can be summarized as follows: the exposed class contributed to a significant secondary infection, which is the first term of
where a proportion of
individuals entered the infected class. The infected class generated a significant secondary infection, which is the second term of
The quarantined class contributed to a significant secondary infection, which is the third term of
The isolated class contributed to a significant secondary infection, which is the last term of
We consider the dynamics of vaccination and infectivity reductions on
as follows in Eqs
9–
11respectively.
• In the absence of vaccination, the reproduction number, is given as
• In the absence of infectivity reductions, the reproduction number, is given as
• In the absence of vaccination and infectivity reductions, the reproduction number, is given as
3.3 Stability analysis
Next, since the system has unique equilibrium points, we check its stability.
3.4 Local stability of disease-free equilibrium
The Jacobian matrix, , is given in Equation 7. It can be seen from (6) that are the three eigenvalues of . The other eigenvalues are derived from the reduced matrix in Eq. (12):where and
Therefore, the remaining eigenvalues are the roots of the following characteristic polynomial:where
Therefore, using the Routh–Hurwitz stability conditions, the roots of in Eq. (13) of the reduced matrix have negative real parts if the following conditions hold: and It is obvious that the condition holds if . This proves that in model (2), disease-free equilibrium, , is locally stable if and unstable if .
3.5 Existence of endemic equilibrium of the model
Let us consider as the endemic equilibrium for system (2) and equate the derivative to zero (0). The following results are obtained:where for
Then, and solving Eq. 14 simultaneously gives the following
Substituting the expressions for in Eq. 15 into the fifth equation of Eq. 14 and simplifying, we obtain the following equation for :
Simplifying the above equation gives Eq. (16)where
Note that, in evaluating and
Substituting and and simplifying gives
Let .
It follows that
The expressions for and can be rewritten as follows
Let us consider the denominators of
and
in the following cases:
(a) If , then This results in the endemic equilibrium point since the expressions for and are nonnegative, which implies .
(b) If , then since which implies , endemic equilibrium does not exist.
Using Eq. (17), for when it implies that the endemic equilibrium has a positive and unique equilibrium point when
3.6 Global stability of endemic equilibrium
Let us consider a globally positively definite and unbounded function with a globally negative time derivative. For for all then at equilibrium, is globally stable for the autonomous system and is the Lyapunov function.
The system (14) has an endemic equilibrium that satisfies . If the endemic equilibrium point globally asymptotically stable in a positive region on with
Proof. We define the Lyapunov function and its endemic equilibrium as follows:from Eq. (19), we have
Substituting (2) and (18) into (20) gives
Solving the (21) value gives
Let
are obtained by setting the coefficients of equal to 0, which after solving gives
Hence, the result is
From the above Eqs (22, 23), it can be realized that if the arithmetic mean is greater than or equal to their geometric mean, then only holds if It implies that ; therefore, the largest invariant set is The endemic equilibrium is globally asymptotically stable in the positive region when based on the Lyapunov–LaSalle stability theorem.
4 Model parameterization and sensitivity analysis
We evaluate the reproduction number to find the transmission routes of the pandemic and the proportion of individuals. A sensitivity analysis was carried out to evaluate the contribution of each parameter to the reproduction number.
4.1 Estimation of model parameters
We evaluate the reproduction numbers using the values in Table 3. The effective reproduction number and a proportion of 0.19998 of the exposed individuals move to the infected class at any given time. The contribution from each transmission route is shown in the table below.
Table 2 reveals that the infectious and exposed classes contribute significantly to the transmission of the disease. This implies that there must be admissible measures to halt the transmission of the disease.
TABLE 2
| Total | Exposed () | Infected () | Quarantined ( | Isolated | |
|---|---|---|---|---|---|
| 0.85575 | 0.21686 | 0.41773 | 0.09027 | 0.13088 | |
| 1.07268 | 0.27184 | 0.52362 | 0.11316 | 0.16406 | |
| 4.20958 | 0.72285 | 0.41773 | 0.45137 | 2.61762 | |
| 5.27683 | 0.90612 | 0.52364 | 0.56581 | 3.28126 |
Contribution of transmission routes to the effective reproduction number .
4.2 Sensitivity analysis
We illustrate the arithmetic behavior of model (2) parameters’ sensitivity. The sensitivity of the parameters in for the model is defined in Eq. 24 aswhere
From the above analysis, it is observed that will increase proportionally as the following parameters increase; on the other hand, decreases proportionally as the following parameters increase since sensitivity indices are sign determined. This can be used to show the numerical importance of the various parameters in . Table 3 illustrates the numerical results of the sensitivity index of the various parameters in .
TABLE 3
| Parameter | Value | Reference | Sensitivity index |
|---|---|---|---|
| 0.4531 | [12] | +0.8192 | |
| 0.3531 | Fitted | +0.7104 | |
| 0.3 | [13] | +0.4754 | |
| 0.2 | Fitted | +0.0743 | |
| 0.05 | Fitted | +0.00035 | |
| 0.2 | Fitted | −0.4345 | |
| 0.1 | Estimated | −0.2012 | |
| 0.000043 | Estimated | −0.8784 | |
| 0.15 | Fitted | −0.0674 | |
| 0.2 | [14] | −0.3605 | |
| 0.023 | Fitted | −0.1546 | |
| 0.015 | Fitted | −0.5721 | |
| 0.05 | [15] | −0.0814 | |
| 0.5482 | [16] | −0.2347 | |
| 0.15 | Assumed | −0.1413 | |
| 0.11 | Fitted | −0.1203 |
Values of the sensitivity index of parameters in .
Based on the numerical illustrations shown in Table 3, the effective reproduction number will change as the parameters change. Sensitivity analysis is interpreted based on the sign associated with the particular parameter. The effective reproduction number will decrease as the parameter values with the negative sign increase, while it increases when they decrease. Considering that it indicates that will increase by 81% whenever there is a 1% increment in the transmission rate of the susceptible individuals. This depicts the explanation for all the parameter values with a positive sign. On the other hand, indicates that will decrease by 15% whenever there is a 1% increment in the recovery rate of the infected individuals. This depicts the explanation for all the parameter values with a negative sign.
Now, we examine the relationship between the parameters of effective reproduction number by checking the partial rank correlation coefficients (PRCCs) of the parameters. The following figure illustrates the behavioral pattern of each parameter in the transmission dynamics of the disease.
From Figure 2, one could realize that have a high positive effect on while have a high negative effect on One could realize that all these parameters are associated with vaccination and interventions to mitigate the infectivity rate in the population.
FIGURE 2
4.3 Effects of vaccination and infectivity reductions on the effective reproduction number
This subsection explores the dynamics of with respect to vaccination and infectivity reduction parameters . The following figure demonstrates the changing effects of vaccination and infectivity reduction parameters on
Figure 3A demonstrates the influence of vaccination rate and vaccine inefficacy One can realize that the effective reproduction number increases as the vaccination rate decreases, while decreases as vaccine inefficacy decreases. Figure 3B–g show the effects of vaccination rate , vaccine inefficacy , and infectivity reductions. It is observed that to mitigate the disease’s transmission, there should be effective control measures to minimize the infectivity levels and maximize the vaccination rate to halt the disease’s transmission. Figures 3H–J show the effect of the infectivity reductions on It could be observed that an increase in at any given time results in an increase in (Figure 3H; Figure 3I). It can be seen that and have a significant effect on , as demonstrated in Figure 2. Based on the above graphical representation, it is important to minimize the infectivity levels of the exposed class and quarantined class to halt the disease’s transmission.
FIGURE 3
5 Numerical simulations and discussions
In this section, the real data are observed to depict the situation of the transmission. The real data are compared with the model’s solution accuracy. The proposed model is applied to explore the transmission of COVID-19 disease in Ghana using data from the WHO [17]. This includes the daily number of infections and the total number of daily infections after the introduction of vaccines. The data are analyzed and compared to other literature works for the simulations. The daily infections and the total number of infections are simulated to analyze the behavioral pattern of the transmission of the disease. In particular, we illustrate the changing effects of the vaccination rate , vaccine inefficacy , and infectivity reductions on the size of infectious individuals.
5.1 Application of the model
Once the model is formulated, it is necessary to compare it with data to check its validity. Here, we want to check the accuracy and authenticity of the model by verifying the extent to which the model can represent the real situation, as described in [18]. The simulation consists of the application of the data from Ghana to illustrate the transmission of COVID-19 for the period 1st January 2022 to 1st March 2022 when the individuals were vaccinated. The following state variables are considered using data from Ghana for the period 1st January 2022 to 1st March 2022 [17]. The estimated total population of Ghana is [19]; hence, and the assumed initial values are as follows: All the parameters used for the simulations are shown in Table 3.
The SVEQIMR model (fitted) depicts the pattern of the real situation (real data), as shown in Figure 4A, which represents the plot of the daily number of infections, and Figure 4B represents the total number of infections for the period of 1st January 2022 to 1st March 2022. The plot of the model results indicates the pattern of the real situation. That is, the fitted model and real situation agree with each other and illustrate the transmission dynamics of the spread of the pandemic.
FIGURE 4
5.2 Effects of vaccination and infectivity reductions on COVID-19 incidences
Here, we analyze the influence of vaccination rate and vaccine inefficacy on the spread of COVID-19 disease using the same dataset. Let us consider the following parameters: vaccination rates and vaccine inefficacy for the infected class.
Again, we demonstrate the changing effects of the different infectivity reductions on the behavioral pattern of the transmission of COVID-19 disease by considering the following arbitrary values for the infectivity reduction of exposed individuals , quarantined individuals , and isolated individuals . The results are presented in Figure 5.
FIGURE 5
From Figure 5A, it is observed that an increase in the vaccination rate results in a significant decrease in the daily number of reported cases. On the contrary, an increase in vaccine inefficacy results in an increase in the daily number of reported cases, as shown in Figure 5B. It could also be observed that the disease dies out gradually with time as vaccine inefficacy approaches zero (0); however, a high vaccination rate reduces the number of infections. Therefore, there is a significant effect of different values of vaccination rate and vaccine inefficacy on the disease’s spread. The above illustration shows that vaccination against the spread of the coronavirus is very important and must be adhered to.
From Figure 5C, it is realized that an increase in the infectivity reduction of exposed individuals results in a significant increase and higher peaks in the daily reported cases. Furthermore, an increase in the infectivity reduction of the quarantined class results in a moderate increase in the daily reported cases compared to the peaks of , as shown in Figure 5D. The peaks of daily reported cases are lower than those of the exposed class because the individuals in this group are confined and monitored. In Figure 5E, there is a significant decline in the daily reported cases as the infectivity reduction of the isolated individuals decreases, which has fewer peaks than the others. From the above graphical representation, the infectivity reduction of the exposed individuals should be controlled carefully because the higher the infectivity reduction, the higher the number of daily reported cases and vice versa. This graphical representation can be confirmed by the numerical illustration of the contribution of the transmission route in Table 2. The exposed compartment contributed significantly to the effective reproduction number , with a proportion of 0.19998 of the individuals moving to the infected class. This means that approximately 20% of individuals in the exposed compartment become infected at a given time. The above diagram depicts the order of contribution of the transmission route to the effective reproduction number
6 Optimal control problem and cost–benefit analysis
We modify Equation 2 with the following optimal control variables: representing the public awareness of the prevalence of COVID-19 and related infections where , which reduces the forces of infection and , by ; denoting the control effort to intensify COVID-19 vaccination; and denoting the control effort for infectivity reduction, which reduces infectivity by . All these efforts denote the control interventions in minimizing the transmission of COVID-19 and its reinfections. The modified Eq. (2) is given in Eq. (25) as follows:where and remain as in Section 2 with the initial conditions given in the model (2).
6.1 Objective functional
We now formulate the optimal trajectories that show the effect of the control efforts subjected to (36); the objective functional is given as
We focus on minimizing the cost function (26), and the total cost of implementing the optimal control is given as
The parameters and in Eq. 27 are the balancing cost factors for , respectively. All the control efforts are assumed to be bounded by Lebesgue measurable time-dependent functions on the interval , where is the final time with the control effort set defined in Eq. 28 as
Now, we establish point-wise Hamiltonian through Pontryagin’s maximum principle to transform the optimal control system (25) and its associated objective functional (26). The following optimal solution is achieved.
where in Eq. 29 are the co-state variables with respect to the state variables,
Given as the optimal controls and the corresponding solutions of the system (25), which minimizes over then there exist co-state variables that satisfy
with conditions where Then, the optimality conditions that minimize the Hamiltonian, H, of (29) with respect to the controls are given as
Proof: We formulate the adjoint equation for the optimal system by taking the partial derivative of Eq. 30 as follows in Eq. (32):
The control set illustrates the co-state system with the optimal conditions.
We solve for , and , and of Eq. (33), and the results confirms the expression in Eq. (31) are as follows:
Therefore, using the bounds of the controls , and the control efforts are in the compact form given by the optimal condition of the system in Eq. (34); hence, the proof is complete.
6.2 Optimal control strategies
Here, our aim is to determine the number of infections after deploying the optimal control interventions. We explore the effects of implementing the interventions; therefore, the optimality system (36) is solved forward in time and the adjoint system backward in time with the corresponding lower and upper bounds of the controls. We used the population of Ghana to study the behavioral pattern of COVID-19. The estimated total population of Ghana is [19]; hence, and the assumed initial values are as follows: together with Λ = 1,364, , and and parameter values illustrated in Table 3. The balance costs associated with the objective functional are assumed to be and weight , where The lower bound (LB) and upper bound (UB) are assumed to be The results are illustrated according to the strategies to implement the control strategies.
6.2.1 Strategy 1: implementation of public awareness ()
The optimal solutions illustrated in Figure 6 account for the observations when the control effort is applied accordingly.
FIGURE 6
The optimal solutions illustrated above depict the following observations when public education is only applied:
(a) Figure 6A represents the effect of the control effort on the infectious individuals. It implies that the number of individuals will decrease if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 6B represents the profile of the control effort for public awareness of COVID-19. It implies that education on COVID-19 should reach 50% of the population throughout the implementation to halt COVID-19 transmission.
6.2.2 Strategy 2: implementation of vaccination ()
The optimal solutions illustrated in Figure 7 account for the observations when the control effort is applied accordingly.
FIGURE 7
The optimal solutions illustrated above depict the following observations when vaccination is only applied:
(a) Figure 7A represents the effect of the control effort on the infectious individuals. It implies that the number of individuals will decrease if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 7B represents the profile of the control effort for vaccination to prevent COVID-19. It implies that approximately 25% of the population should be vaccinated within 80 days and intensified further up to 75% in the subsequent days throughout the implementation to halt COVID-19 transmission.
6.2.3 Strategy 3: implementation of infectivity treatment ()
The optimal solutions illustrated in Figure 8 account for the observations when the control effort is applied accordingly.
FIGURE 8
The optimal solutions illustrated above depict the following observations when infectivity treatment is only applied:
(a) Figure 8A represents the effect of the control effort on the infectious individuals. It implies that the number of individuals will decrease if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 8B represents the profile of the control effort for COVID-19 reinfection. It implies that approximately 75% of the individuals suspected to be carriers of the virus should be treated/monitored throughout the implementation to halt COVID-19 transmission.
6.2.4 Strategy 4: implementation of public awareness and vaccination ()
The optimal solutions illustrated in Figure 9 account for the observations when the control efforts are applied accordingly.
FIGURE 9
The optimal solutions illustrated above depict the following observations when control efforts for public awareness and vaccination are applied:
(a) Figure 9A represents the effect of the control efforts on the infectious individuals. It implies that the number of individuals will decrease if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 9B represents the profile of the control efforts for public awareness and vaccination against COVID-19. It implies that approximately 80% of the population should be educated on COVID-19 within 22 days, which can be relaxed to approximately 25%, and 50% of the population should be vaccinated throughout the implementation of these interventions to halt COVID-19 transmission.
6.2.5 Strategy 5: implementation of public awareness and infectivity treatment ()
The optimal solutions illustrated in Figure 10 account for the observations when all the control efforts are applied accordingly.
FIGURE 10
The optimal solutions illustrated above depict the following observations when all the control efforts for public awareness and infectivity treatment are applied:
(a) Figure 10A represents the effect of the control efforts on the infectious individuals. It implies that the number of individuals will decrease to a minimum within 30 days if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 10B represents the profile of the control efforts for public awareness and infectivity treatment. It implies that approximately 80% of the population should be educated on COVID-19, and 50% of suspected carriers of the virus should be treated/monitored throughout the implementation of these interventions to halt COVID-19 transmission.
6.2.6 Strategy 6: implementation of vaccination and infectivity treatment ()
The optimal solutions illustrated in Figure 11 account for the observations when the control efforts are applied accordingly.
FIGURE 11
The optimal solutions illustrated above depict the following observations when all the control efforts for vaccination infectivity reduction are applied:
(a) Figure 11A represents the effect of the control efforts on the infectious individuals. It implies that the number of individuals will decrease to a minimum within 20 days if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 11B represents the profile of the control efforts for public education and vaccination against COVID-19. It implies that all the interventions should be implemented at levels higher than 25% from the start of implementation throughout the subsequent days to halt COVID-19 transmission.
6.2.7 Strategy 7: implementation of all controls ()
The optimal solutions illustrated in Figure 12 account for the observations when the control efforts are applied accordingly.
FIGURE 12
The optimal solutions illustrated above depict the following observations when all the control efforts are applied:
(a) Figure 12A represents the effect of the control efforts on the infectious individuals. It implies that the number of individuals will decrease to a minimum within 10 days if the control effort is optimally implemented in halting the disease’s transmission. Conversely, it will increase significantly.
(b) Figure 12B represents the profile of all control efforts. It implies that approximately 80% of the population should be educated on COVID-19, approximately 50% of the population should be vaccinated, and 50% of suspected carriers of the virus should also be treated/monitored throughout the implementation period of these interventions to halt COVID-19 transmission.
6.3 Cost–benefit analysis
Once the strategies are given, it is imperative to know the cost associated with implementing such intervention(s). Therefore, we explore the cost associated with each control strategy to check their effectiveness. We outline some cost-effectiveness approaches to further understand the control strategies.
We consider two procedures, namely, average cost-effectiveness ratio (ACER) and incremental cost-effectiveness ratio (ICER), which have been explained in [20–23], to carry out epidemiological studies.
6.4 Average cost-effectiveness ratio
We define the ACER of implementing a strategy as
The overall cost stated in (27) would be used to evaluate the total cost that the intervention would generate in Eq. 36. We then compare the ACER values of each strategy, and the one with the lowest value is the most cost-effective, saving costs. Therefore, the cost-effective intervention is considered the strategy with the least ACER value. The expression in Eq. 35 is illustrated as follows.
From Table 4, control strategy 1, which involves the implementation of public education only, has the least value of ACER, indicating cost savings. However, relying solely on this metric is not enough to choose a strategy; hence, we further explore other approaches.
TABLE 4
| Strategy | Overall infection averted | Overall cost involved | ACER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 2 | |||
| Strategy 3 | |||
| Strategy 4 | |||
| Strategy 5 | |||
| Strategy 6 | |||
| Strategy 7 |
Strategies’ ACER values with their overall infection averted and cost incurred.
6.5 Incremental cost-effectiveness ratio
We define the ICER of implementing a strategy as
The total cost function stated in (27) would be used to estimate the overall cost that the intervention would generate in Eq. (36). It is worth knowing that the averted total number of infections is the difference between the initial values of , without control(s) and with controls. The outcomes are tabulated in increasing order of infection averted.
The ICER in Table 5 is calculated as
TABLE 5
| Strategy | Overall infection averted | Overall cost incurred | ICER value |
|---|---|---|---|
| Strategy 2 | |||
| Strategy 1 | |||
| Strategy 4 | |||
| Strategy 3 | |||
| Strategy 5 | |||
| Strategy 6 | |||
| Strategy 7 | 0.0017 |
Strategies’ ICER values with their overall infection averted and cost incurred.
Assessing strategies 2 and 1 in Table 5, it is noticed from the ICER that strategy 2 is expensive to deploy in a resource-limited setting; therefore, strategy 2 is removed from the list of possible controls, and the ICER is calculated again. This is presented in Table 6.
TABLE 6
| Strategy | Overall infection averted | Overall cost incurred | ICER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 4 | |||
| Strategy 3 | |||
| Strategy 5 | |||
| Strategy 6 | |||
| Strategy 7 | 0.0017 |
Strategies’ ICER values with their overall infection averted and cost incurred.
Assessing strategies 1 and 4 in Table 6, it is noticed from the ICER that strategy 4 is expensive to deploy in a resource-limited setting; therefore, strategy 4 is removed from the list of possible controls, and the ICER is calculated again. This is presented in Table 7.
TABLE 7
| Strategy | Overall infection averted | Overall cost incurred | ICER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 3 | |||
| Strategy 5 | |||
| Strategy 6 | |||
| Strategy 7 | 0.0017 |
Strategies’ ICER values with their overall infection averted and cost incurred.
Assessing strategies 1 and 3 in Table 7, it is noticed from the ICER that strategy 3 is expensive to deploy in a resource-limited setting; therefore, strategy 3 is removed from the list of possible controls, and the ICER is calculated again. This is presented in Table 8.
TABLE 8
| Strategy | Overall infection averted | Overall cost incurred | ICER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 5 | |||
| Strategy 6 | |||
| Strategy 7 | 0.0017 |
Strategies’ ICER values with their overall infection averted and cost incurred.
Assessing strategies 1 and 5 in Table 8, it is noticed from the ICER that strategy 5 is expensive to deploy in a resource-limited setting; therefore, strategy 5 is removed from the list of possible controls, and the ICER is calculated again. This is presented in Table 9.
TABLE 9
| Strategy | Overall infection averted | Total cost involved | ICER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 6 | |||
| Strategy 7 | 0.0017 |
Strategies’ ICER values with their overall infection averted and cost involved.
Assessing strategies 1 and 6 in Table 9, it is noticed from the ICER that strategy 6 is expensive to deploy in a resource-limited setting; therefore, strategy 6 is removed from the list of possible controls, and the ICER is calculated again. This is presented in Table 10.
TABLE 10
| Strategy | Overall infection averted | Overall cost incurred | ICER value |
|---|---|---|---|
| Strategy 1 | |||
| Strategy 7 |
Strategies’ ICER values with their overall infection averted and cost incurred.
Finally, assessing strategies 1 and 7 in Table 10, it is noticed from the ICER that strategy 7 is expensive to deploy in a resource-limited setting; therefore, strategy 7 is removed from the list of possible controls. Therefore, we conclude that strategy 1 is the most cost-effective strategy to use among the several strategies under study here. From the above analysis, it is obvious that strategy 1, which involves public education, is the intervention that saves cost.
7 Conclusion
We have presented a work that analyzes the changing effects of vaccination and infectivity reductions on the transmission of COVID-19 using data from Ghana. We have estimated the model’s parameters and analyzed their effects on disease transmission through numerical and graphical illustrations. Again, we have exhibited the threshold dynamics of the effective reproduction number together with the contributions from the transmission routes (Table 2). We have demonstrated the sensitivity of the model’s parameters to study their effects on the effective reproduction number (Table 3; Figure 2; Figure 3).
The aim of this work is to study the effect of vaccination and infectivity reductions in controlling COVID-19 transmission and devise control interventions that save cost to mitigate the transmission; therefore, we have formulated optimal control strategies together with the cost–benefit analysis that consider control measures involving both pharmaceutical and non-pharmaceutical interventions to control COVID-19. We implemented the strategies (Figure 6A–Figure 12A), and it was realized that public education, vaccination, and infectivity reductions to prevent COVID-19 should be intensified and reach approximately 25% of the population from the beginning and intensified in the subsequent days (Figure 6B–Figure 12B).
It is also worth knowing that public education saves cost as per the cost–benefit analysis compared to the other strategies raised in this work. This intervention can minimize COVID-19, as illustrated in Figure 6A. This intervention should reach approximately 50% of the population throughout the period of its implementation in order to realize the results of strategy 1 (Figure 6). Although strategy 1 saves cost, other strategies elaborated in this work can also be applied, but one has to consider the cost involved in implementing the strategy. The cost involved in applying the optimal control strategies is presented in Table 4–Table 10.
The outcomes of the findings imply that both pharmaceutical and non-pharmaceutical measures are very important in controlling the transmission of COVID-19. These control measures should always be vigorously implemented to create public awareness on COVID-19 and its reinfection, as illustrated in Figure 6B–Figure 12B, in order to reduce the effective contact rates and rates of acquiring COVID-19, as illustrated in Figure 3.
Although we have demonstrated the dynamics of COVID-19 transmission with vaccination and different infectivity reductions, this work is focused on the homogeneity of the population, and we hope to extend this study to explore the transmission dynamics of COVID-19 reinfection by considering heterogeneity of the population such as age and sex. We encourage individuals to adhere to personal hygiene and be aware of COVID-19 reinfection.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Author contributions
RA: conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, and writing–review and editing. ZJ: funding acquisition, supervision, and writing–review and editing. JY: formal analysis, supervision, and writing–review and editing. JA: formal analysis, investigation, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the National Natural Science Foundation of China grants 12231012 and 61873154.
Acknowledgments
The authors appreciate the support of Complex Systems Research Center, Shanxi University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
World Health Organization. Coronavirus disease 2019(COVID-19), situation report -51, data as reported by 8 March 2020 (2019). Available at: https://www.who.int/emergencies/disease/novel-coronavirus-2019/situation-reports (Accessed January 10, 2024).
2.
ChenYWangYFlemingJYuYGuYLiuCet alActive or latent tuberculosis increases susceptibility to COVID-19 and disease severity (2020). Available at: https://www.medrxiv.org/ (Accessed January 10, 2024).
3.
SalmanAMAhmedIMohdMHJamiluddinMSDheyabMA. Scenario analysis of COVID-19 transmission dynamics in Malaysia with the possibility of reinfection and limited medical resources scenarios. Comput Biol Med (2021) 133:104372. 10.1016/j.compbiomed.2021.104372
4.
ZamirMNadeemFAlqudahMAAbdeljawadT. Future implications of COVID-19 through Mathematical modeling. Results Phys (2022) 33:105097. 10.1016/j.rinp.2021.105097
5.
Worldometers. Reported cases and deaths by country or territory (2022). Available at: https://www.worldometers.info/coronavirus/#countries (Accessed April 30, 2022).
6.
ZegarraMAAInfanteSDCarrascoDBLiceagaDO. COVID-19 optimal vaccination policies: a modeling study on efficacy, natural and vaccine-induced immunity responses. Math Biosci (2021) 337:108614. 10.1016/j.mbs.2021.108614
7.
AdachiA. Grand challenge in human/animal virology: unseen, smallest replicative entities shape thewholeglobe. FrontMicrobiol (2019) 11:431. 10.3389/fmicb.2020.00431
8.
1mg. Coronavirus No-panic help guide (2020). Available at: https://smef.org.uk/wpcontent/uploads/2020/03/Corona-Ebook.pdf. (Accessed January 10, 2024)
9.
AnnasSPratamaMIRifandiMSanusiWSideS, Stability analysis and numerical simulation of SEIR model for pandemic COVID-19 spread in Indonesia. Chaos, Solitons and Fractal (2020) 139:110072. 10.1016/j.chaos.2020.110072
10.
PangLLiuSZhangXTianTZhaoZ, Transmission dynamics and control strategies of COVID-19 in Wuhan, China. J Biol Syst (2020) 28:543–60. 10.1142/s0218339020500096
11.
AsamoahJKKOwusuMAJinZOduroFAbidemiAGyasiEO, Global stability and cost-effectiveness analysis of COVID-19 considering the impact of the environment: using data from Ghana. Chaos, Solitons and Fractals (2020) 140:110103. 10.1016/j.chaos.2020.110103
12.
TchoumiSDiagneMRwezauraHTchuencheJ, Malaria and COVID-19 co-dynamics: a mathematical model and optimal control. Appl Math Model (2021) 99:294–327. e327. 10.1016/j.apm.2021.06.016
13.
YanXZouY, Optimal and sub-optimal quarantine and isolation control in SARS epidemics. Math Comp Model (2008) 47(1-2):235–45. 10.1016/j.mcm.2007.04.003
14.
KhanMAAtanganaA, Mathematical modeling and analysis of COVID-19: a study of new variant omicron. Physica A: Stat Mech Its Appl (2022) 599:127452. Article 127452. 10.1016/j.physa.2022.127452
15.
RwezauraHDiagneMOmameAde EspindolaATchuencheJ, Mathematical modeling and optimal control of SARS-CoV-2 and tuberculosis co-infection: a case study of Indonesia. Model Earth Syst Environ (2022) 8(4):5493–520. e5520. 10.1007/s40808-022-01430-6
16.
OmameAAbbasMOnyenegechaC, A fractional-order model for COVID-19 and tuberculosis co-infection using Atangana–Baleanu derivative. Chaos, Solitons and Fractals (2021) 153:111486. Article 111486. 10.1016/j.chaos.2021.111486
17.
World Health Organization. World health organization covid– 19 report (2024). Available at: https://covid19.who.int/region/afro/country/gh. (Accessed January 10, 2024).
18.
DietzK. The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res (1993) 2:23–41. 10.1177/096228029300200103
19.
OseiEAmuHKye-DuoduGKwablaMPDansoEBinkaFNet alImpact of COVID-19 pandemic on Tuberculosis and HIV services in Ghana: an interrupted time series analysis. PLoS ONE (2023) 18(9):e0291808. 10.1371/journal.pone.0291808
20.
AsamoahJKKOwusuMAJinZOduroFTAbidemiAGyasiEO, Global stability and cost-effectiveness analysis of COVID-19 considering the impact of the environment: using data from Ghana. Chaos, Solitons and Fractals (2020) 140:110103. 10.1016/j.chaos.2020.110103
21.
AgustoFLeiteM, Optimal control and cost-effective analysis of the 2017 meningitis outbreak in Nigeria. Infect Dis Model (2017) 4(2019):161–87. 10.1016/j.idm.2019.05.003
22.
AsamoahJKKOkyereEAbidemiAMooreSESunGJinZet alOptimal control and comprehensive cost-effectiveness analysis for COVID-19. Results Phys (2022) 33:105177. 10.1016/j.rinp.2022.105177
23.
AsamoahJKKJinZSunG, Non-seasonal and seasonal relapse model for Q fever disease with comprehensive cost-effectiveness analysis. Results Phys (2021) 22:103889. 10.1016/j.rinp.2021.103889
Summary
Keywords
COVID-19, vaccination, infectivity reductions, optimal control, cost–benefit analysis
Citation
Appiah RF, Jin Z, Yang J and Asamoah JKK (2024) Cost–benefit analysis of the COVID-19 vaccination model incorporating different infectivity reductions. Front. Phys. 12:1383357. doi: 10.3389/fphy.2024.1383357
Received
07 February 2024
Accepted
19 March 2024
Published
08 May 2024
Volume
12 - 2024
Edited by
Xiaoke Xu, Dalian Nationalities University, China
Reviewed by
Xiaofeng Luo, North University of China, China
Segun Oke, Ohio University, United States
Updates
Copyright
© 2024 Appiah, Jin, Yang and Asamoah.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhen Jin, jinzhn@263.net; Junyuan Yang, yjyang66@sxu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.