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

Real-time Nowcasting and Forecasting of COVID-19 Dynamics in England: the first wave?

View ORCID ProfilePaul Birrell, Joshua Blake, View ORCID ProfileEdwin van Leeuwen, PHE Joint Modelling Cell, View ORCID ProfileNick Gent, View ORCID ProfileDaniela De Angelis
doi: https://doi.org/10.1101/2020.08.24.20180737
Paul Birrell
1National Infection Service, Public Health England, 61 Colindale Avenue, NW9 5HT, UK
2MRC Biostatistics Unit, University of Cambridge Institute of Public Health, Cambridge Biomedical Campus, CB2 0SR, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Paul Birrell
  • For correspondence: paul.birrell{at}phe.gov.uk
Joshua Blake
2MRC Biostatistics Unit, University of Cambridge Institute of Public Health, Cambridge Biomedical Campus, CB2 0SR, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edwin van Leeuwen
1National Infection Service, Public Health England, 61 Colindale Avenue, NW9 5HT, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Edwin van Leeuwen
Nick Gent
3Emergency Response Department, Public Health England, Porton Down, SP4 0JG, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Nick Gent
Daniela De Angelis
2MRC Biostatistics Unit, University of Cambridge Institute of Public Health, Cambridge Biomedical Campus, CB2 0SR, UK
1National Infection Service, Public Health England, 61 Colindale Avenue, NW9 5HT, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniela De Angelis
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

England has been heavily affected by the SARS-CoV-2 pandemic, with severe ‘lock-down’ mitigation measures now gradually being lifted. The real-time pandemic monitoring presented here has contributed to the evidence informing this pandemic management. Estimates on the 10th May showed lock-down had reduced transmission by 75%, the reproduction number falling from 2.6 to 0.61. This regionally-varying impact was largest in London of 81% (95% CrI: 77%–84%). Reproduction numbers have since slowly increased, and on 19th June the probability that the epidemic is growing was greater than 5% in two regions, South West and London. An estimated 8% of the population had been infected, with a higher proportion in London (17%). The infection-to-fatality ratio is 1.1% (0.9%–1.4%) overall but 17% (14%–22%) among the over-75s. This ongoing work will be key to quantifying any widespread resurgence should accrued immunity and effective contact tracing be insufficient to preclude a second wave.

1 Introduction

Since the 31st of December 2019, when the government of Wuhan reported treating cases of pneumonia of unknown cause (1), more than 13 million individuals have been reported as being infected by SARS-COV-2 globally and over 500,000 have died (2). Spread of SARS-COV-2 from Wuhan to other Chinese provinces, Thailand, Japan and the Republic of Korea rapidly occurred in the first 3 weeks of January 2020 (3), with an early incursion into Europe centred around clusters in Bavaria, Germany and Haute-Savoie, France, both linked to subsequent cases in Spain (4). This international spread eventually led the World Health Organisation to declare a pandemic on 11th March, 2020 (5). In the UK, pandemic preparedness plans, developed since 2009 A/H1N1pdm, were rapidly activated with governmental emergency bodies and advisory groups convening before the end of January (6).

The UK Response to the COVID-19 pandemic escalated from an initial containment effort to the suppression or ‘lock-down’ strategy introduced on the 23rd of March (7). Over this period, through participation in governmental advisory groups, scientists from a number of research institutions fed into the pandemic decision-making processes. Quantitative understanding of the pandemic was based on the work of SPI-M (Scientific Pandemic Influenza Sub-Group on Modelling) (8). This work has informed the various phases of the UK response, from constructing planning scenarios for the health system, to monitoring these scenarios through regular now-casting of new infections and forecasting of severe disease and health service demand, see (9,10) and papers available from (8).

Here we report the work of one of the participating groups, the Public Health England (PHE)/University of Cambridge modelling group. This collaboration, funded to develop modelling methodology for real-time pandemic influenza monitoring (11), was re-activated for the COVID-19 pandemic (12). The age and spatially structured transmission model developed for influenza (13,14), was re-purposed and, as data have increasingly become available, continuously adapted to SARS-COV-2 epidemiology, see Supplementary Materials (SM) ‘Transmission Model’.

The model is implemented through a Bayesian statistical analysis of pandemic surveillance data, incorporating information on the natural history of infection from emerging literature. The surveillance data used are age- and region-specific counts of deaths of people with a lab-confirmed CoVID-19 diagnosis (see Fig. 1A-B); and, from 21st April onwards, weekly batches of serological data, indicating the fraction of the population carrying CoVID-19 antibodies, from NHS Blood and Transplant (NHSBT) samples (see Fig. 1C) (15). The Wuhan outbreak additionally informs estimation of model parameters: the duration of infectiousness, the mean time from infection to symptom onset (16); and mean time from symptoms onset to death (17). Central to age-specific epidemic modelling, contact patterns between age groups have been derived from the POLYMOD study (18) stratified by setting (school, workplace, leisure etc), with these matrices sequentially updated using the Google mobility study and the UK time-use survey (19), to quantify the change in population mobility and access to these contact settings over time.

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

(A) Data on deaths by region, and (B) age; (C) Serological positivity by region and sampling date.

This transmission model has been used consistently throughout the pandemic. Through this regular monitoring we have been able to: anticipate and understand the impact of the lock-down; provide sequential updates of the pandemic dynamics, by estimating the basic (R0) and effective (Rt) reproduction numbers, (i.e. the average number of individuals infected by an infectious individual in: a totally susceptible population with no control measures; and a partially susceptible population under some control measures respectively); and inform the gradual relaxation of the lock-down.

Here we focus on the contribution of our work at specific dates spanning three key periods: pre-lock-down (before 23rd of March), lock-down (until 11th May); and two subsequent dates that allow the assessment of the gradual easing of the lockdown.

Pre-lock-down

After initial attempts to contain the pandemic through trace and test strategies (20) and to mitigate the burden on the National Health Service (NHS) through combinations of non-pharmaceutical interventions (e.g. case isolation, restrictions on foreign travel, shielding of vulnerable groups, cancellation of mass gatherings), the pressing question became: what level of stringent social distancing measures would be necessary to suppress transmission? At this stage, infection was not sufficiently widespread in each of the seven NHS regions for the data to inform a fully stratified model, by this time there had only been four deaths in the whole of the South-West. Therefore, no age structure was used and the country was stratified into two regions: London, where the number of deaths was significantly higher (Fig 1A), and Outside London. Assuming a pandemic intervention is imposed on 23rd March, the model was fitted to data on CoVID-19 confirmed deaths to 15th March, and then used to project epidemic curves forward a further eight weeks. These projections assumed differing reductions in contact rates (24%, 48%, 64%) (see Equation (6), SM) and consequently transmission.

Fig. 2 shows the projections for different levels of this reduction. The dashed red vertical line shows the date of the most recent data included in the analysis and the dashed purple line represents the timing of the intervention. Each column shows the projected epidemic infection and death rates under three assumed intervention impacts.

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

Estimated and projected COVID-19 infections and deaths by efficacy of social restric-tion measures.

The most optimistic scenario in Fig. 2 corresponded to an immediate 64% reduction in trans-mission. Under this assumption, Rt was estimated to be 1.2 (95% CrI: 0.83–2.1) in London and 1.2 (95% CrI: 0.84–2.2) elsewhere. In this scenario, the probability that the imposed measures were successful in reducing Rt to the threshold of 1 required for declining transmission was only 19% and 17% in the two regions respectively. To be 95% certain that the intervention would lead to a sustained decline in infection, the intervention would need to induce an 81% reduction in transmission. Such a reduction could only be achieved through the implementation of extreme mitigation measures.

The lock-down period

The number of deaths continued to rise until April 8th, particularly in older age-groups (Fig. 1(B)), permitting stratification of the model by both age (8 groups) and region (7 NHS regions). Also new information from serological studies (Fig. 1(C)) started to become available and weekly data downloaded from the Google mobility survey could be used to update contact matrices.

A rhythm for pandemic monitoring was established. The model was run daily, with results feeding into local planning tools as well as the SPI-M consensus view on the state of the pandemic, and with periodic publication of web-reports summarising the latest results (21). These outputs included a number of key indicators: regional estimates of Rt and epidemic growth rates r, indicating whether transmission is increasing (Rt > 1) and the rate at which it is increasing (22); region and age-specific attack rates (i.e. the proportion of the population already infected); and predictions of the burden due to mortality, both in terms of age-specific infection-fatality ratios and number of CoVID19 deaths. Public attention has focused on Rt as a headline figure for the state of the pandemic, but a more complete assessment requires all these indicators. Table 1 presents estimates of a selection of these indicators, giving snapshots of the pandemic state at the three chosen times.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1:

Table of estimates (with 95% credible intervals attached) for key epidemic parameters and derived quantities.

The 10th May section of the table shows the success of the lock-down at curtailing transmission: the Rt in England is now estimated to be 0.75 (95% CrI: 0.72–0.77) having dropped from 2.6 (95% CrI: 2.4–2.9) to 0.61 (95% CrI: 0.57–0.67) at the time of the lock-down, a reduction of 75% (95% CrI: 73%–77%), in line with the anticipation of the pre-lock-down modelling. The ‘growth’ rate for England indicates the daily number of infections were halving every log(2)/r = 11.5 days. London stands out as the region with the highest estimated attack rate (20% of people infected, CrI: 16%–26%), largely due to pre-lock-down levels of infection; the largest reduction in transmission, a drop of 81% (95% CrI: 77%–84%) to Rt = 0.4, 95% CrI: 0.36–0.43, and the steepest rates of decline in both the number of infections (halving every 4 days) and the observed fall in the number of deaths (Fig. 1A). The temporal patterns in infection are disrupted at the lock-down date, with the top row of Fig. 3(A) illustrating the size of this effect for both the London and North West Regions, alongside the estimated Rt.

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

(A) Sequentially obtained estimates infection trajectories; Rt; and three week forecast of the number of deaths by date of analysis (rows). (B) Probability that Rt exceeds 1 as of 3rd June, by region. See SM, Fig. S1 for 10th May8and 19th June. In both plots, vertical blue lines indicate the timings of significant policy changes, vertical red lines indicate the time of analysis.

Launched on the 10th May, the UK COVID 4-level Alert System, co-ordinated by the newly established Joint Biosecurity Centre (JBC) and based on both estimated Rt and current infec-tions, highlighted the key role of these indicators (7) in guiding the relaxation of lock-down measures without re-igniting transmission.

Lock-down Relaxation

The first tranche of relaxations were announced on the 10th May. From this time, accounting for likely changes in behaviour became crucial. In addition, it was clear, from the over-precision of the estimates of both incidence and Rt in the top row of

Fig. 3(A), that the model required a greater flexibility to capture such changing behaviours. From an appropriately adapted model (see SM ‘Transmission model’ for details), we estimated that at the 3rd June, Rt for England reached 0.99 (95% CrI: 0.91–1.09), with the probability of exceeding the value of 1 rapidly increasing over time (see Fig. 3(B)). This figure masks regional heterogeneity in transmission. The North West and the South West were characterised by Rt values above 1 (with a probability Rt > 1 of over 50%, Fig. 3(B)) and growth rate estimates encompassing positive values. For the North West, we estimated 4,170 (95% CrI:1,580–9,840) daily infections, the highest number in the country (see Table 1). The step changes in the plots of Rt over time in Fig. 3 for 10th May are entirely due to changes in the Google mobility matrices. Looking at the equivalent plot for the 3rd June analysis, over the same interval, the step changes are larger. This difference in the level of fluctuation over time suggests that the increases in Rt are too large to be solely attributed to mobility-driven changes in the contact matrices. Furthermore, for the North West, the drop in Rt around the lock-down is not as sharp as in London, but rather staggered over three weeks. This might suggest a different response to the lock-down in these two regions, which we had previously not been able to identify. The estimated steady resurgence of Rt in the North West ultimately led to a policy change, delaying the staged re-opening of schools (23).

Continuing to monitor the pandemic evolution in the post-lock-down era, we adapted the model to incorporate new evidence on differential susceptibility to infection by age (24) (see SM, around equation (7)). Results from 19th June data (Table 1, Fig. 3(A)) show lower estimates for Rt, negative growth rates and the estimated number of infections in England decreasing to 4,300 (95% CrI: 2,400 - 7,300). There is still regional heterogeneity, with two regions for which the credible intervals for Rt exclude 1 (North East & Yorkshire, and the South East); and the probability that the epidemic is growing is 30% in the South West and below 15% in each of the other regions (Fig S1(B), SM). Throughout we have been estimating age-specific infection-fatality ratios (see Tables 1 and S4). Allowing for differential age susceptibility, the age-specific estimates of the infection-fatality ratio fall to 17% (95% CrI: 14%–22%) in the over-75s (from 23%, 95% CrI: 20%–27%) with a rise to 2.9% in the 65–74ys, (see Table S4) and to 1.1% (95% CrI: 0.9%–1.4%) from 0.9% (95% CrI: 0.8%–1.0%) overall. These less severe estimates (in comparison to the 3rd June analysis) led to the UK Chief Medical Officers agreeing with a JBC recommendation that the alert level should be downgraded to level three.

Discussion

Since the 19th of June, weekly updates have continued to be posted online (21). Recent estimates of the IFR have varied in the region of 0.9% to 1.4% overall and 14% to 19% in the over-75s. Rt is approaching the value of 1 in most regions without exceeding it, which, together with a decreasing number of daily infections indicates an epidemic still in decline, although local outbreaks are being increasingly detected (25). The estimates presented here are consistent with the SPI-M consensus on the values of Rt both nationally and regionally (26). Incidence estimates can be contrasted with estimates from community cross-sectional studies. In a report of the 9th July, the Office for National Statistics (ONS) estimates 1,700 (range 700–3,700) new daily infections over the two weeks leading up to 4th July (27), while the ZOE study (28) reports an average 1,470 infections per day over a similar period. In our work incidence ranges from a central estimate of 4,200 daily infections down to 3,500 over the period. These are not incongruous to the ZOE app estimates, which only records symptomatic infection and require some scaling to derive an estimate for total infections. The ONS estimate is likely to be an under-estimate as the survey does not include individuals in institutionalised settings e.g. care homes, where incidence may be far higher than in the community. The degree to which this is an under-estimate, however, is unclear.

The ONS also reports (29) that 6.3% (95% CrI: 4.7%–8.1%) of individuals showed the presence of antibodies to the COVID-19 in blood sera samples (as of 19th June), a little lower than the 8% estimated attack rate for England. This discrepancy may well be due to the timing of the samples, due to the waning of the antibody response over time. Estimates of attack rates (see Table 1) show that our belief on the proportion of the population that has been infected is being revised downwards at each sequential analysis. This, together with the emerging evidence of waning immunity (30), paint a muddied picture in term of potential for a population-level resurgence in infection.

On the 4th July there was a significant step change in the gradual relaxation of pandemic mitigation measures as leisure facilities, tourist attractions, pubs and cafes all became accessible to the public once again (31). The impact of this is still largely to be observed and assessed. If the concerning worse-case scenario presented for the Winter in the recent publication of the Academy of Medical Science (32) comes through, we will be facing a very challenging second wave. Our monitoring tool, continuously adapted to incorporate the accumulating surveillance data, will be key to providing quantification of all aspects of any resurgence.

Data Availability

The data used in the manuscript are routinely collected surveillance data on the ongoing CoVID-19 pandemic. We are in the process of agreeing and appropriate level of aggregation under which these data can be made available in a public repository.

Supplementary Materials

Methods

Surveillance data

The UK confirmed its first case of CoVID infection, imported from China, on the 30th January, with the first CoVID-linked death reported on the 6th March. Initially deaths were monitored by date of report, (33) with a full line-listing of deaths being made available by PHE from the 21st March, enabling tracking of the pandemic using actual dates of death.

Prior to this time, the most complete data were the time series of reported cases, individuals who have a returned a PCR-positive swab sample. (33) Initial analyses of the pandemic were based on either or both of these datasets, attempting to quantify the emerging pandemic but also to monitor the degree of coherence between the datasets. However, rates of testing and the country’s testing capacity were significantly increasing, (34), whilst government policy was changing the testing threshold. These led to a significant ascertainment bias to the confirmed case data which in turn would lead to estimated incidence curves based on these data to be significantly skewed.

A third data source was the testing of residual blood sera from blood donors in England via National Health Service Blood and Transplant (NHSBT). These data inform the population prevalence of antibodies, informing the fraction of the population that have ever been infected. From repeated follow-up of known infections, 74.7% (35) of infections return a positive result for the presence of antibodies within 21 days of the onset of symptoms

Transmission Model

An age stratified transmission model is fitted to each region under study simultaneously with the sharing of some global parameters. Within each region, the transmission dynamics are governed by a system of ordinary differential equations, discretised to give the following set of first order difference equations: Embedded Image where: Embedded Image, Embedded Image, Embedded Image, d =1, 2 represent the time tk,k = 1,…, K, partitioning of the population of individuals in a region r,r = 1, …, nr, in age-group i,i = 1,…, nA, into S (susceptible), E (exposed) and I (infectious) disease states. The mean latent and infectious periods are dL and dI respectively; and Embedded Image is the time- and age-varying rate with which sus-ceptible individuals become infected. Time steps of δt = 0.5 days are chosen to be sufficiently small relative to the anticipated latent and infectious periods. New infections are generated as Embedded Image where Embedded Image

Here, Embedded Image is the probability of a susceptible individual in region r of age group i being infected by an infectious individual in age group j at time tk. It is a function of:

  • a set of time-varying contact matrices Embedded Image, the entries of which describe the expected number of contacts between any two individuals of the different strata within a single time unit.

  • Embedded Image, a region-specific matrix, the (i, j)th element of which gives the relative susceptibility of someone in age-group i to an infection from an infectious individual in age-group j assuming contact between the two.

  • βtk,r, a time-varying parameter encapsulating further temporal fluctuation in transmission that applies to all ages.

  • R0,r are the initial reproduction numbers for the pandemic in each region at time t0

  • Embedded Image are the dominant eigenvalues of the initial next-generation matrices, Λ0,r such that Embedded Image where Nr,t is the population size in region r and age-group i, and crk are a set of matrices defined by Embedded Image where the ⊙ notation indicates element-wise multiplication, s.t A = B ⊙ C if Aij = BijCij.

As described in the following sub-section, the Ctk matrices encode the information con-tained about contact rates between different age groups derived from the POLYMOD study (18), Google mobility and the time-use survey. The Embedded Image matrices capture any misspecification of these matrices in terms of the changing pattern of infection between the age groups, whereas the Embedded Image parameters account for misspecification of the changing scale of transmission over time as described by the matrices.

In the main manuscript, a number of model developments are identified that were imple-mented in between the analyses presented. The general expression of Embedded Image is: Embedded Image but each of the models implemented constitute some simplification of this expression. These model developments are itemised here, starting with the earliest, and most simple, analysis:

  • 20th March: No age dimension to the analysis, so transmission was simply broken up with a changepoint at tk = tlock. The dynamics of the two regions (London, Outside London) differed only in that they had a differing initial seeding of infection. Here Embedded Image and Embedded Image, giving m the interpretation as a global change in the susceptibility to infection given contact with an infectious individual. From (4), Embedded Image. Simplifying (5): Embedded Image

  • 10th May: Analysis is now age-specific. Reproduction number R0,r and susceptibility parameter mr now region-specific. The matrices Embedded Image are a matrix extension of the 20th March analysis: Embedded Image where the susceptibility parameter mr now has region dependence, Embedded Image is a vector of ones of length nA, and Embedded Image is a corresponding matrix of ones. This has no age dependence, and so it is assumed that the POLYMOD and Google-derived matrices, Embedded Image, accurately account for changing patterns of infection over time). Again (4) can be simplified to: Embedded Image

  • 3rd June: Introduction of the time-varying transmissibility parameter, Embedded Image It was decided to change these piecewise with weekly changepoints. Denote wk = w(tk) giving the week in which time tk falls. Over time these are not allowed to vary unconstrained and a smoothing is imposed by assuming, a priori that they develop according to a random-walk. Embedded Image

    Where Embedded Image is the rate that applies in all weeks up to and including the first week of the lock-down. As all other components of the model remain unchanged, the form for Embedded Image is: Embedded Image

  • 19th June: Separate levels of susceptibility are introduced for the over-75s, both before and after the lock-down to account for some lack of fit. To characterise this, we need to define two matrices, Embedded Image and Embedded Image that apply before and after tlock respectively, Embedded Image such that> Embedded Image

    We would then feed back into equation (3): Embedded Image

Calculating Rt

We estimate a region-specific parameter describing the initial exponential growth in infections, denoted ψr. Using the formula of Wearing et al. (2005), this is related to the initial reproduction number R0, r via the formula

Embedded Image

Over time the value of the reproduction number will change as contact patterns shift and the supply of susceptible individuals deplete. The formula for the time-t reproduction number is Embedded Image where Embedded Image is the dominant eigenvalue of the time tk next-generation matrix, Λk,r, with elements:

Embedded Image

To get an ‘all England’ value for Embedded Image a weighted average of the regional Embedded Image is calculated, where the weights are given by the sum of the infections in each region:

Embedded Image

Contact matrices

Time dependent contact matrices are based on location-specific POLYMOD matrices (where locations include “at work”, “at home”, “on transport” etc), combined with the time-use survey (19) to identify 18 different activities, including school, work, social visits, shopping etc.. The traditional POLYMOD matrices are used until 23rd March, the time of the lock-down (18, 36). From this point on, the Google mobility and time-use survey data were used to calculate proportionate reductions in the location-specific POLYMOD matrices, which are summed together to give a weekly-varying contact matrix, Embedded Image.

Following (19) we identified 18 different activities (Table S1). Some of these activities were not allowed during lock down and, therefore, assumed to have been stopped. For other activities we used the most relevant data source. For some activities no suitable data source was available. In that case we used the retail and recreation mobility data provided by Google, because this data was assumed to best represent the general adherence level in the UK. For example, visits are unlikely to have stopped completely. Instead the retail and recreation level is used (Table S1). Note that the activities which in a well mixed model have the most effect on the base reproductive number are school, work, visits and unspecified (19).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S1:

The data sources used for the different activities identified by the time use survey. Values represent the assumed activity level. Where direct data was available we used other (live) data sources. Here GM represents the relevant Google mobility category.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S2:

The fraction of contacts at home attributed to members of the household.

School attendance

School size and age range is publicly available for England. We used this to calculate the number of students for each school in the modelled age groups, based on the assumption that the students were evenly distributed across different school years. We also had access to attendance levels over time for some schools. This data was then used to calculate the average weekly attendance level in each local authority, weighted by the size of the school in that age group. Finally, we combined these values to calculate attendance in England, weighted by local authority population size. Not all schools reported attendance every week, as a result the attendance by local authority was, in rare cases, only based on the report of 1 or more small schools. To ensure that these did not skew our results, we ignored any attendance estimate based on schools with less than 100 students in total.

Visits at home

(19) used the time usage data for visits at home to estimate the number of visit related contacts at home. We found that this underestimates the number of visitors and instead used the POLYMOD data on household size to estimate the fraction of contacts at home with other household members versus the contacts with others (e.g. visitors). First, we extracted the mean number of contacts at home (Table S2). Next, we limited numbers of contacts at home (ci) to household size (hi) minus one Embedded Image, i.e. the maximum number of contacts any participant (i) can have with just household members and calculated the mean contacts based on that value. Note that this provides a conservative estimate of the fraction of contacts from visits, because some participants will not have met all their household members during the day.

Google mobility data

During the pandemic Google provided aggregated mobility data from Android phones for many countries (https://www.google.com/covid19/mobility/). The mobility data gives an indication of the activity level for 5 different activities: retail and recreation, grocery and pharmacy, parks, transit stations, workplaces and residential. The UK data was further subdivided into activity by local authority, by matching the Google provided locations to local authority districts in England. This data was then combined, weighted by population size. Finally, the daily values were averaged by week, to produce weekly activity levels.

Likelihood

Counts of deaths are assumed to have a negative binomial distribution. If Embedded Image is the number of observed deaths on day tk in age group i in region r, then the expected number of deaths are derived using (2), an assumed-known distribution of the time from infection to death from COVID-19, f and an estimated age-specific infection-fatality ratio pi: Embedded Image where fl gives the probability the death occurs on the lth day after infection. Then we assume Embedded Image where η is a dispersion parameter such that Embedded Image and Embedded Image.

The serological data is dependent on two parameters, the sensitivity and the specificity of the serological testing process, ksens and kspec respectively. If, on day tk, Embedded Image blood samples are taken from individuals in region r and age-group i, and the observed number of positive tests is Embedded Image, then

Embedded Image

Inference

Table S3 contains a full list of parameters. The top portion refers to the unknown parameters that are estimated and the bottom portion contain parameters that have been fixed. Wherever possible, a brief justification of the chosen prior/value is given.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S3:

Model parameters with assumed prior distributions or fixed values Name Prior source

Estimation is conducted in a Bayesian framework by combining the prior information and likelihoods specified above, to derive posterior distributions for the unknown parameters and any functions of such parameters. Sampling from these posterior distribution is carried out through Markov chain Monte Carlo (13,14). The most recent analysis featured in this paper was based on 900,000 iterations, with an initial adaptive phase of 45,000 iterations within a burn-in period of 90,000 iterations. Parameter estimates are based on the full sample following burn-in thinned to retain every 125 iterations, and projections are based on a sample thinned to every 250 iterations. All central estimates are pointwise medians of quantities calculated on the basis of this sample, and uncertainty is expressed through 95% credible intervals (CrI) derived from the 2.5% and 97.5% quantiles.

Additional Results

In Fig. S1, we present the region-specific probabilities that Rt > 1. In the 10th May analysis, the plot is binary - Rt is definitely above the threshold before the lockdown, and we estimated the reduction in transmission to be sufficient to reduce it to below the threshold with near certainty. However, by the time of the 19th June analysis, there is much greater flexibility and uncertainty in the evolution of Rt over time and correspondingly, after 3–4 weeks beyond the lock-down, Rt can no longer be said with certainty to be below 1. From this time there is an initial resurgence in the probability Rt > 1 which then plateaus from mid-May onwards. By the end of the period, there is the possibility that Rt > 1 in both the South West and in London, but it remains unlikely.

Table S4 gives the full estimation of the infection-fatality ratios only partially presented in Table 1.

Fig. S2 shows the goodness-of-fit to the death data by age, for each of the three post-lockdown analyses considered. In the 3rd June analysis, there is a consistent under-estimation of the number of deaths from the beginning of May in the over-75s, with corresponding over-estimation in the younger age-groups. The model adaptation to allow for a different susceptibility to infection for the over-75s appears to rectify this lack of fit.

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

Equivalent plots to Fig. 3(B) for the analysis on (A) 10th May, and (B) 19th June.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S4:

Estimates for the age-specific fatality ratio as they change over time

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

Goodness-of-fit to the deaths data by age on the 10th May, 3th June and 19th June

Acknowledgements

This work was supported by the Medical Research Council (Unit programme number MC UU 00002/11) in partnership with Public Health England. Prior to the pandemic, the project was de-veloped under a grant from the National Institute for Health Research (HTA Project: 11/46/03). We gratefully acknowledge the access to the data from the United Kingdom Time Use Survey through the UK Data Service (http://doi.org/10.5255/UKDA-SN-8128-1). We ac-knowledge the support of the PHE Epidemiology Cell in consistently providing the data streams used. We would also like to thank Dr Nigel Gay for the initial development of the model and constructive comments on the manuscript. This work has also benefited from insightful discus-sions with Dr Petra Klepac.

The core membership of the PHE Modelling Cell consists of: Abbygail Jaccard, Andre Charlett, Anna Rance, Anne Presanis, Archana Purohit, Brian Ferguson, Brodie Walker, David Mustard, Declan Bays, Dianne Addei, Emilia Vynnycky, Emily Agnew, Emma Bennett, Emma Gillingham, Hannah Williams, Ian Hall, James Lewis, Jonathan Carruthers, Joseph Shingleton, Joshua Blake, Judith Field, Martin Grunnill, Matt Edmunds, Matt Hennessey, Nick Gent, Peter White, Simona Baracaia, Stephanie Shadwell, Steven Dyke, Thomas Finnie, Virginia Cox and Xu-Sheng Zhang.

References

  1. 1.↵
    World Health Organization, Pneumonia of unknown cause - China, https://www.who.int/csr/don/05-january-2020-pneumonia-of-unkown-cause-china/en/ (2020). WHO, Accessed 25 June, 2020.
  2. 2.↵
    Center for Systems Science and Engineering (CSSE), Johns Hopkins University, COVID-19 Case Tracker, https://coronavirus.jhu.edu/map.html (2020). Johns Hopkins University, Accessed 25 June, 2020.
  3. 3.↵
    J. M. Read, J. R. Bridgen, D. A. Cummings, A. Ho, C. P. Jewell, medRxiv p. 2020.01.23.20018549 (2020).
  4. 4.↵
    G. Spiteri, et al., Eurosurveillance 25, 2000178 (2020).
    OpenUrl
  5. 5.↵
    World Health Organization, WHO Director-General’s opening remarks at the media briefing on COVID-19 - 11 March 2020, https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19---11-march-2020 (2020). WHO, Accessed 25 June, 2020.
  6. 6.↵
    UK Government, Precautionary SAGE meeting on Wuhan Coronavirus (WN-CoV), https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/888767/S0369_Precautionary_SAGE_meeting_on_Wuhan_Coronavirus_WN-CoV.pdf (2020). UK Government, Accessed 9 July, 2020.
  7. 7.↵
    UK Government, Prime Minister’s statment on coronavirus (COVID-19): 23 March 2020, https://www.gov.uk/government/speeches/pm-address-to-the-nation-on-coronavirus-23-march-2020 (2020). UK Government, Accessed 9 July, 2020.
  8. 8.↵
    UK Government, Scientific Advisory Group for Emergences (SAGE): Coronavirus (COVID-19) response, https://www.gov.uk/government/groups/scientific-advisory-group-for-emergencies-sage-coronavirus-covid-19-response (2020). UK Government, Accessed 9 July, 2020.
  9. 9.↵
    S. Abbott, et al., Wellcome Open Research 5, 112 (2020).
    OpenUrl
  10. 10.↵
    M. J. Keeling, et al., medRxiv p. 2020.05.10.20083683 (2020).
  11. 11.↵
    C. Simpson, D. Beever, K. Challen, et al., Lancet Infect Dis 19, 295 (2019).
    OpenUrl
  12. 12.↵
    C. R. Simpson, et al., The UK hibernated pandemic influenza research portfolio: triggered for COVID-19 (2020).
  13. 13.↵
    P. J. Birrell, X.-S. Zhang, R. G. Pebody, N. J. Gay, D. De Angelis, Scientific Reports 6, 29004 (2016).
    OpenUrl
  14. 14.↵
    P. Birrell, R. Pebody, A. Charlett, X.-S. Zhang, D. De Angelis, Health Technology Assess-ment 21 (2017).
  15. 15.↵
    Public Health England Surveillance Cell, Serological Surveillance: Summary Report 3, https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/890194/s0337-serological-surveillance-summary-3-060520-sage34.pdf (2020). PHE, Accessed 12 July, 2020.
  16. 16.↵
    Q. Li, et al., New England Journal of Medicine 382, 1199 (2020).
    OpenUrlCrossRefPubMed
  17. 17.↵
    R. Verity, et al., The Lancet Infectious Diseases (2020).
  18. 18.↵
    J. Mossong, et al., PLoS Medicine 5, 0381 (2008).
    OpenUrl
  19. 19.↵
    E. van Leeuwen, {PHE Joint modelling cell}, F. Sandmann, medRxiv p. 2020.06.03.20067793 (2020).
  20. 20.↵
    N. L. Boddington, et al., medRxiv p. 2020.05.18.20086157 (2020).
  21. 21.↵
    P. J. Birrell, J. Blake, E. van Leeuwen, D. De Angelis, Nowcasting and fore-casting of COVID-19, https://www.mrc-bsu.cam.ac.uk/tackling-covid-19/nowcasting-and-forecasting-of-covid-19/ (2020). MRC-BSU, Ac-cessed 13 July, 2020.
  22. 22.↵
    Gog, Julia and Thomas, Rachel and Freiberger, Marianne, The growth rate of COVID-19, https://plus.maths.org/content/epidemic-growth-rate (2020). +plus magazine, Accessed 14 July, 2020.
  23. 23.↵
    The Guardian, https://www.theguardian.com/education/2020/jun/06/schools-north-west-england-postpone-reopening-coronavirus (2020). The Guardian, Accessed 14 July, 2020.
  24. 24.↵
    R. M. Viner, et al., medRxiv p. 2020.05.20.20108126 (2020).
  25. 25.↵
    E. Mahase, BMJ (Clinical research ed.) 369, m2635 (2020).
    OpenUrlFREE Full Text
  26. 26.↵
    UK Government, The R number and growth rate in the UK, https://www.gov.uk/guidance/the-r-number-in-the-uk (2020). UK Government, Accessed 14 July, 2020.
  27. 27.↵
    Office for National Statistics: COVID-19 Infection Survey, https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/bulletins/coronaviruscovid19infectionsurveypilot/england9july2020#incidence-rate (2020). Office for National Statistics.
  28. 28.↵
    COVID Symptom Study, https://covid.joinzoe.com/post/covid-incidence-uk (2020). COVID Symptom Study, Accessed 13 July, 2020.
  29. 29.↵
    Office for National Statistics: COVID-19 Infection Survey, https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/bulletins/coronaviruscovid19infectionsurveypilot/2july2020#antibody-data (2020). Office for National Statistics.
  30. 30.↵
    J. Seow, et al., medRxiv p. 2020.07.09.20148429 (2020).
  31. 31.↵
    Dunn, Phoebe and Allen, Lucinda and Cameron, Genevieve and Alderwick, Hugh, COVID-19 policy tracker, https://www.health.org.uk/news-and-comment/charts-and-infographics/covid-19-policy-tracker (2020). The Health Foundation, Accessed 13 July, 2020.
  32. 32.↵
    The Academy of Medical Sciences, https://acmedsci.ac.uk/file-download/51353957 (2020). The Academy of Medical Sciences, Accessed 14 July, 2020.
  33. 33.↵
    UK Government, Coronavirus (COVID-19) in the UK, https://coronavirus.data.gov.uk/ (2020). UK Government, Accessed 12 June, 2020.
  34. 34.↵
    UK Government, Number of coronavirus (COVID-19) cases and risk in the UK, https://www.gov.uk/guidance/coronavirus-covid-19-information-for-the-public (2020). UK Government, Accessed 12 June, 2020.
  35. 35.↵
    J. Duggan, Rare and Imported Pathogens Laboratory, Public Health England, Evaluation of the Euroimmun Anti-SARS-CoV-2 ELISA (IgG) serology assay for the detection of anti-SARS-CoV-2 antibodies, https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/893433/Evaluation_of_Euroimmun_SARS_CoV_2_ELISA_IgG__1_.pdf (2020). Public Health England, Accessed 20 July, 2020.
  36. 36.↵
    S. Funk, socialmixr: Social mixing matrices for infectious disease modelling, The Compre-hensive R Archive Network, http://datacompass.lshtm.ac.uk/646/ (2018).
Back to top
PreviousNext
Posted August 30, 2020.
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.
Real-time Nowcasting and Forecasting of COVID-19 Dynamics in England: the first wave?
(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
Real-time Nowcasting and Forecasting of COVID-19 Dynamics in England: the first wave?
Paul Birrell, Joshua Blake, Edwin van Leeuwen, PHE Joint Modelling Cell, Nick Gent, Daniela De Angelis
medRxiv 2020.08.24.20180737; doi: https://doi.org/10.1101/2020.08.24.20180737
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Real-time Nowcasting and Forecasting of COVID-19 Dynamics in England: the first wave?
Paul Birrell, Joshua Blake, Edwin van Leeuwen, PHE Joint Modelling Cell, Nick Gent, Daniela De Angelis
medRxiv 2020.08.24.20180737; doi: https://doi.org/10.1101/2020.08.24.20180737

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

  • Public and Global Health
Subject Areas
All Articles
  • Addiction Medicine (430)
  • Allergy and Immunology (756)
  • Anesthesia (221)
  • Cardiovascular Medicine (3294)
  • Dentistry and Oral Medicine (364)
  • Dermatology (279)
  • Emergency Medicine (479)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1171)
  • Epidemiology (13376)
  • Forensic Medicine (19)
  • Gastroenterology (899)
  • Genetic and Genomic Medicine (5153)
  • Geriatric Medicine (482)
  • Health Economics (783)
  • Health Informatics (3269)
  • Health Policy (1140)
  • Health Systems and Quality Improvement (1190)
  • Hematology (431)
  • HIV/AIDS (1017)
  • Infectious Diseases (except HIV/AIDS) (14629)
  • Intensive Care and Critical Care Medicine (913)
  • Medical Education (477)
  • Medical Ethics (127)
  • Nephrology (523)
  • Neurology (4925)
  • Nursing (262)
  • Nutrition (730)
  • Obstetrics and Gynecology (883)
  • Occupational and Environmental Health (795)
  • Oncology (2524)
  • Ophthalmology (724)
  • Orthopedics (281)
  • Otolaryngology (347)
  • Pain Medicine (323)
  • Palliative Medicine (90)
  • Pathology (543)
  • Pediatrics (1302)
  • Pharmacology and Therapeutics (550)
  • Primary Care Research (557)
  • Psychiatry and Clinical Psychology (4212)
  • Public and Global Health (7504)
  • Radiology and Imaging (1706)
  • Rehabilitation Medicine and Physical Therapy (1013)
  • Respiratory Medicine (980)
  • Rheumatology (480)
  • Sexual and Reproductive Health (497)
  • Sports Medicine (424)
  • Surgery (548)
  • Toxicology (72)
  • Transplantation (236)
  • Urology (205)