 Research Article
 Open access
 Published:
Impact of vaccine supplies and delays on optimal control of the COVID19 pandemic: mapping interventions for the Philippines
Infectious Diseases of Poverty volumeÂ 10, ArticleÂ number:Â 107 (2021)
Abstract
Background
Around the world, controlling the COVID19 pandemic requires national coordination of multiple intervention strategies. As vaccinations are globally introduced into the repertoire of available interventions, it is important to consider how changes in the local supply of vaccines, including delays in administration, may be addressed through existing policy levers. This study aims to identify the optimal level of interventions for COVID19 from 2021 to 2022 in the Philippines, which as a developing country is particularly vulnerable to shifting assumptions around vaccine availability. Furthermore, we explore optimal strategies in scenarios featuring delays in vaccine administration, expansions of vaccine supply, and limited combinations of interventions.
Methods
Embedding our work within the local policy landscape, we apply optimal control theory to the compartmental model of COVID19 used by the Philippine governmentâ€™s pandemic surveillance platform and introduce four controls: (a) precautionary measures like community quarantines, (b) detection of asymptomatic cases, (c) detection of symptomatic cases, and (d) vaccinations. The model is fitted to local data using an LBFGS minimization procedure. Optimality conditions are identified using Pontryaginâ€™s minimum principle and numerically solved using the forwardâ€“backward sweep method.
Results
Simulation results indicate that early and effective implementation of both precautionary measures and symptomatic case detection is vital for averting the most infections at an efficient cost, resulting in \(>99\%\) reduction of infections compared to the nocontrol scenario. Expanding vaccine administration capacity to 440,000 full immunizations daily will reduce the overall cost of optimal strategy by \(25\%\), while allowing for a faster relaxation of more resourceintensive interventions. Furthermore, delays in vaccine administration require compensatory increases in the remaining policy levers to maintain a minimal number of infections. For example, delaying the vaccines by 180 days (6 months) will result in an \(18\%\) increase in the cost of the optimal strategy.
Conclusion
We conclude with practical insights regarding policy priorities particularly attuned to the Philippine context, but also applicable more broadly in similar resourceconstrained settings. We emphasize three key takeaways of (a) sustaining efficient case detection, isolation, and treatment strategies; (b) expanding not only vaccine supply but also the capacity to administer them, and; (c) timeliness and consistency in adopting policy measures.
Graphic Abstract
Background
In the year since the emergence of the global coronavirus disease 2019 (COVID19) pandemic, national policies have had to decisively manage diverse issues of resource availability, institutional capacity, and collective behavioral change [1,2,3]. Striking the right balance of multiple strategies at the right time has been vital for implementing successful pandemic responses [4, 5]. Mathematical modelling has helped scientists and policymakers incorporate emergent discoveries about COVID19 with existing knowledge to design effective interventions [6, 7].
In early 2021, the global introduction of vaccination as a viable counter to the disease prompts new analytical efforts. Regional inequalities introduce challenges to the global vaccine supply chain which may disrupt a straightforward vaccine rollout for a significant proportion of various national populations [8, 9]. Important questions emerge with respect to how governments may adequately adjust existing policies available for pandemic control in relation to multiple scenarios.
This paper undertakes an optimal control study of policies to control the COVID19 outbreak in the Philippines, a developing country that may be particularly vulnerable to experiencing challenges to vaccine rollouts. Amidst largescale preparations for the evaluation, selection, and distribution of vaccines, ongoing policies to respond to the pandemic continue to inform the Philippine governmentâ€™s strategies for pandemic management [10]. Questions around their optimal implementation are particularly salient for developing countries that face heavier burdens from both the pandemic and overly restrictive quarantine measures [11]. This study therefore asks: How should the Philippine government implement existing strategies of community quarantine and case detection in conjunction with the introduction of vaccine rollouts?
Related work
Mathematical modelling for forecasting COVID19 outbreaks
Since the beginning of the COVID19 pandemic, the academic literature has witnessed a vast surge of modelling studies. Existing reviews highlight the importance of compartmental models of COVID19, in connection with other models based on time series forecasting and machine learning [12, 13]. Compartmental models mathematically encode known and emerging information about the transmission dynamics of the disease and have been locally applied across numerous contexts around the world, including major sites of COVID19 transmission like China, India, Brazil, the United States, and the United Kingdom [14,15,16,17,18].
Mathematical modelling efforts have been beneficial for forecasting and intervention assessment [19,20,21]. For instance, in the United Kingdom, a stochastic, agestructured transmission model was used to quantify the costs and mortalities of unmitigated outbreaks without interventions, highlighting the need for sustaining combined control efforts [22]. In another example, an agestructured model with social contact matrices was used to compare the impacts of different reopening strategies on the relative reduction of cases in different regions in China [23].
Optimal control theory for modelling pandemic response
In this work, we utilize optimal control theory to model effective pandemic response. Optimal control theory refers to a field of study that deals with finding optimal solutions to a problem expressed in the form of a nonlinear dynamical system [24, 25]. This helps identify efficient methods of achieving desired outcomes, such as costeffective infection control [15, 26].
Numerous studies have implemented optimal control theory toward similar end goals. In the absence of vaccines, most early studies focused on nonpharmaceutical interventions, including various combinations of rapid testing, contact tracing, and awareness campaigns [27,28,29,30,31]. Ullah and colleagues sought to disentangle the impacts of quarantine and case detection rates on exposed, critical, and hospitalized COVID19 patients [32]. Other research modelled the effect of limited total testing resources, through the addition of an isoperimetric constraint to the optimal control problem [33].
Eventually, however, newer research was further able to consider the impacts of eventual vaccine availability. In an agestructured model, Bonnans and Gianatti studied minimization of the death toll, cost of confinement, and hospitalization peaks discussing a possible extension of their model when a vaccine becomes available [34]. Libotte and colleagues likewise explored programs for vaccine administration within a multiobjective setup, determining a set of Pareto optimal strategies that would minimize infections while also minimizing the number of vaccines needed [35].
Responding to COVID19 in the Philippines
The present work specifically investigates the dynamics of COVID19 in the Philippines and aims to identify optimal strategies for efficiently controlling infections. We draw on existing modelling efforts by the national pandemic surveillance system [36] to derive realistic parameters which match existing epidemic trends and available intervention strategies in the country [37]. By deploying models informed by local parameters of the disease, we therefore aim to provide both theoretically optimal and contextually practical recommendations for policymakers [38].
In the Philippines, nonpharmaceutical interventions have primarily included phased community quarantines and mandated wearing of face masks [39]. Enhancements to the capacity of the health system to efficiently detect asymptomatic and symptomatic infectious individuals have also been key [40, 41]. In early 2021, imminent vaccine rollouts posed a salient new factor for pandemic control. We therefore sought to design optimal strategies for their distribution, and consider appropriate responses to potential obstacles which may arise in resourceconstrained settings.
Aims of the current study
Burgeoning scholarship points to rich global knowledge of the effectiveness and efficiency of various policy tools against the pandemic. However, both nationally specific impacts of the pandemic and the limitations faced by intervening bodies highlight the importance of grounding optimal control analysis in the local context.
In this view, the present work therefore aims to achieve the following goals. First, we frame pandemic interventions with vaccinations as an optimal control problem to identify scenarios for effective pandemic control. Second, we explore various vaccination scenarios featuring both delays and expansions of vaccine administration. This enables a futureoriented analysis of how local policymakers may compensate for unforeseen developments in the global supply chain. Finally, we perform a systematic ablation analysis, whereby we restrict various combinations of available controls to model more limited control scenarios.
Methods
Model formulation
To capture the local dynamics of COVID19 transmission in the country, we form a model that utilizes the local incidence data from the Department of HealthEpidemiology Bureau (DOHEB) [42]. The COVID19 model utilizes six compartments to subdivide the population: susceptible (S), exposed (E), infectious but asymptomatic \((I_a)\), infectious and symptomatic \((I_s)\), confirmed (C), and removed (R). These compartments are governed by the epidemic flow as illustrated in Fig. 1. Compartment S consists of individuals who have not been infected with COVID19 but may contract the disease once exposed to the virus. Compartment E consists of individuals who have been infected but are still within the latency stage of the disease. These individuals will eventually become infectious, and categorized into two compartments depending on the presence (\(I_s\)) or absence (\(I_a\)) of symptoms. Once detected, these infectious individuals will move to compartment C, where they will be included among the active cases. Individuals in this compartment are assumed to be isolated, and hence not capable of infecting the susceptible population, and receiving treatment while in isolation. Lastly, compartment R consists of individuals who have acquired immunity from the disease. We assume that those who have recovered from the disease acquire permanent immunity and therefore move to compartment R.
The movements of individuals toward and out of the different compartments are governed by several parameters of the model. The transmission rate \(\beta\) is a function of the disease transmission rate \(\beta _0\), based on an assumed basic reproduction number \(R_0\), and a reduction factor \((1  \lambda )\). The parameter \(\lambda\) reflects the effect of community quarantine imposed by the government, as well as the degree of compliance to minimum health standards, which includes practicing proper hygiene, social distancing, and wearing protective face coverings. Moreover, the parameter \(\psi\) accounts for the infectiousness of asymptomatic individuals relative to those who have symptoms. The rates of transfer to the two infectious compartments, \(\alpha _a\) for asymptomatic and \(\alpha _s\) for symptomatic, are both dependent on the incubation period \(\tau\) of the virus.
Other parameters in the model include the constant recruitment rate A into the S compartment, which is driven by new births in the population. To account for deaths by natural causes, a constant rate of \(\mu\) per unit time is applied to all compartments in the model. In addition, deaths due to the disease are included through the parameters \(\epsilon _I\) and \(\epsilon _T\), affecting the infectious symptomatic and confirmed compartments, respectively.
By taking into account the above assumptions, a mathematical model is developed, which can be described by the following system of six ordinary differential equations:
where \(\beta = \beta _0 (1\lambda )\), \(\alpha _a = \frac{c}{\tau }\), \(\alpha _s = \frac{1c}{\tau }\), \(N = S + E + I_a + I_s + C + R\). The functions S, E, \(I_a\), \(I_s\), C, and R are differentiable realvalued functions on \({\mathbb {R}}\). Moreover, all parameters are nonnegative constants.
Parameter values
The values of various parameters were determined using several sources and methods. Local COVID19 data [42] were used in calculating detection rate (\(\delta _a,\delta _s\)), postdetection recovery rate (r), and death rate of COVID19 cases (\(\epsilon _I\), \(\epsilon _T\)). The recruitment rate (A) and natural death rate (\(\mu\)) are calculated from population data [43,44,45]. For the other parameters which cannot be computed directly from data, we rely on references that estimate their values. The basic reproduction number of COVID19 (\(R_0\)) and relative infectiousness of asymptomatic cases (\(\psi\)) are based on estimates by the US Centers for Disease Control [46]. The incubation rate (\(\tau\)) and symptomatic transition (\(\omega\)) are obtained from reports by the World Health Organization [47, 48]. The proportion of asymptomatic cases (c) is from Mizumoto et al. [49] that studied the COVID19 outbreak at the Diamond Princess cruise ship.
We fit the model output to data by employing a curvefitting algorithm to estimate the value of the transmission reduction rate (\(\lambda\)), the initial values for the exposed (E(0)), infectious asymptomatic (\(I_a(0)\)), and infectious symptomatic (\(I_s(0)\)). In particular, the constrained LBFGS optimization procedure [50] was utilized to minimize the sum of squared errors between the model output and the empirical time series. The parameter \(\lambda\) is fitted on a permonth basis starting from March 2020, to coincide with the changes in the disease dynamics and the corresponding transmission reduction policies implemented by the government that tends to be updated monthly [51]. The output is a vector of bestfit transmission reduction parameters \([\lambda _1,\lambda _2,...,\lambda _n]\) where n is the number of months since March 2020. This fitting procedure is utilized to produce forecasts for the Philippine COVID19 epidemic [36]. Following the above parametrization, our model fits well to the Philippine data for cumulative cases of COVID19 (Fig. 2). Table 1 summarizes the parameter values used in the model.
Model with optimal control
We explore four control strategies to mitigate the COVID19 epidemicprecautionary measures, detection of asymptomatic cases, detection of symptomatic cases, and vaccination. The definitions of each of these controls and how they are incorporated in the model are given below:

1.
Precautionary measures (\(u_1(t)\)) refer to governmentled efforts to inhibit possible contacts between susceptible and infectious individuals by regulating public gatherings, closing schools, suspending office work, enforcing adherence to health protocols such as social distancing, maskwearing, handwashing, etc. This control affects the transmission rate \(\beta\) and is incorporated in the model as a factor (\(1u_1(t)\)), replacing (\(1  \lambda\)). The value of \(u_1(t)\) represents the effort of precaution at time t. A value of 0 indicates that no precautionary measure is being practiced, and a value of 1 indicates full effort on precaution prohibiting any form of infection.

2.
Detection of asymptomatic cases (\(u_2(t)\)) entails identifying and isolating infectious individuals who do not have symptoms of COVID19. This may be done through laboratory tests such as reverse transcriptionpolymerase chain reaction (RTPCR) to determine whether an individual is infectious or not. A positive case is taken to be immediately followed by isolation at home or in a dedicated quarantine facility to prevent transmission. It is assumed, therefore, that after an individual is confirmed to have COVID19, s/he is not able to infect susceptible individuals. We incorporate this control to the model by replacing \(\delta _a\) with a timevarying control function \(u_2(t)\). The value of \(u_2(t)\) represents the effort of testing and isolation at a given time t. A value of 0 indicates the absence of testing and isolating, and a value of 1 indicates testing and isolating all infectious asymptomatic individuals on a given unit of time.

3.
Detection of symptomatic cases (\(u_3(t)\)) follows the same definition as the detection of asymptomatic cases but applied to individuals that exhibit symptoms of COVID19. We replace \(\delta _s\) by \(u_3(t)\) to incorporate this control to the model. Similarly, a value of 0 of this control indicates the absence of effort to test and isolate symptomatic individuals while a value of 1 indicates full testing and detection of all symptomatic individuals on a given unit of time.

4.
Vaccination (\(u_4(t)\)) refers to the full inoculation of susceptible individuals for them to acquire protection against COVID19 infection or protection against a severe case of the disease. We assume in this paper that vaccines give protection against infection, that is, an individual who is fully vaccinated gets immunity to COVID19 over the period considered. Multiple vaccines with varying effectiveness rates have been identified for use against COVID19 such as those developed and manufactured by PfizerBioNTech, Moderna, Sinovac, etc. We consider the average effectiveness rate of the vaccines weighted by the usage, denoted by \(\sigma\), where \(0\le \sigma \le 1\). To incorporate vaccination in the model, we add a rate of transfer from compartment S to R equal to \(\sigma u_4(t)\). The value of \(u_4(t)\) represents the effort of vaccination for the susceptible population. A value of 0 represents no vaccination efforts while a value of 1 represents vaccination of all susceptible individuals on a single unit of time.
Our goal is to identify the optimal strategy for limiting the spread of SARSCoV2 in a population using minimal cost of controls. In this study, the optimal control problem minimizes the number of asymptomatic (\(I_a\)) and symptomatic individuals (\(I_s\)) and the control costs. The controls are expressed in quadratic forms to incorporate nonlinear costs for the implementation of each control and to ensure the convexity of the cost function. This is a common form of an objective functional in optimal control problems [25, 52]. The objective functional is represented by:
where \(t_0\) and \(t_f\) represent January 1, 2021 and December 31, 2022 respectively, reflecting a 2year period. The parameters \(w_i, i=1,2,3,4,\) account for the relative costs of implementing controls \(u_i\). They represent the weights of corresponding terms in the integrand and their importance in the optimal control problem.
We aim to identify \(u_i^*(t), i=1,2,3,4,\) such that
where for Lebesgue integrable \(u_i\),
Here, \(u_{i}^{\mathrm{{min}}}\) and \(u_{i}^{\mathrm{{max}}}\) are the lower and upper bounds of the control \(u_i\), representing minimum and maximum implementation efforts.
The constraints of the optimal control problem are given by:
where
The existence of the optimal solution can be shown using standard results in optimal control theory [25, 52]. The necessary convexity of the integrand of the objective functional, positive definiteness of system (4), and the linear dependence of the state differential equations to the controls are satisfied in our model.
We apply Pontryaginâ€™s minimum principle [24] to determine the necessary conditions using the optimality system for our problem (see Additional file 1: Appendix). This system is a twopoint boundary problem with initial conditions for the state variables and terminal conditions for the adjoint variables. The solutions are solved numerically using a Rungeâ€“Kutta fourthorder scheme. The state variables are solved forward in time while the adjoint variables are solved backwards, referred to as Forwardâ€“Backward Sweep Method [25]. We update the controls using a convex combination of the latest and previous values. This process is iterated until the updates in the control values are very small or less than the machine epsilon.
The initial state values are computed using a combination of data and model fitting. We relied on model fitting to data to get the values for \(E(0),I_a(0)\) and \(I_s(0)\) [36]. The initial value for confirmed cases (C(0)) is based on data from the Department of Health [42]. The initial number of removed individuals (R(0)) is assumed to be higher than the detected recoveries on January 1, 2021, to include recoveries from undetected asymptomatic cases. We estimate that this is equal to 450,000, consistent with the output of our model [36]. Lastly, the initial susceptible population is estimated to be equal to the whole population minus the assumed values for the other compartments. Table 3 lists the initial state values used in our simulations.
We fix \(u_{i}^{\mathrm{min}}=0\) for \(i\in \{1,2,3,4\}\) while upper bounds of the controls are varied to reflect the realistic maximum efforts that can be achieved with each control. Results from model fitting to data (see Table 1) show that the highest value for precaution is 0.85. Direct computations from epidemiological data provided by the Department of Health [42] show that the minimum monthly average duration of detection, from symptom onset to confirmation of test results, may take 5 days. We take the inverse of this duration as our upper bound for both detection controls, hence \(u_{2}^{\mathrm{{max}}} = u_{3}^{\mathrm{{max}}} = 0.2\). Lastly, the upper bound for vaccination (\(u_4\)) is based on government proclamations [10]. The boundaries for the control values are summarized in Table 2.
The weight parameters \(w_i, i=1,2,3,4,\) are adjusted to balance the terms in the integrand of the cost function (2). These parameters reflect the total costs and payoffs of implementing each control strategy including the cost of the products used (test kits, vaccines, masks, etc.), operational costs (personnel salary, rent, procurement, refrigeration units, etc.), opportunity costs for the economy due to lockdowns, and so forth. To determine the values of the weight parameters, we consider the fact that the upper bounds for the controls already reflect the realistic and achievable efforts that the Philippine government can exert, given the historical and prospective cost and availability of each control. Recall that the upper bounds for precaution, detection of asymptomatic cases, and detection of symptomatic cases are based on data, and the upper bound for the vaccination control is based on government targets. Based on this, we assume that a lower \(u_i^{max}\) signifies a relatively higher implementation cost for the ith control as implementing the control beyond this upper bound is not readily available to the government. Given this, we rescale the terms in the cost function (2) by equating the weight parameters to the inverse of the maximum allowable effort for each control (\(w_1=1/0.85, w_2=w_3=1/0.2, w_4 = 1/0.002\)). This improves the balance of the terms in the cost function and reduces the bias to implement controls that have lower upper bounds.
Limited information is available as of writing to estimate the vaccine effectiveness parameter (\(\sigma\)) for COVID19. Moreover, various vaccines with different effectiveness rates will be deployed in the Philippines as they become available [10], making it more difficult to give a realistic estimate of this parameter. The best alternative is to equate this to the pooled effectiveness of vaccines for a similar disease such as influenza, which is at \(\sigma =0.7\) [53].
Results
Optimal control strategies for COVID19 in the Philippines
Solving the stated optimal control problem in Eqs. 2â€“4 generates the optimal levels of precaution (\(u_1\)), asymptomatic detection (\(u_2\)), symptomatic detection (\(u_3\)), and vaccination (\(u_4\)) over the 2year period from January 1, 2021 to December 31, 2022 (Fig. 3). Notably, the maximum feasible vaccination rate must be sustained throughout the entire 2year period. Symptomatic detection must likewise maintain a high value close to the maximum feasible value, lowering slightly in the first 6 months of 2021. Asymptomatic detection and precaution must likewise be implemented at their maximum respective values early in 2021. But asymptomatic detection may be eased to nearly zero by the second quarter of 2021, while precaution eases a bit by the second half of 2021, then is gradually reduced throughout the remainder of the 2 years under consideration.
We also simulate a nocontrol scenario by setting the controls to 0 throughout the 2year period. A dramatic difference is observed between the withcontrol and withoutcontrol conditions (Fig. 4). Without controls, a peak number of 100 million infectious individuals is achieved within the first quarter of 2021. Meanwhile, with the optimal implementation of all controls, the total number of running infections is driven down quickly in early 2021, without ever breaching the 10,000mark. After the fullthrottle implementation of all controls in early 2021, sustained efforts at detecting symptomatic individuals and proactively vaccinating susceptible populations may thus be sufficient to prevent the infected population from rising. So long as the latter strategies are maintained, the majority of the population may slowly ease stringent distancing rules and fewer resources need to be urgently allocated to asymptomatic detection.
Policy impacts of vaccine delays
With the bottleneck in the global supply of vaccines, it is of chief concern when a country can start vaccinating its population. It is not farfetched for countries to experience delays in vaccination which in turn, would have an effect on policy. Here, we look into the impact of vaccine delay on the optimal control strategy. To achieve this, we add the following constraint to the optimal control problem 2â€“4:
where \(t_d\) is the vaccine delay in days. We first solve for the optimal control profiles given vaccine delays of 180, 360, and 540 days (Fig. 5). Results reveal that increased efforts on the other controls become necessary given longer delays in vaccination. Primarily, precautionary measures should compensate when vaccines are delayed. Detection of symptomatic infectious individuals should also be strengthened for mitigation if vaccine rollout is slowed down.
To further evaluate the effect of vaccine delay, we compute the cost of the optimal strategy in each scenario. The cost of the control strategy (\({\mathscr {C}}\)) is defined as the integral of the last four terms in the cost function of the optimal control problem over the time period, specifically:
We observe that the optimal strategy in the nodelay scenario has the least cost and will also result in the least number of total infections. Given this, we decided to compare the cost of the optimal strategy in the scenarios with vaccine delay relative to the nodelay scenario. Specifically, given vaccine delays of 30k days, where \(k\in \{0,1,2,3,...,24\}\), we examine the resulting relative cost and total infections of the optimal control strategy (Fig. 6).
Based on our simulations, longer delays in vaccination result in higher relative costs of the optimal strategy. For example, delaying the vaccines by 180 days will result in an \(18\%\) increase in cost, and delaying the vaccine by 360 days will increase the cost by \(32\%\), due to the compensation of the other controls. We also observe that delaying the vaccine by 60 days or more will only increase the total infections in the optimal strategy by a relatively small amount (\(<1000\) additional infections). These two findings suggest that while the number of infections can still be effectively managed when vaccines are delayed, vaccine delay may pose more deleterious effects on the economy than on the overall health status of the population, even if the optimal strategy is implemented and new cases are minimized.
Policy impacts of expanding vaccine supply
Recall that the upper bound for the vaccination control (\(u_{4}^{\mathrm{{max}}}\)) was fixed based on the vaccination plan by the local government. Note however that the actual vaccine capacity is unknown and is dependent on negotiations and supply. Here, we want to look into whether increasing vaccine supply will have a significant effect on the optimal control strategy. To discern this relationship, we modify the value of \(u_{4}^{\mathrm{{max}}}\) to double and triple the initial value. The value of the weight parameter \(w_4\) is equal to \(1/u_{4}^{\mathrm{{max}}}\) in each scenario.
We solve for the optimal control profiles and the resulting number of vaccinations if \(u_{4}^{\mathrm{{max}}}=0.002,0.004\) or 0.006 (Figs. 7, 8). We observe that increasing the vaccine capacity will have a significant impact on the optimal control strategy. For the three scenarios considered, the maximum vaccination effort must be utilized for almost the entire period, but vaccination effort is eased earlier if the vaccine capacity is larger. Another important consequence of increasing vaccine capacity is the earlier relaxation of the other controls, mainly precautionary measures and detection of symptomatic cases.
Comparing the relative cost and the resulting total infections reveals that increasing the vaccine capacity by double or triple the initial amount will reduce the cost of the optimal strategy (Table 4). We observe a 25% cost reduction when vaccine supply is doubled, and 37% cost reduction when vaccine supply is tripled, coupled with a slight reduction in the total number of infections. This reinforces the proposition that dedicating more resources to vaccinations is more favorable in the long run owing to the reduced efforts necessary for implementing the other interventions.
Managing cost and impacts of pandemic control strategies
Finally, to integratively consider the dynamics of all interventions, we analyzed the results of a control setup featuring all interventions in conjunction with ablated control scenarios featuring various subsets of the proposed controls. We compared outcomes for optimized single control, dual control, and triple control strategies to simulate scenarios when the other controls are not available, as well as to shed light on their contributions to controlling the epidemic, and highlight the significance of implementing all four in concert. To do this, controls that are not being implemented in each scenario are fixed at 0. Full details on various control profiles are available in Additional file 1: Appendix.
We compare all possible combinations of the controls in terms of both their cost and the infections averted relative to the nocontrol scenario (Fig. 9). Intuitively, an ideal scenario entails low cost and high infections averted. Strikingly, the control scenarios appear to cluster into three major categories. First are low cost, low impact strategies, which do not entail high costs, but also do not effectively curb infections. These correspond to intervention programs that do not mobilize sufficient resources to address the health crisis and subsequently do not achieve the desired impact. We observe here that this cluster of intervention combinations primarily exclude precautionary measures like community quarantines, meaning various scenarios implementing only vaccinations and case detection strategies. This entails that, even if these strategies might be less costlyespecially from an economic perspectivethan prolonged lockdowns, they may not be sufficient on their own to control outbreak trajectories. Especially in the early months of 2021, it will be vital for local governments to limit unnecessary contact between individuals, and enforce such procedures reliably and consistently. Otherwise, even maximally implemented case detection and vaccination strategies will not be able to protect a significant proportion of the population from infection.
The second category represents the worstcase scenario: high cost, low impact strategies. These indicate attempts by governing entities to invest resources in public health interventions, which ultimately still do not effectively control outbreaks. This therefore presents a severe misuse of resources without achieving desired outcomes. Note here that these intervention combinations primarily exclude efficient detection of symptomatic cases. This means that without efficient detection of symptomatic casesâ€”even at full implementation of precautionary quarantine measures, vaccinations, and overall high costs to the economy at largefew infections will be averted.
The final category represents the most favorable category of interventions. Here, medium cost, high impact strategies pertain to scenarios involving some investment of resources, directed towards the most efficient policy levers. This results in an effective minimization of infections, thereby constituting welltargeted policy decisions that achieve the objective of controlling the pandemic. Now we see that both precautionary measures and efficient symptomatic case detection are vital to achieving this set of outcomes. Even if these interventions do introduce higher costs, they can effectively quell outbreak trajectories early on. Moreover, in these setups, their implementation is even given an allowance for relaxation over time if executed effectively and consistently in the early months. This therefore reduces costs from a broader perspective, as overall fewer infections arise nationally, and fewer resources are demanded to address them.
Discussion
In this study, we studied optimal strategies for controlling the spread of COVID19 in the Philippines. We considered existing policy interventions as well as the introduction of vaccination rollouts to quell outbreak trajectories. Dramatic differences were detected in simulated infections depending on which controls were prioritized. In particular, we observed the importance of early and effective implementation of precautionary measures like community quarantines, coupled with efficient detection of symptomatic cases. Furthermore, we found that even if vaccinations alone do not constitute an efficient response to the pandemic, expanding vaccine supply relaxes the need for these more resourceintensive interventions. Meanwhile, although less than ideal, delays in vaccine administration may also be compensated through the remaining policy levers.
These insights bear particular consequences for policies in developing countries like the Philippines [1, 2]. Here, we highlight three key takeaways. First, more than a year into the pandemic, it remains crucial to sustain efficient case detection, isolation, and treatment strategies, particularly for symptomatic cases. In the Philippines, where longterm states of community lockdown have prevailed as the governmentâ€™s response to shortterm fluctuations in COVID19 cases [39], our findings suggest that an optimal, costeffective strategy would actually entail relaxations to such measuresâ€”but only under the condition that symptomatic case detection is properly implemented. Hence, improving the capacity of the local health system to identify, process, and manage these cases efficiently should be a top priority beyond cyclically adjusting quarantine levels [54].
Second, policymakers need to consider how to expand not just vaccine supply, but also the capacity to administer them. This includes both logistical concerns regarding the strategic use of facilities to vaccinate individuals, inform the public regarding vaccine availability and eligibility, as well as reducing vaccine hesitancy through culturally sensitive health promotion programs that strengthen public trust [3, 55, 56]. Only when such health communication objectives are accomplished and collective behavioral change is initiated can the vaccination strategies posited in this work be made feasible. Otherwise, even procuring sufficient supplies of vaccines will not achieve its intended effects to stop local outbreaks. Initial efforts along these lines are underway by the Philippine National Vaccine Operations Center, for which our model results strongly reaffirm the urgency.
Third, timeliness and consistency must be emphasized in adopting policy measures [4, 5]. Across all favorable scenarios simulated, high levels of key interventions were needed in the early months of 2021, with relaxations projected only mid2021 or in 2022. Systems for detecting existing cases, while preventing new ones, are needed for vaccinations to meaningfully impact outbreak trajectories and reduce overall costsâ€”especially as these systems need to be robustly sustained in the event of potential delays in acquiring sufficient vaccines for the entire population. This ensures that even if developing countries like the Philippines do not hold sway over the global supply chain of vaccines, the pandemic may still be kept under control through means over which local policymakers do wield authority.
It is worth noting that the conclusions of this work rest on several assumptions. These assumptions constrain the interpretation of our findings but likewise point to promising avenues for future work [6, 7]. A number of limitations pertain to the realism of our model. For instance, we assume a total population that is unaffected by immigration and emigration flows. This can be a source of confounding given that, despite additional precautions in travel protocols, relevant susceptible, exposed, and infectious populations, in reality, include individuals who may leave or arrive within national borders. Furthermore, ordinary differential equation models of diseases, such as the one utilized in this paper, assume homogeneous mixing of individuals in the population. That is, our model assumes that a susceptible individual has a uniform chance of being infected by any infectious individual, regardless of their geographic proximity. However, in an archipelago such as the Philippines, this assumption of freemixing does not necessarily hold due to different patterns of movement within the country and the tendency of the outbreaks to be concentrated within more urban areas. In relation to this heterogeneity, we further hypothesize that a multiregion approach to mitigation such as in [57] will further lower the total cost of the optimal control strategies. Hence, though we do not include such migratory flows and heterogeneity in our analysis, extensions may valuably consider these factors as well.
We additionally assume that vaccines work by transitioning individuals to a removed compartment. However, it is not the case that all vaccines guarantee 100% immunity for all individuals. Additionally, for parsimony, our model does not incorporate a number of wideranging issues which remain pressing to address, yet are beyond the scope of this work. These include: prioritized vaccinations of various segments of the population, lesserknown dynamics of reinfection with COVID19, the variability in the economic cost of the controls, the impact of emerging variants of the pathogen, documented distinctions between being protected from symptoms while being able to transmit the virus to others, or the practical circumstances of administering vaccines requiring multiple doses [58, 59]. With regard to this latter limitation, we specifically do not model potential logistical impediments in vaccine scheduling or temporary states of partial protection resulting from initial doses [60]. Such considerations therefore place caveats on the implications of our results, further highlighting the need for robust investment in these strategies when translating them into realworld policies. These factors may likewise be modelled with greater precision in succeeding work as growing knowledge continues to accumulate in line with close monitoring by the scientific community [61].
Conclusions
This study applies optimal control theory to an epidemiological model to calculate the optimal efforts required for precautionary measures, asymptomatic case detection, symptomatic case detection, and vaccination to mitigate the impact of the COVID19 pandemic. Using parameter values suitable to the Philippines, we show that precautionary measures and symptomatic case detection are essential interventions to minimize infections at an efficient cost. Furthermore, relaxation of measures is feasible after an early and maximal implementation of all controls. Our results also highlight that increasing vaccination capacity and timely acquisition of vaccines are key to reducing the total implementation cost, leading to earlier relaxation of the other nonpharmaceutical interventions in the optimal strategy. This work provides a quantitative reference for drafting policies designed to control the pandemic in the most efficient manner.
Availability of data and materials
The data analysed for this study are available in the Department of Health COVID19 tracker https://www.doh.gov.ph/covid19tracker.
Abbreviations
 COVID19:

Coronavirus disease 2019
References
Ahmed F, Ahmed N, Pissarides C, Stiglitz J. Why inequality could spread COVID19. Lancet Public Health. 2020;5(5):e240.
Chiriboga D, Garay J, Buss P, Madrigal RS, Rispel LC. Health inequity during the COVID19 pandemic: a cry for ethical global leadership. Lancet (Lond, Engl). 2020;395(10238):1690.
Van Bavel JJ, Baicker K, Boggio PS, Capraro V, Cichocka A, Cikara M, et al. Using social and behavioural science to support COVID19 pandemic response. Nat Hum Behav. 2020;4(5):460â€“71.
Haug N, Geyrhofer L, Londei A, Dervic E, DesvarsLarrive A, Loreto V, et al. Ranking the effectiveness of worldwide COVID19 government interventions. Nat Hum Behav. 2020;4(12):1303â€“12.
Ruktanonchai NW, Floyd J, Lai S, Ruktanonchai CW, Sadilek A, RenteLourenco P, et al. Assessing the impact of coordinated COVID19 exit strategies across Europe. Science. 2020;369(6510):1465â€“70.
Becker AD, Grantz KH, Hegde ST, BÃ©rubÃ© S, Cummings DA, Wesolowski A. Development and dissemination of infectious disease dynamic transmission models during the COVID19 pandemic: what can we learn from other pathogens and how can we move forward? Lancet Digit Health. 2020;3:e41â€“50.
Vespignani A, Tian H, Dye C, LloydSmith JO, Eggo RM, Shrestha M, et al. Modelling COVID19. Nat Rev Phys. 2020;2(6):279â€“81.
Nelson R. COVID19 disrupts vaccine delivery. Lancet Infect Dis. 2020;20(5):546.
Usher AD. COVID19 vaccines for all? Lancet. 2020;395(10240):1822â€“3.
Department of Health. The Philippine National Deployment and vaccination plan for COVID19 vaccines. 2021. https://doh.gov.ph/node/27220. Accessed 16 Feb 2021.Â
Loayza NV. Costs and tradeoffs in the fight against the COVID19 pandemic: a developing country perspective. Washington, DC: World Bank; 2020.
Mohamadou Y, Halidou A, Kapen PT. A review of mathematical modeling, artificial intelligence and datasets used in the study, prediction and management of COVID19. Appl Intell. 2020;50(11):3913â€“25.
Rahimi I, Chen F, Gandomi AH. A review on COVID19 forecasting models. Neural Comput Appl. 2021. https://doi.org/10.1007/s00521020056268.
IHME COVID19 Forecasting Team. Modeling COVID19 scenarios for the United States. Nat Med. 2020;27(1):94â€“105.
Lin F, Muthuraman K, Lawley M. An optimal control theory approach to nonpharmaceutical interventions. BMC Infect Dis. 2010;10(1):1â€“13.
Neto OP, Kennedy DM, Reis JC, Wang Y, Brizzi ACB, Zambrano GJ, et al. Mathematical model of COVID19 intervention scenarios for SÃ£o PauloBrazil. Nat Commun. 2021;12(1):1â€“13.
Samui P, Mondal J, Khajanchi S. A mathematical model for COVID19 transmission dynamics with a case study of India. Chaos Solitons Fract. 2020;140:110173.
Jia J, Ding J, Liu S, Liao G, Li J, Duan B, et al. Modeling the control of COVID19: impact of policy interventions and meteorological factors. Electron J Differ Eq. 2020;23(23):1â€“24.
Chatterjee R, Bajwa S, Dwivedi D, Kanji R, Ahammed M, Shaw R. COVID19 risk assessment tool: dual application of risk communication and risk governance. Prog Disaster Sci. 2020;7:100109.
Pluchino A, Biondo A, Giuffrida N, Inturri G, Latora V, Le Moli R, et al. A novel methodology for epidemic risk assessment of COVID19 outbreak. Sci Rep. 2021;11(1):1â€“20.
Sangiorgio V, Parisi F. A multicriteria approach for risk assessment of Covid19 in urban district lockdown. Saf Sci. 2020;130:104862.
Davies NG, Kucharski AJ, Eggo RM, Gimma A, Edmunds WJ, Jombart T, et al. Effects of nonpharmaceutical interventions on COVID19 cases, deaths, and demand for hospital services in the UK: a modelling study. Lancet Public Health. 2020;5(7):e375â€“85.
Liu Y, Gu Z, Xia S, Shi B, Zhou XN, Shi Y, et al. What are the underlying transmission patterns of COVID19 outbreak? An agespecific social contact characterization. EClinicalMedicine. 2020;22:100354.
Pontryagin LS, Boltyanskii VG, Gamkrelize RV, Mishchenko EF. The mathematical theory of optimal processes. New York: Wiley; 1962.
Lenhart S, Workman JT. Optimal control applied to biological models. New York: Chapman and Hall/CRC; 2007.
Rowthorn RE, Laxminarayan R, Gilligan CA. Optimal control of epidemics in metapopulations. J R Soc Interface. 2009;6(41):1135â€“44.
Tsay C, Lejarza F, Stadtherr MA, Baldea M. Modeling, state estimation, and optimal control for the US COVID19 outbreak. Sci Rep. 2020;10(1):10711.
Perkins TA, EspaÃ±a G. Optimal control of the COVID19 pandemic with nonpharmaceutical interventions. Bull Math Biol. 2020;82(9):118.
Sasmita NR, Ikhwan M, Suyanto S, Chongsuvivatwong V. Optimal control on a mathematical model to pattern the progression of coronavirus disease 2019 (COVID19) in Indonesia. Glob Health Res Policy. 2020;5(1):38.
Madubueze CE, Dachollom S, Onwubuya IO. Controlling the spread of COVID19: optimal control analysis. Comput Math Methods Med. 2020;2020:1â€“14.
Obsu LL, Balcha SF. Optimal control strategies for the transmission risk of COVID19. J Biol Dyn. 2020;14(1):590â€“607.
Ullah S, Khan MA. Modeling the impact of nonpharmaceutical interventions on the dynamics of novel coronavirus with optimal control analysis with a case study. Chaos Solitons Fract. 2020;139:110075.
Ndii MZ, Adi YA. Modelling the transmission dynamics of COVID19 under limited resources. Commun Math Biol Neurosci. 2020.
Bonnans JF, Gianatti J. Optimal control techniques based on infection age for the study of the COVID19 epidemic. Math Model Nat Phenom. 2020;15:48.
Libotte GB, Lobato FS, Platt GM, Neto AJS. Determination of an optimal control strategy for vaccine administration in COVID19 pandemic treatment. Comput Methods Programs Biomed. 2020;196:105664.
FASSSTER. COVID19 Philippines LGU monitoring platform. 2020. https://fassster.ehealth.ph/covid19/. Accessed 29 May 2021.
Estuar MRJE, Uyheng J, De Leon M, Benito DJ, De LaraTuprio E, Estadilla C, et al. Science and public service during a pandemic: reflections from the scientists of the Philippine Governmentâ€™s COVID19 surveillance platform. Philipp Stud Hist Ethnogr Viewp. 2020;68(3):493â€“504.
Uyheng J, Pulmano CE, Estuar MRJ. Deploying system dynamics models for disease surveillance in the Philippines. In: International conference on social computing, behavioralcultural modeling and prediction and behavior representation in modeling and simulation. Springer; 2020. p. 35â€“44.
Vallejo BM Jr, Ong RAC. Policy responses and government science advice for the COVID 19 pandemic in the Philippines: January to April 2020. Prog Disaster Sci. 2020;7:100115.
Buhat CAH, Duero JCC, Felix EFO, Rabajante JF, Mamplata JB. Optimal allocation of COVID19 test kits among accredited testing centers in the Philippines. J Healthc Inform Res. 2021;5(1):54â€“69.
Caldwell JM, deÂ LaraTuprio E, Teng TRY, Estuar MRJE, Sarmiento RF, Abayawardana M, etÂ al. Understanding COVID19 dynamics and the effects of interventions in the Philippines: a mathematical modelling study. medRxiv. 2021;p. 2021â€“01.
Department of HealthEpidemiology Bureau. COVID19 Tracker Philippines. 2020. https://www.doh.gov.ph/covid19tracker. Accessed 8 May 2020.
macrotrends. Philippines birth rate 1950â€“2020. 2020. https://www.macrotrends.net/countries/PHL/philippines/birthrate. Accessed 8 May 2020.
macrotrends. Philippines life expectancy 1950â€“2020. 2020. https://www.macrotrends.net/countries/PHL/philippines/lifeexpectancy. Accessed 8 May 2020.
Philippine Statistics Authority. Census of population and housing. 2020. https://psa.gov.ph/populationandhousing. Accessed 8 May 2020.
US Centers for Disease and Control. COVID19 pandemic planning scenarios. 2021. https://www.cdc.gov/coronavirus/2019ncov/hcp/planningscenarios.html. Accessed 20 Mar 2021.
World Health Organization. Coronavirus disease 2019 situation report73. 2020. https://www.who.int/docs/defaultsource/coronaviruse/situationreports/20200402sitrep73covid19.pdf. Accessed 8 May 2020.
World Health Organization. Report of the WHOChina Joint Mission on Coronavirus Disease 2019 (COVID19). World Health Organization; 2020. https://www.who.int/docs/defaultsource/coronaviruse/whochinajointmissiononcovid19finalreport.PDF.
Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance. 2020;25(10):2000180.
Varadhan R, et al. Numerical optimization in R: Beyond optim. J Stat Softw. 2014;60(1):1â€“3.
Official Gazette. Interagency task force for the management of emerging infectious diseases resolutions. 2020. https://www.officialgazette.gov.ph/section/laws/otherissuances/interagencytaskforceforthemanagementofemerginginfectiousdiseasesresolutions/. Accessed 15 Nov 2020.
Fleming W. Deterministic and stochastic optimal control. New York: Springer; 1975.
Osterholm MT, Kelley NS, Sommer A, Belongia EA. Efficacy and effectiveness of influenza vaccines: a systematic review and metaanalysis. Lancet Infect Dis. 2012;12(1):36â€“44.
World Health Organization. COVID19 strategy update; 2020. https://www.who.int/emergencies/diseases/novelcoronavirus2019/strategiesandplans. Accessed 8 May 2020.
Dror AA, Eisenbach N, Taiber S, Morozov NG, Mizrachi M, Zigron A, et al. Vaccine hesitancy: the next challenge in the fight against COVID19. Eur J Epidemiol. 2020;35(8):775â€“9.
Volpp KG, Loewenstein G, Buttenheim AM. Behaviorally informed strategies for a national COVID19 vaccine promotion program. JAMA. 2021;325(2):125â€“6.
Carli R, Cavone G, Epicoco N, Scarabaggio P, Dotoli M. Model predictive control to mitigate the COVID19 outbreak in a multiregion scenario. Annu Rev Control. 2020;50:373â€“93.
Logunov DY, Dolzhikova IV, Shcheblyakov DV, Tukhvatulin AI, Zubkova OV, Dzharullaeva AS, et al. Safety and efficacy of an rAd26 and rAd5 vectorbased heterologous primeboost COVID19 vaccine: an interim analysis of a randomised controlled phase 3 trial in Russia. Lancet. 2021;397(10275):671â€“81.
Malkov E. Simulation of coronavirus disease 2019 (COVID19) scenarios with possibility of reinfection. Chaos Solitons Fract. 2020;139:110296.
SaadRoy CM, Morris SE, Metcalf CJE, Mina MJ, Baker RE, Farrar J, et al. Epidemiological and evolutionary considerations of SARSCoV2 vaccine dosing regimes. Science. 2021;372:363â€“70.
del Rio C, Malani P. COVID19 in 2021continuing uncertainty. JAMA. 2021;325:1389â€“90.
Acknowledgements
CDSE would like to thank Aurelio A. de los Reyes V, Dr. rer. nat. for his invaluable feedback on the initial design and findings of the study. All authors would like to thank the anonymous reviewers for their constructive comments.
Funding
This work was supported by the Department of Science and TechnologyPhilippine Council for Health Research and Development (DOSTPCHRD) and the United Nations Development Programme (UNDP) Pintig Lab through the FASSSTER project.Â This work is also supported by the Rizal Library Open Access Journal Publication Grant of the Ateneo de Manila University.
Author information
Authors and Affiliations
Contributions
CDSE contributed to the design of the work and coding of the optimal control algorithms. CDSE and JU contributed to the presentation of the results. EDLT, CDSE, JMRM, JU, and TRT contributed to the design of the model. All authors contributed to the drafting of the work. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Supplementary Information
Additional file 1: Appendix S1.
The Appendix is included as a separate file as part of this manuscript submission and contains supplementary information on: (1) Optimality conditions, (2) Control profiles for (ablated) single, dual, and triple control scenarios.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Estadilla, C.D.S., Uyheng, J., de LaraTuprio, E.P. et al. Impact of vaccine supplies and delays on optimal control of the COVID19 pandemic: mapping interventions for the Philippines. Infect Dis Poverty 10, 107 (2021). https://doi.org/10.1186/s40249021008865
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40249021008865