The Impact of COVID-19 Vaccination Delay: A Modelling Study for Chicago and NYC Data
Vinicius V. L. Albani, Jennifer Loria, Eduardo Massad, Jorge P. Zubelli
11 The Impact of COVID-19 Vaccination Delay: A Modelling Study for Chicago and NYC Data
Vinicius V. L. Albani , Jennifer Loria , Eduardo Massad , and Jorge P. Zubelli February 5, 2020
Abstract Background:
By the beginning of December 2020, some vaccines against COVID-19 already presented efficacy and security, which qualify them to be used in mass vaccination campaigns. Thus, setting up strategies of vaccination became crucial to control the COVID-19 pandemic.
Methods:
We use daily COVID-19 reports from Chicago and NYC from 01-Mar2020 to 28-Nov-2020 to estimate the parameters of an SEIR-like epidemiological model that accounts for different severity levels. To achieve data adherent predictions, we let the model parameters to be time-dependent. The model is used to forecast different vaccination scenarios, where the campaign starts at different dates, from 01-Oct-2020 to 01-Apr-2021. Federal University of Santa Catarina, 88.040-900 Florianopolis, Brazil, [email protected] Instituto de Matemática Pura e Aplicada, Rio de Janeiro, Brazil, [email protected] Universidad de Costa Rica, San José, Costa Rica School of Medicine, University of São Paulo and LIM01-HCFMUSP, São Paulo, Brazil School of Applied Mathematics, Fundação Getúlio Vargas, Rio de Janeiro, Brazil, [email protected] Mathematics Department, Khalifa University, Abu Dhabi, UAE, [email protected]
To generate realistic scenarios, disease control strategies are implemented whenever the number of predicted daily hospitalizations reaches a preset threshold.
Results:
The model reproduces the empirical data with remarkable accuracy. Delaying the vaccination severely affects the mortality, hospitalization, and recovery projections. In Chicago, the disease spread was under control, reducing the mortality increment as the start of the vaccination was postponed. In NYC, the number of cases was increasing, thus, the estimated model predicted a much larger impact, despite the implementation of contention measures. The earlier the vaccination campaign begins, the larger is its potential impact in reducing the COVID-19 cases, as well as in the hospitalizations and deaths. Moreover, the rate at which cases, hospitalizations and deaths increase with the delay in the vaccination beginning strongly depends on the shape of the incidence of infection in each city.
Keywords:
Vaccination; Epidemiological Models; COVID-19; Public Health Strategies; SEIR-type models Introduction
Previous pandemics have demonstrated that, as a general rule, pharmaceutical interventions are less important than non-pharmaceutical intervention in controlling the infection, however, there is a possibility that this will not be the case with the vaccines against COVID-19 [1–3].
Some few months after the emergence of SARS-CoV-2 in China, several academic laboratories and pharmaceutical industries around the world started the development of more than 100 types of different vaccines, short-circuiting in less than one year the usual time frame of new vaccines development and testing of around ten years [3,4].
There is, therefore, an enormous variety of COVID-19 vaccines being developed. As of November 2020, there were 48 vaccines in clinical trials and 146 candidate vaccines in pre-clinical evaluation [5]. Of these, 12 vaccines were in the pipeline, of which ten were in Phase 3 of clinical trials (four have already completed this phase) and two were in Phase 2 [5]. In the US, three vaccines completed Phase 3 trials, namely, Moderna, Pfizer, and AstraZeneca, and two were still in Phase 3 [5]. In order to have a significant impact on the course of the pandemic, however, safe and effective vaccines have to emerge in less time it would take the affected populations to reach natural herd immunity [3]. Therefore, an unprecedented time-schedule to roll out any effective vaccine is urgently needed. In December 2020, the CDC proposed the Phase 1 allocation schedule of vaccination, covering an estimated 264 million people in about 25 weeks from the beginning of vaccination. Phase 1a would cover 21 million health personnel and three million nursing residents. Phase 1b would cover 87 million essential workers, 100 million persons with risky medical conditions and 53 million adults older than 65 years of age [6]. This ambitious rolling out plan, however, is way behind schedule. By 7-Jan-2021, only about five million people have been vaccinated [7].
We quantify the delay impact in vaccination deployment under different scenarios using publicly available data. This is done by implementing a novel Susceptible-Exposed-Infective-Recovered-like (SEIR) model that accounts for the different levels of disease severity, asymptomatic infection, age range, and regime changes in disease spread. The model captures well the time evolution of the outbreak leading to the forecast of realistic scenarios. It is tested with publicly available data from Chicago and NYC confirming adherence to historical data. We observe that according to the disease-spread control level, the impact of postponing a mass vaccination campaign is considerable. Reopening strategies after lockdown are also accounted for in our study. It is worth mentioning that the politicization of the vaccination in many countries, the polemic around safety and efficacy of the candidate vaccines and the anti-vaccination groups campaigning against the vaccine are all contributing to hesitancy [8] and an inevitable delay in vaccination in many places around the world (not to mention the technical hurdles to roll out billions of doses necessary to control the SARS-CoV-2 pandemic). All these issues make the estimation of the death toll caused by vaccination delay of utmost importance. Materials and Methods
This section presents the epidemiological model as well as the estimation techniques used to calibrate the model parameters from observed cases of COVID-19.
The Epidemiological SEIR-like Model
The proposed epidemiological model accounts for the distribution of the population into n age ranges. For the i th age range, the corresponding group of individuals is further classified into nine compartments, namely, susceptible (S), vaccinated (V), exposed (E), asymptomatic and infective ( I A ), mildly infective ( I M ), severely infective ( I S ), critically infective ( I C ), recovered ( R ), and deceased ( D ). All infective individuals that are symptomatic but do not need to be hospitalized are considered mildly infective. By severely infective we mean those individuals that were admitted to a regular hospital bed. Those individuals that were admitted to an intensive care unit (ICU) are considered as critically infective. Only susceptible individuals are considered vaccinated, which means that if someone is vaccinated after being exposed, then he or she will pass to the asymptomatic or mildly infective compartments. Before presenting the model, let us introduce the vector notation, i.e., 𝐒 = [𝑆 , ⋯ , 𝑆 𝑛 ] 𝑇 , where 𝑆 𝑖 (𝑖 = 1, ⋯ , 𝑛) represents the susceptible individuals in the i th age range. 𝐄 , 𝐕 , 𝐈 𝐴 , 𝐈 𝑀 , 𝐈 𝑆 , 𝐈 𝐶 , 𝐑 , and 𝐃 are defined similarly. Also consider the tensor product between two n -dimensional vectors, defined as 𝐗: 𝐘 = [𝑥 𝑦 , ⋯ , 𝑥 𝑛 𝑦 𝑛 ] 𝑇 . Thus, the movement between the model compartments is determined by the following system of ordinary differential equations:
𝐒̇ = −𝐒: (𝛽 𝐴 𝐈 𝐴 + 𝛽 𝑀 𝐈 𝑀 + 𝛽 𝑆 𝐈 𝑆 + 𝛽 𝐶 𝐈 𝐶 ) − 𝜈: 𝐒 , (1) 𝐕̇ = ν: 𝐒, (2)
𝐄̇ = 𝐒: (𝛽 𝐴 𝐈 𝐴 + 𝛽 𝑀 𝐈 𝑀 + 𝛽 𝑆 𝐈 𝑆 + 𝛽 𝐶 𝐈 𝐶 ) − σ𝐄, (3) 𝐈̇ 𝐴 = (1 − 𝑝)𝜎𝐄 − 𝛾 𝑅,𝐴 : 𝐈 𝐴 , (4) 𝐈̇ 𝑀 = 𝑝𝜎𝐄 − (𝛾 𝑅,𝑀 + 𝛼 𝑆 ): 𝐈 𝑀 , (5) 𝐈̇ 𝑆 = 𝛼 𝑆 : 𝐈 𝑀 − (𝛾 𝑅,𝑆 + 𝛼 𝐶 ): 𝐈 𝑆 , (6) 𝐈̇ 𝐶 = 𝛼 𝐶 : 𝐈 𝑆 − (𝛾 𝑅,𝐶 + 𝛿 𝐷 ): 𝐈 𝐶 , (7) 𝐑̇ = 𝛾
𝑅,𝐴 : 𝐈 𝐴 + 𝛾 𝑅,𝑀 : 𝐈 𝑀 + 𝛾 𝑅,𝑆 : 𝐈 𝑆 + 𝛾 𝑅,𝐶 : 𝐈 𝐶 , (8) 𝐃̇ = 𝛿 𝐷 : 𝐈 𝐶 . (9) The schematic representation of the model can be found in Figure 1. Figure 1: Schematic representation of the epidemiological model of Eqs. (1)-(9).
The time-dependent transmission parameters for asymptomatic, mildly, severely, and critically infective individuals are denoted, respectively, by 𝛽 𝐴 , 𝛽 𝑀 , 𝛽 𝑆 , and 𝛽 𝐶 . The rate of vaccination is 𝜈 , which is given by the product of the daily rate of vaccination of susceptible individuals by the effectiveness of the used vaccine. The meantime from contagion to become infective is the inverse of the parameter 𝜎 . The recovery rate of mildly, severely, critical and asymptomatic infective individuals are denoted by 𝛾 𝑅,𝑀 , 𝛾 𝑅,𝑆 , 𝛾 𝑅,𝐶 , and 𝛾 𝑅,𝐴 , respectively. The rates of hospitalization and ICU admission are denoted by 𝛼 𝑆 and 𝛼 𝐶 , respectively. According to the World Health Organization, only people in severe conditions generally die by COVID-19, thus, the corresponding death rate is 𝛿 𝐷 [9]. The unknown parameters are 𝛽 𝐴 , 𝛽 𝑀 , 𝛽 𝑆 , and 𝛽 𝐶 , as well as the initial number of mildly and asymptomatic infective individuals, that shall be estimated from the daily numbers of infections. In order to reduce the number of unknowns, we assume that 𝛽 𝑆 = 𝑎𝛽 𝑀 , 𝛽 𝐶 = 𝑏𝛽 𝑀 , and 𝛽 𝐴 = 𝑐𝛽 𝑀 , (10) with 𝑎 = 0.1 , 𝑏 = 0.01 , and 𝑐 = 0.58 , which means that the infection rate of hospitalized, in ICU and asymptomatic individuals are 10%, 1% and 58%, respectively, of the transmission rate of those ones in the mildly infective compartment [10,11]. The mean time between infection and becoming infective is set to 5.1 [12]. The proportion of exposed individuals becoming mildly infective is 𝑝 , which is set to 0.83 [11]. The recovery rates of mildly, severely, and critically infective individuals are simply set as one minus the rates of hospitalization, ICU admission, and death, respectively. All the asymptomatic individuals will recover in 14 days, so, 𝛾 𝑅,𝐴 = 14 −1 . The rate of ICU admission is set as 𝛼 𝐶 = 0.4 [13]. The hospitalization and death rates are time-dependent and defined as follows: 𝛼 𝑆 (𝑡) = 𝐻(𝑡)𝐼(𝑡−1) and 𝛿 𝐷 (𝑡) = 𝐷(𝑡)𝛼 𝐶 𝐻(𝑡−1) (11) where 𝐼 , 𝐻 , and 𝐷 represent the time series of seven-day moving averages of daily numbers of infections, hospitalizations and deaths, respectively. If the number of age ranges n is larger than one, the entries of the matrix 𝛽 𝑀 are defined as: [𝛽 𝑀 ] 𝑖𝑖 = 𝛽 𝑖 (𝑡)𝑏 𝑖 , [𝛽 𝑀 ] 𝑖𝑗 = 𝑝 𝑗 (𝛽 𝑖 (𝑡)𝑏 𝑖 + 𝛽 𝑗 (𝑡)𝑏 𝑗 ), 𝑖 ≠ 𝑗 (𝑖, 𝑗 = 1, ⋯ , 𝑛), where 𝛽 𝑖 (𝑡) (𝑖 = 1, . . . , 𝑛) are time-dependent scalar coefficients, and 𝑏 𝑖 , as well as 𝑝 𝑖 (𝑖 =1, . . . , 𝑛) are time-independent [10]. 𝛽 𝑖 (𝑡) represents the time-dependent part of the transmission coefficient for the i th age range, whereas 𝑏 𝑖 is the time-independent part. They, respectively, capture the short-term and the long-term pattern of the disease spread. The parameter 𝑝 𝑖 represents the probability of any individual from the j th age range to interact with an individual from the ( 𝑖 + 𝑗 – 1 )th age range. In this case, we set 𝑝 = 1 . Under these assumptions, the number of unknown parameters in 𝛽 𝑀 , for each 𝑡 , is , instead of 𝑛 − 𝑛 . Estimation Techniques
As aforementioned, the set of unknown parameters is composed by the initial number of individuals at the mildly infective individuals 𝐈 𝑀 (0) = [𝐼 𝑀1 (0), ⋯ , 𝐼 𝑀𝑛 (0)] 𝑇 , the time-dependent components of the matrix 𝛽 𝑀 , 𝛽⃗(𝑡) = [𝛽 (𝑡), ⋯ , 𝛽 𝑛 (𝑡)] 𝑇 , and the time-independent components of 𝛽 𝑀 , 𝐛 = [𝑏 , ⋯ , 𝑏 𝑛 ] 𝑇 and 𝐩 = [1, 𝑝 , ⋯ , 𝑝 𝑛 ] 𝑻 . Let ℑ 𝑖 = {ℐ 𝑖 (𝑡 ), . . . , ℐ 𝑖 (𝑡 𝑁 )} represent the time series of seven-day moving averaged version of daily reports of COVID-19 infections for the i th age range. The estimation procedure is performed in two steps. In the first one, the time-independent parameters, 𝐈 𝑀 (0) and 𝐛 , are estimated and, in the second one, 𝛽⃗(𝑡) is calibrated. The estimated 𝐈 𝑀 (0) and 𝐛 are minimizers of the objective function below, which is closely related to the so-called posterior density of the model parameters, given the set of reports [14]. Π(𝐈 𝑀 (0), 𝐛, 𝐩|ℑ , ⋯ , ℑ 𝑛 ) = 𝐿(ℑ , ⋯ , ℑ 𝑛 |𝐈 𝑀 (0), 𝐛, 𝐩) + 𝛼𝑃𝑟(𝐈 𝑀 (0), 𝐛, 𝐩) , (12) where 𝛼 > 0 is a penalty parameter, 𝐿(ℑ , ⋯ , ℑ 𝑛 |𝐈 𝑀 (0), 𝐛, 𝐩) = ∑ ∑[ℐ 𝑖 (𝑡 𝑙 )𝑙𝑜𝑔(𝜎𝐸 𝑖 (𝑡 𝑙 )) − 𝜎𝐸 𝑖 (𝑡 𝑙 ) − 𝑙𝑜𝑔(ℐ 𝑖 (𝑡 𝑙 )!)] 𝑁𝑙=1𝑛𝑖=1 is the data misfit function with 𝑙𝑜𝑔(ℐ 𝑖 (𝑡 𝑙 )!) approximated by the Stirling’s formula 𝑙𝑜𝑔(ℐ 𝑖 (𝑡 𝑙 )!) ≈ 𝑙𝑜𝑔(2𝜋ℐ 𝑖 (𝑡 𝑙 )) + ℐ 𝑖 (𝑡 𝑙 )𝑙𝑜𝑔(ℐ 𝑖 (𝑡 𝑙 )) − ℐ 𝑖 (𝑡 𝑙 ) , and 𝑃𝑟(𝐈 𝑀 (0), 𝐛, 𝐩) = ‖(𝐈 𝑀 (0), 𝐛, 𝐩) − (𝐈 𝑀0 (0), 𝐛 , 𝐩 )‖ ℓ is the penalty term, with (𝐈 𝑀0 (0), 𝐛 , 𝐩 ) a set of a priori parameters. The set of time-dependent components of 𝛽 𝑀 , 𝛽⃗(𝑡) is estimated by minimizing the function below, for each 𝑡 𝑙 (𝑙 = 1, ⋯ , 𝑁) : Π(𝛽⃗(𝑡 𝑙 )|ℐ (𝑡 𝑙 ), . . . , ℐ 𝑛 (𝑡 𝑙 ), 𝐈 𝑀 (0), 𝐛, 𝐩)= ∑[ℐ 𝑖 (𝑡 𝑙 )𝑙𝑜𝑔(𝜎𝐸 𝑖 (𝑡 𝑙 )) − 𝜎𝐸 𝑖 (𝑡 𝑙 ) − 𝑙𝑜𝑔(ℐ 𝑖 (𝑡 𝑙 )!)] 𝑛𝑖=1 + α‖𝛽⃗(𝑡 𝑙 ) − 𝛽⃗(𝑡 𝑙−1 )‖ ℓ . The minimization of the objective functions above is performed by a gradient-based technique [15]. Confidence intervals (CIs) are generated by bootstrapping, based on a set of 200 samples [16]. Results
This section presents the comparison between model predictions and reported data after calibration, as well as vaccination scenarios created with the calibrated model using data from Chicago and New York City (NYC). Estimation Results
The parameters of the epidemiological model of Eqs. (1)-(9) are estimated from seven-day moving average time-series of daily new infections from Chicago and NYC. The time-series of daily reports of COVID-19 infections, as well as related hospitalizations and deaths for Chicago and NYC, are available from public resources [17,18]. Recent census data containing the total population of the considered cities and their distributions by age ranges were also used [19]. During the model estimation, the time series were divided into sets of consecutive 20 days. Besides the set corresponding to the beginning of the COVID-19 outbreak in these cities, for each 20-day dataset, 𝐛 and 𝛽⃗(𝑡) are estimated. We start by not distributing the population into age ranges, which means that n is set to 1 in the model of Eqs. (1)-(9). The time-dependent effective reproduction number denoted by ℛ(𝑡) is evaluated through the next-generation matrix technique [20,21]. Model prediction using estimated parameters of daily new cases, hospitalizations and deaths, as well as the corresponding reported numbers for Chicago and NYC for the period 01-Mar-2020 to 28-Nov-2020 can be found in Figures 2-3, respectively. The corresponding effective reproduction numbers are also presented. Figure 2: Model predictions of infections, hospitalizations and deaths (solid lines), using data from Chicago. The bars represent the reports and the envelopes are 90% CIs. The corresponding seven-day moving average of the time-dependent basic reproduction rate is depicted in the bottom right panel.
For both cities, the model predictions of daily new cases, hospitalizations and deaths (in Figure 2, top left and right, as well as bottom left panels, respectively) are adherent to the reports. It is explained by the effectiveness of the calibration procedure, and the use of the hospitalization and death rates defined in Eq. (11). We decided to present the seven-day moving average of
ℛ(𝑡) since it is less fuzzy, allowing to see the qualitative trend of the spread dynamics, such as the effectiveness of control measures. The periods when control measures effectively reduced the number of new COVID-19 infections are illustrated by the graph of
ℛ(𝑡) , where, during such dates, its value remained below one (solid horizontal line). Vaccination Scenarios
Let us consider that a vaccination campaign is implemented in Chicago and NYC. The vaccine is 95% effective. Firstly, during the campaign, on each day, 1% of the susceptible population is immunized, until the number of susceptible individuals is less than 40%. This threshold was chosen based on an estimate of the proportion of US citizens that accept to get a vaccine against COVID-19 [22].
Figure 3: Model predictions of infections, hospitalizations and deaths (solid lines), using data from NYC. The bars represent the reports and the envelopes are 90% CIs. The corresponding seven-day moving average of the time-dependent basic reproduction rate is depicted in the bottom right panel.
In order to forecast scenarios, the time-dependent parameters are extended to the forecast period by repeating the average of the values estimated in the last ten days of the calibration period. To avoid unrealistic numbers, whenever the predicted number of daily hospitalizations reaches the value 300, the time-dependent transmission coefficient 𝛽(𝑡) is set to the average of the values estimated in the period 07-Sept-2020 to 16-Sept-2020, when the disease spread was controlled. In this period, the effective reproduction numbers in Figures 2-3 were close to the value one, indicating that the disease spread was under control. The vaccination campaign is set in the period 01-Oct-2020 to 31-May-2021, starting on different dates, but finishing at 31-May-2020. Table 1 shows the accumulated numbers of infections, hospitalizations and deaths corresponding to the different starting dates. The evolution of the number of accumulated deaths with respect to the starting date of the vaccination campaign can be found in Figure 4. The increasing number of deaths, as the beginning of the campaign is postponed, also illustrate that vaccination must start as soon as possible. Starting Date Cases Hospitalizations Deaths Total Vaccinated
Table 1: Accumulated numbers of infections, hospitalizations and deaths in Chicago, when the vaccination campaign starts at different dates. Such values are based on model predictions using the estimated parameters . Starting Date Cases Hospitalizations Deaths Total Vaccinated
Table 2: Accumulated numbers of infections, hospitalizations and deaths in NYC, when the vaccination campaign starts at different dates. Such values are based on model predictions using the estimated parameters. Figure 4: Evolution of the model predicted accumulated deaths in Chicago (left) and in NYC (right) with respect to the starting date of the vaccination campaign. The envelopes are 90% CIs.
Figure 5: Model predictions for the increment in the accumulated deaths by postponing the starting date of the vaccination campaign in Chicago (left) and in NYC (right). The envelopes are 90% CIs.
An example using a vaccination strategy that accounts for age range can be found in the supplement. The corresponding conclusions and results are similar to the ones above. Discussion
In this paper, to generate vaccination scenarios, we propose an SEIR-like model that accounts for the different levels of disease severity, asymptomatic infection, age range, and regime changes in disease spread as times goes by. The model parameters are calibrated from reports of daily COVID-19 infections, as well as published reports. We end-up with a modeling tool that captures well the time evolution of the outbreak, reproducing the empirical data with remarkable accuracy, helping to forecast realistic scenarios. Such features are illustrated using publicly available data from Chicago and NYC. Depending on, whether the disease spread is under control or not, that is, whether the daily incidence curve of infection is increasing or decreasing, the impact of postponing the beginning of a mass vaccination campaign is considerable. As expected, such impact is more serious in regions where the incidence curve is increasing than in cities where the infection is controlled. We use different strategies and consider the implementation of contention measures, as the daily reports of hospitalizations reach a threshold. Reopening strategies after lockdown are also accounted for in our study. The model has some important limitation worth mentioning. First, it assumes that 60% of susceptible are vaccinated with a 95% efficacy vaccine in a short period of time at a rate of 1% per day. Although this scenario is logistically feasible, it is a daunting task. The current scenario of the pandemic, in which new variants of SARS-CoV-2 are emerging in some countries, should be considered in the simulation of future vaccination models [23]. However, there is not enough empirical evidence of the repercussion of these new variants of the vaccine efficacy. Finally, we should point out that the SARS-CoV-2, like any other viruses, is evolving, with new strains showing increased transmissibility. It should be expected, however, that its case fatality rate (or virulence) should be decreasing with time. This is a general rule in the evolution of pathogens which helps them to increase their basic reproduction number [24]. If this is the case, then it is possible to predict that in few years, COVID-19 tends to be a mild disease as other coronaviruses, like OC-43 which probably caused the so-called “Russian flu” in 1889 and nowadays is responsible for about 10% of the common cold [25]. The future of vaccines against SARS-CoV-2, therefore, will very much depend on the virulence the virus will eventually evolve. ICMJE Statement
All authors attest they meet the ICMJE criteria for authorship.
Contributors
VA, EM and JZ proposed the mathematical model. VA and JL performed numerical simulations. VA and JL analyzed the data. All the authors contributed to the writing of the article. EM and JZ critically revised the manuscript. All the authors had full access to all data used in the study and take responsibility for the accuracy of the data analysis. All authors revised and approved the final version of the article.
Declaration of interests
All authors declare no competing interests.
Data sharing
The data that support the findings of this study are available from publicly available sources [17,18]. The numerical scripts used to generate the simulated scenarios can be found in the GitHub repository https://github.com/JennySorio/Vaccination_Scenarios . Funding
This work was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) [grant numbers 305544/2011-0 and 307873/2013-7], the Fundação Butantan [grant number 01/2020], the Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro [grant number E-26/202.927/2017], and the Universidad de Costa Rica (UCR) [grant number OAICE-CAB-02-022-2016].
Appendix A
A.1 One Additional Example
We now simulate vaccination campaigns, starting on different dates, where people over 80 years old are immunized one month earlier than individuals in other age ranges. In addition, people under 18 years old are not vaccinated. Again, we assume 95% of effectiveness of the vaccine and the rate of vaccination of the population in the i th age range is 1%. Starting Date Cases Hospitalizations Deaths Total Vaccinated
Table 3: Accumulated numbers of infections, hospitalizations and deaths in NYC, when the vaccination campaign starts at different dates. Such values are based on model predictions using the estimated parameters obtained in Section 3.1.
Table 3 presents the accumulated numbers of COVID-19 cases, hospitalizations and deaths, as well as of immunized individuals during the period 01-Nov-2020 to 31-May-2021 for NYC. As the starting date of the campaign is delayed, there is a remarkable increase in the accumulated numbers. Figure 6: Right: Evolution of the model predicted accumulated deaths in NYC with respect to the starting date of the vaccination campaign. Left: increment in the accumulated deaths by postponing the starting date of the vaccination campaign. The envelopes are 90% CIs.
The left panel in Figure 6 shows the evolution of the accumulated number of deaths as a function of the vaccination campaign starting date. The right panel shows the increment in the number of deaths for each starting date, in comparison to starting vaccination one month earlier.
References [1] McKeown T, Brown RG. Medical evidence related to english population changes in the eighteenth century. Popul Stud (NY) 1955;9:119–41. https://doi.org/10.1080/00324728.1955.10404688. [2] Mckeown T, Lowe CR. An Introduction to Social Medicine. Oxford and Edinburgh: Blackwell; 1966. [3] Christakis N. Appolo’s Arrow: The Profound and Enduring Impact of Coronavirus on the Way we Live. New York: 2020. [4] Pronker ES, Weenen TC, Commandeur H, Claassen EHJHM, Osterhaus ADME. Risk in Vaccine Research and Development Quantified. PLoS One 2013;8. https://doi.org/10.1371/journal.pone.0057755.