Network Perspectives on Infectious Disease DynamicsView this Special Issue
Research Article | Open Access
Hiroshi Nishiura, Alex R. Cook, Benjamin J. Cowling, "Assortativity and the Probability of Epidemic Extinction: A Case Study of Pandemic Influenza A (H1N1-2009)", Interdisciplinary Perspectives on Infectious Diseases, vol. 2011, Article ID 194507, 9 pages, 2011. https://doi.org/10.1155/2011/194507
Assortativity and the Probability of Epidemic Extinction: A Case Study of Pandemic Influenza A (H1N1-2009)
Unlike local transmission of pandemic influenza A (H1N1-2009), which was frequently driven by school children, most cases identified in long-distance intranational and international travelers have been adults. The present study examines the relationship between the probability of temporary extinction and the age-dependent next-generation matrix, focusing on the impact of assortativity. Preferred mixing captures as a good approximation the assortativity of a heterogeneously mixing population. We show that the contribution of a nonmaintenance host (i.e., a host type which cannot sustain transmission on its own) to the risk of a major epidemic is greatly diminished as mixing patterns become more assortative, and in such a scenario, a higher proportion of non-maintenance hosts among index cases elevates the probability of extinction. Despite the presence of various other epidemiological factors that undoubtedly influenced the delay between first importations and the subsequent epidemic, these results suggest that the dominance of adults among imported cases represents one of the possible factors explaining the delays in geographic spread observed during the recent pandemic.
Since it was first identified in early 2009, a novel strain of influenza A (H1N1-2009) has caused a global pandemic. Although the rapid international spread created various epidemiological challenges, such as quantifying the strain's transmission potential and virulence during the very early stages of the pandemic [1, 2], many key insights have been obtained to date . Prior to the pandemic, the importance of contact networks in elucidating the epidemiological dynamics of infectious diseases has been emphasized with applications to severe acute respiratory syndrome (SARS), sexually transmitted infections and other directly transmitted diseases [4, 5]. The age specificity in the transmission of the H1N1-2009 indicates the relevance of contact heterogeneity [6–9]. Although the differential attack rates in different age groups by H1N1-2009 have multiple explanatory factors, including age-specific susceptibility and pre-existing immunity [10–12], age-dependent contact is also thought to be associated with a higher susceptibility to infection and greater infectiousness once infected in children . The consequences of this are that community-wide epidemics have been frequently driven by school outbreaks , while attack rates of H1N1-2009 were highest among school-age children in various parts of the world [8, 11]. A network model was used to describe the temporal variations in the age-specific composition of cases during the course of the pandemic and found that attack rates for a novel strain of influenza tend to be initially biased towards children and then shift towards adults .
A parsimonious simplification of the complexity of an age-structured contact network can be obtained by approximating the network by an appropriately quantified age-dependent next-generation matrix. This is accomplished by using the next-generation matrix, the square matrix with generic entry , the average number of secondary cases in age-group generated by a single primary case in age-group in a fully susceptible population. Discretizing chronological age into a small number of age groups, the matrix, , for H1N1-2009 has been quantified using age-stratified epidemic data [2, 6, 7, 15]. The age-dependent has two important properties in understanding epidemiological dynamics. First, the dominant eigenvalue of corresponds to the basic reproduction number, , which frequently yields a threshold condition, that is, a major epidemic is possible if and only if . Second, the final proportion of hosts of type that become infected, , is given by the solution to . One important use of is to the development of optimal vaccination strategies before a pandemic [18, 19] and during a pandemic .
The present study investigates the relationship between the age-dependent next-generation matrix, , and the invasion of a novel virus into a large population, which has not been well clarified to date. Whereas this subject has been partly explored via percolation theory , it is important from an epidemiological perspective to address this using a simpler model that can more readily be fitted to observational data from the outbreak in question. More specifically, we examine the impact of assortativity, that is, preferential mixing of different host types (here, age), on the probability of epidemic extinction, because age-dependent human contact networks have been shown to be highly assortative in contact surveys [21, 22]. Such an assortative network is known to allow disease percolation more easily than disassortative ones , echoed by the finding that increasing the preferential mixing component of simpler models such as ours tends to allow an epidemic to grow more easily .
In addition to these issues, the present study investigates the role of the age of cases importing infection to a local area by long-distance travel either intranationally or internationally on the resulting growth of a local epidemic. As a practical example, the age-dependent transmission of the H1N1-2009 pandemic is considered, and we first present our study motivations in the next section.
2. Materials and Methods
2.1. Study Motivation
Figure 1(a) shows the time delay from the introduction of first imported case on 1 May 2009 to the subsequent increase in local transmission in Hong Kong . Despite the number of imported cases, it took 39 days to observe the first locally acquired case. The interpretation of Figure 1(a) is affected by several factors, such as case ascertainment, ecological factors such as seasonality, and disease control , but as with many other countries exponential growth in the local epidemic did not start for some time after the first imported case. Hong Kong instigated particularly stringent quarantine measures, but a recent study comparing the time delay in local transmission between countries with and without entry screening has shown that the entry screening measures were not associated with a substantial delay in the start of local transmission .
Figure 1(b) shows weekly hospitalization rates due to H1N1-2009 in three coastal areas in the Netherlands . A surge in hospitalizations is first seen in Amsterdam followed by Rotterdam. The peak hospitalization rate in Zeeland occurs three weeks later than that in Amsterdam. There may be various interpretations for the delay before exponential growth, and, in particular, the spatial heterogeneity in Figure 1(b) is likely to have been associated with differing inflows of infected individuals and intrinsically differing patterns of spread within each region. Despite the presence of various possible factors explaining Figure 1(b), it is clear that the spatio-temporal dynamics are not synchronized even in this geographically limited country, and thus, Figure 1(b) at least indicates that stochastic effects may not have been insignificant for the intranational spread. A similar substantial delay in interregional spread has also been seen in the results of seroepidemiological study in England .
Both Figures 1(a) and 1(b) indicate a delay in causing international or interregional spread, but from a sufficiently high number of homogeneous index cases that repeated stochastic extinction is unlikely as an explanation. A more plausible reason is the contrast in age distributions between local and imported cases: whereas imported cases have been predominantly adults , local transmissions are frequently driven by school children. That is, adults were more likely to travel than children, and those aged 25 years and older accounted for more than half of the imported cases in Japan . Similarly, adults may also more likely be the source of spread within a country, especially as the movement distance becomes longer. However, adults are less likely to cause secondary transmissions than children in a local setting [2, 6, 13, 15], making it critically important to understand the differential probability of extinction of the infection tree emerging from a typical child index and a typical adult index case. Because assortativity regulates the frequencies of within- and between-group transmissions, examining the effects of assortativity provides a natural avenue for assessing this. Accordingly, in this paper we use a simple stochastic model to clarify the different roles of children and adults in causing a major epidemic and its relevance to assortativity.
2.2. A Model for Clade Extinction
We employ a multitype branching process to approximate the probability of extinction of the clade of infection emanating from a single index case . Consider a large population which is fully susceptible, and let there be two subpopulations, that is, children and adults. For simplicity, we ignore pre-existing immunity among adults. Throughout this paper, we label children as type 1 and adults as type 2. Let () be the recovery rate of infectious individuals of type and be the birth rate (i.e., the rate of new infection) of type infected individuals caused by a single type infected individual during the initial stage of an epidemic. We consider the case when a small number of infected individuals of type invades a fully susceptible large population. Given the large and (assumed) fully susceptible population, and the small initial number of infectives, depletion of the susceptible stock can be ignored and the initial stages of the outbreak interpreted as a multivariate birth-and-death process . For mathematical convenience, we assume that the generation time is exponentially distributed, and thus, . This assumption is common to many compartmental models, although its realism is dubious. The proposed approach considers a linearized system for the early epidemic period with a crude approximation of host types, but the similar approach of mapping next generation with the use of a square matrix can be employed for explicit network models . In addition to the aforementioned assumptions, we assume that the age-specificity of is fully attributable to the infection rate , and thus the infectious period is assumed to be a constant , independent of host type. Consequently, if we further ignore the age-specific susceptibility and infectiousness, is determined only by the frequency of contact within- and between-age groups.
Letting the random vector represent the number of child and adult infected individuals in the population in the th generation, we consider the process as a multitype branching process. Assuming that an individual of type has probability of infecting, in the next generation, children and adults, we define the probability generating function as Following foregoing studies [32, 34], the generating function with an exponentially distributed generation time is known to be given by for . Since is assumed to be independent of host type , , (2) simplifies to The clade of infections, , emanating from the initial index cases becomes extinct with probability 1 if and only if the dominant eigenvalue of is less than or equal to unity, that is, .
Let be the probability of extinction given that a single infected individual of type is introduced to the population. The extinction probability is the nonnegative root of the equations As is standard in branching process models, each of the secondary cases of type generated by a primary case becomes an ancestor of an independent subprocess (which restarts with a type individual) behaving identically among the same type [31, 35]. Because of this multiplicative nature, we have the probability of extinction of the entire clade with initial vector .
In the two-host population, that is, a population consisting of children and adults, the probabilities of extinction given a single child or adult infected individual, and , satisfy In other words, given that the next-generation matrix is known, the problem of calculating the probability of extinction given a certain number of infected individuals of host and in the zero generation is replaced by the problem of solving two quadratic equations with two unknown parameters. There are four possible combinations of the solutions for (6) including complex numbers, but we iteratively find the only nonnegative real numbers in the range of (see [31, page 18]), except for a combination .
2.3. Quantitative Illustrations
The probability of extinction is investigated for the following three different scenarios. First, to gain an overview of the extinction probabilities and for the H1N1-2009, (6) are solved using published estimates of from Mexico  and Japan . Approximating the original into a two-host population, we use for Mexico, and for Japan. The originally estimated dominant eigenvalues are 1.58 and 1.22, respectively. It should be noted that the child group in Mexico is assumed to be up to age 14 years while that in Japan is up to age 19 years. Assuming that the reproduction number, possibly ranges from 1.2–1.6 [2, 36], we rescale the next-generation matrices by where is the reproduction number to be examined.
Second, in Mexico is further examined in relation to the assortativity. The element of has been parameterized as where is the relative size of the subpopulation (i.e., ). and are originally described as the age-specific susceptibility and infectiousness , and these can also be regarded as the so-called proportionate mixing components. The biological interpretation of proportionate mixing is that irrespective of its own type, an individual can acquire infection from any given individual (i.e., the secondary transmission from host to is determined by host ). Introduction of the most important parameter in the present study, is classically referred to as “preferred” or “preferential” mixing [37, 38]. Although the original definition of the term preferred mixing has a broader meaning, in (10) represents the proportion of contacts reserved for within-group mixing, and () represents the proportion of contacts subject to proportionate mixing. If , the mixing is referred to as fully assortative (Figure 2). If , the mixing corresponds to random mixing (though it should be noted that the mixing matrix still includes a proportionate mixing component). An empirical estimate of from Mexico is 0.50, although the 95% confidence interval is broad: 0–0.72 . Therefore, we examine the sensitivity of the probabilities of extinction, and , to different in the range of 0-1 and in the range of 1.2–1.6. Other parameters are fixed at , , .
Third, to clarify the practical implications of the predomination of adults among travelers, we examine the sensitivity of the probability of extinction to the proportion of adult travelers over various and . Specifically, we consider the probability of extinction given a small importation of ten cases in the zero generation independently entering a large susceptible population at their infection-age 0 (i.e., immediately after their own infections: for simplicity, we ignore the infection-age distribution of imported cases at the time of invasion in the present study, because its realistic incorporation enforces us to account for the epidemic dynamics in exporting countries and thus, the exporting country and travel distance for each imported case would be required ). Among the 10 cases, we vary the number of adult cases from 0 to 10, and examine the probability of extinction given by (5).
Figure 3 shows the probabilities of extinction, and , using published estimates of and in (7) and (8). In both panels, using estimates from Mexico and Japan, , the probability of extinction given a single adult case, always appeared to be higher than , and thus the clade of infections resulting from the introduction of an adult index case is more likely to be self-limiting than from a child index case. The estimates of and using the published estimates of were 61 and 73%, respectively, for Mexico (with ) and 81 and 93%, respectively, for Japan (with ), indicating that the reproduction number in the range of 1.2–1.6 is not far from the critical level and the impact of variations in on epidemic extinction is large. The reader should note the crudeness of the dichotomization of the population into two subpopulations, and that incorporating more detailed network structure (e.g., by dividing the population into many more types of host) tends to yield higher probability of extinction [2, 5]. Moreover, whereas the present study assumes an exponentially distributed generation time, a more realistic depiction, for example, gamma-distributed generation time, tends to capture overdispersion of the offspring distribution more appropriately [40, 41], and thus, again yields a higher probability of extinction than is shown herein.
Figure 4 examines the probabilities of extinction, and , as a function of , the proportion of within-group mixing and . As expected from the randomly mixing interpretation, and were equal to with . However, for populations with more within-group mixing, clades from adults were more likely to go extinct, reaching 100% with . This is attributable to the next-generation matrices (7) and (8) involving the typical reservoir dynamics : children act as a maintenance host (), among whom transmission can be maintained by themselves, while adults constitute a nonmaintenance host group (), and thus, with little relative mixing between the two groups, an adult index case would never lead to a major epidemic. The probability of extinction, , given a single child index case, reached a minimum with in the range of 0.4-0.5, although this probability was not very sensitive to . Such may lead to a “well-mixed” population, thereby allowing the child index case to involve both child and adult secondary cases effectively.
Figure 5 examines the probability of extinction given 10 index cases as a function of and the proportion of adult index cases. With random mixing, the probability of extinction given a single index case is ()% (e.g., 71% with ). The extinction probability given 10 index cases in the randomly mixing population was independent of the proportion of adults, (with ), indicating that a major epidemic is almost unavoidable without any intervention. However, as assortativity increased, the increase in the proportion of adult index cases promoted extinction. The vertical reference line of 63% indicates the empirically observed proportion of adults (i.e., those aged ≥15 years in (7) as defined by ) among all imported cases in Japan . At around that proportion, the probability of extinction was estimated to be 1%–58% over the full range of from 0 to 1, for 10 index cases. Given 10 index cases with , the results were independent of the proportion of adult index cases, because and (for ), whatever the number of adults, , the extinction probability is That is, as mixing becomes more and more assortative, the contribution of the initial number of adults, to the risk of a major epidemic becomes less important, and moreover, an increase in the proportion of adults indirectly reduces , leading to an increase in the probability of extinction. Even when we divide the entire population into many more subpopulations, this argument holds as long as the host type of interest is incapable of maintaining disease by itself, that is, when the host-specific reproduction number, . Of course, if , for .
The present study investigated the relationship between the next-generation matrix and the probability of extinction, employing a simple model that may be viewed as an approximation to a full network model. The modelled heterogeneous mixing accounted for assortativity via an assumption of preferred mixing, and the probability of extinction was derived from a multidimentional branching process model. As a practical example, the age dependency in the transmission of pandemic influenza A (H1N1-2009) was considered, dividing the population into children and adults. Through quantitative illustrations, it has been shown that the probability of extinction given an adult index case increases with , at least for diseases with similar transmissibility as influenza. Although this exercise employed several simplifying assumptions, a formal hypothesis can be developed for explaining a slow interregional and international spread of the H1N1-2009 even in today's highly mobile world population. That is, whereas empirically observed delays in local transmission can be influenced by a large number of factors including pre-existing immunity, public health interventions and seasonality, the dominance of adults among travelers is one possible explanation for the high probability of extinction, and may play an important role in describing the underlying reason (Figure 1). Since the present study adopted three simplifying assumptions (i.e., (1) the crude dichotomization of hosts into two different types, (2) the adoption of exponentially distributed generation times, and (3) ignorance of infection-age among imported cases), the probability of extinction is likely to have been underestimated. The extinction probabilities become higher with more precise network structure (e.g., due to localized burnout of susceptible individuals) and more detailed natural history of infection [5, 40, 41].
Three practical implications are drawn from our exercise. First, the importance of assortativity in appropriately capturing the probability of extinction highlights a critical need to account for this aspect when quantifying the next-generation matrix in an approximately modelled heterogeneous population. Whenever the statistical inference of the next-generation matrix is made for directly transmitted diseases, the estimation framework should ideally account for assortative mixing. Whereas the social contact survey revealed that the age-dependent contact pattern is highly assortative [20, 21], the definition of a contact can be too broad to be practical for all diseases, and more realistic incorporation of assortative mixing and its precise estimation should be the subject of future studies.
Second, as was highlighted with an application, accounting for the age specificity in the surveillance of international and interregional mobility patterns and its use for statistical inference of epidemic dynamics are of utmost importance. For example, global airline transportation is one of the most well-studied networks, and this has been analyzed for H1N1-2009 , but a full description of global dynamics should better account for age-specific travel patterns. In addition, whereas imported cases from Mexico have been utilized to make statistical inference (e.g., spatial backcalculation) of the incidence in Mexico [6, 44, 45], the present study emphasizes a critical need to examine age-specificity in relevant frameworks, so that ultimately, the global dynamics can be described by a multihost metapopulation model [46, 47].
Third, as a disease control implication, although adults dominate imported cases, it should be remembered that the more important target host is still children. If stringent border control measures, for example, travel reduction and movement restrictions among all incoming passengers , are adopted as containment strategies against a highly virulent novel virus, the target host to promote radical reductions in travel-induced illness would be children, at least for diseases with a similar next-generation matrix to that of the recent pandemic.
Although the role of heterogeneously mixing population in the spread of infectious diseases has been examined using stochastic modelling approaches, past studies tended to focus on final epidemic size and its relevance to disease control policy [49, 50]. Moreover, a limited number of studies examining the probability of extinction took an average of the probabilities over different types of host (e.g., by weighting the relative population size to the type specific probability of extinction) [24, 49]. The present study emphasized the importance of capturing type-specificity of index cases in estimating the probability of extinction and examining the impact of assortativity on extinction. In conclusion, we believe our simple exercise successfully illustrated the diminished role of nonmaintenance hosts in causing a major epidemic when assortativity is high, indicating a critical need to capture the assortativity in modelling the initial invasion of an epidemic disease.
Unlike local transmission of the H1N1-2009 which was frequently driven by school children, imported cases were predominantly adults. This study examined the relationship between the age-dependent next-generation matrix and the probability of extinction, focusing on the role of nonmaintenance hosts and the impact of assortativity on the epidemic extinction. The preferred mixing assumption captures assortativity in a much simpler way than full contact network models, allowing analysis in place of Monte Carlo calculations. The contribution of nonmaintenance hosts to the risk of a major epidemic is diminished as the mixing pattern becomes more assortative, so that an increase in the proportion of nonmaintenance hosts among index cases increases the probability of extinction, if temporary in the face of repeat importations. These results helped us to formulate a hypothesis that the dominance of adults in imported cases was one of the possible causes of observing substantial delay in interregional and international spread of the 2009 influenza pandemic. The importance of capturing the assortativity in estimating the next-generation matrix was highlighted.
The work of H. Nishiura is supported by the Japan Science and Technology Agency (JST) PRESTO program. A. R. Cook is grateful to the National University of Singapore for supporting his research. B. J. Cowling is supported by the Harvard Center for Communicable Disease Dynamics from the US National Institutes of Health Models of Infectious Disease Agent Study program (grant no. 1 U54 GM088558).
- Anonymous, “Mathematical modelling of the pandemic H1N1 2009,” Weekly Epidemiological Record, vol. 84, no. 34, pp. 341–348, 2009.
- H. Nishiura, G. Chowell, M. Safan, and C. Castillo-Chavez, “Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009,” Theoretical Biology and Medical Modelling, vol. 7, no. 1, article 1, 2010.
- D. Butler, “Portrait of a year-old pandemic,” Nature, vol. 464, no. 7292, pp. 1112–1113, 2010.
- L. A. Meyers, B. Pourbohloul, M. E. J. Newman, D. M. Skowronski, and R. C. Brunham, “Network theory and SARS: predicting outbreak diversity,” Journal of Theoretical Biology, vol. 232, no. 1, pp. 71–81, 2005.
- M. E. J. Newman, “Spread of epidemic disease on networks,” Physical Review E, vol. 66, no. 1, Article ID 016128, pp. 1–11, 2002.
- C. Fraser, C. A. Donnelly, S. Cauchemez et al., “Pandemic potential of a strain of influenza A (H1N1): early findings,” Science, vol. 324, no. 5934, pp. 1557–1561, 2009.
- H. Nishiura, C. Castillo-Chavez, M. Safan, and G. Chowell, “Transmission potential of the new influenza A(H1N1) virus and its age-specificity in Japan,” Eurosurveillance, vol. 14, no. 22, p. 19227, 2009.
- D. N. Fisman, R. Savage, J. Gubbay et al., “Older age and a reduced likelihood of 2009 H1N1 virus infection,” New England Journal of Medicine, vol. 361, no. 20, pp. 2000–2001, 2009.
- A. L. Greer, A. Tuite, and D. N. Fisman, “Age, influenza pandemics and disease dynamics,” Epidemiology and Infection, vol. 138, no. 11, pp. 1542–1549, 2010.
- Y. Itoh, K. Shinya, M. Kiso et al., “In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses,” Nature, vol. 460, no. 7258, pp. 1021–1025, 2009.
- T. Reichert, G. Chowell, H. Nishiura, R. A. Christensen, and J. A. McCullers, “Does Glycosylation as a modifier of Original Antigenic Sin explain the case age distribution and unusual toxicity in pandemic novel H1N1 influenza?” BMC Infectious Diseases, vol. 10, article 5, 2010.
- M. I. C. Chen, V. J. M. Lee, W. Y. Lim et al., “2009 Influenza A(H1N1) seroconversion rates and risk factors among distinct adult cohorts in Singapore,” Journal of the American Medical Association, vol. 303, no. 14, pp. 1383–1391, 2010.
- B. J. Cowling, K. H. Chan, V. J. Fang et al., “Comparative epidemiology of pandemic and seasonal influenza A in households,” New England Journal of Medicine, vol. 362, no. 23, pp. 2175–2184, 2010.
- S. Bansal, B. Pourbohloul, N. Hupert, B. Grenfell, and L. A. Meyers, “The shifting demographic landscape of pandemic influenza,” PLoS ONE, vol. 5, no. 2, Article ID e9360, 2010.
- Y. Yang, J. D. Sugimoto, M. Elizabeth Halloran et al., “The transmissibility and control of pandemic influenza a (H1N1) virus,” Science, vol. 326, no. 5953, pp. 729–733, 2009.
- O. Diekmann, J. A. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.
- F. Ball and D. Clancy, “The final size and severity of a generalised stochastic multitype epidemic model,” Advances in Applied Probability, vol. 25, no. 4, pp. 721–736, 1993.
- I. M. Longini, E. Ackerman, and L. R. Elveback, “An optimization model for influenza A epidemics,” Mathematical Biosciences, vol. 38, no. 1-2, pp. 141–157, 1978.
- J. Medlock and A. P. Galvani, “Optimizing influenza vaccine distribution,” Science, vol. 325, no. 5948, pp. 1705–1708, 2009.
- J. Wallinga, M. Van Boven, and M. Lipsitch, “Optimizing infectious disease interventions during an emerging epidemic,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 2, pp. 923–928, 2010.
- J. Mossong, N. Hens, M. Jit et al., “Social contacts and mixing patterns relevant to the spread of infectious diseases,” PLoS Medicine, vol. 5, no. 3, article e74, pp. 0381–0391, 2008.
- J. Wallinga, P. Teunis, and M. Kretzschmar, “Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents,” American Journal of Epidemiology, vol. 164, no. 10, pp. 936–944, 2006.
- M. E. J. Newman, “Assortative mixing in networks,” Physical Review Letters, vol. 89, no. 20, Article ID 208701, pp. 1–4, 2002.
- I. C. Marschner, “The effect of preferential mixing on the growth of an epidemic,” Mathematical Biosciences, vol. 109, no. 1, pp. 39–67, 1992.
- The Government of the Hong Kong Special Administrative Region, “Update on human swine influenza (July 2009),” The Government of the Hong Kong Special Administrative Region, Hong Kong, 2009.
- Rijksinstituut voor volksgezondheid en milieu (RIVM), “The Netherlands, Ziekenhuisopnamen Nieuwe Influenza A (H1N1),” RIVM, Bilthoven, 2010.
- J. T. Wu, B. J. Cowling, E. H. Y. Lau et al., “School closure and mitigation of pandemic (H1N1) 2009, Hong Kong,” Emerging Infectious Diseases, vol. 16, no. 3, pp. 538–541, 2010.
- B. J. Cowling, L. L. H. Lau, P. Wu et al., “Entry screening to delay local transmission of 2009 pandemic influenza A (H1N1),” BMC Infectious Diseases, vol. 10, article 82, 2010.
- E. Miller, K. Hoschler, P. Hardelid, E. Stanford, N. Andrews, and M. Zambon, “Incidence of 2009 pandemic influenza A H1N1 infection in England: a cross-sectional serological study,” The Lancet, vol. 375, no. 9720, pp. 1100–1108, 2010.
- H. Nishiura, “Travel and age of influenza A (H1N1) 2009 virus infection,” Journal of Travel Medicine, vol. 17, no. 4, pp. 269–270, 2010.
- C. J. Mode, “Multitype branching processes: theory and applications,” in Modern Analytic and Computational Methods in Science and Mathematics, R. Bellman, Ed., American Elsevier, New York, NY, USA, 1971.
- D. A. Griffiths, “Multivariate birth-and-death processes as approximations to epidemic processes,” Journal of Applied Probability, vol. 10, no. 1, pp. 15–26, 1973.
- M. Boguñá and R. Pastor-Satorras, “Epidemic spreading in correlated complex networks,” Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, vol. 66, no. 4, Article ID 047104, pp. 1–4, 2002.
- F. Ball, “The threshold behaviour of epidemic models,” Journal of Applied Probability, vol. 20, no. 2, pp. 227–241, 1983.
- T. E. Harris, The Theory of Branching Processes, Springer, Berlin, Germany, 1963.
- Anonymous, “Transmission dynamics and impact of pandemic influenza A (H1N1) 2009 virus,” Weekly Epidemiological Record, vol. 84, no. 46, pp. 481–484, 2009.
- J. A. Jacquez, C. P. Simon, J. Koopman, L. Sattenspiel, and T. Perry, “Modeling and analyzing HIV transmission: the effect of contact patterns,” Mathematical Biosciences, vol. 92, no. 2, pp. 119–199, 1988.
- H. W. Hethcote, “Modelling heterogeneous mixing in infectious disease dynamics,” in Models for Infectious Human Diseases: Their Structure and Relation to Data, V. S. Isham and G. F. Medley, Eds., pp. 215–238, Cambridge University Press, Cambridge, UK, 1996.
- P. Mukherjee, P. L. Lim, A. Chow et al., “Epidemiology of travel-associated pandemic (H1N1) 2009 infection in 116 patients, Singapore,” Emerging Infectious Diseases, vol. 16, no. 1, pp. 21–26, 2010.
- T. Garske and C. J. Rhodes, “The effect of superspreading on epidemic outbreak size distributions,” Journal of Theoretical Biology, vol. 253, no. 2, pp. 228–237, 2008.
- J. O. Lloyd-Smith, S. J. Schreiber, P. E. Kopp, and W. M. Getz, “Superspreading and the effect of individual variation on disease emergence,” Nature, vol. 438, no. 7066, pp. 355–359, 2005.
- H. Nishiura, B. Hoye, M. Klaassen, S. Bauer, and H. Heesterbeek, “How to find natural reservoir hosts from endemic prevalence in a multi-host population: a case study of influenza in waterfowl,” Epidemics, vol. 1, no. 2, pp. 118–128, 2009.
- K. Khan, J. Arino, W. Hu et al., “Spread of a novel influenza A (H1N1) virus via global airline transportation,” New England Journal of Medicine, vol. 361, no. 2, pp. 212–214, 2009.
- M. Lipsitch, M. Lajous, J. J. O'Hagan 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, vol. 4, no. 9, Article ID e6895, 2009.
- D. Balcan, H. Hu, B. Goncalves et al., “Seasonal transmission potential and activity peaks of the new influenza A(H1N1): a Monte Carlo likelihood analysis based on human mobility,” BMC Medicine, vol. 7, article 45, 2009.
- J. Arino, J. R. Davis, D. Hartley, R. Jordan, J. M. Miller, and P. van den Driessche, “A multi-species epidemic model with spatial dynamics,” Mathematical Medicine and Biology, vol. 22, no. 2, pp. 129–142, 2005.
- J. Arino, R. Jordan, and P. van den Driessche, “Quarantine in a multi-species epidemic model with spatial dynamics,” Mathematical Biosciences, vol. 206, no. 1, pp. 46–60, 2007.
- H. Nishiura, N. Wilson, and M. G. Baker, “Quarantine for pandemic influenza control at the borders of small island nations,” BMC Infectious Diseases, vol. 9, article 27, 2009.
- N. Becker and I. Marschner, “The effect of heterogeneity on the spread of disease,” in Stochastic Processes in Epidemic Theory, J. P. Gabriel, C. Lefevre, and P. Picard, Eds., vol. 86 of Lecture Notes in Biomathematics, pp. 90–103, Springer, New York, NY, USA, 1990.
- H. Andersson and T. Britton, “Heterogeneity in epidemic models and its effect on the spread of infection,” Journal of Applied Probability, vol. 35, no. 3, pp. 651–661, 1998.
Copyright © 2011 Hiroshi Nishiura et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.