Real-time monitoring the transmission potential of COVID-19 in Singapore, February 2020 ======================================================================================= * Amna Tariq * Yiseul Lee * Kimberlyn Roosa * Seth Blumberg * Ping Yan * Stefan Ma * Gerardo Chowell ## Abstract **Background** The ongoing COVID-19 epidemic that has spread widely in China since December 2019 is now generating local transmission in several countries including Singapore of February 19, 2020. This highlights the need to monitor transmission potential of unfolding SARS-CoV-2 outbreaks in real time. In particular, five major COVID-19 clusters have emerged in Singapore. **Methods** Here we estimate the effective reproduction number, Rt, for Singapore from the daily case series of imported and autochthonous cases by date of symptoms onset, after adjusting the local cases for reporting delays, using a generalized growth model and employing a renewal equation. We also derive the reproduction number from the distribution of cluster sizes using a branching process analysis. **Results** As of February 19th, 2020, our results indicate that the mean effective reproduction number Rt is at 1.1 (95% CI: 1.1,1.3) while an estimate of Rt based on cluster size distribution is at 0.63 (95% CI: 0.4, 1.05). **Conclusion** These estimates underscore positive impact of containment efforts in Singapore while at the same time suggest the need to sustain social distancing and active case finding efforts to stomp out all active chains of transmission. Keywords * SARS-CoV-2 * COVID-19 * Singapore * Transmission heterogeneity * Reproduction number * Cluster distribution * Reporting delay ## Background The ongoing COVID-19 epidemic started with a cluster of pneumonia cases of unknown etiology in Wuhan, China back in December 2019 (1, 2). The initial cases have been linked to a wet market in the city of Wuhan, pointing to an animal source of the epidemic (3). Subsequently, rapid human-to-human transmission of the disease was confirmed in January 2020, and the etiological agent was identified as severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) due to its genetic similarity to the SARS-CoV discovered in 2003 (1, 4). The total case tally has reached 75,700 infections including 2128 deaths as of February 19, 2020 (1, 2). As the virus continues to spread in the human population, obtaining an accurate “real-time” picture of the epidemic’s trajectory is complicated by several factors including reporting delays and changes in the case definition (5, 6). Although the epidemic appears to be on a downward trend in China, cases of the novel coronavirus have been reported in 26 countries around the globe, and local transmission is ongoing in Singapore, Japan, Korea and Thailand (7). Detailed epidemiological data from these countries is needed to monitor transmission potential of SARS-CoV-2 in near real-time. Outside of China, Singapore is one of the hardest hit countries with the COVID-19 epidemic (8), reporting its first symptomatic imported case (a 66 years old Chinese male) on January 23, 2020 (9). Up until February 19, 2020, Singapore has reported 84 laboratory confirmed cases, including 22 reported case importations from Wuhan China (8). Imported cases include six evacuees who departed China between January 30 and February 9, 2020 (10, 11). On February 4, 2020, the Ministry of Health of Singapore reported its first local cluster of COVID-19, which was linked to the Yong Thai Hang shop (12). A total of 9 clusters with 1 or more COVID-19 local cases have been reported thus far. Table 1 summarizes the characteristics of the clusters. View this table: [Table 1:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/T1) Table 1: Characteristics of the largest COVID-19 outbreak in Singapore as of February 19, 2020. Although large-scale community transmission has not been reported in Singapore, the novel coronavirus can rapidly spread in confined and crowded places, as illustrated by a large cluster of COVID-19 cases that has been linked to the Grace Assembly of God Church, which has an average attendance of 4,800 (13). In China, substantial hospital-based transmission of SARS-CoV-2 has been reported, with a total of 1700 cases reported among healthcare workers (14, 15). This pattern aligns well with past outbreaks of SARS and MERS (16), including substantial nosocomial transmission during the 2003 SARS outbreak in Singapore (17). Thus, the distributions of cluster sizes for SARS and MERS outbreaks are characterized by high heterogeneity (16). To minimize the risk of hospital-based transmission of SARS-CoV-2, the Ministry of Health of Singapore has restricted the movement of patients and staff across hospitals (18). Also, because multiple unlinked community cases of COVID-19 have been reported (19), various social distancing measures have been put in place including advising the public against large social gatherings in order to mitigate the risk of community transmission (20, 21). These social distancing measures reduce the risk of onward transmission not only within Singapore, but also beyond the borders of this highly connected nation (22). The reproduction number is a key threshold quantity to assess the transmission potential of an emerging disease such as COVID-19 (23, 24). It quantifies the average number of secondary cases generated per case. If the reproduction number is below 1.0, infections occur in isolated clusters as self-limited chains of transmission, and persistence of the disease would require continued undetected importations. On the other hand, reproduction numbers above 1.0 indicate sustained community transmission (16, 24). Using epidemiological data and mathematical modeling tools, we are monitoring the transmission potential of SARS-CoV-2 in Singapore in real-time, and here we report results as of February 19, 2020. Specifically, we assess the effective reproduction number from the daily case series of imported and autochthonous cases by date of symptoms onset after adjusting for reporting delays, and we also derive an estimate of the reproduction number based on the characteristics of the clusters of COVID-19 in Singapore. ## Methods ### Data We obtained the daily series of confirmed cases of COVID-19 in Singapore as of February 19, 2020. Individual-level case details including the dates of symptoms onset, the date of reporting, and whether the case is autochthonous (local transmission) or imported are publicly available from the Ministry of Health, Singapore (25). Clusters consisting of one or more cases according to the infection source were also assembled from case descriptions obtained from field investigations conducted by the Ministry of Health, Singapore (25). Single imported cases or unlinked cases are analyzed as clusters of size 1. ### Adjusting for reporting delays As an outbreak progresses in real time, epidemiological curves can be distorted by reporting delays arising from several factors that include: (i) delays in case detection during field investigations, (ii) delays in symptom onset after infection, (iii) delays in seeking medical care, (iv) delays in diagnostics and (v) delays in processing data in surveillance systems (26). However, it is possible to adjust epidemic curves by reporting delays using statistical methods (27). The reporting delay for a case is defined as the time lag in days between the date of onset and date of reporting. Here we adjusted the COVID-19 epidemic curve of local cases by reporting delays using a non-parametric method that employs survival analysis techniques (28-30). The 95% prediction limits are derived according to Lawless et al. (31). For this analysis, we exclude 6 imported cases and 3 local cases without dates of symptoms onset. ### Effective reproduction number from case incidence We assess the effective reproduction number over time, Rt, which quantifies the temporal variation in the average number of secondary cases generated per case during the course of an outbreak after considering multiple factors including behavior changes, cultural factors, and the implementation of public health measures (16, 24, 32). Estimates of Rt>1 indicate sustained transmission; whereas, Rt <1 implies that the outbreak is slowing down and the incidence trend is declining. Hence, maintaining Rt <1 is required to bring an outbreak under control. Using the reporting delay adjusted incidence curve, we estimate the most recent estimate of Rt for COVID-19 in Singapore. We first characterize daily incidence of local cases using the generalized growth model (GGM) after adjusting for imported cases. This model characterizes the growth profile via two parameters: the growth rate (r) and the scaling of the growth parameter (p). The GGM can reproduce a range of growth dynamics, including constant growth (p=0), sub-exponential or polynomial growth (0 In this equation the numerator represents the new cases *I**i*, and the denominator represents the total number of cases that contribute to the new cases *I**i* at time *t**i*. Parameter 0≤*α*≤1 represents the relative contribution of imported cases to the secondary disease transmission. We perform a sensitivity analyses by setting *α* = 0.15 and *α* =1.0(39). Next, in order to derive the uncertainty bounds around the curve of *R*t,we estimate *R*t for 300 simulated curves assuming a Poisson error structure (40). ### Reproduction number (Rt) from the analysis of cluster sizes In the branching process analysis, the number of transmissions caused by each new infection is modeled as a negative binomial distribution. This is parameterized by the effective reproduction number, Rt, and the dispersion parameter, k. The reproduction number provides the average number of secondary cases per index case, and the dispersion parameter varies inversely with the heterogeneity of the infectious disease. This parameterization allows variation in the both the strength and heterogeneity of disease transmission, where a lower dispersion parameter indicates higher transmission heterogeneity. Assuming each cluster of infection is initiated by a single infection, branching process theory provides an analytic representation of the size distribution of cluster sizes (41). This permits direct inference of the maximum likelihood estimate and confidence interval for Rt and k. ## Results ### Incidence data and reporting delays The COVID-19 epidemic curve by the date of reporting and stratified for local and imported incidence case counts is shown in Figure 1. It is worth noting that the majority of the imported cases occurred at the beginning of the outbreak in Singapore. The cumulative curve shows an upward trajectory with an average of 3.4 new cases reported per day (Figure 1). Out of 22 imported cases, only 4 cases have been linked to secondary cases. Meanwhile, a total of 61 autochthonous cases have been reported as of February 19, 2020 including 8 cases that remain unlinked to any known transmission chains. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F1) Figure 1: Unadjusted local and imported Incidence cases by date of reporting as of February 19, 2020. The solid line represents the cumulative cases count for the COVID-19 cases in Singapore by the date of reporting. The reporting-delay adjusted epidemic curve of local cases by date of symptoms onset is shown in Figure 2, and the distribution for reporting delays is displayed in Figure 3. The gamma distribution provided a reasonable fit to the distribution of overall reporting delays and estimated a mean reporting delay at 6.9 days (95% CI: 5.9, 7.9) (Figure 3). We also found that imported cases tend to have shorter mean reporting delays compared to local cases (mean 3.5 vs. 8.3 days, Wilcoxon test, p<0.001) as imported cases tend to be identified more quickly. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F2) Figure 2: Reporting delay adjusted incidence for the COVID-19 outbreak in Singapore accounting for local cases as of February 19, 2020. Blue bars represent the raw incidence, red solid line represents the adjusted incidence, red dotted lines represent the 95% lower and upper bound of the adjusted incidence. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F3) Figure 3: Distribution of reporting delays for all cases as of February 19, 2020 represented as a histogram. The red circle represents the mean of reporting delays and the horizontal line represents the 95% CI assuming reporting delays have a gamma distribution. ### Transmission clusters As of February 19 2020, nine different clusters of COVID-19 cases with 1-22 cases per cluster have been reported in Singapore. A schematic diagram and characteristics of the COVID-19 clusters in Singapore are given in Figure 4 and Table 1. The geographic location of the five largest clusters accounting for 53.6% of the total cases is shown in Figure 5 whereas the corresponding distribution of cluster sizes is shown in Figure 6. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F4) Figure 4: Cluster network of the cases in Singapore for the COVID-19 global epidemic as of February 19, 2020. The pink circles represent the imported cases where as the blue circles represent cases with no travel history to China. The larger dotted circles represent the COVID-19 disease clusters. Each blue arrow represents the direction in which the disease was transmitted. Pink arrows represent immediate family. Dates below the circles are the dates of case reporting. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F5) Figure 5: Map depicting the spatial distribution of the 5 largest COVID-19 clusters in Singapore; Grand Hayatt cluster, Yong Thai Hang cluster, Grace Assembly of God Church, Seletar Aerospace cluster and The Life Church and Missions cluster as of February 19, 2020. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F6) Figure 6: Distribution of COVID-19 cluster sizes in Singapore as of February 19, 2020. ### Yong Thai Hang cluster This cluster with 9 cases was the first to be reported in Singapore. It has nine traceable links, including eight Chinese and one Indonesian national associated with the visit on January 23, 2020 of Chinese tourists to the Yong Thai Hang health products store, a shop that primarily serves the Chinese population. Four shop employees and the tour guide were first identified as a cluster on February 4, 2020 (12, 42, 43). The tour guide subsequently infected her husband, a newborn and the domestic helper (43). No further cases have been added to this cluster as of February 8, 2020. ### Grand Hayatt hotel This cluster with 3 local cases was the second cluster to receive international attention, as it originated from a business meeting held at the Grand Hayatt hotel (44). Four international cases associated with this cluster had left Singapore before the local cluster was identified. All Singaporean residents associated with this cluster have recovered as of February 19, 2020 (44). No additional cases have been added to this cluster as of February 8, 2020. ### The Life Church and Missions cluster This cluster comprising of 6 cases was associated with The Life Church and Missions in Paya Lebar on February 7, 2020. This cluster is composed of two imported cases who are related to each other and four local cases. No further cases have been added to this cluster as of February 19, 2020. ### Seletar Aerospace Heights cluster This cluster with 5 Bangladeshi work pass holders was identified on February 9, 2020. No further cases have been added to this cluster as of February 15, 2020. ### Grace Assembly of God Church cluster The fifth and largest cluster composed of 22 local cases has been linked to the two branches of Grace Assembly of God Church at Tanglin and Bukit Batok, which was identified on February 12, 2020 (42, 45). This church has about 4800 people in attendance on average over the weekend. While the church has momentarily closed, field investigations have not led to conclusive evidence regarding super-spreading transmission (46). No further cases have been added to this cluster as of February 19, 2020. ### Cluster A The first case of this cluster comprising 2 imported and one local case was identified on January 30, 2020. The two imported cases are related to each other. No further cases have been added in this cluster as of February 10, 2020. ### Cluster B The first case of this cluster composed of 3 local cases was identified on February 14, 2020. No further cases have been added in this cluster since February 18, 2020. ### Cluster C This first case of this cluster composed of 4 local cases was identified on February 12, 2020. No further cases have been added to this cluster since February 17, 2020. ### Cluster D The first case of this cluster composed of two local cases who are related to each other was identified on February 9, 2020. No further cases have been added in this cluster since February 16, 2020. ## Reproduction Numbers Using the delay-adjusted incidence curve of local cases, we estimated the intrinsic growth rate r at 0.9 (95% CI: 0.4,1.7) and the scaling of growth parameter p at 0.4 (95%CI: 0.2,0.6), indicating sub-exponential growth dynamics of COVID-19 (Figure 7). Our R estimate was not sensitive to changes in parameter *α*. Specifically, our mean effective reproduction number Rt was estimated at 1.1 (95% CI:1.1, 1.3) when *α* =0.15 and 1.1 (95% CI:1.1, 1.2) when *α* =1(Figure 7). It is worth noting that ignoring the presence of imported cases, our estimates would have been slightly underestimated (by about 0.1 in absolute value). ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/25/2020.02.21.20026435.1/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/02/25/2020.02.21.20026435.1/F7) Figure 7: Individual level reproduction number with 95% CI estimated by adjusting for the imported cases *α* =0.15. Estimates for growth rate (r) and the scaling of the growth rate parameter (p) are also provided. The plot at the bottom depicts the fit of the Generalized Growth Model to the Singapore adjusted incidence data assuming Poisson error structure. Based on the entire distribution of cluster sizes, we jointly estimated the reproduction number R and the dispersion parameter k. Fitting the negative binomial distribution to the cluster data in the empirical distributions of the realizations during the early stages of the outbreak in Singapore, the reproduction number is estimated at 0.6 (95% CI: 0.4, 1.1) after adjusting for the independent and identically distributed observation probability of 50%. The dispersion parameter is estimated at 0.4 (95% CI: 0.1, Inf.) indicating higher transmission heterogeneity. ## Discussion Overall, current estimates of transmission potential in Singapore, based on two different data sources and different methods, suggest the potential for further local transmission of SARS-CoV-2 in Singapore. While large-scale transmission has not been reported, the fact that asymptomatic and subclinical cases are now well documented for COVID-19 (47) suggests that our estimates may be underestimated (48). On the other hand, it is not clear if asymptomatic or subclinical cases are as infectious as symptomatic cases. Indeed, we have reported that multiple local cases have yet to be traced to existing transmission chains. Additional data collected during the course of the outbreak will help obtain an improved picture of the transmission dynamics (8). These findings emphasize the need to strengthen public health interventions including active case contact tracing activities in countries with emerging transmission of SARS-CoV-2. It is worth noting that current data suggests that imported cases have contributed little to secondary cases in Singapore, with most of the imported cases dating back to the early phase of the epidemic. However, there are examples such as the Grand Hyatt Singapore conference cluster and the Yong Thai Hang cluster that were linked to imported sources, and the original sources had left Singapore before these local clusters emerged (44, 49). Our Rt estimates for Singapore are substantially lower than mean estimates reported for the COVID-19 epidemic in China, which range from: 2.2-7.1 (50-53). This indicates that containment efforts are having a significant impact in Singapore. However, some differences in the reproduction numbers reported for the epidemic in China may result from different methods, differences in data sources, and time periods used to estimate the reproduction number. Similarly, a recent study has shown an average reporting delay of 6.1 days in China (54) which agrees with our mean estimate for local cases (6.9 days). Moreover, the scaling parameter for growth rate (p) indicates a sub-exponential growth pattern in Singapore, reflecting the effective isolation and control strategies in the region. This is consistent with a sub-exponential growth pattern for Chinese provinces excluding Hubei (p=0.67), as estimated by a recent study (55); whereas, an exponential growth pattern was estimated for Hubei (p∼1.0) (55). A previous study (56) on the 2015 MERS outbreak in South Korea reported a substantial potential for superspreading despite a subcritical Rt. The lower estimate of the dispersion parameter in our study also indicate significant transmission heterogeneity in Singapore. Super-spreading events of MERS-CoV and SARS-CoV associated with nosocomial outbreaks have been observed in the past and attributed to substantial diagnostic delays (16, 30). Although the average time from onset of symptoms to diagnosis for COVID-19 patients in Singapore is 7.9 days and no super-spreading event has been observed yet, the dispersion parameter, k<1, indicates the probability of observing large clusters and the potential for super-spreading (56, 57). Therefore, public health measures enacted by public health authorities in Singapore that advise the public to avoid mass gatherings and confined places are crucial to prevent disease amplification events. However, the presence of asymptomatic cases in the community represent an ongoing threat (58, 59) although it is not currently known if subclinical cases are less infectious. This highlights the need for rapid testing suspected cases for rapid detection and isolation of confirmed cases. To achieve this goal, public health authorities in Singapore are reactivating 900 GP clinics (60). Investigations are also underway to identify additional cases linked to the “Grace Assembly of God” cluster in Singapore. While most clusters have stabilized (no recent additional cases in most clusters), this cluster is still being consolidated (49, 61). Beyond Singapore, COVID-19 clusters have been reported in Germany and on the Diamond Cruise ship in Japan (62-65). Moreover, Singapore has also produced secondary chains of disease transmission beyond its borders (22). Although Singapore has been detecting and isolating cases with diligence, our findings from the analysis underscore the need for continued and sustained containment efforts to prevent large-scale community transmission including nosocomial outbreaks. Overall, the current situation in Singapore highlights the need to investigate the unlinked and asymptomatic cases that could be a potential source of secondary cases. Although Singapore has a highly efficient contact tracing mechanism in place (22), continued epidemiological investigations and active case finding efforts are needed to control the outbreak. Our study is not exempt from limitations. First, the outbreak is still ongoing and we continue to monitor the transmission potential of COVID-19 in Singapore. Second, onset dates are missing in nine cases, which were excluded from our analyses. Third, we cannot rule out that additional cases will be added to existing clusters, which may lead to underestimating the reproduction number based on the cluster size distribution. Lastly, the cluster size analysis is modeled under the assumption that there is one primary case (i.e. all cases can be connected eventually to one person), rather than being multiple imported cases. Additional work could be conducted to account for the uncertainty regarding the primary number of cases per cluster and the possibility of ongoing transmission in some clusters. ## Conclusion This is the first real-time study to estimate the evolving transmission potential of SARS-CoV-2 in Singapore. Our current findings point to ongoing transmission of SARS-CoV-2, although the reproduction number appears to be slightly above the epidemic threshold of 1.0. These estimates underscore positive impact of containment efforts in Singapore while at the same time suggest the need to enhance social distancing and active case finding efforts to stomp out all active chains of transmission. ## Data Availability All data are publicly available. ## List of abbreviations COVID-19 SARS-CoV-2 ## Ethics approval and consent to participate Not applicable ## Consent for publication Not applicable ## Conflict of Interest The authors declare no conflicts of interest. ## Funding G.C. is supported by NSF grants 1610429 and 1633381. ## Data declaration All data are publicly available. ## Author Contributions A.T, S.B., P.Y. and G.C. analyzed the data. A.T., Y. L, P.Y and S.M. retrieved and managed data; A.T and G.C wrote the first draft of the manuscript. All authors contributed to writing and revising subsequent versions of the manuscript. All authors read and approved the final manuscript. * Received February 21, 2020. * Revision received February 25, 2020. * Accepted February 25, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.WHO, W.H.O., Novel Coronavirus (2019-nCoV), Situation report - 1, 21 January 2020. 2020: WHO. 2. 2.Worldometer. COVID-19 coronavirus / cases. 2020 [cited 2020 February 19,2020]; Available from: [https://www.worldometers.info/coronavirus/coronavirus-cases/](https://www.worldometers.info/coronavirus/coronavirus-cases/). 3. 3.Nishiura, H., et al., The Extent of Transmission of Novel Coronavirus in Wuhan, China, 2020. Journal of Clinical Medicine, 2020. 9(2): p. 330. 4. 4.Tian, H., et al., Early evaluation of Wuhan City travel restrictions in response to the 2019 novel coronavirus outbreak. medRxiv, 2020: p. 2020.01.30.20019844. 5. 5.Olowski, L.J., Lew Olowski: Coronavirus worse than reported – here’s how China is catching up to reality, in Fox News. 2020: Fox News. 6. 6.Kottasová, I., China’s coronavirus numbers highlight the challenges of an evolving epidemic, in CNN. 2020: CNN. 7. 7.WHO, W.H.O., Coronavirus disease 2019 (COVID-19) Situation Report – 27. 2020: World Health Organization. 8. 8.Singapore, M.o.H., One more case discharged; five new cases of COVID-19 infection confirmed. 2020: Ministry of Health, Singapore. 9. 9.Singapore, M.o.H., Confirmed imported case of novel coronavirus infection in Singapore; multi-ministry taskforce ramps up precautionary measures. 2020: Ministry of Health Singapore. 10. 10.Ang, D. The Big Story: 92 Singaporeans evacuated from Wuhan; Hikikomori in Singapore. 2020 [cited 2020 February 13]; Available from: [https://www.straitstimes.com/multimedia/the-big-story-92-singaporeans-evacuated-from-wuhan-hikikomori-in-singapore](https://www.straitstimes.com/multimedia/the-big-story-92-singaporeans-evacuated-from-wuhan-hikikomori-in-singapore). 11. 11.Tan, A., Coronavirus: Flight with 174 Singaporeans evacuated from Wuhan lands at Changi Airport on Sunday morning, in The Straits Times. 2020: ST. 12. 12.Singapore, M.o.H., Confirmed cases of local transmission of novel coronavirus infection in Singapore. 2020: Minstry of Health Singapore. 13. 13.Singapore, M.o.H., Five more cases discharged, three new cases of COVID-19 infection confirmed. 2020: Ministry of Health. 14. 14.Gan, N., N. Thomas, and D. Culver, Hundreds of frontline medics likely infected with coronavirus in China presenting new crisis for the government, in CNN. 2020: CNN. 15. 15.Boseley, S., More than 1,700 health workers infected by coronavirus in China, in The Guardian. 2020: The Guardian. 16. 16.Chowell, G., et al., Transmission characteristics of MERS and SARS in the healthcare setting: a comparative study. BMC Medicine, 2015. 13(1): p. 210. 17. 17.Gopalakrishna, G., et al., SARS transmission and hospital containment. Emerging infectious diseases, 2004. 10(3): p. 395–400. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15109403&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F25%2F2020.02.21.20026435.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220079400002&link_type=ISI) 18. 18.Lim, J., Coronavirus: Doctors, staff and patients to restrict movements to within one hospital, in The Straits Times Singapore. 2020: SPH Website. 19. 19.Aravindan, A. and J. Geddie, Singapore lifts virus alert to SARS level, sparking panic buying, in Rueters. 2020: Rueters. 20. 20.Sim, D., Coronavirus: why did Singapore have more cases than Hong Kong – until now?, in This Week in Asia. 2020: Health and Environment. 21. 21.Rettner, R., Can the new coronavirus spread through building pipes?, in Live Science. 2020: Live Science. 22. 22.Vaswani, K., Coronavirus: Why Singapore is so vulnerable to coronavirus spread, in BBC News. 2020: BBC. 23. 23.Diekmann, O. and J.A.P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation 1st ed. 2000: Wiley. 24. 24.Anderson, R.M. and R.M. May, Infectious Diseases of Humans, ed. Oxford. 1991, Oxford Univeristy Press. 25. 25.Singapore, M.o.H., NEWS HIGHLIGHTS, in Press Releases. 2020. 26. 26.Yan, P. and G. Chowell, Quantitative methods for infectious disease outbreak investigations. 2019: Springer. 27. 27.Yan, P. Estimation for the infection curves for the spread of Severe Acute Respiratory Syndrome (SARS) from a back-calculation approach. 2018 [cited 2019 February 8]; Available from: [http://archive.dimacs.rutgers.edu/Workshops/Modeling/slides/Yan.pdf](http://archive.dimacs.rutgers.edu/Workshops/Modeling/slides/Yan.pdf). 28. 28.Lawless, J.F., Adjustments for reporting delays and the prediction of occurred but not reported events. 1994. 22(1): p. 15–31. 29. 29. G.C. Taylor, Claims Reserving In Non Life Insurance. Insurance Series. 1985: Elsevier. 30. 30.Tariq, A., et al., Assessing reporting delays and the effective reproduction number: The Ebola epidemic in DRC, May 2018-January 2019. Epidemics, 2019. 26: p. 128–133. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00046318&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F25%2F2020.02.21.20026435.1.atom) 31. 31.Kalbfleisch, J.D. and J.F. Lawless, Estimating the incubation time distribution and expected number of cases of transfusion-associated acquired immune deficiency syndrome. Transfusion, 1989. 29(8): p. 672–6. 32. 32.Nishiura, H., et al., The ideal reporting interval for an epidemic to objectively interpret the epidemiological time course. Journal of the Royal Society, Interface, 2010. 7(43): p. 297–307. 33. 33.Viboud, C., L. Simonsen, and G. Chowell, A generalized-growth model to characterize the early ascending phase of infectious disease outbreaks. Epidemics, 2016. 15: p. 27–37. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2016.01.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27266847&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F25%2F2020.02.21.20026435.1.atom) 34. 34.Nishiura, H., N.M. Linton, and A.R. Akhmetzhanov, Serial interval of novel coronavirus (2019-nCoV) infections. medRxiv, 2020: p. 2020.02.03.20019497. 35. 35.You, C., et al., Estimation of the Time-Varying Reproduction Number of COVID-19 Outbreak in China. medRxiv, 2020: p. 2020.02.08.20021253. 36. 36.Nishiura, H. and G. Chowell, Early transmission dynamics of Ebola virus disease (EVD), West Africa, March to August 2014. Euro Surveill, 2014. 19(36). 37. 37.H. Nishiura, G.C., The Effective Reproduction Number as a Prelude to Statistical Estimation of Time-Dependent Epidemic Trends. Mathematical and Statistical Estimation Approaches in Epidemiology, ed. D. Springer. 2009. 103–112. 38. 38.Paine, S., et al., Transmissibility of 2009 pandemic influenza A(H1N1) in New Zealand: effective reproduction number and influence of age, ethnicity and importations. Euro Surveill, 2010. 15(24). 39. 39.Nishiura, H. and M.G. Roberts, Estimation of the reproduction number for 2009 pandemic influenza A(H1N1) in the presence of imported cases. Eurosurveillance, 2010. 15(29): p. 19622. 40. 40.Chowell, G., Fitting dynamic models to epidemic outbreaks with quantified uncertainty: A primer for parameter uncertainty, identifiability, and forecasts. Infectious Disease Modelling, 2017. 2(3): p. 379–398. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.idm.2017.08.001&link_type=DOI) 41. 41.Blumberg, S. and J.O. Lloyd-Smith, Inference of R(0) and transmission heterogeneity from the size distribution of stuttering chains. PLoS computational biology, 2013. 9(5): p. e1002993–e1002993. 42. 42.Khalik, S., Coronavirus in Singapore: Of the 50 cases so far, local cases exceed imported ones, in The Straits Times. 2020: ST Singapore. 43. 43.CNA, Coronavirus cases in Singapore: Trends, clusters and key numbers to watch, in CNA. 2020. 44. 44.Wong, L., Coronavirus: Full recovery for all three local cases linked to Grand Hyatt event, in The Straits Times. 2020: The Straits Times. 45. 45.MalayMail, COVID-19: SAF regular among three new victims in Singapore as Grace Assembly of God cluster grows to 18 cases, in Malaymail. 2020. 46. 46.Lee, J., COVID-19 outbreak: Grace Assembly of God cluster grows to 21 cases, in Mothership. 2020: Mothership. 47. 47.Team, T.N.C.P.E.R.E., The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China. 2020. 41(2): p. 145–151. 48. 48.Belluz, J., Why the coronavirus outbreak might be much bigger than we know, in Vox. 2020: Vox. 49. 49.Kurohi, R., Coronavirus: 3 new cases in Singapore, including one first warded as dengue patient; five more patients discharged, in The Straits Times. The Straits Times. 50. 50.Shen, M., et al., Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China. bioRxiv, 2020: p. 2020.01.23.916726. 51. 51.Mizumoto, K., K. Kagaya, and G. Chowell, Early epidemiological assessment of the transmission potential and virulence of 2019 Novel Coronavirus in Wuhan City: China, 2019-2020. medRxiv, 2020: p. 2020.02.12.20022434. 52. 52.Zhao, S., et al., Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak. International Journal of Infectious Diseases, 2020. 92: p. 214–217. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F25%2F2020.02.21.20026435.1.atom) 53. 53.Liu, T., et al., Transmission dynamics of 2019 novel coronavirus (2019-nCoV). bioRxiv, 2020: p. 2020.01.25.919787. 54. 54.Kucharski, A.J., et al., Early dynamics of transmission and control of 2019-nCoV: a mathematical modelling study. medRxiv, 2020: p. 2020.01.31.20019901. 55. 55.Roosa, K., et al., Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020. Infectious Disease Modelling, 2020. 5: p. 256–263. 56. 56.Kucharski, A.J. and C.L. Althaus, The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission. Eurosurveillance, 2015. 20(25): p. 21167. 57. 57.Lloyd-Smith, J.O., et al., Superspreading and the effect of individual variation on disease emergence. Nature, 2005. 438(7066): p. 355–359. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F25%2F2020.02.21.20026435.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233300200048&link_type=ISI) 58. 58.Sim, D. and K. Xinghui, Singapore: nine new coronavirus cases, no plans to raise outbreak alert to red, in This Week in Asia. 2020. 59. 59.Baker, J.A. and A.H. Mamud, 9 new COVID-19 cases in Singapore, including 6 linked to Grace Assembly of God cluster, in Channel New Asia. 2020: CNA. 60. 60.Chang, N. and M. Yong, Public Health Preparedness Clinics reactivated to reduce risk of COVID-19 spread, in CNA. 2020: CNA. 61. 61.Hui, M., A cluster of coronavirus cases in South Korea has been traced to a cult, in Quartz. 2020: Quartz. 62. 62.Marcus, I., Chronology: Germany and COVID 19 (‘Coronavirus’), in The Berlin Spectator. 2020: The Berlin Spectator. 63. 63.Rueters, Germany Confirms Seventh Coronavirus Case, in Rueters. 2020. 64. 64.McCurry, J., Japan reports first coronavirus death as 44 more cases confirmed on cruise ship, in The Guardian. 2020: The Guardian. 65. 65.Woods, A., 44 new coronavirus cases reported on quarantined Diamond Princess cruise, in Newyork Post. 2020. [1]: /embed/graphic-3.gif