 Research article
 Open Access
 Open Peer Review
 Published:
A statistical method utilizing information of imported cases to estimate the transmissibility for an influenza pandemic
BMC Medical Research Methodologyvolume 17, Article number: 31 (2017)
Abstract
Background
In a new influenza pandemic, travel data such as arrival times of cases seeded by the originating country can be regarded as a combination of the epidemic size and the mobility networks of infections connecting the originating country with other regions. It can be a complete and timely source for estimating the basic reproduction number (R _{ 0 }), a key indicator of disease transmissibility.
Method
In this study, we developed a likelihoodbased method using arrival times of infected cases in different countries to estimate R _{ 0 } for influenza pandemics. A simulation was conducted to assess the performance of the proposed method. We further applied the method to the outbreak of the influenza pandemic A/H1N1 in Mexico.
Results
In the numerical application, the estimated R _{ 0 } was equal to 1.69 with a 95% confidence interval (1.65, 1.73). For the simulation results, the estimations were robust to the decline of travel rate and other parameter assumptions. Nevertheless, the estimates were moderately sensitive to the assumption of infectious duration. Generally, the findings were in line with other relevant studies.
Conclusions
Our approach as well as the estimate is potential to assist officials in planning control and prevention measures. Improved coordination to streamline or even centralize surveillance of imported cases among countries will thus be beneficial to public health.
Background
In responding to infectious disease outbreaks in the 21st century, a reliable method to estimate the transmission intensity during the early phases of a new influenza pandemic is critical. The basic reproduction number (R _{ 0 }), defined as the average number of secondary infections produced by a typical infectious individual in a whole susceptible population, is a common measure of disease transmissibility. Valid and reliable estimate of the R _{ 0 } can assist officials in planning control and prevention measures for an influenza pandemic. On the contrary, underestimation of the disease transmissibility will induce an insufficient public awareness on the risk of infection and therefore lower the general public’s incentive to protect themselves against infectious diseases through vaccines, hand washing or other protective measures.
Many kinds of data such as sentinel surveillance of influenzalikeillness (ILI), serological survey, and syndromic data are able to infer the transmissibility of infectious diseases. Common surveillance data can be fitted into an exponential growth model or KermackMcKendricktype models to estimate R _{ 0 } [1, 2]. Likelihoodbased methods are an alternative [3, 4]. Chowell et al. [5] showed that these types of methods are not sensitive on R _{ 0 } estimation given an acceptable goodness of fit. Nevertheless, one of the caveats of using ILI surveillance data is the underreporting although several approaches were developed to adjust this problem [6].
Serological data is another source for inferring transmissibility of an influenza virus [7, 8]. This kind of data helps detecting subclinical infections and the sampling is not affected by the reporting practices. However, it is comparatively costly and requires laboratory resources. A longer time will be taken to ascertain a disease and thus it cannot provide initial estimate for a pandemic outbreak [9]. Seroprevalence surveys are usually suggested to monitor the level of disease spread after a local outbreak in community [8].
In a recent decade, syndromic data such as web search queries [10] are alternative sources to infer estimates of key epidemiological parameters. Nevertheless, the human behaviors of web searching can be highly affected by intensive media coverage and changes of public’s perceptions on disease severity.
Apart from the abovementioned data sources, travel data related to imported cases could be an alternative information to provide an estimate of transmissibility for a new influenza outbreak. During the initial outbreak of a new influenza pandemic, control measures will be initiated at the border points of entry by various concerned countries. On preventing infected persons from seeding local epidemics, local governments may identify the symptomatic passengers by implementing thermal screening, monitor the suspected cases having a travel history to the originating country and even take quarantine measures [11]. Some travel information such as arrival times of infected cases seeded by the originating country is a combination of the epidemic size and the transport network between the originating country and other atrisk countries. By adapting this information, the transmissibility of an originating country during the initial phase of a pandemic can be backcalculated in order to reduce errors from the underestimation of routine surveillance data [9, 12, 13]. The data about importation events is usually timelier than other available data sources such as seroprevalence data. In the H1N1 pandemic, empirical evidence showed that the flow of airline passenger from Mexico had a significant correlation to the detection of exported cases [14]. Because of intense screening for ILI travelers returning from Mexico, ascertainment of early cases in atrisk countries was timelier and complete compared with routine surveillance in Mexico. Hence, it is believed that the data type is able to provide a reliable estimate on disease transmissibility.
In this paper, we developed a likelihoodbased method to estimate the R _{ 0 } for a new influenza pandemic using the information of imported cases i.e. time of the first introduction of infected individuals in different countries seeded by the originating country. The method was demonstrated to the influenza pandemic A/H1N1 (pH1N1) in midMarch 2009.
Methods
Mathematical model
SusceptibleInfectiousRecovered (SIR) model is a common model to describe the dynamic system of the infectious disease [1]. For each time point t (t = 0, 1, 2, 3…), a closed population is divided into three groups (‘compartments’), namely the susceptible (S(t)), the infectious (I(t)) and the recovered (R(t)) populations. In this compartmental model, susceptible individuals (S(t)) are infected at a transmission rate β; infectious individuals (I(t)) recover at a rate γ. Given an exponential assumption, the length of the infectious period is equal to 1/γ. Using S, I, and R to represent each compartment, the system dynamics can be described by the differential equations:
S _{ L } is the susceptible size of local population and S _{ M } is the number of visitors in the originating country. Because S _{ M } is far smaller than S _{ L } i.e. S _{ L } ≫ S _{ M }, the number of susceptible individuals is approximately equal to the population N i.e. S _{ L } ≈ N in the initial phase of a new epidemic. We assumed that there is no prior immunity for all subjects. We approximated
and
In the disease transmission process, the basic reproduction number R _{ 0 } is a measure of disease transmission intensity. By using the linearization method [2], R _{ 0 } is always equal to β/γ by assuming that the whole population is susceptible at first. Given only one infection at time 0, the prevalence I(t) in Eq. (2) is
Distribution of arrival times of infected individuals in countries seeded by originating country
Supposed the passengers counts between the seed country and a particular country ith at day t is m _{ i } (t) and d is the average length of stay the travellers in the seed country i.e. S _{ M(t)} = d∑_{ i } m _{ i }(t), the product (βd/N)m _{ i } (t)I(t) is the average number of infected cases arisen in country ith at time t. Same daily exposure risk to the domestic cases I(t) and uniformity of duration of stay are assumed. By assuming Poisson generations, the probability of importing at least an infected individual (Y _{ i } ) from the originating country at time t is
where θ is the proportion of asymptomatic individuals that cannot be discovered by surveillance measures such as border screening [11, 15]. By substituting I(t) obtained from the Eq. (4), the probability p _{ i,t } is
The first introduction of an infection that arrives on day t in country ith follows a geometric distribution with a timevarying parameter p _{ i,t }. The probability distribution is
where j = 1, 2, 3…t1. Due to different reporting details for case detections in various countries, the times of arrival are sometimes unavailable and only the days of illness onset will be reported during the early period of pandemic. When the time of arrival is unknown, we have assumed that it follows a multinomial distribution with a probability vector u = {u _{ k }}, k = k _{ min },…1, 0, 1,…k _{ max } (the minus numbers of days represent previous days with a minimum k _{ min } number of days, and the positive days represent following days with a maximum k _{ max } number of days) as the days from arrival to illness onset can be varied [15]. Given that the distribution is known, the adjusted probability is
Likelihood formulation
Given fixed parameters of γ, d, and θ, the likelihood for the arrival times of the first infected individuals for all countries is
where m = {m _{ 1 } (t), m _{ 2 } (t), m _{ 3 } (t), …, m _{ n } (t)} and t = {t _{ 1 } , t _{ 2 } , t _{ 3 } , …, t _{ n }} for n total number of countries that have reported their importing cases from the originating country. The estimate can be obtained by the maximum likelihood estimation (MLE).
Simulation testing
The simulation testing aimed to assess the performance of the proposed method by varying the levels of disease transmissibility and travel rates. The simulation tested the robustness of decreasing trend of travels and other parameter assumptions on the estimation. In each setting, 10,000 datasets were simulated for each settings of samples size (n) i.e. numbers of countries having imported cases. The population N was fixed at 1,000,000 people. The mean, standard deviation (SD), average standard errors of means (SE), and 95% credible intervals (Crl) of an estimate in a simulation were obtained to evaluate the performances.
Performance at different levels of reproduction numbers and travel rates
The estimation approach was evaluated at different plausible values of R _{ 0 } (1.2, 1.7, and 2.2) [16, 17]. We defined the low and high travel volumes by setting the daily rates of travel of each country (m) which followed uniform distributions (100, 1000) and (1000, 2000) respectively. The average duration of individual stay was fixed at 3 days. The length of infectious period was set as 3 days and the θ was fixed at 30%. Based on these settings, the occurrence of infections at day t could be simulated by Eq. (6) assuming a Bernoulli distribution, and the arrival times could be further obtained. By fitting the simulated data, the estimation performances at different levels of R _{ 0 } and travel rates were assessed.
Robustness against a decreasing trend of travel
In a pandemic outbreak, the travel rate could be reduced by time due to the risk perception of publics from mass media and international travel advice [18]. A scenario was designed to test the robustness of constant travelling rate assumption against the decreasing pattern of travels. Using the same scenario settings, we considered a linear decreasing trend of daily travel rates by assuming m _{ i } (t) = m _{ max } (1rt), where r is the daily dropping rate and m _{ max } is the passenger rate in the initial day. With reference to the information of pH1N1 pandemic, we assumed m _{ max } during an epidemic period followed a uniform distribution (200, 2500) with fixing r to be 0.4% and 0.8%.
Sensitivity analysis of parameter assumptions
In order to examine the sensitivity of parameter assumptions, the following distributional assumptions were considered:

1/γ ~ Gamma distribution with a mean of 3 days and 1 day SD

θ ~ Uniform distribution ranged from 10 to 50%

d ~ Uniform distribution ranged from 1 day to 7 days
The datasets were simulated based on these assumptions. Their impacts on the precision levels were drawn by fixing the values of parameters (1/γ = 3 days, θ = 30%, d = 3 days) in the estimation.
Numerical study: 2009 influenza pandemic a/H1N1
We applied the estimation approach to the Mexico pH1N1 in midMarch 2009. According to the World Health Organization (WHO), there had been an unexpected increase of influenzalike illness cases shown by the routine influenza surveillance in midApril, which was not the peak of an influenza season for normal outbreaks to take place, of 2009. Subsequent to a report revealing that two cases of an acute respiratory illness were discovered in two children living in the Southern California of the United States and further confirmed as infections of a new strain of H1N1 virus, additional cases were soon discovered in the country. By the end of April, WHO had further raised the pH1N1 alert to Phase Five, although both the general public and the governments still lacked adequate knowledge on the early stages of pH1N1 at that time. Several years later, there was an estimate of around 284 thousands deaths associated with the 2009 pH1N1 [19].
According to the National Council for Population of Mexico, we fixed the population of Mexico (N) as 106,682,518 in 2009 [20]. The start date of the pandemic outbreak was set as March 14, 2009 [21]. The arrival time and illness onset time of infected cases seeded by Mexico, and the number of passengers were retrieved from the studies of Fraser et al. [9], Balcan et al. [12], other articles [22–24] and some news from press [25, 26] with evidence that cases were imported from Mexico. Given the number of passengers from Mexico [9], we assumed the daily rate of passengers: 1. distributed uniformity (m _{ i } (t) = m _{ i }) and, 2. linearly decreased (m _{ i } (t) = m _{ max } (1rt), where r is the daily dropping rate and m _{ max } is the passenger rate in the initial day with the total number of passengers during the epidemic period held fixed. The dropping rate was assumed to be 40% divided by the length of the study period [18]. For those cases with missing arrival dates, the probability of arrival was adjusted with the illness onset date adapting the information from Sakaguchi et al. [15] in Eq. (8) i.e. u = {u _{ −6 } ,…, u _{ 4 } } = {0.007,…,0.014}. The information for the estimation is listed in Table 1.
In addition, the length of the infectious period was fixed at 3 days and the proportion of asymptomatic infections was fixed at 30% [13]. The average duration of stay of travellers was fixed as 3 days. A sensitivity analysis for the parameters was also conducted. The sensitivity of the estimate was tested by setting the infectious duration equal to 2 days and 4 days, θ as 10% and 50%, and d as 3 days and 5 days respectively.
The simulation was performed on R 2.15.2 [27] and the code can be provided upon requested.
Results
Simulations
Table 2 summarizes the estimation performances at different levels of R _{ 0 } and travel rates. The estimates of the proposed method were close to the true values of different R _{ 0 } settings. Given n = 4, the 95% CrI ranged from 2.08 to 2.43 for an expected R _{ 0 } of 2.2, whereas the 95% CrI ranged from 1.18 to 1.25 for an expected R _{ 0 } of 1.2. The results showed that with a low transmissibility, small samples could still achieve a stable precision level, for example, the 95% CrIs were around 1.19 to 1.22 when n ≥ 8 in the setting of R _{ 0 } = 1.2. Higher values of R _{ 0 } were associated with a slightly lower precision of the estimates. The estimation performance did not show any obvious difference between travel volumes. Apart from that, the estimates from our proposed method were robust to plausible drops of travel rates (Table 3). It was because most cases probably arose in the early phase of pandemic in which changes of travel may not subject to any impact in such a short period of time. The robustness of parameter assumptions was also assessed (Fig. 1). As with other studies, R _{ 0 } estimation was moderately sensitive to the infectious duration (or the generation interval). Given a plausible distributional assumption, the bias increased and the precision levels decreased with the increase of R _{ 0 }. For example, the estimate was around 2.4 and the 95% CrI ranged from 1.7 to 3.7 for an expected R _{ 0 } of 2.2. Nevertheless, the estimation performance was not sensitive in the scenarios with lower values of R _{ 0 }. Besides, the simulation results were robust to the assumptions of average stay duration and the proportion of asymptomatic infections. Like most simulation scenarios, the performances were better for a low transmissibility setting.
Numerical study
By employing the estimation method, we analyzed the data of the initial outbreak of pH1N1 in Mexico. The estimation results are summarized in Table 4. Showed in the table, the estimation results only have slightly difference in the settings between the constant and linearly decreasing daily rate of passengers. The estimated \( {\widehat{R}}_0 \) was equal to 1.69 with a 95% confidence intervals (CI) (1.65, 1.73) in a general situation i.e. 1/γ = 3 days, d = 3 days, and θ = 30%. The results were robust to the parameter assumptions (Table 4). Overall the estimated \( {\widehat{R}}_0 \) s were in range from 1.41 to 2.00 for our tested ranges of parameters. In line with the simulation testing, the effects from d and θ were comparatively small. The impacts on the estimate were minor for θ = 10% (\( {\widehat{R}}_0=1.67 \)) and θ = 50% (\( {\widehat{R}}_0=1.73 \)). Apart from that, the length of infectious period contributed more impact on the \( {\widehat{R}}_0 \). The estimate \( {\widehat{R}}_0 \) was equal to 1.45 (95% CI: 1.42, 1.48) when the infectious duration was equal to 2 days, whereas \( {\widehat{R}}_0 \) was equal to 1.94 (95% CI: 1.88, 1.99) when the infectious duration was equal to 4 days.
Discussion
Apart from influenzalikeillness, serological survey, and syndromic data, travel data can be an alternative source to timely infer the transmissibility for a new influenza pandemic. Arrival times of cases seeded by the originating country can be regarded as an integration of the number of infected individuals in the originating country and their mobility networks connecting with other regions [28, 29]. By adapting such travel data as a proxy of the size of epidemic expansion, the transmissibility of the influenza virus can be estimated. The estimation based on this information can prevent the bias induced by underestimations of different surveillance data. Undoubtedly, numbers of cases reported by surveillance systems are usually underestimated due to sampling methods, underdetection, and changes of capacities and requirements over time [9, 13, 30]. Apart from that, routine seroprevalence studies that require laboratory resources and a considerably long sample collection time may no longer serve as a suitable monitoring approach during the initial outbreak of a pandemic [31]. In consideration that serological data can refine parameter estimates, a reliable way of using serial crosssectional serological data alongside with surveillance data to estimate infection rate can be adopted [8]. To address possible errors in the estimation, multifaceted surveillance measures shall be adopted, especially during the early stages of new influenza outbreaks.
In this study, we developed a likelihoodbased approach that employed the arrival times of cases to estimate R _{ 0 } for a new influenza pandemic. We also conducted a simulation to assess the performance of proposed method. The method was applied to the initial outbreak of influenza pH1N1 in Mexico. We showed that the estimated \( {\widehat{R}}_0 \) was equal to 1.69 with a 95% CI (1.65, 1.73). The estimated \( {\widehat{R}}_0 \) and the corresponding range of sensitivity were consistent to other findings of pH1N1, mentioned in the review of R _{ 0 } from Boëlle et al., the R _{ 0 }s of the 2009 pH1N1 were in range from 1.2 to 2.3 [16]. Compared with other similar studies, Fraser et al. developed a Bayesian method to draw a posterior distribution of R _{ 0 } for the pH1N1 in Mexico by assuming that tourists infected at a rate proportional to the density of tourists per local resident [12]. The approach adopted the number of infected travelers within a fixed period rather than using the times of exporting infections to countries. Although their estimate (posterior median R _{ 0 } = 1.4) was similar to us, time to event data combining the mobility and pattern of epidemic invasion is usually preferred to the count data. Balcan et al. used a MLE approach to fit the arrival time data to the simulated epidemics through a highly parameterized metapopulation model in which the R _{ 0 } was their parameter of interest [12]. Their estimate (\( {\widehat{R}}_0=1.75 \)) was also closer to our finding. Although their transmission model required many details for parameterization and was computational intensive, it provides an estimate closer to the biological realism over a large population. Nevertheless, we believe our method provides a comparatively simpler way for the transmissibility estimation.
The reliability of the proposed method depends on the quality and quantity of travel data available during the early stage of a influenza pandemic. On estimating the transmissibility, data can be closely aligned with the proposed estimation method when provided with a timely collection of information during the early outbreak. Yet, with this proposed method based on a variety of data sources from the originating country and other regions, differences and incompatibility of surveillance systems across countries, in addition to disparate policies on international reporting and collaboration, have posed challenges to the acquisition of largescale samples. In pH1N1 of 2009, for instance, not many countries have reported their confirmed cases with known travel history to Mexico [9, 12, 13]. For the betterment of public health in the future, improvements shall be made on the coordination and technical innovation of the streamline or even the centralized surveillance systems for monitoring the imported cases of infections among countries.
One of the advantages of using likelihoodbased approach lies in the flexibility in incorporating other available data source. For example, the likelihood can potentially adapt the distribution of days from arrival to illness onset given that the information is available. Moreover, incorporation of demographic data in the transmission model will make the likelihood more informative, as, for instance, younger individuals are more likely to be infected by pH1N1. Yet, the reliability of such extension of the method shall be further justified by sufficient data.
There are several caveats for our study. We did not capture the effects from multistep journeys in the model although previous articles revealed a single and multistep travels did not differ much [32]. Moreover, no adjustment for disease transmissions on aircrafts was incorporated in the estimation. Nevertheless, a retrospective cohort study indicated inflight transmissions were unlikely [33]. As with many other approaches, we assumed homogeneous mixings at a country level and same risks of infections between travelers and the domestic population. In the early stage of a pandemic, stochastic effects usually induce spatial variations [9]. Apart from that, since airports are the main hubs of connection between countries, passengers will be more likely to be infected in the areas near the airport than areas far from the airports especially for an infectious disease with a short generation interval. Some cases would even recover before their travels. Further developments can take these issues into a consideration.
Conclusions
The study presents how to use travel data in an influenza pandemic to estimate the basic reproduction number (R _{ 0 }), a key parameter to determine what level of control measures should be used. Compared with other data sources, travel data is relatively more complete and timely for a new pandemic outbreak. Our approach as well as the estimate is potential to assist officials in planning control and prevention measures. Improved coordination as well as centralizing surveillance of imported cases among various regions would thus be beneficial to global health.
Abbreviations
 CI:

Confidence interval
 CrI:

Credible intervals
 MLE:

Maximum likelihood estimation
 pH1N1:

Pandemic A/H1N1
 SD:

Standard deviation
 SE:

Mean standard error
 SIR:

SusceptibleInfectiousRecovered model
 WHO:

World Health Organization
References
 1.
Kermack WO, McKendrick AG. Contributions to the mathematical theory of epidemics. Part I. Proc R Soc London Ser A. 1927;115:700–21.
 2.
Diekmann O, Heesterbeek JAP. Mathematical epidemiology of infectious diseases: model building, analysis and interpretation. New York: Wiley; 2000.
 3.
Wallinga J, Teunis P. Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures. Am J Epidemiol. 2004;160:509–16.
 4.
White LF, Pagano M. A likelihoodbased method for realtime estimation of the serial interval and reproductive number of an epidemic. Stat Med. 2008;27:2999–3016.
 5.
Chowell G, Nishiura H, Bettencourt LM. Comparative estimation of the reproduction number for pandemic influenza from daily case notification data. J R Soc Interface. 2007;4(12):155–66.
 6.
Hens N, Van Ranst M, Aerts M, Robesyn E, Van Damme P, Beutels P. Estimating the effective reproduction number for pandemic influenza from notification data made publicly available in real time: a multicountry analysis for influenza A/H1N1v 2009. Vaccine. 2011;29:896–904.
 7.
Glass K, Kelly H, Mercer GN. Pandemic influenza H1N1: reconciling serosurvey data with estimates of the reproduction number. Epidemiol. 2012;23:86–94.
 8.
Wu JT, Ho A, Ma ESK, Lee CK, Chu DKW, et al. Estimating infection attack rates and severity in real time during an influenza pandemic, analysis of serial crosssectional serologic surveillance data. PLoS Med. 2011;8:e1001103.
 9.
Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, et al. Pandemic potential of a strain of influenza a (H1N1), early findings. Science. 2009;324:1557–61.
 10.
Yang W, Lipsitch M, Shaman J. Inference of seasonal and pandemic influenza transmission dynamics. Proc Natl Acad Sci. 2015;112(9):2723–8.
 11.
Fielding R. Pandemic influenza A/H1N1 (pH1N1) in Hong Kong: anatomy of a response. NTS Working Paper Series.2010;3.
 12.
Balcan D, Hu H, Goncalves B, Bajardi P, Poletto C, Ramasco JJ, et al. Seasonal transmission potential and activity peaks of the new influenza A (H1N1): a Monte Carlo likelihood analysis based on human mobility. BMC Med. 2009;7:45.
 13.
Tizzoni M, Bajardi P, Poletto C, Ramasco JJ, Balcan D, Gonçalves B, et al. Realtime numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm. BMC Med. 2012;10:165.
 14.
Khan K, Arino J, Hu W, Raposo P, Sears J, Calderon F, Heidebrecht C, Macdonald M, Liauw J, Chan A, Gardam M. Spread of a novel influenza A (H1N1) virus via global airline transportation. N Engl J Med. 2009;361(2):212–4.
 15.
Sakaguchi H, Tsunoda M, Wada K, Ohta H, Kawashima M, Yoshino Y, et al. Assessment of border control measures and community containment measures used in Japan during the early stages of Pandemic (H1N1) 2009. PLoS One. 2012;7:e31289.
 16.
Boëlle PY, Ansart S, Cori A, Valleron AJ. Transmission parameters of the A/H1N1 (2009) influenza virus pandemic: a review. Influenza Other Respir Viruse. 2011;5:306–16.
 17.
Mills C, Robins J, Lipsitch M. Transmissibility of 1918 pandemic influenza. Nature. 2004;432:904–6.
 18.
Bajardi P, Poletto C, Ramasco JJ, Tizzoni M, Colizza V, Vespignani A. Human mobility networks, travel restrictions, and the global spread of 2009 H1N1 pandemic. PLoS One. 2009;6:e16591.
 19.
Dawood FS, Iuliano AD, Reed C, Meltzer MI, Shay DK, Cheng PY, et al. Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study. Lancet Infect Dis. 2012;12:687–95.
 20.
National Council for Population of Mexico. National Survey of Demographic Dynamics 2009. http://www.conapo.gob.mx/es/CONAPO/ENADID_2009. Accessed 27 Jan 2017.
 21.
The Ministry of Health of Mexico. Distribution of confirmed cases according to date of onset of symptoms. In: Current status of the epidemic. 2009. http:// portal.salud.gob.mx/descargas/pdf/influenza/situacion_ actual_epidemia_121009.pdf. Accessed 29 Oct 10.
 22.
CastroJiménez MÁ, CastilloPabón JO, ReyBenito GJ, PulidoDomínguez PA, BarbosaRamírez J, VelandiaRodriguez DA, AnguloMartínez ES, on behalf of the Virology Group and the Communicable Diseases Surveillance Group. Epidemiologic analysis of the laboratoryconfirmed cases of influenza A(H1N1)v in Colombia. Euro Surveill. 2009;14:19284.
 23.
New influenza A(H1N1) investigation teams. New influenza A(H1N1) virus infections in France, April  May 2009. Euro Surveill. 2009;14:19221.
 24.
Lisena F, Bordi L, Carletti F, Castilletti C, Ferraro F, Lalle E, et al. Influenza A (H1N1) in Rome, Italy in family: three case reports. Cases J. 2009;2:9123.
 25.
Bustamante J. Argentina confirms first H1N1 flu case. 2009. Available from, http://www.reuters.com/article/2009/05/07/usfluargentinaidUSTRE5468RC20090507. Accessed 29 Oct 2010.
 26.
Information Services Department. 1st H1N1 case confirmed in HK. 2009. http://archive.news.gov.hk/isd/ebulletin/en/category/healthandcommunity/090501/html/090501en05004.htm. Accessed 29 Oct 2010.
 27.
R Development Core Team. R, A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2006.
 28.
Chong KC, Zee BCY. Modeling the impact of air, sea, and land travel restrictions supplemented by other interventions on the emergence of a new influenza pandemic virus. BMC Infect Diseases. 2012;12:309.
 29.
Chong KC, Fong HF, Zee CY. Estimating the incidence reporting rates of new influenza pandemics at an early stage using travel data from the source country. Epidemiol Infect. 2014;142:955–63.
 30.
Lipsitch M, Hayden FG, Cowling BJ, Leung GM. How to maintain surveillance for novel influenza A H1N1 when there are too many cases to count. Lancet. 2009;374:1209–11.
 31.
Achonu C, Rosella L, Gubbay JB, Deeks S, Rebbapragada A, Mazzulli T, et al. Seroprevalence of pandemic influenza H1N1 in Ontario from January 2009May 2010. PLoS One. 2011;6:e26427.
 32.
Epstein JM, Goedecke DM, Yu F, Morris RJ, Wagener DK, Bobashev GV. Controlling pandemic flu: the value of international air travel restrictions. PLoS One. 2007;2:e401.
 33.
Baker MG, Thornley CN, Mills C, Roberts S, Perera S, Peters J, et al. Transmission of pandemic A/H1N1 2009 influenza on passenger aircraft: retrospective cohort study. BMJ. 2010;340:c2424.
Acknowledgements
The authors would like to thank Katherine Min Jia for the revision on the final proof. The authors are grateful to the financial support from National Natural Science Foundation of China, Research Grants Council of the University Grants Committee, Funding of New Drug Clinical Evaluation Research Technology Platform, Health and Medical Research Fund (Project reference no: CU16C11) of Food and Health Bureau of the Hong Kong Special Administrative Region Government, Shenzhen Municipal Science and Technology and the Chinese University of Hong Kong. All the abovementioned funders have no role in study design, collection, analysis, interpretation of data and in the writing of the manuscript and the decision to submit the manuscript for publication.
Funding
The research is supported by National Natural Science Foundation of China [31401124, 81201817, 81473035]; General Research Fund of Research Grants Council of the Hong Kong Special Administrative Region Government [476013]; The Chinese University of Hong Kong Direct grant [2014.1]; Health and Medical Research Fund (Project reference no: CU16C11) of Food and Health Bureau of the Hong Kong Special Administrative Region Government; Funding of New Drug Clinical Evaluation Research Technology Platform [2012ZX090312002]; Shenzhen Municipal Science and Technology Program  Basic research project (JCYJ20150630165236958). The abovementioned funders had no role in the study design, data collection, data analysis, interpretation of results, and manuscript preparation.
Availability of data and materials
Study data that support the results are available from the corresponding author upon request, while datasets generated as the central part of conclusion are included in this article.
Authors’ contributions
KCC designed the study, analyzed the data, and wrote the first draft of the manuscript. BCYZ and MHW designed the study. All authors reviewed and approved the final version.
Competing interests
All authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Basic reproduction number
 Epidemic models
 Travel data
 SIR model
 Influenza pandemic