Skip to main content
Advertisement
  • Loading metrics

The importance of non-pharmaceutical interventions during the COVID-19 vaccine rollout

Abstract

The promise of efficacious vaccines against SARS-CoV-2 is fulfilled and vaccination campaigns have started worldwide. However, the fight against the pandemic is far from over. Here, we propose an age-structured compartmental model to study the interplay of disease transmission, vaccines rollout, and behavioural dynamics. We investigate, via in-silico simulations, individual and societal behavioural changes, possibly induced by the start of the vaccination campaigns, and manifested as a relaxation in the adoption of non-pharmaceutical interventions. We explore different vaccination rollout speeds, prioritization strategies, vaccine efficacy, as well as multiple behavioural responses. We apply our model to six countries worldwide (Egypt, Peru, Serbia, Ukraine, Canada, and Italy), selected to sample diverse socio-demographic and socio-economic contexts. To isolate the effects of age-structures and contacts patterns from the particular pandemic history of each location, we first study the model considering the same hypothetical initial epidemic scenario in all countries. We then calibrate the model using real epidemiological and mobility data for the different countries. Our findings suggest that early relaxation of safe behaviours can jeopardize the benefits brought by the vaccine in the short term: a fast vaccine distribution and policies aimed at keeping high compliance of individual safe behaviours are key to mitigate disease resurgence.

Author summary

The start of vaccination campaigns is a decisive turning point in the global effort against COVID-19. Nonetheless, at least in the short and medium-term, vaccine availability and the logistical issues associated with an unprecedented mass vaccination suggest that non-pharmaceutical interventions will still play an important role in virus containment. Here, we propose an epidemic model to study the possible effects induced by a relaxation of COVID-safe behaviours in response to the vaccine rollout. Individuals may see this milestone as the end of the emergency and thus give up preventive measures potentially exposing themselves to higher infection risk. We explore the interplay between such behavioural changes and different population pyramids, contact patterns, epidemic conditions, vaccine allocation strategies, rollout speed, and vaccine efficacy. We show that early relaxation of COVID-safe behaviours can jeopardize and even nullify the benefit brought by the vaccine in the short and medium-term. Our results indicate that a high level of compliance to NPIs during vaccines rollout is crucial to avoid hindering the gigantic effort of the vaccination campaigns.

Introduction

The COVID-19 pandemic has been largely fought with non-pharmaceutical interventions (NPIs). Bans of events and social gatherings, limitations in national and international travels, school closures, shifts towards remote working, curfews, closure of pubs and restaurants, cordon sanitaires, national and regional lockdowns are examples of governmental interventions implemented around the world to curb the spreading of SARS-CoV-2 [15]. While extremely effective, such top-down NPIs induce profound behavioural changes, bring many social activities to a halt, and thus have huge socio-economic costs. Hence, alongside these measures, governments nudged and/or mandated populations to adopt another set of NPIs. Social distancing, face masks, and increased hygiene are examples [6]. Although far from being cost-free, they are more feasible, sustainable, and allow for higher levels of socio-economic activity. As such, they have been the leitmotif of the post first COVID-19 wave in many countries. Unfortunately, awareness, adoption, and compliance with these NPIs have been spotty [6]. Furthermore, they have not been complemented with sufficiently aggressive test and trace programs. As result, many countries experienced marked disease resurgences after the summer 2020 and some had to resort to new lockdowns [7].

As we write, we have turned a crucial corner in hampering the resurgence, diffusion and severe outcomes of the disease. Several vaccines have shown great results of their phase 3 trials, and have been authorized for emergency use by several regulating agencies and tens of others are in the pipeline [810]. As result, we have witnessed the start of vaccination campaigns around the world. However, the logistical issues linked to the production, delivery, and administration of billions of doses on a global scale impose caution when evaluating the impact vaccines will have on the pandemic in the short and medium term. They will be a scarce resource and it will take time to vaccinate the fraction of the population necessary for herd immunity [11, 12]. Furthermore, vaccines are not perfect and there are many other unknowns [13]. The extent to which vaccines limit further transmission, how long the immunity lasts, and the levels of protection offered against new variants are key features currently under investigation [13]. So far, the results show high levels, around 95% [14, 15], of direct protection against the disease. Despite these first signs of the vaccine impact in the real world are extremely encouraging [16, 17], these figures might end up to be lower. Finally, vaccines’ acceptance is a complex challenge. A recent survey among 13, 426 participants in 19 countries shows that, while 71.5% of the sample is very or somewhat likely to take the vaccine, there are large heterogeneities [18]. Acceptance rates vary from 90% in China to less than 55% in Russia. Furthermore, they are linked to socio-economic, socio-demographic features, and education attainment [18]. Arguably, vaccines alone will not be able to contain the spreading of the virus, at least in the short term [19, 20]. Social distancing, face masks, hygiene measures, and other NPIs will be still key during the delivery of such vaccination programs, especially considering the emergence of variants of concern as well as the challenges in vaccines procurement and administration in the developing world [21].

In this context, an important question emerges. What will happen to adoption and compliance to NPIs as vaccination campaigns progress? Their arrival and delivery might induce individual and collective behavioural changes. Some might see this milestone as the official end of the emergency and as result relax their COVID-safe behaviours. Somehow paradoxically, vaccines might have, at least initially, a net negative effect. According to the health-belief model, one of the most commonly used psychological theories to characterize health-related behaviours, beliefs, perceptions, barriers to take action, and other modifying variables such as socio-demographic and socio-economic factors are key ingredients driving behavioural changes [2224]. A recent study may offer an empirical evidence of this tendency [25]. By means of surveys delivered before the vaccination rollout in December 2020, authors found that vaccine information reduce adherence to social distancing, hygiene measures, and the willingness to stay at home. Several surveys conducted during the COVID-19 pandemic, well before any concrete hope for a vaccine, confirm this picture, provide hints of how the arrival of vaccines might corrode even more adoption, highlight how compliance is a complex multi-faced problem [5] and that risk-perception as well as NPIs adoption are indeed associated to several socio-economic determinants such as age, gender, wealth, urban-rural divide [2641].

The literature aimed at estimating the epidemiological and societal impact of COVID-19 vaccines has been focused mainly on two very important points. The first line of research has been devoted to quantifying the effects of a vaccine on the evolution of the pandemic, considering different efficacy and coverage levels [42]. The second instead tackled the issue of vaccine allocation investigating strategies that target first different groups (i.e., age brackets, high-risk individuals) or particular occupations (i.e., doctors, nurses) [10, 4345]. Very recently, the intuition that social distancing remains key during vaccination rollout stimulated few studies on the effects of a vaccine on the adoption of NPIs in specific settings [4650]. An agent-based model applied to North Carolina showed that lifting NPIs during vaccines distribution would imply a substantial increase in infections and deaths [51]. Similarly, a data-driven model of SARS-CoV-2 transmission for China estimated that NPIs need to remain in place at least one year after the start of vaccination to avoid a generalized disease resurgence [52]. Despite these examples, the impact on real-world scenarios of the interplay of disease dynamics, behaviour change, and vaccines rollout is still largely unexplored. Furthermore, because of the novelty of the problem, also modeling frameworks able to characterize the possible behavioural dynamics linked to vaccinations are scarce.

To tackle such limitations, we introduce an age-structured compartmental epidemic model capturing the possible relaxation of NPIs adoption in response to the vaccine rollout. We model different compliance levels as distinct compartments and consider different behavioural dynamics driving the relaxation of NPIs. We test extensively the effects of behaviour change on disease spreading for different prioritization strategies, vaccine efficacy, and vaccination rollout speeds, using real demographic data and contacts matrices for six countries: Egypt, Peru, Serbia, Ukraine, Canada, and Italy. We choose these countries sampling levels of economic development. Indeed, in the World Economic Situation and Prospects 2020 issued by the United Nations, Egypt and Peru are classified as Developing Economies [53], Serbia and Ukraine as Economies in Transition, Canada and Italy as Developed Economies. Furthermore, considering dissimilar countries allows us to explore the interplay between vaccination and behaviours also as a function of population pyramids and intra/inter-generational mixing observed around the world. High income countries are typically characterized by higher average age, but lower inter-generational interactions respects to mid/low income countries [54]. These observations, together with the dependence of COVID-19 fatality rates on age, point to the possibility of non-trivial dependencies which we aim to explore here. As a way to realistically account for the different epidemic trajectories, we also explore the model after calibrating it on COVID-19 weekly deaths in the period 2020/09/01–2020/12/31 in these countries. To this end, we incorporate the timeline and effects of governmental restrictions on social contacts.

Our results provide quantitative insights on the interaction between sustained NPIs and an effective vaccination program. We show that an early relaxation of COVID-safe behaviours may lower, and even nullify, the advantages brought by the vaccine in the short term. Overall, the picture that emerges from the analysis of the different countries is consistent: a high level of compliance towards NPIs such as mask-wearing, social distancing, and avoidance of large gatherings, is needed in order to avoid spoiling the great effort of the vaccination campaigns.

Results

Vaccine-behaviour model

We consider an age-structured epidemic model based on a Susceptible-Latent-Infectious-Recovered compartmentalization with the addition of pre-symptomatic and asymptomatic infection stages and deaths. On top of the disease dynamics, we model both the vaccination rollout and the behavioural responses linked to it. After the start of the vaccination rollout, on each day a fraction of the population receives a vaccine that decreases both the probability of infection (with efficacy VES) and of developing symptoms (with efficacy VESymp) and enters in the compartment V [44]. As result, in our simulations, the overall efficacy of vaccine against severe outcomes such as death is VE = 1 − (1 − VES)(1 − VESymp). We introduce the rollout speed rV as the number of daily administered vaccine doses expressed as a percentage of the total population.

We study and compare three different strategies of vaccine prioritization. The first considers distributing vaccines in decreasing order of age. Previous research on COVID-19 vaccines allocation has shown that this strategy is the most effective in reducing the number of deaths and overall severity [10, 43, 44]. It is also the main strategy currently deployed in the different vaccination campaigns across the globe. In the second and third strategy vaccines are either distributed homogeneously or first to the individuals in age brackets 20–49 [55] and then to the rest of the population. The last strategy seeks to reduce symptomatic transmission targeting the most social active part of the population [10]. For simplicity, we will refer to the three rollout approaches as vaccine strategy 1, vaccine strategy 2, and vaccine strategy 3. In parallel to the vaccination, individuals (both susceptible S and vaccinated V) may start giving up safe behaviours and expose themselves to higher infection risks. We introduce the compartments SNC and VNC for susceptible and vaccinated individuals respectively, to describe different behavioural classes (for convenience “NC” refers to “non-compliant” to COVID-safe behaviours). The increased infection risk for these individuals is captured by the parameter r > 1 (for example, r = 1.3 indicates an increased risk of 30%). Sensitivity analysis on parameter r has been performed for values in the range of current estimates, as reported in Materials and methods section). We model the transition from S and V towards riskier behavioural classes as a function of the fraction of the vaccinated population and the parameter α. In turn, we imagine that a worsening of the epidemiological conditions may push non-compliant individuals back to safer behaviours. Here, we consider the number of fatalities per 100, 000 individuals in the previous time step (i.e., day) and the parameter γ to control the second behavioural transition. We refer the reader to the Materials and Methods section and the S1 Text for more details on the model. We refer to the the S1 Text, for results related to different behavioural mechanisms: a simpler one with constant rate, and a more complex one with different rates for vaccinated vs. unvaccinated individuals.

Interplay between NPIs adoption and vaccination campaign

We use data for six different countries: Egypt, Peru, Serbia, Ukraine, Canada, and Italy. As mentioned above, these countries have been selected to sample a range of demographics and socio-economic contexts. In Fig 1A, we represent some key characteristics of the demographic and the mixing patterns between age groups of these countries. First, we consider the fraction of people aged over 65. This is a measure of the epidemic fragility of a country: indeed 65+ individuals are at particularly high risk of mortality from COVID-19 disease [56]. Italy is the country showing the highest fraction of 65+ people (24%), while the developing economies show the lowest fraction (Egypt 5.3% and Peru 8.7%). Indeed, low/middle-income countries tend to have a younger population with respect to the high-income ones. Second, we present also the contacts intensity of different age groups. This is defined as the total number of contacts that an individual in a certain age group has, on average, with all the others in a day. We observe a typical decreasing trend, with younger people that tend to have more contacts. In particular, developing economies show a much higher number of daily contacts for individuals aged under 30. This can play an important role in the spreading even in developed economies. According to a recent study the resurgence of the COVID-19 Pandemic in the US after Summer 2020 was mainly sustained by younger people [55]. Finally, we represent a measure of inter-generational mixing. We define it simply as the number of daily contacts that an individual in the age groups at high mortality risk from COVID-19 (65+) receives from individuals in the age brackets 0–49. We observe that Egypt is the country showing the highest inter-generational mixing, followed closely by Peru.

thumbnail
Fig 1. Demographic, mixing patterns, and relative deaths difference for different rollout speeds and prioritization strategies.

A) Demographic characteristics and mixing patterns of the six countries considered are shown: percentage of the population aged over 65, contacts intensity for different age groups, and a measure of inter-generational mixing. B) Relative deaths difference is computed as the fraction of deaths that are avoided with a vaccine with respect to a baseline simulation without vaccine. We display results of the simulations for three vaccine rollout speed and prioritization strategies. Other parameters used are γ = 0.5, R0 = 1.15, r = 1.3, VES = 70% and VESymp such that VE = 90%, 0.5% of initially infected, 10% of initially immune individuals, and simulations length is set to 1 year.

https://doi.org/10.1371/journal.pcbi.1009346.g001

It is important to notice how heterogeneities in health infrastructures, access to health care, and comorbities are expected across the six countries under examination [54]. These features might induce variations in the IFR. For simplicity however, in the following, we use the same fatality rates for all countries. We leave the exploration of this important aspect to future work.

Prioritization strategies, rollout speed, and vaccine efficacy.

The proposed model allows to investigate the impact of the behavioural response under different conditions. We consider six populations matching the characteristics of the countries under consideration and exploring the model with the same epidemic initial conditions. The experimental setup considers that each population has already experienced a previous wave of infections and that restrictive measures are in place to mitigate the spreading. Therefore, we set the basic reproductive number R0 = 1.15, 0.5% of initially infected individuals, and 10% of immune individuals. In line with estimates of vaccine efficacy against COVID-19 we set VES = 70% and we choose VESymp such that VE = 90% [57], while we let vary the vaccine rollout speed rV between 0.1% and 1% to cover the spectrum of real vaccination rollout speeds of the vaccination campaigns across the globe. Peru, for example, administered on average 0.05 daily doses per 100 people in the week commencing on the 8th of March, Italy administered 0.30, and Serbia 0.70 [58]. In Fig 1B we show the relative deaths difference for the three vaccine prioritization strategies and the three vaccine rollout speed rV = 0.1%, 0.5%, 1.0%. The relative deaths difference (RDD) is defined as the fraction of deaths averted thanks to the vaccine with respect to a baseline simulation without vaccine and, therefore, no behavioural response. For example, RDD = 0.2 indicates that 20% of deaths are averted. Note how this quantity may become negative if the behavioural response causes more deaths than those averted thanks to the vaccine. We refer the reader to the Materials and Methods section for further details. Moreover, we explore a range of behavioural responses running the simulations for a range of values of the parameter α. Starting from α = 0 (which implies no behavioural response), we perform simulations with increasing α values (implying stronger reactions), while keeping constant the other behavioural parameter (γ = 0.5).

As a first observation, across the different countries and rollout speeds considered, the strategy aiming to curb the severity of the pandemic (i.e., vaccine strategy 1) is indeed the most effective in reducing the number of deaths.

As illustrative and concrete example let us consider the case of Canada. When α = 0 (i.e., no behaviour response) and rV = 1% with this strategy the fraction of averted deaths is 0.77 with respect to a baseline without vaccine. This fraction reduces to 0.72 with the homogeneous vaccination strategy, and to 0.70 with the strategy prioritizing younger individuals. When rV = 0.5%, these fractions reduce to, respectively, 0.67, 0.64, 0.62. The ordering of the strategies slightly changes when rV = 0.1%, in this case we obtain 0.28. 0.35, 0.36. Overall, the picture that emerges is consistent with previous studies on the allocation of COVID-19 vaccines among age groups showing that the strategy minimizing deaths is generally the one targeting the elderly, even though for some combination of the parameters the strategy targeting the younger might be preferable [43, 44].

When the behavioural response to the vaccine rollout is considered (i.e., α > 0), we note a consistent decreasing trend of the relative deaths difference for increasing values of α. This shows that the behavioural response impacts the unfolding of the epidemic, and that a relaxation of NPIs leads to a smaller fraction of averted deaths. Interestingly, a concerning effect also emerges: in some conditions, as non-compliance becomes larger, the benefit brought by the vaccine is nullified and the number of observed deaths increases with respect to the no-vaccine and no-behaviour change scenario (i.e., the relative deaths difference become negative). This is solely attributable to the behavioural reaction to the vaccination campaign which in turn is not efficient enough to balance behaviour relaxation. Indeed, this phenomenon is observed in particular when strategies that do not target a reduction in severity are employed, and when the vaccine rollout speed is low. In the case of Serbia, with rV = 1% and a strategy aimed at reducing severity, the fraction of averted deaths goes from 0.84 when α = 0 to 0.73 when α = 10, with a potential loss of 0.11 attributable to the NPIs relaxation. This effect is more pronounced in the case of the vaccination strategy 2 (or 3): the fraction of averted deaths goes in this case from 0.76 (0.73) when α = 0 to 0.59 (0.50) when α = 10, with a potential loss of 0.17 (0.23) attributable to the NPIs relaxation. Similarly, lower vaccine rollout speed are impacted more significantly by behavioural responses. In the previous example, an immunization campaign with a rV of 1%, 0.5%, and 0.1% would induce a relative deaths difference, in the case of the vaccination strategy 1, of 0.84, 0.75, 0.40 when α = 0. By setting instead α = 10, these figures would drop to 0.73, 0.44, −0.67 with a loss of 0.11 in the first, of 0.31 in the second, and of 1.07 in the third case.

In the results presented in Fig 1B, we observe that the most disadvantaged countries are the developing economies, Egypt and Peru. At first, this may seem counter-intuitive. Indeed, these countries have a lower fraction of individuals aged over 65. A possible explanation is given by the high activity of young individuals combined with the high inter-generational mixing. Moreover, it is worth stressing that these results are to be intended in relative terms: a relative worst performance in averting deaths does not necessarily imply a worst absolute performance. In other words, the relative impact of behavioral responses might be stronger in those countries, but their age pyramid might induce still a smaller number of deaths in absolute terms. On this point, it is important to remember that we are not considering possible modulations of the IFRs due to the heterogeneity in health care infrastructure, prevalence of co-morbidities and access to health care. In the S1 Text we provide an example highlighting this point by comparing Italy and Egypt. We observe also that the difference between countries is more profound for the vaccination strategy aimed at reducing transmission (i.e., strategy 3). Furthermore, behavioural relaxation widens the distance between the countries. With rV = 0.1% and vaccination strategy 3, the gap between the relative deaths difference for Italy and Egypt is 0.24 when α = 0. This figure increases to 0.87 when α = 10.

In the S1 Text we repeat the analysis considering the fraction of averted infections instead of deaths. We find that the ranking of allocation strategies is inverted. Indeed, in line with previous findings [43], when considering averted infections the most efficient strategy is the one targeting the younger population first, followed by the homogeneous, and the strategy aimed at curbing severity. In the S1 Text, we also repeat the results presented in Fig 1B exploring different vaccine efficacy VE (50%, 70%, 90%) rather than rollout speed. Similarly to the case presented in the main text, also in these additional analyses we find that early behavioural relaxation reduces the fraction of averted infections, influences heterogeneously the prioritization strategies, and impacts more significantly lower vaccine efficacy or slower rollout.

To further explore how to mitigate the impact of behavioural relaxation once that the vaccination campaign started, in Fig 2 we systematically explored the interplay between vaccine efficacy VE and vaccination rollout speed rV. The black dashed lines highlight the combinations of VE and rV that achieve a 30% drop in observed deaths (taken as reference value), in the absence of a behavioural response (i.e., α = 0). In most countries, it can be achieved with rV smaller than 0.2%. On the contrary, when even a mild behavioural response is active (red dash-dotted lines, α = 1), the rollout speed has to increase greatly when vaccine efficacy diminishes: for the case of Italy, a 30% drop in deaths would be achieved with a rV of 0.4% when the vaccine efficacy is 90%, but in the case of VE = 60% the rollout speed has to grow up to 0.8% to achieve the same result.

thumbnail
Fig 2. Interplay between vaccine efficacy and rollout speed.

Contourplot of the relative deaths difference in the scenario with a mild behavioural response (α = 1) for different combinations of VE and rV. We let VES vary between 30% and 70% and we choose different VESymp such that the overall efficacy VE vary between 50% and 90%. A 30% reduction of deaths is highlighted with a red dashdotted line. The black dashed line highlight the 30% death drop achieved with a vaccination campaign without the behavioural component (α = 0). Vaccination Strategy 1 is considered, parameters used are γ = 0.5, R0 = 1.15, r = 1.3, 0.5% of initially infected, 10% of initially immune individuals, and simulations length is set to 1 year.

https://doi.org/10.1371/journal.pcbi.1009346.g002

Epidemic trajectories and impact of behavioural relaxation.

As a way to ground the model on more realistic epidemic scenarios, we calibrate it to the real epidemic trajectories of the six countries considered. The calibration is performed via an Approximate Bayesian Computation technique (ABC) [59] on weekly deaths for the period 2020/09/01–2020/12/31. In doing so, we estimate some key epidemiological parameters matching the context of each country, thus defining realistic initial conditions for the model. To capture the effects of non-pharmaceutical interventions on the contacts between people we consider data from the Google Mobility Report [60] and the Oxford Coronavirus Government Response Tracker [61] until week 11 of 2021. After the calibration step, the model evolves from 2021/01/01 up to 2021/06/01 with the vaccination rollout for each of the countries. Note that in some of these countries, the vaccination campaigns officially started in the second half of December, though mostly symbolically.

In Fig 3A we show key indicators resulting from the calibration. In particular, for the six countries we report the boxplot for the calibrated infection parameter β, the projected number of symptomatic infections per 100, 000, and the fraction of recovered individuals on the 2021/01/01, which is the start of the vaccination campaign in our simulations. We acknowledge a significantly high acquired immunity in the case of Peru, where the estimated attack rate as of 1st January, 2021, is around 38%. Nonetheless, this figure is line with other available estimates [58]. We observe how the calibration allows for an heterogeneous representation of the epidemiological conditions of the different countries. In Fig 3A we also show the effects of restrictive measures on contacts. As a proxy, we consider the ratio between the leading eigenvalue of the contacts matrix (considering the restrictions) and the leading eigenvalue of the baseline contacts matrix (without restrictions). The leading eigenvalue of the contacts matrix influences the reproductive number of the disease (see Materials and methods), and more broadly it measures the intensity of contacts among people. By normalizing with respect to the baseline contacts matrix (with no restrictions) we can grasp the strictness of the measures in place. For example, a ratio of 0.3 would imply, in our simulations, a 70% reduction of the reproductive number with respect to the baseline without restrictions. The effect of restrictive measures on contacts varies over time up to week 11, 2021. In the case of Italy, for example, we observe the partial ease of the measures during January and February (i.e., ρ(Ct)/ρ(Cb) increases) followed by the tightening of measures in March aimed at curbing the third wave of infections. Afterwards, we keep the mixing levels as those observed for week 11 (more details in the Materials and methods section). This is a conservative assumption as we can imagine that the seasonality and impact of the vaccination campaign might induce a relaxation of NPIs as the number of cases and deaths, hopefully, will go down. In this scenario, the countries with the strictest measures in place after week 11, are Peru and Italy, and Canada, while Egypt and Ukraine show the most permissive measures. In the S1 Text, we repeat the analysis using real data to compute contacts reduction parameters for the whole period. The results are in line with the observations presented below.

thumbnail
Fig 3. Giving up NPIs during rollout may nullify the benefits brought by the vaccine.

A) We display for the different countries the boxplot of the calibrated infection parameter β, the projected number of symptomatic infectious cases per 100′000 and the fraction of recovered at the 2021/01/01, start of the vaccination campaign in our simulation. We also report the ratio between the leading eigenvalue of the contacts matrix considering restrictions and of the baseline contacts matrix with no restrictions. B) We display the median relative deaths difference for the calibrated model in the different countries. We consider the three vaccination strategies and two possible rollout speed: rV = 1% (faster rollout), and rV = 0.25% (slower rollout). We run the model over the period 2021/01/01–2021/06/01. Other parameters are γ = 0.5, α = 102, r = 1.3, VES = 70% (VE = 90%).

https://doi.org/10.1371/journal.pcbi.1009346.g003

In Fig 3B we study, for the calibrated model in the six countries, the difference between the fraction of averted deaths with respect to a baseline without vaccine in the case of a strong behavioural reaction to vaccine deployment (i.e., α = 102) and in the case of no relaxation of NPIs (i.e., α = 0). This quantity, which we indicate with ΔRDD, represents the additional fraction of deaths that occur due to lower compliance to COVID-safe behaviours in our simulations (see Materials and methods for more details). We study two vaccination rates. A fast deployment that manages to vaccinate rV = 1% of the population daily (aligned with the preparedness plan for Influenza pandemic [62] and vaccination rates achieved by countries like Israel or Chile), and a slower rollout rV = 0.25%. Worryingly, in the different countries, even in a scenario with several restrictions and a successful vaccination campaign in place, it is possible to witness a significant increase in deaths exclusively due to a relaxation of COVID-safe behaviours. As a result of the calibration step, the countries under investigation face different epidemiological conditions (e.g., effective reproductive number, immunity from previous waves, initial number of infected, estimated effects of restrictions). These translate into phenomenological differences with respect to the results of the previous analyses. In particular, Ukraine and Egypt are significantly impacted by behavioural relaxation in the case of the slower rollout rV = 0.25%, while Peru, differently from the previous analysis, is less impacted. From Fig 3A we see that in the first two countries the estimated effect of restrictions on contacts (especially in the last weeks) is smaller with respect to the others. On the other hand, Peru features tougher NPIs, together with Italy and Canada, which are also less impacted by behaviour in Fig 3B. This underlines how, besides the characteristics of the vaccination campaign, also the epidemiological conditions and the measures in place are important factors influencing the behaviour-vaccine interplay. In the case of the faster rollout rV = 1%, the obtained values of ΔRDD are smaller and the differences between countries are less pronounced. A possible explanation is that when a fast rollout is employed the influence of epidemiological conditions and restrictions become less important because of the efficiency of vaccine administration. Furthermore, we observe that, for both rV, the vaccination strategy targeting a reduction of severity (strategy 1) performs generally better in case of NPIs relaxation with respect to the other strategies. In the S1 Text we investigate more in detail the comparison between vaccination strategies in terms of robustness to behaviour relaxation considering both averted deaths and infections. We find that, when considering deaths, a vaccine allocation aimed at reducing disease severity (i.e., strategy 1) is always preferable. Instead, vaccination strategy 2 and 3 perform better when the number of avoided infections is used as evaluation metric. Information about vaccination strategy robustness against behavioural changes might be used to tune and design resilient rollout campaigns. Policy-makers might also consider to optimize vaccination strategies with respect to a combination of multiple metrics. In this direction, our analysis provide additional insights, and suggests that behavioural changes might play an important role possibly modifying the impact of vaccine prioritization strategies.

Discussion

For almost a year, in the midst of a global pandemic, policymakers struggled to implement sustainable restrictions to slow SARS-CoV-2 spreading. Every non-pharmaceutical intervention was aimed at slowing (in few countries stopping) the disease progression buying time for the development, test, production and distribution of vaccines that might ultimately protect the population. With an impressive scientific endeavor, several vaccines have been developed and an early distribution campaign was rolled out by the last days of December 2020. Besides the potential threats emerging from new virus strains [63], the current vaccination campaign represents the beginning of a new normal and a gigantic step towards complete virus suppression. However, we demonstrated that if the growth in vaccination uptake would lead to overconfident conducts inducing relaxation of COVID-safe behaviours, additional avoidable deaths will occur. We extended the literature proposing a mechanistic compartmental model able to simulate the unfolding of COVID-19, the vaccination dynamics and the compliance/non-compliance transition modulated by different behavioural mechanisms. Performing in-silico simulations allowed us to explore, from a theoretical standpoint, the interplay among different vaccination strategies, rollout speeds, vaccine efficacy, and behavioural responses. We found that behaviour impacts non-linearly vaccine effectiveness. Indeed, early NPIs relaxation affects more significantly slower rollouts, lower vaccine efficacy, and allocation strategies that target reduction in transmission rather in severity. We included in our analysis six different countries (Egypt, Peru, Serbia, Ukraine, Canada, and Italy), representatives of various points of the spectrum of world economies. This allowed us to observe the effects of behaviour-vaccine relationships for various population pyramids and mixing patterns. We observed that the developing economies, characterized by a younger population, but higher contacts activity and inter-generational mixing, are generally more affected by behaviour change with respect to developed economies. Then, as a way to ground the model on more realistic epidemiological conditions, we calibrated it using real epidemic and mobility data for the six countries considered and simulated the unfolding of the first months of the vaccination campaign. In such realistic scenario, we observed that even with restrictive measures in place and a successful vaccination campaign, it is possible to witness to non-negligible increases in deaths attributable to an early relaxation of COVID-safe behaviours. The calibration step allowed us to highlight that also the epidemiological conditions related to the country-specific unfolding of the disease are an important factor influencing the interplay behaviour-vaccine.

We acknowledge some limitations in the present study. First, we considered the vaccines fully working immediately after the first dose and we neglected that vaccination campaigns are using a portfolio of vaccines rather than a single one. We have also studied three simple vaccination strategies that neglect the complexities of an unprecedented mass vaccination. As result, both the vaccination priorities, vaccine effects and vaccination rates are an approximation of reality. In the S1 Text, we studied a data-driven vaccination rollout for Italy, where vaccines are distributed to the various age brackets following the real daily administration data. The results we obtained are qualitatively similar to those presented in the main text. Second, while the model calibration suggests that our approach can nicely capture national trends, the model is not meant to provide accurate forecasts of the local unfolding of the disease, but rather to test what-if scenarios in a comparative fashion. We have considered a simple age-structure compartmental model that does not capture spatio-temporal heterogeneity both in terms of spreading and of NPIs implementation which have instead been observed in the countries under investigation. Third, our model does not include new, more transmissible, variants of concern and assumes the same IFR across the six countries. Finally, we propose and model two potential mechanisms leading to behavioural changes, but data are not available to perform a quantitative validation of the behavioural components of the model.

Implementation of individual protective behaviours and adherence to NPIs have been vital in order to reduce the transmission of SARS-CoV-2 leading to substantial population-level effects [5, 6475]. Behavioural science can provide valuable insights for managing policies, incentives, communication strategies and can help coordinate efforts to control threats and evaluate such interventions [76]. As during the first waves of COVID-19, when NPIs were the only available mitigation measures [77], the results of our paper call for adequate strategies to keep high the attention and compliance towards individual COVID-safe behaviours, such as mask-wearing, social distancing, and avoidance of large gatherings now that vaccines are finally available. Communication strategies and policies should keep targeting such non-pharmaceutical intervention to avoid frustrating the immense effort of the vaccination campaigns.

Materials and methods

Epidemic model

Model definition.

We propose an age-stratified compartmental model that incorporates both the vaccination process and the behaviour dynamics possibly linked to it. See Fig 4 for a schematic representation. Individuals are divided into 16 five-year age groups (last group is 75+). The virus transmission is modeled using the following approach. Susceptible individuals (S compartment) in contact with infectious compartments become latent (L compartment). After the latent period (ϵ−1), L individuals enter the pre-symptomatic stage of the infection (P). Then, they can transition either in the asymptomatic (A) or the symptomatic stage (I) at rate ω (the length of time spent in the L and P compartments is the incubation period). The probability of being asymptomatic is f. After the infectious period (μ−1), both I and A individuals enter the Recovered compartment (RI, RA). Alternatively, I individuals can also die and transit to the D compartment. Note how we include a delay of Δ days from the time individuals enter the compartment D and die (Do compartment, the superscript stands for “observed”). The infectious compartments are P, A, I. The transmission rate is β and the force of infection is dependent on the age-stratified contact matrix , whose element Cij represents the average number of contacts that an individual in age group i make with individuals in j per day. The matrix C has four location-specific contributions: contacts at home, workplace, school, and other locations. We adopt country-specific contacts matrices provided in Ref. [78]. We assume that P and A have lower infectiousness with respect to symptomatic I (βχ, with χ < 1). Similar approaches have been used in several modeling studies in the context of the current pandemic [70, 79].

thumbnail
Fig 4. Compartmental model.

We consider an extension of the classic SLIR model adding compartments for presymptomatic (P) and asymptomatic (A). We also design compartments for vaccinated (V), dead (D describes individuals that will die with a delay of Δ−1 entering the compartment Do), susceptible (SNC) and vaccinated (VNC) individuals that do not comply with COVID-safe behaviours. The vaccine offers a protection VES against infection and VESymp against symptoms that can lead to severe outcomes such as death. The transmission rate for susceptible is β and for susceptible non-compliant (r > 1). The parameter α regulates the transition from compliant to non-compliant behaviours, while γ regulates the opposite flow. Arrows describe the transitions between compartments. For simplicity of visualization, we do not display the compartments mediating the different transitions. For example, the transition SL is mediated by the infectious compartments (P, A, I, PV, IV, AV). The compartmentalization is then extended to account for empirical age-structure and contact matrices.

https://doi.org/10.1371/journal.pcbi.1009346.g004

On top of the disease dynamics, we model both the vaccination process and the behavioural change that is possibly coupled to the vaccination. More in detail, after the start of the vaccination campaign at tV, at each time step, a fraction of the susceptible population receive a vaccine and transit to compartment V. We introduce the rollout speed rV as the number of daily available doses expressed as a percentage of the total population. We adjust the number of doses available per day considering the fraction of susceptibles among all those that can receive the vaccine (S, L, P, A, and RA) and we assume that the remaining doses are wasted. We consider a “leaky” vaccine that reduces susceptibility with a certain efficacy VES and the probability of developing symptoms by VESymp. In other words, the infection rate for V individuals is β(1 − VES) and the probability of entering the infectious symptomatic compartment IV from PV is (1 − f)(1 − VESymp). Similar approaches have been used in previous works in the context of mathematical modeling of COVID-19 immunization campaigns [44, 80]. As mentioned above, we consider three vaccination strategies: one in which the vaccine is given in decreasing order of age (vaccination strategy 1 aimed at reducing the severity), one in which it is given homogeneously to the population (vaccination strategy 2), and one in which it is first given to individuals in age brackets 20–49 and then to the rest of the population (vaccination strategy 3 aimed at reducing the transmission of the virus) [10]. In parallel, we imagine a behavioural dynamics triggered by the presence of the vaccine. Indeed, susceptible individuals (both vaccinated and not) may start adopting less safe behaviours because reassured by the presence of an effective vaccine. This is encoded in the model with a transition from the compartment S (V) to a new compartment SNC (VNC)—NC stands for non-compliant—of individuals that protect less themselves and as a result get infected at a higher rate. The parameter r > 1 captures the increased risk of contagion for NC individuals. It modulates the contribution to the force of infection from the SNC and VNC compartments. In the simulations presented in the main text we used r = 1.3 imagining a scenario where the relaxation of COVID safe behaviors increases the transmission rate of 30%. This choice is informed by the estimated effects of NPIs such as face masks and social distancing on COVID-19 spreading, as reported in the literature [4, 20, 81, 82]. Nonetheless, in the S1 Text we perform a sensitivity analysis considering also r = 1.1, 1.5. We propose two mechanisms to model the behavioural transition. In the first mechanism, the transition towards non-compliance happens at a rate α and it is catalyzed by the cumulative fraction of individuals that received a vaccine (vt, including both compliant and non-compliant). The opposite transition from SNC (VNC) to S (V) happens at rate γ and is catalyzed by the number of deaths per 100, 000 observed in the previous time step (, including both compliant and non-compliant). Indeed, an increase in deaths is frequently used—especially by media—as an indicator of the severity of the current epidemiological situation. Existing literature suggests that risk perception (in the form of number of infected individuals or deaths) and communication of such risk significantly affect adherence to personal mitigation strategies such as social distancing and wearing face masks [5, 83]. Expressing the number of deaths in proportion to the population allows us to compare countries of different size. Overall, this approach aims to depict the adaptive nature of human behaviour where individual choices are influenced by vaccination and pandemic progression [84]. In the second mechanism, S (V) individuals transit to the non-compliant compartment SNC (VNC) at a constant rate α. We also account for the possibility of non-compliant going back to safer behaviours again at a constant rate γ. To simplify the narrative, we present the results considering this second mechanism only in the S1 Text. The overall picture discussed above does not change significantly. Note that in order to avoid issues with transition probabilities large than one we model the rates as λXXNC = 1 − expg(α), λXNCX = 1 − exph(γ) with X = [S, V] and where the specific expression of the exponent depend on the two mechanisms described above: g(α) = αvt, in the first one, and g(α) = α, h(γ) = γ in the second mechanism. Note how for small values of g(α), h(γ) the rate converges to the usual mass-action law. The model just described can be written down as the following system of differential equations for individuals in age group k: (1) Where is the force of infection for age group k.

The basic reproduction number is , where is the contacts matrix weighted by the relative population in different age groups. In the S1 Text we provide details on its derivation. We adopt the model integrating numerically the equations, thus it is deterministic. However, it is worth stressing that when we calibrate the model to the real epidemic trajectory in the six countries, we use a probabilistic framework through an Approximate Bayesian Computation technique. Said differently, the calibrated parameters are characterized by posterior probability distributions rather than exact values. For this reason, the results presented in Fig 3 (i.e., the median of model’s projections) are to be intended as an ensemble of multiple trajectories generated sampling from the posterior distribution of the fitted parameters. The model is implemented in the programming language python with the use of the libraries scipy [85], numpy [86], and numba [87]. The visualization of the results is realized with the library matplotlib [88].

In the Results section, we kept for simplicity γ = 0.5 and we let vary α. The choice of γ is informed by the maximum number of deaths observed in the countries of focus. In Italy, for example, the maximum number of deaths reported on a single day is around 1, 000. Therefore, this value of γ is such that, in a similar situation, non-compliant individuals would likely return to COVID-safe behaviours. As described in the previous paragraph, in the dynamic-rate behavioural mechanism we model the transition rate as (X = [S, V]), where is the number of deaths per 100, 000 individuals observed at time t − 1. Hence at the peak of deaths the transition rate towards compliance is . The sensitivity to this choice is discussed in the S1 Text, where we also report additional analyses to better understand the effects of the behavioural parameters. More in detail, we plot the rates and together with the actual number of NC individuals in time for several values of α, γ, vt, and different epidemiological conditions.

In Table 1 we report a list of the model’s parameters together with their values used in the simulations and the related sources. We also indicate which values are optimized during the calibration step (described in the next sections).

Relative deaths difference.

In the previous analyses, we compared the vaccination strategies in the six countries using different evaluation metrics. The relative deaths difference (RDD) represents the fraction of deaths that are avoided in a simulation with vaccines (and behavioural response) with respect to a baseline simulation without vaccines (and no behavioural response). We compute it as: Where we made explicit the dependence from α. By definition, a value of RDD = 0.2 indicates that 20% less deaths are observed in the simulation with vaccines administered. Negative values of RDD capture instead scenarios where the behavioral response induces more deaths that in the case without vaccines. An analogous quantity can be easily computed for infections. In Fig 3B we have also considered the difference between the fraction of averted deaths with respect to a baseline without vaccine in the case of behavioural reaction (i.e., α > 0) and in the case of no relaxation of NPIs (i.e., α = 0). We indicated this quantity with ΔRDD(α), and we compute it as follows:

With respect to the simple RDD, the ΔRDD provides a slightly different piece of information. Indeed, by subtracting the RDD obtained with α = 0 we isolate only the effects of behaviour. As a clarifying example, ΔRDD(α) ∼ 0 refers to a negligible effect of behaviour on the fraction of averted deaths. Instead, ΔRDD(α) < 0, indicates that due to behaviour relaxation more deaths are occurring. Also in the case of ΔRDD(α), an analogous for infections can be easily derived.

Vaccination strategies

We consider three vaccination strategies. In the first one, the vaccine is given in decreasing order of age. Since the IFR of COVID-19 strongly correlates with age, many countries worldwide are adopting similar strategies prioritizing the vaccination among the elderly. Previous modeling works in the context of COVID-19 showed that this is the preferable strategy when considering the number of averted deaths [10, 43, 44]. In practice, this means that in our simulations we start giving the vaccine to the 75+ age bracket and we proceed in decreasing order only when everyone in this group is vaccinated. In the second strategy the vaccine is given homogeneously to the population respecting the age distribution. This means that, of the X vaccines available at step t, the fraction Nk/N is given to age group k (where Nk is the number of individuals in age bracket k, and N is the total number of individuals). In the third strategy, we prioritize age groups 20–49. Indeed, individuals at high mortality risk from COVID-19 disease may be protected indirectly by vaccinating age brackets that sustain the transmission [98, 99]. In practice, in our simulations we start giving the vaccine homogeneously to the 20–49 age brackets, and when they are all vaccinated we give it homogeneously to all other groups. In the previous two strategies, if, on a given time step, the number of people remained in the age group which is currently being vaccinated is smaller than the total number of vaccines available, the exceeding part is given to the next age group in line. This implies that the number of vaccines given in different time steps is always constant. Since doses may be administered to non-susceptible, on each day we compute the fraction of S individuals among all those who may have received a vaccine (latent, presymptomatic, infectious asymptomatic, and recovered asymptomatic individuals) and we adjust accordingly the number of doses available. A similar approach has been used in Refs. [44, 80]. Though a simplification of reality this approach allows to avoid making assumptions about the possible modulations of vaccine efficacy for individuals vaccinated when not susceptible. In the S1 Text, we repeat part of the analyses considering that some individuals may decide not to get vaccinated. We find that, despite shifts due to vaccine hesitancy, the main findings remain unchanged.

Model calibration

In this section we describe the methods used to calibrate the model to the real epidemic trajectories in the countries considered. We use an Approximate Bayesian Computation technique to calibrate the model on weekly deaths during the period 2020/09/01–2020/12/31. We account for government-mandated restrictions and their effect on the contacts between individuals modifying the contacts matrices using data from the Google Mobility Report [60] and the Oxford Coronavirus Government Response Tracker [61].

Approximate Bayesian computation.

We calibrate the model for each country using the Approximate Bayesian Computation (ABC) rejection method [59, 100]. At each step of the rejection algorithm, a set of parameters θ is sampled from a prior distribution and an instance of the model is generated using θ. Then, an output quantity E′ of the model is compared to the corresponding real quantity E using a distance measure s(E′, E). If this distance is greater (smaller) than a predefined tolerance ψ, then the sampled set of parameters is discarded (retained). After accepting N sets, the iteration stops and the distribution of accepted parameters is an approximation of the real posterior distribution P(θ, E). The free parameters of our model and the related prior distributions are:

  • The transmission rate β. We explore uniformly values of β such that the related R0 is between 0.8 and 2.2. The basic reproduction number of SARS-CoV-2 is higher [101], but we consider lower values since our calibration starts on 2020/09/01 when restrictions were in place to mitigate the spreading.
  • The delay in deaths Δ ∼ [14, 25]. Indeed, for COVID-19 the average time between symptoms onset and death is about 2 weeks [97] and we also account for possible additional delays in death reporting.
  • The initial number of infected individuals. We explore uniformly values between 0.5 and 15 times the number of reported cases in the week before the start of the simulation. We then assign these individuals to the infected compartments (L, P, A, I) proportionally to the time spent there by individuals (ϵ−1 for L, ω−1 for P, and μ−1 for I, A), and we split between I and A individuals considering the fraction of asymptomatic f.

The model is calibrated on the period 2020/09/01–2020/12/31 using the weighted mean absolute percentage error on weekly deaths as an error metric with a tolerance ψ = 0.3 (ψ = 0.4 for Egypt for convergence issues) and 5, 000 accepted parameters set.

Model initialization.

The number of individuals in different age groups is initialized considering the 2019 United Nation World Population Prospects [102]. We consider 16 age brackets of five-years, except for the last one that includes 75+ individuals. As specified above, the initial number of infected individuals is calibrated considering the total number of confirmed cases in the week before the start of the simulation according to the data issued from the European Centre for Disease Prevention and Control [103]. To initialize the number of non-susceptible individuals (placed in the R compartment) we compute the average of several publicly available projections of total COVID-19 infections up to 2020/09/01 (i.e., start of the simulation) from different modeling approaches [58]. Both the initial number of infected and of non-susceptibles is assigned homogeneously across age groups. Other parameters used are ϵ−1 = 3.7 days, μ−1 = 2.5 days, ω−1 = 1.5 days χ = 0.55, f = 0.35 in line with current estimates of COVID-19 infection dynamics parameters [9092, 95, 96]. We use the age-stratified Infection Fatality Rate (IFR) from Ref. [56].

Modeling the effects of NPIs on contact matrices.

In our model, we incorporate the implementation of top-down NPIs by changing the contacts patterns defined by the contacts matrix C. As mentioned, we use the country-specific contacts matrices provided in Ref. [78]. These are made by four contributions: contacts that happen at school, workplace, home, and other locations. In a baseline scenario, the overall contacts matrix C is simply the sum of these four contributions. Following Ref. [104], we implement the reductions in contacts, due to the restrictions, multiplying the single contribution by a reduction factor ηi(t). Thus, in general the overall contacts matrix at time t become: (2)

For simplicity, we assume no changes to contacts at home, though lockdowns tend to increase them [64]. For the contacts locations work and other locations we use data from the Google Mobility Report [60]. In detail, we consider the fields workplaces percent change from baseline to model reduction in the contacts location work and the average of the fields retail and recreation percent change from baseline and transit stations percent change from baseline for other locations. A general entry of the report pl(t) represents the percentage change on day t of total visitors to a specific location l with respect to a pre-pandemic baseline. From pl(t) we derive a contacts reduction coefficient as follows: ηl(t) = (1 + pl(t)/100)2. Indeed, the number of potential contacts in a specific location is proportional to the square of the the number of visitors. We model contacts reduction in school considering instead the Oxford Coronavirus Government Response Tracker [61]. More in detail, we use the ordinal index C1 School closing. The index ranges from a minimum of 0 (no measures) to a maximum of 3 (require closing all levels). We turn this quantity into the contacts reduction coefficient in school as follows: ηs(t) = (3 − C1 School closing)/3. We use these datasets to inform contacts reductions up to week 11, 2021. After, we assume that contacts remain at the same level of week 11. Note how we let vary both the transmission rate β and the contact matrix. The former describes the risk of infection given contacts with infectious individuals. This is function of the disease (which is assumed to be the same, we don’t consider multiple or emergent new strains possibly more transmissible) and of the protective behaviours such as social distancing and use of face masks. The latter describes variations to the number and types of contacts induced by top-down NPIs as for example remote working, schools closure and lockdowns. By splitting the contributions to the force of infection of transmission rate and contact matrix we are able to take into consideration different behavioural attitudes which, given the same number of contacts, might lead to higher or lower risks of infection. This allows us to consider explicitly both top-down and bottom-up NPIs.

Calibration results.

In Fig 5, we report the results of the calibration. It is important to stress how our goal is not to develop a predictive model aimed at forecasting the pandemic trajectory. The fit is used to ground the model and to define the epidemic conditions at the start of the vaccination campaign in the six countries. In fact, our aim is to understand the possible interplay between behaviours and vaccine rollout which is also function of the epidemic progression. In the Figure we report the official and simulated weekly number of deaths (median, 50% and 95% confidence interval). Despite its simplicity and approximations, the model is able to reproduce the evolution of the pandemic in the six countries capturing well its progression after the summer. In the S1 Text we also report the posterior distributions for the parameters.

thumbnail
Fig 5. Calibration results.

For each country we represent the observed and simulated weekly deaths (median, 50% and 95% confidence intervals).

https://doi.org/10.1371/journal.pcbi.1009346.g005

Supporting information

S1 Text. Supplementary and sensitivity analysis to the modeling framework.

In this supplementary material we present additional analyses and extensions of our model. We also run extensive sensitivity analysis on the parameters used in the main text.

https://doi.org/10.1371/journal.pcbi.1009346.s001

(PDF)

Acknowledgments

All authors thank the High Performance Computing facilities at Greenwich University.

References

  1. 1. Desvars-Larrive A, Dervic E, Haug N, Niederkrotenthaler T, Chen J, Di Natale A, et al. A structured open dataset of government interventions in response to COVID-19. Scientific Data. 2020;7(1):285. pmid:32855430
  2. 2. Cowling BJ, Ali ST, Ng TW, Tsang TK, Li JC, Fong MW, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study. The Lancet Public Health. 2020;.
  3. 3. Flaxman S, Mishra S, Gandy A, Unwin HJT, Mellan TA, Coupland H, et al. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature. 2020;584(7820):257–261. pmid:32512579
  4. 4. Haug N, Geyrhofer L, Londei A, Dervic E, Desvars-Larrive A, Loreto V, et al. Ranking the effectiveness of worldwide COVID-19 government interventions. Nature Human Behavior. 2020;. pmid:33199859
  5. 5. Perra N. Non-pharmaceutical interventions during the COVID-19 pandemic: A review. Physics Reports. 2021; pmid:33612922
  6. 6. Seale H, Dyer CE, Abdi I, Rahman KM, Sun Y, Qureshi MO, et al. Improving the impact of non-pharmaceutical interventions during COVID-19: examining the factors that influence engagement and the impact on individuals. BMC Infectious Diseases. 2020;20(1):1–13. pmid:32807087
  7. 7. Pullano G, Di Domenico L, Sabbatini CE, Valdano E, Turbelin C, Debin M, et al. Underdetection of COVID-19 cases in France threatens epidemic control. Nature. 2020;. pmid:33348340
  8. 8. Callaway E. What Pfizer’s landmark COVID vaccine results mean for the pandemic. Nature. 2020;. pmid:33169027
  9. 9. Callaway E. COVID vaccine excitement builds as Moderna reports third positive result. Nature. 2020;. pmid:33199901
  10. 10. Wang W, Wu Q, Yang J, Dong K, Chen X, Bai X, et al. Global, regional, and national estimates of target population sizes for covid-19 vaccination: descriptive study. bmj. 2020;371. pmid:33323388
  11. 11. Britton T, Ball F, Trapman P. A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2. Science. 2020;369(6505):846–849. pmid:32576668
  12. 12. Read MC. EID: High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerg Infect Dis. 2020;26.
  13. 13. Lipsitch M, Dean NE. Understanding COVID-19 vaccine efficacy. Science. 2020;370(6518):763–765. pmid:33087460
  14. 14. Polack FP, Thomas SJ, Kitchin N, Absalon J, Gurtman A, Lockhart S, et al. Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine. New England Journal of Medicine. 2020;. pmid:33301246
  15. 15. Baden LR, El Sahly HM, Essink B, Kotloff K, Frey S, Novak R, et al. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. New England Journal of Medicine. 0;0(0):null. pmid:33378609
  16. 16. Dagan N, Barda N, Kepten E, Miron O, Perchik S, Katz MA, et al. BNT162b2 mRNA Covid-19 vaccine in a nationwide mass vaccination setting. New England Journal of Medicine. 2021;.
  17. 17. PHE monitoring of the early impact and effectiveness of COVID-19 vaccination in England; 2021. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/963532/COVID-19_vaccine_effectiveness_surveillance_report_February_2021_FINAL.pdf.
  18. 18. Lazarus JV, Ratzan SC, Palayew A, Gostin LO, Larson HJ, Rabin K, et al. A global survey of potential acceptance of a COVID-19 vaccine. Nature medicine. 2020; p. 1–4.
  19. 19. Peiris M, Leung GM. What can we expect from first-generation COVID-19 vaccines? Lancet (London, England). 2020;. pmid:32971042
  20. 20. Yang B, Huang AT, Garcia-Carreras B, Hart WE, Staid A, Hitchings MDT, et al. Effect of specific non-pharmaceutical intervention policies on SARS-CoV-2 transmission in the counties of the United States. Nature Communications. 2021;12(1):3560. pmid:34117244
  21. 21. Walensky RP, Walke HT, Fauci AS. SARS-CoV-2 Variants of Concern in the United States—Challenges and Opportunities. JAMA. 2021;. pmid:33595644
  22. 22. Rosenstock IM. The Health Belief Model and Preventive Health Behavior. Health Education Monographs. 1974;2(4):354–386.
  23. 23. Hochbaum GM. Public Participation in Medical Screening Programs: A Socio-psychological Study. Public Health Service publication. U.S. Department of Health, Education, and Welfare, Public Health Service, Bureau of State Services, Division of Special Health Services, Tuberculosis Program; 1958. Available from: https://books.google.it/books?id=YXEmHAAACAAJ.
  24. 24. Hayden J. Introduction to Health Behavior Theory. Jones & Bartlett Learning; 2013. Available from: https://books.google.it/books?id=9YZSAAAAQBAJ.
  25. 25. Andersson O, Campos-Mercade P, Meier A, Wengström E. Anticipation of COVID-19 Vaccines Reduces Social Distancing. Lund University, Department of Economics; 2020. 2020:29.
  26. 26. Kasting ML, Head KJ, Hartsock JA, Sturm L, Zimet GD. Public perceptions of the effectiveness of recommended non-pharmaceutical intervention behaviors to mitigate the spread of SARS-CoV-2. PloS one. 2020;15(11):e0241662. pmid:33147261
  27. 27. Wise T, Zbozinek TD, Michelini G, Hagan CC, Mobbs D. Changes in risk perception and self-reported protective behaviour during the first week of the COVID-19 pandemic in the United States. Royal Society open science. 2020;7(9):200742. pmid:33047037
  28. 28. Galasso V, Pons V, Profeta P, Becher M, Brouard S, Foucault M. Gender differences in COVID-19 attitudes and behavior: Panel evidence from eight countries. Proceedings of the National Academy of Sciences. 2020;117(44):27285–27291. pmid:33060298
  29. 29. Cvetković VM, Nikolić N, Nenadić UR, Öcal A, Noji EK, Zečević M. Preparedness and Preventive Behaviors for a Pandemic Disaster Caused by COVID-19 in Serbia. IJERPH. 2020;17(11):4124. pmid:32527056
  30. 30. Barber SJ, Kim H. COVID-19 Worries and Behavior Changes in Older and Younger Men and Women. The Journals of Gerontology: Series B. 2020;.
  31. 31. Lüdecke D, von dem Knesebeck O. Protective Behavior in Course of the COVID-19 Outbreak. Survey Results From Germany. Front Public Health. 2020;8. pmid:33072712
  32. 32. Muto K, Yamamoto I, Nagasu M, Tanaka M, Wada K. Japanese citizens behavioral changes and preparedness against COVID-19: An online survey during the early phase of the pandemic. PLoS ONE. 2020;15(6):e0234292. pmid:32525881
  33. 33. ul Haq S, Shahbaz P, Boz I. Knowledge, behavior and precautionary measures related to COVID-19 pandemic among the general public of Punjab province, Pakistan. The Journal of Infection in Developing Countries. 2020;14(08):823–835.
  34. 34. Wright AL, Sonin K, Driscoll J, Wilson J. Poverty and economic dislocation reduce compliance with COVID-19 shelter-in-place protocols. Journal of Economic Behavior and Organization. 2020;180:544–554. pmid:33100443
  35. 35. Lee M, Zhao J, Sun Q, Pan Y, Zhou W, Xiong C, et al. Human mobility trends during the early stage of the COVID-19 pandemic in the United States. PLoS One. 2020;15(11):e0241468. pmid:33166301
  36. 36. Weill JA, Stigler M, Deschenes O, Springborn MR. Social distancing responses to COVID-19 emergency declarations strongly differentiated by income. Proc Natl Acad Sci USA. 2020;117(33):19658–19660. pmid:32727905
  37. 37. Bonaccorsi G, Pierri F, Cinelli M, Flori A, Galeazzi A, Porcelli F, et al. Economic and social consequences of human mobility restrictions under COVID-19. Proc Natl Acad Sci USA. 2020;117(27):15530–15535. pmid:32554604
  38. 38. Gozzi N, Tizzoni M, Chinazzi M, Ferres L, Vespignani A, Perra N. Estimating the effect of social inequalities on the mitigation of COVID-19 across communities in Santiago de Chile. Nature Communications. 2021;12(1):2429. pmid:33893279
  39. 39. Chang S, Pierson E, Koh PW, Gerardin J, Redbird B, Grusky D, et al. Mobility network models of COVID-19 explain inequities and inform reopening. Nature. 2021;589(7840):82–87. pmid:33171481
  40. 40. Gauvin L, Bajardi P, Pepe E, Lake B, Privitera F, Tizzoni M. Socio-economic determinants of mobility responses during the first wave of COVID-19 in Italy: from provinces to neighbourhoods. Journal of The Royal Society Interface. 2021;18(181):20210092. pmid:34343450
  41. 41. Garnier R, Benetka JR, Kraemer J, Bansal S. Socio-economic disparities in social distancing during the COVID-19 pandemic in the United States: an observational study. Journal of medical Internet research. 2020;.
  42. 42. Mukandavire Z, Nyabadza F, Malunguza NJ, Cuadros DF, Shiri T, Musuka G. Quantifying early COVID-19 outbreak transmission in South Africa and exploring vaccine efficacy scenarios. PloS one. 2020;15(7):e0236003. pmid:32706790
  43. 43. Bubar KM, Reinholt K, Kissler SM, Lipsitch M, Cobey S, Grad YH, et al. Model-informed COVID-19 vaccine prioritization strategies by age and serostatus. Science. 2021;371(6532):916–921. pmid:33479118
  44. 44. Matrajt L, Eaton J, Leung T, Brown ER. Vaccine optimization for COVID-19: Who to vaccinate first? Science Advances. 2020;7(6):eabf1374. pmid:33536223
  45. 45. National Academies of Sciences, Engineering, and Medicine. Framework for Equitable Allocation of COVID-19 Vaccine. Gayle H, Foege W, Brown L, Kahn B, editors. Washington, DC: The National Academies Press; 2020. Available from: https://www.nap.edu/catalog/25917/framework-for-equitable-allocation-of-covid-19-vaccine.
  46. 46. Love J, Keegan LT, Angulo FJ, McLaughlin J, Shea KM, Swerdlow DL, et al. Continued need for non-pharmaceutical interventions after COVID-19 vaccination in long-term-care facilities. medRxiv. 2021;.
  47. 47. Gumel AB, Iboi EA, Ngonghala CN, Ngwa GA. Mathematical assessment of the roles of vaccination and non-pharmaceutical interventions on COVID-19 dynamics: a multigroup modeling approach. medRxiv. 2021; p. 2020–12.
  48. 48. Galanti M, Pei S, Yamana TK, Angulo FJ, Charos A, Swerdlow DL, et al. Social distancing remains key during vaccinations. Science. 2021;371(6528):473–474. pmid:33510018
  49. 49. Moore S, Hill EM, Tildesley MJ, Dyson L, Keeling MJ. Vaccination and non-pharmaceutical interventions for COVID-19: a mathematical modelling study. The Lancet Infectious Diseases. 2021; pmid:33743847
  50. 50. Moghadas SM, Vilches TN, Zhang K, Wells CR, Shoukat A, Singer BH, et al. The Impact of Vaccination on Coronavirus Disease 2019 (COVID-19) Outbreaks in the United States. Clinical Infectious Diseases. 2021;. pmid:33515252
  51. 51. Patel MD, Rosenstrom E, Ivy JS, Mayorga ME, Keskinocak P, Boyce RM, et al. Association of Simulated COVID-19 Vaccination and Nonpharmaceutical Interventions With Infections, Hospitalizations, and Mortality. JAMA Network Open. 2021;4(6):e2110782–e2110782. pmid:34061203
  52. 52. Yang J, Marziano V, Deng X, Guzzetta G, Zhang J, Trentini F, et al. Despite vaccination, China needs non-pharmaceutical interventions to prevent widespread outbreaks of COVID-19 in 2021. Nature Human Behaviour. 2021;. pmid:34158650
  53. 53. World Economic Situation and Prospects 2020; 2020. https://www.un.org/development/desa/dpad/wp-content/uploads/sites/45/WESP2020_Annex.pdf.
  54. 54. Walker PG, Whittaker C, Watson OJ, Baguelin M, Winskill P, Hamlet A, et al. The impact of COVID-19 and strategies for mitigation and suppression in low-and middle-income countries. Science. 2020;369(6502):413–422. pmid:32532802
  55. 55. Monod M, Blenkinsop A, Xi X, Hebert D, Bershan S, Tietze S, et al. Age groups that sustain resurging COVID-19 epidemics in the United States. Science. 2021;371 (6536). pmid:33531384
  56. 56. Verity R, Okell L, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet Infectious Diseases. 2020;20. pmid:32240634
  57. 57. Lopez Bernal J, Andrews N, Gower C, Robertson C, Stowe J, Tessier E, et al. Effectiveness of the Pfizer-BioNTech and Oxford-AstraZeneca vaccines on covid-19 related symptoms, hospital admissions, and mortality in older adults in England: test negative case-control study. BMJ. 2021;373. pmid:33985964
  58. 58. Coronavirus Pandemic (COVID-19); 2020. https://ourworldindata.org/coronavirus.
  59. 59. Sunnåker M, Busetto AG, Numminen E, Corander J, Foll M, Dessimoz C. Approximate Bayesian Computation. PLOS Computational Biology. 2013;9(1):1–10. pmid:23341757
  60. 60. Google LLC “Google COVID-19 Community Mobility Reports”; 2020. https://www.google.com/covid19/mobility/.
  61. 61. Oxford COVID-19 Government Response Tracker; 2020. https://www.bsg.ox.ac.uk/research/research-projects/coronavirus-government-response-tracker#data.
  62. 62. Ferguson NM, Cummings DA, Fraser C, Cajka JC, Cooley PC, Burke DS. Strategies for mitigating an influenza pandemic. Nature. 2006;442(7101):448–452. pmid:16642006
  63. 63. Davies NG, Abbott S, Barnard RC, Jarvis CI, Kucharski AJ, Munday JD, et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science. 2021;372 (6538). pmid:33658326
  64. 64. Sun K, Wang W, Gao L, Wang Y, Luo K, Ren L, et al. Transmission heterogeneities, kinetics, and controllability of SARS-CoV-2. Science. 2020; p. eabe2424.
  65. 65. Tian H, Liu Y, Li Y, Wu CH, Chen B, Kraemer MUG, et al. An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China. Science. 2020;368(6491):638–642. pmid:32234804
  66. 66. Reiner RC, Barber RM, Collins JK, Zheng P, Adolph C, Albright J, et al. Modeling COVID-19 scenarios for the United States. Nat Med. 2020;. pmid:33020651
  67. 67. Schlosser F, Maier BF, Jack O, Hinrichs D, Zachariae A, Brockmann D. COVID-19 lockdown induces disease-mitigating structural changes in mobility networks. Proc Natl Acad Sci USA. 2020; p. 202012326. pmid:33273120
  68. 68. Kraemer MUG, Yang CH, Gutierrez B, Wu CH, Klein B, Pigott DM, et al. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science. 2020;368(6490):493–497. pmid:32213647
  69. 69. Maier BF, Brockmann D. Effective containment explains subexponential growth in recent confirmed COVID-19 cases in China. Science. 2020;368(6492):742–746. pmid:32269067
  70. 70. Di Domenico L, Pullano G, Sabbatini CE, Boëlle PY, Colizza V. Impact of lockdown on COVID-19 epidemic in Île-de-France and possible exit strategies. BMC Medicine. 2020;18(1):240. pmid:32727547
  71. 71. Walker PGT, Whittaker C, Watson OJ, Baguelin M, Winskill P, Hamlet A, et al. The impact of COVID-19 and strategies for mitigation and suppression in low- and middle-income countries. Science. 2020; p. eabc0035. pmid:32532802
  72. 72. Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. The Lancet Public Health. 2020;5(5):e261–e270. pmid:32220655
  73. 73. Jarvis CI, Zandvoort KV, Gimma A, Prem K, Klepac P, et al. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Med. 2020;18(1). pmid:32375776
  74. 74. Davies NG, Kucharski AJ, Eggo RM, Gimma A, Edmunds WJ, Jombart T, et al. Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. The Lancet Public Health. 2020;5(7):e375–e385. pmid:32502389
  75. 75. Sniehotta FF, Araújo-Soares V, Brown J, Kelly MP, Michie S, West R. Complex systems and individual-level approaches to population health: a false dichotomy? The Lancet Public Health. 2017;2(9):e396–e397. pmid:29253408
  76. 76. Van Bavel JJ, Baicker K, Boggio PS, Capraro V, Cichocka A, Cikara M, et al. Using social and behavioural science to support COVID-19 pandemic response. Nature Human Behaviour. 2020; p. 1–12.
  77. 77. West R, Michie S, Rubin GJ, Amlôt R. Applying principles of behaviour change to reduce SARS-CoV-2 transmission. Nature Human Behaviour. 2020; p. 1–9. pmid:32377018
  78. 78. Prem K, Cook AR, Jit M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLOS Computational Biology. 2017;13(9):1–21. pmid:28898249
  79. 79. Hao X, Cheng S, Wu D, Wu T, Lin X, Wang C. Reconstruction of the full transmission dynamics of COVID-19 in Wuhan. Nature. 2020;584 (7821). pmid:32674112
  80. 80. Matrajt L, Eaton J, Leung T, Dimitrov D, Schiffer JT, Swan DA, et al. Optimizing vaccine allocation for COVID-19 vaccines shows the potential role of single-dose vaccination. Nature Communications. 2021;12(1):3449. pmid:34103510
  81. 81. Liang M, Gao L, Cheng C, Zhou Q, Uy JP, Heiner K, et al. Efficacy of face mask in preventing respiratory virus transmission: A systematic review and meta-analysis. Travel Medicine and Infectious Disease. 2020;36:101751. pmid:32473312
  82. 82. Mitze T, Kosfeld R, Rode J, Wälde K. Face masks considerably reduce COVID-19 cases in Germany. Proceedings of the National Academy of Sciences. 2020;117(51):32293–32301. pmid:33273115
  83. 83. Voysey M, Clemens SAC, Madhi SA, Weckx LY, Folegatti PM, Aley PK, 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. The Lancet. 2020;. pmid:33306989
  84. 84. Fenichel EP, Castillo-Chavez C, Ceddia MG, Chowell G, Parra PAG, Hickling GJ, et al. Adaptive human behavior in epidemiological models. Proceedings of the National Academy of Sciences. 2011;108(15):6306–6311. pmid:21444809
  85. 85. Virtanen P, Gommers R, Oliphant TE, Haberland M, Reddy T, Cournapeau D, et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nature Methods. 2020;17:261–272. pmid:32015543
  86. 86. Harris CR, Millman KJ, van der Walt SJ, Gommers R, Virtanen P, Cournapeau D, et al. Array programming with NumPy. Nature. 2020;585(7825):357–362. pmid:32939066
  87. 87. Lam SK, Pitrou A, Seibert S. Numba: A llvm-based python jit compiler. In: Proceedings of the Second Workshop on the LLVM Compiler Infrastructure in HPC; 2015. p. 1–6.
  88. 88. Hunter JD. Matplotlib: A 2D graphics environment. Computing in Science & Engineering. 2007;9(3):90–95.
  89. 89. Temporal variation in transmission during the COVID-19 outbreak; 2021. https://epiforecasts.io/covid/.
  90. 90. Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith H, et al. The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Annals of Internal Medicine. 2020;172(9):577–582. pmid:32150748
  91. 91. Ferretti L, Wymant C, Kendall M, Zhao L, Nurtay A, Abeler-Dörner L, et al. Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science. 2020;368 (6491). pmid:32234805
  92. 92. Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance. 2020;25(10). pmid:32183930
  93. 93. Lavezzo E, Franchin E, Ciavarella C, Cuomo-Dannenburg G, Barzon L, Del Vecchio C, et al. Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo’. Nature. 2020;584. pmid:32604404
  94. 94. Li R, Pei S, Chen B, Song Y, Zhang T, Yang W, et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science. 2020;368(6490):489–493. pmid:32179701
  95. 95. Backer JA, Klinkenberg D, Wallinga J. Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20–28 January 2020. Eurosurveillance. 2020;25(5). pmid:32046819
  96. 96. Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science. 2020;368(6493):860–868. pmid:32291278
  97. 97. Centers for Disease Control and Prevention, COVID-19 Pandemic Planning Scenarios; 2021. https://www.cdc.gov/coronavirus/2019-ncov/hcp/planning-scenarios.html#table-1.
  98. 98. Weycker D, Edelsberg J, Elizabeth Halloran M, Longini IM, Nizam A, Ciuryla V, et al. Population-wide benefits of routine vaccination of children against influenza. Vaccine. 2005;23(10):1284–1293. pmid:15652671
  99. 99. Medlock J, Galvani AP. Optimizing Influenza Vaccine Distribution. Science. 2009;325(5948):1705–1708. pmid:19696313
  100. 100. Minter A, Retkute R. Approximate Bayesian Computation for infectious disease modelling. Epidemics. 2019;29:100368. pmid:31563466
  101. 101. Liu Y, Gayle AA, Wilder-Smith A, Rocklöv J. The reproductive number of COVID-19 is higher compared to SARS coronavirus. Journal of Travel Medicine. 2020;27(2). pmid:32052846
  102. 102. United Nations, Department of Economic and Social Affairs, Population Division. World Population Prospects: The 2019 Revision; 2020. https://population.un.org/wpp/Download/Metadata/Documentation/.
  103. 103. Daily number of new reported cases of COVID-19 by country worldwide; 2020. https://www.ecdc.europa.eu/en/publications-data/download-todays-data-geographic-distribution-covid-19-cases-worldwide.
  104. 104. Gozzi N, Chinazzi M, Davis JT, Mu K, y Piontti AP, Ajelli M, et al. Estimating the spreading and dominance of SARS-CoV-2 VOC 202012/01 (lineage B. 1.1. 7) across Europe. medRxiv. 2021;.