PMC:7107546 / 0-7 JSON TXT

Quantifying the association between domestic travel and the exportation of novel coronavirus (2019-nCoV) cases from Wuhan, China in 2020: a correlational analysis In the end of 2019, a novel coronavirus (2019-nCoV) emerged in Wuhan, China and was causing a serious outbreak of acute respiratory illness.1 Wuhan locates in the centre of mainland China with a population of 14 million and is very conveniently connected to other parts of China through airlines and high-speed rails.2 As of 31 January 2020 (5:00 p.m., GMT + 8), there were 9809 confirmed 2019-nCoV cases in mainland China, including 213 deaths and 180 discharges.3,4 Cases infected in Wuhan were also detected in many foreign countries or regions including Thailand, Japan, Republic of Korea, the United States, Canada and some European countries.4 The World Health Organization (WHO) has declared that the novel coronavirus outbreak is a public health emergency of international concern. Official reports on the newly confirmed cases are released very rapidly (several times a day) after January 16,3–5 as the official diagnosis protocol was released by the WHO on January 17.6 Recent studies indicated the likelihood of travel-related risks of 2019-nCoV spreading both domestically and internationally.7, 8 Many major cities in mainland China reported the finding of imported cases, including Beijing, Shanghai, Guangzhou and Shenzhen. The outbreak is still on-going with an increasing trend in daily new cases.3,4 Before the Wuhan lockdown (official travel restriction) on January 23, virtually all cases found in other major cities were exported cases from Wuhan. Population flow data between major cities in mainland China are available online due to the rapid development of internet in recent decades, see https://qianxi.baidu.com/ (in Chinese). In this work, we quantified the association between the domestic travel load and the number of cases exported from Wuhan to other city-clusters in mainland China. Our city-clusters are the pool of top five cities in the top 10 provinces (in number of cumulative cases). Thus, we included 10 city-clusters in the analysis, and the details of the selected city-clusters can be found in Supplementary Data S1. We examined the association between the load of domestic passengers departed from Wuhan and the number of confirmed cases to the 10 city-clusters (including the three municipalities, Beijing, Shanghai and Chongqing). Data sets of the daily numbers of domestic passengers were obtained from the location-based services database of Baidu company from January 1 to 20. We selected the top 10 provincial regions (except Hubei) with the largest cumulative number of cases, which accounts for 68% of all cases reported outside Hubei, before the implementation of the city lockdown, on 23 January 2020. Cases from other provinces are scattered and showed no clear pattern. The daily numbers of passengers from Wuhan to the city-clusters of each province were adopted to measure the load of domestic passengers departed from Wuhan to that cluster. From now on, province means the city-clusters in that province. The daily number of cases time series are obtained via the online outbreak situation reports.3,4 Daily cases for each cluster are scaled from daily province total. The detail is given in Supplementary Data. The association was formulated as follows:(1) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}\begin{equation*} \log \left[\mathbf{E}\left({c}_{i,t}\right)\right]={\alpha}_i\cdot{\mathrm{province}}_i+\beta \cdot{\xi}_{i,t-\tau}\cdot{\varepsilon}_{t-\tau }. \end{equation*}\end{document} Here, ci,t represents the daily number of new cases in the ith provincial region on day t. The function E(\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}) is the expectation; the ‘provincei’ denotes the dummy variable for the ith provincial region in accounting for the heterogenicity among different provinces, and thus, αi is a locality-varying interception term. Term ξi is the daily number of passengers from Wuhan to the ith province. The time index was denoted by t, and the term τ modelled the delay from exposure to being detected. The ε represented the daily number of new cases in Wuhan. Hence, the product term (ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ε) was proportional to the rate of cases exported from Wuhan to other places, and this setting was consistent with the framework in Imai et al.,9 which estimated the outbreak size through exported 2019-nCoV cases overseas, as well as in more complex frameworks.10 The β is the regression coefficient to quantify the association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan. Hence, the change rate, denoted by Δ = [exp(β × 100) − 1] × 100%, may be interpreted as the expected percentage change in the daily number of cases offsite found associated with per 100 increase in the daily number of passengers departed from Wuhan where there was one new case daily. We estimated β for three means of transportation, and the P-value less than 0.05 was considered as statistical significance. Table 1 Summary of the estimated percentage change rate (%) modelled in Equation (1) between daily passengers and number of exported 2019-nCoV cases The highlighted estimates were treated as the main results. Note: the ‘LR test’ is the likelihood-ratio (LR) test of the model in Equation (1) against the model in Equation (2). The ‘pseudo-R2’ is the McFadden’s pseudo-R-squared. We also considered a baseline version of the model by replacing term ‘β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ξi,t−τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}εt−τ’ in Equation (1) with ‘β0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ξi,t−τ’, which ignored the variation in the force of infection in Wuhan, thus a constant, as in Equation (2)(2) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}\begin{equation*} \log \left[\mathbf{E}\left({c}_{i,t}\right)\right]={\alpha^{(0)}}_i\cdot{\mathrm{province}}_i+{\beta}_0\cdot{\xi}_{i,t-\tau }. \end{equation*}\end{document} The likelihood ratio (LR) test was adopted to justify the model structure in Equation (1) against the baseline form in Equation (2). For the delay term (τ), it was expected to be equivalent to the incubation period, which was reported to be 5.2 days (95% CI: 4.1–7.0) in.1 Finally, we examined the association with τ varying at 3, 4, 5, 6 and 7 days in sensitivity analysis. The LR test yielded statistically significant outcomes, which suggests that the model in Equation (1) is more reasonable than the baseline form in Equation (2). Through the goodness-of-fit in terms of the McFadden’s pseudo-R-squared, we found that τ = 5 days attains the best fitting performance to explain the patterns of the cases offsite detected. This well matched the estimate of the mean incubation period of the infection at 5.2 in.1 We found a statistically significant positive association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan, see Table 1. We estimated that per 100 cases increase in the daily number of newly reported cases in Wuhan together with per 100 persons increase in the daily number of passengers departed Wuhan were likely to cause a 16.25% (95% CI: 14.86–17.66%) increase in the daily number of cases offsite detected on average. The sensitivity analysis by varying τ suggested this fundamental relationship holds, and the details can be found in Table 1 and the Supplementary Data S2. This analysis has limitations. The data used in this study are in the early phase of the outbreak and may not represent subsequent waves. We aim at quantifying the association between domestic travel and 2019-nCoV exportation in China. Although the differences in case ascertainment in different cities clusters are considered in the model by the dummy term ‘provincei’, the association could also bear heterogeneity. Temporal and spatial correlations were not addressed in this simple modelling analysis due to lack of data. The 2019-nCoV surveillance data were too scattered and short at this early stage to consider temporal and spatial correlation. The correlation between population flow, number of cases offsite detected, and source infection prevalence was addressed in this work. Our modelling framework would be extended to a more complex and realistic form for exploring the potential spatial correlations, and benefit from more detailed disease surveillance and travel population flow data. Supplementary Material dataset_taaa022 Click here for additional data file. supp_taaa022 Click here for additional data file.

Document structure show

Title	Quantifying the association between domestic travel and the exportation of novel coronavirus (2019-nCoV) cases from Wuhan, China in 2020: a correlational analysis
Body	In the end of 2019, a novel coronavirus (2019-nCoV) emerged in Wuhan, China and was causing a serious outbreak of acute respiratory illness.1 Wuhan locates in the centre of mainland China with a population of 14 million and is very conveniently connected to other parts of China through airlines and high-speed rails.2 As of 31 January 2020 (5:00 p.m., GMT + 8), there were 9809 confirmed 2019-nCoV cases in mainland China, including 213 deaths and 180 discharges.3,4 Cases infected in Wuhan were also detected in many foreign countries or regions including Thailand, Japan, Republic of Korea, the United States, Canada and some European countries.4 The World Health Organization (WHO) has declared that the novel coronavirus outbreak is a public health emergency of international concern. Official reports on the newly confirmed cases are released very rapidly (several times a day) after January 16,3–5 as the official diagnosis protocol was released by the WHO on January 17.6 Recent studies indicated the likelihood of travel-related risks of 2019-nCoV spreading both domestically and internationally.7, 8 Many major cities in mainland China reported the finding of imported cases, including Beijing, Shanghai, Guangzhou and Shenzhen. The outbreak is still on-going with an increasing trend in daily new cases.3,4 Before the Wuhan lockdown (official travel restriction) on January 23, virtually all cases found in other major cities were exported cases from Wuhan. Population flow data between major cities in mainland China are available online due to the rapid development of internet in recent decades, see https://qianxi.baidu.com/ (in Chinese). In this work, we quantified the association between the domestic travel load and the number of cases exported from Wuhan to other city-clusters in mainland China. Our city-clusters are the pool of top five cities in the top 10 provinces (in number of cumulative cases). Thus, we included 10 city-clusters in the analysis, and the details of the selected city-clusters can be found in Supplementary Data S1. We examined the association between the load of domestic passengers departed from Wuhan and the number of confirmed cases to the 10 city-clusters (including the three municipalities, Beijing, Shanghai and Chongqing). Data sets of the daily numbers of domestic passengers were obtained from the location-based services database of Baidu company from January 1 to 20. We selected the top 10 provincial regions (except Hubei) with the largest cumulative number of cases, which accounts for 68% of all cases reported outside Hubei, before the implementation of the city lockdown, on 23 January 2020. Cases from other provinces are scattered and showed no clear pattern. The daily numbers of passengers from Wuhan to the city-clusters of each province were adopted to measure the load of domestic passengers departed from Wuhan to that cluster. From now on, province means the city-clusters in that province. The daily number of cases time series are obtained via the online outbreak situation reports.3,4 Daily cases for each cluster are scaled from daily province total. The detail is given in Supplementary Data. The association was formulated as follows:(1) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}\begin{equation} \log \left[\mathbf{E}\left({c}_{i,t}\right)\right]={\alpha}_i\cdot{\mathrm{province}}_i+\beta \cdot{\xi}_{i,t-\tau}\cdot{\varepsilon}_{t-\tau }. \end{equation}\end{document} Here, ci,t represents the daily number of new cases in the ith provincial region on day t. The function E(\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}) is the expectation; the ‘provincei’ denotes the dummy variable for the ith provincial region in accounting for the heterogenicity among different provinces, and thus, αi is a locality-varying interception term. Term ξi is the daily number of passengers from Wuhan to the ith province. The time index was denoted by t, and the term τ modelled the delay from exposure to being detected. The ε represented the daily number of new cases in Wuhan. Hence, the product term (ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ε) was proportional to the rate of cases exported from Wuhan to other places, and this setting was consistent with the framework in Imai et al.,9 which estimated the outbreak size through exported 2019-nCoV cases overseas, as well as in more complex frameworks.10 The β is the regression coefficient to quantify the association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan. Hence, the change rate, denoted by Δ = [exp(β × 100) − 1] × 100%, may be interpreted as the expected percentage change in the daily number of cases offsite found associated with per 100 increase in the daily number of passengers departed from Wuhan where there was one new case daily. We estimated β for three means of transportation, and the P-value less than 0.05 was considered as statistical significance. Table 1 Summary of the estimated percentage change rate (%) modelled in Equation (1) between daily passengers and number of exported 2019-nCoV cases The highlighted estimates were treated as the main results. Note: the ‘LR test’ is the likelihood-ratio (LR) test of the model in Equation (1) against the model in Equation (2). The ‘pseudo-R2’ is the McFadden’s pseudo-R-squared. We also considered a baseline version of the model by replacing term ‘β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ξi,t−τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}εt−τ’ in Equation (1) with ‘β0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\cdot}$\end{document}ξi,t−τ’, which ignored the variation in the force of infection in Wuhan, thus a constant, as in Equation (2)(2) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}\begin{equation} \log \left[\mathbf{E}\left({c}_{i,t}\right)\right]={\alpha^{(0)}}_i\cdot{\mathrm{province}}_i+{\beta}_0\cdot{\xi}_{i,t-\tau }. \end{equation}\end{document} The likelihood ratio (LR) test was adopted to justify the model structure in Equation (1) against the baseline form in Equation (2). For the delay term (τ), it was expected to be equivalent to the incubation period, which was reported to be 5.2 days (95% CI: 4.1–7.0) in.1 Finally, we examined the association with τ varying at 3, 4, 5, 6 and 7 days in sensitivity analysis. The LR test yielded statistically significant outcomes, which suggests that the model in Equation (1) is more reasonable than the baseline form in Equation (2). Through the goodness-of-fit in terms of the McFadden’s pseudo-R-squared, we found that τ = 5 days attains the best fitting performance to explain the patterns of the cases offsite detected. This well matched the estimate of the mean incubation period of the infection at 5.2 in.1 We found a statistically significant positive association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan, see Table 1. We estimated that per 100 cases increase in the daily number of newly reported cases in Wuhan together with per 100 persons increase in the daily number of passengers departed Wuhan were likely to cause a 16.25% (95% CI: 14.86–17.66%) increase in the daily number of cases offsite detected on average. The sensitivity analysis by varying τ suggested this fundamental relationship holds, and the details can be found in Table 1 and the Supplementary Data S2. This analysis has limitations. The data used in this study are in the early phase of the outbreak and may not represent subsequent waves. We aim at quantifying the association between domestic travel and 2019-nCoV exportation in China. Although the differences in case ascertainment in different cities clusters are considered in the model by the dummy term ‘provincei’, the association could also bear heterogeneity. Temporal and spatial correlations were not addressed in this simple modelling analysis due to lack of data. The 2019-nCoV surveillance data were too scattered and short at this early stage to consider temporal and spatial correlation. The correlation between population flow, number of cases offsite detected, and source infection prevalence was addressed in this work. Our modelling framework would be extended to a more complex and realistic form for exploring the potential spatial correlations, and benefit from more detailed disease surveillance and travel population flow data. Supplementary Material dataset_taaa022 Click here for additional data file. supp_taaa022 Click here for additional data file.
Table caption	Table 1 Summary of the estimated percentage change rate (%) modelled in Equation (1) between daily passengers and number of exported 2019-nCoV cases The highlighted estimates were treated as the main results. Note: the ‘LR test’ is the likelihood-ratio (LR) test of the model in Equation (1) against the model in Equation (2). The ‘pseudo-R2’ is the McFadden’s pseudo-R-squared.
Section	Supplementary Material dataset_taaa022 Click here for additional data file. supp_taaa022 Click here for additional data file.
Title	Supplementary Material

projects that have annotations to this span

There is no project

TAB JSON ListView MergeView

PMC:7107546 / 0-7 JSONTXT

Document structure show

projects that have annotations to this span

PMC:7107546 / 0-7 JSON TXT