Efficacy of stay-at-home policy and transmission of COVID-19 in Toronto, Canada: a mathematical modeling study

BackgroundIn many parts of the world, restrictive non-pharmaceutical interventions (NPI) that aim to reduce contact rates, including stay-at-home orders, limitations on gatherings, and closure of public places, are being lifted, with the possibility that the epidemic resurges if alternative measures are not strong enough. Here we aim to capture the combination of use of NPIs and reopening measures which will prevent an infection rebound. MethodsWe employ an SEAIR model with household structure able to capture the stay-at-home policy (SAHP). To reflect the changes in the SAHP over the course of the epidemic, we vary the SAHP compliance rate, assuming that the time to compliance of all the people requested to stay-at-home follows a Gamma distribution. Using confirmed case data for the City of Toronto, we evaluate basic and instantaneous reproduction numbers and simulate how the average household size, the stay-at-home rate, the efficiency and duration of SAHP implementation, affect the outbreak trajectory. FindingsThe estimated basic reproduction number R_0 was 2.36 (95% CI: 2.28, 2.45) in Toronto. After the implementation of the SAHP, the contact rate outside the household fell by 39%. When people properly respect the SAHP, the outbreak can be quickly controlled, but extending its duration beyond two months (65 days) had little effect. Our findings also suggest that to avoid a large rebound of the epidemic, the average number of contacts per person per day should be kept below nine. This study suggests that fully reopening schools, offices, and other activities, is possible if the use of other NPIs is strictly adhered to. InterpretationOur model confirmed that the SAHP implemented in Toronto had a great impact in controlling the spread of COVID-19. Given the lifting of restrictive NPIs, we estimated the thresholds values of maximum number of contacts, probability of transmission and testing needed to ensure that the reopening will be safe, i.e. maintaining an Rt < 1. Research in contextO_ST_ABSEvidence before this studyC_ST_ABSA survey on published articles was made through PubMed and Google Scholar searches. The search was conducted from March 1 to August 13, 2020 and all papers published until the end of this research were considered. The following terms were used to screen articles on mathematical models: "household structure", "epidemic model", "SARS-CoV-2", "COVID-19", "household SIR epidemic", "household SIS epidemic", "household SEIR epidemic", "quarantine, isolation model", "quarantine model dynamics", "structured model isolation". Any article showing, in the title, application of epidemic models in a specific country/region or infectious diseases rather than SARS-CoV-2 were excluded. Articles in English were considered. Added value of this studyWe develop an epidemic model with household structure to study the effects of SAHP on the infection within households and transmission of COVID-19 in Toronto. The complex model provides interesting insights into the effectiveness of SAHP, if the average number of individuals in a household changes. We found that the SAHP might not be adequate if the size of households is relatively large. We also introduce a new quantity called symptomatic diagnosis completion ratio (d_c). This indicator is defined as the ratio of cumulative reported cases and the cumulative cases by episode date at time t, and it is used in the model to inform the implementation of SAHP. If cases are diagnosed at the time of symptom onset, isolation will be enforced immediately. A delay in detecting cases will lead to a delay in isolation, with subsequent increase in the transmission of the infection. Comparing different scenarios (before and after reopening phases), we were able to identify thresholds of these factors which mainly affect the spread of the infection: the number of daily tests, average number of contacts per individual, and probability of transmission of the virus. Our results show that if any of the three above mentioned factors is reduced, then the other two need to be adjusted to keep a reproduction number below 1. Lifting restrictive closures will require the average number of contacts a person has each day to be less than pre-COVID-19, and a high rate of case detection and tracing of contacts. The thresholds found will inform public health decisions on reopening. Implications of all the available evidenceOur findings provide important information for policymakers when planning the full reopening phase. Our results confirm that prompt implementation of SAHP was crucial in reducing the spread of COVID-19. Also, based on our analyses, we propose public health alternatives to consider in view of a full reopening. For example, for different post-reopening scenarios, the average number of contacts per person needs to be reduced if the symptomatic diagnosis completion ratio is low and the probability of transmission increases. Namely, if fewer tests are completed and the usage of NPIs decreases, then the epidemic can be controlled only if individuals can maintain contact with a maximum average number of 4-5 people per person per day. Different recommendations can be provided by relaxing/strengthening one of the above-mentioned factors.