SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, January 1-July 6, 2020 ================================================================================================= * Isaac Chun-Hai Fung * Yuen Wai Hung * Sylvia K. Ofori * Kamalich Muniz-Rodriguez * Po-Ying Lai * Gerardo Chowell ## Abstract **Objective** To investigate COVID-19 epidemiology in Alberta, British Columbia and Ontario, Canada. **Methods** We calculated the incidence rate ratio (January 1—July 6, 2020) between the 3 provinces, and estimated time-varying reproduction number, *Rt*, starting from March 1, using EpiEstim package in R. **Results** Using British Columbia as a reference, the incidence rate ratios in Alberta and Ontario are 3.1 and 4.3 among females, and 3.4 and 4.0 among males. In Ontario, *Rt* fluctuated ~1 in March, reached values >1 in early and mid-April, then dropped <1 in late April and early May. *Rt* rose to ~1 in mid-May and then remained <1 from late May through early July. In British Columbia, *Rt* dropped <1 in early April, but it increased towards the end of April. *Rt* <1 in May while it fluctuated around 1.0 in June and early July. In Alberta, *Rt* > 1 in March; *Rt* dropped in early April and rose again in late April. In much of May, *Rt* <1, but *Rt* increases in early June and fluctuates ~1 since mid-June. **Conclusions** *Rt* wavering around 1.0 indicated that three provinces of Canada have managed to achieve limited onward transmission of SARS-CoV-2 as of early July 2020. Keywords * Canada * Coronavirus * COVID-19 * Epidemiology * Infectious Disease * Transmission * SARS-CoV-2 ## Introduction In 2020, the pandemic of coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome virus 2 (SARS-CoV-2), spread across Canada. The first imported case was presented in a Toronto hospital on January 23. 2020.1,2 On March 13, Quebec was the first province to declare a public health emergency;3 four days later, Alberta, British Columbia and Ontario also declared public health emergencies (Table S1).4 As of July 10, a cumulative total of 107,126 cases, including 8759 deaths, have been reported in Canada.5 As seen in the Canadian government’s outbreak update,6,7 the national epidemic trajectory appears to have stabilized in early July, but spatial-temporal pattern displays substantial heterogeneity across provinces. Epidemiologists have explored various aspects of the pandemic in Canada, including data availability,8 syndromic surveillance,9 disease burden and mortality,10 as well as the epidemiology in specific settings or subpopulations, such as nursing homes,11 and intensive care units.12 Mathematical modelers have also made projections of epidemic trajectories to assist Ottawa and the provincial governments in their decision-making.7,13-16 Against the backdrop of recent literature, our study provides a unique perspective by estimating and interpreting the time-varying reproduction number, *Rt*, of SARS-CoV-2 at both provincial and subprovincial levels in 3 Canadian provinces. This study described the COVID-19 pandemic as it unfolded in Alberta, British Columbia, and Ontario. Our objectives are to compare COVID-19 incidence rates from January 1 to July 6, 2020, and to estimate the time-varying reproduction number, *Rt*, of SARS-CoV-2 since the beginning of March 2020. ## Methods ### Scope and data sources We investigated the descriptive epidemiology and the *Rt* of SARS-CoV-2 in Alberta, British Columbia, and Ontario, and their public health subdivisions, using publicly available line lists of COVID-19 case data downloaded from provincial government websites. Alberta, British Columbia, and Ontario were chosen for this analysis given their line list data availability. These datasets contained no personal identifiers to protect patients’ privacy. Population data by age and sex for Alberta, British Columbia and Ontario was downloaded from Alberta government website.17 At the time of our study, Quebec has not yet made the line list of cases publicly available and was not included herein. As of early July, the cumulative number of cases in the other provinces were relatively small and therefore those provinces were not included. The Canadian government’s national line list dataset provided the accurate episode time in epidemiologic weeks instead of dates, and thus could not be analyzed using our method here. All the subprovincial incidence curves and *Rt* trajectories were presented in the Supplementary Materials. ### Timeframe The timeframe of this study started on January 1, 2020 (the accurate episode date of the earliest confirmed case in Ontario) and ended on July 5 (Alberta and Ontario) or 6 (British Columbia), the last episode date or reported date in each of our datasets. This timeframe is characterized by importation of cases from overseas followed by community transmission while social distancing measures were gradually introduced. March 17 was the day when the public health emergency was declared in Alberta, British Columbia and Ontario.4 By early July, some mandatory social distancing measures remain implemented therein. The timeframe for *Rt* estimation began with March 1. We did not estimate *Rt* in January and February because the small number of cases in Ontario and British Columbia during that time would lead to very uncertain *Rt* estimates. No cases were reported in Alberta before March. ### Ontario dataset The Ontario dataset included cases until the end of July 5, 2020 and was downloaded on July 6, 2020, at 5.34pm (EDT).18 In Ontario, the subprovincial unit is a Public Health Unit (PHU). The Ontario dataset contained the following variables: accurate episode date, case reported date, test reported date, specimen date, age group, gender, case acquisition information, clinical outcome, outbreak-related (yes or no), reporting PHU, and the address, city, postal code, website, latitude and longitude of the reporting PHU. In contrast to the Alberta and British Columbia datasets, which did not contain information on how the case-patients acquired the virus, the Ontario dataset provided information on how a case-patient acquired the virus (“Case Acquisition Info”): 1761 travel-related cases, 10839 “close contact” cases, 6687 cases with “no epi-link”, 14259 outbreak-related cases, 1120 cases with “no info-missing” and 1280 cases with “no info-unknown”. In our *Rt* estimation, we categorized all travel-related cases as “imported” cases and all others as “local” as per user instructions for the EpiEstim package. Please note that to run the EpiEstim package, the first case(s) on the first day of a time series must be rendered as “imported” cases. Thus, if the first case(s) of a time series was not travel-related, we manually denoted them as “imported” cases so that *Rt* can be estimated using the EpiEstim package. ### Alberta dataset The Alberta dataset included cases until the end of July 5, 2020 and was downloaded on July 7, 2020 at 7.59am (Eastern Daylight Time, EDT).19 In Alberta, each subprovincial unit is a Health Services Zone. The Alberta dataset contained the following variables: date of report, health services zone, gender, age group, case status (clinical outcome) and case type (confirmed or probable). The Alberta dataset contained 7678 confirmed cases and 711 probable cases; both types are included in the main analysis. In the sensitivity analysis, only confirmed cases were analyzed. ### British Columbia dataset The British Columbia dataset included cases until July 6, 2020 and was downloaded on July 7 at 8.01am (EDT).20 In British Columbia, each subprovincial unit is a Health Authority. The British Columbia dataset contained the following variables: date of report, health authorities, sex, age group and classification reported (“Epi-linked” versus laboratory-diagnosed). The British Columbia dataset contained 2970 laboratory-diagnosed cases (aka confirmed cases) and 8 “epi-linked” cases (aka probable cases); both are included in the main analysis. In the sensitivity analysis, only laboratory-diagnosed cases were analyzed. ### Incidence rate ratio The incidence rate was calculated by dividing the cumulative number of cases in a province by its total population in 2019. Given that the timeframe of our datasets is essentially the same (with only a difference in 1 day), we did not use person-time as the denominator. We used British Columbia as the reference, and calculate the incidence rate ratio (IRR) of Alberta and that of Ontario. We did not calculate the 95% confidence interval for the estimates, because the estimates were true population estimates and not sample estimates, as we used data from province-wide surveillance systems that cover everyone in the province. ### Case fatality ratio We limited our analysis to individuals whose clinical outcomes are known (either deaths or recovered), and used the following simple estimator to obtain an approximation of the case fatality ratio: case fatality ratio = number of deaths / (number of deaths + number of recovered).21 ### Time-varying reproduction number In contrast with the basic reproduction number, *R*, which represents the average number of secondary cases generated by an infectious individual in a totally susceptible population in the absence of interventions or behavioral changes, *Rt* is a time-varying indicator which represents the average number of secondary cases per infectious individual in a population as the epidemic unfolds in the presence of interventions and behavioral changes. *Rt* >1 indicates sustained transmission and epidemic growth; *Rt* <1 indicates unsustainable transmission and epidemic decline. As the epidemic runs its course, it is possible to quantify *Rt* over time after accounting for reductions in susceptibility in the population as more people acquire immunity through natural infection, behavior changes in the population, and interventions that mitigate the transmission rate.22 ### Instantaneous reproduction number method implemented in R package EpiEstim There are multiple statistical methods available for *Rt* estimation,23,24 one of which is the instantaneous reproduction number method implemented in the R package EpiEstim.25,26 This method is considered to be well-suited to near-real-time estimation of *Rt*, and is sensitive to signals of changes in transmission potential given recent implementation and cessation of interventions or behavioral changes.23 This method has been applied in various studies to estimate the *Rt* of SARS-CoV-2 in different locations globally.27-34 A summary of the instantaneous reproduction number method is provided in the Online Supplementary Materials. ### Software version Data management and analysis (except for maps) was performed using Microsoft Excel and R version 3.6.2 (R Core Team, Vienna, Austria). Version 2.2-3 of the EpiEstim package was used. Maps were made in R version 3.5.1 (R Core Team, Vienna, Austria). ### Ethics The Georgia Southern University Institutional Review Board made a non-human subjects determination for this project (H20364), under the G8 exemption category. ## Results The dates of the first cases and the total number of cases in our datasets by each subprovincial public health unit in Alberta, British Columbia and Ontario are presented in Table S2. ### Incidence rate ratio and case-fatality ratio at the provincial level Table S3 presents the number of cases per 100,000 population, stratified by sex/gender and age groups (see Tables S4-S6 for details). Interestingly, the incidence rates in Alberta and Ontario were higher than those in British Columbia, which served as the reference group. Among females, the respective IRR in Alberta and Ontario were 3.1 (187.3/59.9) and 4.3 (259.9/59.9). Among males, the respective IRR in Alberta and Ontario were 3.4 (196.2/57.2) and 4.0 (196.2/57.2). This may reflect a more aggressive testing strategy in Alberta and Ontario. As of July 11, 2020, Alberta and Ontario had the highest rate of testing in Canada (102,299 and 110,612 tests per million population, respectively) whereas in British Columbia the rate was only 36,059 tests per million population.6 In both provinces, anyone who feels needing a test, including asymptomatic individuals, can get a test without referral. For example, Ontario waived the requirement to get a referral to qualify for testing on May 24, 2020, and increased testing capacity to over 20,000 a day.35 Likewise, Alberta started testing asymptomatic cases around mid-May. The alternative hypothesis was that this might reflect a genuine difference across the 3 provinces in their epidemiology of COVID-19. Tables S7 and S8 present the clinical outcomes of cases in Alberta and Ontario. As expected, the elderly has a high case-fatality ratio than other age groups, as described in a recent study by Bignami and Van Assche,10 and this is similar to observations in other parts of the world, such as England.36 Likewise, as in other countries, there had been outbreaks in long-term care facilities in Canada.11,37,38 ### Epidemic curve and Rt The accurate episode date of the first two cases in Ontario was January 1, 2020, one without epi-link and one being a close contact. The “no epi-link” case was a male in his 80s whose specimen date was April 23 and whose test result was reported on April 24 (case reported on the same date). The “close contact” case was a male in his 50s whose specimen date was May 15 and whose test result was reported on May 17 (case reported on the same date). A third case was a male in his 50s whose accurate episode date was January 10 and was a “close contact” case but whose specimen was taken on June 9 (test reported and case reported on the next day). This cluster of 3 cases had episode dates before the first two reported cases in Ontario (date of report: January 23, 2020) whose accurate episode dates were January 21 and 22 respectively.1,2 In Ontario (Figure 1), *Rt* fluctuated around 1 in March, rose to >1 but <1.25 in early and mid-April, then dropped <1 in late April and early May. *Rt* rose to ~1 in mid-May and then remained <1 from late May through early July. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F1) Figure 1. The epidemic trajectory of coronavirus disease 2019 in Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). A map showing the cumulative number of cases over time by PHU in Ontario is presented in Figure 2, and a similar map showing the cumulative number of deaths is in Figure S1. Epidemic trajectories in selected PHUs are as follows. In City of Toronto PHU (Figure S2) and its neighbors, namely, Durham Regional PHU (Figure S3), Peel Regional PHU (Figure S4), and York Regional PHU (Figure S5), we observed similar tragectories: *Rt* ≥ 1 in sustained period of time from late March to mid-April. *Rt* subsided to <1 in late April and early May but rose to ~1 in mid-May. *Rt* decreased to <1 consistently in early June but it was on an increasing trend towards early July. Slightly further away from Toronto, in Halton Regional PHU (Figure S6), *Rt* fluctuates around 1 through the study period. It was <1 in much of June. In City of Hamilton PHU (Figure S7), *Rt* fluctuates around 1 through the timeframe of this study, except for mid-May when *Rt* reached a peak of 2 (1-window estimate). In Region of Waterloo PHU (Figure S8), *Rt* fluctuates around 1 in March but increases steadily to >1 in mid-April, followed by a continuous drop to <1, until late May when *Rt* rebounded to around 1 in June. In Niagara Region PHU (Figure S9), *Rt* fluctuates around 1 across the timeframe of this study, except for late May when there was a surge in cases leading to *Rt* >1 in late May followed by a trough <1 in June. In Ottawa PHU (Figure S10), *Rt* was <1 in late March but it increased to slightly >1 in April. It maintained <1 in much of May and fluctuates around 1 in June. For the other PHUs from Algoma PHU to Windsor-Essex County Health PHU, see Figures S11 to S35. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F2) Figure 2. Map of cumulative number of cases in Ontario by public health unit by the end of month and on July 5, 2020. The first case in British Columbia was reported by Vancouver Coastal Health Authority on January 26, 2020. In British Columbia, *Rt* dropped below 1 the first time in early April, but it increased towards the end of April. For most of May, *Rt* <1 while it fluctuated around 1 in June and early July (Figure 3). Figure 4 is a map showing the cumulative number of cases in British Columbia by Health Authority over time. Below are epidemic trajectories in selected Health Authorities. In Vancouver Coastal Health Authority (Figure S36), *Rt* was ~2 in mid-March and dropped to 1 by April and below <1 in mid-April but increased in late April to ~1. It fluctuates <1 for much for May but rose to ~2 in mid-June, followed by fluctuation around 1 in late June and early July. In Fraser Health Authority (Figure S37), *Rt* fluctuated between 1.5 and 2 in mid-March before dropping to ~1 through April. In May it fluctuated <1 but increased to 1 and fluctuated around 1 in June and early July. For Interior, Northern and Vancouver Island Health Authorities, see Figure S38 to S40. As a sensitivity analysis we analyzed the data excluding 8 probable cases (“epi-linked cases”) (Table S9); focusing at the 2970 confirmed cases, the results in *Rt* are similar (Figure S41). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F3) Figure 3. The epidemic trajectory of coronavirus disease 2019 in British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F4) Figure 4. Map of cumulative number of cases in British Columbia by health authority by the end of month and on July 6, 2020. The first case in Alberta was reported in Calgary Zone on March 6, 2020. In Alberta, *Rt* > 1 in March, corresponding to reaching a peak in incidence in late March and early April. A small dip in *Rt* happened in early April and *Rt* rose again in late April, as the incidence increased and reached a peak in late April. In much of May, *Rt* <1, but *Rt* increases in early June and fluctuates around 1 (slightly above 1) since mid-June (Figure 5). Figure 6 is a map showing the cumulative number of cases in Alberta by Health Services Zone over time, while Figure S42 is a similar map showing the cumulative number of deaths. Below are epidemic trajectories in selected Health Services Zones. In Calgary Zone (Figure S43), *Rt* began at ~2 in mid-March and dropped to ~1 in early April and rose again to ~1.5 and maintained >1 in the second half of April. In May and early July, *Rt* was <1, but increased to ~1 from mid-June through early July. In Edmonton Zone (Figure S44), it fluctuates around 1.5 in late March, dropped to 1 in early April and <1 in mid-April. *Rt* rose to ~1 in late April and then maintained at a <1 level for much of May. *Rt* increased to ~2 in late May and early June, and then gradually dropped to ~1 in early July. In South Zone (Figure S45), Rt fluctuated around 1 in late March and early April and rose to ~3 in the second half of April as there was an outbreak from mid-April to early May. Rt below <1 in much of May and fluctuated around 1 in much of June and made a sharp increase to ~2 in early July. For Central and North Zones, see Figures S46 and S47. As a sensitivity analysis we analyzed the data excluding 711 probable cases; focusing at 7678 confirmed cases only, the results in *Rt* are similar (Figure S48). ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F5) Figure 5. The epidemic trajectory of coronavirus disease 2019 in Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F6) Figure 6. Map of cumulative number of cases in Alberta by health service zone by the end of month and on July 5, 2020. ### Specimen collection date in Ontario dataset Among the 35948 case-patients in Ontario, the median delay in specimen collection for testing is 1 day after symptom onset (accurate episode date), with an interquartile range from 0 to 5 days. A decrease over time in the average time difference between date of symptom onset and the date of specimen taken, and the variance of that time difference (Figure 7 Panel a) represents the improvement in the testing capacity allowing individuals to get tested promptly after symptom onset. As seen in Figure 4a, the delay on March 1 would be 7.4 days and the delay was close to 0 around June 21. See Appendix for details of the regression model. The longest time difference (an outlier) was a delay of 151 days (that might suggest that these were serological tests and not polymerase chain reaction tests,39 but the dataset do not provide further details). A total of 305 (0.8%) individuals do not have specimen dates. There were 2391 (6.7%) whose specimens were collected prior to date of symptom onset, i.e. (accurate episode date - specimen date) < 0. Of these 2391 individuals, 926, 766, 397 and 302 individuals developed symptom 1, 2, 3 or ≥4 days, respectively, *after* having specimen taken for testing (Table S10). Most of these 2391 case-patients were outbreak-related (n=1473) or close contacts of cases (n=619), while 28 were travel-related. However, 211 of them did not have any epi-link, and 24 and 36 were without case acquisition info (“missing” or “unknown” respectively) (Table S11). This is not surprising, since Ontario enhanced testing capacity for testing asymptomatic people at risk of exposure. Of these 2391 individuals whose specimens taken prior to symptom onset, except for 8 who did not have a test report date, 1463 (61%) had their tests reported the same day as their symptom onset; 307 (13%), 1 day after symptom onset, 385 (16%), ≥2 days later, and 228 (9.5%) had their test reported *before* symptom onset (Table S10). Of these 2391 individuals, 1641 (69%) had their case reported the day they developed symptoms, 337 (14%) on the next day of symptom onset, 413 (17%), ≥2 days after symptom onset (Table S10). ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F7/graphic-7.medium.gif) [](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F7/graphic-7) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F7/graphic-8.medium.gif) [](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F7/graphic-8) Figure 7. **(a)** Scatter plot of the time difference between symptom onset (accurate episode date) and testing (specimen date) of each case-patient in Ontario over time (accurate episode date), Jan 1-July 6, 2020. **(b)** Scatter plot of the time difference between symptom onset (accurate episode date) and case report date of each case-patient in Ontario over time (accurate episode date), Jan 1-July 6, 2020. ### Delay in test report and case report in Ontario dataset The delay in test report decreased over time. On March 1, the delay was approximately 11 days, while on July 5, it was about 1 day (Figure S49). The delay in case report has also decreased over time. On March 1, the delay was approximately 10 days, while on July 5, it was about 1 day (Figure 7 panel b). See Appendix for details of the regression models. ## Discussion Alberta, British Columbia and Ontario, Canada, appeared to have managed to achieve limited transmission of SARS-CoV-2 with Rt~1 at the time of writing in early July 2020 (Figures 1, 3, 5), through a combination of multiple interventions, including community social distancing, increased testing of individuals at higher risk or suspected of SARS-CoV-2 infection, and isolation of COVID-19 cases. Canada’s experience with SARS epidemic in 2003 may have contributed to the coordinated government’s response in COVID-19. The 2003 SARS epidemic in Canada resulted in a set of recommendations to federal, provincial, and territorial leader, and in particular, led to the creation of the Public Health Agency of Canada which has been leading the response to COVID-19.40 Our results show that since mid-March in the 3 Canadian provinces under study, SARS-CoV-2 has been spreading with an *Rt* <2 that is smaller than the estimated R of >241,42 For the majority of the time period under study, *Rt* fluctuated around 1. Our analysis of the COVID-19 epidemic by Canadian subprovincial units present a similar picture as our province-level analysis with few exceptions. Nowhere in Canada experienced the trauma that residents in Wuhan, Milan and New York City had experienced. While Canadian provinces have substantial administrative authority over COVID-19 response, the federal and provincial government have collaborated to produce a coordinated response to the pandemic, as evidenced by the similar COVID-19 response measures in the three provinces, as well as the united politicians’ communications on the nature and severity of the COVID-19 pandemic, and the necessity of social distancing across the different Canadian political parties.43 As the three provinces began reopening businesses, services, and public spaces, all three provinces have been gradually relaunching by stages which depended on local infection rate and the capacity of the healthcare system. 44-46 Our study has several limitations. First, our *Rt* estimation for Alberta and British Columbia relied on data by reporting date and not date of onset, as the latter was unavailable to us. Only Ontario provided the accurate episode date. Future research can study the epidemic data by date of onset when it is available to researchers to better estimate *Rt*. Second, the number of reported cases per day in the Alberta and British Columbia datasets would be influenced by the testing rate, since limited testing capacity could lead to underdiagnosis. We did not attempt to do a nationwide analysis for Canada because the testing criteria has not been consistent across Canadian provinces. Third, our estimates are contingent that the testing ratio does not change significantly over the study period. However, provincial governments might have been ramping out testing capacity during the early epidemic phase in their jurisdiction. Fourth, reporting delays are known to vary over time and could influence the observed trajectory of the epidemics by date of report in Alberta and British Columbia. For Ontario, our results showed that the delay to testing has been reduced over time and Ontario achieved testing on the same or next date since symptom onset in June. It is possible that similar improvement in reduction in testing delay may have happened in Alberta and British Columbia, but we do not have such data. Fifth, improvement in testing over time could potentially reduce generation time over time. As we do not have data on infector-infectee pairs, we cannot estimate if serial interval has changed over time. In our *Rt* estimation, our serial interval distribution is parametrically determined *a priori* based on the literature.47 Sixth, the Ontario dataset is by date of symptom onset and therefore the most recent data points may be lower than what they should have been, because there could be cases who had not yet been reported to the provincial government. Thus, the low *Rt* estimate for Ontario in the 1st week of July may be an artifact due to right-censoring of the data. Seventh, we obtained an approximation of the case-fatality ratio using a simple estimator by restricting the analysis to cases with known outcomes. There are estimation methods that are more statistically advanced,21,48,49 but they are less accessible to practitioners of applied epidemiology and they are beyond the scope of this paper. ## Conclusions In conclusion, Alberta, British Columbia, and Ontario managed to keep the *Rt* of SARS-CoV-2 wavering around 1, with our recent estimates indicating sporadic episodes of sustained transmission in Alberta and British Columbia where *Rt* fluctuates around 1 or slightly >1 in the first week of July. In Ontario, *Rt* may have been suppressed <1 in early July. Sustained social distancing interventions and further limiting case importations could potentially suppress the epidemic by driving *Rt* <1 for a sustained period of time. ## Data Availability All the data used are publicly available data downloaded from websites of the provincial governments of Alberta, British Columbia and Ontario. The URL links to those websites are provided in the reference list of the manuscript. ## Disclaimers The opinions expressed in this article do not represent the official positions of the Centers for Disease Control and Prevention or the United States Government. ## Author Contributions Dr. Fung is the first and corresponding author. Dr. Hung, Ms. Ofori, Ms. Muniz-Rodriguez and Ms. Lai contributed equally as co-second authors. Dr. Chowell serves as the senior author. ### Instantaneous reproduction number method implemented in R package EpiEstim We summarized the instantaneous reproduction number method below, of which detailed description has been made available elsewhere.4,5 The instantaneous reproduction number is the ratio between the number of incident cases at time t, and the total infectiousness of all the infected individuals at time t: *Rt* = *It*/*Λt*, where ![Graphic][1]. The total infectiousness is the sum of infection incidence up to the previous time point (*t-1*) weighted by an infectivity function *ws*. The infectivity function is a probability distribution that describes the average infectiousness profile of an individual after being infected, and it is usually represented by the serial interval distribution as an approximation to the generation time distribution. The number of new cases at time *t* is therefore a Poisson distribution with a mean of ![Graphic][2]. Given the reproduction number *Rt*, and conditional on the time series of daily new cases, *I, …, It-1*, the likelihood of the incidence *It* is therefore, ![Graphic][3] where ![Graphic][4]. Therefore, ![Graphic][5]. We can interpret the instantaneous reproduction number as the average number of secondary cases that an infectious person will infect at time *t*, if conditions in previous time points remain the same at time *t*. Given this formulation of the instantaneous reproduction number, *Rt* is highly sensitive to fluctuation in the incidence data. To minimize *Rt* variation, the EpiEstim package allows user to specify a time window over which the *Rt* estimate is assumed to be constant. The window is of size *t* and ends at time t. If transmissibility is assumed constant over time from (*t*-*τ*+1) to *t*, and is denoted by a reproduction number, *Rt*,*τ* the likelihood of the incidence during the time period, from *I*(*t-τ+*1) to *It*, conditional on incidence prior to the time period, from *I*, to *I(t-τ)*, is found to be:5 ![Formula][6] The EpiEstim package utilized a Bayesian framework with a Gamma-distributed prior for *Rt,τ*, and the analytical expression of the posterior distribution of *Rt* and estimated its median, the variance and the 95% credible interval (CrI).5 Given the small number of cases in January and February and thus wide uncertainty bounds of the *Rt* estimates, we presented *Rt* estimates from March to early July. *Rt* estimates of two respective time windows of 7 days and 14 days, over which the *Rt* estimates is assumed to be constant, were presented with their 95% CrI, as in Muniz-Rodriguez et al.14 As a sensitivity analysis, the 14-day time window provides a more stable estimate as it would not be as responsive to sudden fluctuations in incident case count. We did not use a time window shorter than 7 days because such *Rt* estimate would be very sensitive to daily incidence fluctuation. Furthermore, for surveillance data by date of report (Alberta and British Columbia), weekend effect is common given the 5-day work week, i.e., daily number of cases reported over Saturday and Sunday would be lower than those from Monday to Friday. Therefore, it is reasonable to use time-windows of 7 or 14 days. Serial interval is assumed to follow a gamma distribution (mean: 4.60 days; standard deviation: 5.55).15 #### Map creation The health region boundaries shapefile was downloaded at Statistics Canada at [https://www150.statcan.gc.ca/n1/pub/82-402-x/2018001/hrbf-flrs-eng.htm](https://www150.statcan.gc.ca/n1/pub/82-402-x/2018001/hrbf-flrs-eng.htm). The health regions were based on 2016 Census geographic units. Maps were made in R version 3.5.1 (R Core Team, Vienna, Austria). The boundaries shapefiles are transformed by R packages rnaturalearth (version 0.1.0) and sf (version 0. 9.0), merged with the study data by health regions, and presented by R package ggplot2 (version 3.3.0). The maps show the cumulative number of cases and fatal cases on the last day of each month between January 2020 to June 2020, and the last date of the data (5 or 6 July 2020) with non-missing values. ### Time delay in specimen collection, test report and case report in Ontario, January 1-July 5, 2020 We ran linear regression models to obtain the trend line over time. 1. Time delay in specimen collection With (accurate episode date - specimen date) as the outcome variable (y) and accurate episode date as the predictor variable (x), the linear model is: y = −0.06627x + 1221.67026 Given that in R, date is counted from Jan 1, 1970 (when it is “day 0”), for January 1, 2020, x = 18262, y = 11.41 days; for March 1, 2020, x = 18322, y = 7.43 days; for June 21, 2020, x = 18434, y = 0.0062 day(s). 2. Time delay in test report With (accurate episode date - test report date) as y and accurate episode date as x, the linear model is: y = −0.08383x+1547.21994 For January 1, 2020, x = 18262, y = 16.31 days; for March 1, 2020, x = 18322, y = 11.29 days; for July 5, 2020, x = 18448, y = 0.72 day(s). 3. Time delay in case report With (accurate episode date - case report date) as y and accurate episode date as x, the linear model is: y = −0.07254+1339.42505 For January 1, 2020, x = 18262, y = 14.67 days; for March 1, 2020, x = 18322, y = 10.32 days; for July 5, 2020, x = 18448, y = 1.18 day(s). ## Supplementary Tables View this table: [Table S1.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T1) Table S1. Timeline of COVID-19-related events in Alberta, British Columbia, and Ontario in 2020. View this table: [Table S2.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T2) Table S2. Date of first cases and total number of cases as of July 5, 2020 (Alberta and Ontario) or July 6, 2020 (British Columbia), by provincial and subprovincial jurisdiction. View this table: [Table S3.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T3) Table S3. Cases per 100,000 and incidence rate ratio between Alberta, British Columbia and Ontario by sex (or gender) and age group. View this table: [Table S4.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T4) Table S4. Cumulative number of cases and cumulative incidence rate by age and gender in Alberta (N=8,389), as of July 5, 2020. View this table: [Table S5.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T5) Table S5. Cumulative number of cases and cumulative incidence rate by age and sex in British Columbia (N=2,978), as of July 6, 2020. View this table: [Table S6.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T6) Table S6. Cumulative number of cases and cumulative incidence rate by age and gender in Ontario, as of July 5, 2020. View this table: [Table S7.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T7) Table S7. Cases by age, type and clinical outcome in Alberta, as of July 5, 2020. View this table: [Table S8.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T8) Table S8. Cases by age and clinical outcome in Ontario, as of July 5, 2020. View this table: [Table S9.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T9) Table S9. Cases by age and classification reported in British Columbia, as of July 6, 2020. View this table: [Table S10.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T10) Table S10. Among 2391 Ontario cases whose specimen collection (specimen date) happened prior to symptom onset (accurate episode date), as of July 5, 2020. View this table: [Table S11.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/T11) Table S11. Among 2391 Ontario cases whose specimen collection (specimen date) happened prior to symptom onset (accurate episode date), as of July 5, 2020. ## Supplementary Figure legends ![Figure S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F8.medium.gif) [Figure S1.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F8) Figure S1. Map of cumulative number of deaths in Ontario by public health unit by the end of month and on July 6, 2020. ![Figure S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F9.medium.gif) [Figure S2.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F9) Figure S2. The epidemic trajectory of coronavirus disease 2019 in Toronto Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F10.medium.gif) [Figure S3.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F10) Figure S3. The epidemic trajectory of coronavirus disease 2019 in Durham Region Health Department, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F11.medium.gif) [Figure S4.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F11) Figure S4. The epidemic trajectory of coronavirus disease 2019 in Peel Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F12.medium.gif) [Figure S5.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F12) Figure S5. The epidemic trajectory of coronavirus disease 2019 in York Region Public Health Services, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F13.medium.gif) [Figure S6.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F13) Figure S6. The epidemic trajectory of coronavirus disease 2019 in Halton Public Health Services, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F14.medium.gif) [Figure S7.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F14) Figure S7. The epidemic trajectory of coronavirus disease 2019 in Hamilton Public Health Services, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F15.medium.gif) [Figure S8.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F15) Figure S8. The epidemic trajectory of coronavirus disease 2019 in Region of Waterloo, Public Health Services, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F16.medium.gif) [Figure S9.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F16) Figure S9. The epidemic trajectory of coronavirus disease 2019 in Niagara Region Public Health Department, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). North Bay Parry Sound District Health Unit ![Figure S10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F17.medium.gif) [Figure S10.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F17) Figure S10. The epidemic trajectory of coronavirus disease 2019 in Ottawa Public Health Services, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F18.medium.gif) [Figure S11.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F18) Figure S11. The epidemic trajectory of coronavirus disease 2019 in Algoma Public Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S12.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F19.medium.gif) [Figure S12.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F19) Figure S12. The epidemic trajectory of coronavirus disease 2019 in Brant County Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S13.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F20.medium.gif) [Figure S13.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F20) Figure S13. The epidemic trajectory of coronavirus disease 2019 in Chatham-Kent Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S14.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F21.medium.gif) [Figure S14.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F21) Figure S14. The epidemic trajectory of coronavirus disease 2019 in Eastern Ontario Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S15.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F22.medium.gif) [Figure S15.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F22) Figure S15. The epidemic trajectory of coronavirus disease 2019 in Grey Bruce Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S16.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F23.medium.gif) [Figure S16.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F23) Figure S16. The epidemic trajectory of coronavirus disease 2019 in Haldimand-Norfolk Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S17.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F24.medium.gif) [Figure S17.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F24) Figure S17. The epidemic trajectory of coronavirus disease 2019 in Haliburton, Kawartha, Pine Ridge District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S18.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F25.medium.gif) [Figure S18.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F25) Figure S18. The epidemic trajectory of coronavirus disease 2019 in Hastings and Prince Edward Counties Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S19.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F26.medium.gif) [Figure S19.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F26) Figure S19. The epidemic trajectory of coronavirus disease 2019 in Huron Perth District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S20.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F27.medium.gif) [Figure S20.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F27) Figure S20. The epidemic trajectory of coronavirus disease 2019 in Kingston, Frontenac and Lennox & Addington Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S21.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F28.medium.gif) [Figure S21.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F28) Figure S21. The epidemic trajectory of coronavirus disease 2019 in Lambton Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S22.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F29.medium.gif) [Figure S22.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F29) Figure S22. The epidemic trajectory of coronavirus disease 2019 in Leeds, Grenville and Lanark District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S23.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F30.medium.gif) [Figure S23.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F30) Figure S23. The epidemic trajectory of coronavirus disease 2019 in Middlesex-London Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S24.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F31.medium.gif) [Figure S24.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F31) Figure S24. The epidemic trajectory of coronavirus disease 2019 in North Bay Parry Sound District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S25.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F32.medium.gif) [Figure S25.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F32) Figure S25. The epidemic trajectory of coronavirus disease 2019 in Northwestern Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S26.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F33.medium.gif) [Figure S26.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F33) Figure S26. The epidemic trajectory of coronavirus disease 2019 in Peterborough Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S27.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F34.medium.gif) [Figure S27.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F34) Figure S27. The epidemic trajectory of coronavirus disease 2019 in Porcupine Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S28.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F35.medium.gif) [Figure S28.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F35) Figure S28. The epidemic trajectory of coronavirus disease 2019 in Renfrew County and District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S29.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F36.medium.gif) [Figure S29.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F36) Figure S29. The epidemic trajectory of coronavirus disease 2019 in Simcoe Muskoka District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S30.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F37.medium.gif) [Figure S30.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F37) Figure S30. The epidemic trajectory of coronavirus disease 2019 in Southwestern Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S31.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F38.medium.gif) [Figure S31.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F38) Figure S31. The epidemic trajectory of coronavirus disease 2019 in Sudbury & Districts Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S32.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F39.medium.gif) [Figure S32.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F39) Figure S32. The epidemic trajectory of coronavirus disease 2019 in Thunder Bay District Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S33.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F40.medium.gif) [Figure S33.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F40) Figure S33. The epidemic trajectory of coronavirus disease 2019 in Timiskaming Health Unit, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S34.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F41.medium.gif) [Figure S34.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F41) Figure S34. The epidemic trajectory of coronavirus disease 2019 in Windsor-Dufferin-Guelph Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S35.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F42.medium.gif) [Figure S35.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F42) Figure S35. The epidemic trajectory of coronavirus disease 2019 in Windsor-Dufferin-Guelph Public Health, Ontario, January 1-July 5, 2020: Upper panel: Daily number of new imported and local cases by date of symptom onset (aka accurate episode date) (upper panel); Rt with a 1-week window (middle panel); Rt with a 2-week window (lower panel). ![Figure S36.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F43.medium.gif) [Figure S36.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F43) Figure S36. The epidemic trajectory of coronavirus disease 2019 in Vancouver Coastal Health Authority, British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure S37.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F44.medium.gif) [Figure S37.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F44) Figure S37. The epidemic trajectory of coronavirus disease 2019 in Fraser Health Authority, British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure S38.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F45.medium.gif) [Figure S38.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F45) Figure S38. The epidemic trajectory of coronavirus disease 2019 in Interior Health Authority, British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure S39.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F46.medium.gif) [Figure S39.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F46) Figure S39. The epidemic trajectory of coronavirus disease 2019 in Northern Health Authority, British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure S40.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F47.medium.gif) [Figure S40.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F47) Figure S40. The epidemic trajectory of coronavirus disease 2019 in Vancouver Island Health Authority, British Columbia, January 1-July 6, 2020: Daily number of new cases by date of report (upper panel); *Rt* with a 1-week window (middle panel); *Rt* with a 2-week window (lower panel). ![Figure S41.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F48.medium.gif) [Figure S41.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F48) Figure S41. Sensitivity analysis (lab-diagnosed cases only) The epidemic trajectory of coronavirus disease 2019 in British Columbia, 2020: Daily number of new lab-diagnosed cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S42.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F49.medium.gif) [Figure S42.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F49) Figure S42. Map of cumulative number of deaths in Alberta by health service zone by the end of month and on July 6, 2020. ![Figure S43.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F50.medium.gif) [Figure S43.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F50) Figure S43. The epidemic trajectory of coronavirus disease 2019 in Calgary Zone, Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S44.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F51.medium.gif) [Figure S44.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F51) Figure S44. The epidemic trajectory of coronavirus disease 2019 in Edmonton Zone, Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S45.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F52.medium.gif) [Figure S45.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F52) Figure S45. The epidemic trajectory of coronavirus disease 2019 in South Zone, Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S46.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F53.medium.gif) [Figure S46.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F53) Figure S46. The epidemic trajectory of coronavirus disease 2019 in Central Zone, Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S47.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F54.medium.gif) [Figure S47.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F54) Figure S47. The epidemic trajectory of coronavirus disease 2019 in North Zone, Alberta, January 1-July 5, 2020: Daily number of new cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S48.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F55.medium.gif) [Figure S48.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F55) Figure S48. Sensitivity analysis (confirmed cases only) The epidemic trajectory of coronavirus disease 2019 in Alberta, 2020: Daily number of new confirmed cases by date of report (upper panel), *Rt* with a 1-week window (middle panel) and *Rt* with a 2-week window (lower panel). ![Figure S49.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/26/2020.07.18.20156992/F56.medium.gif) [Figure S49.](http://medrxiv.org/content/early/2020/08/26/2020.07.18.20156992/F56) Figure S49. Scatter plot of the time difference between symptom onset (accurate episode date) and test report date of each case-patient in Ontario over time (accurate episode date), Jan 1-July 6, 2020. ## Acknowledgements G.C. received support from NSF grant 1414374 as part of the joint NSF-NIH-USDA Ecology and Evolution of Infectious Diseases program. I.C.-H.F. received salary support from the Centers for Disease Control and Prevention (CDC) (19IPA1908208). This article is not part of I.C.-H.F’s CDC-sponsored projects. ## Footnotes * Isaac Chun-Hai Fung: [https://orcid.org/0000-0001-5496-2529](https://orcid.org/0000-0001-5496-2529) * Yuen Wai Hung: [https://orcid.org/0000-0002-1171-109X](https://orcid.org/0000-0002-1171-109X) * Sylvia Ofori: [https://orcid.org/0000-0002-9961-9975](https://orcid.org/0000-0002-9961-9975) * Kamalich Muniz-Rodriguez: [https://orcid.org/0000-0002-4786-9438](https://orcid.org/0000-0002-4786-9438) * Po Ying Lai: [https://orcid.org/0000-0002-6156-1878](https://orcid.org/0000-0002-6156-1878) * Gerardo Chowell: [https://orcid.org/0000-0003-2194-2251](https://orcid.org/0000-0003-2194-2251) * The labels of the y-axes in Figure 7A, Figure 7B and Figure S49, respectively, have been corrected. * Received July 18, 2020. * Revision received August 25, 2020. * Accepted August 26, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Silverstein WK, Stroud L, Cleghorn GE, Leis JA. First imported case of 2019 novel coronavirus in Canada, presenting as mild pneumonia. Lancet. 2020;395(10225):734. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30370-6&link_type=DOI) 2. 2.Marchand-Senecal X, Kozak R, Mubareka S, et al. Diagnosis and Management of First Case of COVID-19 in Canada: Lessons applied from SARS. Clin Infect Dis. 2020. 3. 3.National Post. As the COVID-19 pandemic hit, provinces declared states of emergency. Now many are up for renewal. April 15, 2020. 2020; [https://nationalpost.com/news/provincial-states-of-emergencies-were-issued-a-month-ago-most-are-coming-up-for-renewal](https://nationalpost.com/news/provincial-states-of-emergencies-were-issued-a-month-ago-most-are-coming-up-for-renewal). Accessed May 2, 2020. 4. 4.The Global and Mail. Coronavirus March 17: Ontario, Alberta, B.C. declare states of emergency; B.C. closes schools to combat virus spread. 2020;[https://www.theglobeandmail.com/canada/article-the-latest-on-the-coronavirus-cineplex-closes-its-doors-a-quieter-st/](https://www.theglobeandmail.com/canada/article-the-latest-on-the-coronavirus-cineplex-closes-its-doors-a-quieter-st/). Accessed April 7, 2020. 5. 5.Government of Canada. Coronavirus disease 2019 (COVID-19): Epidemiology update (as of July 10, 2020, 7pm EDT). 2020; [https://health-infobase.canada.ca/covid-19/epidemiological-summary-covid-19-cases.html](https://health-infobase.canada.ca/covid-19/epidemiological-summary-covid-19-cases.html). Accessed July 11, 2020. 6. 6.Government of Canada. Coronavirus disease (COVID-19): Outbreak update (as of July 12, 2020). 2020; [https://www.canada.ca/en/public-health/services/diseases/2019-novel-coronavirus-infection.html](https://www.canada.ca/en/public-health/services/diseases/2019-novel-coronavirus-infection.html). Accessed July 12, 2020. 7. 7.Government of Canada. Update on COVID-19 in Canada: Epidemiology and Modelling (July 8, 2020). 2020; [https://www.canada.ca/content/dam/phac-aspc/documents/services/diseases-maladies/coronavirus-disease-covid-19/epidemiological-economic-research-data/updatecovid\_19\_Canada\_epidemiology\_modelling\_20200708.pdf](https://www.canada.ca/content/dam/phac-aspc/documents/services/diseases-maladies/coronavirus-disease-covid-19/epidemiological-economic-research-data/updatecovid\_19\_Canada_epidemiology_modelling_20200708.pdf). Accessed July 13, 2020. 8. 8.Berry I, Soucy JR, Tuite A, Fisman D, Group C-CODW. Open access epidemiologic data and an interactive dashboard to monitor the COVID-19 outbreak in Canada. CMAJ. 2020;192(15):E420. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMToiMTkyLzE1L0U0MjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNi8yMDIwLjA3LjE4LjIwMTU2OTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 9. 9.Lapointe-Shaw L, Rader B, Astley CM, et al. Syndromic Surveillance for COVID-19 in Canada (Version 1). *medRxiv*. 2020;2020.05.19.20107391. 10. 10.Bignami S, Van Assche A. Assessing the burden of COVID-19 in Canada (Version 1). *medRxiv* 2020; [https://www.medrxiv.org/content/medrxiv/early/2020/06/16/2020.06.14.20130815.full.pdf](https://www.medrxiv.org/content/medrxiv/early/2020/06/16/2020.06.14.20130815.full.pdf), 2020.06.14.20130815. 11. 11.Brown KA, Jones A, Daneman N, et al. Association Between Nursing Home Crowding and COVID-19 Infection and Mortality in Ontario, Canada (Version 1). *medRxiv* 2020; [https://www.medrxiv.org/content/medrxiv/early/2020/06/23/2020.06.23.20137729.full.pdf](https://www.medrxiv.org/content/medrxiv/early/2020/06/23/2020.06.23.20137729.full.pdf), 2020.06.23.20137729. 12. 12.Mitra AR, Fergusson NA, Lloyd-Smith E, et al. Baseline characteristics and outcomes of patients with COVID-19 admitted to intensive care units in Vancouver, Canada: a case series. CMAJ. 2020;192(26):E694–E701. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMToiMTkyLzI2L0U2OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNi8yMDIwLjA3LjE4LjIwMTU2OTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 13. 13.Tuite AR, Fisman DN, Greer AL. Mathematical modelling of COVID-19 transmission and mitigation strategies in the population of Ontario, Canada. CMAJ. 2020;192(19):E497–E505. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMToiMTkyLzE5L0U0OTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNi8yMDIwLjA3LjE4LjIwMTU2OTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 14. 14.Institut national de santé publique du Québec. Mesures de distanciation et de confinement au Québec: impact et projections. 2020; [https://www.inspq.qc.ca/covid-19/donnees/projections/distanciation](https://www.inspq.qc.ca/covid-19/donnees/projections/distanciation). Accessed July 13, 2020. 15. 15.British Columbia Centre for Disease Control. Modelling & Projections. 2020; [http://www.bccdc.ca/health-info/diseases-conditions/covid-19/modelling-projections](http://www.bccdc.ca/health-info/diseases-conditions/covid-19/modelling-projections). Accessed July 13, 2020. 16. 16.Government of Alberta. COVID-19 Modelling (Update, April 8, 2020). 2020; [https://www.alberta.ca/assets/documents/covid-19-case-modelling-projection.pdf](https://www.alberta.ca/assets/documents/covid-19-case-modelling-projection.pdf). Accessed July 13, 2020. 17. 17.Government of Alberta. Population by age and sex: 1971 - 2019. 2020; [https://open.alberta.ca/dataset/alberta-population-estimates-data-tables/resource/7e96a710-3725-45d1-ba58-36ea9cfe144c](https://open.alberta.ca/dataset/alberta-population-estimates-data-tables/resource/7e96a710-3725-45d1-ba58-36ea9cfe144c). Accessed July 8, 2020. 18. 18.Government of Ontario. Confirmed positive cases of COVID-19 in Ontario. 2020; [https://data.ontario.ca/dataset/confirmed-positive-cases-of-covid-19-in-ontario](https://data.ontario.ca/dataset/confirmed-positive-cases-of-covid-19-in-ontario). Accessed July 6, 2020. 19. 19.Government of Alberta. COVID-19 Alberta statistics. 2020; [https://www.alberta.ca/stats/covid-19-alberta-statistics.htm](https://www.alberta.ca/stats/covid-19-alberta-statistics.htm). Accessed July 7, 2020. 20. 20.British Columbia Centre for Disease Control. BC COVID-19 Data. 2020; [http://www.bccdc.ca/health-info/diseases-conditions/covid-19/data](http://www.bccdc.ca/health-info/diseases-conditions/covid-19/data). Accessed July 7, 2020. 21. 21.Ghani AC, Donnelly CA, Cox DR, et al. Methods for estimating the case fatality ratio for a novel, emerging infectious disease. Am J Epidemiol. 2005;162(5):479–486. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwi230&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16076827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231363800012&link_type=ISI) 22. 22.Vynnycky E, White RG. An Introduction to Infectious Disease Modelling. Oxford: Oxford University Press; 2010. 23. 23.Gostic KM, McGough L, Baskerville E, et al. Practical considerations for measuring the effective reproductive number, Rt. *medRxiv*. 2020;2020.06.18.20134858. 24. 24.O’Driscoll M, Harry C, Donnelly CA, Cori A, Dorigatti I. A comparative analysis of statistical methods to estimate the reproduction number in emerging epidemics with implications for the current COVID-19 pandemic. *medRxiv*. 2020;2020.05.13.20101121. 25. 25.Thompson RN, Stockwin JE, van Gaalen RD, et al. Improved inference of time-varying reproduction numbers during infectious disease outbreaks. Epidemics. 2019;29:100356. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100356&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 26. 26.Cori A, Ferguson NM, Fraser C, Cauchemez S. A new framework and software to estimate time-varying reproduction numbers during epidemics. Am J Epidemiol. 2013; 178(9): 1505–1512. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 27. 27.Abbott S, Hellewell J, Thompson R, et al. Estimating the time-varying reproduction number of SARS-CoV-2 using national and subnational case counts [version 1; peer review: awaiting peer review]. 2020;5(112). 28. 28.Leung K, Wu JT, Liu D, Leung GM. First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment. Lancet. 2020;395(10233): 1382–1393. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30746-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32277878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 29. 29.Cowling BJ, Ali ST, Ng TWY, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study. Lancet Public Health. 2020;5(5):e279–e288. 30. 30.Najafi F, Izadi N, Hashemi-Nazari SS, Khosravi-Shadmani F, Nikbakht R, Shakiba E. Serial interval and time-varying reproduction number estimation for COVID-19 in western Iran. New Microbes New Infect. 2020;36:100715. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 31. 31.Zhuang Z, Zhao S, Lin Q, et al. Preliminary estimates of the reproduction number of the coronavirus disease (COVID-19) outbreak in Republic of Korea and Italy by 5 March 2020. Int J Infect Dis. 2020;95:308–310. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.04.044&link_type=DOI) 32. 32.Moirano G, Schmid M, Barone-Adesi F. Short-Term Effects of Mitigation Measures for the Containment of the COVID-19 Outbreak: An Experience From Northern Italy. Disaster Med Public Health Prep. 2020:1–2. 33. 33.Adegboye OA, Adekunle AI, Gayawan E. Early Transmission Dynamics of Novel Coronavirus (COVID-19) in Nigeria. Int J Environ Res Public Health. 2020;17(9):3054. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph17093054&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32353991&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 34. 34.Scire J, Nadeau S, Vaughan T, et al. Reproductive number of the COVID-19 epidemic in Switzerland with a focus on the Cantons of Basel-Stadt and Basel-Landschaft. Swiss Med Wkly. 2020;150:w20271. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4414/smw.2020.20271&link_type=DOI) 35. 35.Government of Ontario. Protecting Ontarians Through Enhanced Testing. 2020; [https://files.ontario.ca/moh-covid-19-protecting-ontarians-through-enhanced-testing-en-2020-05-29.pdf](https://files.ontario.ca/moh-covid-19-protecting-ontarians-through-enhanced-testing-en-2020-05-29.pdf). Accessed July 12, 2020. 36. 36.Williamson EJ, Walker AJ, Bhaskaran K, et al. OpenSAFELY: factors associated with COVID-19 death in 17 million patients. Nature. 2020;DOI:[https://doi.org/10.1038/s41586-020-2521-4](https://doi.org/10.1038/s41586-020-2521-4). 37. 37.Bensadoun E. Nearly half of Canada’s COVID-19 deaths linked to long-term care facilities: Tam. Global News 2020; [https://globalnews.ca/news/6811726/coronavirus-long-term-care-deaths-canada/](https://globalnews.ca/news/6811726/coronavirus-long-term-care-deaths-canada/). Accessed July 13, 2020. 38. 38.Hsu AT, Lane N. Impact of COVID-19 on residents of Canada’s long-term care homes - ongoing challenges and policy response. Report in LTCcovid.org, International Long-Term Care Policy Network, CPEC-LSE, 23 April 2020. 2020; [https://ltccovid.org/wp-content/uploads/2020/04/LTC-COVID19-situation-in-Canada-22-April-2020-1.pdf](https://ltccovid.org/wp-content/uploads/2020/04/LTC-COVID19-situation-in-Canada-22-April-2020-1.pdf). Accessed July 13, 2020. 39. 39.Government of Canada. Testing devices for COVID-19: Serological testing devices. 2020; [https://www.canada.ca/en/health-canada/services/drugs-health-products/covid19-industry/medical-devices/testing/serological.html](https://www.canada.ca/en/health-canada/services/drugs-health-products/covid19-industry/medical-devices/testing/serological.html). Accessed July 13, 2020. 40. 40.Webster P. Canada and COVID-19: learning from SARS. Lancet. 2020;395(10228):936–937. 41. 41.Tsang TK, Wu P, Lin Y, Lau EHY, Leung KSM, Cowling BJ. Effect of changing case definitions for COVID-19 on the epidemic curve and transmission parameters in mainland China: a modelling study. Lancet Public Health. 2020:DOI:10.1016/S2468-2667(1020)30089-X. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-2667(1020)30089-X&link_type=DOI) 42. 42.Muniz-Rodriguez K, Fung IC-H, Ferdosi SR, et al. Severe Acute Respiratory Syndrome Coronavirus 2 Transmission Potential, Iran, 2020. Emerg Infect Dis. 2020;26(8):DOI:[https://doi.org/10.3201/eid2608.200536](https://doi.org/10.3201/eid2608.200536). 43. 43.Merkley E, Bridgman A, Loewen PJ, Owen T, Ruths D, Zhilin O. A Rare Moment of CrossPartisan Consensus: Elite and Public Response to the COVID-19 Pandemic in Canada. Canadian Journal of Political Science. 2020:1–8. 44. 44.Government of Alberta. Opening soon: Alberta’s relaunch strategy (Updated June 2020). 2020; [https://open.alberta.ca/dataset/61f54c09-d6d7-4a12-a5be-0bc663a02c31/resource/e158ff14-eab7-4f24-94f4-b67c3639d0d5/download/covid-19-alberta-relaunch-strategy-2020-06.pdf](https://open.alberta.ca/dataset/61f54c09-d6d7-4a12-a5be-0bc663a02c31/resource/e158ff14-eab7-4f24-94f4-b67c3639d0d5/download/covid-19-alberta-relaunch-strategy-2020-06.pdf) Accessed July 13, 2020. 45. 45.Government of Ontario. Reopening Ontario. 2020; [https://www.ontario.ca/page/reopening-ontario](https://www.ontario.ca/page/reopening-ontario). Accessed July 13, 2020. 46. 46.Government of British Columbia. BC’s Restart Plan. 2020; [https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/bc-restart-plan](https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/bc-restart-plan). Accessed July 13, 2020. 47. 47.You C, Deng Y, Hu W, et al. Estimation of the time-varying reproduction number of COVID-19 outbreak in China. Int J Hyg Environ Health. 2020;228:113555. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijheh.2020.113555&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 48. 48.Garske T, Legrand J, Donnelly CA, et al. Assessing the severity of the novel influenza A/H1N1 pandemic. BMJ. 2009;339:b2840. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzMzkvanVsMTRfMy9iMjg0MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzI2LzIwMjAuMDcuMTguMjAxNTY5OTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 49. 49.Mizumoto K, Chowell G. Estimating Risk for Death from Coronavirus Disease, China, January-February 2020. Emerg Infect Dis. 2020;26(6):1251–1256. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/EID2606.200233&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) ## References 1. 1.Vynnycky E, White RG. An Introduction to Infectious Disease Modelling. Oxford: Oxford University Press; 2010. 2. 2.Gostic KM, McGough L, Baskerville E, et al. Practical considerations for measuring the effective reproductive number, Rt. *medRxiv*. 2020;2020.06.18.20134858. 3. 3.O’Driscoll M, Harry C, Donnelly CA, Cori A, Dorigatti I. A comparative analysis of statistical methods to estimate the reproduction number in emerging epidemics with implications for the current COVID-19 pandemic. *medRxiv*. 2020;2020.05.13.20101121. 4. 4.Thompson RN, Stockwin JE, van Gaalen RD, et al. Improved inference of time-varying reproduction numbers during infectious disease outbreaks. Epidemics. 2019;29:100356. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100356&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 5. 5.Cori A, Ferguson NM, Fraser C, Cauchemez S. A new framework and software to estimate time-varying reproduction numbers during epidemics. Am J Epidemiol. 2013; 178(9): 1505–1512. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 6. 6.Abbott S, Hellewell J, Thompson R, et al. Estimating the time-varying reproduction number of SARS-CoV-2 using national and subnational case counts [version 1; peer review: awaiting peer review]. 2020;5(112). 7. 7.Leung K, Wu JT, Liu D, Leung GM. First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment. Lancet. 2020;395(10233): 1382–1393. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30746-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32277878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 8. 8.Cowling BJ, Ali ST, Ng TWY, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study. Lancet Public Health. 2020;5(5):e279–e288. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s2468-2667(20)30090-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 9. 9.Najafi F, Izadi N, Hashemi-Nazari SS, Khosravi-Shadmani F, Nikbakht R, Shakiba E. Serial interval and time-varying reproduction number estimation for COVID-19 in western Iran. New Microbes New Infect. 2020;36:100715. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 10. 10.Zhuang Z, Zhao S, Lin Q, et al. Preliminary estimates of the reproduction number of the coronavirus disease (COVID-19) outbreak in Republic of Korea and Italy by 5 March 2020. Int J Infect Dis. 2020;95:308–310. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.04.044&link_type=DOI) 11. 11.Moirano G, Schmid M, Barone-Adesi F. Short-Term Effects of Mitigation Measures for the Containment of the COVID-19 Outbreak: An Experience From Northern Italy. Disaster Med Public Health Prep. 2020:1–2. 12. 12.Adegboye OA, Adekunle AI, Gayawan E. Early Transmission Dynamics of Novel Coronavirus (COVID-19) in Nigeria. Int J Environ Res Public Health. 2020;17(9):3054. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph17093054&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32353991&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 13. 13.Scire J, Nadeau S, Vaughan T, et al. Reproductive number of the COVID-19 epidemic in Switzerland with a focus on the Cantons of Basel-Stadt and Basel-Landschaft. Swiss Med Wkly. 2020;150:w20271. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4414/smw.2020.20271&link_type=DOI) 14. 14.Muniz-Rodriguez K, Chowell G, Schwind JS, et al. Time-varying reproduction numbers of COVID-19 in Georgia, USA, March 2-June 14, 2020: comparison of two estimation methods. 2020; [https://www.medrxiv.org/content/10.1101/2020.05.19.20107219v2](https://www.medrxiv.org/content/10.1101/2020.05.19.20107219v2), 2020.05.19.20107219. 15. 15.You C, Deng Y, Hu W, et al. Estimation of the time-varying reproduction number of COVID-19 outbreak in China. Int J Hyg Environ Health. 2020;228:113555. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijheh.2020.113555&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F26%2F2020.07.18.20156992.atom) 16. 16.Government of Canada. Coronavirus disease (COVID-19): Outbreak update (as of July 12, 2020). 2020; [https://www.canada.ca/en/public-health/services/diseases/2019-novel-coronavirus-infection.html](https://www.canada.ca/en/public-health/services/diseases/2019-novel-coronavirus-infection.html). Accessed July 12, 2020. 17. 17.Government of Ontario. Protecting Ontarians Through Enhanced Testing. 2020; [https://files.ontario.ca/moh-covid-19-protecting-ontarians-through-enhanced-testing-en-2020-05-29.pdf](https://files.ontario.ca/moh-covid-19-protecting-ontarians-through-enhanced-testing-en-2020-05-29.pdf). Accessed July 12, 2020. 18. 18.Bignami S, Van Assche A. Assessing the burden of COVID-19 in Canada (Version 1). *medRxiv* 2020; [https://www.medrxiv.org/content/medrxiv/early/2020/06/16/2020.06.14.20130815.full.pdf](https://www.medrxiv.org/content/medrxiv/early/2020/06/16/2020.06.14.20130815.full.pdf), 2020.06.14.20130815. 19. 19.Williamson EJ, Walker AJ, Bhaskaran K, et al. OpenSAFELY: factors associated with COVID-19 death in 17 million patients. Nature. 2020;DOI: [https://doi.org/10.1038/s41586-020-2521-4](https://doi.org/10.1038/s41586-020-2521-4). 20. 20.Bensadoun E. Nearly half of Canada’s COVID-19 deaths linked to long-term care facilities: Tam. Global News 2020; [https://globalnews.ca/news/6811726/coronavirus-long-term-care-deaths-canada/](https://globalnews.ca/news/6811726/coronavirus-long-term-care-deaths-canada/).Accessed July 13, 2020. 21. 21.Hsu AT, Lane N. Impact of COVID-19 on residents of Canada’s long-term care homes - ongoing challenges and policy response. Report in LTCcovid.org, International Long-Term Care Policy Network, CPEC-LSE, 23 April 2020. 2020; [https://ltccovid.org/wp-content/uploads/2020/04/LTC-CQVID19-situation-in-Canada-22-April-2020-1.pdf](https://ltccovid.org/wp-content/uploads/2020/04/LTC-CQVID19-situation-in-Canada-22-April-2020-1.pdf). Accessed July 13, 2020. 22. 22.Brown KA, Jones A, Daneman N, et al. Association Between Nursing Home Crowding and COVID-19 Infection and Mortality in Ontario, Canada (Version 1). *medRxiv* 2020; [https://www.medrxiv.org/content/medrxiv/early/2020/06/23/2020.06.23.20137729.full.pdf](https://www.medrxiv.org/content/medrxiv/early/2020/06/23/2020.06.23.20137729.full.pdf), 2020.06.23.20137729. 23. 23.Silverstein WK, Stroud L, Cleghorn GE, Leis JA. First imported case of 2019 novel coronavirus in Canada, presenting as mild pneumonia. Lancet. 2020;395(10225):734. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30370-6&link_type=DOI) 24. 24.Marchand-Senecal X, Kozak R, Mubareka S, et al. Diagnosis and Management of First Case of COVID-19 in Canada: Lessons applied from SARS. Clin Infect Dis. 2020. 25. 25.Government of Canada. Testing devices for COVID-19: Serological testing devices. 2020; [https://www.canada.ca/en/health-canada/services/drugs-health-products/covid19-industry/medical-devices/testing/serological.html](https://www.canada.ca/en/health-canada/services/drugs-health-products/covid19-industry/medical-devices/testing/serological.html). Accessed July 13, 2020. 26. 26.Webster P. Canada and COVID-19: learning from SARS. Lancet. 2020;395(10228):936–937. 27. 27.Tsang TK, Wu P, Lin Y, Lau EHY, Leung KSM, Cowling BJ. Effect of changing case definitions for COVID-19 on the epidemic curve and transmission parameters in mainland China: a modelling study. Lancet Public Health. 2020:DOI:10.1016/S2468-2667(1020)30089-X. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-2667(1020)30089-X&link_type=DOI) 28. 28.Muniz-Rodriguez K, Fung IC-H, Ferdosi SR, et al. Severe Acute Respiratory Syndrome Coronavirus 2 Transmission Potential, Iran, 2020. Emerg Infect Dis. 2020;26(8):DOI:[https://doi.org/10.3201/eid2608.200536](https://doi.org/10.3201/eid2608.200536). 29. 29.Merkley E, Bridgman A, Loewen PJ, Owen T, Ruths D, Zhilin O. A Rare Moment of CrossPartisan Consensus: Elite and Public Response to the COVID-19 Pandemic in Canada. Canadian Journal of Political Science. 2020:1–8. 30. 30.Government of Alberta. Opening soon: Alberta’s relaunch strategy (Updated June 2020). 2020; [https://open.alberta.ca/dataset/61f54c09-d6d7-4a12-a5be-0bc663a02c31/resource/e158ff14-eab7-4f24-94f4-b67c3639d0d5/download/covid-19-alberta-relaunch-strategy-2020-06.pdf](https://open.alberta.ca/dataset/61f54c09-d6d7-4a12-a5be-0bc663a02c31/resource/e158ff14-eab7-4f24-94f4-b67c3639d0d5/download/covid-19-alberta-relaunch-strategy-2020-06.pdf). Accessed July 13, 2020. 31. 31.Government of Ontario. Reopening Ontario. 2020; [https://www.ontario.ca/page/reopening-ontario](https://www.ontario.ca/page/reopening-ontario). Accessed July 13, 2020. 32. 32.Government of British Columbia. BC’s Restart Plan. 2020; [https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/bc-restart-plan](https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/bc-restart-plan). Accessed July 13, 2020. 33. 33.CTV News. Tracking every case of COVID-19 in Canada. 2020; [https://www.ctvnews.ca/health/coronavirus/tracking-every-case-of-covid-19-in-canada-1.4852102](https://www.ctvnews.ca/health/coronavirus/tracking-every-case-of-covid-19-in-canada-1.4852102). Accessed May 2, 2020. 34. 34.Global News. Timeline: How Canada has changed since coronavirus was declared a pandemic. April 11 2020.. 2020; [https://globalnews.ca/news/6800118/pandemic-one-month-timeline/](https://globalnews.ca/news/6800118/pandemic-one-month-timeline/). Accessed May 2, 2020. 35. 35.Public Health Ontario. Coronavirus Disease 2019 (COVID-19) Testing. 2020; [https://www.publichealthontario.ca/en/laboratory-services/test-information-index/wuhan-novel-coronavirus](https://www.publichealthontario.ca/en/laboratory-services/test-information-index/wuhan-novel-coronavirus). Accessed May 2, 2020. 36. 36.The Global and Mail. Coronavirus March 17: Ontario, Alberta, B.C. declare states of emergency; B.C. closes schools to combat virus spread. 2020;[https://www.theglobeandmail.com/canada/article-the-latest-on-the-coronavirus-cineplex-closes-its-doors-a-quieter-st/](https://www.theglobeandmail.com/canada/article-the-latest-on-the-coronavirus-cineplex-closes-its-doors-a-quieter-st/). Accessed April 7, 2020. 37. 37.Global News. Alberta closes some non-essential businesses, prevents evictions as 542 COVID-19 cases confirmed. March 27, 2020. 2020; [https://globalnews.ca/news/6742251/alberta-health-coronavirus-covid-19-march-27/](https://globalnews.ca/news/6742251/alberta-health-coronavirus-covid-19-march-27/).Accessed May 2, 2020. 38. 38.Vogel L. COVID-19: A timeline of Canada’s first-wave response. CMAJ News 2020;[https://cmajnews.com/2020/06/12/coronavirus-1095847/](https://cmajnews.com/2020/06/12/coronavirus-1095847/). Accessed July 13, 2020. 39. 39.Government of Ontario. Reopening Ontario in stages: What can open in Stage 1. 2020; [https://www.ontario.ca/page/reopening-ontario-stagestfsection-5](https://www.ontario.ca/page/reopening-ontario-stagestfsection-5). Accessed July 13, 2020. 40. 40.Government of Ontario. Reopening Ontario: what’s in each stage: What can open in Stage 2. 2020; [https://www.ontario.ca/page/reopening-ontario-whats-each-stagetfsection-2](https://www.ontario.ca/page/reopening-ontario-whats-each-stagetfsection-2). Accessed July 13, 2020. 41. 41.Government of British Columbia. Phase 3 - BC’s Restart Plan. 2020; [https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/phase-3](https://www2.gov.bc.ca/gov/content/safety/emergency-preparedness-response-recovery/covid-19-provincial-support/phase-3). Accessed July 13, 2020. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: /embed/graphic-9.gif