Age groups and spread of influenza: implications for vaccination strategy
© Hsieh; licensee BioMed Central Ltd. 2010
Received: 21 September 2009
Accepted: 30 April 2010
Published: 30 April 2010
The unpredictable nature of the potentially devastating impact of 2009 pH1N1 influenza pandemic highlights the need for pandemic preparedness planning, where modeling studies could be most useful for simulations of possible future scenarios.
A compartmental model with pre-symptomatic and asymptomatic influenza infections is proposed which incorporates age groups as well as intervention measures such as age-specific vaccination, in order to study spread of influenza in a community.
We derive the basic reproduction number and other effective reproduction numbers under various intervention measures. For illustration, we make use of the Pneumonia and Influenza (P&I) mortality data and vaccination data of the very young (age 0-2) and the very old (age >64) during 2004-2005 Taiwan winter influenza season to fit our model and to compute the relevant reproduction numbers. The reproduction number for this winter flu season is estimated to be slightly above one (~1.0001).
Comparatively large errors in fitting the P&I mortality data of the elderly (>64) were observed shortly after winter school closings in January, which may indicate the impact of younger, more active age groups transmitting influenza to other age groups outside of the school settings; in particular, to the elderly in the households. Pre-symptomatic infections seemed to have little effect on the model fit, while asymptomatic infection by asymptomatic infectives has a more pronounced impact on the model fit for the elderly mortality, perhaps indicating a larger role in disease transmission by asymptomatic infection. Simulations indicate that the impact of vaccination on the disease incidence might not be fully revealed in the change (or the lack thereof) in the effective reproduction number with interventions, but could still be substantial. The estimated per contact transmission probability for susceptible elderly is significantly higher than that of any other age group, perhaps highlighting the vulnerability of the elderly due to close contacts with their caretakers from other age groups. The relative impact of targeting the very young and the very old for vaccination was weakened by their relative inactivity, thus giving evidence of the lack of impact of vaccinating these two groups on the overall transmissibility of the disease in the community. This further underscores the need for morbidity-based strategy to prevent elderly mortality.
In the spring of 2009, the novel H1N1 influenza virus first emerged in Mexico and later spread widely throughout the world within just a few months. The World Health Organization (WHO) announced on June 11 the start of 2009 influenza pandemic , and further issued an advisory on August 28 for countries in the northern hemisphere to prepare for a second wave of pandemic spread in the coming fall/winter . As of November 8, more than 206 countries and overseas territories or communities worldwide have reported laboratory confirmed cases of the pandemic pH1N1 virus, including over 6250 deaths .
To lessen the severity of this pandemic, developing an effective flu vaccine and a global vaccination strategy is considered to be among the most important medical interventions . However, to have the greatest impact, pandemic vaccines need to be available quickly and in large quantities, and to be delivered to the population optimally. Moreover, vaccines against a novel pandemic strain might take up to six months to manufacture and deliver, even in developed countries . Given the potential threat of drug-resistance resulting from widespread use of antiviral treatment against pandemic flu, vaccine appears to be our primary weapon to prevent and to mitigate a pandemic. However, in addition to the need to consider the logistics of implementing large-scale vaccination, distinctly different age-specific mortalities had also been observed during some past flu pandemics (e.g., in 1918 ), which require different priorities when large-scale vaccination is to be implemented.
Moreover, vaccine for influenza is known to have different efficacy (i.e., reduction in the number of infectives) and effectiveness (i.e., reduction in symptomatic case number) for different age groups, see e.g., [7–9]. Setting priority for vaccination by targeting age groups most vulnerable (the elderly, infants, etc.) to prevent mortalities is commonly employed in most countries. However, when vaccinating those at greatest risk of mortality becomes impractical (if, e.g., medical care is relatively inaccessible) or inefficient (if, e.g., immune response is deficient), targeting those most likely to expose them to infection might be more preferable . In this way, the very young and the old might be better protected by vaccinating those who are most likely to be in contact with them (thereby reducing their risk of exposure), rather than by being vaccinated. Comparison of influenza mortality among elderly Japanese during time periods when schoolchildren were and were not vaccinated suggests that the infected children pose a risk to others , including the elderly. Moreover, several past US experiences (as summarized in ) also are consistent with this conclusion. Nonetheless, influenza policymakers have typically advocated protecting those individuals of ages 6-24 months and >65 years directly.
Bansal et al.  recently carried out a comparative analysis of two classes of suggested vaccination strategies, namely, the mortality-based strategies that target the high-risk populations and the morbidity-based strategies that target the high-prevalence populations, by applying the methods of contact network epidemiology to a model of disease transmission in a large urban population. Using a range of mortality rates reported previously for past influenza epidemics and pandemics, they concluded that the optimal strategy depends critically on the viral transmission level (or reproduction number) of the virus. That is, the morbidity-based strategies outperform the mortality-based strategies for moderately transmissible strains, while the reverse is true for highly transmissible strains. However, they also cautioned that when information pertaining to viral transmission rate of a particular disease and the frequency of new introductions into the community prior to an outbreak is unreliable or not available, a mortality-based vaccination priority is recommended. This further demonstrates the importance of targeting and, moreover, the uncertainty surrounding this issue.
To further the uncertainties regarding influenza pandemic preparedness planning, it is widely believed that asymptomatic cases (i.e., individuals who had been infected but showed little or no symptoms) and asymptomatic infection of influenza (i.e., infection caused by an asymptomatic case) do indeed occur regularly (e.g., [14–17]).
Model with only asymptomatic infections, by either asymptomatic or subclinical infectives, during their infectivity period had been recently studied in . In this current study, we will consider a traditional compartmental model which incorporates both pre-symptomatic and asymptomatic infections, in order to explore the role which they could play in the overall spread of disease, if any. Moreover, the age-group structure of the model, by dividing the population into seven groups of the very young, preschool children, younger and older schoolchildren, young adults, adults, and the elderly, allows us to study targeted public health policies (e.g., immunization) aimed at different age groups. Our model also allows for inclusion of immunity and other age-dependent intervention measures such as quarantine and voluntary home withdrawal (see e.g., [15, 16]). A full model will be proposed to take into account of the above-mentioned factors that may be important in determining the best vaccine strategy.
S i (t): number of susceptible individuals of the ith age group at time t;
V i (t): number of vaccinated individuals of the ith age group at time t;
E i (t): number of exposed (infected) individuals of the ith age group at time t;
Q i (t): number of quarantined infected individuals of the ith age group at time t;
I i (t): number of infective individuals of the ith age group at time t;
I i v (t): number of vaccinated infective individuals of the ith age group at time t;
A i (t): number of asymptomatic (subclinical) infective individuals of the ith age group at time t;
H i (t): number of hospitalized (treated) individuals of the ith age group at time t;
Z i (t): number of recovered and immune individuals of the ith age group at time t;
D i (t): cumulative number of influenza deaths of the ith age group at time t;
cij: contact rate of an individual of ith group with an individual of jth group;
β ij : per contact transmission probability of a susceptible individual of ith group by an infective of jth group.
π i : age-specific vaccine efficacy for age group i.
λ i (t) and : disease incidence rates for the susceptible and vaccinated individuals of age group i. See [Additional file 1] for detailed formulae.
mean vaccine waning rate
mean incubation period of exposed individuals
mean time to onset of those who had been quarantined
mean recovery rate of untreated infectives
mean recovery rate of vaccinated infectives with ρ V ≥ ρ
mean recovery rate of asymptomatic infectives
mean immune waning rate
migration rate of the population
migration rate of symptomatic infectives
immigration rate of the population
immigration rate of symptomatic infectives
quarantine rate of unvaccinated exposed individuals
quarantine rate of vaccinated exposed individuals
home withdrawal rate of untreated symptomatic infectives
home withdrawal rate of all "well" individuals
reduction in infectivity of unvaccinated pre-symptomatic infectives
reduction in infectivity of asymptomatic infectives
reduction factor in contact due to hospital isolation
reduction in infectivity of vaccinated infectives
time-dependent vaccination rate ν i (t)
υ 1 (t)
υ 7 (t)
mortality rate of untreated infectives δ i
6 × 10-4
5 × 10-4
5 × 10-5
5 × 10-5
5 × 10-5
1 × 10-4
6 × 10-4
5 × 10-4
5 × 10-5
5 × 10-5
5 × 10-5
1 × 10-4
5 × 10-4
5 × 10-4
5 × 10-4
mean recovery rate of the hospitalized infectives ρ 1
hospitalization rate of unvaccinated symptomatic infectives h i
per contact transmission probability of infectees in i th group β ij ≡ β i
per contact transmission probability of infectors in i th group β ij ≡ β j
(1 )Exposed individuals are infective during the incubation period. It is commonly known (e.g., ) that the pre-symptomatic (exposed) individuals cannot transmit the disease in the noninfectious latent period, during which the viral titres gradually increase to detectable and transmissible levels when they became infective for only a short period (0.25 days in ) before the onset of symptoms. To avoid adding an extra compartment to account for the (short) period of infectivity after the end of the latent period and before the end of incubation period (see e.g., ), we assume that the individuals are infective during the incubation period with the infectivity averaged out over the whole incubation period. This simplification is reasonable since the incubation period for influenza is typically very short, e.g., 1.48 days in .
(2) Influenza vaccine is efficacious in preventing influenza infection and effective against influenza-like illness, albeit at different levels of efficacy and effectiveness for different age groups [8, 9]. Moreover, the vaccinated individuals are less infectious, once they become infected, when compared to those who had not been vaccinated.
(4) A fraction of the infectives has no symptoms or only subclinical symptoms, and is classified as asymptomatic infectives with reduced infectivity .
(5) A hospitalized person is removed from isolation either by death or discharged due to recovery from illness.
(6) Homogeneous mixing within subpopulations is assumed.
(7) Negligible births and deaths (excluding disease deaths) during the course of the disease outbreak are assumed.
The two blue arrows denotes within-group infections (R 11 and R 22 ) and the black arrows denote inter-group infection cycle (R 12 R 21 ). The term R 11 R 22 subtracted in (3) accounts for the redundancy that resulted when adding R 11 and R 22 . Similar results for the basic reproduction number of a multi-group model were also obtained in [26, 27].
Moreover, we have the following effective reproduction numbers due to interventions:
i, j = 1, 2, ..., n.
Detailed derivations of the reproduction numbers are also given in [Additional file 1]. Model fit using 2005-2006 Taiwan winter seasonal influenza data and age-specific vaccination data as it was implemented during that flu season, as well as simulation studies of hypothetical scenarios, will be carried out.
Simulations with Taiwan Seasonal Influenza
With the issue of morbidity-based vs. mortality-based vaccine strategy for pandemic influenza still open to debate, the Taiwan Centers for Disease Control (TCDC) launched a new program of free flu vaccination for 1st and 2nd grades elementary school students (age 6-7) prior to the 2007-2008 winter flu season which was expanded further to include grades 1-4 in the fall of 2008. The aim of this vaccine program is hopefully to lower the seasonal influenza incidence across all age groups of the population. In anticipation of future investigation on public health impact of this program, we carry out model simulations by dividing the Taiwan population into 7 age groups (see Table 2), and by making use of the weekly Taiwan influenza vaccination data for young children of age 2 or less (age group 1, with free flu vaccination since 2004) and the elderly of 65 or older (age group 7, with free flu vaccination since 2001) during the 2004-2005 winter flu season. The average vaccine coverages of age groups 1 and 7 during the 2004-2005 winter season in Taiwan are 63.4% and 58.2%, respectively.
There were no other nonpharmaceutical intervention measures during this winter flu season. That is, all parameters pertaining to quarantine and home withdrawal in the model are set to be 0 in Table 1. Moreover, for the sake of simplicity, we assume no noticeable level of migration and no waning of immunity during the flu season. The rest of the parameter values used are given in Tables 1-2. We also assume a conservative 20% pre-epidemic immunity in our simulation based on a recent sero-epidemiological survey conducted during 2005-2006 winter flu season in Taiwan .
For the contact rates between different age groups, we make use of the age-specific contact matrix obtained by Wallinga et al.  for Utrecht, the Netherlands, 1986. We adjust for the discrepancy in the population age distributions of Taiwan in 2005 and the Netherlands in 1986, by considering the ratios of demographic age structures of Netherlands in 1987 (Appendix Table 1 in ) and of Taiwan in 2005 [Additional file 1: Table A1]. The resulting contact matrix, of the average daily number of contacts for each individual in a certain age group with individuals in another age group, is given in [Additional file 1: Table A2].
Fitting with Seasonal Influenza Data
The age-dependent hospitalization rates and per contact transmission probability in the last four rows (in bold) in Table 2 were obtained by least-squared curve-fitting with the 2004-2005 Taiwan winter P&I (Pneumonia and Influenza) mortality data from October 9, 2004 to March 5, 2005 for age groups 6 (age 22-64) and 7 (>64) using MATLAB software. ν 7(τ) and ν 1(τ) are piecewise linear (by week) approximations of the respective weekly vaccination data for elderly (>64) and young children in Taiwan during this time period. To simply our data fitting, we first fitted the data by assuming the per contact transmission probability of the infectees in ith group are averaged, i.e., β ij ≡ β i (see next to the last row in bold in Table 2).
Discussion and Conclusion
By fitting only the hospitalization rates and the per contact transmission probabilities to the 2004-2005 Taiwan winter season P&I mortality for age groups 6 (age 22-64) and 7 (>64), we are able to obtain satisfactory model fit (see Fig. 3), although we note that fitting P&I deaths might conceivably lead to an overestimate of hospitalizations rates. However, simulation studies showed that the model fit is less sensitive to the hospitalization rates than to the transmission probabilities, and intuitively, most sensitive to changes in these rates for the elderly group.
It is interesting to note that comparatively significant errors in both theoretical curves occur in late January roughly after day 110. We note that day 105 was January 22, 2005, when all schools in Taiwan closed for the winter vacation which lasted until after the traditional lunar New Year holiday, in mid-February of that year. The school closure, and subsequent shutdown of all non-essential venues during the week-long New Year holiday, surely had a significant impact on the contact rates which was not reflected in our constant contact matrix that implicitly assumes that inter-age group contact patterns remain the same during the whole season. More precisely, the P&I deaths for the elderly exhibits a slight increase for about 3 weeks after the closing of schools (shortly after the New Year) when compared to the model predicted values, while the P&I deaths for the adults of age 22-64 dropped substantially below the theoretical curve. These results indicate that the school closure and the subsequent New Year holiday led to more frequent contacts by the elderly in the households with family members who spent more time at home during the holidays At the same time, there were less contacts for the adults at the workplace (especially for those working in the educational facilities who had longer holidays) during this time period. Since most individuals in the elder group (and the very young) are not directly affected by the school closings and the holidays, one may speculate that the difference in the actual mortality and the theoretical mortality of the elderly, as averaged over the whole time period, is due at least in parts to the impact of interaction between the elderly and younger children with the school-age children and adults, when the activity levels (in terms of frequency as well as whom to make contact with) of the latter groups were significantly changed by the holidays.
Furthermore, the elderly had a higher per contact transmission probability (Table 2) but lower contact frequency [Additional file 1: Table A2]. This study showed that, with the combination of these two factors, morbidity-based vaccination strategy still could be effective for the prevention of elderly mortality. Although a recent study to quantify the effect of school closures during the 2008 winter influenza season in Hong Kong did not find the school closures as having a substantial effect on the community transmission , our results indicate the need for further studies on this topic. We also note that, normally, one would like to use averaged excess winter season P&I mortality data over other seasons for the data fitting of seasonal influenza. However, Taiwan often experiences summer influenza epidemics which would offset any attempt to obtain a reasonable "excess" P&I mortality for the winter season.
The Reproduction Number
Using Equations (2-3) we obtain R0 = 1.0001, just above unity. However, we note that it is more appropriately the reproduction number R of the winter flu epidemic, given the effect of pre-epidemic immunity that must exist. Chowell et al.  used several weekly seasonal flu mortality data, derived from P&I excess deaths and influenza-specific deaths from US, France, and Australia during 1972-1997 (1972-2002 for US), to estimate the (mean) reproduction number Rp over 3 decades of seasonal flu. They found that the mean of Rp to be around 1.3, with year-to-year variability of 0.9-2.1. Our estimate is lower, but within their range.
Furthermore, in our simulations we had assumed vaccine efficacy (proportion of infection prevented) of 40% for elderly and 70% for young children, to be in line with current literature [8, 9]. However, efficacy depends very much on matching of the circulating strains with vaccine strains each year, where mismatch often causes low efficacy and might have affected our resulting data fit. We also assumed vaccine effectiveness (proportion of reduction in symptomatic cases) to be 60% for elderly and 40% for young children (also see [8, 9]).
The per contact transmission probability for the susceptible elderly (>64) infectees (see next to the last row in Table 2) is estimated to be 0.155, more than three-fold of any other age group (with the young children of 0-2 being the next highest), which may be due to the common need for very close contact while the elderly (or younger children of age <3) are being cared for, typically by individuals from adults of ages 22-64 in age group 6, even though the frequency of contact might be less than that of with other age groups. This further highlights the high vulnerability of the elderly (or the younger children) to exposure from other age groups, and demonstrates the need for morbidity-based strategy to prevent the elderly (or younger children) influenza mortality. On the other hand, the estimated per contact transmission probabilities for the younger (<3) and the elderly (>64) infectors are also both higher than those of the other groups (see the last row in Table 2), but not significantly so except when compared with those between ages 8-14. The higher transmission probabilities of young children and elderly infectives could reflect, again, the fact that contacts with these individuals are usually of a more intimate nature, although these probabilities are not as drastically different as their vulnerability to be infected. The less likelihood of schoolchildren of 8-14 to infect others might also be attributable to signs of their less intimate contacts with others as they grow into adolescence.
We know from Equations (2) and (7). We also know that κ ∈ [0, 1] and ν i π i for all i, j. Subsequently, is a nonincreasing function of ν i π i for all i, j. However, since R V = (-1) n+1 det R V + 1, where is the matrix with its ijth element, the effective reproduction number with vaccination only, R V , does not necessarily decrease as the effective vaccination rate ν i π i increases. In other words, vaccination is not always beneficial in reducing incidence and the design of an effective vaccination program in multi-group model is highly nontrivial, unlike in simple epidemic models where there is a simple formula for the critical vaccination coverage level necessary for eradication (pc) . In fact, it has been shown mathematically by Hadeler and Castillo-Chavez , using a model with a core group, that partially effective vaccination program may actually increase the total number of cases. Explicit quantification of an optimal vaccination policy () in a multi-group population scenario remains to be a challenge for mathematical modelers and is beyond the scope of this work. Moreover, the model fit was carried out with vaccination (of the elderly and small children of age 2 or less) as the only intervention since vaccination was the only intervention that was implemented during the flu seasonal for which the fitted P&I data was collected. Further simulations with quarantine and home withdrawal also can be easily carried out.
A related issue is that of pre-epidemic immunity. Further sensitivity analysis using pre-immunity in the range of 10%-30% has shown that the results of curve-fittings are not sensitive to small changes in the pre-epidemic immunity. One would expect that pre-epidemic immunity does impact the outbreak, but nonetheless it is not reflected, at least not in the reproduction numbers.
Finally, as the WHO Strategic Advisory Group of Experts (SAGE) recently made its recommendations for priorities in vaccination for the H1N1 pandemic in terms of the social groups (e.g., healthcare workers those with chronic medical conditions), health groups (pregnant women), and age groups , the model also can be used to divide population into social/health groups. For example, to study vaccine policy for the elderly, we could divide the elderly into those living in households and those living in old age homes, since they mix differently in these two distinct settings. The model is also useful for simulations of the cost-effectiveness of vaccine and other intervention measures, such as prophylaxis treatment, perhaps in future work.
YHH was supported by Taiwan CDC Pandemic Influenza Vaccine Research and Development Program (DOH95-DC-1407), National Science Council of Taiwan (NSC 97-2118-M-039-004), and China Medical University (CMU 97 323). The author is grateful to Leng Leng Lim and Yun-Shih Wang for their help with the numerical simulations, and to Jen-Hsiang Chuang for providing Taiwan age-specific P&I mortality data. The author is grateful to two reviewers for their constructive and insightful comments which significantly improved this paper.
- World Health Organization: World now at the start of 2009 influenza pandemic, WHO, Director-General statement to the press. Accessed September 18, 2009, [http://www.who.int/mediacentre/news/statements/2009/h1n1_pandemic_phase6_20090611/en/index.html]
- World Health Organization: Preparing for the second wave: lessons from current outbreaks. Pandemic (H1N1) 2009 briefing note 9. accessed August 31, 2009, [http://www.who.int/csr/disease/swineflu/notes/h1n1_second_wave_20090828/en/index.html]
- World Health Organization: Pandemic (H1N1). accessed November 17, 2009, [http://www.who.int/csr/don/2009_11_13/en/index.html]
- World Health Organization: Safety of pandemic vaccines. Accessed September 18, 2009, [http://www.who.int/csr/disease/swineflu/notes/h1n1_safety_vaccines_20090805/en/index.html]
- Yamada T, Dautry A, Walport M: Ready for avian flu?. Nature. 2008, 454: 162-10.1038/454162a.View ArticlePubMedGoogle Scholar
- Taubenberger JK, Morens DM: 1918 influenza: the mother of all pandemics. Emerg Infect Dis. 2006, 12: 15-22.View ArticlePubMedPubMed CentralGoogle Scholar
- Nicholson KG, Wood J, Zambon M: Influenza. Lancet. 2003, 362: 1733-1745. 10.1016/S0140-6736(03)14854-4.View ArticlePubMedGoogle Scholar
- Jefferson T, Smith S, Demicheli V, Harnden A, Rivetti A, Di Pietrantonj C: Assessment of the efficacy and effectiveness of influenza vaccines in healthy children: systematic review. Lancet. 2005, 365: 773-780.View ArticlePubMedGoogle Scholar
- Jefferson T, Rivetti D, Rivetti A, Rudin M, Di Pietrantonj C, Demicheli V: Efficacy and effectiveness of influenza vaccines in elderly people: a systematic review. Lancet. 2005, 366: 1165-1174. 10.1016/S0140-6736(05)67339-4.View ArticlePubMedGoogle Scholar
- Lloyd-Smith JO, Schreiber SJ, Kopp PE, Getz WM: Superspreading and the effect of individual variation on disease emergence. Nature. 2005, 438: 355-359. 10.1038/nature04153.View ArticlePubMedGoogle Scholar
- Reichert TA, Sugaya N, Fedson DS, Glezen WP, Simonsen L, Tashiro T: The Japanese experience with vaccinating schoolchildren against influenza. N Engl J Med. 2001, 44: 889-896. 10.1056/NEJM200103223441204.View ArticleGoogle Scholar
- Glezen WP: Emerging infections: pandemic influenza. Epidemiol Rev. 1996, 18: 64-76.View ArticlePubMedGoogle Scholar
- Bansal S, Pourbohloul B, Meyers LA: A Comparative Analysis of Influenza Vaccination Programs. PLoS Med. 2006, 3 (10): e387-10.1371/journal.pmed.0030387.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu SB, Hsieh YH: On the Role of Asymptomatic Infection in Transmission Dynamics of Infectious Diseases. Bull Math Biology. 2008, 70: 134-155. 10.1007/s11538-007-9245-6.View ArticleGoogle Scholar
- Longini IM, Halloran ME, Nizam A, Yang Y: Containing pandemic influenza with antiviral agents. Am J Epidemiol. 2004, 159: 623-633. 10.1093/aje/kwh092.View ArticlePubMedGoogle Scholar
- Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings DA, et al: Containing pandemic influenza at the source. Science. 2005, 309: 1083-1087. 10.1126/science.1115717.View ArticlePubMedGoogle Scholar
- Germann TC, Kadau K, Longini IM, Macken CA: Mitigation strategies for pandemic influenza in the United States. PNAS. 2006, 103 (15): 5935-5940. 10.1073/pnas.0601266103.View ArticlePubMedPubMed CentralGoogle Scholar
- Baccam P, Beauchemin C, Macken CA, Hayden FG, Perelson AS: Kinetics of influenza A virus infection in humans. J Virol. 2006, 80: 7590-7599. 10.1128/JVI.01623-05.View ArticlePubMedPubMed CentralGoogle Scholar
- Alexander ME, Bowman CS, Feng Z, Gardam M, Moghadas SM, Röst G, Wu J, Yan P: Emergence of drug resistance: implications for antiviral control of pandemic influenza. Proc Biol Sci. 2007, 274 (1619): 1675-1684. 10.1098/rspb.2007.0422.View ArticlePubMedPubMed CentralGoogle Scholar
- Anderson R, May R: Dynamics and Control. 1991, Oxford: Oxford University PressGoogle Scholar
- Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, Lamsirithaworn S, Burke DS: Strategies for containing an emerging influenza pandemic in southeast Asia. Nature. 2005, 437: 209-214. 10.1038/nature04017.View ArticlePubMedGoogle Scholar
- Hsieh YH, King CC, Ho MS, Chen CWS, Lee JY, Liu FC, Wu YC, Wu JSJ: Quarantine for SARS, Taiwan. Emerg Infect Dis. 2005, 11 (2): 278-282.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsieh YH, King CC, Chen CWS, Ho MS, Hsu SB, Wu YC: Impact of Quarantine on the 2003 SARS Outbreak: a retrospective modeling study. J Theoretical Biology. 2007, 244: 729-736. 10.1016/j.jtbi.2006.09.015.View ArticleGoogle Scholar
- Diekmann O, Heesterbeek JAP, Metz JAJ: On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990, 28: 365-382. 10.1007/BF00178324.View ArticlePubMedGoogle Scholar
- Driessche van den P, Watmough J: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002, 180: 29-48. 10.1016/S0025-5564(02)00108-6.View ArticlePubMedGoogle Scholar
- Hsieh YH, Chen CH: Modeling the social dynamics of the sex industry in Thailand: Its implications for spread of HIV. Bull Math Biol. 2004, 66 (1): 143-166. 10.1016/j.bulm.2003.08.004.View ArticlePubMedGoogle Scholar
- Hsieh YH, Wang YS: Basic reproduction number HIV model incorporating commercial sex and behavior change. Bull Math Biol. 2006, 68: 551-575. 10.1007/s11538-005-9050-z.View ArticlePubMedGoogle Scholar
- Lin CY: Serological surveillance of human influenza viruses among elementary schoolchildren in Taiwan during the 2005-2006 flu season. Master thesis. 2007, National Taiwan University, Graduate Institute of EpidemiologyGoogle Scholar
- Wallinga J, Teunis P, Kretzschmar M: Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. Am J Epidemiol. 2006, 164: 936-944. 10.1093/aje/kwj317.View ArticlePubMedGoogle Scholar
- Cowling BJ, Lau EHY, Lam CLH, Cheng CKY, Kovar J, Chan KH, Peiris JSM, Leung GM: Effects of school closures, 2008 winter influenza season, Hong Kong. Emerg Infect Dis. 2008, 14 (10): 1660-1662. 10.3201/eid1410.080646.View ArticlePubMedPubMed CentralGoogle Scholar
- Chowell G, Miller MA, Viboud C: Seasonal influenza in the United States, France, and Australia: transmission and prospects for control. Epidemiol Infect. 2008, 136 (6): 852-864. 10.1017/S0950268807009144.View ArticlePubMedGoogle Scholar
- Hadeler KP, Castillo-Chavez C: A core group model for disease transmission. Math Biosci. 1995, 128: 41-55. 10.1016/0025-5564(94)00066-9.View ArticlePubMedGoogle Scholar
- Hill AN, Longini IM: The critical vaccination fraction for heterogeneous epidemic models. Math Biosci. 2003, 181 (1): 85-106. 10.1016/S0025-5564(02)00129-3.View ArticlePubMedGoogle Scholar
- World Health Organization: WHO recommendations on pandemic (H1N1) 2009 vaccines. Accessed September 18, 2009, [http://www.who.int/csr/disease/swineflu/notes/h1n1_vaccine_20090713/en/index.html]
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/106/prepub