Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

The 2020 SARS-CoV-2 epidemic in England: key epidemiological drivers and impact of interventions

Edward S. Knock, Lilith K. Whittles, John A. Lees, Pablo N. Perez-Guzman, Robert Verity, Richard G. FitzJohn, Katy AM Gaythorpe, Natsuko Imai, Wes Hinsley, Lucy C. Okell, Alicia Rosello, Nikolas Kantas, Caroline E. Walters, Sangeeta Bhatia, Oliver J Watson, Charlie Whittaker, Lorenzo Cattarino, Adhiratha Boonyasiri, Bimandra A. Djaafara, Keith Fraser, Han Fu, Haowei Wang, Xiaoyue Xi, Christl A. Donnelly, Elita Jauneikaite, Daniel J. Laydon, Peter J White, Azra C. Ghani, Neil M. Ferguson, Anne Cori, View ORCID ProfileMarc Baguelin
doi: https://doi.org/10.1101/2021.01.11.21249564
Edward S. Knock
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lilith K. Whittles
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John A. Lees
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pablo N. Perez-Guzman
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert Verity
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard G. FitzJohn
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katy AM Gaythorpe
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Natsuko Imai
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wes Hinsley
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lucy C. Okell
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alicia Rosello
4Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nikolas Kantas
5Faculty of Natural Sciences, Department of Mathematics, Imperial College London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Caroline E. Walters
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sangeeta Bhatia
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Oliver J Watson
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Charlie Whittaker
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lorenzo Cattarino
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Adhiratha Boonyasiri
3Department of Infectious Disease, School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bimandra A. Djaafara
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keith Fraser
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Han Fu
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haowei Wang
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiaoyue Xi
5Faculty of Natural Sciences, Department of Mathematics, Imperial College London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christl A. Donnelly
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
6Department of Statistics, University of Oxford, Oxford, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elita Jauneikaite
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel J. Laydon
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter J White
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
2National Institute for Health Research Health Protection Research Unit in Modelling and Health Economics, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Azra C. Ghani
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Neil M. Ferguson
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: m.baguelin@imperial.ac.uk neil.ferguson@imperial.ac.uk
Anne Cori
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marc Baguelin
1MRC Centre for Global Infectious Disease Analysis, Abdul Latif Jameel Institute for Disease and Emergency Analytics (J-IDEA), School of Public Health, Imperial College London; UK
4Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marc Baguelin
  • For correspondence: m.baguelin@imperial.ac.uk neil.ferguson@imperial.ac.uk
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

We fitted a model of SARS-CoV-2 transmission in care homes and the community to regional surveillance data for England. Among control measures implemented, only national lockdown brought the reproduction number below 1 consistently; introduced one week earlier it could have reduced first wave deaths from 36,700 to 15,700 (95%CrI: 8,900–26,800). Improved clinical care reduced the infection fatality ratio from 1.25% (95%CrI: 1.18%–1.33%) to 0.77% (95%CrI: 0.71%–0.84%). The infection fatality ratio was higher in the elderly residing in care homes (35.9%, 95%CrI: 29.1%–43.4%) than those residing in the community (10.4%, 95%CrI: 9.1%–11.5%). England is still far from herd immunity, with regional cumulative infection incidence to 1st December 2020 between 4.8% (95%CrI: 4.4%–5.1%) and 15.4% (95%CrI: 14.9%–15.9%) of the population.

One-sentence summary We fit a mathematical model of SARS-CoV-2 transmission to surveillance data from England, to estimate transmissibility, severity, and the impact of interventions

1 Introduction

England is among the countries worst-affected by the global pandemic of COVID-19, caused by the novel Betacoronavirus SARS-CoV-2. As of 2nd December 2020, over 51,000 deaths have been reported nationally, or 91 deaths per 100,000 people (1). The impact of the epidemic has varied across the country, with regional epidemics differing in their severity and timing. A key feature in all regions is the burden suffered by older adults living in care homes, where mortality has been high.

We use a mathematical model of SARS-CoV-2 transmission to reproduce the first two waves of the epidemic across England’s seven NHS regions and assess the impact of interventions implemented by the UK government. We analyse the epidemic from importation of SARS-CoV-2 into each region to the 2nd December 2020: encompassing the first national lockdown from March –May, the interventions implemented as COVID-19 deaths increased again in the autumn, eventually leading to the second national lockdown in November.

We built an age-structured stochastic transmission model of SARS-CoV-2, representing care homes, hospital clinical pathways and the wider community (Materials and Methods). We developed a Bayesian evidence-synthesis approach to estimate model parameters and to reconstruct regional epidemics using data from daily recorded deaths, PCR testing, hospital admissions, hospital bed occupancy, individual patient outcomes, contact surveys, and serological surveys. We evaluated temporal changes in transmission as new control measures were implemented and then relaxed, and population immunity accrued. Inclusion of serological data allowed us to robustly estimate region- and age-specific disease severity, to compare severity in care home residents to elderly individuals in the community, and estimate the total epidemic size, by calculating the proportion of individuals infected over time in each region. Finally, we examined counterfactual epidemic scenarios, varying the date and duration of the first national lockdown and the effectiveness of restricting care home visits, to quantify the resulting impact on mortality.

Our analysis, which synthesises multiple data sources and parametrically accounts for their biases, provides a comprehensive overview of transmission, hospitalisation, and mortality patterns of SARS-CoV-2 in the first and second waves (up to 2nd December) in all regions of England. Our results provide crucial insights for controlling the epidemic in the future, emphasising the importance of acting fast to save lives.

2 Results

2.1 Epidemic trajectory

We used our evidence-synthesis approach, to infer the COVID-19 epidemic start date in each NHS England region (Figure 1A), then reconstructed epidemic trajectories for hospitalisations (Figure S7) and deaths in care homes and hospitals (Figure 1B-H). We estimated the basic reproduction number, R0, defined as the expected number of onward infections from an infectious individual in a fully susceptible population to be 2.9 (95% CrI: 2.8-3.1) nationally. Figure 1I shows how the effective reproduction number Embedded Image (the expected number of onward infections from an individual infected at time t) changed in each region over time, in relation to government control measures and accrual of population immunity.

Figure 1:
  • Download figure
  • Open in new tab
Figure 1:

Trajectory of the England COVID-19 epidemic. A, The inferred epidemic start date in each NHS England region. B-H, The model fit to reported daily deaths from COVID-19 in care homes and hospitals for each NHS England region. The points show the daily data, solid lines the median posterior and the shaded area shows the 95% CrI. I, The mean effective reproduction number within the general community (i.e. excluding care homes) in each region from March to December. Vertical lines and labels represent dates of key policy changes, defining the breaking points of the underlying piecewise linear transmission rate. Dashed horizontal line depicts reproduction number of 1.

The first COVID-19 death in England occurred on 5th March 2020 (2). Seven days later, in response to the growing epidemic, the government began to introduce control measures, initially requiring individuals with a dry persistent cough and/or fever to self-isolate (3). On 23rd March this escalated to a full national lockdown (4, 5). Irrespective of initial differences, the level of transmission during lockdown was similar across all regions (Figure 1I), consistent with mobility data showing movement during lockdown reduced to a consistent level nationally (6).

The epidemic in London began 15 days before (95% CrI: 28 days before, 3 days after) the rest of the country (Figure 1A), meaning the lockdown occurred at a later stage of its epidemic. London experienced a mortality of 88.5 (95% CrI: 79.9–95.3) per 100,000 during the first wave, compared to the national average of 70.7 (95% CrI: 64.6–77.1), despite having a younger population and a smaller care home population than other regions (296 vs 603 per 100,000 nationally).

A key feature of the first epidemic wave in England, in common with other European countries, was the high death toll within care homes, which accounted for 22.6% laboratory-confirmed COVID-19 deaths in England as of 1st August 2020. Although community transmission rates fell during lockdown, transmission within care homes continued to rise, with infection risk peaking in care home residents, between 26th March in London and 12th April in North East and Yorkshire (Figure 2A). Deaths in care homes peaked on average 13 days later than hospital deaths (Figure 1B-H).

Figure 2:
  • Download figure
  • Open in new tab
Figure 2:

A, Inferred daily SARS-CoV-2 infections in England care home residents (right axis) and the wider community (left axis). B-H, Comparison of modelled (shaded bands) and observed (solid line) proportion of PCR tests that are positive, under pillar-2 testing (community swab testing for symptomatic individuals) in >25 year olds. Shaded bands depict 95% CrI, 50% CrI and median model outputs.

The first lockdown in England continued until 11th May, when people unable to work remotely were permitted to resume their jobs. Over the summer restrictions were successively eased, with non-essential shops, pubs and restaurants opening, followed by the government’s ‘Eat Out to Help Out’ restaurant subsidy scheme in August (7). This led to a steady increase in transmission, with Embedded Image rising above 1 in all regions by mid-August (Figure 1I).

Increasing PCR test positivity marked the beginning of a second epidemic wave (Figure 2B-H, S6). The accompanying introduction of non-pharmaceutical interventions (NPIs) began with the “Rule of Six” (limiting social gatherings to 6 persons maximum) on 14th September (8), followed by the localised tiered restrictions on 14th October (9). These measures limited transmission in most regions but were not sufficient to reduce Embedded Image below 1 (Figure 1I). Consequently, on 31st October, the government announced a second national lockdown, which lasted from 5th November to 1st December (10).

Restrictions during the second lockdown were less stringent than the first, with schools and some workplaces remaining open. This was reflected in Embedded Image estimates of 0.83 (95% CrI: 0.81–0.85) at the start of the second lockdown, compared to Embedded Image (95% CrI: 0.50–0.59) at the start of the first. We estimate that without the population immunity accrued during the first wave, contact rates during the second lockdown would have resulted in a reproduction number of Rt = 0.95 (95% CrI: 0.93–0.98). Hence, population immunity helped to reduce transmission further below the critical threshold of Embedded Image.

2.2 Severity and hospitalisation

COVID-19 manifests a broad spectrum of severity, from asymptomatic infection to life-threatening illness requiring intensive care. We estimated age-patterns of clinical progression in people admitted to hospital using individual-level data from 17,702 patients admitted between 18th March and 31st May 2020 (inclusive) in the COVID-19 Hospitalisation in England Surveillance System (CHESS, (11)) (Materials and Methods). We derived estimates of the time spent in each stage of the hospital pathway (including general wards, ICU and post-ICU stepdown care), as well as age-stratified probabilities of progression through that pathway (Figure 3 and Figure S8). We accounted for differing length of stays given different outcomes; there were marked differences in average length of ICU stay for those who died in ICU, those who later died in stepdown care and those who were discharged following stepdown care (Figure 3F). Among patients over 65, we found the probability of admission to ICU decreased with increasing age. Severity of COVID-19 increases with age, but for older patients and those with most severe illness, the benefit of ICU admission, ventilation and the corresponding prognosis may not be better than with oxygen therapy in a general ward (12). Thus, older and more severely infected patients may be directed to care on a general ward rather than admitted to ICU.

Figure 3:
  • Download figure
  • Open in new tab
Figure 3:

Age-dependent probabilities of progression through hospital pathways. A, Probability of admission to ICU. B, Probability of death in a general ward. C, Probability of death in ICU. D, Probability of death in stepdown care. E, Probability of death through all hospital pathways. Black circles and vertical segments show posterior mean and 95% credible intervals of splines fitted to data, blue circles and vertical segments show raw mean values and 95% confidence intervals (exact binomial) for each 5-year age group. F, Average length of stay in each ward (posterior mean and 95% credible intervals).

We used estimates of clinical progression to parametrise the transmission model, enabling us to infer temporal and regional differences in disease severity, informed by local demography, observed daily hospital admissions, bed occupancy and deaths. We measured severity of disease by the infection fatality ratio (IFR) and the infection hospitalisation ratio (IHR).

The severity of disease increased with age in all regions with the steepest increase above 65 years (Figure 4A-C), in line with observations worldwide (5). Regional estimates of age-aggregated disease severity depend on the population age distribution, which is similar in most regions of the country, except London, where the median age is 34.6 years (vs 39.5 years nationally). At the start of the first wave, London experienced an IFR (respectively IHR) of 0.91% (95% CrI: 0.82%–1.00%) (resp. 3.02%; 95% CrI: 2.82%–3.19%) compared to the national average of 1.25% (95% CrI: 1.18%–1.33%) (resp. 3.52%; 95% CrI: 3.29%–3.72%) (Figure 4D-E).

Figure 4:
  • Download figure
  • Open in new tab
Figure 4:

Relative severity of disease by age group and region. A, and B, Variation in the Infection fatality ratio (IFR) and Infection Hospitalisation Ratio (IHR) by age group in each region. Ages 80+ were modelled as a single risk group, care home residents are not included. C, The England IFR and IHR by age group and in care home residents (estimates denoted CHR at the right-hand side of the panel). National severity estimates are produced by aggregating regional estimates based on infection incidence. D, The regional IHR, aggregated over age and risk group by infection incidence. Plots a-d use parameter estimates, and incidence weightings calculated as of 1st December 2020. E, The England IFR over time, coloured dots show regional estimates of IFR at the start of the epidemic and on 1st December 2020 (clusters each correspond to one time-point, LON: London). In plots C-E Shaded bands depict 95% CrI and interquartile ranges, points depict medians.

Regional variation in the population age distribution did not fully account for differences in severity, with London still experiencing lower mortality when stratified by age (Figure 4A-B). The oldest age group (80+) in London had an IFR of 6.1% (95% CrI: 5.2%–6.8%) compared to 12.7% (95% CrI: 10.8%–14.3%) in the North West.

We estimated temporal trends in the IFR for England, by weighting regional estimates by incidence and population demographics. At the start of the first wave, the national IFR was 1.25% (95% CrI: 1.18%–1.33%) (Figure 4E), consistent with earlier reports from serology data alone (13). The national IFR initially appeared to increase, as transmission widened from London to regions with older populations and greater disease severity. Over the first wave, the proportion of hospital admissions resulting in death decreased, due to improvements in clinical management and alleviation of capacity constraints (14), leading to a national IFR of 0.77% (95% CrI: 0.71%–0.84%) by the end of the first wave. The magnitude of the relative reduction in IFR over time varied between regions, from 36.5% (95% CrI: 26.5%–47.5%) in the North West to 64.6% (95% CrI: 58.6%–68.8%) in London.

The IFR was greater among care home residents (35.9%, 95% CrI: 29.1%–43.4%) than in the 80+ in the community (10.4%, 95% CrI: 9.1%–11.5%, Figure 4C). Many care home residents did not transfer into hospital, and instead died in the facilities where they lived, so conversely the IHR was lower in care home residents (19.1%, 95% CrI: 11.5%–26.8%) than in those aged 80+ in the community (51.1%, 95% CrI: 47.6%–54.3%). We present national estimates of severity at the end of the second wave, stratified by age and care home residency in Table S9.

2.3 Epidemic size

Data from repeated serological surveys of blood donors aged 15-65 informed our estimation of the total regional epidemic size (Figure 5A-G), accounting for imperfect sensitivity and specificity of serological tests (Materials and Methods) (15). The cumulative proportion of the population ever infected with SARS-CoV-2 ranged from 4.8% (95% CrI: 4.4%–5.1%) in the South West to 15.4% (95% CrI: 14.9%–15.9%) in London (Figure 5H). Predicted seropositivity was initially greater than cumulative incidence, due to imperfect test specificity. The increase in seropositivity lagged cumulative infections by two weeks, reflecting the time from infection to seroconversion.

Figure 5:
  • Download figure
  • Open in new tab
Figure 5:

Cumulative incidence and seropositivity by region. A-G, Comparison of the estimated proportion of the population testing seropositive with observations from serological surveys. Vertical grey shaded bands show serological survey timings, black points the observed seroprevalence (bars: 95% exact confidence intervals), blue and purple lines show estimated proportion of the population infected and seropositive respectively (shaded bands the 95% CrI, 50% CrI and median). H, Comparison by region of the estimated cumulative attack rate in care home residents vs in the 80+ age group in the community (median, 95% CrI). The final epidemic size in each England NHS region I) in total and J) in care home residents.

Seropositivity notably declined following the first wave in some regions (Figure 5A-G). This may reflect antibody waning (16), or temporal trends in the composition of the surveyed population. Lockdown restrictions made attending blood donation centres difficult for all except key workers, who were more likely to have been infected (17), and may therefore be overrepresented in the sample of blood donors during the two lockdowns.

The proportion of care home residents ever infected with SARS-CoV-2 was 13.7% (95% CrI: 10.7%–16.7%), much higher than the 4.2% (95% CrI: 4.0%–4.4%) estimated in >80-year olds living in the community. This difference was consistently observed across all regions (Figure 5H). Regional differences in care home attack rates mirrored the patterns seen in the general community, with regions with larger community epidemics also experiencing larger care home epidemics (Figure 5I,J).

2.4 Impact of non-pharmaceutical interventions (NPIs)

We explored counterfactual intervention scenarios and examined the potential impact on mortality of initiating the first national lockdown one week earlier or later; ending that lockdown two weeks earlier or later; and 50% more or less restricted care home visits throughout the epidemic (Figure 6).

Figure 6:
  • Download figure
  • Open in new tab
Figure 6:

Counterfactual analysis of the impact on mortality aggregated across NHS England regions of A, D, initiating lockdown one week earlier / later, B, E Relaxing lockdown two weeks earlier / later, and C, F 50% more / less restricted care home visits. Panels A, B, D and E all present counterfactual outcomes for daily deaths in England but have different y-axis scales to better highlight differences between the observed data and each alternative lockdown scenario.

The timing of the initial national lockdown was crucial in determining the eventual epidemic size in England. Locking down a week earlier could have reduced the first wave death toll (up to 1st July 2020) from 36,700 to 15,700 (95% CrI: 8,900–26,800) while delaying lockdown by a week would have increased the deaths to 102,600 (95% CrI: 66,400–154,800) (Figure 6A, D). The impact varied by region, with regions with less established epidemics at the time of the first lockdown more sensitive to the timing of the intervention (Figure S10). Locking down a week later may have increased deaths, with large variability by region, from 105% in London to 274% in the Midlands but with very large uncertainty (Figure S9). Initiating a lockdown to interrupt the exponential growth phase of an epidemic has a much greater impact on reducing total mortality than extending an existing lockdown. Due to this asymmetry, relaxing the lockdown measures two weeks earlier (respectively later), could have increased deaths by 9,300 (95% CrI: 700–17,000) (respectively prevented 9,800 (95% CrI: 7,400–12,100) deaths) prior to 2nd December (Figure 6B, D).

We also explored counterfactual scenarios varying the level of visit restriction in care homes and estimated that reducing contact between the general population and care home residents by 50% could have reduced care home deaths by 44% (95% CrI: 17%–64%) (Figure 6C).

3 Discussion

We present a comprehensive overview of SARS-CoV-2 transmission, hospitalisation, mortality and intervention impact in the first two epidemic waves across all regions of England between March and December 2020. We successfully reproduce the transmission dynamics of the two epidemic waves, in terms of cases, PCR prevalence, seroprevalence, hospitalised cases (general wards and ICU), and deaths in hospitals and in care homes. We estimate intense transmission in care homes even during the first national lockdown when Embedded Image in the community was well below one in all regions (Figure 2) (18–20). Combined with our counterfactual analysis of restricting visits (Figure 6) this suggests that reducing infection levels in care home residents is challenging. This highlights the difficulty of protecting care home residents from COVID-19: due to the necessarily close contact between staff and residents within a care home, once a care home outbreak has begun it is very difficult to reduce transmission, which overrides any impact of reducing the number of introductions (21, 22).

We find that, consistent with existing literature (23), disease severity increases with age. Assessment of severity is complicated by the broad clinical spectrum of COVID-19 (24–26) hence, recent published estimates are still based on data from early in the pandemic (27). Here we provide updated severity estimates based on multiple contemporary data streams. We estimate considerable regional heterogeneity in severity, broadly consistent in the general population and in care homes for IFR and IHR. London experienced the lowest severity even after adjusting for its younger population. The estimated two-fold reduction over time in IFR (Figure 4) cannot be explained solely by the introduction of dexamethasone which reduces mortality amongst ICU patients (28), but rather a combination of factors including improvements in clinical management, greater experience in treating patients in ICU, and alleviation of capacity constraints (14, 29).

Our analysis shows large regional variation in burden, especially in the first wave. This is likely due to the pattern of seeding and the timing of lockdown relative to how advanced each region’s epidemic was (Figure 1A). Our counterfactual scenarios of initiating the first national lockdown one week earlier or later highlight the importance of early interventions to reduce overall mortality (Figure 6).

The extent and duration of infection-induced immunity to SARS-CoV-2 and its relationship to seropositivity remains unclear. Related seasonal coronaviruses induce immunity that wanes in one or two years (30), though antibody titres following SARS-CoV-1 infection appear to decay more slowly (31). Our estimated cumulative incidence over time (Figure 5), strongly supports the hypothesis that the epidemic decline after the first national lockdown was due to NPIs, with immunity playing a minimal role (32). Population-level immunity was insufficient to prevent a second wave of infection in any region (Figure 1), illustrated by the increase in reported cases and deaths which prompted the second national lockdown (33).

With the authorisation of the first SARS-CoV-2 vaccines in December 2020, we are now entering a new phase in the control of the COVID-19 pandemic. However, our estimates of current population immunity are low, with regional cumulative attack rates ranging from 4.8% to 15.4%, therefore any vaccination campaign will need to achieve high coverage and high levels of protection in vaccinated individuals to allow NPIs to be lifted without a resurgence of transmission. While vaccinating the most vulnerable age and risk groups will considerably reduce the burden of COVID-19, a large proportion of younger age groups may also need to be vaccinated to reach the immunity threshold for control. Our high estimates of transmission in care homes imply that vaccine uptake there will need to be especially high, particularly if vaccine efficacy is lower amongst older age groups.

We make a number of simplifying assumptions in our analysis. First, due to the compartmental nature of our model, we do not explicitly model individual care homes, rather the regional care home sector as a whole. However, as care home workers may work across multiple facilities leading to within and between care home transmission, we do not expect the simplification to substantially affect our conclusions. Similarly, we do not model individual households or transmission within and between them. When assessing the impact of NPIs on transmission we therefore capture population averages, rather than the contribution of household and non-household contacts. Second, hospital-acquired infections may have contributed to overall transmission, especially around the peak of the epidemic, and to persistence of infection in England over the summer months (34, 35). Our model does not explicitly represent nosocomial transmission; therefore such effects will be encompassed within our regional Embedded Image estimates. Third, each data stream was subject to competing biases, which we statistically accounted for as far as possible (supplement section 1.1.2). A key strength of our evidence-synthesis approach is that we do not rely on any single data source, combining multiple perspectives to provide a robust overall picture of the epidemic. Finally, we model the epidemics in each NHS region in England independently without accounting for spatial effects across regional boundaries, or spatial heterogeneity within regions. This spatial scale was determined by the data and reflects limited movement between regions due to travel restrictions but does allow for movement within regions.

Our analysis provides a comprehensive overview of transmission, hospitalisation, and mortality patterns of COVID-19 in the first and second waves of the epidemic in all regions of England, one of the countries worst-affected by the pandemic. Integration of multiple data streams into a single cohesive modelling framework, enables us to disentangle transmission and severity from features of the surveillance system and provide robust estimates of the epidemiological characteristics of the COVID-19 epidemic in England. As nationwide vaccination programmes are rolled out, our results will help to inform how NPIs are applied in the future.

Data Availability

All code and de-identified regionally aggregated data (see supplementary materials for full details) required to reproduce this analysis are available at https://github.com/mrc-ide/sarscov2-transmission-england (https://zenodo.org/record/4384864)

https://github.com/mrc-ide/sarscov2-transmission-england

https://zenodo.org/record/4384864

Disclaimer

The views expressed are those of the authors and not necessarily those of the United Kingdom (UK) Department of Health and Social Care, the National Health Service, the National Institute for Health Research (NIHR), Public Health England (PHE), UK MRC, UKRI or European Union.

List of Supplementary Materials

Supplementary materials (Materials and Methods, Supplementary Results) Supplementary data files: data_rtm.csv, data_serology.csv, support_progression.csv, support_severity.csv

Data availability statement

All code and de-identified regionally aggregated data (see supplementary materials for full details) required to reproduce this analysis are available at https://github.com/mrc-ide/sarscov2-transmission-england (https://zenodo.org/record/4384864)

Supplementary Materials

1 Materials and Methods

Understanding the transmission of SARS-CoV-2 is challenging. The available data are subject to competing biases, such as dependence on case definition for testing and reporting, as well as being influenced by capacity and logistical constraints. These factors are further complicated by the nature of SARS-CoV-2 transmission, whereby a substantial proportion of infected individuals develop very mild symptoms, or remain asymptomatic, but are nonetheless able to infect others (1). In this section, we describe the data used in our analyses, give details on the dynamic transmission model, and present the methods used for fitting the model to the various data sources, accounting for the inherent biases in those data.

1.1 Data sources

Here we detail the datasets used to calibrate the model to the regional epidemics. We fitted our model to time series data spanning 16th March 2020 to 2nd December 2020 (inclusive), using the data available to us on 14th December 2020, by which point the effect of remaining reporting lags would be minimal.

1.1.1 Hospital admissions and bed occupancy

We use healthcare data for each NHS region from the UK Government Dashboard (supplementary data files: data_rtm.csv, columns: phe_admissions, phe_occupied, phe_patients) (2).

For admissions data, we use the daily number of confirmed COVID-19 patients admitted to hospital, which includes people admitted to hospital who tested positive for COVID-19 in the 14 days prior to admission and inpatients who tested positive in hospital after admission, with the latter being reported as admitted on the day prior to their diagnosis.

For ICU bed occupancy, we use the daily number of (confirmed) COVID-19 patients in beds which are capable of delivering mechanical ventilation.

For the occupancy in general (i.e. non-ICU) hospital beds, we use the daily number of confirmed COVID-19 patients in hospital beds with ICU occupancy subtracted.

1.1.2 Deaths

We use the number of deaths by date of death for people who had a positive COVID-19 test result and died within 28 days of their first positive test provided Public Health England.

These can be found on (2). We also use the number among these deaths occurring in hospital (as reported by NHS England) and consider the remainder to have occurred in care homes. While non-hospital deaths may include deaths in other settings, such as in private residences, comparison with ONS data suggests that care home deaths from COVID-19 may also have been under-reported. As such we consider non-hospital deaths to be an appropriate proxy for care home deaths, and do not expect the margin for under or over-ascertainment to affect our conclusions. These data were provided by PHE and the data we have been using is provided as a supplementary file (supplementary data file: data_rtm.csv, columns: death2, death3) to allow reproducibility of our analysis.

1.1.3 Pillar 2 testing

We use pillar 2 testing data (see supplementary data files), which covers PCR testing for the general population (as compared with pillar 1 testing, which mainly occurred in hospitals). Since such testing was not available to the whole population for much of the spring wave of the pandemic, we only use this data from June 1st onwards.

We use the daily number of positives and negative tests by specimen date. Each individual who tested positive was only counted once in the number of positives, on the specimen date of their first positive test. Multiple negatives were allowed per individual, but the negatives of all individuals who ever tested positive had been removed. We only consider PCR tests and thus exclude lateral flow tests, which have been introduced recently in trials of population mass testing. We also only use pillar 2 data for those aged 25 or over, to avoid bias resulting from increased testing of university students around the reopening of universities (supplementary data file: data_rtm.csv, columns: pillar2_negatives_non_lft_over25, pillar2_positives_over25).

1.1.4 Serology surveys

Serological survey data come from antibody testing by Public Health England of samples from healthy adult blood donors, supplied by NHS Blood and Transplant (NHSBT) (supplementary data file: data_serology.csv).

1.1.5 REACT-1 prevalence survey

We use the daily number of positives and negatives by specimen date from the first 7 rounds of the REACT-1 (Real-time Assessment of Community Transmission) infection prevalence survey (supplementary data file: data_rtm.csv, columns: react_positive, react_samples) (3). Note that results published in REACT preprints use data aggregated using the administrative regions of England, whereas for the purposes of this study the data has been aggregated using NHS regions. Additionally, small changes can occur in the aggregated datasets that were published in real time because of participant withdrawals and additional data cleaning.

1.1.6 Summary of the data used for calibration

Table S 1 details the datasets used to calibrate the model to the regional epidemics.

View this table:
  • View inline
  • View popup
Table S 1:

Data sources and definitions.

1.1.7 Other data sources
1.1.7.1 Patient progression in hospital

The COVID-19 Hospitalisation in England Surveillance System (CHESS) data consists of a line list of daily individual patient-level data on COVID-19 infection in persons requiring hospitalisation, including demographic and clinical information on severity and outcomes. We use the individual dates of progression through hospital wards, from admission to eventual death or discharge, to produce age-stratified estimates of hospital progression parameters to be passed to the wider transmission model (see Section 1.9.2 and (supplementary data file: support_progression.csv, support_severity.csv).

1.1.7.2 Demographic data

We use data from the Office for National Statistics (ONS (6)) to get the number of individuals in each of the 17 age-groups, i.e. 16 five-year age bands (0-4, 5-9, …, 75-79) and an 80+ group. We get the number of care-home beds in England from (7) giving us the number of care-home beds for each NHS regions. We then got an estimate of the total population of care-home residents in the UK from (8) that we scaled down to the England population size, combined with the estimate of the total number of beds in England, we derived a value of the total occupancy of care-homes of 74.2%. We assumed that the occupancy is the same in all the NHS regions. Care-home residents are subtracted from the 4 oldest age group (5% from age 65-69, 5% age 70-74, 15% age 75-79 and 75% age 80+ (9)). We then assume a 1:1 ratio of care-home residents to care-home workers and assume that the care-home workers population is homogeneously distributed among the 25-65 population in the region.

The contact matrix between the 17 age-groups is based on the POLYMOD contact survey. See parameterisation for more details (10).

1.2 Evidence synthesis

Figure S 1 shows the functional relationships between data sources, modelled outputs and parameters in our study.

Figure S 1:
  • Download figure
  • Open in new tab
Figure S 1:

Graph showing the functional relationships between data sources (rectangles), modelled outputs (ovals) and parameters (hexagons).

1.3 Model description

We developed a stochastic compartmental transmission-dynamic model incorporating hospital care pathways to reconstruct the course of the COVID-19 epidemic in the seven NHS regions of England (Figure S 2). All analyses were done by regions, and then aggregated somehow if needed (e.g. for national IFR, or cumulative incidence). In the following description we do not mention any index denoting the region and thus all notations refer to the same NHS region.

1.3.1 Stratification of population into groups

We divided each regional population into 19 strata, denoted by the superscript i, 17 strata representing age groups within the general population, and two separate risk groups comprising care home workers (CHW) and care home residents (CHR). The 17 age groups consisted of 16 five-year age bands (0-4, 5-9, …, 75-79) and an 80+ group. The total size of the care home worker and resident groups were calculated assuming that 74.2% of available care home beds are occupied and there is a 1:1 carer to resident ratio (11). The care home workers were then split equally between all 8 age categories in the range 25 –64-year-old and removed from the corresponding age categories in the general population. Despite the care-home workers being removed from all age categories in the range 25 –64-year-old, they care-home workers are assumed to constitute one single group in our model for simplicity. The care home residents were drawn from the 65+ year old general population, such that 5% were aged 65-69, 5% aged 70-74, 15% aged 75-79 and 75% aged 80+ (9) and similarly removed from the corresponding age groups in the general population. Again, similarly to care-home workers they do constitute a single group in our model. We thus do not capture specific transmission dynamics within each care home, but rather an average mixing between residents and workers in the regional care home sector as a whole.

1.3.2 Progression of infection and hospitalisation

Prior to the importation of COVID-19, all individuals were assumed equally susceptible to infection (S). Upon infection, individuals pass through a latent period (E) before becoming infectious. A proportion (pC) of infectious individuals develop symptoms (IC) while the rest remain asymptomatic (IA). All asymptomatic individuals are assumed to recover naturally. Those with symptoms may also recover naturally (R), however a proportion (Embedded Image, age/care home-dependant as indicated by the i superscript) develop severe disease requiring hospitalisation. Of these, a proportion Embedded Image die at home without receiving hospital care. In practice this proportion is set to zero except among care home residents. Of the patients who are admitted to hospital, a proportion (p∗(t)) have their COVID-19 diagnoses confirmed prior to admission, while the remainder may be diagnosed during their inpatient stay. All hospital compartments are divided between suspected (but not yet confirmed) and confirmed diagnoses (indicated by superscript ∗). A proportion Embedded Image of new hospital admissions are triaged (ICUpre) before admission to the intensive care unit (ICU), where a fraction Embedded Image die; those who do not die get out of ICU to a ward (W) where a proportion Embedded Image die, while the remainder recover, following an inpatient care stepdown period. Inpatients not triaged to the ICU are assigned to general hospital beds (H), where a proportion Embedded Image die, while the remainder recover. Recovered individuals are assumed to be immune to reinfection for at least the duration of the simulation.

In addition, there are two parallel flows which we use for fitting to testing data: (i) for PCR positivity and (ii) for seropositivity. Upon infection, an individual enters the PCR flow in a pre-positivity compartment Embedded Image before moving into the PCR positivity compartment Embedded Image and then ultimately into the PCR negativity compartment Embedded Image. Meanwhile, individuals move into the seropositivity flow upon becoming infectious, entering first into a pre-seropositivity compartment Embedded Image. A proportion of individuals Embedded Image then seroconvert and move into the seropositivity compartment Embedded Image, while the remainder move into the seronegativity compartment Embedded Image.

We calibrated the duration distributions for each hospital compartment, and the age-stratified probabilities of moving between compartments, using the analysis of individual-level patient data (presented below in Section 1.9.2). The required Erlang distributional form was achieved within the constraints of the modelling framework by splitting each model compartment into k sequential sub-compartments (Table S 2).

Figure S 2:
  • Download figure
  • Open in new tab
Figure S 2:

Model structure flow diagram with rates of transition between infection states. Variable names defined in text.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 2:

Description of model compartments and distribution of time spent in each. For each named compartment, we give the associated duration. Due to the Markovian structure these are model Erlang-like distributions with kj the number of exponential-like compartments and γj the rate of the exponential-like compartment. 𝔼 [δj] gives the mean duration in days spent in the corresponding compartment. The structure and duration of each stage was assumed to be the same for unconfirmed and confirmed cases in hospital (see Figure S2). For length of stays related to hospital pathways, more detail is given in section 1.9.2.

1.3.3 Progression of infection and hospitalisation

The force of infection, λi(t), for individuals in group i ∈ {[0,5), …, [75,80), [80 +), CHW, CHR} depends on time-varying social mixing between age groups and prevalence in all age/care home groups: Embedded Image where mi,j (t) is the (symmetric) time-varying person-to-person transmission rate from group j to group i, and Πj (t) is the number of infectious individuals in group j, given by: Embedded Image Broadly, to parameterise mi,j (t), we informed mixing in the general population, and between the general population and care home workers using POLYMOD (10) via the R package socialmixr using age-structured regional demography (18).

Transmission between different age groups (i,j) ∈ {[0,5), …, [75,80), [80 +)}2 was parameterised as follows: Embedded Image Here ci,j is the (symmetric) person-to-person contact rate between age group i and j, derived from pre-pandemic data (10). β (t) is the time-varying transmission rate Dhich encompasses both changes over time in transmission efficiency (e.g. due to temperature) and temporal changes in the overall level of contacts in the population (due to changes in policy and behaviours).

We assumed β (t) to be piecewise linear: Embedded Image with 12 change points t i corresponding to major announcements and changes in COVID-19 related policy, as detailed in Table S 3.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 3:

Changepoints for β (t)

The contact matrix ci,j between different age groups (i,j) ∈ {[0,5), …, [75,80), [80 +)}2 is derived from the POLYMOD survey (10) for the United Kingdom using the socialmixr package (18,30), scaling by the local population demography to yield the required person-to-person daily contact rate matrix.

We defined parameters representing transmission rates within care homes (between and among workers and residents), which were assumed to be constant over time. Parameter mCHW represents the person-to-person transmission rate among care home workers and between care home workers and residents; mCHR represents the person-to-person transmission rate among care home residents. Hence, Embedded Image Embedded Image Transmission between the general population and care home workers was assumed to be similar to that within the general population, accounting for the average age of care home workers, with, for i ∈ {[0,5), …, [75,80), [80 +)}, Embedded Image where ci,CHW is the mean of ci,[25,30), ci,[30,35), …, ci,[60,65) (i.e. of the age groups that the care home workers are drawn from).

Transmission between the general population and care home residents was assumed to be similar to that between the general population and the 80+ age group, adjusted by a reduction factor (∈, which was estimated), such that, for i ∈ {[0,5), …, [75,80), [80 +)}, Embedded Image These represent contact between visitors from the general community and care home residents. This might involve a slightly different age profile than the age profile of the contact made by people in the 80+ age group.

1.3.4 Age-varying and time-varying infection progression probabilities

Various probabilities of clinical progression within the model are assumed to vary across age groups to account for severity of infection varying with age, and some are assumed to vary in time in order to model improvements in clinical outcomes, such as those achieved through the use of dexamethasone (31).

Two probabilities are age-varying but not time-varying, the probability of admission to hospital for symptomatic cases, and the probability of death for severe symptomatic cases in care homes. These were modelled as follows: Embedded Image Embedded Image where for probability Embedded Image is the maximum across all groups and Embedded Image is the age scaling such that Embedded Image for the group corresponding to the maximum, against which all other groups are scaled.

As well as varying with age, four probabilities also vary with time: the probability of admission to ICU for hospitalised cases, the probability of death in ICU, the probability of death for hospitalised cases not admitted to ICU, and the probability of death in hospital after discharge from ICU: Embedded Image Embedded Image Embedded Image Embedded Image where here for probability Embedded Image gives the maximum initial value across groups and h(μ, t) = 1 before April 1st, h(μ,, t) = μ < 1 after June 1st, with a linear reduction in between.

Care home residents with severe disease leading to death are assumed to remain in compartment GD for 5 days on average before dying (modelled with Embedded Image and Embedded Image), 95% range 0.6-13.9 days broadly consistent with durations in (15) and with duration about half the length observed in hospital streams (see Figure S 5).

For care home workers, the age scaling Embedded Image is taken as the mean of the age scalings Embedded Image for i ∈ {[25,30), [30,35), …, [60,65)}. For care home residents, we assume that Embedded Image, with the exception of the probability of individual with severe disease requiring hospitalisation dying at home (without receiving hospital care), where we assume Embedded Image and Embedded Image for all other groups, to effectively allow death outside hospital only for care home residents.

Reproduction number YW and effective reproduction number Embedded Image

We calculated the reproduction number over time, Rt, and effective reproduction number over time, Embedded Image, using next generation matrix methods (32). The reproduction numbers are calculated for the general population, i.e. excluding care home workers and residents. We define Rt as the average number of secondary infections a case infected at time t would generate in a large entirely susceptible population, and Embedded Image as the average number of secondary infections generated by a case infected at time t would accounting for the finite population size and potential immunity in the population.

To compute the next generation matrix, we calculated the mean duration of infectiousness ΔI, as Embedded Image where parameter and model compartment notations are defined in Table S 2 - Table S 8. For this calculation, the expected durations of stay in compartments were adjusted to account for the discrete-time nature of the model, via calculating the expected number of time-steps (of length dt) spent in a given compartment. Note that if in continuous-time a compartment duration is τ ∼Erlang (k, λ), then the corresponding discrete-time mean duration is: Embedded Image The next generation matrix was calculated as, for (i, j) ∈ {[0,5), …, [75,80), [80 +)}2, Embedded Image where Ni is the total population of group i and Rt is taken to be the dominant eigenvalue of NGM(t), while the effective next generation matrix was calculated as: Embedded Image with Embedded Image taken to be the dominant eigenvalue of NGMeff(t).

1.5 Infection severity

Posterior estimates of severity, namely the infection hospitalisation and infection fatality ratios, were calculated in each group i as follows: Embedded Image Embedded Image Note that for simplicity the notation we use do refer explicitly to the NHS region of interest. We calculated age-aggregated estimates for each region by weighting the age-specific severity estimates by the cumulative incidence in that age group. Aggregate estimates for England were then calculated by weighting the region-specific estimates by the regional attack rates.

1.6 Compartmental model equations

To clearly illustrate the model dynamics, we describe a deterministic version of the model in differential equations (22)-(57), followed by the stochastic implementation used in the analysis. Each compartment is stratified by mixing category i ∈ {[0,5), …, [75,80), [80 +), CHW, CHR}. Full definitions of compartments and model parameters are set out in Table S 2 - Table S 8. Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image We used the tau-leap method (33) to create a stochastic, time-discretised version of the model described in equations (58-162), taking four update steps per day. The process was initialised with ten asymptomatic infectious individuals aged 15-19 on the epidemic start date t0, a parameter we estimate. For each time step, the model iterated through the procedure described below. In the following, we introduce a small abuse of notation: for transitions involving multiple onward compartments (e.g transition from compartment E to compartments IA or IC), for conciseness, we write Embedded Image instead of Embedded Image where Embedded Image is a dummy variable counting the number of individuals remaining in compartment E i,2.

Using this convention, transition variables are drawn from the following distributions, with probabilities defined below: Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Model compartments were then updated as follows: Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image

1.7 Observation process

To describe the epidemic in each NHS region, we fitted our model to time series data on hospital admissions, hospital ward occupancy (both in general beds and in ICU beds), deaths in hospitals, deaths in care homes, population serological surveys and PCR testing data (section 1.1 and Table S 1).

1.7.1 Notation for distributions used in this section

If Y ∼ NegBinom(m, k), then Y follows a negative binomial distribution with mean m and shape κ, such that Embedded Image where Γ(x) is the gamma function. The variance of Y is Embedded Image.

If Z ∼ BetaBinom(n, ω, ρ), then Z follows a beta-binomial distribution with size n, mean probability ω and overdispersion parameter ρ, such that Embedded Image where Embedded Image and B(a, b) is the beta function. The mean of Z is nω and the variance is nω (1 − ω)[1 + (n − 1) ρ].

1.7.2 Hospital admissions and new diagnoses in hospital

We represented the daily number of confirmed COVID-19 hospital admissions and new diagnoses for existing hospitalised cases, Yadm (t), as the observed realisations of an underlying hidden Markov process, Xadm (t), defined as: Embedded Image Which was related to the data via a reporting distribution: Embedded Image We allow for overdispersion in the observation process to account for noise in the underlying data streams, for example due to day-of-week effects on data collection. We adopt k = 2 for all NHSE data streams, so that they contribute equal weight to the overall likelihood.

1.7.3 Hospital bed occupancy by confirmed COVID-19 cases

The model predicted general hospital bed occupancy by confirmed COVID-19 cases, Xhosp (t) as: Embedded Image Which was related to the observed daily general bed-occupancy via a reporting distribution: Embedded Image with κhosp = 2 as above.

Similarly, the model predicted ICU bed occupancy by confirmed COVID-19 cases, XICU (t) as: Embedded Image Which was related to the observed daily ICU bed-occupancy via a reporting distribution: Embedded Image with κICU = 2.

1.7.4 Hospital and care homes COVID-19 deaths

The reported number of daily COVID-19 deaths in hospitals, Embedded Image was considered as the observed realisation of an underlying hidden Markov process, Embedded Image, defined as: Embedded Image Which was related to the data via a reporting distribution: Embedded Image with Embedded Image.

Similarly, we represented the reported number of daily COVID-19 deaths in care homes, Embedded Image, as the observed realisations of an underlying hidden Markov process, Embedded Image, defined as: Embedded Image Which was related to the data via a reporting distribution: Embedded Image with Embedded Image.

1.7.5 Serosurveys

We model serological testing of all individuals aged 15-65, and define the resulting number of seropositive and seronegative individuals (were all individuals aged 15-65 to be tested), as: Embedded Image Embedded Image We compared the observed number of seropositive results, Embedded Image, with that predicted by our model, allowing for i) the sample size of each serological survey, Embedded Imageand ii) imperfect sensitivity Embedded Image and specificity Embedded Image of the serological assay: Embedded Image Where: Embedded Image

1.7.6 PCR testing

As described in the data section (1.1), we fitted the model to PCR testing data from two separate sources:

  • pillar 2: the government testing programme, which recommends that individuals with COVID-19 symptoms are tested (34),

  • the REACT-1 study, which aims to quantify the prevalence of SARS-CoV-2 in a random sample of the England population on an ongoing basis (35).

We only use Pillar 2 PCR test results for individuals aged 25 and over (we assume this includes all care home workers and residents). We assume that individuals who get tested through Pillar 2 PCR testing are either newly symptomatic SARS-CoV-2 cases (who will test positive): Embedded Image or non-SARS-CoV-2 cases who have symptoms consistent with COVID-19 (who will test negative): Embedded Image where pNC is the probability of non SARS-CoV-2 cases having symptoms consistent with COVID-19 leading them to seek a PCR test.

We compared the observed number of positive PCR tests, Embedded Image with that predicted by our model, accounting for the number of PCR tests conducted each day under pillar 2, Embedded Image, by calculating the probability of a positive PCR result (assuming perfect sensitivity and specificity of the PCR test): Embedded Image People may seek PCR tests for many reasons and thus the pillar 2 data are subject to competing biases. We therefore allowed for an over-dispersion parameter Embedded Image, which we fitted separately for each region in the modelling framework: Embedded Image We incorporated the REACT-1 PCR testing data into the likelihood analogously to the serology data, by considering the model-predicted number of PCR-positives Embedded Image and PCR-negatives Embedded Image, were all individuals aged over five and not resident in a care home to be tested: Embedded Image Embedded Image We compared the daily number of positive results observed in REACT-1, Embedded Image, given the number of people tested on that day, Embedded Image, to our model predictions, by calculating the probability of a positive result, assuming perfect sensitivity and specificity of the REACT-1 assay: Embedded Image Embedded Image The overall likelihood function was then calculated as the product of the likelihoods of the individual observations.

1.8 Bayesian inference and model fitting

A closed-form expression of the likelihood of the observed data given the model and its parameters was not analytically tractable, so we used particle filtering methods to obtain an unbiased estimate of the likelihood which can be efficiently sampled from (36). Where appropriate, we used estimates from the literature to set model parameters at fixed values. We limited the parameters being inferred to just those with particular epidemiological interest, or with large uncertainty in existing literature.

The model was fitted independently to each NHS region. For each NHS region, we aimed to infer the values of 26 model parameters:

  • the local epidemic start-date, t0;

  • thirteen transmission rates at different time points β1, …, β12;

  • three parameters governing transmission to and within care homes m mCHW, mCHR, ϵ;

  • the probability of symptomatic individuals developing serious disease requiring hospitalisation, Embedded Image, for the group with the largest probability;

  • the probability of a care home resident dying in a care home if they have severe disease requiring hospitalisation, Embedded Image;

  • the probability, at the start of the pandemic, of a patient being admitted to ICU after hospitalisation, Embedded Image, for the group with the largest probability;

  • the probabilities, at the start of the pandemic, of dying in a hospital general ward, Embedded Image, in the ICU, Embedded Image, and in a stepdown ward following ICU, Embedded Image, for the groups with the largest probability;

  • the multiplier for hospital mortality after improvement in care, μH;

  • the multiplier for probability of admission to ICU after improvement in care, μICU;

  • the daily proportion pNC, of the population seeking to get tested for an infection of SARS-Cov-2 following COVID-19 like symptoms and the overdispersion of the corresponding observation distribution Embedded Image.

We used particle Monte Carlo Markov Chain (pMCMC) methods (37), implementing a particle marginal Metropolis-Hastings algorithm with a bootstrap particle filter (38) with 96 particles (for sufficient variance in likelihood and a convenient multiple of number of available CPU cores for efficiency), to obtain a sample from the posterior distribution of the model parameters given the observed data. If the expected values of count distributions are zero when observed values are non-zero, this results in particles of zero weight, which can lead to the particle filter estimating the marginal likelihood to be 0. Therefore, to get a small but non-zero weight for each particle at every observation, within our particle filter likelihood we add a small amount of noise (exponentially distributed with mean 10−6) to count values from the model.

Within our particle filter we add small amounts of exponentially-distributed noise (with mean 10−6) to model outputs prior to calculating likelihood weights to avoid particles of zero weight, instead resulting in small but non-zero weights.

We implemented our model and parameter inference in an R package, sircovid (39), available at https://mrc-ide.github.io/sircovid, which uses two further R packages, dust to run the model in efficient compiled code and mcstate to implement the pMCMC sampler using Metropolis-Hastings sampling (40).

At each iteration, the sampler proposes an update to the joint distribution of parameters. These proposals are generated from multivariate Gaussian densities centred on the current parameter values, and with a covariance structure chosen to facilitate efficient mixing of the Markov chain. We specified reflecting boundaries for the proposal kernel to ensure that the proposed parameters are both epidemiologically and mathematically plausible and retain symmetry in the proposals.

For each regional fit, eight parallel chains of the pMCMC were run for 11,000 iterations, with the first 1,000 discarded as burn-in, and a 1/80 thinning. We assessed convergence visually.

1.9 Prior distributions and parameter calibration

1.9.1 Risk of hospital admission

In our Bayesian inference framework, we estimate Embedded Image, the probability of hospital admission for symptomatic cases in the group (across all ages and CHW and CHR) with the largest probability of hospital admission. However, we fix the relative probability of hospital admission for the other age groups, Embedded Image, defined so that Embedded Image, with Embedded Image in the group with largest probability of hospital admission.

In this section we explain how the values of Embedded Image were chosen. We used two sources of information, an individual-level and an aggregated dataset. On the one hand, the COVID-19 Hospitalisation in England Surveillance System (CHESS) is a daily, confidential line list containing highly detailed information on patients admitted to hospital with confirmed COVID-19 (see following section 1.9.2 for further details). On the other hand, the Government’s Coronavirus Dashboard is an aggregated, publicly available situation report updated daily. Amongst other data, it provides updates on the number of daily admissions and hospital occupancy by devolved nation and, for England, by NHS region. We found the demography of hospitalisation in CHESS to be biased toward older patients compared to Dashboard data (Figure S3). We thus undertook a two-step approach to infer the demographic composition of COVID-19 hospitalisations across England.

Firstly, we derived an initial approximation of Embedded Image by dividing the total number of hospital admissions for age group i in CHESS over the total number of positive PCR tests (Pillar 2) for i. Both data sources were censored to include patients admitted to hospital or with a specimen data (i.e. the date the test was taken), respectively, between March 1 and December 2, 2020. We ran our full inference framework using this initial approximation for Embedded Image and observed its fit to the demographic composition of admissions from the data.

As a second step, we refined our initial approximations of Embedded Image over a series of iterations of our inference framework, by drawing the modelled Embedded Image and observed Embedded Image proportion of admissions for each age group (i.e. admissions in age group i divided by all admissions) and using it to derive a re-scaling factor for a new proposal for Embedded Image as follows: Embedded Image This process was repeated to obtain a close approximation to the observed proportion of admissions by age and region (Figure S3). A key strength of our approach is that we did not overfitted demography by individual regions. Rather, by assuming Embedded Image to be independent of geographic region, we allowed our inference framework to derive the number of admissions for each five-year age band i solely based on Embedded Image, the demographic composition of the NHS region and inferred epidemic parameters, such as Rt.

Figure S 3:
  • Download figure
  • Open in new tab
Figure S 3:

Proportion of admissions by age. a) Comparison of model outputs to data from the Government’s Coronavirus Dashboard, aggregated to five broad age categories. b) Age spline fitted (red) to Government’s Coronavirus Dashboard, with age categories disaggregated to five-year bands. The fitted spline (red) was used as input parameters for the probability of hospitalisation by age.

1.9.2 Severity and hospital progression

We also performed extensive preliminary analysis to inform the age-structure of progression parameters within hospital. Data from the COVID-19 Hospitalisation in England Surveillance System (CHESS) were used to fit a simple model of patient clinical progression in hospital. The model structure was designed to mirror the within-hospital component of the wider mechanistic transmission model, but without the complexities arising from unknown admission dates and with greater detail on trends with age (Figure S 4).

Figure S 4:
  • Download figure
  • Open in new tab
Figure S 4:

Directed Acyclic Graph of the hospital pathways fitted to CHESS data, which mirror the model structure described in Figure S 2, but with all parameters varying with age and not over time.

CHESS data consists of a line list of daily individual patient-level data on COVID-19 infection in persons requiring hospitalisation, including demographic and clinical information on severity and outcomes. Data were filtered to patients admitted between 18th March and 31st May 2020 (inclusive), with subsequent progression events possible up until 25th Nov 2020.

This gave >5 months for outcomes to complete, and hence justified filtering to patients with resolved outcomes only. The length of stay in each state was taken as the difference between the registered dates of entering and leaving each hospital ward. Lengths of stay were assumed to follow Erlang distributions, as in the wider model, with a distinct mean and shape parameter for each state. Specifically, the probability of being in state Embedded Image for n ∈ ℕ0 days was taken as the integral over day n of the Erlang distribution with mean mX and shape sX: Embedded Image For a patient of age a, this was combined with the probability of their path through the hospital progression model, taken as the product of the individual transition probabilities at each bifurcation, i.e. values taken from pz(a) for Z ∈ {ICU, HD, ICUD, WD}. Transition probabilities were modelled as functions of age using logistic-transformed cubic splines. Knots were defined at coordinates Embedded Image, where xi values were fixed at {0, 20, 40, 60, 80, 100, 120} and Embedded Image were free parameters to be estimated. The complete spline, yz(a) for a ∈ 0: 120, was obtained from these knots using standard expressions for cubic spline interpolation. Finally, transition probabilities were obtained from the raw yz(a) values using the logistic transformation:.Embedded Image.

In total there were 44 free parameters in the within-hospital progression model: 8 mean length of stay parameters, 8 length of stay shape parameters and 4 × 7 transition probability spline nodes (Figure S 4, Table S 4).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 4:

Descriptions of all states and transitions in the simplified hospital progression model fitted to CHESS data.

All parameters of the hospital progression model were given priors (Table S 5) and estimated within a Bayesian framework. All length of stay parameters were given uniform priors over a plausible range of values. For transition probabilities, the first spline node Embedded Image was given a prior that corresponded to a uniform distribution after logistic transformation, and subsequent spline nodes were given a multivariate normal prior to apply a smoothing constraint to the spline. Parameters were estimated jointly via MCMC using the custom package markovid v1.5.0 (41), which uses the random-walk Metropolis-Hastings algorithm to draw from the joint posterior distribution. MCMC was run for 1000 burn-in iterations and 100,000 sampling iterations replicated over 10 independent chains. Convergence was assessed via the Gelman-Rubin diagnostic (all parameters had potential scale reduction factor <1.1) and sampling sufficiency was assessed by visualising posterior distributions and by effective sample size (ESS) calculations (all parameters had ESS >100,000).

View this table:
  • View inline
  • View popup
Table S 5:

Priors on all length of stay distributions and transition probability splines. Uniformcont(a, b) denotes the continuous uniform distribution, and Uniformdisc(a, b) the discrete uniform distribution between a and b (inclusive).

Parameter estimates (posterior medians) were passed to the wider mechanistic transmission model as fixed values (Figure S 5). For transition probabilities, the full age-spline (Figure 3, main text) was aggregated to 5-year age groups and normalised by the largest value to define the relative risk with age. The absolute risk in the mechanistic transmission model was obtained by multiplying the relative risk by region-specific scaling factors that were fitted as free parameters in the pMCMC. Hence, the preliminary analysis of CHESS data was used to inform trends of severity with age, but not the absolute probability of progression through the hospital states, which was informed by the Government’s Coronavirus Dashboard data.

For the wider mechanistic transmission model, we used Beta distributions for the priors of the various fitted probabilities regarding hospitalisation. The priors for Embedded Image and Embedded Image were all informed by the fitting to CHESS data by taking the median fitted value for the prior mean, which we halve in the case of Embedded Image to account for CHESS being biased to more severe patients. The prior distributions are then calibrated so that the lower bound of the 95% confidence interval is 0.1 lower than the prior mean. For Embedded Image and Embedded Image, we assume prior means of 0.75 and calibrate the prior so that the lower bound of the 95% confidence interval is 0.2 lower than the mean. For the multipliers for hospital mortality after improvement in care, μH, and for probability of admission to ICU after improvement in care, μICU we used uninformative, [0,1] priors.

Figure S 5:
  • Download figure
  • Open in new tab
Figure S 5:

Posterior 95% credible intervals of length of stay mean (left) and shape parameters (right).

1.9.3 Serosurveys

To keep serology parameters consistent between all regions we used estimates from the literature to fix the parameters of the seroconversion process. An alternative would have been to use these estimates as priors within a hierarchical model where some parameters would be shared between regions, but this would be much more involved computationally.

As described in section 1.3.2, the time to seroconversion from leaving the Ei compartment is modelled by an exponential distribution time spent in Embedded Image with a proportion Embedded Image ultimately seroconverting and moving to Embedded Image and the remaining staying negative and moving to Embedded Image.

We fixed Embedded Image to 0.85 based on the estimate of 85% of infections becoming detectably seropositive with the EUROIMMUN assay used in the NHSBT serological surveys (42). The specificity of the serology test is Embedded Image fixed to 0.99 also from (42). Finally, the sensitivity of serology test Embedded Image is assumed to be 1 as it is non-distinguishable from the time varying seroconversion process (Table S7).

1.9.4 PCR positivity

As for other compartments, we modelled the duration of SARS-CoV-2 PCR-positivity after symptom onset using an Erlang distribution τ ∼Erlang(k, γ), with k successive compartments and a total mean time spent of Embedded Image and variance Embedded Image.

We estimated the parameters of this distribution from Omar et al. (16), which reported the cumulative distribution of duration of PCR positivity in 523 individuals with mild COVID-19 disease in home quarantine in a German region. We performed a survival analysis using a gamma-accelerated failure time model fitted to their data, from which we estimated the mean and variance of the time from symptom onset to PCR negativity. This was used to derive values of k and U shown in Table S 2.

1.9.5 Local start date of the epidemic

The start date of the epidemic for each region is assumed to have a uniform prior on the dates from 1st January 2020 to 15th March 2020, inclusive –with the latter date corresponding to the last date before the data begin.

1.9.6 Time-varying transmission rates

We set priors for the transmission rates β1, …, β12 to reflect a Gamma distribution for the reproduction number Rt with a reasonable 95% confidence interval a priori. To obtain a prior for the corresponding βk, we then scale by a factor of 0.0241 (given other parameter values, βk = 0.0241 would correspond approximately to Rt = 1). The 95% ranges for R0 we used are (i) (2.5, 3.5) at the onset of the epidemic (corresponding to β1); and then Rt (ii) (0.4, 3.5) at announcement of the first lockdown (corresponding to β2); and (iii) (0.4, 3) from the implementation of the first lockdown onwards (corresponding to β3, …, β12). The values are consistent with the values of the COMIX study (43).

1.9.7 Transmission within care homes

For the transmission between care home workers and residents, mCHW, and transmission among care home residents, mCHR, we used a prior distributions reflecting that these are person-to-person infectious contact rates and thus should be scaled according to regional care home demography. We then used a Gamma distribution with shape 5 and mean Embedded Image for both of these parameters (recall that we assume there is a 1-to-1 ratio of care home workers to residents in each region, so NCHW = NCHR).

For the parameter governing the reduction in contacts between the general population and care home residents, ϵ, we used an uninformative U[0,1] prior.

1.9.8 Parameters relating to Pillar 2 testing

For both the parameters pNC and Embedded Image, we used uninformative U[0,1] priors.

View this table:
  • View inline
  • View popup
Table S 6:

Inferred model parameter notations, prior and posterior distributions. Note that Γ (a, b) here refers to a Gamma distribution with shape a and scale b (such that the mean is ab), and ‘(a, b) refers to a Beta distribution with shape parameters a and b (such that the mean is a/(a + b)).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 7:

Fixed model parameters (age / care home scaling factors are shown separately in Table S 8).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 8:

Age / care-home scaling factors

2 Supplementary Results

2.1 Model fitting

Figure S 6:
  • Download figure
  • Open in new tab
Figure S 6:

Model fits by region to PCR positivity for individuals aged >25 years (top row) and PCR positivity from the REACT-1 study (bottom row). The points show the data and bars the 95% CI. The solid line the median model fit and the shaded area the 95% CrI.

Figure S 7:
  • Download figure
  • Open in new tab
Figure S 7:

Model fits to daily hospital deaths (top row), daily care home deaths (second row), ICU bed occupancy (third row), general bed occupancy (fourth row), all hospital beds (fifth row), and all daily admissions (bottom row) by region (columns). The points show the data, the solid line the median model fit and the shaded area the 95% CrI.

2.2 Severity estimates

Figure S 8:
  • Download figure
  • Open in new tab
Figure S 8:

Fits to CHESS data broken down into one-year age bands. Blue ribbons show the 95% CrI of the fitted spline, black circles and vertical segments give the raw mean and 95% CI from the data (exact binomial).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S 9:

Age-stratified estimates of disease severity (*to 2sf, ^to 3dp)

Supplementary counterfactual analysis

Figure S 9:
  • Download figure
  • Open in new tab
Figure S 9:

Counterfactual intervention scenarios in each England NHS Region: Panel A1-7 impact of locking down one-week earlier Panel B1-7 impact of locking down one week later; Panel C1-7 impact of relaxing lockdown restrictions two weeks earlier. Panel D1-7 impact of relaxing lockdown restrictions two weeks later; Panel E1-7 impact of 50% less contact between care home residents and the general population; Panel F1-7 impact of 50% more contact between care home residents and the general population.

Acknowledgements

We thank all the colleagues at Public Health England (PHE) and frontline health professionals who have driven and continue to drive the daily response to the epidemic, but also for providing the necessary data to inform this study. This work would not have been possible without their dedication and expertise. The use of pillar 2 PCR testing data was made possible thanks to PHE colleagues and we extend our thanks to Dr Nick Gent and Dr André Charlett for facilitation and their insights into these data. The use of serological data was made possible by colleagues at PHE Porton Down, Colindale, and the NHS Blood Transfusion Service. We are particularly grateful to Dr Gayatri Amirthalingam and Prof Nick Andrews for helpful discussions around these data. We also thank the entire Imperial College London Covid-19 Response Team for their support and feedback throughout. This work was supported by the NIHR HPRU in Modelling and Health Economics, a partnership between PHE, Imperial College London and LSHTM (grant code NIHR200908). We acknowledge funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 programme supported by the European Union.

References

  1. 1.↵
    GOV.UK, Coronavirus (COVID-19) in the UK (2020), (available at https://coronavirus.data.gov.uk/details/download).
  2. 2.↵
    E. Mahase, Covid-19: UK records first death, as world’s cases exceed 100 000. BMJ. 368 (2020), doi:10.1136/bmj.m943.
    OpenUrlFREE Full Text
  3. 3.↵
    GOV.UK, Prime Minister’s statement on coronavirus (COVID-19): 12 March 2020 - GOV.UK, (available at https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-12-march-2020).
  4. 4.↵
    GOV.UK, Prime Minister’s statement on coronavirus (COVID-19): 20 March 2020 - GOV.UK, (available at https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-20-march-2020).
  5. 5.↵
    GOV.UK, Prime Minister’s statement on coronavirus (COVID-19): 22 March 2020 - GOV.UK, (available at https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-22-march-2020).
  6. 6.↵
    B. Jeffrey, C. E. Walters, K. E. C. Ainslie, O. Eales, C. Ciavarella, S. Bhatia, S. Hayes, M. Baguelin, A. Boonyasiri, N. F. Brazeau, G. Cuomo-Dannenburg, R. G. FitzJohn, K. Gaythorpe, W. Green, N. Imai, T. A. Mellan, S. Mishra, P. Nouvellet, H. J. T. Unwin, R. Verity, M. Vollmer, C. Whittaker, N. M. Ferguson, C. A. Donnelly, S. Riley, Anonymised and aggregated crowd level mobility data from mobile phones suggests that initial compliance with covid-19 social distancing interventions was high and geographically consistent across the UK. Wellcome Open Res. 5, 1–10 (2020).
    OpenUrl
  7. 7.↵
    HM Revenue & Customs, Get a discount with the Eat Out to Help Out Scheme. www.gov.uk (2020).
  8. 8.↵
    GOV.UK, Rule of six comes into effect to tackle coronavirus - GOV.UK, (available at https://www.gov.uk/government/news/rule-of-six-comes-into-effect-to-tackle-coronavirus).
  9. 9.↵
    GOV.UK, Prime Minister announces new local COVID Alert Levels - GOV.UK, (available at https://www.gov.uk/government/news/prime-minister-announces-new-local-covid-alert-levels).
  10. 10.↵
    GOV.UK, Prime Minister announces new national restrictions - GOV.UK, (available at https://www.gov.uk/government/news/prime-minister-announces-new-national-restrictions).
  11. 11.↵
    NHS Digital, SGSS and CHESS data - NHS Digital, (available at https://digital.nhs.uk/about-nhs-digital/corporate-information-and-documents/directions-and-data-provision-notices/data-provision-notices-dpns/sgss-and-chess-data).
  12. 12.↵
    NHS, Overview | COVID-19 rapid guideline: critical care in adults | Guidance | NICE.
  13. 13.↵
    N. F. Brazeau, R. Verity, S. Jenks, H. Fu, C. Whittaker, P. Winskill, I. Dorigatti, P. Walker, S. Riley, R. P. Schnekenberg, H. Hoeltgebaum, T. A. Mellan, S. Mishra, H. T. Juliette Unwin, O. J. Watson, Z. M. Cucunubá, M. Baguelin, L. Whittles, S. Bhatt, A. C. Ghani, N. M. Ferguson, L. C. Okell, Infection Fatality Ratio: Estimates from Seroprevalence, doi:10.25561/83545.
    OpenUrlCrossRef
  14. 14.↵
    R. A. Armstrong, A. D. Kane, T. M. Cook, Outcomes from intensive care in patients with COVID-19: a systematic review and meta-analysis of observational studies. Anaesthesia. 75, 1340–1349 (2020).
    OpenUrlCrossRefPubMed
  15. 15.↵
    Public Health England, Sero-surveillance of COVID-19 - GOV.UK, (available at https://www.gov.uk/government/publications/national-covid-19-surveillance-reports/sero-surveillance-of-covid-19).
  16. 16.↵
    F. J. Ibarrondo, J. A. Fulcher, D. Goodman-Meza, J. Elliott, C. Hofmann, M. A. Hausner, K. G. Ferbas, N. H. Tobin, G. M. Aldrovandi, O. O. Yang, Rapid Decay of Anti–SARS-CoV-2 Antibodies in Persons with Mild Covid-19. N. Engl. J. Med., 1–2 (2020).
  17. 17.↵
    L. H. Nguyen, D. A. Drew, M. S. Graham, A. D. Joshi, C. G. Guo, W. Ma, R. S. Mehta, E. T. Warner, D. R. Sikavi, C. H. Lo, S. Kwon, M. Song, L. A. Mucci, M. J. Stampfer, W. C. Willett, A. H. Eliassen, J. E. Hart, J. E. Chavarro, J. W. Rich-Edwards, R. Davies, J. Capdevila, K. A. Lee, M. N. Lochlainn, T. Varsavsky, C. H. Sudre, M. J. Cardoso, J. Wolf, T. D. Spector, S. Ourselin, C. J. Steves, A. T. Chan, C. M. Albert, G. Andreotti, B. Bala, B. A. Balasubramanian, L. E. Beane-Freeman, J. S. Brownstein, F. J. Bruinsma, J. Coresh, R. Costa, A. N. Cowan, A. Deka, S. L. Deming-Halverson, M. Elena Martinez, M. E. Ernst, J. C. Figueiredo, P. Fortuna, P. W. Franks, L. B. Freeman, C. D. Gardner, I. M. Ghobrial, C. A. Haiman, J. E. Hall, J. H. Kang, B. Kirpach, K. C. Koenen, L. D. Kubzansky, J. V Lacey, L. Le Marchand, X. Lin, P. Lutsey, C. R. Marinac, M. E. Martinez, R. L. Milne, A. M. Murray, D. Nash, J. R. Palmer, A. V Patel, E. Pierce, M. M. Robertson, L. Rosenberg, D. P. Sandler, S. H. Schurman, K. Sewalk, S. V Sharma, C. J. Sidey-Gibbons, L. Slevin, J. W. Smoller, C. J. Steves, M. I. Tiirikainen, S. T. Weiss, L. R. Wilkens, F. Zhang, Risk of COVID-19 among front-line health-care workers and the general community: a prospective cohort study. Lancet Public Heal. 5, e475.-e483 (2020).
    OpenUrl
  18. 18.↵
    GOV.UK, COVID-19: number of outbreaks in care homes - management information, (available at https://www.gov.uk/government/statistical-data-sets/covid-19-number-of-outbreaks-in-care-homes-management-information#history).
  19. 19.
    S. N. Ladhani, J. Y. Chow, R. Janarthanan, J. Fok, E. Crawley-Boevey, A. Vusirikala, E. Fernandez, M. S. Perez, S. Tang, K. Dun-Campbell, E. W. Evans, A. Bell, B. Patel, Z. Amin-Chowdhury, F. Aiano, K. Paranthaman, T. Ma, M. Saavedra-Campos, R. Myers, J. Ellis, A. Lackenby, R. Gopal, M. Patel, C. Brown, M. Chand, K. Brown, M. E. Ramsay, S. Hopkins, N. Shetty, M. Zambon, Investigation of SARS-CoV-2 outbreaks in six care homes in London, April 2020. EClinicalMedicine. 26, 100533 (2020).
    OpenUrl
  20. 20.↵
    J. K. Burton, G. Bayne, C. Evans, F. Garbe, D. Gorman, N. Honhold, D. McCormick, R. Othieno, J. E. Stevenson, S. Swietlik, K. E. Templeton, M. Tranter, L. Willocks, B. Guthrie, Evolution and effects of COVID-19 outbreaks in care homes: a population analysis in 189 care homes in one geographical region of the UK. Lancet Heal. Longev. 1, e21–e31 (2020).
    OpenUrl
  21. 21.↵
    L. J. Strausbaugh, S. R. Sukumar, C. L. Joseph, Infectious disease outbreaks in nursing homes: An unappreciated hazard for frail elderly persons. Clin. Infect. Dis. 36, 870–876 (2003).
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    T. Inns, D. Wilson, P. Manley, J. P. Harris, S. J. O’Brien, R. Vivancos, What proportion of care home outbreaks are caused by norovirus? An analysis of viral causes of gastroenteritis outbreaks in care homes, North East England, 2016-2018. BMC Infect. Dis. 20, 2 (2019).
    OpenUrl
  23. 23.↵
    P. N. Perez-Guzman, A. Daunt, S. Mukherjee, P. Crook, R. Forlano, M. D. Kont, A. Løchen, M. Vollmer, P. Middleton, R. Judge, C. Harlow, A. Soubieres, G. Cooke, P. J. White, T. B. Hallett, P. Aylin, N. Ferguson, K. Hauck, M. R. Thursz, S. Nayagam, Clinical characteristics and predictors of outcomes of hospitalized patients with COVID-19 in a multi-ethnic London NHS Trust: a retrospective cohort study. Clin. Infect. Dis., 1–11 (2020).
  24. 24.↵
    A. B. Docherty, E. M. Harrison, C. A. Green, H. E. Hardwick, R. Pius, L. Norman, K. Holden, J. M. Read, F. Dondelinger, G. Carson, L. Merson, J. Lee, D. Plotkin, L. Sigfrid, S. Halpin, C. Jackson, C. Gamble, P. W. Horby, J. S. Nguyen-Van-Tam, A. Ho, C. D. Russell, J. Dunning, P. J. M. Openshaw, J. K. Baillie, M. G. Semple, Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: Prospective observational cohort study. BMJ. 369, 1–12 (2020).
    OpenUrlCrossRefPubMed
  25. 25.
    S. Tabata, K. Imai, S. Kawano, M. Ikeda, T. Kodama, K. Miyoshi, H. Obinata, S. Mimura, T. Kodera, M. Kitagaki, M. Sato, S. Suzuki, T. Ito, Y. Uwabe, K. Tamura, Clinical characteristics of COVID-19 in 104 people with SARS-CoV-2 infection on the Diamond Princess cruise ship: a retrospective analysis. Lancet Infect. Dis. 20, 1043–1050 (2020).
    OpenUrlPubMed
  26. 26.↵
    P. Vetter, D. L. Vu, A. G. L’Huillier, M. Schibler, L. Kaiser, F. Jacquerioz, Clinical features of covid-19. BMJ. 369, 1–2 (2020).
    OpenUrlCrossRefPubMed
  27. 27.↵
    R. Verity, L. C. Okell, I. Dorigatti, P. Winskill, C. Whittaker, N. Imai, G. Cuomo-Dannenburg, H. Thompson, P. G. T. Walker, H. Fu, A. Dighe, J. T. Griffin, M. Baguelin, S. Bhatia, A. Boonyasiri, A. Cori, Z. Cucunubá, R. FitzJohn, K. Gaythorpe, W. Green, A. Hamlet, W. Hinsley, D. Laydon, G. Nedjati-Gilani, S. Riley, S. van Elsland, E. Volz, H. Wang, Y. Wang, X. Xi, C. A. Donnelly, A. C. Ghani, N. M. Ferguson, Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect. Dis. (2020), doi:10.1016/S1473-3099(20)30243-7.
    OpenUrlCrossRefPubMed
  28. 28.↵
    The RECOVERY Collaborative Group, Dexamethasone in Hospitalized Patients with Covid-19 — Preliminary Report. N. Engl. J. Med., 1–11 (2020).
  29. 29.↵
    L. I. Horwitz, S. A. Jones, R. J. Cerfolio, F. Francois, J. Greco, B. Rudy, C. M. Petrilli, Trends in COVID-19 Risk-Adjusted Mortality Rates. J. Hosp. Med. 23, 2020 (2020).
    OpenUrl
  30. 30.↵
    A. W. D. Edridge, J. Kaczorowska, A. C. R. Hoste, M. Bakker, M. Klein, K. Loens, M. F. Jebbink, A. Matser, C. M. Kinsella, P. Rueda, M. Ieven, H. Goossens, M. Prins, P. Sastre, M. Deijs, L. van der Hoek, Seasonal coronavirus protective immunity is shortlasting. Nat. Med. (2020), doi:10.1038/s41591-020-1083-1.
    OpenUrlCrossRefPubMed
  31. 31.↵
    S. M. Kissler, C. Tedijanto, E. Goldstein, Y. H. Grad, M. Lipsitch, Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science (80-.). 368, 860–868 (2020).
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    L. C. Okell, R. Verity, O. J. Watson, S. Mishra, P. Walker, C. Whittaker, A. Katzourakis, C. A. Donnelly, S. Riley, A. C. Ghani, A. Gandy, S. Flaxman, N. M. Ferguson, S. Bhatt, Correspondence Have deaths from COVID-19 in Europe plateaued due to herd. Lancet. 395, e110.-e111 (2020).
    OpenUrl
  33. 33.↵
    Academy of Medical Science, Preparing for a challenging winter 2020/21, 79 (2020).
  34. 34.↵
    NHS England and NHS Improvement, “Hospital Onset Covid-19: IPC evidence from recent survey and next steps.”
  35. 35.↵
    G. Iacobucci, Covid-19: Doctors sound alarm over hospital transmissions. BMJ. 369, m2013 (2020).
    OpenUrlFREE Full Text
  36. 36.
    D. Buitrago-Garcia, D. Egli-Gany, M. J. Counotte, S. Hossmann, H. Imeri, A. M. Ipekci, G. Salanti, N. Low, Occurrence and transmission potential of asymptomatic and presymptomatic SARS-CoV-2 infections: A living systematic review and metaanalysis. PLoS Med. 17, e1003346 (2020).
    OpenUrlCrossRefPubMed
  37. 37.
    S. Riley, C. E. Walters, H. Wang, O. Eales, K. E. C. Ainslie, C. Atchison, C. Fronterre, P. J. Diggle, D. Ashby, C. A. Donnelly, G. Cooke, W. Barclay, H. Ward, A. Darzi, P. Elliott, medRxiv, in press, doi:10.1101/2020.12.15.20248244.
    OpenUrlAbstract/FREE Full Text
  38. 38.
    Office for National Statistics, Office for National Statistics, (available at https://www.ons.gov.uk/).
  39. 39.
    Care Quality Commission, [ARCHIVED CONTENT] UK Government Web Archive - The National Archives, (available at https://webarchive.nationalarchives.gov.uk/20200605160439/ https://www.cqc.org.uk/files/cqc-care-directory-filters-1-june-2020).
  40. 40.
    GOV.UK, “Care Homes Analysis Background” (2020).
  41. 41.
    Age UK, Later Life in the United Kingdom 2019, (available at https://www.ageuk.org.uk/globalassets/age-uk/documents/reports-and-publications/later_life_uk_factsheet.pdf).
  42. 42.
    J. Mossong, N. Hens, M. Jit, P. Beutels, K. Auranen, R. Mikolajczyk, M. Massari, S. Salmaso, G. S. Tomba, J. Wallinga, J. Heijne, M. Sadkowska-Todys, M. Rosinska, W. J. Edmunds, Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 5, 381–391 (2008).
    OpenUrlCrossRefWeb of Science
  43. 43.
    S. A. Lauer, K. H. Grantz, Q. Bi, F. K. Jones, Q. Zheng, H. R. Meredith, A. S. Azman, N. G. Reich, J. Lessler, The incubation period of coronavirus disease 2019 (CoVID-19) from publicly reported confirmed cases: Estimation and application. Ann. Intern. Med. 172, 577–582 (2020).
    OpenUrlCrossRefPubMed
  44. 44.
    Q. Bi, Y. Wu, S. Mei, C. Ye, X. Zou, Z. Zhang, X. Liu, L. Wei, S. A. Truelove, T. Zhang, W. Gao, C. Cheng, X. Tang, X. Wu, Y. Wu, B. Sun, S. Huang, Y. Sun, J. Zhang, T. Ma, J. Lessler, T. Feng, Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. Lancet Infect Dis (2020), doi:10.1016/s1473-3099(20)30287-5.
    OpenUrlCrossRefPubMed
  45. 45.
    M. Bernabeu-Wittel, J. E. Ternero-Vega, P. Díaz-Jiménez, C. Conde-Guzmán, M. D. Nieto-Martín, L. Moreno-Gaviño, J. Delgado-Cuesta, M. Rincón-Gómez, L. Giménez-Miranda, D. Navarro-Amuedo, M. M. Muñoz-García, S. Calzón-Fernández, M. Ollero-Baturone, Death risk stratification in elderly patients with covid-19. A comparative cohort study in nursing homes outbreaks. Arch. Gerontol. Geriatr. 91, 104240 (2020).
    OpenUrl
  46. 46.
    S. Omar, C. Bartz, S. Becker, S. Basenach, S. Pfeifer, C. Trapp, H. Hamm, H. C. Schlichting, M. Friederichs, U. Koch, C. Jestrabek, E. Hilger, M. Vogt, K. Jahn, S. Chen, T. Barnighausen, P. Zanger, Duration of SARS-CoV-2 RNA detection in COVID-19 patients in home isolation, Rhineland-Palatinate, Germany, 2020 - an interval-censored survival analysis. Eurosurveillance. 25, 1–8 (2020).
    OpenUrl
  47. 47.
    B. Benny, G. Amandine, P. Kc, H. Sarah, M. Abby, C. Caitlin, S. Van, L.-S. James, O. Affiliations, Quantifying antibody kinetics and RNA shedding during early-phase SARS-CoV-2 infection, doi:10.1101/2020.05.15.20103275.
    OpenUrlAbstract/FREE Full Text
  48. 48.
    S. Funk, Socialmixr: Social Mixing Matrices for Infectious Disease Modelling (2018).
  49. 49.
    GOV.UK, Prime Minister’s statement on coronavirus (COVID-19): 25 March 2020 - GOV.UK, (available at https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-25-march-2020).
  50. 50.
    GOV.UK, Prime Minister’s statement on coronavirus (COVID-19): 11 May 2020 - GOV.UK, (available at https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-11-may-2020).
  51. 51.
    GOV.UK, Prime Minister sets out timeline for retail to reopen in June - GOV.UK, (available at https://www.gov.uk/government/news/prime-minister-sets-out-timeline-for-retail-to-reopen-in-june).
  52. 52.
    GOV.UK, Pubs, restaurants and hairdressers to reopen from 4 July - GOV.UK, (available at https://www.gov.uk/government/news/pubs-restaurants-and-hairdressers-to-reopen-from-4-july).
  53. 53.
    GOV.UK, Eat Out to Help Out launches today –with government paying half on restaurant bills - GOV.UK, (available at https://www.gov.uk/government/news/eat-out-to-help-out-launches-today-with-government-paying-half-on-restaurant-bills).
  54. 54.
    GOV.UK, Schools and colleges to reopen in full in September - GOV.UK, (available at https://www.gov.uk/government/news/schools-and-colleges-to-reopen-in-full-in-september).
  55. 55.
    O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 28, 365–382 (1990).
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.
    D. T. Gillespie, Approximate accelerated stochastic simulation of chemically reacting systems. J. Chem. Phys. 115, 1716–1733 (2001).
    OpenUrlCrossRefWeb of Science
  57. 57.
    Department of Health and Social Care, COVID-19 testing data: methodology note. www.gov.uk (2020).
  58. 58.
    I. M. C. Martin, C. A. Ison, D. M. Aanensen, K. A. Fenton, B. G. Spratt, Rapid Sequence-Based Identification of Gonococcal Transmission Clusters in a Large Metropolitan Area. J. Infect. Dis. 189, 1497–1505 (2004).
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.
    P. Del Moral, A. Doucet, A. Jasra, Sequential Monte Carlo samplers. J. R. Stat. Soc. Ser. B Stat. Methodol. 68, 411–436 (2006).
    OpenUrlCrossRef
  60. 60.
    C. Andrieu, A. Doucet, R. Holenstein, Particle Markov chain Monte Carlo methods. J. R. Stat. Soc. Ser. B Stat. Methodol. 72, 269–342 (2010).
    OpenUrlCrossRef
  61. 61.
    N. J. Gordon, D. J. Salmond, A. F. M. Smith, Novel approach to nonlinear/non-gaussian Bayesian state estimation. IEE Proceedings, Part F Radar Signal Process. 140, 107–113 (1993).
    OpenUrl
  62. 62.
    M. Baguelin, E. Knock, L. K. Whittles, R. FitzJohn, J. Lees, sircovid (2020).
  63. 63.
    E. S. Knock, L. K. Whittles, P. N. Perez-Guzman, S. Bhatia, F. Guntoro, O. J. Watson, C. Whittaker, N. M. Ferguson, A. Cori, M. Baguelin, R. G. FitzJohn, J. A. Lees, Reproducible parallel inference and simulation of stochastic state space models using odin, dust, and mcstate. Wellcome Open Res. 5, 288 (2020).
    OpenUrl
  64. 64.
    M. Plummer, N. Best, K. Cowles, K. Vines, CODA: Convergence Diagnosis and Output Analysis for MCMC. R News. 6, 7–11 (2006).
    OpenUrlCrossRefPubMed
  65. 65.
    A. Gelman, D. B. Rubin, Inference from iterative simulation using multiple sequences. Stat. Sci. 7, 457–472 (1992).
    OpenUrlCrossRefPubMed
  66. 66.
    R. Verity, R. FitzJohn, mrc-ide/markovid at version 1.5, (available at https://github.com/mrc-ide/markovid/tree/version1.5).
  67. 67.
    C. I. Jarvis, K. Van Zandvoort, A. Gimma, K. Prem, M. Auzenbergs, K. O’Reilly, G. Medley, J. C. Emery, R. M. G. J. Houben, N. Davies, E. S. Nightingale, S. Flasche, T. Jombart, J. Hellewell, S. Abbott, J. D. Munday, N. I. Bosse, S. Funk, F. Sun, A. Endo, Rosello, S. R. Procter, A. J. Kucharski, T. W. Russell, G. Knight, H. Gibbs, Q. Leclerc, B. J. Quilty, C. Diamond, Y. Liu, M. Jit, S. Clifford, C. A. B. Pearson, R. M. Eggo, A. K. Deol, P. Klepac, G. J. Rubin, W. J. Edmunds, Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Med. 18, 1–10 (2020).
    OpenUrlCrossRef
  68. 68.
    E. Lavezzo, E. Franchin, C. Ciavarella, G. Cuomo-Dannenburg, L. Barzon, C. Del Vecchio, L. Rossi, R. Manganelli, A. Loregian, N. Navarin, D. Abate, M. Sciro, S. Merigliano, E. De Canale, M. C. Vanuzzo, V. Besutti, F. Saluzzo, F. Onelia, M. Pacenti, S. G. Parisi, G. Carretta, D. Donato, L. Flor, S. Cocchio, G. Masi, A. Sperduti, L. Cattarino, R. Salvador, M. Nicoletti, F. Caldart, G. Castelli, E. Nieddu, B. Labella, L. Fava, M. Drigo, K. A. M. Gaythorpe, K. E. C. Ainslie, M. Baguelin, S. Bhatt, A. Boonyasiri, O. Boyd, L. Cattarino, C. Ciavarella, H. L. Coupland, Z. Cucunubá, G. Cuomo-Dannenburg, B. A. Djafaara, C. A. Donnelly, I. Dorigatti, S. L. van Elsland, R. FitzJohn, S. Flaxman, K. A. M. Gaythorpe, W. D. Green, T. Hallett, A. Hamlet, D. Haw, N. Imai, B. Jeffrey, E. Knock, D. J. Laydon, T. Mellan, S. Mishra, G. Nedjati-Gilani, P. Nouvellet, L. C. Okell, K. V. Parag, S. Riley, H. A. Thompson, H. J. T. Unwin, R. Verity, M. A. C. Vollmer, P. G. T. Walker, C. E. Walters, H. Wang, Y. Wang, O. J. Watson, C. Whittaker, L. K. Whittles, X. Xi, N. M. Ferguson, A. R. Brazzale, S. Toppo, M. Trevisan, V. Baldo, C. A. Donnelly, N. M. Ferguson, I. Dorigatti, A. Crisanti, Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo’. Nature (2020), doi:10.1038/s41586-020-2488-1.
    OpenUrlCrossRefPubMed
  69. 69.
    NHS England and NHS Improvement, Statistics » COVID-19 Hospital Activity, (available at https://www.england.nhs.uk/statistics/statistical-work-areas/covid-19-hospital-activity/).

References

  1. 1.↵
    Buitrago-Garcia D, Egli-Gany D, Counotte MJ, Hossmann S, Imeri H, Ipekci AM, et al. Occurrence and transmission potential of asymptomatic and presymptomatic SARS-CoV-2 infections: A living systematic review and meta-analysis. PLoS Med. 2020;17(9):e1003346.
    OpenUrlCrossRefPubMed
  2. 2.↵
    GOV.UK. Coronavirus (COVID-19) in the UK [Internet]. 2020 [cited 2020 Dec 3]. Available from: https://coronavirus.data.gov.uk/details/download
  3. 3.↵
    Riley S, Walters CE, Wang H, Eales O, Ainslie KEC, Atchison C, et al. REACT-1 round 7 updated report: regional heterogeneity in changes in prevalence of SARS-CoV-2 infection during the second national COVID-19 lockdown in England. medRxiv. 2020 Dec 16;2020.12.15.20248244.
  4. 4.
    Public Health England. Sero-surveillance of COVID-19 - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/publications/national-covid-19-surveillance-reports/sero-surveillance-of-covid-19
  5. 5.
    NHS Digital. SGSS and CHESS data - NHS Digital [Internet]. [cited 2020 Dec 4]. Available from: https://digital.nhs.uk/about-nhs-digital/corporate-information-and-documents/directions-and-data-provision-notices/data-provision-notices-dpns/sgss-and-chess-data
  6. 6.↵
    Office for National Statistics. Office for National Statistics [Internet]. [cited 2020 Dec 16]. Available from: https://www.ons.gov.uk/
  7. 7.↵
    Care Quality Commission. [ARCHIVED CONTENT] UK Government Web Archive - The National Archives [Internet]. [cited 2020 Dec 16]. Available from: https://webarchive.nationalarchives.gov.uk/20200605160439/ https://www.cqc.org.uk/files/cqc-care-directory-filters-1-june-2020
  8. 8.↵
    GOV.UK. Care Homes Analysis Background. 2020.
  9. 9.↵
    Age UK. Later Life in the United Kingdom 2019 [Internet]. [cited 2020 Dec 16]. Available from: https://www.ageuk.org.uk/globalassets/age-uk/documents/reports-and-publications/later_life_uk_factsheet.pdf
  10. 10.↵
    Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5(3):381–91.
    OpenUrlCrossRefWeb of Science
  11. 11.↵
    Ladhani SN, Chow JY, Janarthanan R, Fok J, Crawley-Boevey E, Vusirikala A, et al. Investigation of SARS-CoV-2 outbreaks in six care homes in London, April 2020. EClinicalMedicine. 2020 Sep 1;26:100533.
    OpenUrl
  12. 12.
    Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, et al. The incubation period of coronavirus disease 2019 (CoVID-19) from publicly reported confirmed cases: Estimation and application. Ann Intern Med. 2020;172(9):577–82.
    OpenUrlCrossRefPubMed
  13. 13.
    Bi Q, Wu Y, Mei S, Ye C, Zou X, Zhang Z, et al. Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. Lancet Infect Dis. 2020;(PG-).
  14. 14.
    Docherty AB, Harrison EM, Green CA, Hardwick HE, Pius R, Norman L, et al. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: Prospective observational cohort study. BMJ. 2020;369(March):1–12.
    OpenUrlCrossRefPubMed
  15. 15.↵
    Bernabeu-Wittel M, Ternero-Vega JE, Díaz-Jiménez P, Conde-Guzmán C, Nieto-Martín MD, Moreno-Gaviño L, et al. Death risk stratification in elderly patients with covid-19. A comparative cohort study in nursing homes outbreaks. Arch Gerontol Geriatr. 2020;91:104240.
    OpenUrl
  16. 16.↵
    Omar S, Bartz C, Becker S, Basenach S, Pfeifer S, Trapp C, et al. Duration of SARS-CoV-2 RNA detection in COVID-19 patients in home isolation, Rhineland-Palatinate, Germany, 2020 - an interval-censored survival analysis. Eurosurveillance. 2020;25(30):1–8.
    OpenUrl
  17. 17.
    Benny B, Amandine G, Kc P, Sarah H, Abby M, Caitlin C, et al. Quantifying antibody kinetics and RNA shedding during early-phase SARS-CoV-2 infection.
  18. 18.↵
    Funk S. Socialmixr: Social Mixing Matrices for Infectious Disease Modelling. The Comprehensive R Archive Network; 2018.
  19. 19.
    GOV.UK. Prime Minister’s statement on coronavirus (COVID-19): 12 March 2020 - GOV.UK [Internet]. [cited 2020 Dec 3]. Available from: https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-12-march-2020
  20. 20.
    GOV.UK. Prime Minister’s statement on coronavirus (COVID-19): 22 March 2020 - GOV.UK [Internet]. [cited 2020 Dec 3]. Available from: https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-22-march-2020
  21. 21.
    GOV.UK. Prime Minister’s statement on coronavirus (COVID-19): 25 March 2020 - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-25-march-2020
  22. 22.
    GOV.UK. Prime Minister’s statement on coronavirus (COVID-19): 11 May 2020 - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/speeches/pm-statement-on-coronavirus-11-may-2020
  23. 23.
    GOV.UK. Prime Minister sets out timeline for retail to reopen in June - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/news/prime-minister-sets-out-timeline-for-retail-to-reopen-in-june
  24. 24.
    GOV.UK. Pubs, restaurants and hairdressers to reopen from 4 July - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/news/pubs-restaurants-and-hairdressers-to-reopen-from-4-july
  25. 25.
    GOV.UK. Eat Out to Help Out launches today –with government paying half on restaurant bills - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/news/eat-out-to-help-out-launches-today-with-government-paying-half-on-restaurant-bills
  26. 26.
    GOV.UK. Schools and colleges to reopen in full in September - GOV.UK [Internet]. [cited 2020 Dec 16]. Available from: https://www.gov.uk/government/news/schools-and-colleges-to-reopen-in-full-in-september
  27. 27.
    GOV.UK. Rule of six comes into effect to tackle coronavirus - GOV.UK [Internet]. [cited 2020 Dec 7]. Available from: https://www.gov.uk/government/news/rule-of-six-comes-into-effect-to-tackle-coronavirus
  28. 28.
    GOV.UK. Prime Minister announces new local COVID Alert Levels - GOV.UK [Internet]. [cited 2020 Dec 7]. Available from: https://www.gov.uk/government/news/prime-minister-announces-new-local-covid-alert-levels
  29. 29.
    GOV.UK. Prime Minister announces new national restrictions - GOV.UK [Internet]. [cited 2020 Dec 7]. Available from: https://www.gov.uk/government/news/prime-minister-announces-new-national-restrictions
  30. 30.↵
    Funk S. socialmixr @ github.com.
  31. 31.↵
    The RECOVERY Collaborative Group. Dexamethasone in Hospitalized Patients with Covid-19 — Preliminary Report. N Engl J Med. 2020;1–11.
  32. 32.↵
    Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28(4):365–82.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    Gillespie DT. Approximate accelerated stochastic simulation of chemically reacting systems. J Chem Phys. 2001;115(4):1716–33.
    OpenUrlCrossRefWeb of Science
  34. 34.↵
    Department of Health and Social Care. COVID-19 testing data: methodology note. www.gov.uk. 2020.
  35. 35.↵
    Martin IMC, Ison CA, Aanensen DM, Fenton KA, Spratt BG. Rapid Sequence-Based Identification of Gonococcal Transmission Clusters in a Large Metropolitan Area. J Infect Dis. 2004;189(8):1497–505.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    Del Moral P, Doucet A, Jasra A. Sequential Monte Carlo samplers. J R Stat Soc Ser B Stat Methodol. 2006;68(3):411–36.
    OpenUrlCrossRef
  37. 37.↵
    Andrieu C, Doucet A, Holenstein R. Particle Markov chain Monte Carlo methods. J R Stat Soc Ser B Stat Methodol. 2010;72(3):269–342.
    OpenUrlCrossRef
  38. 38.↵
    Gordon NJ, Salmond DJ, Smith AFM. Novel approach to nonlinear/non-gaussian Bayesian state estimation. IEE Proceedings, Part F Radar Signal Process. 1993;140(2):107–13.
    OpenUrl
  39. 39.↵
    Baguelin M, Knock E, Whittles LK, FitzJohn R, Lees J. sircovid. 2020.
  40. 40.↵
    Knock ES, Whittles LK, Perez-Guzman PN, Bhatia S, Guntoro F, Watson OJ, et al. Reproducible parallel inference and simulation of stochastic state space models using odin, dust, and mcstate. Wellcome Open Res. 2020 Dec 11;5:288.
    OpenUrl
  41. 41.↵
    Verity R, FitzJohn R. mrc-ide/markovid at version1.5 [Internet]. [cited 2020 Dec 4]. Available from: https://github.com/mrc-ide/markovid/tree/version1.5
  42. 42.↵
    Brazeau NF, Verity R, Jenks S, Fu H, Whittaker C, Winskill P, et al. Infection Fatality Ratio: Estimates from Seroprevalence.
  43. 43.↵
    Jarvis CI, Van Zandvoort K, Gimma A, Prem K, Auzenbergs M, O’Reilly K, et al. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Med. 2020;18(1):1–10.
    OpenUrlCrossRef
  44. 44.
    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;
  45. 45.
    NHS England and NHS Improvement. Statistics » COVID-19 Hospital Activity [Internet]. [cited 2020 Dec 17]. Available from: https://www.england.nhs.uk/statistics/statistical-work-areas/covid-19-hospital-activity/
Back to top
PreviousNext
Posted January 13, 2021.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
The 2020 SARS-CoV-2 epidemic in England: key epidemiological drivers and impact of interventions
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
The 2020 SARS-CoV-2 epidemic in England: key epidemiological drivers and impact of interventions
Edward S. Knock, Lilith K. Whittles, John A. Lees, Pablo N. Perez-Guzman, Robert Verity, Richard G. FitzJohn, Katy AM Gaythorpe, Natsuko Imai, Wes Hinsley, Lucy C. Okell, Alicia Rosello, Nikolas Kantas, Caroline E. Walters, Sangeeta Bhatia, Oliver J Watson, Charlie Whittaker, Lorenzo Cattarino, Adhiratha Boonyasiri, Bimandra A. Djaafara, Keith Fraser, Han Fu, Haowei Wang, Xiaoyue Xi, Christl A. Donnelly, Elita Jauneikaite, Daniel J. Laydon, Peter J White, Azra C. Ghani, Neil M. Ferguson, Anne Cori, Marc Baguelin
medRxiv 2021.01.11.21249564; doi: https://doi.org/10.1101/2021.01.11.21249564
Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Citation Tools
The 2020 SARS-CoV-2 epidemic in England: key epidemiological drivers and impact of interventions
Edward S. Knock, Lilith K. Whittles, John A. Lees, Pablo N. Perez-Guzman, Robert Verity, Richard G. FitzJohn, Katy AM Gaythorpe, Natsuko Imai, Wes Hinsley, Lucy C. Okell, Alicia Rosello, Nikolas Kantas, Caroline E. Walters, Sangeeta Bhatia, Oliver J Watson, Charlie Whittaker, Lorenzo Cattarino, Adhiratha Boonyasiri, Bimandra A. Djaafara, Keith Fraser, Han Fu, Haowei Wang, Xiaoyue Xi, Christl A. Donnelly, Elita Jauneikaite, Daniel J. Laydon, Peter J White, Azra C. Ghani, Neil M. Ferguson, Anne Cori, Marc Baguelin
medRxiv 2021.01.11.21249564; doi: https://doi.org/10.1101/2021.01.11.21249564

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (70)
  • Allergy and Immunology (172)
  • Anesthesia (51)
  • Cardiovascular Medicine (456)
  • Dentistry and Oral Medicine (83)
  • Dermatology (55)
  • Emergency Medicine (160)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (192)
  • Epidemiology (5317)
  • Forensic Medicine (3)
  • Gastroenterology (199)
  • Genetic and Genomic Medicine (767)
  • Geriatric Medicine (81)
  • Health Economics (216)
  • Health Informatics (704)
  • Health Policy (364)
  • Health Systems and Quality Improvement (224)
  • Hematology (100)
  • HIV/AIDS (166)
  • Infectious Diseases (except HIV/AIDS) (5971)
  • Intensive Care and Critical Care Medicine (368)
  • Medical Education (107)
  • Medical Ethics (25)
  • Nephrology (83)
  • Neurology (781)
  • Nursing (43)
  • Nutrition (135)
  • Obstetrics and Gynecology (147)
  • Occupational and Environmental Health (237)
  • Oncology (483)
  • Ophthalmology (155)
  • Orthopedics (40)
  • Otolaryngology (98)
  • Pain Medicine (39)
  • Palliative Medicine (20)
  • Pathology (141)
  • Pediatrics (224)
  • Pharmacology and Therapeutics (138)
  • Primary Care Research (99)
  • Psychiatry and Clinical Psychology (868)
  • Public and Global Health (2051)
  • Radiology and Imaging (356)
  • Rehabilitation Medicine and Physical Therapy (160)
  • Respiratory Medicine (288)
  • Rheumatology (96)
  • Sexual and Reproductive Health (75)
  • Sports Medicine (77)
  • Surgery (110)
  • Toxicology (25)
  • Transplantation (30)
  • Urology (40)