Comparison of methods to Estimate Basic Reproduction Number (R_{0}) of influenza, Using Canada 2009 and 201718 A (H1N1) Data
Roya Nikbakht^{1}, Mohammad Reza Baneshi^{2}, Abbas Bahrampour^{2}, Abolfazl Hosseinnataj^{2}
^{1} HIV/STI Surveillance Research Center, and WHO Collaborating Center for HIV Surveillance, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Department of Biostatistics and Epidemiology, Faculty of Health Kerman, Iran ^{2} Department of Biostatistics and Epidemiology, Faculty of Health, Modeling in Health Research Center, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman, Iran
Date of Submission  18Nov2018 
Date of Decision  13Mar2019 
Date of Acceptance  17May2019 
Date of Web Publication  24Jul2019 
Correspondence Address: Prof. Abbas Bahrampour Department of Biostatistics and Epidemiology, Faculty of Health, Modeling in Health Research Center, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman Iran
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/jrms.JRMS_888_18
Background: The basic reproduction number (R_{0}) has a key role in epidemics and can be utilized for preventing epidemics. In this study, different methods are used for estimating R_{0}'s and their vaccination coverage to find the formula with the best performance. Materials and Methods: We estimated R_{0}for cumulative cases count data from April 18 to July 6, 2009 and 352017 to 342018 weeks in Canada: maximum likelihood (ML), exponential growth rate (EG), timedependent reproduction numbers (TD), attack rate (AR), gammadistributed generation time (GT), and the final size of the epidemic. Gamma distribution with mean and standard deviation 3.6 ± 1.4 is used as GT. Results: The AR method obtained a R_{0 (}95% confidence interval [CI]) value of 1.116 (1.1163, 1.1165) and an EG (95%CI) value of 1.46 (1.41, 1.52). The R_{0}(95%CI) estimate was 1.42 (1.27, 1.57) for the obtained ML, 1.71 (1.12, 2.03) for the obtained TD, 1.49 (1.0, 1.97) for the gammadistributed GT, and 1.00 (0.91, 1.09) for the final size of the epidemic. The minimum and maximum vaccination coverage were related to AR and TD methods, respectively, where the TD method has minimum mean squared error (MSE). Finally, the R_{0}(95%CI) for 2018 data was 1.52 (1.11, 1.94) by TD method, and vaccination coverage was estimated as 34.2%. Conclusion: For the purposes of our study, the estimation of TD was the most useful tool for computing the R_{0}, because it has the minimum MSE. The estimation R_{0}>_{}1 indicating that the epidemic has occurred. Thus, it is required to vaccinate at least 41.5% to prevent and control the next epidemic.
Keywords: Basic reproduction number, influenza A virus, vaccination coverage
How to cite this article: Nikbakht R, Baneshi MR, Bahrampour A, Hosseinnataj A. Comparison of methods to Estimate Basic Reproduction Number (R_{0}) of influenza, Using Canada 2009 and 201718 A (H1N1) Data. J Res Med Sci 2019;24:67 
How to cite this URL: Nikbakht R, Baneshi MR, Bahrampour A, Hosseinnataj A. Comparison of methods to Estimate Basic Reproduction Number (R_{0}) of influenza, Using Canada 2009 and 201718 A (H1N1) Data. J Res Med Sci [serial online] 2019 [cited 2019 Oct 20];24:67. Available from: http://www.jmsjournal.net/text.asp?2019/24/1/67/263366 
Introduction   
Pandemic influenza, a global outbreak, defines as spreading influenza virus between peoples (with little or lack of immunity) over a wide geographic field.^{[1]} In the 20^{th} century, three pandemics of influenza happened which were “Spanish flu,” “Asian flu,” and “Hong Kong flu” in the years “1918–1919,” “1957–1958,” and “1968–1969,” respectively.^{[2]} In early 2009, H1N1 influenza at first occurred in Mexico and the United States and speared rapidly worldwide (>200 countries involved).^{[3],[4]} The influenza virus can spread among people by direct contact (a cough, sneeze or talk), inhalation of virusladen aerosols, and touch fomites (contaminated objects) that has the flu virus.^{[5],[6]} The most affected groups for developing flurelated complications are children, pregnant women, elders (adults older than 64yearold), and persons with a specific disease (chronic pulmonary disease, chronic heart disease, diabetes, etc.,).^{[7],[8]} The mortality and morbidity related to the annual influenza in the worldwide estimated approximately one million people, a considerable number.^{[9]} For example, the number of deaths for “United States flu (2009)” reported 12,469 and for “Asian flu” was 1–4 million.^{[10],[11]} Therefore, controlling and preventing the epidemic of influenza is an important issue. The basic reproduction number (R_{0}) is an important metric that used for measuring the vaccination coverage (to prevent epidemic), eradicating an infectious disease, controlling and immunizing the disease which is defined the mean number of secondary infections generated by a single infectious individual in a fully susceptible population without immunity and interventions.^{[12]} In particular, the R_{0} determines whether an infection spreads through a population.^{[13]} The basic reproduction number or threshold parameter applied for determining the critical immunity coverage can be a real number greater than, less than, or equal to one. The disease will fade out when R_{0}<1 and an epidemic will occur (the infection will grow) if R_{0}≥1, showing an endemic in the population.^{[13],[14]}
Since the R_{0} has a key role in measuring the transmission of diseases and is crucial in preventing epidemics, thus it is important to know which methods and formulas to apply to estimate R_{0} and have better performance. We estimate the R_{0} and its related vaccination coverage for Canadian influenza data during 2009 and 2017–2018.
Materials and Methods   
Objectives
In this study, we reviewed the investigated methods and formulas used for estimating the R_{0} of influenza in various published research papers from 1954 to 2017. After a scientific systematic review on R_{0}, we found out that there are many basic reproduction formula which are applied for determining the vaccination coverage so it is necessary to characterize a formula which gives more accurate result to use in vaccination strategies which leads to optimize the costs. We extracted more commonlyutilized formulas [Appendix Table 1]. We considered six common formulas and applied them to real data to determine which formula most closely approximates the real epidemic threshold parameter with high efficacy.
Then, R_{0}s and related vaccination coverage of these methods was estimated for a secondary real data of Canadian influenza (2009). The calculated R_{0} was compared with R_{0} of the Canadian paper^{[15]} and also simulations were performed. Finally, the best method was chosen based on mean squared error (MSE), then R_{0} calculated by selected method for the H1N1 Canadian data in the 35^{th} week in 2017–34^{th} week in 2018.
Data
In Canada, circulating of influenza A virus is very common. The data sets in this study were obtained from the Public Health Agency of Canada (PHAC) website^{[16]} and the last FluWatch weekly report of the 2017–2018 influenza surveillance season achieved from the Respiratory Virus Detections in Canada Report website.^{[17]}
The total number of patients was 927 during the 2009 influenza season which were based on month/day and the number of new cases was 1280 for Canada 2017–2018 H1N1 data which report every Thursday in Canada. We fitted all the six models to Canadian 2009 pH 1N1 cumulative cases data.^{[16]} Then, the best model was applied to the data of Canada (34^{th} week in 2017 to 34^{th} week in 2018).^{[17]}
Statistical analysis
The models used in this article included the Richard model, attack rate (AR), exponential growth rate (EG), maximum likelihood (ML), timedependent reproduction numbers (TD), gammadistributed generation time (GT), and R_{0} using the final size of the epidemic. The above mentioned methods were applied for estimating R_{0} using R software (R_{0} package and programming). R software was created by Ross Ihaka and Robert Gentleman at the University of Auckland, New Zealand, and is currently developed by the R Development Core Team (of which Chambers is a member).
Generation time
The timegap between infection of a primary case and infection of a secondary case that is generated by the primary case.^{[18]}
The attack rate
The R_{0} can be described by the AR with the following formula:
where AR defines the ratio of the people generating an infection disease and S_{0} show the initial susceptible ratio.^{[19]}
The exponential growth rate
The following formula was applied for computing the R:
In this formula, M is the momentgenerating function of the GT.^{[20]} The parameter r is determined by the Poisson regression. Furthermore, the parameter w is GT.
The maximum likelihood
Let N_{0}, N_{1},...., N_{T} identify incident cases over sequential time. The loglikelihood function is:
where
and R is the maximum value of the loglikelihood function.^{[21]} Furthermore, the parameter w is estimated by maximizing loglikelihood is GT.
Timedependent reproduction numbers
In this method, R_{t} is computed by averaging R_{j}, which is the mean of all transmission networks corresponding to the cases observed.^{[22]}
where
And
Consider that person i and person j are in times t_{i} and t_{j}, respectively, then displays the probability of infection transmission from person j to person i so R_{t} compute by averaging all R_{j} which is the mean of all transmission networks correspondent with the cases that observed.
The gammadistributed generation time
The number of cases on the day “t,” denoted by n_{t} in (t_{1}, t_{2}) grows exponentially where
n_{t}= n_{t1} exp (r[t – t 1]) (9)
And
The EG denotes by r. The mean and standard deviation of the GT are μ and σ, respectively, where a = μ^{2}/σ^{2} and b = μ/σ^{2}.^{[23]}
R_{0} using the final size of the epidemic
The R_{0} can be estimated with the below formula:
where the total population at risk and total number of infections are denoted by N and C, respectively.^{[24]}
Vaccination coverage
The vaccination coverage is computed by the basic reproduction number with formula:
which shows the proportion of peoples who should be received the vaccine.^{[7]}
Comparison of methods
For exploring the closeness of the estimation of the mentioned methods to the actual R_{0}s and comparing them with each other, we applied 10000 times simulation for each formula based on the Canada data. The epidemics were simulated with the following properties. The distribution of the GT was considered gamma with the mean of 3.6 and standard deviation of 1.4. According to real data (the Canada data), the length of the epidemic was 80 days. Moreover, the peak value (the threshold value for the incidence before epidemics begin decreasing) for the Canada data occurred in the day 54. Therefore, we applied the value equal to 54 for the peak value in the simulation command [For details, see the simulation command under [Table 1] in the results section]. Simulation of the basic reproduction number was made with above characteristics and the MSE was calculated for evaluating the performance of models with below formula. The lowest MSE value corresponds to the method which fitted the data best.  Table 1: The simulated R_{0}s and their 95% confidence interval for each method
Click here to view 
Results   
Canadian 2009 H1N1 influenza data
We fitted the six models to the daily dataset of Canada, throughout the 80day period of the studies. All dates of the Canada data were based on month/day form 18 April, 2009 to 6 July, 2009. Moreover, the number of infected people was plotted as frequency [Figure 1].  Figure 1: The incidence case counts influenza data of Canada during 18 April, 2009–6 July, 2009
Click here to view 
In order to demonstrate the difference in modeling with various formulas, the result of the Richard model (presented in Hsieh's study)^{[15]} as well as the results of the other six models are presented in [Table 2]. The reported R_{0}(95% confidence interval [CI]) (vaccination coverage%) using the Richard model was 1.68 (1.45, 1.91) (40.47) that means every person infected 1.68 other people on average during the infection period. Note that, R_{0}(95%CI) (vaccination coverage%) for the estimation of TD (1.71 [1.12, 2.03] (41.52)) was clearly close to R_{0} for the Richard model. The second method with the closest R_{0}(95%CI) to that of the Richard model was the gammadistributed GT (1.49 [1.0, 1.97] (32.88)). On the other hand, the computed R_{0}(95%CI) using the EG was 1.46 (1.41, 1.52) (31.51). The ML method revealed that the calculated R_{0}(95%s CI) for this model was different from that for the Richard model (1.42 [1.27, 1.57] [29.58]). In addition, the estimated R_{0}(95% CI) (vaccination coverage%) by the AR with two approaches was 1.000388 (1.000383, 1.000392) (0.04) and 1.1164 (1.1163, 1.1165) (10.43). The minimum computed R_{0}(95% CI) was related to the estimation of the final size of the epidemic obtained as 1.0 (0.91, 1.09). The estimates of vaccination coverage for the six methods were vary. The lowest and highest vaccination coverage values in this setting were associated with AR and TD methods, respectively.  Table 2: The Reproduction number estimation by the different methods for the Canada data (2009)
Click here to view 
In order to compare the mentioned models to find the formula with better fit to the actual values, we conducted a simulation with R software and calculated R_{0} based on the six models reported in [Table 2]. We used gamma distribution for the GT with the mean of 3.6 and standard deviation of 1.4. The peak value determined right over the original data were equal to 54. Then, using the above parameters, the simulation was implemented and R_{0} was computed for each method. The simulation results for comparing the quality of the six methods are represented in [Table 1] and [Figure 2]. In order to carry out the simulation, the number of runs to achieve the R_{0} was 10000.  Figure 2: The plots of the actual and simulated R_{0}compared for each method
Click here to view 
The results, given in [Table 1], indicated that there were differences between the actual and simulated R_{0}; however, the TD method had the closest value to the R_{0} calculated from the simulation compared to the other methods. Surprisingly, some variation was considered for the ML estimations when the actual values were equal to one, between one and two and greater than two. In the ML method, we found that the simulated R_{0} for small values was very close to that for the actual values when the actual values were between 1.42 and 1.71; while the simulated R_{0} for large values was very different from that for the actual values. For the gammadistributed GT approach, the simulated R_{0} grew out of the actual values for values close to one. In contrast, the results showed that the computed values for R_{0} in the simulated system were slightly greater than the actual values when we applied R_{0} between 1.42 and 2. By following the same interpretation, we can infer that the EG method had a small variation for small R_{0} values (1.4 < R_{0}<2). On the other hand, the R_{0} estimations using the EG diverged from the actual R_{0} but was not significant. Finally, the computed R_{0} by the AR and final size of the epidemic methods seemed likely to reflect stability for all R_{0}s. In particular, for the latest assumed R_{0}s, the estimated R_{0} was equal to one.
We also plotted [Figure 2] the actual R_{0} and simulated R_{0} based on six methods with the parameters described in [Table 1]. For evaluating the performance of models, we computed MSE for all methods [Table 3]. The TD method had the lowest MSE value in comparison to other methods. The MSE of AR and final size of the epidemic methods was very varied. In addition, MSE of ML, EG, and gammadistributed GT methods were also calculated. For ML, EG, and gammadistributed GT, the mean of MSE of all points were 4.85, 3.81, and 3.31, respectively. As noted above, the TD introduced the approach with the nearest estimation to the actual R_{0} based on MSE criterion.  Table 3: Mean squared error of reproduction number estimation for each method
Click here to view 
We also performed a sensitivity analysis with the incidence data of Canada on the GT with the gamma distribution [Figure 3]. The sensitivity analysis demonstrated that R_{0}(95% CI) for the mean GT (days) of 3.6 and 4.9 was estimated as 1.47 (1.41, 1.53) and 1.67 (1.58, 1.76). Thus, the computed R_{0} was approximately near that of the Richard and TD methods when the mean GT was equal to 4.9.  Figure 3: Sensitivity of R_{0}to mean generation time to select the generation time
Click here to view 
Canadian 2017–2018 H1N1 influenza data
The incidence data are reported based on week/year from the 35^{th} week in 2017 to the 34^{th} week in 2018. Peak value for this data has occurred in the 12^{th} week in 2018 after starting the epidemics. The number of infected cases is plotted in [Figure 4].  Figure 4: The incidence case counts influenza data of Canada from the 35^{th} week in 2017 to the 34^{th} week in 2018
Click here to view 
For the given data, R_{0}(95% CI) and vaccination coverage based on TD method was computed. Indeed, we found that the estimated R_{0} by TD method was (1.52 95% CI: 1.11, 1.94). In addition, the estimates of vaccination coverage were 34.2% for 2017–2018.
Discussion   
We implemented six methods (the ML, EG, TD, AR, gammadistributed and final size of the epidemic), which permitted the estimation of the R_{0} as key parameters of the epidemic based on the A/H1N1 Influenza cumulative case counts data in Canada (2009). The R_{0} for the ML, EG, TD, AR, gammadistributed and final size of the epidemic methods were estimated 1.42, 1.46, 1.71, 1.116, 1.49, and 1.0, respectively. In most cases, the R_{0} was greater than unity; hence, the epidemic outbreak was observed. In addition, the computed R_{0} for Canadian data (2018) by TD method was greater than one indicating that an epidemic occurred in Canada (R_{0}>1). Thus, it seems necessary to consider appropriate solutions in order to control, decrease and prevent the epidemic or pandemic of influenza. One of the most effective methods to protect people against influenza is vaccination that can be determined by using R_{0}(vaccination coverage = 1 − 1/R_{0}). On the other hand, annual influenza vaccination in the highrisk groups such as elderly people, ill person, pregnant woman, and children can reduce mortality rate. In addition, vaccination can also reduce the incidence of disease, cost, exacerbations of the disease, and hospitalizations. The vaccination coverage for Canada (2009) ranged between 10.43 and 41.52 using various methods and this value was 34.2% for 2017–2018 influenza Canada data.
Moreover, we performed a simulation using R software for several R_{0} and obtained their estimates based on the epidemic data of Canada (2009) for the six methods. The computed R_{0} in the TD method was nearly the same as the actual R_{0} based on MSE criterion. Comparing the simulation results from the ML, gammadistributed GT and EG methods showed variation for different values of the actual R_{0}; however, some of the calculated R_{0}s applying the simulation were close to the actual values. For the most actual R_{0}, the simulated R_{0} by the AR and final size of the epidemic methods was equal to one. Whereas these type of modeling approaches are not able to differentiate between various R_{0}. We believe that this may correspond to the small number of the infected cases compared to the susceptible cases.
Note that, our basic reproduction number estimated using the TD method was consistent with that derived from the Richard model in the Canadian papers.^{[15]} Not only the simulated R_{0} for the value 1.68 almost agreed with that of the TD approach but also the other simulated R_{0} by the TD method was nearly consistent with the actual R_{0}. In other words, the lowest MSE values were obtained for TD method.
From the methods reviewed in [Appendix Table 2], which can be applied to estimate the R_{0}, the approaches presented in [Table 1] fitted to the cumulative cases data. All the methods reviewed in this paper, as any modeling techniques, had advantageous and disadvantageous. One of the strengths of this study is to review all studies done related to influenza and then selected some of the frequently used model and determine their strengths and weaknesses; seven of them used for the R_{0} estimation in the Canada data, as shown in [Table 1], are explained in details in [Table 4].  Table 4: Limitation and power of the methods used for the cumulative case counts data
Click here to view 
Regarding [Table 4], it seemed that the TD, ML and EG methods had superiority compared to the other methods. These models were used by researchers to estimate R_{0} of influenza.
Some studies estimated the R_{0} from influenza data using different models and compared the results. Obadia et al. obtained estimates of R_{0} from the “Germany 1918” epidemic data based on five approaches which including the AR, ML, sequential Bayesian and TD methods. In addition, comparing results from different methods showed that the biased ML and TD methods were least.^{[30]} Another study applied four different methods (the EG, simple susceptibleexposedinfectiousrecovered [SEIR], more complex SEIRtype model, and ML model) in order to compare these estimation approaches. The EG had large uncertainty while ML had a consistent estimate with the estimate of the autumn wave.^{[20]} In general, the TD had a good fit on the data as confirmed with the Richard model and MSE criterion.
A weakness of this study is that the 2009 Canada data have been used for comparing methods, which looks old. The reason for this, is comparing R_{0} with pervious article^{[15]} and comparing the methods with the actual values which are exist on this data in the mentioned paper. Finally, a more comprehensive study for influenza as an annual national disaster using new method such as Bayesian is needed that we are going to do in the future research.
Conclusion   
Awareness of the basic reproduction number of influenza is useful for calculating vaccination coverage and then applying vaccine strategy. Therefore, it is necessary to know the method which has better performance for influenza data that our results showed the TD method is preferred. One advantage of the TD method in compared to the other methods was that it was useful for computing the R_{0} regarding the real cumulative case count data. Another advantage of the mentioned modeling was that it did not require extensive, detailed data as well as more parameters to calculate the basic reproduction number. Therefore, we recommend using this method in order to estimate the basic reproduction number.
Acknowledgments
This research is part of Roya Nikbakht's PhD dissertation. We are very grateful to the Public Health Agency of Canada (PHAC) website and the Respiratory Virus Detections in Canada Report website to provide online database and access to use their data.
Financial support and sponsorship
Nil.
Conflicts of interest
There are no conflicts of interest.
Appendix   
Search Strategy   
In order to review the literature on basic reproduction number of influenza, we searched in the electronic databases such as the web of knowledge, PubMed, EMBASE, and Google Scholar to find published papers between 1954 and 2017. The medical subject heading was applied to find a wide range of keywords that had a maximum sensitivity. The following keywords were searched: influenza, human, and reproduction number. In detail searched keywords were (“influenza, human”[MeSH Terms] OR (“influenza”[All Fields] AND “human”[All Fields]) OR “human influenza”[All Fields] OR “influenza”[All Fields]) AND ((“reproduction”[MeSH Terms] OR “reproduction”[All Fields]) AND number [All Fields]).
Study Selection   
Two reviewers independently extracted relevant studies from the keywords search. All types of original articles were investigated. The studies which included “influenza reproduction number” in their titles or abstracts were included. The irrelevant articles, based on the title and abstract evaluation, were excluded. Moreover, we eliminated the duplicated articles to determine unique studies. Animal studies and human studies that included special populations such as pregnant women and schizophrenia were excluded. We then extracted data and formulas from the full text of the included studies.
[Figure 1] shows the search strategy, through which 1213 papers were obtained in the initial round. The number of the retained papers was 910, which estimated R_{0} for epidemic or pandemic influenza with A/H1N1, A/H1N5, H1N2, H1N3, H5N1, pH 1N1, A/H3N2, influenza B, A (H7N9), Spanish flu, H2N2, H3N2, AH1, AH3, A (H5N1), and Asian flu. The number of papers identified through other sources was 5. Overall, 89 papers presented the basic reproduction number estimation and its formula, as summarized in [Table 1].
In addition, detailed information of the study characteristics provided in the systematic review is given in [Table 2], of which 10 studies were taken into consideration. In some of the studies, pH1N1, A (H1N1), A (H3N2), type B, and A (H7N9) were reported as types of influenza. The models used for estimating R_{0} in these 8 studies were the multicontrol measure, growth rate of exponential, and multiphase Richards. In several of the studies, laboratoryconfirmed cases were investigated for determining the reproduction number of influenza. Maximum, minimum and median of the reproduction number were 10.03 (in Mainland China), 0.08 (in China) and 1.39, respectively. The reproduction number of the influenza type A (H1N1) in Taiwan (2013) and Mexico was reported 1.54 (95% confidence interval [CI]: 0.22–8.88) and 1.69 (95% CI: 1.65–1.73), respectively. For A (H7N9), the reproduction number and its 95% CI in China for the first wave was estimated 0.27 (0.14, 0.44).
References   
 Ajelli M, Merler S. Transmission potential and design of adequate control measures for marburg hemorrhagic fever. PLoS One 2012;7:e50948.
 Ajelli M, Poletti P, Melegaro A, Merler S. The role of different social contexts in shaping influenza transmission during the 2009 pandemic. Sci Rep 2014;4:7218.
 Andreasen V, Viboud C, Simonsen L. Epidemiologic characterization of the 1918 influenza pandemic summer wave in Copenhagen: Implications for pandemic control strategies. J Infect Dis 2008;197:2708.
 Barakat A, Ihazmad H, El Falaki F, Tempia S, Cherkaoui I, El Aouad R, et al. 2009 pandemic influenza A virus subtype H1N1 in Morocco, 20092010: Epidemiology, transmissibility, and factors associated with fatal cases. J Infect Dis 2012;206 Suppl 1:S94100.
 Buckley D, Bulger D. Estimation of the reproductive number for the 2009 pandemic H1N1 influenza in rural and metropolitan New South Wales. Aust J Rural Health 2011;19:5963.
 Caley P, Philp DJ, McCracken K. Quantifying social distancing arising from pandemic influenza. J R Soc Interface 2008;5:6319.
 Cauchemez S, Bhattarai A, Marchbanks TL, Fagan RP, Ostroff S, Ferguson NM, et al. Role of social networks in shaping disease transmission during a community outbreak of 2009 H1N1 pandemic influenza. Proc Natl Acad Sci U S A 2011;108:282530.
 Chen SC, Liao CM. Probabilistic indoor transmission modeling for influenza (sub) type viruses. J Infect 2010;60:2635.
 Chen SC, Liao CM. Costeffectiveness of influenza control measures: A dynamic transmission modelbased analysis. Epidemiol Infect 2013;141:258194.
 Chen SC, Chio CP, Jou LJ, Liao CM. Viral kinetics and exhaled droplet size affect indoor transmission dynamics of influenza infection. Indoor Air 2009;19:40113.
 Chen SC, Hsieh NH, You SH, Wang CH, Liao CM. Behavioural response in educated young adults towards influenza A (H1N1) pdm09. Epidemiol Infect 2015;143:184657.
 Chen TM, Chen QP, Liu RC, Szot A, Chen SL, Zhao J, et al. The transmissibility estimation of influenza with early stage data of smallscale outbreaks in Changsha, China, 20052013. Epidemiol Infect 2017;145:42433.
 Chen T, Chen T, Liu R, Xu C, Wang D, Chen F, et al. Transmissibility of the influenza virus during influenza outbreaks and related asymptomatic infection in Mainland China, 20052013. PLoS One 2016;11:e0166180.
 Cheng YH, Liao CM. Modeling control measure effects to reduce indoor transmission of pandemic H1N1 2009 virus. Build Environ 2013;63:119.
 Cheng YH, Wang CH, You SH, Hsieh NH, Chen WY, Chio CP, et al. Assessing coughinginduced influenza droplet transmission and implications for infection risk control. Epidemiol Infect 2016;144:33345.
 Chong KC, Wang X, Liu S, Cai J, Su X, Zee BC, et al. Interpreting the transmissibility of the avian influenza A (H7N9) infection from 2013 to 2015 in Zhejiang province, China. Epidemiol Infect 2016;144:158491.
 Chong KC, Zee BC, Wang MH. A statistical method utilizing information of imported cases to estimate the transmissibility for an influenza pandemic. BMC Med Res Methodol 2017;17:31.
 Chowell G, Ammon CE, Hengartner NW, Hyman JM. Estimating the reproduction number from the initial phase of the Spanish flu pandemic waves in Geneva, Switzerland. Math Biosci Eng 2007;4:45770.
 Chowell G, Bettencourt LM, Johnson N, Alonso WJ, Viboud C. The 19181919 influenza pandemic in England and wales: Spatial patterns in transmissibility and mortality impact. Proc Biol Sci 2008;275:5019.
 Chowell G, DiazDueñas P, Miller JC, AlcazarVelazco A, Hyman JM, Fenimore PW, et al. Estimation of the reproduction number of dengue fever from spatial epidemic data. Math Biosci 2007;208:57189.
 Chowell G, Nishiura H, Bettencourt LM. Comparative estimation of the reproduction number for pandemic influenza from daily case notification data. J R Soc Interface 2007;4:15566.
 Chowell G, Simonsen L, Towers S, Miller MA, Viboud C. Transmission potential of influenza A/H7N9, February to May 2013, China. BMC Med 2013;11:214.
 Chowell G, Towers S, Viboud C, Fuentes R, Sotomayor V, Simonsen L, et al. The influence of climatic conditions on the transmission dynamics of the 2009 A/H1N1 influenza pandemic in Chile. BMC Infect Dis 2012;12:298.
 Chowell G, Viboud C, Munayco CV, Gómez J, Simonsen L, Miller MA, et al. Spatial and temporal characteristics of the 2009 A/H1N1 influenza pandemic in Peru. PLoS One 2011;6:e21287.
 Chowell G, Viboud C, Simonsen L, Miller M, Alonso WJ. The reproduction number of seasonal influenza epidemics in Brazil, 19962006. Proc Biol Sci 2010;277:185766.
 Cowling BJ, Lau MS, Ho LM, Chuang SK, Tsang T, Liu SH, et al. The effective reproduction number of pandemic influenza: Prospective estimation. Epidemiology 2010;21:8426.
 Dorigatti I, Cauchemez S, Pugliese A, Ferguson NM. A new approach to characterising infectious disease transmission dynamics from sentinel surveillance: Application to the Italian 20092010 A/H1N1 influenza pandemic. Epidemics 2012;4:921.
 Earn DJ, Andrews PW, Bolker BM. Populationlevel effects of suppressing fever. Proc Biol Sci 2014;281:20132570.
 Ejima K, Aihara K, Nishiura H. The impact of model building on the transmission dynamics under vaccination: Observable (symptombased) versus unobservable (contagiousnessdependent) approaches. PLoS One 2013;8:e62062.
 Fielding JE, Kelly HA, Glass K. Transmission of the first influenza A (H1N1) pdm09 pandemic wave in Australia was driven by undetected infections: Pandemic response implications. PLoS One 2015;10:e0144331.
 Fierro A, Liccardo A. A simple stochastic lattice gas model for H1N1 pandemic. Application to the Italian epidemiological data. Eur Phys J E Soft Matter 2011;34:11.
 Fraser C, Cummings DA, Klinkenberg D, Burke DS, Ferguson NM. Influenza transmission in households during the 1918 pandemic. Am J Epidemiol 2011;174:50514.
 Furuya H. Risk of transmission of airborne infection during train commute based on mathematical model. Environ Health Prev Med 2007;12:7883.
 Glass K, Mercer GN, Nishiura H, McBryde ES, Becker NG. Estimating reproduction numbers for adults and children from case data. J R Soc Interface 2011;8:124859.
 Gran JM, Iversen B, Hungnes O, Aalen OO. Estimating influenzarelated excess mortality and reproduction numbers for seasonal influenza in Norway, 19752004. Epidemiol Infect 2010;138:155968.
 Gurav YK, Chadha MS, Tandale BV, Potdar VA, Pawar SD, Shil P, et al. Influenza A (H1N1) pdm09 outbreak detected in interseasonal months during the surveillance of influenzalike illness in Pune, India, 20122015. Epidemiol Infect 2017;145:1898909.
 Haghdoost A, Baneshi MR, Zolala F, Farvahari S, Safizadeh H. Estimation of basic reproductive number of flulike syndrome in a primary school in Iran. Int J Prev Med 2012;3:40813.
 Hens N, Calatayud L, Kurkela S, Tamme T, Wallinga J. Robust reconstruction and analysis of outbreak data: Influenza A (H1N1) v transmission in a schoolbased population. Am J Epidemiol 2012;176:196203.
 Hens N, Van Ranst M, Aerts M, Robesyn E, Van Damme P, Beutels P, et al. Estimating the effective reproduction number for pandemic influenza from notification data made publicly available in real time: A multicountry analysis for influenza A/H1N1v 2009. Vaccine 2011;29:896904.
 Hsieh YH. Age groups and spread of influenza: Implications for vaccination strategy. BMC Infect Dis 2010;10:106.
 Hsieh YH, Cheng KF, Wu TN, Li TC, Chen CY, Chen JH, et al. Transmissibility and temporal changes of 2009 pH 1N1 pandemic during summer and fall/winter waves. BMC Infect Dis 2011;11:332.
 Hsieh YH, Fisman DN, Wu J. On epidemic modeling in real time: An application to the 2009 novel A (H1N1) influenza outbreak in Canada. BMC Res Notes 2010;3:283.
 Hsieh YH, Huang HM, Lan YC. On temporal patterns and circulation of influenza virus strains in Taiwan, 20082014: Implications of 2009 pH 1N1 pandemic. PLoS One 2016;11:e0154695.
 Hsieh YH, Ma S, Velasco Hernandez JX, Lee VJ, Lim WY. Early outbreak of 2009 influenza A (H1N1) in Mexico prior to identification of pH 1N1 virus. PLoS One 2011;6:e23853.
 Hsieh YH. Pandemic influenza A (H1N1) during winter influenza season in the Southern hemisphere. Influenza Other Respir Viruses 2010;4:18797.
 Inaba H, Nishiura H. The statereproduction number for a multistate class age structured epidemic system and its application to the asymptomatic transmission model. Math Biosci 2008;216:7789.
 Jackson C, Vynnycky E, Mangtani P. Estimates of the transmissibility of the 1968 (Hong Kong) influenza pandemic: Evidence of increased transmissibility between successive waves. Am J Epidemiol 2010;171:46578.
 Kadi AS, Avaradi SR. A bayesian inferential approach to quantify the transmission intensity of disease outbreak. Comput Math Methods Med 2015;2015:256319.
 Kelly HA, Mercer GN, Fielding JE, Dowse GK, Glass K, Carcione D, et al. Pandemic (H1N1) 2009 influenza community transmission was established in one Australian state when the virus was first identified in North America. PLoS One 2010;5:e11341.
 Kelso JK, Halder N, Postma MJ, Milne GJ. Economic analysis of pandemic influenza mitigation strategies for five pandemic severity categories. BMC Public Health 2013;13:211.
 Kim K, Omori R, Ito K. Inferring epidemiological dynamics of infectious diseases using Tajima's D statistic on nucleotide sequences of pathogens. Epidemics 2017;21:219.
 Kucharski AJ, Edmunds WJ. Characterizing the transmission potential of zoonotic infections from minor outbreaks. PLoS Comput Biol 2015;11:e1004154.
 Lin CJ, Deger KA, Tien JH. Modeling the tradeoff between transmissibility and contact in infectious disease dynamics. Math Biosci 2016;277:1524.
 Liu R, Leung RK, Chen T, Zhang X, Chen F, Chen S, et al. The effectiveness of agespecific isolation policies on epidemics of influenza A (H1N1) in a large city in central South China. PLoS One 2015;10:e0132588.
 Liu W, Tang S, Xiao Y. Model selection and evaluation based on emerging infectious disease data sets including A/H1N1 and Ebola. Comput Math Methods Med 2015;2015:207105.
 Marquetoux N, Paul M, Wongnarkpet S, Poolkhet C, Thanapongtharm W, Roger F, et al. Estimating spatial and temporal variations of the reproduction number for highly pathogenic avian influenza H5N1 epidemic in Thailand. Prev Vet Med 2012;106:14351.
 Massad E, Burattini MN, Coutinho FA, Lopez LF. The 1918 influenza A epidemic in the city of São Paulo, Brazil. Med Hypotheses 2007;68:4425.
 Mathews JD, McBryde ES, McVernon J, Pallaghy PK, McCaw JM. Prior immunity helps to explain wavelike behaviour of pandemic influenza in 19189. BMC Infect Dis 2010;10:128.
 Meeyai A, Cooper B, Coker R, PanNgum W, Akarasewi P, Iamsirithaworn S, et al. Pandemic influenza H1N1 2009 in Thailand. WHO South East Asia J Public Health 2012;1:5968.
 Modchang C, Iamsirithaworn S, Auewarakul P, Triampo. A modeling study of school closure to reduce influenza transmission: A case study of an influenza A (H1N1) outbreak in a private Thai school. Math Comput Model 2012;55:102133.
 MostaçoGuidolin LC, Bowman CS, Greer AL, Fisman DN, Moghadas SM. Transmissibility of the 2009 H1N1 pandemic in remote and isolated Canadian communities: A modelling study. BMJ Open 2012;2. pii: e001614.
 Nishiura H. Time variations in the transmissibility of pandemic influenza in Prussia, Germany, from 191819. Theor Biol Med Model 2007;4:20.
 Nishiura H, Chowell G, Safan M, CastilloChavez C. Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009. Theor Biol Med Model 2010;7:1.
 Nishiura H, Mizumoto K, Ejima K. How to interpret the transmissibility of novel influenza A (H7N9): An analysis of initial epidemiological data of human cases from China. Theor Biol Med Model 2013;10:30.
 Nyirenda M, Omori R, Tessmer HL, Arimura H, Ito K. Estimating the lineage dynamics of human influenza B viruses. PLoS One 2016;11:e0166107.
 Obadia T, Haneef R, Boëlle PY. The R0 package: A toolbox to estimate reproduction numbers for epidemic outbreaks. BMC Med Inform Decis Mak 2012;12:147.
 Pamaran RR, Kamigaki T, Hewe TT, Flores KM, Mercado ES, Alday PP, et al. Epidemiological characterization of influenza A (H1N1) pdm09 cases from 2009 to 2010 in Baguio city, the Philippines. PLoS One 2013;8:e79916.
 Rizzo C, Ajelli M, Merler S, Pugliese A, Barbetta I, Salmaso S, et al. Epidemiology and transmission dynamics of the 191819 pandemic influenza in florence, Italy. Vaccine 2011;29 Suppl 2:B2732.
 Roberts MG, Nishiura H. Early estimation of the reproduction number in the presence of imported cases: Pandemic influenza H1N12009 in New Zealand. PLoS One 2011;6:e17835.
 Roll U, Yaari R, Katriel G, Barnea O, Stone L, Mendelson E, et al. Onset of a pandemic: Characterizing the initial phase of the swine flu (H1N1) epidemic in Israel. BMC Infect Dis 2011;11:92.
 Sattenspiel L. Regional patterns of mortality during the 1918 influenza pandemic in Newfoundland. Vaccine 2011;29 Suppl 2:B337.
 Sertsou G, Wilson N, Baker M, Nelson P, Roberts MG. Key transmission parameters of an institutional outbreak during the 1918 influenza pandemic estimated by mathematical modelling. Theor Biol Med Model 2006;3:38.
 Ali ST, Kadi AS, Ferguson NM. Transmission dynamics of the 2009 influenza A (H1N1) pandemic in India: The impact of holidayrelated school closure. Epidemics 2013;5:15763.
 Shil P, Gurav YK, Chadha MS, Mishra A. Transmission dynamics of novel influenza A/H1N1 2009 outbreak in a residential school in India. Curr Sci 2011;100:1177.
 Tan X, Yuan L, Zhou J, Zheng Y, Yang F. Modeling the initial transmission dynamics of influenza A H1N1 in Guangdong province, China. Int J Infect Dis 2013;17:e47984.
 Tang S, Xiao Y, Yang Y, Zhou Y, Wu J, Ma Z, et al. Communitybased measures for mitigating the 2009 H1N1 pandemic in China. PLoS One 2010;5:e10911.
 Tang S, Xiao Y, Yuan L, Cheke RA, Wu J. Campus quarantine (Fengxiao) for curbing emergent infectious diseases: Lessons from mitigating A/H1N1 in Xi'an, China. J Theor Biol 2012;295:4758.
 Truscott J, Fraser C, Hinsley W, Cauchemez S, Donnelly C, Ghani A, et al. Quantifying the transmissibility of human influenza and its seasonal variation in temperate regions. PLoS Curr 2009;1:RRN1125.
 Tsukui S. Casebased surveillance of pandemic (H1N1) 2009 in Maebashi city, Japan. Jpn J Infect Dis 2012;65:1327.
 Van Kerckhove K, Hens N, Edmunds WJ, Eames KT. The impact of illness on social networks: Implications for transmission and control of influenza. Am J Epidemiol 2013;178:165562.
 Ward MP, Maftei D, Apostu C, Suru A. Estimation of the basic reproductive number (R0) for epidemic, highly pathogenic avian influenza subtype H5N1 spread. Epidemiol Infect 2009;137:21926.
 White LF, Wallinga J, Finelli L, Reed C, Riley S, Lipsitch M, et al. Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A/H1N1 pandemic in the USA. Influenza Other Respir Viruses 2009;3:26776.
 Xiao Y, Sun X, Tang S, Wu J. Transmission potential of the novel avian influenza A (H7N9) infection in Mainland China. J Theor Biol 2014;352:15.
 Yang Y, Zhang Y, Fang L, Halloran ME, Ma M, Liang S, et al. Household transmissibility of avian influenza A (H7N9) virus, China, February to May 2013 and October 2013 to March 2014. Euro Surveill 2015;20:21056.
 Yu H, Cauchemez S, Donnelly CA, Zhou L, Feng L, Xiang N, et al. Transmission dynamics, border entry screening, and school holidays during the 2009 influenza A (H1N1) pandemic, China. Emerg Infect Dis 2012;18:75866.
 Zhang S, Yan P, Winchester B, Wang J. Transmissibility of the 1918 pandemic influenza in Montreal and Winnipeg of Canada. Influenza Other Respir Viruses 2010;4:2731.
References   
1.  World Health Organization. Pandemic H1N1 2009. WHO Regional Office for SouthEast Asia; 2009. 
2.  Kilbourne ED. Influenza pandemics of the 20 ^{th} century. Emerg Infect Dis 2006;12:914. 
3.  Askarian M, Danaei M, Vakili V. Knowledge, attitudes, and practices regarding pandemic H1N1 influenza among medical and dental residents and fellowships in Shiraz, Iran. Int J Prev Med 2013;4:396403. 
4.  Al Hajjar S, McIntosh K. The first influenza pandemic of the 21 ^{st} century. Ann Saudi Med 2010;30:110. 
5.  RahimiRad S, Alizadeh A, Alizadeh E, Hosseini SM. The avian influenza H9N2 at avianhuman interface: A possible risk for the future pandemics. J Res Med Sci 2016;21:51. [Full text] 
6.  Mubareka S, Lowen AC, Steel J, Coates AL, GarcíaSastre A, Palese P. Transmission of influenza virus via aerosols and fomites in the Guinea pig model. J Infect Dis 2009;199:85865. 
7.  Bridges CB, Harper SA, Fukuda K, Uyeki TM, Cox NJ, Singleton JA. Prevention and control of influenza. Recommendations of the advisory committee on immunization practices (ACIP). MMWR Recomm Rep 2003;52:134. 
8.  PlansRubió P. Prevention and control of influenza in persons with chronic obstructive pulmonary disease. Int J Chron Obstruct Pulmon Dis 2007;2:4153. 
9.  Viboud C, Alonso WJ, Simonsen L. Influenza in tropical regions. PLoS Med 2006;3:e89. 
10.  SaundersHastings PR, Krewski D. Reviewing the history of pandemic influenza: Understanding patterns of emergence and transmission. Pathogens 2016;5. pii: E66. 
11.  Shrestha SS, Swerdlow DL, Borse RH, Prabhu VS, Finelli L, Atkins CY, et al. Estimating the burden of 2009 pandemic influenza A (H1N1) in the United States (April 2009April 2010). Clin Infect Dis 2011;52 Suppl 1:S7582. 
12.  Wu KM, Riley S. Estimation of the basic reproductive number and mean serial interval of a novel pathogen in a small, wellobserved discrete population. PLoS One 2016;11:e0148061. 
13.  Becker NG, Bahrampour A. Preventing epidemics with agespecific vaccination schedules. Math Biosci 1997;142:6377. 
14.  Woolhouse ME, Hasibeder G, Chandiwana SK. On estimating the basic reproduction number for schistosoma haematobium. Trop Med Int Health 1996;1:45663. 
15.  Hsieh YH, Fisman DN, Wu J. On epidemic modeling in real time: An application to the 2009 novel A (H1N1) influenza outbreak in Canada. BMC Res Notes 2010;3:283. 
16.  
17.  
18.  Svensson A. A note on generation times in epidemic models. Math Biosci 2007;208:30011. 
19.  Dietz K. The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res 1993;2:2341. 
20.  Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci 2007;274:599604. 
21.  Arruda AG, Alkhamis MA, VanderWaal K, Morrison RB, Perez AM. Estimation of timedependent reproduction numbers for porcine reproductive and respiratory syndrome across different regions and production systems of the US. Front Vet Sci 2017;4:46. 
22.  Roberts MG, Nishiura H. Early estimation of the reproduction number in the presence of imported cases: Pandemic influenza H1N12009 in New Zealand. PLoS One 2011;6:e17835. 
23.  Jackson C, Vynnycky E, Mangtani P. Estimates of the transmissibility of the 1968 (Hong Kong) influenza pandemic: Evidence of increased transmissibility between successive waves. Am J Epidemiol 2010;171:46578. 
24.  Haghdoost A, Baneshi MR, Zolala F, Farvahari S, Safizadeh H. Estimation of basic reproductive number of flulike syndrome in a primary school in Iran. Int J Prev Med 2012;3:40813. 
25.  White LF, Pagano M. A likelihoodbased method for realtime estimation of the serial interval and reproductive number of an epidemic. Stat Med 2008;27:29993016. 
26.  McBryde E, Bergeri I, van Gemert C, Rotty J, Headley E, Simpson K, et al. Early transmission characteristics of influenza A (H1N1) v in Australia: Victorian state, 16 May – 3 June 2009. Euro Surveill 2009;14. pii: 19363. 
27.  White LF, Wallinga J, Finelli L, Reed C, Riley S, Lipsitch M, et al. Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A/H1N1 pandemic in the USA. Influenza Other Respir Viruses 2009;3:26776. 
28.  Chowell G, Nishiura H, Bettencourt LM. Comparative estimation of the reproduction number for pandemic influenza from daily case notification data. J R Soc Interface 2007;4:15566. 
29.  Lessler J, Cummings DA, Fishman S, Vora A, Burke DS. Transmissibility of swine flu at Fort Dix, 1976. J R Soc Interface 2007;4:75562. 
30.  Obadia T, Haneef R, Boëlle PY. The R0 package: A toolbox to estimate reproduction numbers for epidemic outbreaks. BMC Med Inform Decis Mak 2012;12:147. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3], [Table 4]
