# Expected Evolution of COVID-19 Epidemic in France for Several Combinations of Vaccination Strategies and Barrier Measures

Simon Pageaud, Catherine Pothier, Christophe Rigotti, Anne Eyraud-Loisel, Jean-Pierre Bertoglio, Alexis Bienvenüe, Nicolas Leboisne, Nicolas Ponthus, Romain Gauchon, François Gueyffier, et al.

► To cite this version:

Simon Pageaud, Catherine Pothier, Christophe Rigotti, Anne Eyraud-Loisel, Jean-Pierre Bertoglio, et al. Expected Evolution of COVID-19 Epidemic in France for Several Combinations of Vaccination Strategies and Barrier Measures. Vaccines, 2021, 9 (12), pp.1462. 10.3390/vaccines9121462. hal-03482342

HAL Id: hal-03482342

https://hal.science/hal-03482342

Submitted on 6 Dec 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Expected Evolution of COVID-19 Epidemic in France for Several Combinations of Vaccination Strategies and Barrier Measures

Simon Pageaud<sup>1,2,3,4,5,6</sup>, Catherine Pothier<sup>7,8 ID</sup>, Christophe Rigotti<sup>7,8,9</sup>, Anne Eyraud-Loisel<sup>5</sup>, Jean-Pierre Bertoglio<sup>10,11</sup>, Alexis Bienvenüe<sup>5 ID</sup>, Nicolas Leboisne<sup>5</sup>, Nicolas Ponthus<sup>11,12,13</sup>, Romain Gauchon<sup>5</sup>, François Gueyffier<sup>1,2,3,4</sup>, Philippe Vanhems<sup>14,15</sup>, Jean Iwaz<sup>1,2,3,4</sup>, Stéphane Loisel<sup>5</sup>, Pascal Roy<sup>1,2,3,4,12,* ID</sup> and on behalf of the Group CovDyn (Covid Dynamics) †

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Pageaud, S.; Pothier, C.; Rigotti, C.; Eyraud-Loisel, A.; Bertoglio, J.-P.; Bienvenüe, A.; Leboisne, N.; Ponthus, N.; Gauchon, R.; Gueyffier, F.; et al. Expected Evolution of COVID-19 Epidemic in France for Several Combinations of Vaccination Strategies and Barrier Measures. Vaccines **2021**, *9*, 1462. https://doi.org/10.3390/vaccines9121462

Academic Editor: S. Louise Cosby

Received: 20 October 2021
Accepted: 2 December 2021
Published: 10 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

1 Université de Lyon, F-69000 Lyon, France; simon.pageaud@gmail.com (S.P.); francois.gueyffier@univ-lyon1.fr (F.G.); jean.iwaz@chu-lyon.fr (J.I.)

2 Université Claude Bernard Lyon 1, F-69100 Villeurbanne, France

3 CNRS UMR 5558, Laboratoire de Biométrie et Biologie Évolutive, F-69100 Villeurbanne, France

4 Service de Biostatistique-Bioinformatique, Pôle Santé Publique, Hospices Civils de Lyon, F-69003 Lyon, France

5 Laboratoire de Sciences Actuarielle et Financière (LSAF), Institut de Science Financière et d'Assurances (ISFA), Université Claude Bernard Lyon 1, F-69007 Lyon, France; anne.loisel@univ-lyon1.fr (A.E.-L.); alexis.bienvenue@univ-lyon1.fr (A.B.); nicolas.leboisne@univ-lyon1.fr (N.L.); romain.gauchon@laposte.net (R.G.); stephane.loisel@univ-lyon1.fr (S.L.)

6 Fondation du Risque, Groupe Louis Bachelier, F-75002 Paris, France

7 CNRS UMR 5205, Laboratoire d'InfoRmatique en Image et Systèmes d'Information (LIRIS), F-69621 Villeurbanne, France; catherine.pothier@insa-lyon.fr (C.P.); christophe.rigotti@insa-lyon.fr (C.R.)

8 Institut National des Sciences Appliquées de Lyon (INSA), F-69621 Villeurbanne, France

9 INRIA Grenoble-Rhône-Alpes, F-38334 Montbonnot, France

10 CNRS UMR 5509, Laboratoire de Mécanique des Fluides et d'Acoustique (LMFA), F-69130 Ecully, France; jean-pierre.bertoglio@ec-lyon.fr

11 École Centrale de Lyon, F-69130 Lyon, France; nicolas.ponthus@gmx.fr

12 CNRS UMR 5513, Laboratoire de Tribologie et Dynamique des Systèmes (LTDS), F-69130 Ecully, France

13 École Nationale des Travaux Publics de l'État (ENTPE), F-69120 Vaulx-en-Velin, France

14 Service d'Hygiène, Épidémiologie, Infectiovigilance et Prévention, Hôpital Edouard Herriot, Hospices Civils de Lyon, F-69003 Lyon, France; philippe.vanhems@chu-lyon.fr

15 Centre International de Recherche en Infectiologie (CIRI: Inserm U1111, CNRS UMR 5308, École Nationale Supérieure de Lyon), F-69007 Lyon, France

* Correspondence: pascal.roy@chu-lyon.fr

† Members of the group are listed in Appendix A.

**Abstract:** The outbreak of the SARS-CoV-2 virus, enhanced by rapid spreads of variants, has caused a major international health crisis, with serious public health and economic consequences. An agent-based model was designed to simulate the evolution of the epidemic in France over 2021 and the first six months of 2022. The study compares the efficiencies of four theoretical vaccination campaigns (over 6, 9, 12, and 18 months), combined with various non-pharmaceutical interventions. In France, with the emergence of the Alpha variant, without vaccination and despite strict barrier measures, more than 600,000 deaths would be observed. An efficient vaccination campaign (i.e., total coverage of the French population) over six months would divide the death toll by 10. A vaccination campaign of 12, instead of 6, months would slightly increase the disease-related mortality (+6%) but require a 77% increase in ICU bed-days. A campaign over 18 months would increase the disease-related mortality by 17% and require a 244% increase in ICU bed-days. Thus, it seems mandatory to vaccinate the highest possible percentage of the population within 12, or better yet, 9 months. The race against the epidemic and virus variants is really a matter of vaccination strategy.

**Keywords:** vaccination; COVID-19; agent-based model; decision support techniques
# 1. Introduction

The outbreak of SARS-CoV-2 virus is causing major national, European, and international health crises, with serious public health and economic consequences. The COVID-19 epidemic started in China, in late autumn 2019, and, on 31 December [1], the WHO China office received reports on severe pneumonia cases in the city of Wuhan (Hubei Province). On 14 February 2021, more than 100 million cases and 2.4 million officially recognized deaths were recorded worldwide, of which there were 34 million cases and 760 thousand deaths in Europe and 3.4 million cases and 81.6 thousand deaths in France. In most countries, the dynamics of the epidemic is continuously assessed through monitoring the number of hospital admissions, occupancy rate of intensive care units (ICUs), and cumulative number of confirmed COVID-19-related deaths.

In France, in 2020, two lockdown periods were imposed to avoid saturation of ICUs that would increase mortality. The first period (17 March to 11 May 2020) aimed at limiting the severity of the first wave. On 8 April 2020, the maximum number of daily hospital deaths (a 7-day moving average) was 532 and, between 6 June and 16 September 2020, less than 30 daily hospital deaths being observed to form a plateau. At the beginning and end of the second lockdown (30 October and 15 December 2020, respectively), the numbers of daily hospital deaths were 222 and 276, respectively, and a maximum of 419 deaths was reached on 19 November. Nevertheless, the resulting epidemic situation was less favorable after the second lockdown; it showed a residual plateau with 240 to 340 daily hospital deaths.

On 14 December 2020, the Alpha variant (lineage B.1.1.7) was identified in the United Kingdom [2], whereas, in France, the first case was identified by the end of December. The World Health Organisation designated that variant as a 'variant of concern' (VOC), which would require new appropriate public health actions. In France, the Alpha variant was deemed responsible for 3.3% of diagnosed cases on 8 January 2021, 14% on 27 January, and 36% on 16 February. Preliminary studies suggested that it increased transmissibility (infectiousness) but not disease severity. The prevalence of this variant in France is available from the public data of Santé Publique France (https://www.santepubliquefrance.fr/, accessed on 15 October 2020). Lineage B.1.672 of SARS-CoV-2 was first detected in December 2020 in India [3], but it appeared only months later in France. The Delta variant (lineage B.1.617.2) was labeled as a VOC on 11 May 2021; on that date, in France, it was deemed responsible for 1% of diagnosed cases.

Vaccines against COVID-19 have been rapidly developed to limit the spread of the epidemic through herd immunity after reaching a high vaccination coverage. Their efficacy against the severe forms of the disease, caused by the historical strain, was estimated to lie between 75% [4] and 94% [5,6]. By the time of those vaccine releases, their efficacy in preventing disease transmission was not known.

In France, protective measures have been constantly adapted to the dynamics of the epidemic. Between September 2020 and January 2021, three periods with different levels of non-pharmaceutical interventions (NPIs) took place. The first period (September 2020–October 2020) corresponded to moderate-NPIs, such as mask wearing and social distancing, but restaurants and leisure facilities remained open. The second period (November 2020) was characterized by extended-NPIs, which consisted of very strict measures, such as lockdown during weekends or full lockdown, as well as closing of schools, universities, and shopping facilities. The third period (December 2020–January 2021) was characterized by intensive-NPIs (intermediate between moderate-NPIs and extended-NPIs), combining several strict protective measures, such as overnight curfew, social distancing, mask-wearing, and increased distancing at work, as well as the closing of bars, restaurants, leisure facilities, etc.

This study aimed to predict the course of the COVID-19 epidemic over an 18-month period, starting on 1 January 2021, by applying several vaccination campaign and barrier measure combinations and taking into account the increasing proportion of the Alpha variant. It was based on a previously published scientific report on a non-calibrated
theoretical version [7]. (Note: the proposed displayed results were obtained before the emergence of the Delta variant).

## 2. Results

### 2.1. *Dynamics of the Epidemic with the Historical Strain without Vaccination*

Table 1 presents the results obtained by restricting the simulations to the historical strain without vaccination and setting the prevalence of asymptomatic cases at 20%.

<table><caption>**Table 1.** Historical strain without vaccination. Cumulative numbers of individuals removed and hospital deaths on 30 June 2022.</caption><thead><tr><td>NPIs</td><td>Removed</td><td>Deceased</td></tr></thead><tbody><tr><td>Relaxed-NPIs</td><td>26,397,570</td><td>372,973</td></tr><tr><td>Intensive-NPIs</td><td>16,688,896</td><td>202,411</td></tr><tr><td>Extended-NPIs</td><td>7,906,676</td><td>60,771</td></tr></tbody></table>

After 18 months of relaxed-NPIs (systematic alternation of periods of 45 days between intensive-NPIs and moderate-NPIs), the simulated number of removed individuals was 26.4 millions, and the number of deaths was 373 thousands. In this scenario, the cumulative distribution of deaths plateaued after 500 days (Figure 1a). A maximum daily requirement of 23,732 ICU beds were needed on day 186 (Figure 1b). A slightly delayed effect of changes between intensive and moderate protection measures on ICU needs was observed under relaxed-NPIs, and the control of the epidemic was obtained after 500 days.

After 18 months of intensive-NPIs, the simulated number of removed individuals was much lower (around 16.7 millions), with 202.4 thousand deaths. The cumulative distribution of deaths was still increasing after 500 days (Figure 1a). A maximum daily requirement of 5960 ICU beds were needed on day 163 (Figure 1b). A slow decrease was then observed, indicating a control of the epidemic after 500 days of intensive-NPIs.

After 18 months of extended-NPIs, the simulated number of removed individuals was 7.9 millions, with 60.8 thousand deaths. The cumulative distribution of deaths reached a plateau after 150 days (Figure 1a). A maximum daily requirement of 2700 ICU beds were needed on day 6 (Figure 1b). A rapid decrease was observed, indicating a control of the epidemic by day 150 under extended-NPIs.

Figure 1 shows clearly distinct changes in death (Figure 1a) and ICU beds occupation (Figure 1b), when relaxed-NPIs, intensive-NPIs, and extended-NPIs were applied, respectively.

### 2.2. *Dynamics of the Epidemic with the Introduction of the Alpha Variant without Vaccination*

Table 2 presents the results obtained in a situation of competition between the historical strain and the more infectious Alpha variant in the absence of vaccination. With relaxed-NPIs, the simulated number of removed individuals was around 48.3 millions, with 621.3 thousand deaths. The cumulative distribution of deaths reached a plateau at day 300 (Figure 2a). A maximum daily requirement of 125,824 ICU beds were needed on day 125. A huge single wave occurred, without possible distinction between intensive-NPIs and moderate-NPIs periods of 45 days. With relaxed-NPIs, a control of the epidemic occurred at day 350 (Figure 2b).

With intensive-NPIs, the simulated number of removed individuals was still important (around 45.5 millions), with 570.6 thousand deaths. The cumulative distribution of deaths reached a plateau after 300 days (Figure 2a). A maximum daily requirement of 116,548 ICU beds was needed on day 147. Despite continuous intensive-NPIs, a huge single wave was observed, and the control of the epidemic occurred at day 350 (Figure 2b).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Expected evolution of the COVID-19 epidemic in France for the historical strain, without vaccination, for extended-NPIs, intensive-NPIs, and relaxed-NPIs, as well as for 20% rate of asymptomatic subjects. Day 0 is 1 January 2021; (a) cumulative number of deaths. (b) Prevalence of COVID-19 patients needing intensive care units. The horizontal black lines represent ICU capacities of 5000, 8000, and 12,000 beds.</figcaption></figure>

With extended-NPIs, the simulated number of removed individuals was 20.6 millions, with 213.7 thousand deaths. The cumulative distribution of deaths still increased at 500 days (Figure 2a). A maximum daily requirement of 12,650 ICU beds was needed on day 343. In comparison to the wave simulated under intensive-NPIs, the wave under extended-NPIs was delayed, and the control of the epidemic did not occur at day 500 (Figure 2b).
With relaxed- and intensive-NPIs, close simulated results were obtained, in terms of mortality and ICU overload. Only extended-NPIs seemed able to partially limit the number of deaths without ensuring herd immunity (Figure 2).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Expected evolution of the COVID-19 epidemic in France with the historical strain and Alpha variant, without vaccination, for extended-NPIs, intensive-NPIs, and relaxed-NPIs, as well as for 20% rate of asymptomatic subjects. Day 0 is 1 January 2021; (a) cumulative number of deaths. (b) Prevalence of COVID-19 patients needing Intense Care Units. The horizontal black lines represent an ICU capacity of 5000, 8000, and 12,000 beds.</figcaption></figure>
<table><caption>**Table 2.** Historical strain and Alpha variant development without vaccination. Cumulative numbers of individuals removed and hospital deaths on 30 June 2022.</caption><thead><tr><td>NPIs</td><td>Removed</td><td>Deceased</td></tr></thead><tbody><tr><td>Relaxed-NPIs</td><td>48,293,980</td><td>621,289</td></tr><tr><td>Intensive-NPIs</td><td>45,537,214</td><td>570,555</td></tr><tr><td>Extended-NPIs</td><td>20,551,828</td><td>213,715</td></tr></tbody></table>

## 2.3. Effects of a 6-Month Vaccination Campaign with Intensive-NPIs

Under the assumption of a 90% reduction in virus transmission, due to the vaccine, the number of individuals aged 10 years and older removed by the disease (8.6 millions) was divided by 5, and the number of deaths (56.9 thousands) was divided by 10 (Table 3). The number of deaths reached a plateau at day 200 (Figure 3).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Expected evolution of the COVID-19 epidemic in France for the historical strain and Alpha variant, for different durations of vaccination campaigns and reductions of virus transmission of 50%, 75%, and 90%. Each row displays the prevalence of COVID-19 patients needing intensive care and cumulative number of deaths at hospital for 20% rate of asymptomatic subjects. Day 0 is 1 January 2021.</figcaption></figure>
<table><caption>**Table 3.** Results of a vaccine with 90% transmission reduction, with the historical strain and Alpha variant under different vaccination ages, campaign duration, and non-pharmaceutical interventions (NPIs). Cumulative numbers of individuals removed and deceased on 30 June 2022, cumulative ICU bed-days required, maximum daily ICU bed-days required, and intensive care unit (ICU) overload (exceeding bed-days capacity) over a period of 18 months.</caption><thead><tr><th rowspan="2">NPIs</th><th rowspan="2">Vaccination Ages</th><th rowspan="2">Campaign Duration (Days)</th><th rowspan="2">Removed</th><th rowspan="2">Deceased</th><th rowspan="2">Cumulative ICU Bed-Days</th><th rowspan="2">Max Bed-Days /On Day</th><th colspan="3">Days of ICU Overload If Max Available Beds Is</th></tr><tr><th>1000</th><th>2000</th><th>3000</th></tr></thead><tbody><tr><td>Intensive-NPIs</td><td>10+</td><td>180</td><td>8,625,521</td><td>56,916</td><td>192,728</td><td>2926/30</td><td>69</td><td>48</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>270</td><td>9,446,952</td><td>59,006</td><td>253,715</td><td>2928/30</td><td>86</td><td>57</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>360</td><td>10,715,225</td><td>60,552</td><td>341,083</td><td>3001/27</td><td>109</td><td>69</td><td>1</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>540</td><td>14,160,343</td><td>66,375</td><td>662,125</td><td>3281/80</td><td>270</td><td>155</td><td>72</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>180</td><td>8,537,674</td><td>56,466</td><td>180,477</td><td>2853/26</td><td>65</td><td>46</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>270</td><td>9,190,934</td><td>57,836</td><td>230,283</td><td>2855/30</td><td>79</td><td>52</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>360</td><td>10,074,304</td><td>58,819</td><td>304,086</td><td>2883/29</td><td>100</td><td>61</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>540</td><td>13,321,101</td><td>65,539</td><td>556,832</td><td>3250/36</td><td>256</td><td>107</td><td>27</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>180</td><td>7,681,105</td><td>53,277</td><td>101,474</td><td>2553/3</td><td>42</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>270</td><td>7,702,427</td><td>53,331</td><td>107,054</td><td>2461/7</td><td>44</td><td>22</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>360</td><td>7,725,029</td><td>54,139</td><td>118,005</td><td>2634/7</td><td>47</td><td>26</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>540</td><td>7,816,371</td><td>55,280</td><td>122,126</td><td>2578/4</td><td>47</td><td>25</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>180</td><td>7,695,726</td><td>53,418</td><td>106,390</td><td>2592/3</td><td>44</td><td>24</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>270</td><td>7,785,962</td><td>54,464</td><td>116,136</td><td>2633/7</td><td>47</td><td>26</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>360</td><td>7,734,358</td><td>53,728</td><td>113,134</td><td>2615/7</td><td>45</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>540</td><td>7,749,345</td><td>54,387</td><td>113,459</td><td>2517/7</td><td>46</td><td>22</td><td>0</td></tr></tbody></table>

The ICUs never reached saturation (Figure 3), although ICU bed requirements exceeded 1000 beds for 69 days and 2000 beds for 48 days That need never exceeded 3000 beds. A total of 192,728 ICU bed-days were required, with a maximum daily requirement of 2926 beds on day 30.

In comparison with vaccination of individuals aged 10 years and older, vaccination of people of all ages resulted in an additional 0.8% reduction in the number of deaths.

### 2.4. Effects of a 9-Month Vaccination Campaign with Intensive-NPIs

Under the assumption of a 90% reduction in virus transmission, due to the vaccine, the number of individuals aged 10 years and older removed by the disease (9.4 millions) and number of deaths (59.0 thousands) were quite similar to those obtained with a 6-month vaccination campaign (Table 3). The number of deaths reached a plateau at day 250 (Figure 3).

The pressure on ICUs was increasing (Figure 3); 253,715 ICU bed-days were required.

### 2.5. Effects of a 12-Month Vaccination Campaign with Intensive-NPIs

Under the assumption of a 90% reduction in virus transmission, due to the vaccine, the number of individuals aged 10 years and older removed by the disease (10.7 millions) and the number of deaths (60.6 thousands) were slightly higher than those obtained with 6-month and 9-month vaccination campaigns (Table 3). The number of deaths reached a plateau at day 300 (Figure 3).

The pressure on ICUs was still increasing (Figure 3); 341,083 ICU bed-days were required.

### 2.6. Effects of a 18-Month Vaccination Campaign with Intensive-NPIs

Under the assumption of a 90% reduction in virus transmission due to the vaccine, the number of individuals aged 10 years and older removed by the disease (14.2 millions) and the number of deaths (66.4 thousands) were, respectively, 32% and 10% higher than those obtained with a 12-month vaccination campaigns (Table 3). The number of deaths reached a plateau at day 400 (Figure 3).

The pressure on ICUs was much higher (Figure 3); 662,125 ICU bed-days were required.
## 2.7. Impact of Transmission Reduction on Vaccination Campaign Efficacy

With a vaccine less able to reduce viral transmission, the negative impacts of extending the vaccination campaign on mortality, and on the need for ICU beds, increased (Tables 4 and 5).

With a transmission reduction of 90%, extending the vaccination campaign from 6 to 18 months increased the number of deaths from 56.9 to 66.4 thousands (+17%), and the number of ICU bed-days needed from 192,728 to 662,125 (+244%) (Table 3). With a transmission reduction of 50%, the number of deaths increased from 57.5 to 74.0 thousands (+29%) and the number of ICU bed-days from 216,751 to 1,106,397 (+410%) (Table 5).

<table><caption>**Table 4.** Results of a vaccine with 75% transmission reduction, with the historical strain and Alpha variant under different vaccination ages, campaign duration, and non-pharmaceutical interventions (NPIs). Cumulative numbers of individuals removed and deceased on 30 June 2022, cumulative ICU bed-days required, maximum daily ICU bed-days required, and intensive care unit (ICU) overload (exceeding bed-days capacity), over a period of 18 months.</caption><thead><tr><th rowspan="2">NPIs</th><th rowspan="2">Vaccination Ages</th><th rowspan="2">Campaign Duration (Days)</th><th rowspan="2">Removed</th><th rowspan="2">Deceased</th><th rowspan="2">Cumulative ICU Bed-Days</th><th rowspan="2">Max Bed-Days /On Day</th><th colspan="3">Days of ICU Overload If Max Available Beds Is</th></tr><tr><th>1000</th><th>2000</th><th>3000</th></tr></thead><tbody><tr><td>Intensive-NPIs</td><td>10+</td><td>180</td><td>8,848,757</td><td>57,720</td><td>207,524</td><td>3087/25</td><td>73</td><td>51</td><td>10</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>270</td><td>9,726,806</td><td>58,995</td><td>270,539</td><td>2958/30</td><td>94</td><td>61</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>360</td><td>11,501,532</td><td>61,280</td><td>397,393</td><td>3088/32</td><td>154</td><td>81</td><td>13</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>540</td><td>15,243,733</td><td>67,387</td><td>741,392</td><td>3378/84</td><td>278</td><td>210</td><td>76</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>180</td><td>8,574,155</td><td>55,863</td><td>182,187</td><td>2787/20</td><td>67</td><td>47</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>270</td><td>9,427,345</td><td>58,267</td><td>248,697</td><td>2950/28</td><td>86</td><td>57</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>360</td><td>10,826,562</td><td>60,202</td><td>360,763</td><td>3061/31</td><td>123</td><td>73</td><td>12</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>540</td><td>14,114,090</td><td>65,830</td><td>629,200</td><td>3090/37</td><td>265</td><td>137</td><td>24</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>180</td><td>7,726,318</td><td>53,926</td><td>106,497</td><td>2544/8</td><td>43</td><td>24</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>270</td><td>7,799,063</td><td>54,892</td><td>114,168</td><td>2553/7</td><td>45</td><td>25</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>360</td><td>7,689,234</td><td>53,274</td><td>109,427</td><td>2497/8</td><td>44</td><td>21</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>540</td><td>7,782,863</td><td>54,260</td><td>122,188</td><td>2565/7</td><td>48</td><td>24</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>180</td><td>7,691,260</td><td>52,927</td><td>103,575</td><td>2496/8</td><td>43</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>270</td><td>7,801,188</td><td>54,828</td><td>113,826</td><td>2726/7</td><td>45</td><td>25</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>360</td><td>7,760,065</td><td>54,589</td><td>115,912</td><td>2535/7</td><td>47</td><td>25</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>540</td><td>7,857,610</td><td>56,095</td><td>127,551</td><td>2746/8</td><td>49</td><td>27</td><td>0</td></tr></tbody></table>

<table><caption>**Table 5.** Results of a vaccine with 50% transmission reduction, with the historical strain and Alpha variant under different vaccination ages, campaign duration, and non-pharmaceutical interventions (NPIs). Cumulative numbers of individuals removed and deceased on 30 June 2022, cumulative ICU bed-days required, maximum daily ICU bed-days required, and intensive care unit (ICU) overload (exceeding bed-days capacity) over a period of 18 months.</caption><thead><tr><th rowspan="2">NPIs</th><th rowspan="2">Vaccination Ages</th><th rowspan="2">Campaign Duration (Days)</th><th rowspan="2">Removed</th><th rowspan="2">Deceased</th><th rowspan="2">Cumulative ICU Bed-Days</th><th rowspan="2">Max Bed-Days /On Day</th><th colspan="3">Days of ICU Overload If Max Available Beds Is</th></tr><tr><th>1000</th><th>2000</th><th>3000</th></tr></thead><tbody><tr><td>Intensive-NPIs</td><td>10+</td><td>180</td><td>9,202,882</td><td>57,504</td><td>216,751</td><td>2870/17</td><td>80</td><td>52</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>270</td><td>11,043,587</td><td>59,056</td><td>336,887</td><td>2926/31</td><td>131</td><td>72</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>360</td><td>14,036,762</td><td>63,840</td><td>556,218</td><td>3031/38</td><td>241</td><td>119</td><td>11</td></tr><tr><td>Intensive-NPIs</td><td>10+</td><td>540</td><td>19,133,475</td><td>74,014</td><td>1,106,397</td><td>5219/134</td><td>298</td><td>258</td><td>207</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>180</td><td>9,023,689</td><td>56,264</td><td>202,522</td><td>2830/25</td><td>75</td><td>49</td><td>0</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>270</td><td>10,527,221</td><td>59,584</td><td>310,561</td><td>3071/28</td><td>113</td><td>62</td><td>10</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>360</td><td>13,036,072</td><td>62,685</td><td>481,481</td><td>3117/27</td><td>218</td><td>96</td><td>22</td></tr><tr><td>Intensive-NPIs</td><td>All</td><td>540</td><td>17,927,080</td><td>70,458</td><td>926,777</td><td>3877/145</td><td>291</td><td>249</td><td>181</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>180</td><td>7,771,015</td><td>54,196</td><td>107,426</td><td>2530/2</td><td>44</td><td>24</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>270</td><td>7,772,026</td><td>54,103</td><td>111,120</td><td>2502/8</td><td>45</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>360</td><td>7,782,238</td><td>54,307</td><td>117,712</td><td>2541/9</td><td>46</td><td>22</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>10+</td><td>540</td><td>7,775,685</td><td>54,302</td><td>122,420</td><td>2486/8</td><td>48</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>180</td><td>7,705,102</td><td>53,347</td><td>106,052</td><td>2571/7</td><td>43</td><td>23</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>270</td><td>7,778,149</td><td>54,255</td><td>110,905</td><td>2553/7</td><td>46</td><td>22</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>360</td><td>7,784,818</td><td>54,110</td><td>120,321</td><td>2593/6</td><td>48</td><td>24</td><td>0</td></tr><tr><td>Extended-NPIs</td><td>All</td><td>540</td><td>7,836,271</td><td>55,120</td><td>128,861</td><td>2655/5</td><td>48</td><td>26</td><td>0</td></tr></tbody></table>
## 2.8. Impact of Extended-NPIs on Vaccination Campaign Efficacy

With extended-NPIs, the impacts of the probability of viral transmission and duration of the vaccination campaign on mortality and ICU pressure were much smaller than with intensive-NPIs (Figure 3).

## 2.9. Reduction in Disease Transmission by NPIs and Relative Contagiousness of the Alpha Variant

Compared to moderate-NPIs, the estimated relative reductions in disease transmission were 0.783 for intensive-NPIs and 0.534 for extensive-NPIs ($\beta_{2,Intensive}$ and $\beta_{2,Extended}$, Table 6). With the historical strain as reference, the Alpha variant was estimated to be 1.572 times more infectious ($\beta_{3,Alpha}$, Table 6).

The estimated proportion of removed individuals was 11% on 1 January 2021, and the estimated prevalence of the Alpha variant represented 13.2% of SARS-CoV-2 strains on 27 January 2021 and 38.4% on 16 February 2021.

# 3. Discussion

The modelling proposed in this study provides objective and original estimates of the effects of NPIs at the population level. With estimated reductions in the probability of COVID-19 transmission of 25% ($\beta_{2,Intensive}$ = 0.783) and 50% ($\beta_{2,Extended}$ = 0.534), intensive-NPIs and extended-NPIs have slowed the dynamics of the COVID-19 epidemic. The first part of this work simulated the spread of the epidemic in the absence of vaccination. With the historical strain alone, adopting relaxed-NPIs would have led to the extinction of the epidemic in about 400 days but at the cost of a very high mortality (about 372 thousand deaths) and a serious ICU overcrowding. Adopting extended-NPIs would have led to much less deaths (about 61 thousand), without ICU overloading. Intensive-NPIs would have led to intermediate results.

The estimated relative infectiousness of the Alpha variant ($\beta_{3,Alpha}$ = 1.572) was in the range of the literature estimates [8]. The estimated prevalence values of the Alpha variant on 27 January 2021 (13.2%) and 16 February 2021 (38.4%) were close to the observed values on those dates (i.e., 14% and 36%, respectively). These results support the proposed multiplicative approach (see Methods, Equation (1)).

Strikingly, in the absence of vaccination, as soon as the variant strain was included in the simulations, extended-NPIs had to be imposed to avoid large increases in the number of deaths and important ICU overload as intensive-NPIs became ineffective. Extended-NPIs should be maintained, as long as herd immunity is achieved. The population would not be protected against a new epidemic wave if these protective measures were removed. These conclusions are in agreement with the recommendations of Bosetti et al. [9] that reinforced NPIs should be applied to more contagious variants to achieve an epidemic control, similar to that obtained by applying less stringent measures to the historical virus strain.

Within the current context, vaccination seems to be the sole solution to stop the epidemic; thus, it is important to vaccinate the whole population as soon as possible. Reaching herd immunity through vaccination is mandatory to avoid the resumption of the epidemic after NPIs removal. With intensive-NPIs and a vaccine reduction of virus transmission of 90%, extending the vaccination campaign from 6 to 12 months, would have a relatively moderate impact on the disease-specific mortality (increase by 6%) but a high impact on the need for ICU bed-days (+77%). However, maintaining intensive-NPIs would have deleterious health consequences, such as delayed surgical management of benign diseases, inadequate follow-up of chronic diseases, and disruptions of health prevention programs. Besides, intensive-NPIs have deleterious impacts on education, research, economy, etc. Thus, extending the vaccination campaign to 12 months would lead to significant pressure on ICUs and the entire health system, as well as to social impacts from maintaining intensive-NPIs during that period. A substantial increase in disease-related mortality of 17% would be seen if the vaccination campaign were extended over 18 months, with a strong need for ICU bed-days (+244%). Moreover, maintaining intensive-NPIs for 18 months seems unrealistic for social and economic reasons. All those elements
are in favour of the fastest possible mass vaccination, providing vaccine availability and population compliance.

Despite a significant acceleration of the current French vaccination campaign, the objective of vaccinating the entire French population (67 millions) in 6, or even 9, months has definitely not been achieved. In France, the vaccination campaign started on 11 January 2021. By 1 April 2021, 3.1 million people were fully vaccinated and 9.0 millions had received a first dose and, by 1 June 2021, 12.5 million people were fully vaccinated (16% of the population) and 26.4 millions had received a first dose. On 1 June 2021, the cumulative number of observed hospital deaths was 84 thousands, whereas 65 thousand deaths were predicted by the worst-case scenario (Table 5). Two main reasons explain this difference. The first is the vaccination coverage rate among the elderly. In simulated vaccination campaigns, all people aged 75 years or older received a first dose of vaccine on day 46 of 2021 (15 February) for the 18-month vaccination campaigns, and earlier for shorter campaigns. In contrast, by 15 February 2021, only 22.2% of individuals aged 75 years or older had actually received a first dose of vaccine. The second reason concerns NPIs. A relaxation of NPIs was seen in France during the first half of 2021: the curfew starting at 6 PM, decided in January, was postponed to 7 PM in March, 9 PM in May, and 11 PM in June; the isolation of travelers entering the territory was eased during a part of that period; and some business activities reopened in February, March, and May. An important additional argument in favour of limiting the duration of the vaccination campaign is the difficulty of maintaining acceptance of appropriate NPIs over a long period.

Another argument in favour of rapid and extensive routine vaccination is the race against the emergence of new VOC that is directly related to the number of infected individuals. Within this context, our opinion is that vaccinating the entire population may be necessary after ensuring the safety and efficacy of vaccination of children (see NCT04796896 and NCT04816643 trials on clinicaltrials.gov).

The present study makes the strong hypothesis of a 100% acceptance level of vaccination. Virus extinction will happen if the population is widely vaccinated. This means that more efforts should be made, in terms of communication and vaccination strategy. It would be interesting to tackle the issue of vaccine acceptance by measuring the impacts of partial acceptance levels on the course of the epidemic. Another modelling choice was ignoring the duration of vaccine protection. This might be the object of future investigations.

The model presented in this paper has been calibrated on the historical strain, variants with contagiousness similar to that of the Alpha variant, virulence similar to that of the historical strain, and the absence of vaccine resistance. Once sufficient data will be available to estimate a realistic $\beta_3$ for each variant, this model will allow comparing the effects of various vaccination strategies on mixtures of virus strains. It will also be possible to incorporate a variant-specific vaccine resistance parameter into the model, as soon as meaningful estimates of that parameter will be available. An increase in the virus virulence may be handled by changing the model transition probabilities (see [10]).

The results presented and discussed here do not take the Delta variant into account because the model was designed and developed before the emergence of that variant. The study was, thus, focused on the effect of the Alpha variant and the model calibration used Alpha-variant-specific data. Therefore, the results obviously under-estimate the real impact of the whole range of existing variants. The model used is sufficiently flexible to be adapted to measure the impacts of the Delta variant on the vaccination strategies as soon as sufficient data will be available. This is a perspective for a future work, that would also integrate scenarios about future more dangerous and/or contagious variants.

One advantage of the modelling presented in this paper is that it takes into account the particularity of any vaccination strategy. If it appears that it will not be possible to vaccinate all the population within one year, alternative solutions will have to be examined. This work did not consider the logistics details of the vaccination strategies; nevertheless, the model flexibility allows including those details. In addition, future works may tackle
the delay between the two vaccine doses or the pace of vaccination through increasing the number of vaccination centers, vaccinating personnel, or available vaccine supplies.

The present study considered the side effects of anti-COVID-19 vaccines to be negligi-
ble in comparison with the expected benefit at the population level (control of the epidemic
through herd immunity). In fact, the rare serious side effects can be taken into account
by considering the benefit-risk balance at an individual level, especially when the risk of
severe COVID-19 is very low (e.g., in young children). In a safety analysis of BNT162b2,
the mRNA COVID-19 vaccine in individuals aged 16 or older, the estimated excess of
myocarditis associated with the vaccine was 2.7 events per 100,000 vaccinations (95% CI:
1.0 to 4.6), whereas the estimated excess of myocarditis associated with COVID-19 disease
was 11.0 events per 100,000 cases (95% CI: 5.6 to 15.8) [11]. In a cohort of individuals aged
18–65 years, vaccinated with the ChAdOx1-S, the estimated excess of venous thromboem-
bolic events was 11 events per 100,000 vaccinations (95% CI: 5.6 to 17.0), of which an excess
of cerebral venous thrombosis of 2.5 events per 100,000 vaccinations (95% CI: 0.9 to 5.2) [12].
The reporting of those rare vaccine side effects influenced the French vaccination campaign
by increasing the use of the BNT162b2 mRNA COVID-19 vaccine. Underestimating the
vaccination benefit–risk balance might have a negative impact on the acceptance of a
vaccine and even vaccination.

Other limitations concern taking into account the details of the transmission and the severity of the disease. In particular, the model did not consider the living conditions (e.g., residence in cities vs. villages, holidays, etc.) that can influence the number of contacts and the application of barrier measures; in addition, it did not consider the prevalence of comorbidities [13]. Fortunately, the agent-based model used in this study simulated encounters and transmissions between individuals; thus, it may be extended to include any selection of individual characteristics. Furthermore, as the simulation can be run with a large number of individuals (more than a million), it may represent various comorbidities, living, habits, or other conditions within a given population. Such improvements of the core model may be used to investigate the potential effects of targeted public health decisions.

4. Methods

Section 4.1 presents the agent-based model and Section 4.2 shows the simulated vacci-
nation scenarios obtained from that model. The model calibration is presented Section 4.3.

4.1. *The Agent-Based Model*

The agent-based model (ABM) allowed simulating individual trajectories between specific states of disease occurrence and evolution. This model (Figure 4) extended the set of compartments initially proposed by Di Domenico et al. [14] to an ordinary differential equation (ODE) model. The ABM included two parts. The left-hand side of the model, which corresponds to ‘disease spreading’, reproduced the compartment pattern of Di Domenico et al. [14]. The right-hand part (the dotted frame in Figure 4) describes the ‘outcome of hospitalized cases’ and includes specific states to account explicitly for hospitalization, before and after ICU periods. Each state of the model was divided into nine parts for the nine age groups, as described by Gauchon et al. [15] (only three age groups were considered in the model of Di Domenico et al. [14]). The flexibility of the ABM allowed incorporating vaccination, variant occurrence, and types of NPI measures (three aspects that were not considered by Di Domenico et al. [14]). In the ABM model, the transition probabilities—that are similar to those of a multi-state Markov model—were applied at individual levels. The states and trajectories are summarized in Figure 4; the parameters estimated or extracted from the literature are shown in Tables 6 and 7.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** The states and their connections in the model. The left-hand part of the model represents the disease spreading with the average sojourn time beneath each compartment. The right-hand part of the model (bottom-right dotted frame) represents the hospitalized cases and their outcomes. It is composed of states $H_1$, $H_2$, $H_3$, ICU, $R_2$, and D. Each arrow is labeled with the individual daily transition probability between two states.</figcaption></figure>

<table><caption>**Table 6.** Estimated parameters (the values in italics for the left-hand part of the model are extracted from a previous work [15]).</caption><thead><tr><th>Age</th><th>0-9</th><th>10-19</th><th>20-29</th><th>30-39</th><th>40-49</th><th>50-59</th><th>60-69</th><th>70-79</th><th>80+</th></tr></thead><tbody><tr><td colspan="10" style="text-align:center;">Left-hand part of the model</td></tr><tr><td>β<sub>1,age</sub></td><td>0.01</td><td>0.01</td><td>0.017</td><td>0.008</td><td>0.007</td><td>0.012</td><td>0.016</td><td>0.16</td><td>0.164</td></tr><tr><td>p<sub>1ps</sub></td><td>0.249</td><td>0.249</td><td>0.247</td><td>0.243</td><td>0.242</td><td>0.225</td><td>0.209</td><td>0.189</td><td>0.031</td></tr><tr><td>p<sub>1ms</sub></td><td>0.746</td><td>0.748</td><td>0.741</td><td>0.730</td><td>0.727</td><td>0.674</td><td>0.626</td><td>0.568</td><td>0.092</td></tr><tr><td>p<sub>1ss</sub></td><td>0.006</td><td>0.003</td><td>0.012</td><td>0.026</td><td>0.031</td><td>0.102</td><td>0.166</td><td>0.243</td><td>0.877</td></tr><tr><td>β<sub>2,Intensive</sub></td><td colspan="9">0.783 for all ages</td></tr><tr><td>β<sub>2,Extended</sub></td><td colspan="9">0.534 for all ages</td></tr><tr><td>β<sub>3,Alpha</sub></td><td colspan="9">1.572 for all ages</td></tr><tr><td colspan="10" style="text-align:center;">Right-hand part of the model for<br>historical strain (calibration period September 2020–December 2020)</td></tr><tr><td>p<sub>ICU|H1</sub></td><td>0.151</td><td>0.5</td><td>0.077</td><td>0.078</td><td>0.117</td><td>0.146</td><td>0.228</td><td>0.227</td><td>0.128</td></tr><tr><td>p<sub>H2|H1</sub></td><td>0.849</td><td>0.95</td><td>0.923</td><td>0.922</td><td>0.874</td><td>0.854</td><td>0.771</td><td>0.773</td><td>0.872</td></tr><tr><td>p<sub>H3|ICU</sub></td><td>0.149</td><td>0.055</td><td>0.143</td><td>0.065</td><td>0.058</td><td>0.049</td><td>0.051</td><td>0.062</td><td>0.22</td></tr></tbody></table>
<table><caption>**Table 6.** *Cont.*</caption><thead><tr><th>Age</th><th>0-9</th><th>10-19</th><th>20-29</th><th>30-39</th><th>40-49</th><th>50-59</th><th>60-69</th><th>70-79</th><th>80+</th></tr></thead><tbody><tr><td>p<sub>D|H2</sub></td><td>0</td><td>0</td><td>0</td><td>0.001</td><td>0.001</td><td>0.002</td><td>0.006</td><td>0.012</td><td>0.019</td></tr><tr><td>p<sub>R|H2</sub></td><td>0.526</td><td>0.462</td><td>0.442</td><td>0.265</td><td>0.188</td><td>0.142</td><td>0.098</td><td>0.07</td><td>0.036</td></tr><tr><td>p<sub>D|H3</sub></td><td>3 × 10<sup>-4</sup></td><td>5 × 10<sup>-4</sup></td><td>0.001</td><td>0.002</td><td>0.005</td><td>0.007</td><td>0.009</td><td>0.009</td><td>0.008</td></tr><tr><td>p<sub>R|H3</sub></td><td>0.039</td><td>0.001</td><td>0.018</td><td>0.016</td><td>0.016</td><td>0.01</td><td>0.01</td><td>0.009</td><td>0.067</td></tr><tr><td colspan="10" style="text-align: center;">Right-hand part of the model for<br>Alpha variant (calibration period March 2021-April 2021)</td></tr><tr><td>p<sub>ICU|H1</sub></td><td>0.052</td><td>0.031</td><td>0.077</td><td>0.071</td><td>0.116</td><td>0.157</td><td>0.280</td><td>0.268</td><td>0.098</td></tr><tr><td>p<sub>H2|H1</sub></td><td>0.948</td><td>0.967</td><td>0.902</td><td>0.918</td><td>0.873</td><td>0.841</td><td>0.396</td><td>0.732</td><td>0.895</td></tr><tr><td>p<sub>H3|ICU</sub></td><td>0.079</td><td>0.033</td><td>0.070</td><td>0.033</td><td>0.033</td><td>0.033</td><td>0.075</td><td>0.057</td><td>0.139</td></tr><tr><td>p<sub>D|H2</sub></td><td>0</td><td>0</td><td>0.001</td><td>0.001</td><td>0.002</td><td>0.003</td><td>0.009</td><td>0.012</td><td>0.016</td></tr><tr><td>p<sub>R|H2</sub></td><td>0.288</td><td>0.193</td><td>0.165</td><td>0.151</td><td>0.115</td><td>0.082</td><td>0.040</td><td>0.035</td><td>0.030</td></tr><tr><td>p<sub>D|H3</sub></td><td>0</td><td>0</td><td>0.001</td><td>0.005</td><td>0.004</td><td>0.007</td><td>0.001</td><td>0.004</td><td>0.014</td></tr><tr><td>p<sub>R|H3</sub></td><td>0.084</td><td>0.052</td><td>0.211</td><td>0.057</td><td>0.055</td><td>0.053</td><td>0.102</td><td>0.059</td><td>0.073</td></tr></tbody></table>

<table><caption>**Table 7.** Parameters sojourn time, relative infectiousness, and rate of asymptomatic subjects, as obtained from the literature.</caption><thead><tr><th>Parameter</th><th>Sources</th><th>Value</th></tr></thead><tbody><tr><td colspan="3">Sojourn time</td></tr><tr><td>t<sub>i</sub></td><td>see [15]</td><td>5.1</td></tr><tr><td>t<sub>p</sub></td><td>see [15]</td><td>1.5</td></tr><tr><td>t<sub>s</sub></td><td>see [15]</td><td>7</td></tr><tr><td>t<sub>bh</sub></td><td>see [15]</td><td>3</td></tr><tr><td colspan="3">Relative infectiousness</td></tr><tr><td>i<sub>Im</sub></td><td>[16]</td><td>1</td></tr><tr><td>i<sub>Is</sub></td><td>[16]</td><td>1</td></tr><tr><td>i<sub>A</sub></td><td>[16]</td><td>0.55</td></tr><tr><td>i<sub>Ip</sub></td><td>[16]</td><td>0.55</td></tr><tr><td>i<sub>Ip</sub></td><td>[16]</td><td>0.55</td></tr><tr><td colspan="3">Rate of asymptomatic subjects</td></tr><tr><td>p<sub>A</sub></td><td>[14,17]</td><td>0.20</td></tr><tr><td colspan="3">Proportion new variant (8 January 2021)</td></tr><tr><td>Alpha variant</td><td>Santé Publique France</td><td>3.3%</td></tr></tbody></table>

### 4.1.1. The Model States

The model is a stochastic ABM, with a discrete time step of one day. The states used in the model are the following:

* Individuals who were not infected by the virus begin in state S (susceptible);

* The incubation period includes two states:

- when infected (see next section), an individual moves from state S to state E (exposed), which contains infected individuals who did not develop symptoms yet and are not contagious. The mean sojourn time in E was $(t_i - t_p)$, $t_i$ being the incubation period and $t_p$ the duration of the prodromal state (see next item);

- when contagious, an individual is transferred from state E to state $I_p$ (prodromal phase, which is the short phase that follows contamination without symptoms but
with possible non-specific prodromes). After an average stay $t_p$, the individuals moved to one of the four states A, $I_{ps}$, $I_{ms}$, and $I_{ss}$, defined hereafter.

* A (asymptomatic state, with probability $p_a$): individuals who completed the incubation period, became infectious, but without disease symptoms. The mean sojourn time in A was $t_s$;

* The symptomatic infectious period includes three states for individuals developing symptoms, (with probability $1-p_a$):

- $I_{ps}$ (paucisymptomatic disease): individuals with weak disease symptoms;

- $I_{ms}$ (medium symptoms): individuals with disease symptoms (e.g., fever or cough) who did not require hospitalization. The average sojourn time in states $I_{ps}$ and $I_{ms}$ are the same as in state A;

- $I_{ss}$ (severe symptoms): individuals severely infected requiring hospitalization. They stayed in $I_{ss}$ before being hospitalized. The mean sojourn time in $I_{ss}$ was $t_{bh}$.

Given the onset of symptoms, the probabilities of states $I_{ps}$, $I_{ms}$, and $I_{ss}$ are, respectively, $p_{I_{ps}}$, $p_{I_{ms}}$, and $p_H$ (probabilities summing to one).

An individual leaving states A, $I_{ps}$ or $I_{ms}$ ends in state $R_1$ (removed individuals).

* When leaving $I_{ss}$, individuals enter the hospitalization period, which corresponds to four states:

- $H_1$: individuals hospitalized before it was determined whether they needed intensive care or not;

- $H_2$: individuals hospitalized without the need for ICU;

- ICU: individuals hospitalized in ICU;

- $H_3$: individuals hospitalized after leaving ICU.

From hospital, all individuals end in one of the two absorbing states:

- D (deceased at hospital);

- $R_2$ (removed): individuals who recovered.

During the hospitalization period, each individual follows a Markov chain dynamics (the dotted frame in Figure 4), with daily transition probabilities noted $p_{H2|H1}$, $p_{ICU|H1}$, $p_{H3|ICU}$, $p_{D|H2}$, $p_{R|H2}$, $p_{D|H3}$, and $p_{R|H3}$.

The age distribution of the population came from the INSEE (Institut National de la Statistique et des Études Économiques); nine age groups were considered: 8 ten-year age groups from 0 to 79 years, plus an extra age-group with individuals aged 80 or older. The contact matrix was available for France, on the basis of these age groups [18]. The number of daily contacts per individual was set at the beginning of the simulation and drawn from a Poisson distribution, with a parameter set to the corresponding element of the contact matrix. At each time step, the daily contacts between infectious and susceptible agents were drawn at random.

### 4.1.2. Transition Probabilities

The disease spreads with infectious contacts that move individuals from state S to state E. Daily contacts between agents were based on a contact matrix (C), the element $(C_{i,j})$ of C being the mean number of individuals from age group j encountered per day by individuals from age group i. $C_i = \sum_j C_{i,j}$ is the mean number of individuals encountered per day by individuals from age group i. The daily number of individuals encountered by a new agent was obtained from a Poisson distribution, with parameter $C_i$, the expected proportions of contacts aged j being $C_{i,j}/C_i$.

The probability ($P_{\text{infect}}$) of an individual susceptible to infection because of a contact with an infectious individual was decomposed into a product of four terms:

$$P_{\text{infect}}(\text{age}, \text{NPI}, \text{strain}, Z) = \beta_{1,\text{age}} \beta_{2,\text{NPI}} \beta_{3,\text{strain}} i_Z \quad (1)$$

where:
* $\beta_{1,age}$ was the estimated baseline probability of infection, as a function of age (based on the age of the susceptible individual).

* $\beta_{2,NPI}$, NPI $\in$ {Moderate, Intensive, Extended}, estimated the effects of three NPI levels: moderate, intensive, and extended. The underlying hypotheses were that the NPI effects are independent from the state of the infectious individual, age group, or virus strain. Fixing $\beta_{2,Moderate} = 1$, parameters $\beta_{2,Intensive}$, and $\beta_{2,Extended}$ estimated the relative reductions in disease transmission under intensive- and extensive-NPIs, in comparison to moderate-NPIs. Parameter $\beta_{2,NPI}$ was interpreted as a coefficient, reducing the mean number of daily contacts, mean duration of contacts, and/or their infectiousness ($0 \le \beta_{2,NPI} \le 1$).

* $\beta_{3,strain}$, strain $\in$ {historical, Alpha}, estimated the relative contagiousness of each virus strains, $\beta_{3,Alpha}$, estimating the relative contagiousness of the Alpha variant, compared to that of the historical strain ($\beta_{3,historical} = 1$).

* $Z \in$ {E, $I_p$, A, $I_{ps}$, $I_{ms}$, $I_{ss}$} is the state of the infectious individual, $i_Z$ being the relative infectiousness of individuals in state Z. It was fixed to 1 for $I_{ms}$ and $I_{ss}$ states (reference states), and to 0.55 for $I_p$, $I_A$, and $I_{ps}$ states, according to the literature (Table 7).

Sojourn times, in states (E, A, I_s, I_p, I_m, and I_s), were generated from Weibull distributions. When more than one destination state was possible, the transition was selected at random using the corresponding probabilities (i.e., $p_A$, $p_{Ips}$, $p_{Ims}$, and $p_{Iss}$).

For the right-hand part of the model (dotted frame, Figure 4), each day, for each individual in state $H_1$, $H_2$, $H_3$, or ICU, the state transition was drawn with probabilities $p_{H2|H1}$, $p_{ICU|H1}$, $p_{H3|ICU}$, $p_{D|H2}$, $p_{R|H2}$, $p_{D|H3}$, and $p_{R|H3}$, where staying in the same state (i.e., self-loops) were the complementary events.

The proportion of asymptomatic forms was set to 20% (Table 7). The relative infectiousness of prodromal, asymptomatic, and paucisymptomatic forms of the historical strain was set to $i_{Ip} = i_A = i_{Ips} = 0.55$; medium and severe symptoms forms were taken as references, with $i_{Ims} = i_{Iss} = 1$ (Table 7). The proportions of various symptomatic forms, paucisymptomatic forms, medium symptoms, and severe symptoms were specified, depending on age (Table 6).

## 4.2. Scenarios

In agreement with the three levels of NPIs, three scenarios were defined, the first one being the use of intensive-NPIs and the second being the use of extended-NPIs. The third scenario, termed relaxed-NPIs, was built as an alternation of 45-day periods of intensive and moderate protective measures, the first period being one of intensive measures.

The expected effects of the three protective scenarios were then compared.

### 4.2.1. Vaccination Effects

The expected effects of vaccination were set at their most recently estimated levels of efficacy in preventing severe forms of the disease; that is, 94% [5,6] and 75% [4], depending on the vaccine. When vaccination was effective, a vaccinated individual could develop only an asymptomatic form of the disease. Putative reductions of virus transmission were set to 50, 75, and 90%, the reduction of transmission probability being independent of the vaccine. An age-based priority for being vaccinated was used. Four vaccination priorities were considered and all individuals of a given age group had to receive their first injection before giving the first injections to the next priority age group. The age priorities considered were: 75 and older, 65 and older, 55 and older, and then the rest of the population. In the rest of the population, two vaccination options were compared: vaccination of individuals aged 10 years and over vs. vaccination at all ages (vaccination is currently under investigation in the youngest population). The effects of four vaccination campaigns, over 6, 9, 12, and 18 months, were compared. At the end of each vaccination campaign, all individuals eligible to vaccination should have received their second injection. In all simulations, vaccines with 94% efficacy were attributed to all individuals aged 75 or older and to all individuals aged 54 or younger. For individuals aged between 55 and 74, vaccines with
94% efficacy were attributed to 5/6 of individuals, whereas vaccines with 75% efficacy were attributed to the remaining 1/6 of individuals, aged between 55 and 74. These choices were made to fit vaccine availability for different age groups, as applied in France. An individual was considered to be fully vaccinated after receiving two injections of any vaccine with an interval of 42 days. In all simulations, vaccine supplies were considered sufficient to follow the study vaccination schedule. Furthermore, all vaccines were considered as efficient against the historical strain as against new variants, similar to the Alpha variant [19].

### 4.2.2. Statistical Analyses

Due to the stochasticity of the agent-based simulation, 50 simulations were performed with each set of parameters. Each run involved 1,300,000 agents (i.e., 2% of the French metropolitan population—nearly 64.5 million inhabitants). To obtain estimations for France, the simulation results were multiplied by 50. The mean values and standard errors of the mean across the runs were plotted.

### 4.2.3. Outcome Criteria

The outcome criteria were: the cumulative number of hospital deaths of COVID-19 patients, as well as the cumulative number of removed individuals due to the disease, from the beginning of the epidemic up to 30 June 2022, and the daily prevalence of COVID-19 patients in ICU beds and its saturation indicators during the simulation period (1 January 2021 to 30 June 2022).

The saturation of ICU beds was calculated as the cumulative number of cases requiring the ICU when all beds were already occupied. Three thresholds were considered for scenarios without vaccination: 5000 (the current number of ICU beds), 8000 (adding ICU beds in intermediate care units), and 12,000 (adding ICU beds in conventional care units).

To evaluate the occupation rate of the ICU under the various vaccination campaign scenarios, the numbers of days during which ICU occupancy exceeded the thresholds of 1000, 2000, and 3000 beds were calculated.

### 4.3. Model Calibration

The parameters of the model (Figure 4) obtained from the literature are shown in Table 7, and those estimated by calibration are shown in Table 6. The calibration was performed using public data from Santé Publique France (https://www.data.gouv.fr/fr/datasets/donnees-hospitalieres-relatives-a-lepidemie-de-covid-19/ accessed on 15 October 2020). Comparisons between those reference data and the simulation outputs, obtained from the estimated parameters, are shown in supplementary Figures S1–S5.

The left-hand part of the model was calibrated independently of the right-hand part (the dotted frame in Figure 4). For the left-hand part ('disease spreading'), the parameters were estimated by likelihood maximization using daily hospital admissions from 1 September 2020 to 31 January 2021. The maximization was performed by the Nelder-Mead simplex method [20]. To reduce stochastic error, the calibration used daily hospital admissions averaged over 50 simulations. During each simulation, depending on the period, parameter $\beta_{2,NPI}$ took three values, with $NPI \in \{\text{Moderate, Intensive, Extended}\}$. It was set to $\beta_{2,Moderate}$ from 1 September 2020 to 29 October 2020 (before the second lockdown in France), then to $\beta_{2,Extended}$ from 30 October 2020 to 30 November 2020 (second lockdown) and $\beta_{2,Intensive}$ from 1 December 2020 to 31 January 2021.

Parameter $\beta_{2,Moderate}$ was set to 1, as a reference, whereas $\beta_{2,Intensive}$ and $\beta_{2,Extended}$ were calibrated. The last parameter $\beta_{3,Alpha}$ was estimated as the relative infectiousness of the Alpha variant, in comparison with the historical strain, for which $\beta_{3,historical}$ was set to 1. Santé Publique France has estimated that, in France, on 8 January 2021, the Alpha variant was responsible for 3.3% of all new infections. In the simulations, during calibration, the presence of the Alpha variant was initiated on this date to 3.3%.

The right-hand part of the model ('outcomes of hospitalized cases') was calibrated by likelihood maximization, using the daily prevalence of D (deceased at hospital), ICU,
$R_2$ (recovered after hospitalization), and daily prevalence of hospitalization (sum of $H_1$, $H_2$, and $H_3$). During calibration the input for this part of the model (i.e., the new cases arriving in state $H_1$) was the daily hospital admissions obtained from Santé Publique France, and, thus, did not require using the right-hand part of the model in these simulations. To handle parameter bounds (e.g., make all daily probabilities range between 0 and 1), the likelihood maximization relied on the L-BFGS-B method [21]. Due to the sensitivity to the daily hospital admissions used as input, all Santé Publique France data were smoothed using a sliding average over an interval of three days.

This right-hand part of the model could not be calibrated over a single period ranging from 1 September 2020 to 30 April 2021 (no satisfying convergence). This is likely to be due to changes in protocols and/or some probabilities (e.g., ICU stay or death), relative to the Alpha variant. Thus, it has been calibrated over two different periods. The first one ranged from 1 September to 31 December 2020, to grasp the effect of the historical strain (very small number of cases, related to new variants in France over that period). The second period was from 1 March to 30 April 2021, to capture the dynamics of the new variants (majority of cases of Alpha variant in France over that period). Then, in the different scenarios with the presence of a variant, the values of the parameters obtained for the first period were used in the simulation steps performed before 1 March 2021, while the second set of parameter values was used for the simulation steps from 1 March 2021. In scenarios where there was no variant, only the first set of parameter values was used in all simulation periods.

**5. Conclusions**

The main conclusion of this study is that the key point in the race against the historical COVID-19 strain and its variants is to optimise a combination of vaccination campaign and adequate NPIs. The proposed adaptive model allowed estimating the combined effect of vaccination campaign duration, virus contagiousness, and NPI level on the historical strain and Alpha variant. It appeared mandatory to vaccinate the majority of the French population within one year or less.

The flexibility of the model allows it to follow up the course of the current epidemic, taking into account the progress of the vaccination coverage and emergence of new variants, with the aim of adapting the NPIs until achieving herd immunity.

**Supplementary Materials:** The following are available at https://www.mdpi.com/2076-393X/9/12/1462/s1. Figure S1: The calibrated hospital admission from 1 September 2020 to 31 January 2021. Figure S2: The calibrated prevalence of *H* and *ICU* from 1 September 2020 to 31 December 2020. Figure S3: The calibrated prevalence of $R_2$ and *D* from 1 September 2020 to 31 December 2020. Figure S4: The calibrated prevalence of *H* and *ICU* from 1 March 2021 to 30 April 2021. Figure S5: The calibrated prevalence of $R_2$ and *D* from 1 March 2021 to 30 April 2021.

**Author Contributions:** S.P., C.P., C.R., A.E.-L., J.-P.B., A.B., N.L., N.P., R.G., S.L. and P.R. designed the model; S.P., C.P., C.R., N.P. and R.G. implemented and ran the model; S.P., C.P., C.R., A.E.-L., J.-P.B., A.B., N.P., R.G., F.G., N.L., S.L. and P.R. performed the analyses; S.P., C.P., C.R., A.E.-L., J.-P.B., A.B., N.L. and P.R. drafted and revised the text; F.G. and S.L. edited and commented on the text; P.V. and J.I. revised critically the manuscript for scientific content; P.R. and S.L. coordinated the project. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by Région Auvergne-Rhône-Alpes, Initiative d'excellence (IDEXLYON)—Université de Lyon, Université Claude Bernard Lyon 1, École Centrale de Lyon, and by the PHRCi-2020 COVID-3S, supported by GIRCI Auvergne Rhône-Alpes. This work was also supported by the excellence chair 'DIALog, Digital Insurance And Long-term risks', financed by CNP Assurance and the research initiative 'Actuariat Durable et stabilité du secteur de l'assurance à long terme' financed by Milliman, joined projects with ISFA SAF Lab and Fondation du Risque.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.
**Data Availability Statement:** All data generated during this study are available from the corresponding author on reasonable request. The prevalence values of Alpha variant in France are publicly available from Santé Publique France (https://www.santepubliquefrance.fr/ (accessed on 15 October 2020)). The data used for estimation are also publicly available from Santé Publique France (https://www.data.gouv.fr/fr/datasets/donnees-hospitalieres-relatives-a-lepidemie-de-covid-19/ (accessed on 15 October 2020)). Prevalence of Removed individuals in France in 1 January 2021, are available from the publicly available data of Institut Pasteur (https://modelisation-covid19.pasteur.fr/realtime-analysis/infected-population/ (accessed on 15 October 2020)). The source code is available on gitlab: https://gitlab.com/Alnindra/covdyn_abm/ (accessed on 15 October 2020).

**Conflicts of Interest:** The authors declare no conflicts of interest.

## Appendix A. Members of the CovDyn (Covid Dynamics) Group, by Main Affiliation

### LBBE-HCL

Mathieu Fauvernier, François Gueyffier, Jean Iwaz, Delphine Maucort-Boulch, Simon Pageaud, Muriel Rabilloud, Pascal Roy.

### LSAF

Alexis Bienvenüe, Anne Eyraud-Loisel, Romain Gauchon, Pierre-Olivier Goffard, Nicolas Leboisne, Stéphane Loisel, Xavier Milhaud, Denys Pommeret, Yahia Salhi.

### LAMA

Pierre Vandekerkhove.

### CIRI-HCL

Philippe Vanhems.

### LMFA

Jean-Pierre Bertoglio.

### LTDS

Nicolas Ponthus.

### INL

Christelle Yeromonahos.

### LIRIS

Stéphane Derrode, Véronique Deslandres, Mohand-Saïd Hacid, Salima Hassas, Catherine Pothier, Christophe Rigotti.

### ICJ

Thibault Espinasse, Samuel Bernard, Philippe Michel, Vitaly Volpert.

### INRIA

Mostafa Adimy.

### CHU de Rouen Unité INSERM 1018, CESP

Jacques Benichou.

### CHU de Rouen LIMICS INSERM U1142, Université de Rouen/Sorbonne Université

Stéfan Darmoni.

### CHU de NICE, Université de Nice

Pascal Staccini.
## References

1. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. *N. Engl. J. Med.* **2020**, *382*, 727–733. [CrossRef] [PubMed]

2. Volz, E.; Mishra, S.; Chand, M.; Barrett, J.C.; Johnson, R.; Geidelberg, L.; Hinsley, W.R.; Laydon, D.J.; Dabrera, G.; O'Toole, Á.; et al. Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data. *medRxiv* **2021**, [CrossRef]

3. European Centre for Disease Prevention and Control. In *Emergence of SARS-CoV-2 B.1.617 Variants in India and Situation in the EU/EEA*; European Centre for Disease Prevention and Control: Solna Municipality, Sweden, 2021.

4. Voysey, M.; Clemens, S.A.C.; Madhi, S.A.; Weckx, L.Y.; Folegatti, P.M.; Aley, P.K.; Angus, B.; Baillie, V.L.; Barnabas, S.L.; Bhorat, Q.E.; et al. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: An interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. *Lancet* **2021**, *397*, 99–111. [CrossRef]

5. Baden, L.R.; El Sahly, H.M.; Essink, B.; Kotloff, K.; Frey, S.; Novak, R.; Diemert, D.; Spector, S.A.; Rouphael, N.; Creech, C.B.; et al. Efficacy and safety of the mRNA-1273 SARS-CoV-2 vaccine. *N. Engl. J. Med.* **2021**, *384*, 403–416. [CrossRef] [PubMed]

6. Polack, F.P.; Thomas, S.J.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Perez, J.L.; Pérez Marc, G.; Moreira, E.D.; Zerbini, C.; et al. Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine. *N. Engl. J. Med.* **2020**, *383*, 2603–2615. [CrossRef] [PubMed]

7. Pageaud, S.; Ponthus, N.; Gauchon, R.; Pothier, C.; Rigotti, C.; Eyraud-Loisel, A.; Bertoglio, J.P.; Bienvenüe, A.; Gueyffier, F.; Vanhems, P.; et al. Adapting French COVID-19 vaccination campaign duration to variant dissemination. *medRxiv* **2021**, [CrossRef]

8. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D.; et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. *Science* **2021**, *372*, 141–150. [CrossRef] [PubMed]

9. Bosetti, P.; Kiem, C.T.; Andronico, A.; Paireau, J.; Bruhl, D.L.; Lina, B.; Yazdanpanah, Y.; Fontanet, A.; Colizza, V.; Boëlle, P.Y.; et al. A Race between SARS-CoV-2 Variants and Vaccination: The Case of the B.1.1.7 Variant in France. HAL-pasteur-03149525. 2021. Available online: https://hal-pasteur.archives-ouvertes.fr/pasteur-03149525 (accessed on 2 March 2021).

10. Challen, R.; Brooks-Pollock, E.; Read, J.M.; Dyson, L.; Tsaneva-Atanasova, K.; Danon, L. Increased hazard of mortality in cases compatible with SARS-CoV-2 variant of concern 202012/1-a matched cohort study. *medRxiv* **2021**, [CrossRef]

11. Barda, N.; Dagan, N.; Ben-Shlomo, Y.; Kepten, E.; Waxman, J.; Ohana, R.; Hernán, M.A.; Lipsitch, M.; Kohane, I.; Netzer, D.; et al. Safety of the BNT162b2 mRNA Covid-19 vaccine in a nationwide setting. *N. Engl. J. Med.* **2021**, *385*, 1078–1090. [CrossRef] [PubMed]

12. Pottegård, A.; Lund, L.C.; Karlstad, Ø.; Dahl, J.; Andersen, M.; Hallas, J.; Lidegaard, Ø.; Tapia, G.; Gulseth, H.L.; Ruiz, P.L.D.; et al. Arterial events, venous thromboembolism, thrombocytopenia, and bleeding after vaccination with Oxford-AstraZeneca ChAdOx1-S in Denmark and Norway: Population based cohort study. *BMJ* **2021**, *373*, n1114. [CrossRef] [PubMed]

13. Kiem, C.T.; Massonnaud, C.R.; Levy-Bruhl, D.; Poletto, C.; Colizza, V.; Bosetti, P.; Fontanet, A.; Gabet, A.; Olié, V.; Zanetti, L.; et al. A modelling study investigating short and medium-term challenges for COVID-19 vaccination: From prioritisation to the relaxation of measures. *EClinicalMedicine* **2021**, *38*, 101001. [CrossRef] [PubMed]

14. Di Domenico, L.; Pullano, G.; Sabbatini, C.E.; Boëlle, P.Y.; Colizza, V. Impact of lockdown on COVID-19 epidemic in Île-de-France and possible exit strategies. *BMC Med.* **2020**, *18*, 240. [CrossRef] [PubMed]

15. Gauchon, R.; Ponthus, N.; Pothier, C.; Rigotti, C.; Volpert, V.; Derrode, S.; Bertoglio, J.P.; Bienvenue, A.; Goffard, P.O.; Eyraud-Loisel, A.; et al. Lessons learnt from the use of compartmental models over the COVID-19 induced lockdown in France. *medRxiv* **2021**, [CrossRef]

16. Li, R.; Pei, S.; Chen, B.; Song, Y.; Zhang, T.; Yang, W.; Shaman, J. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). *Science* **2020**, *368*, 489–493. [CrossRef] [PubMed]

17. Beale, S.; Hayward, A.; Shallcross, L.; Aldridge, R.W.; Fragaszy, E. A Rapid Review of the Asymptomatic Proportion of PCR-Confirmed SARS-CoV-2 Infections in Community Settings. 2020. Available online: https://www.gov.uk/government/publications/nervtag-rapid-review-of-the-asymptomatic-proportion-of-pcr-confirmed-sars-cov-2-infections-in-community-settings-9-september-2020 (accessed on 23 October 2020).

18. Prem, K.; Cook, A.R.; Jit, M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. *PLoS Comput. Biol.* **2017**, *13*, e1005697. [CrossRef] [PubMed]

19. Wang, P.; Nair, M.S.; Liu, L.; Iketani, S.; Luo, Y.; Guo, Y.; Wang, M.; Yu, J.; Zhang, B.; Kwong, P.D.; others. Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7. *Nature* **2021**, *593*, 130–135. [CrossRef] [PubMed]

20. Nelder, J.A.; Mead, R. A Simplex Method for Function Minimization. *Comput. J.* **1965**, *7*, 308–313. [CrossRef]

21. Byrd, R.H.; Lu, P.; Nocedal, J.; Zhu, C. A Limited Memory Algorithm for Bound Constrained Optimization. *SIAM J. Sci. Comput.* **1995**, *16*, 1190–1208. [CrossRef]