{'fraction_non_alphanumeric': 0.02681274326077555, 'fraction_numerical': 0.013406371630387775, 'mean_word_length': 4.9604810996563575, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Authors’ response: Estimating the generation interval for COVID-19 based on symptom onset data

Estimating key infectious disease parameters from the COVID-19 outbreak is quintessential for modelling studies and guiding intervention strategies. Whereas different estimates for the incubation period distribution and the serial interval distribution have been reported, estimates of the generation interval for COVID-19 have not been provided.We used outbreak data from clusters in Singapore and Tianjin, China to estimate the generation interval from symptom onset data while acknowledging uncertainty about the incubation period distribution and the underlying transmission network. From those estimates we obtained the proportions pre-symptomatic transmission and reproduction numbers.The mean generation interval was 5.20 (95%CI 3.78-6.78) days for Singapore and 3.95 (95%CI 3.01-4.91) days for Tianjin, China when relying on a previously reported incubation period with mean 5.2 and SD 2.8 days. The proportion of pre-symptomatic transmission was 48% (95%CI 32-67%) for Singapore and 62% (95%CI 50-76%) for Tianjin, China. Estimates of the reproduction number based on the generation interval distribution were slightly higher than those based on the serial interval distribution.Estimating generation and serial interval distributions from outbreak data requires careful investigation of the underlying transmission network. Detailed contact tracing information is essential for correctly estimating these quantities.

By expressing the density of the serial interval Z i as a convolution of X i and Y i , we indeed make the simplifying assumption that the incubation period of the infector, v (i) , is independent of the generation time X i = t i − t v(i) . The possible correlation between the incubation period of the infector and the generation time should be taken into account. Ideally, both should be estimated from the same data. Unfortunately, we did not have these data directly available. Instead we assumed these quantities to be independent and we acknowledge that this assumption may not be realistic. However, to the best of our knowledge, the literature does not yet report clear indications of a strong relation between infectiousness and incubation period for coronavirus disease (COVID-19), with highly varying findings between studies [1]. It should be kept in mind that if our assumption of independence is not valid, our model is mis-specified as the convolution Z = X + Y is defined for independent random variables.
Liu et al. [3] have investigated the impact of correlation between incubation period and serial interval on estimates of presymptomatic transmission. They found that, in the presence of active case finding and assuming a mean serial interval of 4.8 days and mean incubation period of 5.2 days, the percentage of presymptomatic transmission was 48% when assuming no correlation, ranging from 38% under positive correlation to 47% under negative correlation. Our estimate of the proportion of presymptomatic transmission in Singapore (i.e. 48%; where the mean serial interval was 5.21 days when allowing only positive serial intervals and assuming a mean incubation period of 5.2 days) is in line with these estimates, with the credible interval (32-67%) also covering the lower estimate.

Assumption of independent serial intervals
Our likelihood function is indeed an approximation assuming that the serial intervals Z i are independent and identically distributed. Although our estimates are in line with those from other studies [4], we cannot be completely sure how much the assumption of independent serial intervals affects our estimates since those studies might involve the same or other statistical issues as well. This simplifying assumption of independent serial intervals is commonly used (e.g. [5]). An advantage of our method of inferring serial and generation intervals is that we take into account the infectious history of individuals by inferring these quantities based on an epidemic tree. This is in contrast to other studies that have estimated serial intervals using only probable infector-infectee pairs from different settings (e.g. [6]), which could lead to bias by not accounting for infectious histories.
It would be interesting to investigate the effect of this dependency between serial intervals in future work. These may be especially important when considering superspreading events, where there can be long transmission chains of the type i → j, i → k, i → l, i → m, …, when the same individual i generates a large number of secondary cases.

Possibility of cycles in the Singapore data
We have indeed overlooked the fact that the way contacts are defined in the Singapore data may lead to cycles in the network when allowing serial intervals to be negative, when not accounting for directionality (i.e. a case cannot infect any of its ancestors). This could potentially lead to biased estimates of serial and generation intervals. When only allowing missing serial intervals to be positive, the problem of cycles in the network does not occur since infector-infectee pairs are then based on dates of symptom onset.
We have re-analysed the Singapore data to investigate whether this would substantially change our results. In the original Markov chain Monte Carlo algorithm, we implemented a condition that the sampled network will only be accepted if it does not contain cycles. For computational reasons we make the comparison for a scenario where we allow serial intervals to be negative up to −3 days and sample a network every 100th iteration. Table 1 shows our estimates of the generation and serial interval as well as the proportion of presymptomatic transmission p. We can see that the estimates of the serial/generation interval when removing cycles are a bit higher, but the credible intervals largely overlap with those of the original analysis. As expected, the proportion of pre-symptomatic transmission is a bit lower now, but again with overlapping credible intervals. Hence, we do believe our overall conclusions in the original article would not change but for future work, one should make sure that no cycles in the network occur. We have also added an update of the R code on GitHub [7].

License, supplementary material and copyright
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY 4.0) Licence. You may share and adapt the material, but must give appropriate credit to the source, provide a link to the licence and indicate if changes were made.
Any supplementary material referenced in the article can be found in the online version.
This article is copyright of the authors or their affiliated institutions, 2020. The estimates are based on an analysis allowing for negative serial intervals up to −3 days and used an incubation period with mean 5.2 and standard deviation 2.8 days.