Evaluation of Targeted Influenza Vaccination Strategies via Population Modeling

doi:10.1371/journal.pone.0012777

PLoS One. 2010; 5(9): e12777.

Published online 2010 September 17. doi: 10.1371/journal.pone.0012777

PMCID: PMC2941445

Copyright This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.

Evaluation of Targeted Influenza Vaccination Strategies via Population Modeling

John Glasser,¹^* Denis Taneri,¹ Zhilan Feng,² Jen-Hsiang Chuang,³ Peet Tüll,⁴ William Thompson,¹ Mary Mason McCauley,¹ and James Alexander¹

¹National Center for Immunization and Respiratory Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America

²Department of Mathematics, Purdue University, West Lafayette, Indiana, United States of America

³Epidemic Intelligence Center, Centers for Disease Control, Taipei, Taiwan

⁴Scientific Advice Unit, European Centre for Disease Prevention and Control, Stockholm, Sweden

Andy T. Y. Lau, Editor

University of Minnesota, United States of America

* E-mail: jglasser/at/cdc.gov

Conceived and designed the experiments: JWG DT ZF WT. Performed the experiments: JWG. Analyzed the data: JWG JHC PT WT MMM JA. Contributed reagents/materials/analysis tools: DT ZF JHC WT. Wrote the paper: JWG. Revised the manuscript: JHC PT WT MMM JA.

Received June 2, 2010; Accepted August 8, 2010.

This article has been cited by other articles in PMC.

Abstract

Background

Because they can generate comparable predictions, mathematical models are ideal tools for evaluating alternative drug or vaccine allocation strategies. To remain credible, however, results must be consistent. Authors of a recent assessment of possible influenza vaccination strategies conclude that older children, adolescents, and young adults are the optimal targets, no matter the objective, and argue for vaccinating them. Authors of two earlier studies concluded, respectively, that optimal targets depend on objectives and cautioned against changing policy. Which should we believe?

Methods and Findings

In matrices whose elements are contacts between persons by age, the main diagonal always predominates, reflecting contacts between contemporaries. Indirect effects (e.g., impacts of vaccinating one group on morbidity or mortality in others) result from off-diagonal elements. Mixing matrices based on periods in proximity with others have greater sub- and super-diagonals, reflecting contacts between parents and children, and other off-diagonal elements (reflecting, e.g., age-independent contacts among co-workers), than those based on face-to-face conversations. To assess the impact of targeted vaccination, we used a time-usage study's mixing matrix and allowed vaccine efficacy to vary with age. And we derived mortality rates either by dividing observed deaths attributed to pneumonia and influenza by average annual cases from a demographically-realistic SEIRS model or by multiplying those rates by ratios of (versus adding to them differences between) pandemic and pre-pandemic mortalities.

Conclusions

In our simulations, vaccinating older children, adolescents, and young adults averts the most cases, but vaccinating either younger children and older adults or young adults averts the most deaths, depending on the age distribution of mortality. These results are consistent with those of the earlier studies.

Introduction

Seasonal influenza causes an estimated 200,000 hospitalizations and 36,000 deaths on average in the United States, most among the elderly [1]. If a 1918-like pandemic occurred today, 10 million hospitalizations and 1.9 million deaths – many among younger adults – are expected [2]. Vaccination affords the best protection, especially for those at risk of pneumonia and other life-threatening complications [1].

Development and production of influenza vaccines is challenging. In the northern hemisphere, the World Health Organization (WHO) collects relevant information every February for review by experts. Based on which viruses they believe will most likely be circulating, the experts select 3 strains for inclusion in the upcoming season's vaccine. Almost every year, at least one vaccine constituent is replaced, because viral strains drift; i.e., undergo constant genetic change. Even small changes can result in novel strains, and mismatch with circulating strains can reduce vaccine effectiveness, as occurred during the 2007–08 influenza season [3]. Once the experts have identified the strains likely to circulate next season, a vaccine must be manufactured in a slow process that has changed little since its invention. Testing, approval, and distribution also take several months. Problems encountered during production, such as inability to grow sufficient quantities of a viral strain, may cause vaccine shortages or delays in distribution. Such problems have affected vaccine availability in recent influenza seasons in the United States [4].

During influenza pandemics, these challenges are compounded. Pandemic strains may emerge when antigenic shifts – major changes in the genetic makeup of a virus – occur in influenza A, creating new viral subtypes against which populations have little or no immunity [5]. Even when effective vaccines are created, acute shortages are possible, especially in areas with limited production capacity that also have little advance warning, making it difficult or impossible to obtain sufficient vaccine in time to protect at-risk populations. During the recent pandemic, even in wealthy countries that developed and produced an H1N1 vaccine as soon as possible, vaccine supplies were inadequate to accommodate all who sought timely vaccination. The prospect of a shortage motivated health authorities in influenza vaccine-producing countries to devise strategies for ensuring that people who were most likely to suffer complications of influenza were vaccinated first. In the United States, the CDC's Advisory Committee on Immunization Practices (ACIP) determined that pregnant women, caregivers of young infants, health care workers, and people too young to have antibodies to H1N1 had first priority. Next were those most vulnerable to complications of influenza, generally the elderly [6].

In such circumstances, other strategies for using scarce influenza vaccine efficiently also warrant consideration. Among such strategies is indirect protection; that is, immunizing those who might infect vulnerable people. One group whose vaccination might achieve the benefits of indirect protection is schoolchildren. The merits of vaccinating schoolchildren against influenza, partly to protect others, such as the elderly, have been argued from community-intervention trials [7], natural experiments [8], and individual-based models [9]. While trials generally are better controlled than natural experiments, they are relatively expensive and time-consuming. Moreover, only models allow examination of alternative vaccination strategies in exactly the same setting. Models should be evaluated against historical observations to check their predictive ability, but identifying and remedying deficiencies of individual-based models can be prohibitively difficult. Population models are simple enough for evaluation before use to inform public policy making. Analytical results, such as the optimal targets for interventions against infectious diseases, also can be derived.

To identify vaccine allocation strategies with the greatest potential to reduce influenza morbidity and mortality, we studied an age-structured population model whose infection rates we estimated from observed proportions infected [10] and interpersonal contacts weighted by duration [11]. Our model's disease-induced mortality rates were either quotients of deaths attributed to pneumonia or influenza [12] and populations at risk or products of those rates and ratios of 1918 and average 1913–17 mortalities [13]. We refer to the latter as contemporary 1918-like mortality.

Methods

We adapted a demographically-realistic version of a classic population model [14] with 4 disease or immune states: susceptible; infected, but not yet infectious (exposed); infectious; recovered and immune (removed). We added vaccination with age-specific efficacy, based on the work of Govaert et al., who conducted the only randomized, double-blind, placebo-controlled trial of vaccination against morbidity [15], and loss of immunity to circulating strains, via antigenic shifting and drifting [16]. For lack of the requisite information, we ignored transient protection via maternal antibodies, despite how important this protection may be, given seasonal influenza complications among infants aged <6 months [17]. For simplicity, we also ignored immigration and emigration. File S1 and Table S1 describe the system of equations and parameter values.

Age-structured models require multiple infection rates, to which Anderson and May [18] referred collectively as “who-acquires-infection-from-whom.” We derived ours from age-specific proportions of household members infected during the 1957 influenza pandemic [10], commonly called “attack rates,” and from interpersonal contacts weighted by duration in Portland, Oregon [11]. Briefly, the risks of infection, An external file that holds a picture, illustration, etc.
Object name is pone.0012777.e001.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0012777.e001.jpg

, where a_i are average numbers of contacts per person per day; β_i are probabilities of infection upon contact with infectious persons; c_ij are proportions of contacts that members of group i have with those of group j; and y_j

I_j/N_j are probabilities that randomly encountered members of group j are infectious. Using a logistic regression model fitted to the y_i reported by Chin et al. [10] and the duration-weighted contacts of Del Valle et al. [11], together with the relationships An external file that holds a picture, illustration, etc.
Object name is pone.0012777.e002.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0012777.e002.jpg

and

An external file that holds a picture, illustration, etc.
Object name is pone.0012777.e003.jpg

, we estimated λ_i and then β_i.

Figure 7 of Glasser et al. [unpublished manuscript] illustrates β_i obtained using these “attack rates” and the C_ij from several recently published studies of face-to-face conversations or periods in proximity with others during which respiratory diseases might be transmitted. Del Valle et al. [11] not only weighted contacts by duration, but the off-diagonal elements of their contact matrix are relatively large, increasing possible indirect effects (e.g., impact of vaccination on morbidity or mortality in groups not targeted), and they kindly shared their observations. Thus, we could calculate empirical rates of effective contact between members of any age groups simply by averaging (Figure 1).

Figure 1

Effective contact or infection rates derived from attack “rates” during the 1957 pandemic [10] and daily contacts weighted by duration [11].

We calculated age-specific disease-induced mortalities as quotients of 2005 deaths attributed to pneumonia or influenza (Table S1) and simulated infections during an average year. We obtained contemporary 1918-like mortalities by fitting logistic regression models to published 1918 and average 1913–17 rates [13], calculating age-specific ratios for the groups we modeled, <1, 1–4, 5–9, …, 80–84, 85+ years (Figure 2a), and multiplying them by the 2005 rates. Our estimates of mortality conditional on influenza during 2005 and a hypothetical contemporary 1918-like pandemic are illustrated in Figure 2b.

Figure 2

Mortality due to influenza.

In experiments, all else should be equal. We simulated our model without vaccination, with 60% of infants <1 year and adults ≥65 years of age or the same percentage of children aged 1–9, adolescents 10–19, or young adults 20–29 years being vaccinated. These groups are roughly the same size, but coverage actually is <60% among persons <65 years of age [19]. Our hypothetical annual influenza vaccine protected 70% of people 1–64 years of age, but lower proportions of infants and older adults. Efficacy was 35% among infants and declined linearly with age over 64 years (i.e., was 60% among people aged 65–69 years, 50% among those aged 70–74 years, and so on). The resulting efficacies among elderly adults correspond roughly to those reported by Govaert et al. [15]. Annual vaccination occurred November through January; pandemic vaccination began 30 days later and continued for 6 months. Pandemic efficacy was half annual, but as roughly twice as many doses were eventually administered, similar numbers of people were protected.

To assess the impact of these alternative strategies on morbidity and mortality, we averaged daily differences between age-specific cases or deaths with and without vaccination over 365-day periods. Averaging was necessary because our simulation model is stochastic (i.e., we employ Renshaw's discrete event/time method [20]). Finally, in age-structured models, the average number of effective contacts,

₀, may be calculated as the dominant eigenvalue of the next-generation matrix [21] whose associated eigenvector describes the age-specific contributions [22], [23]. We derive these quantities in File S2.

Results

Our matrix of infection rates (Figure 1) illustrates preferential mixing, not only among contemporaries – which is particularly intense among older children, adolescents, and young adults [24] – but also that between parents and children and among co-workers evident in more recent, higher-resolution observations [11], [25], [26]. As indirect effects emanate from off-diagonal matrix elements, such observations increase the accuracy of assessments of intervention impacts via transmission modeling. Similarly, our 1918-like mortalities (Figure 2) resemble those of the 2009 pandemic, although this swine H1N1 was much less virulent than that avian strain.

During simulated pandemic as well as annual influenza outbreaks, vaccinating older children, adolescents, and young adults reduced morbidity the most, especially among target age groups (Figures 3a and b). Despite a contact matrix with relatively large off-diagonal elements, only 20–25% of cases averted were in groups not targeted. By contrast, vaccinating infants and elderly adults reduced mortality most during simulated annual influenza outbreaks (Figure 4a), but vaccinating young adults also reduced mortality during simulated pandemics (Figure 4b).

Figure 3

Cases averted by vaccination.

Figure 4

Deaths averted by vaccination.

While target age groups are similar in size, the numbers of cases averted depend on the vaccine efficacies as well as age distribution of the 2005 U.S. population. Cases averted per efficacy adjusted dose correspond to the proportionate contributions to

₀ (Figure 5), which identifies the optimal target for interventions to reduce transmission.

Figure 5

Normalized age-specific contributions to the reproduction number (File S2).

Discussion

We adapted a classic age-structured population model with parameters chosen to maximize indirect effects due to vaccinating older children, adolescents, and young adults, and to accurately assess direct effects due to vaccinating elderly adults. Comparing the impact of vaccinating these age groups against influenza, we found that vaccinating children, adolescents, and young adults would reduce morbidity the most, with 20–25% of the reduction in other age groups. However, while vaccinating infants and older adults would mitigate mortality most during annual outbreaks, vaccinating young adults also would mitigate mortality during contemporary 1918-like pandemics.

Evidently, which vaccination strategy is superior depends on the objective: mitigating morbidity or mortality, and if mortality, its age-distribution. For many years, U.S. vaccination policy was designed to mitigate mortality, particularly among elderly adults. Relatively recently, it was redesigned to also mitigate morbidity, initially among young children, but then progressively among older children, adolescents, and adults [http://www.cdc.gov/media/pressrel/2010/r100224.htm]. Unlike this policy, in which the 6 month lower age of recommended vaccination has not changed as the upper age has increased, our experimental design maintained similar target group sizes by increasing both lower and upper ages of vaccination simultaneously.

Our findings are comparable to those obtained via other methodologies. The observation that mortality attributed to influenza and pneumonia among elderly Japanese was lower when children were vaccinated routinely [8] suggests that susceptible young people pose a risk to elderly ones, but not necessarily directly. While few such studies are unequivocal, numerous U.S. experiences [27] are consistent with this deduction. Similar conclusions have been reached via community intervention trials [28], [29], [30] as well as individual-based modeling [9]. As our findings support results of these studies using other methodologies, they make a strong case for using relatively simple population models to examine pressing public health issues, and therefore to arrive relatively quickly at sound conclusions about the effectiveness of alternative interventions.

Influenza vaccination strategies have been compared recently using a variety of modeling approaches and perspectives. In 2007, Dushoff et al. [31] explored the same strategies in a 2-group model, one more effective at transmitting the pathogen and other more vulnerable to its effects. These researchers were reluctant to choose among the many interesting scenarios described by various combinations of their parameters, and urged only caution. In 2006, Bansal et al. [32] adopted a more detailed network model with which they also evaluated these strategies, obtaining results qualitatively similar to ours. Three years later, Medlock and Galvani [33] used an age-structured population model with a mixing matrix whose off-diagonal elements are relatively small [25]. Nonetheless, they concluded that vaccinating older children, adolescents, and young adults was the best strategy, regardless of objective.

Impacts of other control measures for pandemic influenza also have been explored recently, by modeling individual members of socially and spatially structured populations [34]–[38]. Our work illustrates several advantages of simpler population models [39]. Insofar as plausible mixing scenarios are modeled, individual behavior is extraneous. Furthermore, systems of equations can be analyzed, whereas computer programs cannot; for example, Areno et al. [40] not only reproduced results with a proportionately-mixed, age-structured population model that had been obtained with a relatively complex individual-based model [9], but also deduced several analytical results. Finally, population models use observations and make predictions familiar to epidemiologists, who group individuals based on characteristics of interest, both in disease surveillance, and to develop and implement interventions. As recently as 2008, for example, Vynnycky and Edmunds used a population model to investigate the impact of school closures on the spread of influenza during a pandemic [41].

Because people of some ages are more active than others, immunizing those potential “super-spreaders” reduces the average number of secondary infections disproportionately. As Figure 5 indicates, adolescents and young adults are the optimal targets for reducing morbidity. Because the main diagonal predominates in all known mixing matrices [11], [23], [25], [26], however, direct effects exceed indirect ones. Unless vaccine efficacy is very low, consequently, the best strategy for reducing mortality will be to vaccinate members of at-risk groups [42]. This analytical result is not limited to vaccination; it may be applied to other interventions that prevent infection or reduce the magnitude or duration of infectiousness. For example, as neuraminidase inhibitors are most effective when administered early [43], timely medication of ill children, adolescents, and young adults could reduce the number needing treatment and possibly the duration of treatments. Treating optimally would be much less costly than widespread prophylaxis, and reduce the risk of drug-resistant strains emerging [44].

Age-specific infection rates are the essence of population models. We calculated risks of infection from Chin et al.'s prospective study of household transmission following illnesses among schoolchildren [10]; households without school-aged children were not represented. Together with clinical observations and individual onset dates, a cross-sectional serological survey would remedy this possible deficiency and might resolve uncertainty about the contribution of asymptomatic infections to transmission. Anderson and May [18] described “who-acquires-infection-from-whom” matrices with as many unique elements as risks of infection, but Nold [45] formulated mixing as a convex combination of age-specific activities (number of contacts per person per day) and constant preference (proportion with others in the same group), and Jacquez et al. [46] allowed preference to vary with age. Recent empirical observations enabled us to include contacts between parents and children and among co-workers [Glasser et al. unpublished manuscript]. Insofar as mixing differs from society to society, if not between rural and urban sub-populations, more diverse subjects would permit continued refinement of methods to permit rapid, robust analysis and interpretation of alternative actions to address public health priorities.

Supporting Information

File S1

A classic population model.

(0.02 MB DOC)

Click here for additional data file.^{(25K, doc)}

File S2

The reproduction number.

(0.06 MB DOC)

Click here for additional data file.^{(60K, doc)}

File S3

Evaluation of the model.

(0.02 MB DOC)

Click here for additional data file.^{(23K, doc)}

Table S1

Demographic parameters, United States, 2005 [12].

(0.07 MB DOC)

Click here for additional data file.^{(70K, doc)}

Acknowledgments

We are grateful to Hsu-Sung Kuo for allowing us to evaluate this model in Taiwan (File S3), to Angus Nicoll for suggesting that we determine which mixing patterns were consistent with historical risks of infection, to Lone Simonsen for guidance about age-specific vaccine efficacy, to Maciej Boni, Gerardo Chowell, Jonathan Dushoff, and Jianhong Wu for constructive reviews of the manuscript, and to Nathaniel Hupert and Michael Washington for their support.

Disclaimer

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention or other institutions with which they are affiliated.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: The authors have no support or funding to report.

References

1. Advisory Committee on Immunization Practices (ACIP) Prevention and Control of Seasonal Influenza. MMWR Recommendations and Reports 2009. 2009;July 31;58(RR-8):1–52. [PubMed]

2. Department of Health and Human Services (HHS) Pandemic Influenza Plan. 2010. Accessed March 23, 2010, at http://www.hhs.gov/pandemicflu/plan/

3. Centers for Disease Control and Prevention. Interim Within-Season Estimate of the Effectiveness of Trivalent Inactivated Influenza Vaccine – Marshfield, Wisconsin, 2007–08 Influenza Season. MMWR Weekly 2008. 2008;April 18; 57(15):393–398. [PubMed]

4. Centers for Disease Control and Prevention (CDC) Seasonal Influenza Vaccine Supply and Distribution in the United States. Accessed March 23, 2010, at http://www.cdc.gov/flu/about/qa/vaxdistribution.htm.

5. Earn DJD, Dushoff J, Levin SA. Ecology and evolution of the flu. Trends Ecol Evol. 2002;17:334–340.

6. ACIP. H1N1 Vaccination Recommendations. 2009. Accessed March 23, 2010 at http://www.cdc.gov/h1n1flu/vaccination/acip.htm.

7. Monto AS, Davenport FM, Napier JA, Francis T. Modification of an outbreak of influenza in Tecumseh, Michigan by vaccinating schoolchildren. J Infect Dis. 1970;122:16–25. [PubMed]

8. Reichert TA, Sugaya N, Fedson DS, Glezen WP, Simonsen L, et al. The Japanese experience with vaccinating schoolchildren against influenza. N Engl J Med. 2001;344:889–896. [PubMed]

9. Longini IM, Halloran ME, Nizam A, Yang Y. Containing pandemic influenza with antiviral agents. Am J Epidemiol. 2004;159:623–633. [PubMed]

10. Chin TDY, Foley JF, Doto IL, Gravelle CR, Weston J. Morbidity and mortality characteristics of Asian strain influenza. Public Health Rep. 1960;75:149–158. [PMC free article] [PubMed]

11. Del Valle SY, Hyman JM, Hethcote HW. Mixing patterns between age groups in social networks. Soc Networks. 2007;29:539–554.

12. CDC. National Vital Statistics System. Public-Use Data Files and Documentation. 2005. Accessed March 23, 2010 at http://www.cdc.gov/nchs/data_access/ftp_data.htm.

13. Luk J, Gross P, Thompson WW. Observations on mortality during the 1918 influenza pandemic. Clin Infect Dis. 2001;33:1375–1378. [PubMed]

14. Bailey NTJ. The Mathematical Theory of Infectious Diseases and Its Applications. London and High Wycombe, England: Charles Griffin and Company; 1975.

15. Govaert TM, Thijs CT, Masurel N, Sprenger MJ, Dinant GJ, et al. The efficacy of influenza vaccination in elderly individuals. A randomized double-blind placebo-controlled trial. JAMA. 1994;272:1661–1665. [PubMed]

16. Koelle K, Cobey S, Grenfell B, Pascual M. Epochal evolution shapes the phylodynamics of interpandemic influenza A (H3N2) in humans. Science. 2006;314:1898–1903. [PubMed]

17. Izurieta HS, Thompson WW, Kramarz P, Shay DK, Davis RL, et al. Influenza and the rates of hospitalization for respiratory disease among infants and young children. N Engl J Med. 2000;342:232–239. [PubMed]

18. Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford, England: Oxford University Press; 1991.

19. CDC, National Center for Health Statistics. NCHS Health E-Stat: Vaccination coverage estimates from the National Health Interview Survey: United States, 2008. 2008. Accessed on March 25, 2010, at http://www.cdc.gov/nchs/data/hestat/vaccine_coverage.htm.

20. Renshaw E. Modeling Biological Populations in Space and Time. Cambridge, England: Cambridge University Press; 1991.

21. Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio R₀ in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28:503–522. [PubMed]

22. Wallinga J, van Boven M, Lipsitch M. Optimizing infectious disease interventions during an emerging epidemic. Proc Natl Acad Sci USA. 2010;107:923–928. [PubMed]

23. 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. [PubMed]

24. Hurford A, Cownden D, Day T. Next-generation tools for evolutionary invasion analyses. J R Soc Interface. 2010;7:561–571. [PMC free article] [PubMed]

25. Mossong J, Hens N, Jit M, Beutels P, Auranen K, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5:381–391. [PMC free article] [PubMed]

26. Zagheni E, Billari FC, Manfredi P, Melegaro A, Mossong J, et al. Using time-use data to parameterize models for the spread of close-contact infectious diseases. Am J Epidemiol. 2008;168:1082–1090. [PubMed]

27. Glezen WP. Emerging infections: pandemic influenza. Epidemiol Rev. 1996;18:64–76. [PubMed]

28. King JC, Jr, Cummings GE, Stoddard J, Readmond BX, Magder LS, et al. A pilot study of the effectiveness of a school-based influenza vaccination program. Pediatrics. 2005;116:e868–873. [PubMed]

29. King JC, Jr, Stoddard JJ, Gaglani MJ, Moore KA, Magderet L, et al. Effectiveness of school-based influenza vaccination. N Engl J Med. 2006;355:2523–2532. [PubMed]

30. Loeb M, Russell ML, Moss L, Fonseca K, Fox J, et al. Effect of influenza vaccination of children on infection rates in Hutterite communities. JAMA. 2010;303:943–950. [PubMed]

31. Dushoff J, Plotkin JB, Viboud C, Simonsen L, Miller M, et al. Vaccinating to protect a vulnerable subpopulation. PLoS Med. 2007;4:921–927. [PMC free article] [PubMed]

32. Bansal S, Pourbohloul B, Meyers LA. A comparative analysis of influenza vaccination programs. PLoS Med. 2006;3:1816–1825.

33. Medlock J, Galvani AP. Optimizing influenza vaccine distribution. Science. 2009;325:1705–1708. [PubMed]

34. Ferguson NM, Cummings DAT, Fraser C, Cajka JC, Cooley PC, et al. Strategies for mitigating an influenza pandemic. Nature. 2006;442:448–452. [PubMed]

35. Germann TC, Kadau K, Longini IM, Macken CA. Mitigation strategies for pandemic influenza in the United States. Proc Natl Acad Sci USA. 2006;103:5935–5940. [PubMed]

36. Glass RJ, Glass LM, Beyeler WE, Min HJ. Targeted social distancing design for pandemic influenza. Emerg Infect Dis. 2006;12:1671–1681. [PMC free article] [PubMed]

37. Haber MJ, Shay DK, Davis XM, Patel R, Jin X, et al. Effectiveness of interventions to reduce contact rates during a simulated influenza pandemic. Emerg Infect Dis. 2007;13:581–589. [PMC free article] [PubMed]

38. Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakulet W, et al. Containing pandemic influenza at the source. Science. 2005;309:1083–1087. [PubMed]

39. May RM. Uses and abuses of mathematics in biology. Science. 2004;303:790–793. [PubMed]

40. Areno J, Brauer F, van den Driessche P, Watmough J, Wu J. Simple models for containment of a pandemic. J R Soc Interface. 2006;3:453–457. [PMC free article] [PubMed]

41. Vynnycky E, Edmunds WJ. Analyses of the 1957 (Asian) influenza pandemic in the United Kingdom and the impact of school closures. Epidemiol Infect. 2008;136:166–179. [PubMed]

42. Nichol KL, Nordin JD, Nelson DB, Mullooly JP, Hak E. Effectiveness of influenza vaccine in the community-dwelling elderly. N Engl J Med. 2007;357:1373–1381. [PubMed]

43. Moscona A. Neuraminidase inhibitors for influenza. N Engl J Med. 2005;353:1363–1373. [PubMed]

44. Lipsitch M, Cohen T, Murray M, Levin BR. Antiviral resistance and the control of pandemic influenza. PLoS Med. 2007;4:111–121. [PMC free article] [PubMed]

45. Nold A. Heterogeneity in disease transmission modeling. Math Biosci. 1980;124:59–82.

46. Jacquez JA, Simon CP, Koopman J, Sattenspiel L, Perry T. Modeling and analyzing HIV transmission: the effect of contact patterns. Math Biosci. 1988;92:119–199.

Articles from PLoS ONE are provided here courtesy of
Public Library of Science