Abstract
Background
Influenza vaccination is vital for reducing H1N1 infectionmediated morbidity and mortality. To reduce transmission and achieve herd immunity during the initial 20092010 pandemic season, the US Centers for Disease Control and Prevention (CDC) recommended that initial priority for H1N1 vaccines be given to individuals under age 25, as these individuals are more likely to spread influenza than older adults. However, due to significant delay in vaccine delivery for the H1N1 influenza pandemic, a large fraction of population was exposed to the H1N1 virus and thereby obtained immunity prior to the wide availability of vaccines. This exposure affects the spread of the disease and needs to be considered when prioritizing vaccine distribution.
Methods
To determine optimal H1N1 vaccine distributions based on individual selfinterest versus population interest, we constructed a game theoretical agestructured model of influenza transmission and considered the impact of delayed vaccination.
Results
Our results indicate that if individuals decide to vaccinate according to selfinterest, the resulting optimal vaccination strategy would prioritize adults of age 25 to 49 followed by either preschoolage children before the pandemic peak or older adults (age 5064) at the pandemic peak. In contrast, the vaccine allocation strategy that is optimal for the population as a whole would prioritize individuals of ages 5 to 64 to curb a growing pandemic regardless of the timing of the vaccination program.
Conclusions
Our results indicate that for a delayed vaccine distribution, the priorities that are optimal at a population level do not align with those that are optimal according to individual selfinterest. Moreover, the discordance between the optimal vaccine distributions based on individual selfinterest and those based on population interest is even more pronounced when vaccine availability is delayed. To determine optimal vaccine allocation for pandemic influenza, public health agencies need to consider both the changes in infection risks among age groups and expected patterns of adherence.
Background
In response to the rapid spread of a pandemic strain of H1N1 influenza A, the World Health Organization (WHO) raised the pandemic alert to its highest phase on June 11, 2009 [1]. The H1N1 pandemic was the first influenza pandemic in over 40 years. Although most H1N1 cases in individuals were mild and the case fatality rate was lower than that of previous influenza pandemics, severe cases frequently occurred in previously healthy, young adults [2].
Vaccines hold considerable promise for reducing the spread of H1N1 influenza A. However, the H1N1 vaccine was not readily available until late October, 2009 [3]. This delayed the US vaccination program until after a large proportion of the population had already been exposed to H1N1.
There is evidence that a substantial proportion of the elderly was protected by crossimmunity from prior infection, resulting in the lowest infection rate in this age group [4]. The 2009 H1N1 influenza disproportionately affected younger patients [5,6]. The median age of hospitalized H1N1 patients was 27 years, which is much lower than the median age of hospitalized seasonalinfluenza cases (between 75 and 79 years) [7,8]. Yet, H1N1 was least likely to turn fatal in patients under age 17 [8]. Such differences in agespecific susceptibility and case fatality for 2009 H1N1 strain posed a challenge to public health agencies that sought to determine optimal vaccine distribution and expected public adherence.
Determining an optimal vaccination policy can be quite challenging. An individual's risk of infection depends not only on his or her decision to be vaccinated, but also on the decisions of others [9,10]. In addition, overwhelming majority of infected people are either asymptomatic or recover without medical attention. Such cases may be unaware that they have been exposed to the virus and still seek vaccination [11]. To calculate the payoff of vaccination to an individual and to the population as a whole, it is important to incorporate the cost of vaccination as well as the benefits of vaccination such as both direct and indirect protection due to herd immunity [10,12,13].
Here, we use game theory to investigate agedependent optimal vaccine distribution against H1N1 influenza A in the US, from both individual and population perspectives. We first model the evolving age distribution of H1N1 cases as the pandemic unfolds, and examine the optimal control strategy assuming that the vaccine becomes available before, at, or after the peak of the influenza pandemic. Then, we find the expected agespecific H1N1 vaccine allocation strategy that would emerge if individuals pursue their own interest, i.e. the Nash strategy, and compare it to a strategy that is optimal to the population as a whole, known as the utilitarian strategy. The personal payoff of vaccination varies among age groups and changes over the course of an outbreak, and we recognize that individuals may not adhere to the utilitarian strategy when acting according to selfinterest.
Our game theoretical analyses of the vaccination program for an influenza A (H1N1) pandemic in the United States show that the utilitarian strategy prioritizes aggressive control among individuals of age 5 and 64 regardless of the timing of vaccination. In contrast, the Nash strategy dictates vaccination of adults, ages 2549, as the first priority group. If the vaccination program implemented before the peak of pandemic wave, then the second priority group to be vaccinated based on the Nash strategy is preschoolage children; however, if vaccination is delayed until the peak of pandemic wave, then the second priority group is older adults (ages 50 to 64).
Methods
To model the transmission of the 2009 H1N1 influenza and vaccination, we developed the agestructured model incorporating six epidemiological compartments (i.e. susceptible, vaccinated, latent, asymptomatic and infectious, symptomatic and infectious, and recovered). Each epidemiological compartment is then subdivided into two depending on an individual’s vaccination decision. The asymptotic dynamics of this model are then used to calculate the probability for individuals to become infected based on their vaccination decision. The expected cost of infection and vaccination associated with vaccine acceptance and refusal are calculated using these infection probabilities. Since the payoff of vaccination depends on both the individual’s decision and the population's average behaviour, we formulate our model as a population game. Monte Carlo methods are employed to determine the optimal vaccination levels driven by selfinterest versus the population interest.
Mathematical model for disease transmission and vaccination
To model H1N1 influenza transmission in the United States, we divide the population into the five age groups (04, 524, 25–50, 50–64, and 65+), according to the age classes used in US CDC case reports [14]. The numbers of people in each age group were set to values estimated for the US 2008 population (Additional File) [15]. In our model, individuals in each age class k are subdivided based on epidemiological status. The dynamics of influenza infection, illness, and infectiousness reflect our current understanding of the natural history of influenza. Here, subscripts U and V represent an unvaccinated and vaccinated population, respectively. We assume that S_{U,k} (t), L_{U,k} (t), A_{U,k} (t), I_{U,k} (t), and R_{U,k} (t) represents the respective number of unvaccinated susceptible, latent, asymptomatic and infectious, symptomatic and infectious, and recovered individuals in age groups k at time t (k = 1, 2, . . . , 5). Similarly, we define S_{V,k} (t), L_{V,k} (t), A_{V,k} (t), I_{V,k} (t), and R_{V,k} (t) as the respective number of vaccinated susceptible, latent, asymptomatic and infectious, symptomatic and infectious, and recovered individuals (k = 1, 2, . . . , 5).
We assume that the vaccine provides partial protection, resulting in vaccinated individuals being less susceptible than unvaccinated ones. Vaccinated individuals become infected at a fraction (1σ_{k}) of the rate at which unvaccinated susceptible individuals become infected, where σ_{k} is the efficacy of the vaccine against infection for individuals of age group k (Additional File). We consider the three vaccination scenarios where vaccines become available before, at, or after the peak of an influenza pandemic. Thus, when vaccines become available, we assume that a proportion, ψ_{k}, of susceptible individuals in age group k is vaccinated. We also assume that the same proportion, ψ_{k}, of individuals in age group k who have been infected asymptomatically still may get vaccinated, because they were not aware of exposure to novel influenza A (H1N1) viruses. However, vaccine doses given to those who were already exposed to H1N1 viruses are assumed to be wasted, because these individuals already gained immunity to the H1N1 strain. Recovered individuals are assumed to be fully protected against further influenza infection for the remainder of the outbreak.
Upon infection, individuals enter a latency period, 1/δ. Latently infected individuals proceed to become infectious, and a proportion, p, of infected individuals becomes symptomatic. Infectious individuals recover after an average period of 1/γ. The inf luenzainduced death rates are α_{U,k} and α_{V,k} for unvaccinated and vaccinated individuals, respectively, for people in age group k. Agespecific influenzarelated death rates are based on estimates of excess pneumonia and on influenza deaths from the H1N1 influenza [16]. The transmission dynamics are thus described by the following differential equations:
for k=1,…, 5.
We used a standardincidence form for the force of infection λ_{k}:
where N is the total population size. Thus, it follows that
where N_{k} is the number of people of age group k, i.e.
Here ϕ_{km} is the number of contacts per day between a person in age group k with people in age group m, and β is the probability of infection for a susceptible person who has contact with an infectious person.
As both epidemiological and serological data are suggestive of residual immunity to H1N1 among adults and seniors, we assume that a proportion (ξ_{k}) of individuals in age group k is immune to H1N1 viruses [4]. The residual immunity incorporates the fact that younger people are more susceptible to the current H1N1 strain than older people due to lack of exposure to a similar virus in the past [4,17]. The demographic effects of aging, birth, and death by causes unrelated to influenza are not included because we only model one influenza season, where these demographic effects are negligible.
The epidemic is initiated with a proportion of each age group assumed to be immune to infection, with one person of each age group assumed infectious, and with the remaining population assumed susceptible. That is,
We assume that an influenza pandemic approaches its peak at time t=ω, and a proportion of the population is vaccinated at time τ=ω±θ where θ=0 or 21 days. We assume that vaccination instantaneously protects people, so that the state variables change discontinuously at t=τ:
with the other state variables remaining the same.
We further assume that the basic reproductive number (R_{0}), defined as the number of secondary cases caused by a single infective case in a completely susceptible population, was 1.4, as estimated for the novel swineorigin H1N1 influenza outbreak [18]. For sensitivity analysis, the basic reproductive ratio was increased from 1.4 to 1.6 (Figure 1). We parameterized agespecific contact rates, ϕ_{km}, using data from a largescale survey of daily contacts [19]. These contact data show strong mixing between people of similar ages and moderately high mixing between children and people of their parents’ ages [20]. Given the contact data and US population size, we reconstructed the contact matrix to match our five age groups [19,20]. Using the relative size of the age group m (N_{m} / N) and the number of contacts per person in age group k with people in age group m, c_{km}, we define the elements of the contact matrix by . To ensure that the number of contacts between age groups is symmetric, N_{m}c_{km} = N_{k}c_{mk}, i.e., ϕ_{km} = ϕ_{mk}, we made further adjustment, , and used ϕ_{km} to be the contact matrix.
Figure 1. (a) Nash and (b) utilitarian strategies when basic reproductive ratio is 1.6. Vaccination is assumed to be offered free of charge. Vaccination is implemented three weeks before, exactly at, or three weeks after the peak of a pandemic influenza.
Cost parameterization
To calculate the average individual net payoff of vaccination strategy, we incorporated the costs associated with infection, vaccination, and the side effects of the vaccine (Table 1). We calculate the cost of infection using weighted average of the costs associated with possible infection outcomes such as mortality, hospitalization, outpatient visits, and cases without medical care. The cost of vaccination includes the value of an individual's time receiving it ($16), and travel cost ($4), resulting in the total cost of vaccination estimated at $20 [21]. The cost of administration is not included in baseline parameters because vaccine for the 2009 novel H1N1 influenza was provided free of charge in the US. However, for sensitivity analysis, we increase the cost of administration from $0 up to $20, in order to examine the elasticity of the Nash and utilitarian strategies to a range of vaccination cost (Figures 2 and 3).
Table 1. Parameterization of infection and vaccination cost
Figure 2. (a) Nash and (b) utilitarian strategies when vaccination administration costs $10
Figure 3. (a) Nash and (b) utilitarian strategies when vaccination administration costs $20
We calculate the cost for vaccine side effects based on the reduction in quality of life and the costs of treating individuals with severe side effects. Mild to moderate side effects are reported to occur at a probability of 5% and to reduce the quality of life by 0.05 for two days on average [21]. To calculate the cost of vaccine side effects, we use the conversion that a qualityadjusted life year (QALY) is monetarily equivalent to $100,000 [2224]. Thus, the cost associated with mild to moderate vaccine side effects can be estimated at
In line with clinical data, severe vaccine side effects are assumed to occur at a probability of 0.001% and result in hospitalization for 7 days, and the cost of hospitalization in ICU to treat severe side effects is taken as $3,739.05 per day [21]. In addition, we assume that severe vaccine adverse effects result in death at 5% of probability [21]. We assume that all individuals value their life equally, irrespective of their age. Thus, the value of life is estimated at $1,045,278 using average expected future lifetime earnings for all ages [25,26]. Estimating the value of life at $1,045,278, the cost of severe side effects is calculated as
Payoff to vaccination strategy
In our vaccination game, the payoff to an individual choosing a particular strategy depends on the average behavior of the population. We considered the two basic strategies, “vaccinator” (obtain vaccination) and “nonvaccinator” (decline vaccination). For both strategies, the payoff to an individual is measured in terms of a monetary cost due to infection and/or vaccination, based on the probability of infections and vaccine risks (Table 1 and Additional File). We also parameterized the payoff calculations with agespecific distributions of vaccine efficacy in reducing influenza morbidity and mortality (Additional File).
The net payoff to vaccinator strategy then is
where x_{k} is the probability of infection among vaccinators, z_{M} is the probability of mild to moderate side effects, and z_{S} is the probability of severe side effects. C_{IV,k} denotes the cost of infection among vaccinators in age group k, C_{V} denotes the cost of vaccination, and C_{M} and C_{S} denote the cost of mild and severe side effects associated with vaccination, respectively.
As the vaccine efficacy is imperfect, the vaccinator may still be infected with reduced probability of infection (x_{k}), which depends on both vaccination probability of age group k (ψ_{k}) and on vaccination probability across all age groups . If infected, vaccinated individuals incur lower infection cost (C_{IV,k}) than unvaccinated ones. The probability of symptomatic infection among vaccinators in age group k who are not yet infected before vaccination is given by
.
Here represents the number of cumulative symptomatic infections until time t = τ. People who had have been symptomatically infected would be aware that they gained immunity against H1N1, thus would not get vaccinated, and therefore the expression, , represents the maximum number of vaccinating people in age group k.
The net payoff to a nonvaccinator is
where C_{IN,k} denotes the cost of infection among nonvaccinators of age group k, and y_{k} is the probability of symptomatic infection among nonvaccinators, given by
.
Here, describes the cumulative number of symptomatic infections among unvaccinated individuals in age group k after vaccination is implemented at time t = τ.
Defining the Nash strategy
For individuals driven by selfinterest, gametheoretic decisions are assumed to settle to a Nash equilibrium at which it is impossible for a few individuals to increase their payoffs by switching to a different strategy [27]. We define these individual decisions at the Nash equilibrium as the Nash strategy. A pure vaccinator strategy cannot be the Nash equilibrium, because when the population vaccine coverage is 100%, an individual who chooses a nonvaccinator strategy reaps the benefits of herd immunity without paying for vaccination and without experiencing possible vaccine side effects. By comparison, a nonvaccinator can result in an individual optimum under certain conditions, such as when the infection risk is sufficiently low when vaccines become available. In our agestructured model, it might be best for some people in an age group to be vaccinated and for others in the same group to choose not to get vaccinated. To allow this scenario, we consider mixed strategies whereby individuals in age group k choose the vaccinator strategy with probability ψ_{k} (0 <ψ_{k} < 1) and the nonvaccinator strategy otherwise. If all individuals play the mixed strategy ψ_{k}, then a proportion ψ_{k} of the population in age group k is vaccinated. The individual optimum can be found by solving for ψ_{k,ind} in the equation U_{vac,k} = U_{nonvac,k} (k=1…5). The individual optimum (ψ_{k,ind}) predicted by this gametheoretical analysis corresponds to the level of coverage ψ_{k,ind} expected under a voluntary program where individuals act in a rational way to maximize their payoffs.
Defining utilitarian strategy
From the perspective of group interest, the objective is to maximize the total payoff of vaccinators and nonvaccinators. If ψ_{k} is the proportion of the population in age group k that is vaccinated, we can express the expected payoff due to vaccination and an influenza pandemic as We now maximize T(ψ_{1}, ψ_{2}, ψ_{3}, ψ_{4}, ψ_{5}) on the parameter space {(ψ_{1}, ψ_{2}, ψ_{3}, ψ_{4}, ψ_{5})  0 ≤ ψ_{k} ≤1} to determine the utilitarian strategy (ψ_{1}*, ψ_{2}*, ψ_{3}*, ψ_{4}*, ψ_{5}*), which is the coverage level that would maximize the total payoff.
Results
Epidemiological impact of the 2009 H1N1 influenza pandemic
Our agestructured model of influenza transmission predicts that 41% of the US population will be infected with pandemic H1N1 influenza in the absence of interventions (Figures 4 and 5). Based on our assumptions that on average 33% of infected people become symptomatic after three day of incubation period [28], we estimate that 13% of the population will be symptomatically infected during the current influenza pandemic, which is consistent with the estimate of previous modeling studies [2931]. However, the agespecific attack rates are predicted to vary considerably between age groups because of agedependent activity patterns and immune profiles. The highest incidence is predicted to occur in individuals of age 524, followed by adult population of age 2549, with symptomatic plus asymptomatic attack rates of 57% and 43%, respectively (Figure 4b). The lowest attack rate (14%) is predicted to occur in the oldest age group (age 65 and older).
Figure 4. Outcomes of influenza A/H1N1 infection. (a) Simulated agestratified daily influenza A/H1N1 infection incidence per 100,000 individuals in the absence of vaccination or other interventions (b) Simulated agespecific attack rates, in the absence of vaccination or other interventions. Both symptomatic and asymptomatic cases are shown.
Figure 5. Cumulative incidence of influenza A/H1N1 when vaccination is guided by the Nash or utilitarian strategies. Vaccination is implemented at free of charge the three weeks before, exactly at, or three weeks after the peak of a pandemic influenza. Solid lines show the cumulative attack rate when vaccination is in alignment with utilitarian strategies when vaccination and is offered free of charge, whereas dotted lines show the cumulative attack rate assuming the vaccination is in alignment with the Nash strategies. For comparison, cumulative incidence without vaccination is also shown.
Our results also suggest that individuals in each age group reach their highest incidence at different times (Figure 4a). That is, schoolage children and young adults (age 524) reach their maximum incidence first, followed by adults (age 2549) and preschoolage children under five years of age. In contrast, the oldest group is the last one to reach maximum incidence.
In the absence of vaccination or other interventions, our model predicts that the 2009 H1N1 influenza pandemic would result in 271 hospitalizations and 13 deaths per 100,000 individuals. With their highest case fatality and hospitalization ratios among all age groups, adults (age 2549) bear the highest case fatality rate (7 out of 13 per 100,000) and hospitalization rate (133 out of 271 per 100,000), followed by schoolage children and young adults of age 524 (Figures 4, 6, 7).
Figure 6. Hospitalizations per 100,000 when vaccination follows the Nash and utilitarian strategies. The number of hospitalizations per 100,000 is shown in age groups with and without vaccination. (a) Vaccination is implemented according to the Nash strategy at different timing of a pandemic wave, i.e. three weeks before, exactly at, or three weeks after the peak of a pandemic influenza. (b) Vaccination is implemented according to the utilitarian strategy at different timing of a pandemic wave, i.e. three weeks before, exactly at, or three weeks after the peak of a pandemic influenza.
Figure 7. Deaths per 100,000 when vaccination follows the Nash and utilitarian strategies. The number of deaths per 100,000 is shown in age groups with and without vaccination. (a) Vaccination is implemented according to the Nash strategy at different timing of a pandemic wave, i.e. three weeks before, exactly at, or three weeks after the peak of a pandemic influenza. (b) Vaccination is implemented according to the utilitarian strategy at different timing of a pandemic wave, i.e. three weeks before, exactly at, or three weeks after the peak of a pandemic influenza.
Optimal H1N1 vaccine distribution based on individual selfinterest
To determine optimal H1N1 vaccine distributions based on individual selfinterest versus population interest, we constructed a game theoretical agestructured model of influenza transmission assuming delayed vaccination. Our calculations show that when vaccination occurs three weeks prior to the peak of a pandemic wave, the Nash (individualbased) strategy prioritizes vaccinating adults (age 25 49) and preschoolage children, followed by schoolage children and young adults (age 524) and then older adults (age 5064) (Figure 8a). The Nash vaccination strategy among senior population of age 65 and older would be to refuse vaccination. With such strategy, the vaccination program is predicted to reduce an overall attack rate from 41% to 15%, averting 8,514 clinical infections, 179 hospitalizations and 9 deaths per 100,000 individuals.
Figure 8. Nash (a) and utilitarian (b) strategies when vaccination is offered free of charge. Vaccination is implemented three weeks before, exactly at, or three weeks after the peak of a pandemic influenza.
The Nash strategy, however, was found to be highly dependent on the timing of vaccine implementation (Figure 8a). If vaccine production is delayed, then the payoff to the vaccinators is diminished because of reduced risk of future infection. For example, if vaccination is implemented at the peak of a pandemic, the Nash vaccination strategy does not include the preschoolage children anymore. Instead, the Nash strategy is to vaccinate 91% of adults (age 2549), 87% of older adults (age 5064), and 63% of schoolage children and young adults. This change in vaccination strategy occurs because preschoolage children have a relatively early pandemic peak and moderate morbidity compared to other age groups. Thus, when a vaccination is significantly delayed, the relative infection risk of this group is low.
The only age groups that are included in the Nash vaccination strategy regardless of vaccine delay are adults (age 2549) and schoolage children/young adults (age 524). In general, the demand for vaccine among adults (age 2549) is the most inelastic to vaccine delay. The Nash level of vaccine coverage for schoolage children and young adults falls rapidly with vaccine delay (Figure 8). For instance, when vaccination is delayed until three weeks after the pandemic peak, the Nash strategy is to vaccinate 88% of adults (age 2549) and 23% of schoolage children/young adults. This vaccination allocation would result in an overall attack rate of 36%, with 229 hospitalizations and 11 deaths per 100,000 individuals (Figures 5, 6, 7).
Our results also demonstrate the dependence of the Nash strategy on basic reproductive ratio of pandemic influenza. At higher transmissibility (R_{0}=1.6), the Nash levels of vaccination are 100% among adults of age 2564 regardless of timing of vaccine implementation. This is, in part, because the case fatality ratio and case hospitalization ratio are highest in these age groups, yielding a high payoff of vaccination (Figure 1). Schoolage children/younger adults (age 524), on the other hand, are expected to seek vaccination only if vaccination is offered on or before the pandemic peak (Figures 9 and 10).
Figure 9. Nash strategy when vaccine is available at the beginning of an influenza pandemic and when vaccination is offered free of charge
Figure 10. Utilitarian strategy when vaccine is available at the beginning of an influenza pandemic and when vaccination is offered free of charge
Finally, we considered the rising cost of vaccination and its impact on the Nash strategy of each age group. We show that the Nash vaccination level among adults (age 2549) is the most inelastic to the changes in vaccination cost (Figures 2 and 3). This inelasticity arises because, in this age group, the infection risk and casefatality rate of H1N1 are high and residual immunity is low. In contrast, the Nash vaccination of older adults (age 50 to 64) is the most elastic to the changes in the vaccination cost, demonstrating the tradeoff between vaccine cost and reduced benefit of vaccine due to vaccine delay. This elasticity is because the infection risk in this age group is relatively low, resulted from their residual immunity against H1N1 and low contact rate. Seniors of age 65 or older are unlikely to seek vaccination if vaccination is voluntary at a wide range of vaccination costs because their risk of infection (and thus their vaccination payoff) is lowest among all age groups.
Optimal H1N1 vaccine distribution based on population interest
The average vaccination level across all age groups for the utilitarian strategy is higher than that for the Nash strategy (Figure 8). For example, if vaccines become available three weeks before the pandemic peak, the overall Nash and utilitarian vaccine coverage are 76% and 82%, respectively. When 93% of young individuals (age under 24), 96% of adults (age 2564) and 1% of seniors (age 65 and older) are vaccinated according to the utilitarian strategy three weeks prior to the peak of the influenza pandemic, 25,767 clinical infections, 182 hospitalizations and 9 deaths would be averted per 100,000 individuals (Figures 5, 6, 7).
The utilitarian strategy is, however, much less effective if vaccination is delayed. For instance, if vaccination is delayed until the peak of influenza pandemic, it is estimated that 15,683 infections, 112 hospitalizations and 6 deaths would be averted per 100,000 individuals, which is considerably fewer than when the vaccination is implemented before the pandemic peak (Figures 5 and 8).
We find that the utilitarian vaccine coverage levels are more inelastic than those under the Nash strategy. For instance, if vaccination is delayed until three weeks after the pandemic peak, the resulting vaccine coverage level according to the Nash strategy falls to 37% whereas the utilitarian strategy is to still vaccinate 73% of population. Thus, the resulting disease incidence and the number of diseaserelated deaths are lower under the utilitarian strategy than under the Nash strategy. The utilitarian strategy includes vaccinating 90% of older adults (age 5064), 88% of adults (age 2549), 83% of schoolage children/young adults (age 524) and 43% of preschoolage children. At a higher basic reproductive ratio of 1.6, the utilitarian strategy also includes the vaccination of preschoolage children (age 04), because the risk of infection increases with transmissibility of influenza virus, increasing the payoff of vaccination (Figure 1). However, the Nash vaccination strategy does not include preschoolage children or older adults, and thus the utilitarian coverage levels may be unachievable under voluntary vaccination.
Conclusions
For pandemic H1N1, we find that the individualbased (Nash) vaccination strategies differ significantly from the utilitarian vaccination strategies. Without vaccination delay, the primary priority group under the utilitarian strategy is schoolage children and young adults (age 524) because of their important role in transmitting disease (Figures 9 and 10). The case hospitalization ratio and the case fatality ratio are the highest, and thus vaccinating these individuals yields high individual and population payoffs. Indeed, regardless of length of the delay and when vaccination is guided by the Nash or utilitarian strategies, younger adults (age 2549) are among the highest priority groups for vaccination.
However, the second priority group changes dramatically under the Nash strategy. If vaccination occurs before the pandemic peak, the second Nash priority group is preschoolage children. If vaccination is delayed, the second Nash priority group is shifted to older adults (age 5064) or to schoolage children/younger adults (age 524). The peak incidence among preschoolage children is relatively early compared to other age groups, thus lowering the benefit of vaccination to these children with time. Because the case fatality ratio is the highest among older adults, and H1N1 morbidity is the highest among schoolage children/younger adults, the benefit of vaccination is relatively inelastic over the course of a pandemic. Therefore, the demand for vaccines among these age groups is high even if vaccination is delayed in a pandemic.
The discordance between the Nash and utilitarian strategies is even more pronounced when vaccine availability is delayed. If vaccination is delayed but implemented near the pandemic peak, the utilitarian vaccination strategy includes individuals of age up to 64, in contrast to the Nash strategy which excludes preschoolage children and older adults (age 5064) (Figure 8). If vaccination is further delayed, the Nash strategy would also exclude adults (age 2549), preschoolage children and older adults (age 5064), whereas the utilitarian vaccination strategy still includes individuals of age up to 64. Therefore, the average vaccination level across all age groups at the utilitarian strategy was found to be higher than that at the Nash strategy.
Overall, our results indicate that a vaccination levels under a voluntary immunization program may not be optimal for the population, regardless of vaccine delay. Such discordance between the Nash and utilitarian strategies is predicted to be robust to the increase in the basic reproductive ratio for pandemic influenza (Figure 1). This finding is consistent with those of previous studies, which demonstrated that, in the context of vaccination against smallpox and seasonal influenza, the vaccination levels driven by selfinterest are likely to be lower than those that are optimal from the population perspective [3235].
There are three primary reasons for the discrepancy between the individualbased and utilitarian agespecific vaccination levels for pandemic H1N1. First, different age groups have different incentives to vaccinate. In particular, an earlier pandemic peak among young individuals results in a relatively low infection risk later in the pandemic compared to that for older adults. Therefore, the young are predicted to undervaccinate under the Nash strategy relative to the utilitarian strategy when vaccination is delayed. Second, the positive externalities of indirect protection by herd immunity also contribute to the differences between utilitarian and Nash vaccination strategies. The benefits of herd immunity contribute to the utilitarian strategy, but also create an incentive for individuals to free ride on the vaccination of others. Consequently, the overall level of population vaccination is lower for the Nash strategy than for the utilitarian strategy. Third, because vaccine delivery was delayed for the H1N1 pandemic, our model predicts that people will be less inclined to vaccinate than if vaccine was available at the beginning of the pandemic. As a consequence, achieving vaccination rates high enough to achieve the utilitarian strategy may be difficult, and the discordance between the Nash and utilitarian strategies is found to increase with vaccine delay.
The guidelines for vaccinating against the 20092010 pandemic H1N1 influenza proposed by the CDC’s Advisory Committee on Immunization Practices (ACIP) prioritize young people aged 6 months to 25 years, who are the most efficient at transmitting influenza viruses [36]. This guideline also reflects the reduced susceptibility among the elderly due to their residual immunity from past exposure [37]. If large stockpiles of vaccines had been available prior to the pandemic, the optimal vaccine distribution strategy would be to vaccinate children in order to reduce transmission and achieve herd immunity [38,39]. However, our analysis suggests that the success of such vaccination strategies depends heavily on the timing of a vaccine’s availability. Nevertheless, our analysis might be limited by the difficulties of knowing the state of the pandemic at the time vaccines become available. In addition, our outcome measure (i.e. cost of infection and vaccination) may oversimplify the vaccination decisions or be incongruous with the consideration of the Advisory Committee on Immunization Practices (ACIP).
We found that, for both the Nash and utilitarian strategies, the optimal vaccination strategies with vaccine delay should prioritize individuals of age 25 to 49. Our results also suggest that a utilitarian vaccine strategy should also include individuals from a wide range of ages, from 5 months to 65 years; and for longer delay in vaccination, vaccination priority should increasingly be given to older individuals. Our results further suggest that agespecific demands for vaccination depend on the risk of infection at the time of vaccine delivery and the severity of the disease. When vaccination is delayed, voluntary adherence to vaccine recommendation might become lower among young individuals. This suggests that influenza pandemic response plans should include efforts to encourage the vaccination of young individuals if vaccine delivery is delayed.
List of abbreviations
H1N1:(H: Hemaglutinin) (N: Neuraminadase); WHO: the World Health Organization; CDC: Centers for Disease Control and Prevention
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
ES developed and analyzed the model, and carried out numerical simulations. ES produced all figures, interpreted results, and wrote the manuscript. AG suggested some of the simulations and helped write the manuscript. LM edited the manuscript. All authors read and approved the final manuscript.
Acknowledgements
We are grateful to funding from NIH grant U01 GM087719. ES was also partially supported by the National Institute of General Medical Sciences MIDAS grant 5U54GM08849102 and by the Vaccine Modeling Initiative (VMI) funded by the Bill and Melinda Gates Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
This article has been published as part of BMC Public Health Volume 11 Supplement 1, 2011: Mathematical Modelling of Influenza. The full contents of the supplement are available online at http://www.biomedcentral.com/14712458/11?issue=S1.
References

Eurosurveillance editorial team: Pandemic alert level 6: scientific criteria for an influenza pandemic fulfilled.
Euro Surveill 2009, 14(23):19237. PubMed Abstract  Publisher Full Text

Maritz J, Maree L, Preiser W: Pandemic influenza A (H1N1) 2009: the experience of the first six months.
Clin Chem Lab Med 48(1):1121. PubMed Abstract  Publisher Full Text

CDC: 38 Million Doses of H1N1 Vaccine Available [http://www.flu.gov/news/blogs/vaccine38million.html] webcite

PerezTrallero E, Pineiro L, Vicente D, Montes M, Cilla G: Residual immunity in older people against the influenza A(H1N1)recent experience in northern Spain.

Bansal S, Pourbohloul B, Hupert N, Grenfell B, Meyers LA: The shifting demographic landscape of pandemic influenza.
PLoS One 5(2):e9360. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Lee BY, Brown ST, Korch GW, Cooley PC, Zimmerman RK, Wheaton WD, Zimmer SM, Grefenstette JJ, Bailey RR, Assi TM, et al.: A computer simulation of vaccine prioritization, allocation, and rationing during the 2009 H1N1 influenza pandemic.
Vaccine 28(31):48754879. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Thompson WW, Shay DK, Weintraub E, Brammer L, Bridges CB, Cox NJ, Fukuda K: Influenzaassociated hospitalizations in the United States.
JAMA 2004, 292(11):13331340. PubMed Abstract  Publisher Full Text

Louie JK, Acosta M, Winter K, Jean C, Gavali S, Schechter R, Vugia D, Harriman K, Matyas B, Glaser CA, et al.: Factors associated with death or hospitalization due to pandemic 2009 influenza A(H1N1) infection in California.
JAMA 2009, 302(17):18961902. PubMed Abstract  Publisher Full Text

Shim E, Chapman GB, Galvani AP: Decision making with regard to antiviral intervention during an influenza pandemic.
Med Decis Making 2010, 30(4):E6481. PubMed Abstract  Publisher Full Text

Shim E, Kochin B, Galvani A: Insights from epidemiological game theory into genderspecific vaccination against rubella.
Math Biosci Eng 2009, 6(4):839854. PubMed Abstract  Publisher Full Text

Influenza A (H1N1): pandemic alert phase 6 declared, of moderate severity [http:/ / www.euro.who.int/ en/ whatwepublish/ informationforthemedia/ sections/ pressreleases/ 2009/ 06/ influenzaah1n1whoannouncespand emicalertphase6,ofmoderatesev erity] webcite

Shim E, Galvani AP: Impact of transmission dynamics on the costeffectiveness of rotavirus vaccination.
Vaccine 2009, 27(30):40254030. PubMed Abstract  Publisher Full Text

Shim E, CastilloChavez C: The epidemiological impact of rotavirus vaccination in the United States and Mexico. In Mathematical and Statistical Estimation Approaches in Epidemiology. Edited by Chowell G, Hyman JM, Bettencourt LMA, CastilloChavez C. Springer; 2009:303323. Publisher Full Text

2009 H1N1 Early Outbreak and Disease Characteristics [http://www.cdc.gov/H1N1FLU/surveillanceqa.htm] webcite

Boni MF, Manh BH, Thai PQ, Farrar J, Hien TT, Hien NT, Van Kinh N, Horby P: Modelling the progression of pandemic influenza A (H1N1) in Vietnam and the opportunities for reassortment with other influenza viruses.
BMC Med 2009, 7:43. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Presanis AM, Lipsitch M, Daniela DA, Hygiene NYCDoHaM, Team SFI, Hagy A, Reed C, Riley S, Cooper B, Biedrzycki P, et al.: The severity of pandemic H1N1 influenza in the United States, April  July 2009.

Kelly HA, Grant KA, Williams S, Fielding J, Smith D: Epidemiological characteristics of pandemic influenza H1N1 2009 and seasonal influenza infection.
Med J Aust 2009, 191(3):146149. PubMed Abstract  Publisher Full Text

Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, et al.: Pandemic potential of a strain of influenza A (H1N1): early findings.
Science 2009, 324(5934):15571561. PubMed Abstract  Publisher Full Text

Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, et al.: Social contacts and mixing patterns relevant to the spread of infectious diseases.
PLoS Med 2008, 5(3):e74. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Medlock J, Galvani AP: Optimizing influenza vaccine distribution.
Science 2009, 325(5948):17051708. PubMed Abstract  Publisher Full Text

Khazeni N, Hutton DW, Garber AM, Hupert N, Owens DK: Effectiveness and costeffectiveness of vaccination against pandemic influenza (H1N1) 2009.

Burstrom K, Johannesson M, Diderichsen F: The value of the change in health in Sweden 1980/81 to 1996/97.
Health Econ 2003, 12(8):637654. PubMed Abstract  Publisher Full Text

Cutler DM, Richardson E: Measuring the health of the US population. Brookings paper on economic activity.

Meltzer MI, Cox NJ, Fukuda K: The economic impact of pandemic influenza in the United States: priorities for intervention.
Emerg Infect Dis 1999, 5(5):659671. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Haddix AC, Teutsch SM, Shaffer PA, Duet DO: Prevention Effectiveness. New York, NY: Oxford University Press; 1996.

Nash JF: Equilibrium points in Nperson games.
PNAS 1950, 36:4849. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Basta NE, Chao DL, Halloran ME, Matrajt L, Longini IM Jr.: Strategies for pandemic and seasonal influenza vaccination of schoolchildren in the United States.
Am J Epidemiol 2009, 170(6):679686. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

CDC Estimates of 2009 H1N1 Influenza Cases, Hospitalizations and Deaths in the United States, April 2009 – March 13, 2010 [http://www.cdc.gov/h1n1flu/estimates_2009_h1n1_April_March13.htm] webcite

Towers S, Feng Z: Pandemic H1N1 influenza: predicting the course of a pandemic and assessing the efficacy of the planned vaccination programme in the United States.
Euro Surveill 2009, 14(41):19358. PubMed Abstract  Publisher Full Text

La Ruche G, Tarantola A, Barboza P, Vaillant L, Gueguen J, GastelluEtchegorry M: The 2009 pandemic H1N1 influenza and indigenous populations of the Americas and the Pacific.

Bauch CT, Galvani AP, Earn DJ: Group interest versus selfinterest in smallpox vaccination policy.
Proc Natl Acad Sci U S A 2003, 100(18):1056410567. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Bauch C, Earn D: Vaccination and the theory of games.
Proc Natl Acad Sci USA 2004, 101:1339113394. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Bauch CT: Imitation dynamics predict vaccinating behaviour.
Proc Biol Sci 2005, 272(1573):16691675. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Galvani AP, Reluga TC, Chapman GB: Longstanding influenza vaccination policy is in accord with individual selfinterest but not with the utilitarian optimum.
PNAS 2007, 104(13):56925697. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Use of influenza A (H1N1) 2009 monovalent vaccine: recommendations of the Advisory Committee on Immunization Practices (ACIP), 2009
MMWR Recomm Rep 2009, 58(RR10):18. PubMed Abstract  Publisher Full Text

Dawood FS, Jain S, Finelli L, Shaw MW, Lindstrom S, Garten RJ, Gubareva LV, Xu X, Bridges CB, Uyeki TM: Emergence of a novel swineorigin influenza A (H1N1) virus in humans.
N Engl J Med 2009, 360(25):26052615. PubMed Abstract  Publisher Full Text

Mathematical modelling of the pandemic H1N1 2009
Wkly Epidemiol Rec 2009, 84(34):341348. PubMed Abstract

Yang Y, Sugimoto JD, Halloran ME, Basta NE, Chao DL, Matrajt L, Potter G, Kenah E, Longini IM Jr.: The transmissibility and control of pandemic influenza A (H1N1) virus.
Science 2009, 326(5953):729733. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Clark TW, Pareek M, Hoschler K, Dillon H, Nicholson KG, Groth N, Stephenson I: Trial of influenza A (H1N1) 2009 monovalent MF59adjuvanted vaccine  preliminary report.
N Engl J Med 2009. PubMed Abstract  Publisher Full Text

Dasta JF, McLaughlin TP, Mody SH, Piech CT: Daily cost of an intensive care unit day: the contribution of mechanical ventilation.
Crit Care Med 2005, 33(6):12661271. PubMed Abstract  Publisher Full Text

CBA Case Study: The Economic Impact of Pandemic Influenza in the United States [http://www.cdc.gov/owcd/EET/CBACase/1.html] webcite

Campbell DS, Rumley MA: Costeffectiveness of the influenza vaccine in a healthy, workingage population.
J Occupational Environ Med 1997, 39:408414. Publisher Full Text

Kavet J: A perspective on the significance of pandemic influenza.
Am J Public Health 1977, 67:10631070. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Serfling RE, Sherman IL, Houseworth WJ: Excess pneumoniainfluenza mortality by age and sex in three major influenza A2 epidemics, United States, 195758, 1960 and 1963.
Am J Epidemiol 1967, 86(2):433441. PubMed Abstract

Coburn BJ, Wagner BG, Blower S: Modeling influenza epidemics and pandemics: insights into the future of swine flu (H1N1).
BMC Med 2009, 7:30. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Loughlin J, Wang FT, El Khoury A, Doherty M, Mast TC, Seeger JD: Adherence to the recommended pentavalent rotavirus vaccine administration schedule in the United States: 20062007.
Pediatr Infect Dis J 2009, 28(7):667668. PubMed Abstract  Publisher Full Text

Presanis AM, Lipsitch M, Angelis DD, The_Swine_Flu_Infestigation_Team, N.Y.C.Department_of_Health_and_Mental_Hygiene, Hagy A, Reed C, Riley S, Cooper B, Biedrzycki P, et al.: The severity of pandemic H1N1 influenza in the United States, April–July 2009.

Meltzer MI, Cox NJ, Fukuda K: Modeling the economic impact of pandemic influenza in the United States: implications for setting priorities for intervention: background paper. Technical Report.
1999.

Neustadt RE, Fineberg HV: The Swine Flu Affair: Decision Making on a Slippery Disease. Washington, DC: National Academies Press Online; 1978.