Analysis of the Tradeoff Between Health and Economic Impacts of the Covid-19 Epidemic

Various measures have been taken in different countries to mitigate the Covid-19 epidemic. But, throughout the world, many citizens don't understand well how these measures are taken and even question the decisions taken by their government. Should the measures be more (or less) restrictive? Are they taken for a too long (or too short) period of time? To provide some quantitative elements of response to these questions, we consider the well-known SEIR model for the Covid-19 epidemic propagation and propose a pragmatic model of the government decision-making operation. Although simple and obviously improvable, the proposed model allows us to study the tradeoff between health and economic aspects in a pragmatic and insightful way. Assuming a given number of phases for the epidemic (namely, 4 in this paper) and a desired tradeoff between health and economic aspects, it is then possible to determine the optimal duration of each phase and the optimal severity level (i.e., the target transmission rate) for each of them. The numerical analysis is performed for the case of France but the adopted approach can be applied to any country. One of the takeaway messages of this analysis is that being able to implement the optimal 4−phase epidemic management strategy in France would have led to 1.05 million of infected people and a GDP loss of 231 billions € instead of 6.88 millions of infected and a loss of 241 billions €. This indicates that, seen from the proposed model perspective, the effectively implemented epidemic management strategy is good economically, whereas substantial improvements might have been obtained in terms of health impact. Our analysis indicates that the lockdown/severe phase should have been more severe but shorter, and the adjustment phase occurred earlier. Due to the natural tendency of people to deviate from the official rules, updating measures every month over the whole epidemic episode seems to be more appropriate.

Various measures have been taken in different countries to mitigate the Covid-19 epidemic. But, throughout the world, many citizens don't understand well how these measures are taken and even question the decisions taken by their government. Should the measures be more (or less) restrictive? Are they taken for a too long (or too short) period of time? To provide some quantitative elements of response to these questions, we consider the well-known SEIR model for the Covid-19 epidemic propagation and propose a pragmatic model of the government decision-making operation. Although simple and obviously improvable, the proposed model allows us to study the tradeoff between health and economic aspects in a pragmatic and insightful way. Assuming a given number of phases for the epidemic (namely, 4 in this paper) and a desired tradeoff between health and economic aspects, it is then possible to determine the optimal duration of each phase and the optimal severity level (i.e., the target transmission rate) for each of them. The numerical analysis is performed for the case of France but the adopted approach can be applied to any country. One of the takeaway messages of this analysis is that being able to implement the optimal 4−phase epidemic management strategy in France would have led to 1.05 million of infected people and a GDP loss of 231 billions e instead of 6.88 millions of infected and a loss of 241 billions e. This indicates that, seen from the proposed model perspective, the effectively implemented epidemic management strategy is good economically, whereas substantial improvements might have been obtained in terms of health impact. Our analysis indicates that the lockdown/severe phase should have been more severe but shorter, and the adjustment phase occurred earlier. Due to the natural tendency of people to deviate from the official rules, updating measures every month over the whole epidemic episode seems to be more appropriate.

INTRODUCTION
One of the goals of this work is to provide a simple but exploitable model to measure the quality of the epidemic management strategy implemented by a government to mitigate the health and macro-economic impacts of the Covid-19 epidemic. The quality is measured in terms of the tradeoff between the total number of infected people over a given period of time and the Gross Domestic Product (GDP) loss, under a constraint of the total number of infected people requiring Intensive Care Units (ICU). To reach this objective, we propose a behavioral model for governmental decision-making operations. Although we assume a simple measure for the quality of the lockdown measures and a simple dynamical model (namely, a classical susceptible-exposed-infected-removed (SEIR) model), the proposed approach is seen to be sufficient to constitute a first step into capturing and quantifying the tradeoff under consideration. In contrast with most studies conducted on the Covid-19 epidemic analysis where the primary goal is to refine the SEIR model [see e.g., (1-5)] or employ the SEIR model by accounting for local variations [by using a given SEIR model per geographical region-see e.g., (6)] or for the impact of class type [by age, sex, risk level-see e.g., (7)] our approach is to use the standard SEIR model for an entire country and choose a simple economic model to focus on the study of the tradeoff between health and economic aspects.
Although there have been many several interesting studies on the economic impact of Covid-19 [see e.g., (8)(9)(10)(11)(12)], the pursued goal of these studies is not to model the behavior of the government. As a consequence, the proposed tradeoff has not been analyzed, at least formally. In fact, the closest contribution to this direction would be given by (13) where generic discretetime epidemics over multiple regions are considered, the particular 4-phase structure is not considered, the focus is not on Covid-19, and the key aspect of the tradeoff analysis is neither developed nor analyzed. Additionally, the numerous studies available on the problem of the transmission rate control generally concern the continuous-time control approach. In this work, the focus is on a multiple phase approach (namely, 4 phases). In the literature dedicated to epidemic control, one can for instance find that some recent studies on how the lockdown strategies and quarantine can be planned in an optimal fashion (14)(15)(16)(17). A common feature to all these works on optimal control and lockdown planning is that the policies under consideration, vary over time in a continuous manner, i.e., the lockdown policy is continuously evolving based on the infected population or just on time. However, from the perspective of a government, implementing such policies is not practical since daily changes of the epidemic control measures are difficult to be implemented and to be followed by people.
Summarizing, compared to the existing literature on epidemics modeling and control of epidemics, the main contribution of our work is four-fold: • A model for capturing the tradeoff between health and economic aspect and therefore for the government decisionmaking operation is proposed and studied; • The focus is on multiple phase control policies and not on general continuous-time control policies (to be precise, 4 phases are assumed, see Figure 1); • The problem of finding the optimal features of the optimal epidemic management policy (i.e., the target severity level for each phase and the switching time instants) is stated and solved exhaustively. Additionally, to refine the analysis, we assume a simple model for the natural time drift in terms of behavior of people; • The numerical analysis of the tradeoff is dedicated to the Covid-19 epidemic and a case study for France.

Epidemic Model
To model the dynamics of the Covid-19 epidemic globally i.e., over an entire country, we assume a standard SEIR model. Let us, respectively, denote by s, e, i, and r the fractions of the population: being susceptible to be infected by the SARS-Cov2 virus, having been exposed to it, being infected, and being removed (including recoveries and deceases). The epidemic is assumed to obey the following continuous-time dynamics: where: • β(t), t ∈ R, represents the time-varying virus transmission rate; • γ denotes the rate at which the exposed subject develops the disease (this includes people presenting symptoms and asymptomatics). The period 1 γ is called the incubation period; • δ denotes the removal rate and 1 δ is called the average recovery period.
We assume that the control action u(t) taken by the decisionmaker (the government or possibly a more local decision-maker) has a linear effect on the transmission virus rate. Additionally, the effectiveness of this action is assumed to undergo a noncontrollable drift or attenuation effect due to the observed fact that people tend to relax their effort over time (18,19), hence the presence of the attenuation factor a(t) yields: where R 0 is the natural reproduction number (namely, without any control or population awareness), u(t) ∈ [0, U] is the control action or severity level of the lockdown measures. Note that U corresponds to the most drastic or severe control action (in theory it could reach the value R 0 δ and make the reproduction number vanishing). In this work, u(t) is a piecewise-constant function. For the numerical analysis, we will assume a(t) to be a linearly decreasing function of time (as detailed in the next FIGURE 1 | One of the goals of this work is to determine numerically, for a given tradeoff between health and economic costs, the best 4−phase epidemic management policy that is, the best values for τ 0 , τ 1 , τ 2 , R 1 , R 2 , R 3 (the epidemic time horizon T and the natural reproduction number R 0 being fixed). section). Therefore, one can define the time-varying effective reproduction number: As illustrated by Figure 1, we are solving an epidemic control problem in which determining the function u(t) or R(t) amounts to jointly determining the switching instants τ 0 , τ 1 ,τ 2 and the targeted reproduction numbers R 1 , R 2 , R 3 ; T is a given period of time for the epidemic analysis. In particular, we will determine the best duration of the lockdown phase τ 1 and the corresponding targeted reproduction number R 1 . Figure 1 shows for instance that if the lockdown measures taken are such the reproduction number is R 1 at the beginning of the lockdown phase, then, because of the drift induced by the typical human behavior, the effective reproduction number increases over time.

Time Drift or People Behavior Model
We propose here a model for the attenuation function a(t), which quantifies the degree to which people relax their effort to implement the government management measures. As the attenuation effect is negligible when a new policy is released, we consider that a(t) = 1 when t ∈ {τ 0 , τ 1 , τ 2 }. The attenuation factor is assumed to increase over time in each phase, and we assume the following piecewise linear behavior between phases: where a 1 , a 2 , and a 3 , respectively, represent the attenuation coefficients during the lockdown phase, after the lockdown phase, and during the adjustment phase.

Decision-Maker Behavior Model
The proposed model for the behavior of the decision-maker is based on the fact that it wants to obtain a given tradeoff between economic and health aspects. For the cost related to the economics loss, we assume the simplest reasonable model. That is, we assume that economic cost is quadratic in the control action. For the health cost, we assume that it is given by the number of infected people over the given period of time. Therefore, the proposed overall cost consists of a convex combination of these two costs. By minimizing the overall cost, one realizes the desired tradeoff between economic and health aspects. On top of this we impose the number of patients requiring intensive care to be under a given threshold N ICU max . Thus, the corresponding minimization is performed under a constraint on the number of people infected at any time t ∈ [0, T]: σ Ni(t) ≤ N ICU max , N being the population size, N ICU max the maximum number of ICU patients, and 0 ≤ σ ≤ 1 is the percentage of infected people requiring intensive care. In France, official records state that the maximum cumulated number of ICU patients has reached 7, 148 (on April 8, 2020) but the capacity over the whole territory has been evaluated to be greater than 15, 000. By denoting α ∈ [0, 1] the weight assigned to the macroeconomic impact of the epidemic and K e > 0, K h > 0, µ > 0 some constants (parameters) defined below, obtaining the desired tradeoff amounts to finding a solution of the following optimization problem (OP) while fixing α to a given value: where: • K e > 0 and K h > 0 are constants that weight the economic and health cost functions (they also act as conversion factors allowing one to obtain appropriate units and orders of magnitude); • τ 0 and τ 1 represent the lockdown starting time and duration, respectively. T min is the minimum lockdown duration to make the lockdown policies effective. The quantity τ 2 represents the duration of the post-lockdown phase; • the parameters µ 1 , µ 2 ≥ 1 accounts for possible differences in terms of economic impact between the lockdown and postlockdown phases; • s(0) and s(T) are, respectively, the fractions of the population susceptible at the beginning and the end of the analysis.
We would like to make additional comments concerning the parameter µ 1 , µ 2 . The motivation for introducing µ 1 , µ 2 is twofold. First, after lockdown, people are more aware and act more responsibly than before lockdown. This means that automatic and costless population distancing typically occurs (8,20,21). Taking µ 1 , µ 2 ≥ 1 precisely amounts to having a smaller reproduction number without any cost for the government. Additionally, as people typically tend to relax their effort to implement the epidemic management measures as time passes, it makes sense to assume that µ 1 ≥ µ 2 in our model. It also allows FIGURE 2 | Comparison between our model (blue curve) and the reported statistics (black curve) from March 1st to September 30th. When the number of tests are sufficiently large, our model matches well with the reported statistics. The reproduction number is linearly increasing during each phase, and is discontinuous during the transition between phases.
Frontiers in Public Health | www.frontiersin.org one to account for the fact that, after lockdown, the economic activity grows after the lockdown and the effects of the pandemic starts vanishing. This means that, in some sense we ignore memory effects due to lockdown measures. Further refinements of the proposed model might be considered to account for the lockdown memory effects. This is out of the scope of the present paper but we believe that, this would correspond to assuming µ i < 1.

4-Phase Optimal Control With Piecewise Constant Control Actions
Solving analytically the optimization problem given by (5) is not trivial. However, since we restrict our attention to a certain class of control policies, the problem turns out to be solvable through exhaustive search. Assuming the attenuation factor a(τ 0 ) = a(τ 1 ) = a(τ 2 ) = 1 (no attenuation at the beginning of each phase) and a constant control action in each phase, by using the relation , the OP (5) can be rewritten under a more convenient form for numerical purposes: (1) and (2). (6) where R i represents the desired or target reproduction number over Phase i ∈ {1, 2, 3} without considering the attenuation factor (also, it is the reproduction number at the start of i-th phase). The second constraint is introduced here as there is a gap between lockdown reproduction number and after lockdown reproduction number. Finally, the conversion factors K e and K h FIGURE 3 | The plots represent the possible tradeoffs between health cost (measured in terms of the total number of infected people) and economic cost (measured in terms of GDP loss) that can be obtained (by choosing the best epidemic management policy). In particular, with the assumed model, it is seen, in retrospect, that it would have been possible to divide the number of infected people by about 6 while saving about 10 billions e in terms of GDP with the optimal 4-phase strategy.
Frontiers in Public Health | www.frontiersin.org FIGURE 4 | When economic losses are assumed to be uniform over time (µ = 1), there is some economic incentive to delay the lockdown, but this delay is seen to be at maximum 4 days. When economic losses are lower after lockdown than during it, it is beneficial to start the lockdown faster, up to 3 days delay.
are chosen as follows. The rationale behind the choice of K e is that when choosing α = 1 the GDP loss should correspond to the best estimations made by economists. The GDP loss over the lockdown period for a given country is denoted by GDP, the conversion factor K e is chosen as follows: For France for example, the GDP loss during the lockdown has been evaluated (on April 20) to be around 120 billions e according to the OFCE (22). At last, the constant K h is merely chosen as K h = N, that is, when α = 0 the cost function corresponds to the number of people infected over the considered period of time.

General Simulation Setup
To perform exhaustive search over the sextuple of variables (τ 0 , τ 1 ,  Figure 5, due to the physical characteristics of the epidemics in France, we set T min = 30 because the lockdown duration is at least 4 weeks or 1 month to make the lockdown effective in real systems (23, 24) 1 . The SEIR model parameters are as follows: 1 δ = 1 0.1857 = 5.4 days, 1 γ = 1 0.16 = 6.25 days; these choices are consistent with many works and in particular the studies performed for France (25) and Italy (16). The population size is set to N = 66.10 6 , the maximum number of patients requiring intensive care is set to N ICU max = 15.10 3 , and σ = 1.5% (23, 24) 1 . Notice that this number is only reached for very small values of α (for which the total number of people infected over the analysis duration would be around 9 millions). The exposed population on March 1 2020 is initialized to Ne(0) = 1.33.10 5 . This number is obtained from analyzing the data provided in (23). We consider the number of reported deaths at a given time to be a more reliable way of tracking the evolution of the pandemic rather than the reported number of infected people. Indeed, as soon as one examines absolute values, they typically become irrelevant. For example, during lockdown, because of the lack of tests and measurements, the real number of infected people was much higher than the official number. Now, even when tests were performed intensively, because of false positives, the absolute number of infected were again completely unreliable. For a prevalence of 1/1, 000 and a test reliability of 5% of false positives we see that the number of infected for 1,000 people is declared to be about 50 whereas the actual number of infected is only 1. Motivated by these critical issues, we have considered figures which are much more reliable such as the number of deceases due to Covid-19 in France (23,25). From the number of deaths and the global average rate (worldwide) of the mortality rate (in the range 0.3-0.5% when averaged over classes of ages and countries), the reconstructed number of infected people turns out to be more accurate. Therefore, by fixing the mortality rate to a given value [in (25) for instance, the mortality rate averaged over all the classes of infected people is evaluated to be around 0.53% for France], one can estimate the exposed and infected population size 3-4 weeks before the measured number of deaths due to Covid-19. For our computation, the initial conditions of the ordinary differential equations (ODE) equations are chosen as s(0) = 1 − e(0), i(0) = r(0) = 0, and the ODE is solved by using the Matlab ode45 solver. Concerning the economic cost for France related to Covid-19, the GDP loss over the lockdown period is estimated by the OFCE (22) to be 120 billions e and we have, as reliable figures, that τ France 1 = 55 days with R France 0 = 3.5 FIGURE 5 | When economic losses are assumed to be uniform over time (µ = 1), it is optimal to have a long lockdown (typically between 60 and 80 days) whatever the tradeoff desired. But, if the economic losses of the third phase are less than during lockdown, the lockdown phase should be shorter.
Frontiers in Public Health | www.frontiersin.org and R France 1 = 0.6. We therefore take K e = 7.379.10 9 e/day. Values of the reproduction number for the first two phases come from past and quite accurate evaluations [see e.g., (24)]. The value R France 2 = 0.9 is less accurate and corresponds to the assumption that the government has been aiming at giving as much as freedom to the population while avoiding a second wave. Based on available statistics on Covid-19 in France 1 , the attenuation coefficients of the drift model have been chosen as follows: a 1 = 0.1%, a 2 = 0.2%, a 3 = 0.2%. Unless stated otherwise, the economic impact parameters are chosen as µ 1 = 1.41 (i.e., µ 2 1 ∼ 2) and µ 2 = 1.3. Also, when α is assumed to be fixed, it is set to 10 −4 .
To justify the choice of the attenuation parameters, that is, a 1 = 0.1%, a 2 = 0.2%, a 3 = 0.2%, we apply the French policy into our model. By comparing the active cases obtained from our model and the official data, it can be illustrated in Figure 2 that our model matches well with the real data, especially in the second-half of the plot, where the number of tests conducted are sufficiently large. This validate the choice of our model and parameters. The mismatch on March and April are mainly due to the lacking number of tests that were taken during the early outbreak of the pandemic, leading to a much lower reported number of active cases.

Optimal Tradeoff Between Economic and Health Impacts
With the proposed government decision-making model, implementing a desired tradeoff between the health cost and economic cost merely amounts to choosing a given value for the parameter α. Figure 3 depicts for various values of α in the interval [10 −7 , 10 −4 ] the total GDP loss and number of infected people that is obtained after choosing the (quantized version of the) sextuple (τ 0 , τ 1 , τ 2 , R 1 , R 2 , R 3 ) that minimizes the combined cost given by Equation (6). At one extreme, when α is relatively large (α = 10 −4 ) (that is, when the government aims at minimizing the economic cost in the first place-always under the ICU capacity constraint) we see that the best epidemic management strategy leads to a GDP loss over the entire study period  30; with our model, the GDP loss over the period of time of interest is 241 billion e and the total number of infected people is about 6.88 million. What the best tradeoff curve indicates is that there were management policies that would allow the French government to have a better "performance" both in terms of GDP loss and the number of infected people. For instance, we indicate a point for which it would have been possible to have about 1.05 million people infected (that is, about 6 times less than what is estimated with the current policy) while ensuring a total GDP loss of 231 billions e. Which type of epidemic management strategy should be used to have such an outcome? The next sections provide a detailed analysis of the features of the optimal strategy.

Optimal Features of the Optimal Epidemic Management Strategy
One of the important features for controlling the Covid-19 epidemic which has been well-commented in newspapers in various countries is the lockdown starting time. To minimize the health cost, the answer is ready: the lockdown phase should always start as soon as possible. But when one wants to realize a tradeoff between health and economic costs, the answer is less immediate. For different values for µ 1 and µ 2 , Figure 4 provides the best day to start locking down the population, for one hundred values of α ranging from 10 −4 to 10 −6 . The main message of this figure is that, even for (relatively) large values for α (that is, when the economic cost dominates the health cost), the optimal lockdown starting day should be before March 4th (i.e., τ 0 ≤ 4). This clearly shows that, once an epidemic has been declared, invoking economic damages to delay the lockdown phase is not acceptable. Note that this conclusion holds when economic losses are assumed to be uniform over time (µ 1 = µ 2 = 1). When the economic cost associated with a given intensity or severity level is lower after lockdown than during it (here µ 1 = 1.41 and µ 2 = 1.3), it is always optimal to start locking down as soon as possible. Note that our model does not capture the possible fact that population needs to be psychologically prepared to follow the lockdown measures. In France, by March 17, there were official figures about the epidemic which were sufficiently critical to make the population accept the measures whereas, starting on March 4 (the optimal starting date for µ 1 = µ 2 = 1) the situation might have not been critical enough to create full adhesion to government measures.
A second key feature of the Covid-19 epidemic control strategy was the lockdown phase duration namely, the value of τ 1 . To better explore the relationship between τ 1 and α (the tradeoff), we relax the lockdown duration constraint here and set T min = 1. For (µ 1 , µ 2 ) = (1, 1) and (µ 1 , µ 2 ) = (1.41, 1.3), FIGURE 7 | Comparison of lockdown policies with different R 0 . When the situation is worse, a sooner and more strict lockdown strategy should be applied.
Frontiers in Public Health | www.frontiersin.org Figure 5 provides the optimal lockdown duration (in days) for values of α ranging from 10 −4 to 10 −6 . For (µ 1 , µ 2 ) = (1, 1) (i.e., when economic losses are uniform over time), the optimal duration ranges from 53 to 83 days for a large fraction of the considered interval for α. Interestingly, we see that these values are relatively close to the lockdown duration effectively imposed in France namely 55 days. For larger values of µ 1 and µ 2 , the optimal lockdown duration is seen to be much smaller. Therefore, if economic losses are uniform over time, the French government policy seems to be very coherent. On the other hand, if the economic impact is smaller after lockdown, our study suggests shorter lockdown durations. In fact, our results show the existence of a critical value for the tradeoff parameter α above which the second phase of the management of the epidemic should not be present. This means that the optimal control consists of three phases instead of four.
To conclude this section, let us consider Figure 6. For the by default scenario studied in this paper [(µ 1 , µ 2 ) = (1.41, 1.3), α = 10 −4 ], the figure represents the evolution of number of infected people, that is, Ni(t), and the transmission rate when the optimal policy is adopted. First, it is seen that for the healtheconomic tradeoff corresponding to α = 10 −4 , there is no interest in delaying the lockdown phase. By "lockdown" phase, the authors mean that it might be any type of phase for which the reproduction number is as low as 0.4 (vs. the estimated 0.6 in France); very efficient digital tracing and intensive use of face masks is also an option which has been successfully adopted in countries such as South Korea [see e.g., (26)]. The optimal lockdown phase duration is seen to be about 1 month (instead of 2 for France). We see that the existence of an adjustment phase is part of the optimal policy. For this point, we see that the adjustment phase should have occurred much before in France (end of June vs. end of September). The next section is precisely dedicated to the impact of the adjustment phase.

Lockdown Policies With Different R 0
Since the natural reproduction number R 0 depends on temperatures, population density, may vary over time due to mutation effects, and in any case is not known perfectly, it is of interest to study the impact of R 0 on the obtained characteristics for the optimal epidemic management policy. This is what Figure 7 represents. It is seen that large variations on R 0 do not involve large variations on the starting day. For instance, moving from R 0 = 2 to R 0 = 3.5 only changes the optimal date by 1 day (Day 2 instead of Day 3), which confirms the need to act fastly even when the transmission is more limited (e.g., thanks to higher temperatures or lower population density). Note that this holds even if the economic impact is accounted for. It is Frontiers in Public Health | www.frontiersin.org also good for economical aspects to react fastly to an epidemics. For the optimal reproduction number during lockdown (namely, R 1 ) it is also seen that moving a scenario in which R 0 = 2 to R 0 = 3.5 does not change very significantly the results: the target severity (or freedom) level would correspond to R 1 = 0.4 instead R 1 = 0.6, which shows that the severity should be high during lockdown even in countries or regions where transmission is more limited.
Alternatively, the impact of R 0 uncertainty can be evaluated by adding a perturbation on R 0 . It is assumed that what is known to determine the optimal parameters is R 0 = R 0 + , being a Gaussian noise ∼ N (0, σ 2 ). The reproduction number have to stay in a given interval of physical relevance of the form [R min , R max ]. Thus, the noise is imposed to stay in the interval [−R 0 + R min , −R 0 + R max ]. With R 0 = 3.5, R min = 1 and R max = 4, Figure 8 depicts the average bias for τ 0 and R 1 induced by uncertainty on R 0 . The average biases for τ 0 and R 1 are where R 1 and τ 0 are obtained with the noisy reproduction number R 0 . Remarkably, the impact of the corresponding noise on the results is seen to be very reasonable and does not affect the main conclusions drawn in the first version of the paper. This indicates that the conducted analysis is robust against some forms of uncertainties. But of course, as mentioned previously, a deeper analysis would be required to state more general conclusions.
In addition, a more global sensitivity analysis to other SEIR model parameters, i.e, δ and γ , have been conducted by simulations. A similar conclusion can be drawn. This shows that our conclusions are quite robust to effects such as parameter uncertainties. This is in part due to the fact that we are mostly seeking discrete parameters and not continuous parameters, which creates a form of robustness.

Impact of the Adjustment Phase
Always for the typical scenario presented in the general setting part, Figure 9 depicts the evolution of the number of infected people (in France) for the policy effectively implemented over the period [March 1, September 30]. Here, only the adjustment phase is assumed to be optimizable. The figure allows one to quantify the impact of the severity level of the adjustment phase. Without the adjustment phase, the fraction of infected people is such that the number of people requiring intensive care exceeds the double of the maximum ICU capacity of France. However, by implementing measures such that R(t) < 1.2 over the adjustment phase, the constraint on the ICU capacity is not violated by the end of 2020. Furthermore, to avoid the overwhelming health service for a longer time, it is better to implement measures such that R(t) < 1 over the adjustment phase since the fraction of infected people will be non-increasing with R(t) < 1.

DISCUSSION
In this work, we propose to model the behavior of a government as far as the epidemic control is concerned. The proposed model, despite its simplicity, has the merit of being able to capture the fundamental tradeoff between economic and health aspects. Obviously, to capture other effects such as the psychological effects of measures on people, a more general model should be considered. The proposed model allows one to provide quantitative answers to issues which have been well-commented in the media. For example, even when a government chooses to assign a high value to the economic aspect, it is seen that the best strategy is almost always to implement a severe phase as soon as possible. This severe phase involves locking the population down, as most countries did, or to make intensive use of digital tracing and face masks as South Korea did. In the latter case, a loss in terms of privacy is the price to be paid for having more freedom. When inspecting the obtained numerical results performed for France, it is seen that the optimal features for the lockdown/severe phase require targeting a reproduction number smaller than the one achieved in France (0.4 vs. 0.6), while having a shorter duration for the severe phase (1 month instead of 2). Note that some countries adopted measures which were more severe than France. For instance, China has been imposing the use of a given food supply system which was very efficient in terms of epidemic mitigation (27)(28)(29). Then, by planning an adjustment phase at the right time, we have seen that as a final result, the number of infected people can be reduced by a factor of 6 when compared to the current French policy (1.05 million of infected instead of 6.88 millions for the current policy) while having similar GDP losses (231 billions e instead of 241 billions e). We have seen that by considering a simple model as we have studied, the need for an adjustment phase could have been anticipated. Such a phase is necessary to avoid the number of patients under intensive care exceeding the capacity of the ICUs. Also, because of the natural tendency of humans to deviate from rules over time, it appears that measures should be updated about every month and not less frequently.

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 author/s.

AUTHOR CONTRIBUTIONS
SL, CZ, VV, and IM designed the proposed model. SL and CZ conducted the review of the related literature. SL collected the data to calibrate the assumed model. CZ and SL performed the simulations. All authors contributed to the writing of the manuscript.

FUNDING
This work was partially supported by ANR via the grant NICETWEET, number ANR-20-CE48-0009.