Hostname: page-component-7c8c6479df-ws8qp Total loading time: 0 Render date: 2024-03-28T14:52:19.866Z Has data issue: false hasContentIssue false

Estimating the incidence reporting rates of new influenza pandemics at an early stage using travel data from the source country

Published online by Cambridge University Press:  10 October 2013

K. C. CHONG*
Affiliation:
Division of Biostatistics, Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong SAR, China
H. F. FONG
Affiliation:
Center for Global Public Health, University of California, Berkeley, CA, USA
C. Y. ZEE
Affiliation:
Division of Biostatistics, Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong SAR, China
*
*Author for correspondence: Dr. K. C. Chong, Division of Biostatistics, Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong SAR, China. (Email: marc@cct.cuhk.edu.hk)
Rights & Permissions [Opens in a new window]

Summary

During the surveillance of influenza pandemics, underreported data are a public health challenge that complicates the understanding of pandemic threats and can undermine mitigation efforts. We propose a method to estimate incidence reporting rates at early stages of new influenza pandemics using 2009 pandemic H1N1 as an example. Routine surveillance data and statistics of travellers arriving from Mexico were used. Our method incorporates changes in reporting rates such as linearly increasing trends due to the enhanced surveillance. From our results, the reporting rate was estimated at 0·46% during early stages of the pandemic in Mexico. We estimated cumulative incidence in the Mexican population to be 0·7% compared to 0·003% reported by officials in Mexico at the end of April. This method could be useful in estimation of actual cases during new influenza pandemics for policy makers to better determine appropriate control measures.

Type
Original Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
The online version of this article is published within an Open Access environment subject to the conditions of the Creative Commons Attribution license .
Copyright
Copyright © Cambridge University Press 2013

INTRODUCTION

During the early outbreak of an influenza pandemic, rapid disease transmission can lead to exponential rises of influenza cases throughout the population. Underreporting of influenza cases in early stages poses problems in estimating both pandemic severity and transmission intensity. Underreporting stems from the short infectious periods of influenza infections; thus, individuals may recover before seeking treatment from their healthcare provider or before being tracked in a surveillance system. Asymptomatic or mild cases may not even be reported at all. A previous study has shown that official surveillance only reveals a small proportion of actual infections during influenza pandemics [Reference Xiaoli1] – in some instances the consultation rate in influenza-like illness case-patients was no more than 50%. Furthermore, cases increase exponentially during the initial stage of an outbreak and the limited capacity of surveillance systems, such as limited serological tests, can also lead to underreporting [Reference Colizza2].

Underreporting has consequential effects on public health response. From a policy perspective, underreporting can lead to officials underestimating public health risk which in turn affects planning and the implementation of systematic control and prevention activities. For example, there may be a delay in implementing entry screening for travellers or inadequate warning to local and national health departments. The underestimation of incidence and pandemic severity can also reduce education and health notices to the general public about the influenza virus, causing the public not to take measures to protect themselves through vaccines, hand washing, or other control measures. On the other hand, the case-fatality rate would be overestimated as being higher than it actually was due to the missing calculation of asymptomatic and mild cases from the rate denominator [Reference Fielding3]. If there is an insufficient system for pandemic control, this situation can place unexpected, unnecessary financial and human resource demands on a healthcare system. Therefore, reliable methods to estimate the reporting rate during early influenza epidemic outbreaks are critical to good public health and infectious disease response systems.

The influenza A(H1N1) pandemic in Mexico in mid-March 2009 is one example of when country officials underestimated influenza incidence rates. Although it was not a peak season for an influenza outbreak, routine influenza surveillance identified an unexpected increase in cases of an influenza-like illness in mid-April 2009 [4]. An acute respiratory illness was discovered in two children and further confirmed as a new strain of H1N1 virus. Subsequently, on 26 April 2009, the World Health Organization (WHO) notified the public of the new H1N1. Additional cases were soon discovered in the USA [5], and the WHO had raised the H1N1 pandemic alert level to phase 5 by the end of April. At the time, governments and the public still lacked sufficient knowledge about the early stages of the outbreak. At this time, according to H1N1 surveillance data from the Ministry of Health in Mexico, cumulative incidence was measured to be as low as 0·003% in Mexico's population [6].

Due to the increasing awareness of H1N1 throughout April 2009, other at-risk countries began control measures at border points of entry to prevent local epidemics. For example, thermal screening was implemented and suspected cases with a travel history to Mexico were monitored and some quarantined [Reference Fielding3]. Because surveillance at the borders was quite thorough for influenza-like illness cases even before the H1N1 virus had spread globally, data on early cases such as time of import from the source country is relatively more complete and timely than other available data. For estimating the size and local expansion of the influenza pandemic, this is a valuable data source that also provides a perspective on how the disease is spread.

Previous studies have demonstrated the usefulness of mathematical modelling in summarizing the epidemiology of infectious illness and in examining impact of the diseases from the external factors [Reference Ferguson7Reference Medlock and Galvani14]. In this study, a mathematical modelling approach was adopted to develop a method to help quantify the spread of infectious disease in the population. The method is able to estimate the incidence reporting rate by using the local routine surveillance data with estimates refined from statistics of travellers from the source country for an influenza pandemic. The approach made use of the 2009 pandemic influenza A(H1N1) (pH1N1) outbreak as an example.

METHODS

Mathematical model

We adopted a susceptible-exposed-infectious-recovered (SEIR) model to describe the dynamic system of the infectious disease [Reference Kermack and McKendrick15]. For each time point t, a whole population is classified into one of four groups (‘compartments’): susceptible [S(t)]; exposed [E(t)]; infectious [I(t)]; or recovered [R(t)]. Using S, E, I, and R to represent each compartment, the SEIR model has four differential equations describing the rates of subject movement for each time step:

$$ \eqalign{& \displaystyle{{{\rm d}S} \over {{\rm d}t}} = - \beta SI = - \beta (S_{\rm L} + S_{\rm T} )I, \cr & \displaystyle{{{\rm d}E} \over {{\rm d}t}} = \beta SI - \alpha E, \cr & \displaystyle{{{\rm d}I} \over {{\rm d}t}} = \alpha E - \gamma I, \cr & \displaystyle{{{\rm d}R} \over {{\rm d}t}} = \gamma I,} $$

In this compartmental model, once a susceptible individual (including local residents and travellers) in compartment S(t) is infected, they move to compartment E(t) and remain there for the latent period. When that latent period is over, they move to compartment I(t) during the infectious period. When the infectious period is over, individuals in compartment I(t) recover and move to compartment R(t). S L is the local susceptible size and S T is the number of travellers from the source country. As S T is far smaller than S L i.e. S LS T, we approximate

$$ - \beta SI = - \beta (S_{\rm L} + S_{\rm T} )I \approx - \beta S_{\rm L} I.$$

In the model, the probability of an individual becoming infected is configured using the basic reproduction number (R 0), the average number of secondary infections produced by a typical infectious individual in a wholly susceptible population. The transmission rate is β, so the force of infection is βI. The total population size N, is equal to S + E + I + R for any time and N = S for time zero. We assumed the lengths of the latent period and the infectious period follow exponential distributions and their averages would be 1/α and 1/γ, respectively. Adopting the linearization method [Reference Diekmann and Heesterbeek16], the basic reproduction number R 0 is equal to βN/γ.

Parameter estimation

In our model, we assumed homogenous mixing between individuals in the system being studied, and that cases reported to officials when infectious. Since the numbers of asymptomatic and non-severe cases may not have been presented for the observed surveillance time-series data U(t), we used f t (.) to represent a functional form of reporting rates. Therefore, f t (.) is defined as (reported cases/actual cases), f t (.) = U(t)/αE. The αE is generated from the SEIR model. We considered two forms of f t (.) in the estimation:

  1. (1) Constant reporting rate: f t (r= r,

  2. (2) Linearly increasing reporting rate:

$$\eqalign{&f_t (r_{\min}, r_{\max} ) \cr&\quad= \left\{ \matrix{r_{\min}, \quad t \lt t_0 \cr r_{\min} + (r_{\max} - r_{\min} )\left( {\displaystyle{{t - t_0} \over {t_1 - t_0}}} \right),\quad t_0 \leqslant t \lt t_1 \cr r_{\max}, \quad t \geqslant t_1 } \right.}$$

For parameter estimation, we first iterated the parameters by fixing their values within a grid search. Thus the reproduction number R 0 can be fitted into the SEIR model using the least-squares method. We then adopted the earliest times of infected cases imported from Mexico (T i ) and the daily rate of travel (m i ) to particular country i. Assuming the travelling cases had the same daily risk of exposure to the influenza virus as local cases, the average imported cases to a country i will be βm i I(t) for time t with m i the daily rate of travel using the fitted SEIR model. Assuming a Poisson event, we assume the probability of importing at least one case from the source country at time t as p i,t  =  (1 − q i )(1−exp[−βm i I(t)]), where q i is the entry screening sensitivity for case detection of country i. Therefore, the estimated time of the first imported case can be simulated as $\hat T_i = \sum\nolimits_{k = 1} {kp_{i,k} \prod\nolimits_{j \lt k} {(1 - p_{i,j} )}} $$ . Iterations were repeated for ranges of fixed values within the grid search. Optimum parameters in f t (.) were obtained with the minimum square root of the sum of standardized squared errors (RSE) between observed data (T i ) and the simulated estimate times of the first cases imported ( $\hat T_i $$ ) from Mexico, i.e.

$${\rm RSE} = \sqrt {\sum\limits_i {\displaystyle{{(\hat T_i - T_i )^2} \over {\hat T_i}}}}. $$

We also developed a bootstrap method to calculate 95% confidence intervals (CIs). Supposing pairs of resample (T i *, m i *) were randomly drawn from the original pairs of T i and m i with replacement, the bootstrapped RSE for the jth iteration of bootstrapping was:

$${\rm RSE}^{\left[ j \right]} = \sqrt {\sum\limits_i {\displaystyle{{(\hat T_i ^* - T_i ^*)^2} \over {\hat T_i ^*}}}}. $$

One thousand bootstrapped RSE[j] was generated with corresponding fitted parameters. The 2·5th and 97·5th percentiles of the fitted parameters were the lower and upper limits, respectively, of the non-parametric 95% CI over the 1000 samples.

Given the estimate R 0, we back-calculated the exponential growth rate of the pandemic [Reference Chowell17]:

$$R_0 = \displaystyle{{\theta ^2 + \theta (\alpha + \gamma )} \over {\alpha \gamma}}, $$

where θ is the exponential growth rate. The date of pH1N1 seeding can be calculated as

$$\eqalign{\left({{\rm date \,of\, first}\atop\, {\rm confirmed\, case}}\right) - {\rm no}.\ {\rm of}\left( {\displaystyle{{{\rm ln}\left\lfloor {U(1)/f_1 (.)} \right\rfloor} \over \theta}} \right){\rm days,}$$

assuming an exponential growth during the early phase of the pandemic.

The estimation method was implemented using SAS v. 9.2.1 software (SAS Institute Inc., USA).

Materials and parameter values

The population of Mexico (N) was 106 682 518 in 2009, a figure provided by the National Council for Population of Mexico [18]. The pH1N1 surveillance data [U(t)], shown in Figure 1, was obtained from the Ministry of Health of Mexico covering the first wave of the pandemic from 14 March 2009 to 27 May 2009 [6]. We assumed the reporting rate remained constant throughout the time period and increased linearly. In the linear increase approach, the start date (t 0) of enhanced surveillance in Mexico was 17 April 2009 [19] and the end date (t 1) was 17 May 2009.

Fig. 1. Confirmed cases in Mexico between 14 March 2009 and 27 May 2009.

Traveller data including earliest dates of infected cases arriving from Mexico into different countries are shown in Table 1. We estimated the daily rates of travel to a particular country i (m i ) by dividing the passenger count in March 2009 and April 2009 by 61 days. We excluded the USA from our study because air travel is not the only means of cross-border transport between the two countries.

Table 1. Number of travellers and earliest date of cases imported from Mexico to a particular country in March and April, 2009

The epidemiological details in the parameter estimation were mostly from the previous findings of pH1N1. The lengths of the latent and infectious period S were set at 1·6 days and 1·4 days, respectively [Reference Yang20Reference Boëlle24].

Sensitivity analysis

Limited entry screening at airports at the initial stage of H1N1 could have led to undetected cases from Mexico in the early stages [Reference Chong and Zee25Reference Nishiura and Kamiya27], especially since Mexico did not implement exit screening. In our model, in order to consider undetected cases, we tested results with a range of entry-screening sensitivities. As exact entry-screening sensitivities would vary for all countries, we varied the screening sensitivities by uniformly choosing from 30% to 100% for each of the countries in every simulation.

We also performed a multivariate sensitivity analysis on the lengths of the latent and infectious periods. The latent period was assumed to follow a gamma distribution with a mean of 1·6 days and a standard deviation of half a day; the infectious period followed a gamma distribution with a mean of 1·4 days and a standard deviation of half a day.

Parameter distributions were drawn from 1000 simulations.

RESULTS

From our model, the value of the constant reporting rate (r) was estimated at 0·46% using a minimum value of RSE. The bootstrapped 95% CI was between 0·28% and 0·69% when the estimated value of R 0 was 1·24 (Fig. 2, Table 2). The value R 0 remained steady when r was >0·1%. The figure demonstrated that an increasing reporting rate was associated with exponential decreases in R 0; thus, only fitting the surveillance data to the epidemic model would provide unreliable findings for the estimation of the r. Using these estimates, there was a 0·7% (95% CI 0·4–1·1) cumulative incidence in the Mexican population at the end of April 2009, which was the time that the pandemic phase 5 alert level was announced by the WHO.

Fig. 2. Values of the minimum square root of the sum of standardized squared errors (RSE) and R 0 given different constant r.

Table 2. Estimates of reporting rates (bootstrapped 95% confidence intervals) given different variations

The reporting rate did not increase after enhanced surveillance in Mexico after mid-April 2009, when officials stepped up surveillance systems (Table 2). The r min was 0·46% (bootstrapped 95% CI 0·27–0·68), whereas the r max was 0·47% (bootstrapped 95% CI 0·28–0·69). Reporting behaviour may not have been significantly affected during this short time-frame.

In the study, we considered the sensitivity of missing ‘detections’ of imported cases in the estimation process. The entry-screening sensitivities of countries were found to be moderately sensitive to our results. If the entry-screening sensitivities were distributed uniformly between 30% and 100%, the constant r was estimated as 0·18% (95% CI 0·09–0·31) (Table 2). The value was relatively lower due to a decrease in the average probability of detection. If a linear trend was assumed, a slight increase of the reporting rate was observed. The rate increased from 0·10% (95% CI 0·03–0·26) to 0·31% (95% CI 0·11–0·80). However, this increasing range was insignificant and did not deviate much from our initial estimates.

The impacts of variation of latent period length and infectious period on our results were also tested. As shown in Figures 3 and 4, the variation of lengths did not make any impact on the reporting rate estimation. The constant r was 0·44% (bootstrapped 95% CI 0·31–0·69) and the values of the linear reporting rates (r min and r max) were both close to this value (Table 2). However, the variations did affect the estimated value of R 0; the estimated median R 0 was 1·24 (range 1·11–1·44) given a constant assumption of r. Insignificant difference was observed for the range of R 0 given a linearly increasing assumption. The range of R 0 was consistent with other studies [Reference Yang20, Reference Fraser28, Reference Hsieh29].

Fig. 3. The effect of variations from the length of the latent period (∼gamma[mean = 1·6, s.d. = 0·5]) and the length of the infectious period (∼gamma[mean = 1·4, s.d. = 0·5]) given a constant reporting rate assumption. Left panel is the box-plot of r and the right panel is the box-plot of R 0.

Fig. 4. The effect of variations in the lengths of the latent and infectious periods given a linearly increasing reporting rate assumption. The impacts to r min, r max, and R 0 are shown by the box-plots from left, middle, and right panels, respectively.

Given the estimates, the date of seeding for pH1N1 was exponentially interpolated from the date of the first confirmed case (i.e. 14 March 2009 with an estimate of about 400 infections). From the results, the date of seeding for pH1N1 was estimated as 24 December 2011 (95% CI 17–29 December 2011) in order to maintain a sufficient large epidemic size for exportation of cases. About 5500 Mexicans were infected by pH1N1 virus before the date of the first confirmed case in the surveillance data. When the minimum and maximum values (1·11 and 1·44, tespectively) of the range of R 0 were adopted in the estimation, the dates of seeding for pH1N1 were 26 September 2011 (95% CI 12 September 2011 to 7 October 2011) and 27 January 2012 (95% CI 23–30 January 2012), respectively.

DISCUSSION

A reliable method to estimate reporting rates during early phases of a new influenza pandemic is critical in addressing infectious disease response in the 21st century, especially with increased travel by air, land, and sea [Reference Chong and Zee25]. The importance of this was highlighted in Mexico's 2009 influenza pandemic, in which the reported incidence by Mexican officials (0·003%) during the early stages of the outbreak was not even close to our estimate during the early outbreak. Even though the strain of pH1N1 virus had been further confirmed in mid-April 2009, the Mexican officials' reporting rates still did not increase. This situation masked the actual growth of pH1N1, leading to a reduction of public awareness and potentially more rapid disease transmission. Inaccurate estimates overstating the risk can provide misleading information to the public and potentially raise levels of anxiety or panic [Reference Leung30].

A reliable estimate can assist officials at local, national, and global levels in planning and implementing prevention and control strategies for a pandemic influenza during the early stages, and better inform policy and protocols for other infectious disease outbreaks. In our study, we introduced such a method using existing information available to countries during a pandemic, the time at which imported cases may be arriving from a source country, to estimate reporting rates. According to our results, the estimated epidemic size was larger than officially reported in 2009; we found an estimate of 0·7% cumulative incidence (about 691 000 individuals) in the Mexican population compared to the 0·003% reported from the Ministry of Health of Mexico [6]. In terms of the epidemic size, our estimates were in line with other studies [Reference Colizza2, Reference Lipsitch31] but were higher than that of Fraser et al. [Reference Fraser28]. The reason for the difference is that our approach adopted time-series data for reported incidence, which can help better validate results when using traveller data. Several studies have employed a cross-sectional set of travel data to estimate actual epidemic size but those approaches did not aim to project the epidemic curve or address trends of reporting behaviour.

Interpolated estimation suggested the date of initiation for pH1N1 was late December 2008, which agrees well with other studies [Reference Fraser28] and, suggests that the pH1N1 virus had the potential to spread to other continents prior to laboratory confirmation of the virus [Reference Hsieh29]. With the use of the SEIR model and the estimates, we were also able to estimate that around 0·005% of the Mexican population was infected prior to the first case being detected by the surveillance system. Therefore, there is a possibility that undetected cases from Mexico, in other countries before the first global case was reported, could have affected our estimates. By using the mathematical model, the probability of having imported case from Mexico for at least one listed country (Table 1) was about 0·21 (results not showed) prior to 14 March 2009. Hence, the early ‘missing’ detection of imported cases from Mexico was not unexpected. This situation has similar potential issues with entry-screening sensitivities mentioned previously in the Results section, and we believe it would have only a minor effect on our findings.

The reliability of our proposed method would depend greatly on the quantity and quality of travel surveillance available at the borders during the early stages of a potential pandemic. If surveillance data from travellers could be collected in a timely way, it could effectively align with the estimation of a new influenza pandemic size and threat. However, there are challenges in acquiring large samples because countries especially those which do not border each other have different and incompatible surveillance systems as well as disparate policies on international reporting and collaboration. In our study, we only found 12 countries that reported their confirmed cases with known travel history in Mexico. Regarding this issue of small sample size, a bootstrap method was the preferred choice. In the future, improved coordination and technical innovations to streamline or even centralize infectious disease surveillance of travellers between countries would be beneficial to public health.

Besides the surveillance data from travellers, routine serological surveys could be another source of estimating incidence. However, compared to the surveillance data at borders, routine samples of seroprevalence may not be suitable during an initial outbreak of a pandemic as it requires laboratory resources and a longer collection time [Reference Achonu32]. Its reliability also relies on the sampling frame of the data [Reference Wu33]. Using serial cross-sectional serological data along with surveillance data could be reliable in estimating infection rates, since serological data could refine parameter estimates [Reference Wu33]. In order to account for possible estimation errors, multi-faceted surveillance measures are recommended, especially for new outbreaks of influenza pandemics during the early stages.

One of the advantages of using our method is its flexibility in adapting to/incorporating other epidemic models. It can be extended using similar concepts which adapt the reporting rate function of incidence in the epidemic models. For example, our approach could potentially be extended to demographic stratified models. As younger age groups were likely to be affected by pH1N1 and to be presented in ascertainment, incorporation of demographic stratified models would make the modelling results more realistic. However, sufficient data is required to support the extension of the method.

One of the caveats for applying the method to pH1N1 in Mexico was the homogenous dispersion of infections throughout the source country [Reference Lipsitch31]. Clearly, the pH1N1 outbreak may have not yet spread to all cities in Mexico at the early stage. Without available infection data at the city level, the resolution of our results would not be high enough and the spatial variation would alter our estimates. Although our method provides further understanding on how to tackle estimates of incidence reporting rates at early stages of an influenza outbreak, future studies could explore further model extensions.

DECLARATION OF INTEREST

None.

References

REFERENCES

1. Xiaoli, , et al. Estimates of the true number of cases of pandemic (H1N1) 2009, Beijing, China. Emerging Infectious Diseases 2010; 16: 17861788.Google Scholar
2. Colizza, V, et al. Estimate of novel influenza A/H1N1 cases in Mexico at the early stage of the pandemic with a spatially structured epidemic model. PLoS Currents 2009; 1: RRN1129.Google Scholar
3. Fielding, R. Pandemic influenza A/H1N1 (pH1N1) in Hong Kong: anatomy of a response. NTS Working Paper Series, 2010; 3.Google Scholar
4. World Health Organization. WHO ad hoc scientific teleconference on the current influenza A(H1N1) situation (http://www.who.int/csr/resources/publications/swineflu/tc_report_2009_04_29/en/index.html). Accessed 19 April 2010.Google Scholar
5. Centers for Disease Control and Prevention. Swine influenza A (H1N1) infection in two children – Southern California, March–April 2009. Morbidity and Mortality Weekly Report 2009; 58 (Dispatch): 13.Google Scholar
6. The Ministry of Health of Mexico. Current status of the epidemic, 12 October 2009. Distribution of confirmed cases according to date of onset of symptoms (http://portal.salud.gob.mx/descargas/pdf/influenza/situacion_actual_epidemia_121009.pdf). Accessed 19 April 2010.Google Scholar
7. Ferguson, N, et al. Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature 2005; 437: 209214.CrossRefGoogle ScholarPubMed
8. Wu, JT, et al. Reducing the impact of the next influenza pandemic using household-based public health interventions. PLoS Medicine 2006; 3: e361.CrossRefGoogle ScholarPubMed
9. Vynnycky, E, Edmunds, W. Analyses of the 1957 (Asian) influenza pandemic in the United Kingdom and the impact of school closures. Epidemiology and Infection 2008; 136: 166179.Google Scholar
10. Wu, JT, et al. School closure and mitigation of pandemic (H1N1) 2009, Hong Kong. Emerging Infectious Diseases 2010; 16: 538541.Google Scholar
11. Gani, R, et al. Potential impact of antiviral drug use during influenza pandemic. Emerging Infectious Diseases 2005; 11: 13551362.Google Scholar
12. Longini, IJ, et al. Containing pandemic influenza with antiviral agents. American Journal of Epidemiology 2004; 159: 623633.CrossRefGoogle ScholarPubMed
13. Riley, S, Wu, JT, Leung, GM. Optimizing the dose of pre-pandemic influenza vaccines to reduce the infection attack rate. PLoS Medicine 2007; 4: e218.Google Scholar
14. Medlock, J, Galvani, AP. Optimizing influenza vaccine distribution. Science 2009; 325: 17051708.CrossRefGoogle ScholarPubMed
15. Kermack, WO, McKendrick, AG. Contributions to the mathematical theory of epidemics, part I. Proceedings of the Royal Society of London, Series A 1927; 115: 700721.Google Scholar
16. Diekmann, O, Heesterbeek, JAP. Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. New York: Wiley, 2000.Google Scholar
17. Chowell, G, et al. The basic reproductive number of Ebola and the effects of public health measures: The cases of Congo and Uganda. Journal of Theoretical Biology 2004; 229: 119126.Google Scholar
18. National Council for Population of Mexico. Estimated population 2009 (http://www.conapo.gob.mx/index.php?option=com_content&view=article&id=125&Itemid=193). Accessed 19 April 2010.Google Scholar
19. Centers for Disease Control and Prevention. Outbreak of swine-origin influenza A (H1N1) virus infection – Mexico, March–April 2009. Morbidity and Mortality Weekly Report 2009; 58 (Dispatch): 13.Google Scholar
20. Yang, Y, et al. The transmissibility and control of pandemic influenza A (H1N1) virus. Science 2009; 326: 729733.CrossRefGoogle ScholarPubMed
21. Nishiura, H, Inaba, H. Estimation of the incubation period of influenza A (H1N1-2009) among imported cases: addressing censoring using outbreak data at the origin of importation. Journal of Theoretical Biology 2011; 272: 123170.Google Scholar
22. Lessler, J, Reich, NG, Cummings, DAT. Outbreak of 2009 pandemic influenza A(H1N1) at a New York City school. New England Journal of Medicine 2009; 361: 26282636.Google Scholar
23. Writing Committee of the WHO Consultation on Clinical Aspects of Pandemic (H1N1) 2009 Influenza. Clinical aspects of pandemic 2009 influenza A (H1N1) virus infection. New England Journal of Medicine 2010; 362: 17081719.Google Scholar
24. Boëlle, PY, et al. Transmission parameters of the A/H1N1 (2009) influenza virus pandemic: a review. Influenza and Other Respiratory Viruses 2011; 5: 306316.Google Scholar
25. 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 Infectious Diseases 2012; 12: 309.Google Scholar
26. Cowling, BJ, et al. Entry screening to delay local transmission of 2009 pandemic influenza A (H1N1). BMC Infectious Diseases 2010; 10: 82.CrossRefGoogle ScholarPubMed
27. Nishiura, H, Kamiya, K. Fever screening during the influenza (H1N1-2009) pandemic at Narita International Airport, Japan. BMC Infectious Diseases 2011; 11: 111.Google Scholar
28. Fraser, C, et al. Pandemic potential of a strain of influenza A(H1N1): early findings. Science 2009; 324: 15571561.Google Scholar
29. Hsieh, YH, et al. Early outbreak of 2009 influenza A(H1N1) in Mexico prior to identification of pH1N1 virus. PLoS ONE 2011; 6: e23853.Google Scholar
30. Leung, GM, et al. The impact of community psychological responses on outbreak control for severe acute respiratory syndrome in Hong Kong. Journal of Epidemiology and Community Health 2003; 57: 857863.Google Scholar
31. Lipsitch, M, et al. Use of cumulative incidence of novel influenza A/H1N1 in foreign travelers to estimate lower bounds on cumulative incidence in Mexico. PLoS ONE 2009; 4: e6895.Google Scholar
32. Achonu, C, et al. Seroprevalence of pandemic influenza H1N1 in Ontario from January 2009–May 2010. PLoS ONE 2011; 6: e26427.Google Scholar
33. Wu, JT, et al. Estimating infection attack rates and severity in real time during an influenza pandemic: analysis of serial cross-sectional serologic surveillance data. PLoS Medicine 2011; 8: e1001103.Google Scholar
34. Health Protection Agency West Midlands H1N1v Investigation Team. Preliminary descriptive epidemiology of a large school outbreak of influenza A(H1N1)v in the West Midlands, United Kingdom, May 2009. Eurosurveillance 2009; 14: pii=19264.Google Scholar
35. Influenza A(H1N1)v Investigation Teams. Modified surveillance of influenza A(H1N1)v virus infections in France. Eurosurveillance 2009; 14: pii=19276.Google Scholar
36. Reuters AlertNet. Colombia reports flu case, Ecuador monitors border. Released 3 May 2009 (http://www.alertnet.org/thenews/newsdesk/N03367692.htm). Accessed 28 January 2010.Google Scholar
37. Channel NewsAsia. El Salvador confirms first two cases of Influenza A(H1N1). Released 4 May 2009 (http://www.channelnewsasia.com/stories/afp_world/view/426728/1/.html). Accessed 28 January 2010.Google Scholar
38. Reuters. Argentina confirms first H1N1 case. Released 12 May 2009 (http://www.reuters.com/article/idUSN07423083). Accessed 28 January 2010.Google Scholar
39. Channel NewsAsia. Brazil confirms first domestically transmitted H1N1 flu case. Released 9 May 2009 (http://www.channelnewsasia.com/stories/afp_world/view/428001/1/.html). Accessed 28 January 2010.Google Scholar
40. Reuters. Cuba says Mexican has island's first H1N1 flu case. Released 12 May 2009 (http://www.reuters.com/article/idUSTRE54B0RZ20090512) Accessed 28 January 2010.Google Scholar
Figure 0

Fig. 1. Confirmed cases in Mexico between 14 March 2009 and 27 May 2009.

Figure 1

Table 1. Number of travellers and earliest date of cases imported from Mexico to a particular country in March and April, 2009

Figure 2

Fig. 2. Values of the minimum square root of the sum of standardized squared errors (RSE) and R0 given different constant r.

Figure 3

Table 2. Estimates of reporting rates (bootstrapped 95% confidence intervals) given different variations

Figure 4

Fig. 3. The effect of variations from the length of the latent period (∼gamma[mean = 1·6, s.d. = 0·5]) and the length of the infectious period (∼gamma[mean = 1·4, s.d. = 0·5]) given a constant reporting rate assumption. Left panel is the box-plot of r and the right panel is the box-plot of R0.

Figure 5

Fig. 4. The effect of variations in the lengths of the latent and infectious periods given a linearly increasing reporting rate assumption. The impacts to rmin, rmax, and R0 are shown by the box-plots from left, middle, and right panels, respectively.