Rapid transmission of coronavirus disease 2019 within a religious sect in South Korea: a mathematical modeling study ==================================================================================================================== * Jong-Hoon Kim * Hyojung Lee * Yong Sul Won * Woo-Sik Son * Justin Im ## Abstract Rapid transmission of coronavirus disease 2019 (COVID-19) was observed in the Shincheonji Church of Jesus, a religious sect in South Korea. The index case was confirmed on February 18, 2020 in Daegu City, and within two weeks, 3,081 connected cases were identified. Doubling times during these initial stages (i.e., February 18 – March 2) of the outbreak were less than 2 days. A stochastic model fitted to the time series of confirmed cases suggests that the basic reproduction number (*R*) of COVID-19 was 8.5 [95% credible interval (CrI): 6.3, 10.9] among the church members, whereas (*R* = 1.9 [95% CrI: 0.4, 4.4]) in the rest of the population of Daegu City. The model also suggests that there were already 4 [95% CrI: 2, 11] undetected cases of COVID-19 on February 7 when the index case reportedly presented symptoms. The Shincheonji Church cluster is likely to be emblematic of other outbreak-prone populations where *R* of COVID-19 is higher. Understanding and subsequently limiting the risk of transmission in such high-risk places is key to effective control. **Highlights** * Basic reproduction number (*R*) of COVID-19 in a religious community of Shincheonji Church of Jesus was estimated to be 8.5 [95% credible interval (CrI): 6.3, 10.9], which is more than 4 times larger than the general population (*R* = 1.9 [95% CrI: 0.4, 4.4]) * There were estimated 4 [95% CrI: 2, 11] undetected cases when the index case from the religious community reported symptom on February 7. * The Shincheonji Church cluster is likely to be emblematic of other outbreak-prone populations where *R* of COVID-19 is higher. Understanding and subsequently limiting the risk of transmission in such high-risk places is key to effective control. Keywords * COVID-19 * Korea * Shincheonji Church * transmission model * reproduction number ## 1. Introduction Coronavirus disease 2019 (COVID-19) has become a global pandemic since it was first reported in Wuhan, China in December 2019 with the name of novel coronavirus disease (Li et al., 2020a). The causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), transmits mainly through human-to-human contact (Chan et al., 2020), which can happen even during the infector is asymptomatic (Rothe et al., 2020; Yu et al., 2020). Infection with the virus causes diseases with varying degree of symptoms including death (Fu et al., 2020; Guan et al., 2020). Infection mortality ratio is lowest among children aged between 5 and 9 years and increases loglinearly with age (O’Driscoll et al., 2021). One key characteristic of COVID-19 pandemic is that transmission events in high-risk settings such as super-spreading events (SSEs) contribute to most transmissions (Adam et al., 2020; Lemieux et al., 2020). The risk of COVID-19 transmission is believed to high in places with high occupancy and poor ventilation (Jones et al., 2020). One extreme example is the outbreak in the Diamond Princess cruise ship, where 17% (619/3711) of the passengers were infected from January 25 to February 20, 2020 (Russell et al., 2020). Other examples include transmission events in bars and wedding (Adam et al., 2020) in Hong Kong, nursing homes in U.S. (Chen et al., 2021), telemarketers working in group in closed places (Park et al., 2020) and fitness classes (Jang et al., 2020) in South Korea, and also religious gatherings, which we describe below. Explosive spread of COVID-19 was observed in the Shincheonji Church of Jesus (“Shincheonji”), a religious sect in South Korea. The index case was confirmed on February 18, 2020 and within two weeks, 3,081 connected cases were identified (Korea Disease Control and Prevention Agency, 2020a). A simple calculation reveals that the outbreak size doubled in less than every 2 days (14/log2(3081) ≈ 1.21), which is smaller than doubling times reported in the early stages of the COVID-19 outbreak in China (2.5 and 3.1 days in Hubei Province and Hunan Province, respectively) (Muniz-Rodriguez et al., 2020b), Spain (2.8 days) (Guirao, 2020), the US (2.7 days) (Lurie et al., 2020), and Korea (2.8 – 10.2 days) (Shim et al., 2020b). A total of 6,684 confirmed cases were reported in Daegu City as of March 31, 2020 of which 4,467 (66.8%) were Shincheonji members, representing close to half (47.9%, 4,467/9,334) of the city’s total Shincheonji membership. Previous studies highlighted that COVID-19 transmissions involve SSEs (Xu et al., 2020), which can play a key role in sustained community transmissions (Adam et al., 2020; Lemieux et al., 2020). However, there have been no attempts to model the dynamics of COVID-19 transmission within high-risk settings and their interaction with the general community. In this study, we modeled the outbreak in the Shincheonji community while accounting for its interaction with the rest of population. We used the stochastic model to account for the stochastic nature of the transmission events. We used the model to explore the differences in the basic reproduction number (*R*) between the high-risk setting and the general community, and quantify uncertainties related to the initial conditions and dynamics of transmission under the dynamic intervention programs. ## 2. Materials and Methods ### 2.1 Backgound on the Shincheonji Church of Jesus Shincheonji was founded by Man-hee Lee in 1984 and has approximately 245,000 members including 30,000 foreigners (Chung and Hill, 2020). At Shincheonji gatherings, worshipers used to sit close together on the floor and facial coverings, such as glasses and face masks, are forbidden. Members were expected to attend services despite illness (Choe, 2020). The index case of the Daegu City outbreak was identified as a Shincheonji member and some 1,000 people were reported to have attended worship together (Yonhap, 2020). Further tracing of church members identified clustering in apartment complexes. Of 142 residents in a particular Daegu apartment block, 94 (66%) were Shincheonji members of whom 46 (38.9%) tested positive for the virus (Myung, 2020). ### 2.2 Data Time series of patients confirmed with COVID-19 in Shincheonji community and the overall Daegu City over the period of 19 February – 31 March 2020 was compiled based on the daily reports from Korea Disease Control and Prevention Agency (KDCA) (Korea Disease Control and Prevention Agency, 2020a) (Figure 1). The reports provide the number of cases confirmed with SARS-CoV 2 based on reverse transcription polymerase chain reaction (RT-PCR) by category (Shincheonji or non-Shincheonji). We made some adjustments to the existing data before we fit the model. First, in the beginning of the outbreak, KDCA provided both daily and cumulative numbers of cases confirmed for SARS-CoV 2, which did not agree always. If there is a discrepancy between these numbers, we prioritized cumulative numbers as this figure was reported continuously throughout the outbreak. Second, data were missing for some days for the number cases for Shincheonji members. We imputed missing values using the cubic spline method (Figure S1 in the Supplementary Material). ### 2.3. Doubling time The epidemic doubling time (*T**d*) represents the duration in which the cumulative incidence doubles. Assuming exponential growth with a constant epidemic growth rate (*r*), the epidemic doubling time can be calculated by the following equation (Anderson et al., 2020; Lurie et al., 2020; Muniz-Rodriguez et al., 2020a) ![Formula][1] Epidemic growth rate (*r*) may be estimated based on the data. For example, *r*(*t*) can be estimated by the following equation: ![Formula][2] where *C*(*t*) indicates the cumulative number of infected people at time *t* and *Δt* is the duration over which *r*(*t*) is assumed to be constant. *r*(*t*) can be calculated over the fixed time interval (e.g., 1 day or 1 week) (Ebell and Bagwell-Adams, 2020; Patel and Patel, 2020) or variable time intervals (e.g., days on which the number of cases doubles, quadruples, etc.) (Muniz-Rodriguez et al., 2020a; Shim et al., 2020b). We calculated doubling based on prior 7 days or 1 day from 18 February to 5 March 2020, when the epidemic peaked and no further doubling of cumulative number of cases occurred onward. The basic reproduction number (*R*) is defined as the average number of secondary cases caused by a single infected case in an entirely susceptible population and it provides sufficient information to produce doubling times in the beginning of an outbreak. However, estimating *R* requires additional information such as generation time or developing a mechanistic model, and its estimates come with higher degree of uncertainty (Anderson et al., 2020). Calculating doubling times requires fewer assumptions and also allows us to compare our results with estimates from different settings where doubling times, but not reproduction numbers, are available. ### 2.4. Mechanistic model of COVID-19 transmission We developed a stochastic model of COVID-19 transmission within the Shincheonji community and the overall population of Daegu City. The model includes six disease states: susceptible (*S*), exposed but not infectious (*E*), pre-symptomatic but infectious (*P*), symptomatic and infectious (*I*), asymptomatic but infectious (*A*), confirmed and isolated (*C*), and recovered (*R*). The model includes two patches to model Shincheonji and non-Shincheonji people, separately. Transmission rates may differ for each patch and person from one patch may infect people from the other patch. (Figure 2). This modeling framework of mixing between two distinct sub-populations has been adopted in previous works, ranging from sexually transmitted diseases (Koopman et al., 1988) to vector-borne diseases such as dengue (Lee and Castillo-Chavez, 2015), where formulations for mixing between patches vary. We adopted the formulation used in the work on modeling transmission of cholera between hotspot and non-hotspot areas (Azman and Lessler, 2015). Mixing between two sub-populations are defined by the 2 × 2 contact matrix, ![Formula][3] where *c**ij* indicates the fraction of time that individuals from patch *i* spends in patch *j*. Next, we impose two conditions on the matrix *C*: 1. Individuals must reside in either of the two patches, i.e., *c**i*1 + *c**i*2 = 1 for *i* =1 (Shincheonji) and 2 (non-Shincheonji). 2. The population in each patch remains constant, i.e., *c*12*N*1 = *c*21*N*2, where *N*1 and *N*2 represent population size for patch 1 and 2, respectively. The above conditions may transform the contact matrix *C* to the following form: ![Formula][4] containing only one unknown parameter *c*12. The force of infection for individuals from patch *i* at time *t, λ**i*(*t*), is defined as follows: ![Formula][5] where *β**j* indicates local transmission rate in patch *j* and *I**k*(*t*) indicates number of infectious individuals from patch *k*. The transitions between states are modeled using an explicit tau-leap algorithm (Gillespie, 2001) to account for stochasticity of the infection transmission process. The number of susceptible people in patch *i* at time *t* + Δ*t, S**i*(*t* + Δ*t*), is written as follows: ![Formula][6] ![Graphic][7] represents the number of people who transit from state *S* to state *E* from *t* to *t* + Δ*t* in patch *i* and is a random variable with binomial distribution: ![Formula][8] That is, it is represented as an integer varying between 0 and *S**i*(*t*). For states from which more than one potential transition exist (e.g., *P* to either *A* or *I*), multinomial distributions were applied. For instance, the number of people transit from *P* to either *I* or *A* are given as follows: ![Formula][9] where *π* is a vector given as ![Formula][10] The first element of *π* indicates a probability of transition from *P* to *I* and the second element indicates the probability of transition from *P* to *A*. The number of people in other states (i.e., *E, A, I, C, R*) at time t can be described similarly. The model was implemented in a combination of R and C++ languages, in which the core transmission model part is expensive and was written in C++. All the computer codes that generate the results in this paper are available at the author’s GitHub repository (Kim, 2021). ### 2.5. Modeling intervention program To account for intensification of the intervention such as case isolation and contact tracing with subsequent testing during the outbreak, we assumed case isolation rate (1 / mean time between symptom onset and case isolation) and transmission rate of the infectious people per unit time change over time. Specifically, we assumed that the case isolation rate, α(*t*), starts increasing on February 20 from the initial value of *α*init when 4,474 out of 9,334 Shincheonji members were identified and were asked to self-isolate. During model fitting, we let data suggest the duration of intervention, *d* in day, which is the time required for the case isolation rate to reach its minimum, *α*(*t*) = *α*final for *t* > February 20 + *d*. We assumed that the mean time between symptom onset and case isolation linearly decreases over the intervention period *d*. In other words, α(*t*) is formulated as follows: ![Formula][11] where *α*final is assumed to be 1 day based on the experiences in Busan City in Korea and α(*t*) is assumed to be zero before February 17 when the index case was detected. Similarly, transmission rate per unit time at time *t, β**i*(*t*) for *i* = 1 (Shincheonji members), 2 (non-Shincheonji people in Daegu City), is assumed to linearly decrease during the intervention period. ![Formula][12] Here, ![Graphic][13] and ![Graphic][14] indicate the transmission rate per unit time before the intervention and after the intervention measures fully take effect, respectively. They can be derived once *R*0,*i* and *R*final are given as: ![Formula][15] and ![Formula][16] ### 2.6. Parameter estimation Our model of COVID-19 transmission requires 15 parameters (Table 1). We divided the model parameters into three classes depending on our belief on their relative certainty. The first class includes parameters related to the natural history of infection and population size and we deemed that available parameter estimates are reliable. For these parameters, we used their point estimates based on analyses of data on COVID-19 transmissions in Korea or China. For the second class, which includes parameters related to intervention programs, we used our best guesses based on supporting evidence but still acknowledged their uncertainty. Therefore, we analyzed the models under various assumptions on their values within some pre-specified ranges. Finally, we defined six parameters that are critical for characterizing dynamics of COVID-19 transmission in Shincheonji members and non-Shincheonji people. We estimated these parameters by fitting the model to daily confirmed COVID-19 cases of Shincheonji members and non-Shincheonji people. Estimation of parameters θ = (*R*0,1, *R*0,2, *I*, *c*12, *d, R*final) was based on Approximate Bayesian Computation Sequential Monte Carlo (ABC-SMC) (Minter and Retkute, 2019). The ABC is a method for approximating posterior distributions given data *D, p*(θ|*D*), by accepting proposed parameter values when the difference between simulated data *D** and *D, d*(*D, D**), is smaller than tolerance ϵ: ![Formula][17] For our model, *d*(*D, D**) is defined as the sum of the squared differences in daily confirmed cases over the outbreak of duration *T* days, that is, ![Formula][18] for Shincheonji (*i* = 1) and non-Shincheonji (*i* = 2). Here, *D**t* and ![Graphic][19] represent observed daily confirmed cases and model predicted values at time day *t*, respectively. ABC-SMC was designed to increase efficiency of the ABC method and ABC is applied in a sequential manner by constructing intermediate distributions, which converge to the posterior distribution. Tolerance ϵ is gradually decreased and each intermediate distribution is obtained as a sample that is drawn with weights from the previous distribution and then perturbed through a kernel *K*(θ|θ*). The kernel helps keep the algorithm from being stuck in local optimum while maintaining the efficiency of the ABC-SMC method. Minimally informative uniform distributions were used as prior distributions and estimation procedure was repeated for ten different random seeds. The resulting distribution was summarized as median, 50% credible intervals (CrI; interval between 25% and 75% percentiles) and 95% CrI (interval between 2.5% and 97.5% percentiles). More details of the algorithm such as prior distribution for each parameter, the number of steps, the tolerance values for each step, perturbation kernel appear in the Supplementary Material. ## 3. Results ### 3.1 Doubling time Over the period of February 18 – March 5, during which doubling of confirmed cases occurred 12 times, doubling times were <1 day in the beginning and increased subsequently with daily doubling time presenting higher variability for both Shincheonji and non-Shincheonji values (Table 2). Doubling times calculated over sliding one-week intervals remained shorter than 3 days for the most part for both Sincheonji and non-Shincheonji population. ### 3.2 Comparison between observations and the mechanistic model Our fitted model projects the trajectory of number of daily and cumulative confirmed cases in Shincheonji and in the rest of the population of Daegu City (Figure 3(a)-(d)). The model correctly projects the decreasing trends in both patches after reaching the peak on around March 3, 2020. However, for the non-Shincheonji, daily new cases are underestimated toward the end of the outbreak. *R* were estimated to be quite different across two patches (Figure 3(e)). The local reproduction number in the patch representing Shincheonji members, *R*0,1, was estimated to be 8.54 [95% credible interval (CrI): 6.30, 10.95] whereas the local reproduction number in the patch representing non-Shincheonji members, *R*0,2, was estimated to be 1.87 [95% CrI: 0.38, 4.40]. The time taken for the intervention program to have exerted highest effect, *d*, is around 9.02 days [95% CrI: 7.85, 10.45], which leads to both reduced transmission rate per unit time and reproduction numbers (*R*final = 0.34 [95% CrI: 0.18, 0.53]). The model also suggests that there were infectious people already when the first cases was symptomatic on February 7 (*I*= 4 [95% CrI: 2, 11]). The proportion of time that a person from Shincheonji members spends mixing with non-Shincheonji people, *c*12, was estimate be around 0.14 [95% CrI: 0.05, 0.22]. Posterior distribution of parameters based on 2,000 samples obtained from 10 different random seeds and two-way correlations appear in Supplementary Material (Figure S2). ## 4. Discussion Rapid transmission of COVID-19 within the Shincheonji community is likely to have been facilitated by high intensity contact between individuals gathering during services and in residential areas. Our mathematical modeling analyses quantify the rapid spread of COVID-19 in Daegu City driven by a community of Shincheonji members. The median *R* among Shincheonji members (*R*0,1) was 8.5, which is over 4-fold higher than what was estimated for the rest of the population in Daegu City (*R*0,2 =1.9). While the *R* in the Shincheonji community is higher than estimates from most transmission hotspots (e.g., in China (Alimohamadi et al., 2020; Imai et al., 2020; Riou and Althaus, 2020; Wu et al., 2020; Zhao et al., 2020b) and Korea (Bae et al., 2020; Choi and Ki, 2020; Ki, 2020; Shim et al., 2020a)), such high *R*is not unusual in particular considering that *R* can be different depending on the local settings with varying contact rates (Temime et al., 2020). Studies do report that *R* estimates of COVID-19 that are comparable or even higher than our estimates for the Shincheonji community. During periods of intensive social contacts near the Chinese New Year in China, *R* was estimated to be 6 (Sanche et al., 2020; Tang et al., 2020). Also, *R* estimates were around 5 among those traveled from Wuhan and were subsequently confirmed in other countries (Zhao et al., 2020a), and around 7 during the initial growth phase in the UK (Dropkin, 2020). In an extreme setting such as the Diamond Princess ship, much higher estimates (*R* = 14.8) were reported (Rocklöv et al., 2020). Although the previous studies that included data on the outbreak in Shincheonji community report smaller *R* estimates (Choi and Ki, 2020; Shim et al., 2020a) than our estimates, the difference might stem from that prior studies did not model the Shincheonji community separately from the rest of the population and therefore measured the *R* averaged across sub-population that are highly heterogeneous. Although estimated daily doubling times show some variability (e.g., 14 days on February 28 and 69.6 days on March 1 for the non-Shincheonji population), they are short overall, which indicates rapid growth of the outbreak, and are compatible with estimates from other settings. Daily doubling times were lower than one day in the beginning of the outbreak and this is similar to the estimates from several regions in China (Muniz-Rodriguez et al., 2020a). The study by Shim *et al*. (Shim et al., 2020b) used the dataset from Daegu City, Korea including Shincheonji population produced the doubling time of 2.8 days [95% CI: 2.5, 4.0]. Our daily doubling time estimates averaged over the period of February 18 – March 5 is 2.9 days and is consistent with the study. The period from February 18 to March 5 is likely to have been used in the study by Shim *et al*. because the authors calculated the doubling times on the days when the reported cases doubled and during the period of February 18 – March 5 the number of cases doubled 12 times and no further doubling occurred since then. The study by Lee *et al*. (Lee et al., 2020) used similar data, but reported seemingly inconsistent findings, doubling time of 2.9 days for the first week and 3.4 over the period around February 18 – March 4 considering that our estimate averaged over the first week is 0.9 day. One likely reason for this difference is that Lee *et al*. calculated the doubling time using the cumulative incidence estimated from a logistic model that used the initial value (i.e., number of infected people on February 18) as a free parameter. Figure 2D from their study indicates that the number of infected people on 18 February is much larger than 1 and this might have led to the higher doubling time than our estimates. This may also explain why Lee *et al*. estimates for a similar period (i.e., February 18 – March 5) is higher than our estimates and those by Shim *et al*. (Shim et al., 2020b). The relationship between the doubling time and the *R* provides two insights on our inferences on *R*. For an *SEIR* model, there exists an algebraic formula that describes the inverse relationship between initial epidemic growth rate and *R* (Ma, 2020; Ma et al., 2014). This inverse relationship suggests that short doubling times during the early phase of the epidemic we calculated using the growth rates are consistent with high *R* for Shincheonji we estimated using the stochastic dynamic transmission model (Table S2 in the Supplementary Material). On the other hand, while doubling times may be reduced and imply high *R* for non-Shincheonji people as well, such short doubling times can arise through mixing (i.e., positive *c*12) with Shincheonji of high *R* even if the *R* for the non-Shincheonji people are not as high. Parameters around asymptomatic infections of COVID-19 are largely unknown (Fox et al., 2020) and we tested the sensitivity of our inferences to our assumptions on two parameters related to asymptomatic infection, namely the proportion of asymptomatic infection, *f*, and relative rate of isolation of asymptomatic people, *ρ* (Figure S3 in Supplementary Material). *R* for Shincheonji people, *R*0,1, and the final reproduction number, *R*final, showed a slight increase with increasing *f* or decreasing *ρ* while other parameter estimates remain relatively constant. We also tested the sensitivity of our parameter estimates to *α*final , maximum rate of isolation near the end of the outbreak. *α*final showed an inverse relationship with other intervention-related parameters such as duration of intervention *d*, and reproduction number at the end of the outbreak, *R*final. Overall, while there are some quantitative differences in our parameter estimates in response to the change in our assumptions on fixed parameters, *R*0,1 was always over 4-fold higher than *R*0,2. While the first case was confirmed on February 18 for the Shincheonji outbreak, it was later revealed that the first case had symptoms on February 7 and even earlier transmission events were also suspected (Korea Disease Control and Prevention Agency, 2020b). This finding is consistent with our model analyses, which suggest there were 4 [95% CrI: 2, 11] infectious people on February 7. These undetected cases are likely to have contributed to the explosive outbreak in the Shincheonji community. Studies suggest that a substantial fraction of all SARS-CoV-2 infections were undetected. For Korea, it was suggested that the number of undetected cases may be larger than the number of detected cases (Lee et al., 2021). A study suggests that > 80% of all infections were undocumented during the initial spread in China (Li et al., 2020b). In France, over the period of 7 weeks since 28 June 2020 after the first lockdown, it was estimated that around 93% of all symptomatic cases were undetected initially and later around 69% of symptomatic cases were undetected by the time when case ascertainment improved (Pullano et al., 2021). We have shown that SARS-CoV 2 has disproportionately affected a religious community generating a large cluster of linked cases in Korea. Similar large clusters of cases in high-risk settings have been observed in Korea and elsewhere. In Korea, many similar outbreaks in high-risk settings have been reported in the news including the outbreak in a dance class (Jang et al., 2020) and a call center (Park et al., 2020). In Singapore, a total of 247 cases were confirmed as of March 17, 2020 and six clusters including the spread in a hotel and in a church accounted for 45.3% of the total cases (Tariq et al., 2020). In Hong Kong. 1,038 cases were confirmed from January 23 to April 28, 2020 and among them, 51.3% of cases were associated with large clusters. Such social settings as bars, restaurants, weddings, and religious sites appeared at increased risk of large outbreaks (Adam et al., 2020). One limitation of our analyses is that the model was fit to date of case confirmation because the date of symptom onset, which is more closely related with the date of infection, was not available. The daily number of confirmed cases can abruptly change depending on the intensity of intervention measures, of which the dynamics may not be consistent with disease transmission process. This means using the data on case confirmation under dynamics intervention measures is challenging. We tried to mitigate this difficulty by incorporating the dynamics of intervention programs by assuming that the start date and duration of enhanced case detection vary while the case detection rate increases over time and let the data suggest the values for those parameters. ## 5. Conclusions The potential for large variations in *R* for COVID-19 has important implications for the design and effectiveness of control strategies. The efficacy of components of intervention programs, such as contact tracing and physical distancing, is dependent on various environmental and societal factors (e.g., large gatherings, physical proximity, high risk behaviors such as singing, etc.) that influence the transmissibility of disease. Our analyses provide important insights that in order to minimize the risk of sudden outbreaks, efforts to identify and preempt high transmission scenarios will be key to controlling the spread of the COVID-19. Understanding and subsequently limiting the risk of transmission in high-risk places such as the Shincheonji Church cluster in Korea is key to effective control of COVID-19 transmission. ## Data Availability All of the data referred to in the manuscript is publicly available. ## Funding This research was partly supported by Government-wide R&D Fund project for infectious disease research (GFID), Republic of Korea (grant number: HG18C0088) and National Institute for Mathematical Sciences (NIMS) grant funded by the Korean Government (NIMS-B21910000). ## Supplementary Material ### Data source Figure S1 shows the cumulative cases before and after imputation. We used the cubic spline method provided in the imputeTS package of R. ![Figure S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/01/2021.08.05.21261683/F1.medium.gif) [Figure S1.](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/F1) Figure S1. Cumulative number of cases in Shincheonji community and the overall Daegu before (left panel) and after imputation (right panel). ### Model equations The model is implemented using a tau-leap algorithm. The number of people who transit from state *x* to state *y* over the time interval from *t* to *t* +Δ*t* in patch ![Graphic][20] , is defined as follows: ![Formula][21] where Bin(*n, p*) represents binomial and with parameters *n* and *p*, respectively. ![Formula][22] where Multi(*n*, π ={π1, π2}) indicates multinomial distribution with parameters *n* and π. π is given as follows: ![Formula][23] The first element of *π* indicates a probability for *x* = *I* (i.e., transition from *P* to *I*) and the second element indicates the probability of transition from *P* to *A* (*x* = *A*). ![Formula][24] The first element of *π* indicates a probability for *x* = C (i.e., transition from *I* to *C*) and the second element indicates the probability of transition from *I* to *R* (*x* = *R*). ![Formula][25] The first element of *π* indicates a probability for *x* = C (i.e., transition from *A* to *C*) and the second element indicates the probability of transition from *A* to *R* (*x* = *R*). The number of people in each state at time *t* + Δ*t* can be described using the terms defined above: ![Formula][26] The model comprises two sets of above equations that describe two patches (i.e., a community of Shincheonji members and the non-Shincheonji people in Daegu City) and these equations are linked through the force of infection function, *λ*(*t*), which is defined in the main text. View this table: [Table S1.](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/T1) Table S1. Delay from onset of symptoms to isolation during the COVID-19 outbreak in Busan City, Korea. ### Model fitting A pseudocode for Approximate Bayesian Computation Sequential Monte Carlo (ABC-SMC) is presented below adopting what was presented in the previous study (Minter and Retkute, 2019): 1. Set the number of generations *G* and the number of particles *N* 2. Set the tolerance schedule *ϵ*1 < *ϵ*2 < *ϵ*3 < ⋯ < *ϵ**G* and set the generation indicator *g* = 1 3. Set the particle indicator *i* = 1 4. If *g* = 1, sample *θ*** from the prior distribution *p*(*θ*). If *g* > 1, sample *θ** from the previous generation {*θ**g*−1} with weights {*w**g*−1} and perturb the particle to obtain *θ***∽*K*(*θ*|*θ**) 5. If *p*(*θ***) = 0, return to Step 4. 6. Generate *n* data sets ![Graphic][27] from the model using *θ*** and calculate ![Formula][28] 7. If![Graphic][29] , return to Step 4 8. Set ![Graphic][30] and calculate the corresponding weight of the accepted particle *i* ![Formula][31] 9. If *i* < *N*, increment *i* = *i* + 1 and go to step 4. 10. Normalize the weights so that ![Graphic][32] 11. If *g* < *G*, set *g* = *g* + 1, go to step 3 *K*(*θ*|*θ**) was assumed to follow a multivariate normal distribution that was truncated to give only positive values. Twenty generations (i.e., *G* = 20) were used with the following tolerance for Shincheonji ![Graphic][33] and non-Shincheonji![Graphic][34] : Initial tolerance values ![Graphic][35] were set as ![Formula][36] where *y**j* represent incidence of confirmed case on day *j* for Shincheonji (*S* = 1) and non-Shincheonji (*S* = 2), respectively. ![Graphic][37] values were determined by setting the minimum values, ![Graphic][38] , as ![Graphic][39] and dividing 20 equidistance pieces. Finally, final two value were manually adjusted to provide good fit between data and the model predictions through trial and error. Below are actual values used but were rounded for presentation. ![Formula][40] Prior distributions for the parameters *θ* = (*R*0,1, *R*0,2*I*, *c*12, *d, R*final) were defined as uniform distribution as follows: *R*0,1 ∽ *U*(1, 20), *R*0,2 ∽ *U*(1, 20), *I* ∽ *U*(1, 20), *c*12 ∽ *U*(0.000001, 1), *d* ∽ *U*(1, 30), *R*final ∽ *U*(1, 20). ![Figure S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/01/2021.08.05.21261683/F2.medium.gif) [Figure S2.](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/F2) Figure S2. Posterior distribution of model parameters (*n*=2000). For each of 10 random seeds, 200 samples were generated. ### Growth rate *r* and basic reproduction number *R* For the differential equation-based *SIR* model, the initial (*i.e*., the entire population is susceptible) epidemic growth rate *r** can be given as *β* − *γ* (Ma, 2020), where *β* and *γ* represent transmission rate and recovery rate, respectively, as we defined in our model. Similarly, for a differential equation-based *SEIR* model *r** is given as below (Ma, 2020; Ma et al., 2014): ![Formula][41] where *ϵ* represent the rate at which the exposed individuals become infectious (*i.e*., ![Graphic][42]) as we defined in the main text. The above equation gives *β* and therefore *R* for given *r**, *γ, ϵ*. Assuming *r** = *r*, which is the growth rate we calculated in the main text, we can see to what value of *R* the doubling times we calculated in the main text are translated and qualitatively see if *R* estimates from the current study are reasonable. View this table: [Table S2.](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/T2) Table S2. Basic reproduction number, *R*, calculated by assuming the empirical daily or weekly growth rate *r* is the same as the *r** calculated for the SEIR model. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/01/2021.08.05.21261683/F3/graphic-32.medium.gif) [](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/F3/graphic-32) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/01/2021.08.05.21261683/F3/graphic-33.medium.gif) [](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/F3/graphic-33) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/01/2021.08.05.21261683/F3/graphic-34.medium.gif) [](http://medrxiv.org/content/early/2021/11/01/2021.08.05.21261683/F3/graphic-34) Figure S3. Sensitivity of our parameter estimates to simplifying assumption of three selected parameters ## Acknowledgments All authors acknowledge discussions with the members of the Research and Development on Integrated Surveillance System Development for Early Warning of Infectious Diseases of Korea. ## Footnotes * Author affiliations updated; Supplemental files updated. * Received August 5, 2021. * Revision received November 1, 2021. * Accepted November 1, 2021. * © 2021, 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/) ## 6. References 1. Adam, D.C., Wu, P., Wong, J.Y., Lau, E.H.Y., Tsang, T.K., Cauchemez, S., Leung, G.M., Cowling, B.J., 2020. Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong. Nat Med 26, 1714–1719. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 2. Alimohamadi, Y., Taghdir, M., Sepandi, M., 2020. The Estimate of the Basic Reproduction Number for Novel Coronavirus disease (COVID-19): A Systematic Review and Meta-Analysis. Korean J Prev Med . 3. Anderson, R.M., Donnelly, C., Hollingsworth, D., Keeling, M., Vegvari, C., Baggaley, R., R, M., 2020. Reproduction number (R) and growth rate (r) of the COVID-19 epidemic in the UK: methods of estimation, data sources, causes of heterogeneity, and use as a guide in policy formulation. The Royal Society. 4. Azman, A.S., Lessler, J., 2015. Reactive vaccination in the presence of disease hotspots. Proc Biol Sci 282, 20141341. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2014.1341&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25392464&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 5. Bae, T.W., Kwon, K.K., Kim, K.H., 2020. Mass Infection Analysis of COVID-19 Using the SEIRD Model in Daegu-Gyeongbuk of Korea from April to May, 2020. J Korean Med Sci 35, e317. 6. Chan, J.F., Yuan, S., Kok, K.H., To, K.K., Chu, H., Yang, J., Xing, F., Liu, J., Yip, C.C., Poon, R.W., Tsoi, H.W., Lo, S.K., Chan, K.H., Poon, V.K., Chan, W.M., Ip, J.D., Cai, J.P., Cheng, V.C., Chen, H., Hui, C.K., Yuen, K.Y., 2020. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet 395, 514–523. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30154-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 7. Chen, M.K., Chevalier, J.A., Long, E.F., 2021. Nursing home staff networks and COVID-19. Proc Natl Acad Sci U S A 118. 8. Choe, S.-H., 2020. Shadowy Church Is at Center of Coronavirus Outbreak in South Korea, The New York Times. 9. Choi, S., Ki, M., 2020. Estimating the reproductive number and the outbreak size of COVID-19 in Korea. Epidemiol Health 42, e2020011. 10. Chung, E., Hill, A., 2020. [DEBRIEFING] What is the Shincheonji Church of Jesus and who are its members? And more importantly, what are its links to the coronavirus?, Korea Joongang Daily. 11. Dropkin, G., 2020. COVID-19 UK Lockdown Forecasts and R0. Front Public Health 8. 12. Ebell, M.H., Bagwell-Adams, G., 2020. Mandatory Social Distancing Associated With Increased Doubling Time: An Example Using Hyperlocal Data. Am J Prev Med 59, 140–142. 13. Fox, S.J., Pasco, R., Tec, M., Du, Z., Lachmann, M., Scott, J., Meyers, L.A., 2020. The impact of asymptomatic COVID-19 infections on future pandemic waves. medRxiv, 2020.2006.2022.20137489. 14. Fu, L., Wang, B., Yuan, T., Chen, X., Ao, Y., Fitzpatrick, T., Li, P., Zhou, Y., Lin, Y.F., Duan, Q., Luo, G., Fan, S., Lu, Y., Feng, A., Zhan, Y., Liang, B., Cai, W., Zhang, L., Du, X., Li, L., Shu, Y., Zou, H., 2020. Clinical characteristics of coronavirus disease 2019 (COVID-19) in China: A systematic review and meta-analysis. J Infect 80, 656–665. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2020.03.041&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 15. Gillespie, D.T., 2001. Approximate accelerated stochastic simulation of chemically reacting systems. J Chem Phys 115, 1716–1733. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1063/1.1378322&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169800000010&link_type=ISI) 16. Guan, W.J., Ni, Z.Y., Hu, Y., Liang, W.H., Ou, C.Q., He, J.X., Liu, L., Shan, H., Lei, C.L., Hui, D.S.C., Du, B., Li, L.J., Zeng, G., Yuen, K.Y., Chen, R.C., Tang, C.L., Wang, T., Chen, P.Y., Xiang, J., Li, S.Y., Wang, J.L., Liang, Z.J., Peng, Y.X., Wei, L., Liu, Y., Hu, Y.H., Peng, P., Wang, J.M., Liu, J.Y., Chen, Z., Li, G., Zheng, Z.J., Qiu, S.Q., Luo, J., Ye, C.J., Zhu, S.Y., Zhong, N.S., 2020. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med 382, 1708–1720. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2002032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 17. Guirao, A., 2020. The Covid-19 outbreak in Spain. A simple dynamics model, some lessons, and a theoretical framework for control response. Infectious Disease Modelling 5, 652–669. 18. Imai, N., Cori, A., Dorigatti, I., Baguelin, M., Donnelly, C., Riley, S., Ferguson, N., 2020. Report 3: Transmissibility of 2019-nCoV. 19. Jang, S., Han, S.H., Rhee, J.Y., 2020. Cluster of Coronavirus Disease Associated with Fitness Dance Classes, South Korea. Emerg Infect Dis 26, 1917–1920. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2608.200633&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32412896&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 20. Jones, N.R., Qureshi, Z.U., Temple, R.J., Larwood, J.P.J., Greenhalgh, T., Bourouiba, L., 2020. Two metres or one: what is the evidence for physical distancing in covid-19? BMJ 370, m3223. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNzAvYXVnMjVfNi9tMzIyMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzExLzAxLzIwMjEuMDguMDUuMjEyNjE2ODMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 21. Ki, M., 2020. Epidemiologic characteristics of early cases with 2019 novel coronavirus (2019-nCoV) disease in Korea. 22. Kim, J.-H., 2021. COVID_Shincheonji, GitHub repository [https://github.com/kimfinale/COVID\_Shincheonji](https://github.com/kimfinale/COVID_Shincheonji). GitHub. 23. Koopman, J., Simon, C., Jacquez, J., Joseph, J., Sattenspiel, L., Park, T., 1988. Sexual partner selectiveness effects on homosexual HIV transmission dynamics. J Acquir Immune Defic Syndr (1988) 1, 486–504. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3221321&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 24. Korea Disease Control and Prevention Agency, 2020a. Press release. 25. Korea Disease Control and Prevention Agency, 2020b. Press release (February 22). 26. Lee, C., Apio, C., Park, T., 2021. Estimation of Undetected Asymptomatic COVID-19 Cases in South Korea Using a Probabilistic Model. Int J Environ Res Public Health 18. 27. Lee, S., Castillo-Chavez, C., 2015. The role of residence times in two-patch dengue transmission dynamics and optimal strategies. J Theor Biol 374, 152–164. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2015.03.005&link_type=DOI) 28. Lee, W., Hwang, S.-S., Song, I., Park, C., Kim, H., Song, I.-K., Choi, H.M., Prifti, K., Kwon, Y., Kim, J., Oh, S., Yang, J., Cha, M., Kim, Y., Bell, M.L., Kim, H., 2020. COVID-19 in South Korea: epidemiological and spatiotemporal patterns of the spread and the role of aggressive diagnostic tests in the early phase. International Journal of Epidemiology 49, 1106–1116. 29. Lemieux, J.E., Siddle, K.J., Shaw, B.M., Loreth, C., Schaffner, S.F., Gladden-Young, A., Adams, G., Fink, T., Tomkins-Tinch, C.H., Krasilnikova, L.A., DeRuff, K.C., Rudy, M., Bauer, M.R., Lagerborg, K.A., Normandin, E., Chapman, S.B., Reilly, S.K., Anahtar, M.N., Lin, A.E., Carter, A., Myhrvold, C., Kemball, M.E., Chaluvadi, S., Cusick, C., Flowers, K., Neumann, A., Cerrato, F., Farhat, M., Slater, D., Harris, J.B., Branda, J., Hooper, D., Gaeta, J.M., Baggett, T.P., O’Connell, J., Gnirke, A., Lieberman, T.D., Philippakis, A., Burns, M., Brown, C.M., Luban, J., Ryan, E.T., Turbett, S.E., LaRocque, R.C., Hanage, W.P., Gallagher, G.R., Madoff, L.C., Smole, S., Pierce, V.M., Rosenberg, E., Sabeti, P.C., Park, D.J., Maclnnis, B.L., 2020. Phylogenetic analysis of SARS-CoV-2 in the Boston area highlights the role of recurrent importation and superspreading events. medRxiv, 2020.2008.2023.20178236. 30. Li, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., Ren, R., Leung, K.S.M., Lau, E.H.Y., Wong, J.Y., Xing, X., Xiang, N., Wu, Y., Li, C., Chen, Q., Li, D., Liu, T., Zhao, J., Liu, M., Tu, W., Chen, C., Jin, L., Yang, R., Wang, Q., Zhou, S., Wang, R., Liu, H., Luo, Y., Liu, Y., Shao, G., Li, H., Tao, Z., Yang, Y., Deng, Z., Liu, B., Ma, Z., Zhang, Y., Shi, G., Lam, T.T.Y., Wu, J.T., Gao, G.F., Cowling, B.J., Yang, B., Leung, G.M., Feng, Z., 2020a. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med 382, 1199–1207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 31. Li, R., Pei, S., Chen, B., Song, Y., Zhang, T., Yang, W., Shaman, J., 2020b. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 368, 489–493. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMS8wMS8yMDIxLjA4LjA1LjIxMjYxNjgzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 32. Lurie, M.N., Silva, J., Yorlets, R.R., Tao, J., Chan, P.A., 2020. Coronavirus Disease 2019 Epidemic Doubling Time in the United States Before and During Stay-at-Home Restrictions. The Journal of Infectious Diseases 222, 1601–1606. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiaa491&link_type=DOI) 33. Ma, J., 2020. Estimating epidemic exponential growth rate and basic reproduction number. Infect Dis Model 5, 129–141. 34. Ma, J., Dushoff, J., Bolker, B.M., Earn, D.J., 2014. Estimating initial epidemic growth rates. Bull Math Biol 76, 245–260. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11538-013-9918-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24272389&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 35. Minter, A., Retkute, R., 2019. Approximate Bayesian Computation for infectious disease modelling. Epidemics 29, 100368. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100368&link_type=DOI) 36. Muniz-Rodriguez, K., Chowell, G., Cheung, C.-H., Jia, D., Lai, P.-Y., Lee, Y., Liu, M., Ofori, S., Roosa, K., Simonsen, L., Viboud, C., Fung, I.C.-H., 2020a. Doubling Time of the COVID-19 Epidemic by Province, China. Emerg Infect Dis 26, 1912. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 37. Muniz-Rodriguez, K., Chowell, G., Cheung, C.-H., Jia, D., Lai, P.-Y., Lee, Y., Liu, M., Ofori, S.K., Roosa, K.M., Simonsen, L., Viboud, C., Fung, I.C.-H., 2020b. Doubling Time of the COVID-19 Epidemic by Province, China. Emerging infectious diseases 26, 1912–1914. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 38. Myung, M.-j., 2020. Apartment block enters cohort isolation for the first time, THE DONG-A ILBO. 39. O’Driscoll, M., Ribeiro Dos Santos, G., Wang, L., Cummings, D.A.T., Azman, A.S., Paireau, J., Fontanet, A., Cauchemez, S., Salje, H., 2021. Age-specific mortality and immunity patterns of SARS-CoV-2. Nature 590, 140–145. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 40. Park, S.Y., Kim, Y.M., Yi, S., Lee, S., Na, B.J., Kim, C.B., Kim, J.I., Kim, H.S., Kim, Y.B., Park, Y., Huh, I.S., Kim, H.K., Yoon, H.J., Jang, H., Kim, K., Chang, Y., Kim, I., Lee, H., Gwack, J., Kim, S.S., Kim, M., Kweon, S., Choe, Y.J., Park, O., Park, Y.J., Jeong, E.K., 2020. Coronavirus Disease Outbreak in Call Center, South Korea. Emerg Infect Dis 26, 1666–1670. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 41. Patel, S.B., Patel, P., 2020. Doubling Time and its Interpretation for COVID 19 Cases. Natl J Community Med 11, 141–143. 42. Pullano, G., Di Domenico, L., Sabbatini, C.E., Valdano, E., Turbelin, C., Debin, M., Guerrisi, C., Kengne-Kuetche, C., Souty, C., Hanslik, T., Blanchon, T., Boëlle, P.Y., Figoni, J., Vaux, S., Campèse, C., Bernard-Stoecklin, S., Colizza, V., 2021. Underdetection of cases of COVID-19 in France threatens epidemic control. Nature 590, 134–139. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 43. Riou, J., Althaus, C.L., 2020. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. Euro Surveill 25. 44. Rocklöv, J., Sjödin, H., Wilder-Smith, A., 2020. COVID-19 outbreak on the Diamond Princess cruise ship: estimating the epidemic potential and effectiveness of public health countermeasures. J Travel Med 27. 45. Rothe, C., Schunk, M., Sothmann, P., Bretzel, G., Froeschl, G., Wallrauch, C., Zimmer, T., Thiel, V., Janke, C., Guggemos, W., Seilmaier, M., Drosten, C., Vollmar, P., Zwirglmaier, K., Zange, S., Wölfel, R., Hoelscher, M., 2020. Transmission of 2019-nCoV Infection from an Asymptomatic Contact in Germany. New England Journal of Medicine 382, 970–971. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2001468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 46. Russell, T.W., Hellewell, J., Jarvis, C.I., van Zandvoort, K., Abbott, S., Ratnayake, R., Cmmid Covid-Working, G., Flasche, S., Eggo, R.M., Edmunds, W.J., Kucharski, A.J., 2020. Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship, February 2020. Euro Surveill 25. 47. Sanche, S., Lin, Y.T., Xu, C., Romero-Severson, E., Hengartner, N., Ke, R., 2020. The Novel Coronavirus, 2019-nCoV, is Highly Contagious and More Infectious Than Initially Estimated. medRxiv, 2020.2002.2007.20021154. 48. Shim, E., Tariq, A., Choi, W., Lee, Y., Chowell, G., 2020a. Transmission potential and severity of COVID-19 in South Korea. Int J Infect Dis 93, 339–344. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.03.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 49. Shim, E., Tariq, A., Chowell, G., 2020b. Spatial variability in reproduction number and doubling time across two waves of the COVID-19 pandemic in South Korea, February to July, 2020. Int J Infect Dis 102, 1–9. 50. Tang, B., Wang, X., Li, Q., Bragazzi, N.L., Tang, S., Xiao, Y., Wu, J., 2020. Estimation of the Transmission Risk of the 2019-nCoV and Its Implication for Public Health Interventions. J Clin Med 9. 51. Tariq, A., Lee, Y., Roosa, K., Blumberg, S., Yan, P., Ma, S., Chowell, G., 2020. Real-time monitoring the transmission potential of COVID-19 in Singapore, March 2020. BMC Med 18, 166. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12916-020-01615-9&link_type=DOI) 52. Temime, L., Gustin, M.-P., Duval, A., Buetti, N., Crépey, P., Guillemot, D., Thiébaut, R., Vanhems, P., Zahar, J.-R., Smith, D.R.M., Opatowski, L., 2020. A Conceptual Discussion About the Basic Reproduction Number of Severe Acute Respiratory Syndrome Coronavirus 2 in Healthcare Settings. Clinical Infectious Diseases. 53. Wu, J.T., Leung, K., Leung, G.M., 2020. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet. 54. Xu, X.K., Liu, X.F., Wu, Y., Ali, S.T., Du, Z., Bosetti, P., Lau, E.H.Y., Cowling, B.J., Wang, L., 2020. Reconstruction of Transmission Pairs for Novel Coronavirus Disease 2019 (COVID-19) in Mainland China: Estimation of Superspreading Events, Serial Interval, and Hazard of Infection. 55. Yonhap, 2020. ‘Shincheonji’ suspected as coronavirus hotbed, The Korea Times. 56. Yu, P., Zhu, J., Zhang, Z., Han, Y., 2020. A Familial Cluster of Infection Associated With the 2019 Novel Coronavirus Indicating Possible Person-to-Person Transmission During the Incubation Period. J Infect Dis 221, 1757–1761. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 57. Zhao, Q., Chen, Y., Small, D.S., 2020a. Analysis of the epidemic growth of the early 2019-nCoV outbreak using internationally confirmed cases. medRxiv, 2020.2002.2006.20020941. 58. Zhao, S., Lin, Q., Ran, J., Musa, S.S., Yang, G., Wang, W., Lou, Y., Gao, D., Yang, L., He, D., Wang, M.H., 2020b. 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. Int J Infect Dis 92, 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%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) ## References 1. Ma, J., 2020. Estimating epidemic exponential growth rate and basic reproduction number. Infect Dis Model 5, 129–141. 2. Ma, J., Dushoff, J., Bolker, B.M., Earn, D.J., 2014. Estimating initial epidemic growth rates. Bull Math Biol 76, 245–260. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11538-013-9918-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24272389&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F01%2F2021.08.05.21261683.atom) 3. Minter, A., Retkute, R., 2019. Approximate Bayesian Computation for infectious disease modelling. Epidemics 29, 100368. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100368&link_type=DOI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-6.gif [7]: /embed/inline-graphic-1.gif [8]: /embed/graphic-7.gif [9]: /embed/graphic-8.gif [10]: /embed/graphic-9.gif [11]: /embed/graphic-10.gif [12]: /embed/graphic-11.gif [13]: /embed/inline-graphic-2.gif [14]: /embed/inline-graphic-3.gif [15]: /embed/graphic-12.gif [16]: /embed/graphic-13.gif [17]: /embed/graphic-14.gif [18]: /embed/graphic-15.gif [19]: /embed/inline-graphic-4.gif [20]: /embed/inline-graphic-5.gif [21]: /embed/graphic-17.gif [22]: /embed/graphic-18.gif [23]: /embed/graphic-19.gif [24]: /embed/graphic-20.gif [25]: /embed/graphic-21.gif [26]: /embed/graphic-22.gif [27]: /embed/inline-graphic-6.gif [28]: /embed/graphic-24.gif [29]: /embed/inline-graphic-7.gif [30]: /embed/inline-graphic-8.gif [31]: /embed/graphic-25.gif [32]: /embed/inline-graphic-9.gif [33]: /embed/inline-graphic-10.gif [34]: /embed/inline-graphic-11.gif [35]: /embed/inline-graphic-12.gif [36]: /embed/graphic-26.gif [37]: /embed/inline-graphic-13.gif [38]: /embed/inline-graphic-14.gif [39]: /embed/inline-graphic-15.gif [40]: /embed/graphic-27.gif [41]: /embed/graphic-29.gif [42]: /embed/inline-graphic-16.gif