Stochastic modelling of infectious diseases for heterogeneous populations
 RuiXing Ming^{1},
 Jiming Liu^{2},
 William K. W. Cheung^{2} and
 Xiang Wan^{3}Email author
DOI: 10.1186/s4024901601995
© The Author(s) 2016
Received: 21 December 2015
Accepted: 10 October 2016
Published: 22 December 2016
Abstract
Background
Infectious diseases such as SARS and H1N1 can significantly impact people’s lives and cause severe social and economic damages. Recent outbreaks have stressed the urgency of effective research on the dynamics of infectious disease spread. However, it is difficult to predict when and where outbreaks may emerge and how infectious diseases spread because many factors affect their transmission, and some of them may be unknown.
Methods
One feasible means to promptly detect an outbreak and track the progress of disease spread is to implement surveillance systems in regional or national health and medical centres. The accumulated surveillance data, including temporal, spatial, clinical, and demographic information can provide valuable information that can be exploited to better understand and model the dynamics of infectious disease spread. The aim of this work is to develop and empirically evaluate a stochastic model that allows the investigation of transmission patterns of infectious diseases in heterogeneous populations.
Results
We test the proposed model on simulation data and apply it to the surveillance data from the 2009 H1N1 pandemic in Hong Kong. In the simulation experiment, our model achieves high accuracy in parameter estimation (less than 10.0 % mean absolute percentage error). In terms of the forward prediction of case incidence, the mean absolute percentage errors are 17.3 % for the simulation experiment and 20.0 % for the experiment on the real surveillance data.
Conclusion
We propose a stochastic model to study the dynamics of infectious disease spread in heterogeneous populations from temporalspatial surveillance data. The proposed model is evaluated using both simulated data and the real data from the 2009 H1N1 epidemic in Hong Kong and achieves acceptable prediction accuracy. We believe that our model can provide valuable insights for public health authorities to predict the effect of disease spread and analyse its underlying factors and to guide new control efforts.
Keywords
Epidemiology Stochastic model Surveillance system Spread patternMultilingual abstracts
Please see Additional file 1 for translations of the abstract into the five official working languages of the United Nations.
Background
Infectious diseases remain a major cause of morbidity and mortality worldwide, triggering immeasurable loss in many societies. Most people may still have a fresh memory of the H1N1 outbreak in 2009, which brought pictures of empty streets and people wearing face masks and collectively caused at least 12799 deaths according to the World Health Organization (WHO) report [1]. The H1N1 pandemic calls for research on accurately modelling the spread dynamics of an infectious disease, which offers a practically useful means for policy makers to evaluate the potential effects of intervention strategies [2–4].
Here, β≥0 is the effective transmission rate and k≥0 is the recovery rate. Because the SIRbased models are well presented in the literature, herein, we omit a verbose introduction of these models. Readers with an interest in such a topic can find the details in [5–7].
The SIRbased models and its variants have proven to be quite useful in the study of the spread dynamics of infectious diseases [8–10]. In [11–13], the progression of disease spread is characterized by tracking the number of S _{ t } with a chain binomial model. The number of susceptible members S _{ t+△t } (△t represents the infectious period of the disease and is always chosen to be 1/k) at time t+△t is a binomial random variable that depends on S _{ t } and I _{ t } α, S _{ t+△t }∼B i n(S _{ t },1−I _{ t } α), which provides a recursive relationship between S _{ t+△t } and S _{ t } and produces a formal stochastic process. However, the power of these models is mainly limited to uniform and homogeneous populations or populations with infinite size and homogeneous interactions. In many cases, the actual spread of infectious diseases occurs in a diverse or dispersed population. To study the spread of infectious diseases in heterogeneous populations, people usually divide a population into subpopulations that differ from each other. Subpopulations can be determined on the basis of social, cultural, economic, demographic, and geographic factors. Next, besides the dynamics of the internal spread within a subpopulation, the transmission dynamics between subpopulations should also be considered in the study of epidemic spreading.
Networkbased epidemic modelling represents a popular approach for heterogeneous populations in which the nodes in the network correspond to subpopulations, and the links indicate the neighboring relationships. Many networkbased models have been proposed, including patch models [14–16], distancetransmission models [17], and multigroup models [18, 19]. However, these models require knowledge of every individual (or host) and all relationships between individuals, which may be not achievable due to information privacyrelated restrictions and the high cost of subject recruitment. To overcome the difficulties of collecting data, researchers have investigated several types of computergenerated networks in the context of disease spread in populationscale studies [20–24]. Grassberger first studied the dynamics of infectious diseases that propagate on regular networks using the percolation theory [25]. Recent studies have revealed that many realworld networks, including social networks in which infectious diseases propagate, are either smallworld [26] or scalefree [27] rather than regular or random, as thought previously [28]. Because the underlying structures of networks will influence the effect that the dynamics of epidemics will have on them, researchers, such as PastorSatorras and Vespignani, have made many contributions to critical value analysis of typical epidemics on different types of complex network [23, 24, 29]. On the basis of the meanfield theory, they found that compared with homogeneous networks, scalefree networks are fragile to the invasion of infectious diseases, computer viruses, or any other type of negative epidemics.
Epidemics have also been studied in various disciplines. Sociologists are concerned with the diffusion of rumors or innovation on social networks [30]; economists have studied viral marketing and recommendation strategies by considering both cascading dynamics and the network effects of vital nodes [31]; and computer scientists are interested in how some topics can quickly cascade in virtual blog spaces and how their propagation trends [32, 33].
Although networkbased studies have contributed to the modelling of disease and/or information dynamics, some models make a strong assumption that the structures of underlying networks over which epidemics spread are known beforehand. In the real world, however, the structures of underlying diffusion networks are not known directly. Many others assume the availability of information about the interactions occurring between individuals [34–37] that are often not valid in the context of disease spread. What may be obtained is only the time at which particular subpopulations become infected, but not how they become infected, nor how they affect their neighboring areas. Moreover, the underlying structures of networks will greatly influence the dynamics of infectious disease spread.
Since the emergence of the H1N1 influenza pandemic in April 2009, its underlying dynamics have been of great public health interest, and many approaches for its study have been proposed [14, 38–41]. Most of them are based on the classic SIR model. For example, Birrell et al. [40] provided an age structurebased compartmental model with a Bayesian synthesis of multiple evidence sources to reveal substantial changes in contact patterns throughout the epidemic. Besides of the compartmental models, other mathematical models are also used to describe the transmission dynamics [3, 42–47]. The chain binomial model was used to calculate the household secondary attack rates to measure the transmissibility of the 2009 H1N1 influenza pandemic by Lessler et al. [44] and Klick et al. [45]. Yang et al. [46] constructed a model based on chains of infections and used the infection hazard function and survival function to study the 2009 H1N1 influenza pandemic. Ferguson et al. [3] and Cauchemez et al. [42, 43] incorporated other factors, such as household risk, withinschool risk, and community risk, in the study of infection spread and found out that younger age groups under 19 years old were more susceptible than older age groups. Jin et al. [47] formulated an epidemic model of influenza A based on networks and calculated the basic reproduction number and studied the effects of various immunization schemes. However, this work required that the individual contact pattern be provided. Nonetheless, none of the aforementioned approaches takes spatial heterogeneity into consideration in the study of disease spread.
Recently, an outbreak of Ebola virus disease (EVD) swept across parts of West Africa from March 2014 to April 2015. By June 10, 2015, WHO had reported 27,237 confirmed, probable, or suspected cases in three countries with 11,158 deaths [48]. This epidemic received extensive research attention on its dynamics of spread [49–57] (for further references in the review article [58]). To name a few, Chowell et al. found that districtlevel Ebola virus disease outbreaks in West Africa follow polynomialbased growth in time instead of the exponential growth that describes the progress of many infectious disease epidemics [52]. Fisman et al. used a simple, two parameter mathematical model to characterize epidemic growth patterns in the 2014 Ebola outbreak [53]. Webb et al. proposed a variant of the classic SIR model with three extra groups, incubating, contaminated and isolated, which can provide a more accurate prediction for the future incidences [56]. Carroll et al. used a deep sequencing approach to gain insight into the evolution of the Ebola virus (EBOV) in Guinea from the ongoing West African outbreak. The viral sequence data can be combined with epidemiological information to retrospectively test the effectiveness of control measures, and provides an unprecedented window into the evolution of an ongoing outbreak of viral haemorrhagic fever [57].
To accurately predict when and where outbreaks will occur, a feasible means is to deploy manual or electronic surveillance systems through regional or national public health and medical organizations [59]. Most of the surveillance data accumulated from such systems contains temporal, spatial, clinical, and demographic information. For instance, Telehealth Ontario is a teletriage helpline that is available free to all Ontario residents, which allows those with suspected infections to connect with experts who can assess their symptoms. The records of such calls provide valuable information on which individual from where was possibly infected and by which type of disease at what time. In this paper, we address the problem of modelling disease spread dynamics in heterogeneous populations from temporalspatial surveillance data. We analyse the role of heterogeneity in a stochastic epidemic model on a twodimensional lattice. Within a particular subpopulation, the speed of spread is controlled by a single parameter, the transmissibility of the pathogen between individuals. Between subpopulations, the transmissibility becomes a random variable drawn from a probability distribution. Our work differs from existing studies in some fundamental ways, in light of the unique nature of infectious disease diffusion dynamics. Our results have practical implications for the analysis of disease control strategies in realistic heterogeneous epidemic systems.
Methods
In this work, we propose a stochastic model to study the dynamics of infectious disease spread in heterogeneous populations from temporalspatial surveillance data. We divide the whole population into m subpopulations on the basis of geographic regions. In the following, we use S _{ i }(t),I _{ i }(t), and R _{ i }(t) to denote the number of susceptible, infected, and recovered people, respectively, at time t in region i, i=1,2,⋯,m and t∈[0,T].
Stochastic model
where α is a parameter that measures the autorecovery rate of one particular infectious disease, which is usually considered as a constant among subpopulations, δ _{ i } is the parameter that measures the different disease transmissibility in different subpopulations, σ _{ i }>0 is the diffusion parameter that measures the disease spread from neighbors, and B _{ i }(t) is a standard Brownian motion. It is worth noting that we assume the parameter δ _{ i }≠0 for technical purposes, and the results in the case of δ _{ i }=0 can be achieved with δ _{ i }→0.
Comparing our model in Eq. (2) with the classic model in Eq. (1), we can see that they both capture the situation in which the change in the number of infected people has a positive relationship with the total number of infected people, which means that the more infected people there are, the more people will get infected. There are two key differences between these two models: first, the key factor (β S _{ i }(t)−k) associated with the disease spread in Eq. (1) is replaced with a single parameter δ _{ i } in Eq. (2), which can be used to analyse the role of heterogeneity in the disease spread; and second, Eq. (2) takes the neighboring relationships into consideration to study the dynamics of the disease spread among different subpopulations.

α>−I _{ i }(0)δ _{ i }
In this case, E[I _{ i }(t)] tends to infinity as t goes to infinity, which implies that all people in that region will be infected if the time is long enough.

α=−I _{ i }(0)δ _{ i }
In this case, the pandemic or epidemic will reach a state of equilibrium.

α<−I _{ i }(0)δ _{ i }
In this case, E[I _{ i }(t)] will reach 0 at some time \(t = \hat {t}\) and go to negative infinity as t goes to infinity, which implies the pandemic or epidemic will end at time \(\hat {t}\).
Parameter estimation
\(\widehat {\alpha }_{i}\) is defined in Eq. (24).
Results and discussion
In this section, we illustrate the performance of our proposed model using both simulated and real data.
Simulation study

Set m=4 (the number of subpopulations) and T=100 (the number of time slots). These two numbers are randomly selected.

Randomly draw α from [0.05,0.09], δ _{ i } from [0.02,0.08], and σ _{ i } from [0.02,0.08].

Initialize I _{ i }(0),1≤i≤m.

Simulate I _{ i }(k+1)=I _{ i }(k)+△I _{ i }(k),k=0,1,⋯,T−1 using Eq. (7) and \(\triangle I_{i}(k)I_{i}(t_{k})\sim N\left ((\alpha +\delta _{i}I_{i}(t_{k}))\triangle t(k),{\sigma _{i}^{2}}\triangle t(k)\right)\).
Real application
Prediction error for real data
District  WAN CHAI  KWUN TONG  KOWLOON CITY  SAI KUNG  WONG TAI SIN  SHA TIN 

X Error  0.12  0.02  0.15  0.31  0.13  0.14 
District  TUEN MUN  CENTRAL & WEST  NORTH  KWAI TSING  EASTER  TSUEN WAN 
Error  0.09  0.09  0.08  0.30  0.25  0.65 
District  SOUTHERN  SHAM SHUI PO  TAI PO  YAU TSIM MONG  ISLANDS  YUEN LONG 
Error  0.38  0.06  0.15  0.29  0.08  0.25 
Conclusions
Epidemic modelling offers a practical means for policy makers to evaluate the potential effects of intervention strategies. To do so, the accuracy of epidemic modelling with respect to the realworld disease transmission dynamics is essential and remains a challenging task due to the inaccessibility of many factors that affect the spread patterns of infectious diseases. In particular, heterogeneity should be taken into consideration when modelling the disease spread in nonrandom mixing populations. Many methods have been proposed to deal with heterogeneity in the study of epidemic dynamics, mostly using networkbased epidemic models in which nodes correspond to spatial locations with reported incidences over time, and the directional links indicate the probability of disease transmission from one node to another over time. However, it is very challenging to determine the network topology. Many studies have used a geographical topology whereas others have used a mobility network inferred from the public transportation network or other sources. How to verify the inferred network topology is another challenging issue because the true epidemic network topology is unknown, and it may vary for different types of infectious diseases for the same population. Furthermore, the neighborhood effect estimation is nontrivial; it involves many parameters (a polynomial of the number of nodes) and requires a large amount of data to avoid overfitting. Such data may not always be available for the inference of network topology. Therefore, in this work, we propose an alternate approach to investigate the spatial heterogeneity from temporalspatial surveillance data without the inference of network topology.
Our proposed model possesses several merits over the previous works. First, it quantifies the role of the heterogeneity in the analysis of the spread dynamics of infectious diseases in heterogeneous populations. Second, parameter estimation can be computed very quickly. Therefore, the prediction and the corresponding intervention policies can be implemented without delay in an outbreak of infectious disease. We apply our model on both the simulated data and the real data from the 2009 H1N1 epidemic in Hong Kong and achieve acceptable prediction accuracy. Based on the study of disease diffusion, the model proposed in this work can be extended to study other propagation patterns such as the Internet and World Wide Web, through which individuals form multiple communities in which information can propagate in a manner similar to that of infectious disease. We believe that our work makes theoretical and empirical contributions in many areas.
There are some limitations in our proposed stochastic model. First, it does not consider the epidemic network topology. However, how to infer such networks is another challenging task. To the best of our knowledge, the best way to do so is to use the contact data among some infected patients to verify the results, but such data are not always available and can be difficult to collect due to many issues (e.g., privacy). This issue may be addressed by using other types of data, such as daily commute data extracted from social networks. Second, our proposed model achieves a prediction accuracy of only around 80 %. We need to further improve it to allow its full use in real applications. Third, the proposed model is only suitable for the situation in which the susceptible population (or subpopulation) maintains a relatively constant size and structure in a region. However, if the number of infected people in an epidemic is large or asymptomatic infection plays a central role (e.g., the malaria epidemic in Africa), the population factor should be taken into consideration in the model. Moreover, for a highly spatially heterogeneous outbreak (e.g., the Ebola epidemic) in which cases may seem to disappear due to reduced transmission in one area while growth may continue or rise in new locales, our proposed model may have problems in capturing these opposite dynamics in different regions. Fourth, because the proposed model is based on the classic SIR model, it only works in the situation in which the number of infected people grows exponentially. We will investigate resolutions to these limitations in our future work.
Notes
Declarations
Acknowledgements
This work is supported by Hong Kong Baptist University Strategic Development Fund and Hong Kong General Research Grant HKBU12202114.
Authors’ contributions
RXM, JML, and XW conceived and designed the experiments. RXM implemented the software. JML and XW analysed the data. All authors were involved in the manuscript preparation. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
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.
Authors’ Affiliations
References
 World Health Organization Pandemic (H1N1). 2009. http://www.who.int/csr/don/2010_01_08/en/index.html. Accessed 1 Aug 2010.
 Cohen J, Enserink M. As swine flu circles globe, scientists grapple with basic questions. Science. 2009; 324(5927):572–3.View ArticlePubMedGoogle Scholar
 Ferguson NM, Cummings DA, Fraser C, Cajka JC, Cooley PC, Burke DS. Strategies for mitigating an influenza pandemic. Nature. 2006; 442(7101):448–52.View ArticlePubMedGoogle Scholar
 Germann TC, Kadau K, Longini IM, Macken CA. Mitigation strategies for pandemic influenza in the United States. Proc Natl Acad Sci. 2006; 103(15):5935–940.View ArticlePubMedPubMed CentralGoogle Scholar
 Bailey NTJ, et al. The Mathematical Theory of Infectious Diseases and Its Applications. High Wycombe: Charles Griffin and Company Ltd; 1975.Google Scholar
 Anderson RM, May RM. Infectious Diseases of Humans vol. 1. Oxford: Oxford University Press; 1991.Google Scholar
 Heesterbeek J. Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. vol 5. New York City: Wiley; 2000.Google Scholar
 Li MY, Muldowney JS. Global stability for the SEIR model in epidemiology. Math Biosci. 1995; 125(2):155–64.View ArticlePubMedGoogle Scholar
 Kuznetsov YA, Piccardi C. Bifurcation analysis of periodic SEIR and SIR epidemic models. J Math Biol. 1994; 32(2):109–21.View ArticlePubMedGoogle Scholar
 Hethcote HW. Qualitative analyses of communicable disease models. Math Biosci. 1976; 28(3):335–56.View ArticleGoogle Scholar
 Becker NG, Britton T. Statistical studies of infectious disease incidence. J R Stat Soc Series B Stat Methodol. 2002; 61(2):287–307.View ArticleGoogle Scholar
 Ferrari MJ, Bjørnstad ON, Dobson AP. Estimation and inference of R0 of an infectious pathogen by a removal method. Math Biosci. 2005; 198(1):14–26.View ArticlePubMedGoogle Scholar
 Allen L. An introduction to stochastic epidemic models. Math Epidemiol. 2008; 1945:81–130.View ArticleGoogle Scholar
 Cooper BS, Pitman RJ, Edmunds WJ, Gay NJ, et al. Delaying the international spread of pandemic influenza. PLoS Med. 2006; 3(6):212.View ArticleGoogle Scholar
 Hufnagel L, Brockmann D, Geisel T. Forecast and control of epidemics in a globalized world. Proc Natl Acad Sci U S A. 2004; 101(42):15124–15129.View ArticlePubMedPubMed CentralGoogle Scholar
 Hollingsworth TD, Ferguson NM, Anderson RM. Will travel restrictions control the international spread of pandemic influenza?Nat Med. 2006; 12(5):497–9.View ArticlePubMedGoogle Scholar
 Keeling MJ, Woolhouse ME, Shaw DJ, Matthews L, ChaseTopping M, Haydon DT, Cornell SJ, Kappey J, Wilesmith J, Grenfell BT. Dynamics of the 2001 UK foot and mouth epidemic: stochastic dispersal in a heterogeneous landscape. Science. 2001; 294(5543):813–7.View ArticlePubMedGoogle Scholar
 Ferguson NM, Cummings DA, Cauchemez S, Fraser C, Riley S, Meeyai A, Iamsirithaworn S, Burke DS. Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005; 437(7056):209–14.View ArticlePubMedGoogle Scholar
 Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings DA, Halloran ME. Containing pandemic influenza at the source. Science. 2005; 309(5737):1083–1087.View ArticlePubMedGoogle Scholar
 PastorSatorras R, Vespignani A. Epidemic dynamics and endemic states in complex networks. Phys Rev E. 2001; 63(6):066117.View ArticleGoogle Scholar
 PastorSatorras R, Vespignani A. Epidemics and immunization in scalefree networks. 2002. arXiv preprint condmat/0205260.
 Kuperman M, Abramson G. Small world effect in an epidemiological model. Phys Rev Lett. 2001; 86(13):2909–912.View ArticlePubMedGoogle Scholar
 Newman ME, Jensen I, Ziff R. Percolation and epidemics in a twodimensional small world. Phys Rev E. 2002; 65(2):021904.View ArticleGoogle Scholar
 Boguná M, PastorSatorras R, Vespignani A. Absence of epidemic threshold in scalefree networks with degree correlations. Phys Rev Lett. 2003; 90(2):028701.View ArticlePubMedGoogle Scholar
 Grassberger P. On the critical behavior of the general epidemic process and dynamical percolation. Math Biosci. 1983; 63(2):157–72.View ArticleGoogle Scholar
 Watts DJ, Strogatz SH. Collective dynamics of smallworld networks. Nature. 1998; 393(6684):440–2.View ArticlePubMedGoogle Scholar
 Barabási AL, Albert R. Emergence of scaling in random networks. Science. 1999; 286(5439):509–12.View ArticlePubMedGoogle Scholar
 Erd, 6s P, Rényi A. On the evolution of random graphs. Publ Math Inst Hungar Acad Sci. 1960; 5:17–61.
 PastorSatorras R, Vespignani A. Epidemics and immunization in scalefree networks. In: Handbook of graphs and networks: from the genome to the internet. Hoboken: Wiley Online Library: 2005. p. 111–130.Google Scholar
 Rogers EM. Diffusion of Innovations. New York: Simon and Schuster; 2010.Google Scholar
 Leskovec J, Adamic LA, Huberman BA. The dynamics of viral marketing. ACM T Web. 2007; 1(1):5.View ArticleGoogle Scholar
 Kumar R, Novak J, Raghavan P, Tomkins A. On the bursty evolution of Blogspace. World Wide Web. 2005; 8(2):159–78.View ArticleGoogle Scholar
 Leskovec J, McGlohon M, Faloutsos C, Glance NS, Hurst M. Patterns of cascading behavior in large blog graphs. In: SDM, vol. 7. SIAM: 2007. p. 551–6.
 Salathé M, Kazandjieva M, Lee JW, Levis P, Feldman MW, Jones JH. A highresolution human contact network for infectious disease transmission. Proc Natl Acad Sci. 2010; 107(51):22020–2025.View ArticlePubMedPubMed CentralGoogle Scholar
 Stehlé J, Voirin N, Barrat A, Cattuto C, Isella L, Pinton JF, Quaggiotto M, Van den Broeck W, Régis C, Lina B, et al. Highresolution measurements of facetoface contact patterns in a primary school. PloS ONE. 2011; 6(8):23176.View ArticleGoogle Scholar
 Dickison M, Havlin S, Stanley HE. Epidemics on interconnected networks. Phys Rev E. 2012; 85(6):066109.View ArticleGoogle Scholar
 Hancock K, Veguilla V, Lu X, Zhong W, Butler EN, Sun H, Liu F, Dong L, DeVos JR, Gargiullo PM, et al. Crossreactive antibody responses to the 2009 pandemic H1N1 influenza virus. N Engl J Med. 2009; 361(20):1945–1952.View ArticlePubMedGoogle Scholar
 Riley S, Fraser C, Donnelly CA, Ghani AC, AbuRaddad LJ, Hedley AJ, Leung GM, Ho LM, Lam TH, Thach TQ, et al. Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions. Science. 2003; 300(5627):1961–1966.View ArticlePubMedGoogle Scholar
 Mills CE, Robins JM, Lipsitch M. Transmissibility of 1918 pandemic influenza. Nature. 2004; 432(7019):904–6.View ArticlePubMedGoogle Scholar
 Birrell PJ, Ketsetzis G, Gay NJ, Cooper BS, Presanis AM, Harris RJ, Charlett A, Zhang XS, White PJ, Pebody RG, et al. Bayesian modelling to unmask and predict influenza A/H1N1pdm dynamics in London. Proc Natl Acad Sci. 2011; 108(45):18238–18243.View ArticlePubMedPubMed CentralGoogle Scholar
 Towers S, Geisse KV, Zheng Y, Feng Z. Antiviral treatment for pandemic influenza: Assessing potential repercussions using a seasonally forced SIR model. J Theor Biol. 2011; 289:259–68.View ArticlePubMedGoogle Scholar
 Cauchemez S, Valleron AJ, Boelle PY, Flahault A, Ferguson NM. Estimating the impact of school closure on influenza transmission from sentinel data. Nature. 2008; 452(7188):750–4.View ArticlePubMedGoogle Scholar
 Cauchemez S, Donnelly CA, Reed C, Ghani AC, Fraser C, Kent CK, Finelli L, Ferguson NM. Household transmission of 2009 pandemic influenza A (H1N1) virus in the United States. N Engl J Med. 2009; 361(27):2619–627.View ArticlePubMedGoogle Scholar
 Lessler J, Reich NG, Cummings DA. Outbreak of 2009 pandemic influenza a (h1n1) at a New York city school. N Engl J Med. 2009; 361(27):2628–636.View ArticlePubMedGoogle Scholar
 Klick B, Nishiura H, Ng S, Fang VJ, Leung GM, Peiris JM, Cowling BJ. Transmissibility of seasonal and pandemic influenza in a cohort of households in Hong Kong in 2009. Epidemiology (Cambridge). 2011; 22(6):793.Google Scholar
 Yang Y, Sugimoto JD, Halloran ME, Basta NE, Chao DL, Matrajt L, Potter G, Kenah E, Longini IM. The transmissibility and control of pandemic influenza a (H1N1) virus. Science. 2009; 326(5953):729–33.View ArticlePubMedPubMed CentralGoogle Scholar
 Jin Z, Zhang J, Song LP, Sun GQ, Kan J, Zhu H. Modelling and analysis of influenza a (H1N1) on networks. BMC Public Health. 2011; 11(Suppl 1):9.Google Scholar
 World Health Organization. Ebola Situation Reports. 2015. Available at: http://apps.who.int/ebola/en/currentsituation/ebolasituationreport. Accessed 29 June 2015.
 Camacho A, Kucharski A, AkiSawyerr Y, White MA, Flasche S, Baguelin M, Pollington T, Carney JR, Glover R, Smout E, et al.Temporal changes in Ebola transmission in sierra leone and implications for control requirements: a realtime modelling study. PLoS Curr. 2015; 7:PMC4339317.Google Scholar
 Chowell D, CastilloChavez C, Krishna S, Qiu X, Anderson KS. Modelling the effect of early detection of Ebola. Lancet Infect Dis. 2015; 15(2):148–9.View ArticlePubMedGoogle Scholar
 Chowell G, Hengartner NW, CastilloChavez C, Fenimore PW, Hyman J. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda. J Theor Biol. 2004; 229(1):119–26.
 Chowell G, Viboud C, Hyman JM, Simonsen L. The Western Africa Ebola virus disease epidemic exhibits both global exponential and local polynomial growth rates. PLoS Curr. 2015. doi:10.1371/currents.outbreaks.8b55f4bad99ac5c5db3663e916803261.
 Fisman D, Khoo E, Tuite A. Early epidemic dynamics of the West African 2014 Ebola outbreak: estimates derived with a simple twoparameter model. PLoS Curr. 2014. doi:10.1371/currents.outbreaks.89c0d3783f36958d96ebbae97348d571.
 Gomes MF, y Piontti AP, Rossi L, Chao D, Longini I, Halloran ME, Vespignani A. Assessing the international spreading risk associated with the 2014 West African Ebola outbreak. PLoS Curr. 2014. doi:10.1371/currents.outbreaks.cd818f63d40e24aef769dda7df9e0da5.
 Althaus CL. Estimating the reproduction number of Ebola virus (EBOV) during the 2014 outbreak in West Africa. PLoS Curr. 2014. doi:10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288.
 Webb G, Browne C, Huo X, Seydi O, Seydi M, Magal P. A model of the 2014 Ebola epidemic in West Africa with contact tracing. PLoS Curr. 2014. doi:10.1371/currents.outbreaks.846b2a31ef37018b7d1126a9c8adf22a.
 Carroll MW, Matthews DA, Hiscox JA, Elmore MJ, Pollakis G, Rambaut A, Hewson R, GarcíaDorival I, Bore JA, Koundouno R, et al.Temporal and spatial analysis of the 20142015 Ebola virus outbreak in West Africa. Nature. 2015; 524(7563):97–101.View ArticlePubMedGoogle Scholar
 Chowell G, Nishiura H. Transmission dynamics and control of Ebola virus disease (EVD): a review. BMC Med. 2014; 12(1):196.View ArticlePubMedPubMed CentralGoogle Scholar
 van Dijk A, Aramini J, Edge G, Moore KM. Realtime surveillance for respiratory disease outbreaks, Ontario, Canada. Emerg Infect Diseases. 2009; 15(5):799.View ArticleGoogle Scholar