Skip to main content

Volume 11 Supplement 1

Mathematical Modelling of Influenza

Reactive strategies for containing developing outbreaks of pandemic influenza

Abstract

Background

In 2009 and the early part of 2010, the northern hemisphere had to cope with the first waves of the new influenza A (H1N1) pandemic. Despite high-profile vaccination campaigns in many countries, delays in administration of vaccination programs were common, and high vaccination coverage levels were not achieved. This experience suggests the need to explore the epidemiological and economic effectiveness of additional, reactive strategies for combating pandemic influenza.

Methods

We use a stochastic model of pandemic influenza to investigate realistic strategies that can be used in reaction to developing outbreaks. The model is calibrated to documented illness attack rates and basic reproductive number (R0) estimates, and constructed to represent a typical mid-sized North American city.

Results

Our model predicts an average illness attack rate of 34.1% in the absence of intervention, with total costs associated with morbidity and mortality of US$81 million for such a city. Attack rates and economic costs can be reduced to 5.4% and US$37 million, respectively, when low-coverage reactive vaccination and limited antiviral use are combined with practical, minimally disruptive social distancing strategies, including short-term, as-needed closure of individual schools, even when vaccine supply-chain-related delays occur. Results improve with increasing vaccination coverage and higher vaccine efficacy.

Conclusions

Such combination strategies can be substantially more effective than vaccination alone from epidemiological and economic standpoints, and warrant strong consideration by public health authorities when reacting to future outbreaks of pandemic influenza.

Background

In April, 2009, the World Health Organization (WHO) announced the emergence of a new influenza A (H1N1) virus, and on June 11, 2009, it declared that the world was at the start of a new influenza pandemic [1]. WHO reported more than 414,000 laboratory-confirmed cases of H1N1 [2] — a gross underestimate, as many countries simply stopped counting individual cases. The US Centers for Disease Control and Prevention reported widespread influenza activity in forty-six states, with influenza-like illness (ILI) activity in October 2009 higher than what is seen during the peak of many regular flu seasons; and further, “Almost all of the influenza viruses identified … are 2009 H1N1 influenza A viruses” [3]. Countries found themselves in the position of having to react to contain already developing Fall outbreaks of influenza due to the new pandemic strain, a position they are likely to find themselves in again if and when future waves of pandemic influenza occur.

Research has suggested that mass vaccination of 60–70% of the population prior to the start of the flu season could effectively contain outbreaks due to pandemic strains [47]; and the public health preparedness plans of most countries have, accordingly, emphasized vaccination intervention strategies. However, the recent experience with H1N1 suggests that high vaccination coverage levels are difficult to achieve. In the case of H1N1, vaccination programs in most northern hemisphere countries started only after the virus was widely circulating. Furthermore, in some countries, supplies of vaccine were limited [8], delivery and administration occurred over a period of several months [9, 10], and there were reports of public skepticism regarding the necessity and safety of vaccination [11, 12], all of which were strong indicators suggesting that high vaccination coverage would be difficult to achieve. While many institutions in the US and elsewhere strongly encouraged and, in some cases, required workers to be vaccinated against seasonal influenza in 2009, H1N1 vaccination guidelines were focused mostly on people in certain age and high-risk groups [13]. Delays, limited and untimely vaccination supplies, and public reluctance to be vaccinated are likely to reduce the effectiveness of vaccination campaigns [4, 5].

The issues outlined above for the recent outbreak of H1N1 are likely to occur again in future outbreaks of pandemic influenza. In this paper, we explore the effectiveness of realistic reactive intervention strategies implemented after the beginning of outbreaks of pandemic influenza. We calibrate our model based on data for the H1N1 pandemic (see Tuite et al. [14]), and we investigate the impacts of (i) the moderate vaccination coverage levels which, based on past experience, are likely to be realized, as well as high levels which would be more ideal; (ii) very limited treatment of cases with antivirals and prophylaxis of cases’ households with antivirals; and (iii) limited and practical social distancing measures such as five-day closure of individual schools on an as-needed basis, encouragement of liberal leave policies in the workplace, and encouragement of self-isolation. Intervention strategies that combine these approaches are also studied (cf. Halloran et al. [15]). For all intervention strategies, we provide cost estimates associated with morbidity and mortality that take into account direct medical costs as well as economic consequences resulting from school closures and work loss.

Methods

The simulation model

We developed a portable and adaptable stochastic, individual-level simulation model of influenza spread within a structured population. The simulator is similar to models developed by Longini et al. [7, 16]. The simulation population of 649,565 people was generated stochastically to represent a typical North American city, namely, Hamilton (Ontario), Canada, which was chosen due to availability of demographic and epidemiological data necessary for constructing and calibrating the simulator. Our population is a collection of heterogeneous individuals with various attributes that impact whom they interact with (and hence whom they may infect or get infected by). More specifically, each individual has the following stochastically generated attributes: age, household, playgroup or daycare attended (for pre-school children), school attended (for school-age children), workgroup (for working adults), household census tract and workplace census subdivision, community, and neighborhood. As in [16], a community consists of approximately 2000 people living within the same census tract, and a neighborhood consists of approximately 500 people living within proximity to each other within the same community; also see the recent papers [17] and [18], which incorporate more-detailed individual-level behavior involving larger populations. Age and household-size distributions, shown in Figures 1 and 2, were matched to 2001 Canadian census data [19, 20]. Household census tract assignments were made so that census tract population sizes were consistent with 2006 census statistics [21]. Workgroups were formed to match 2006 employment statistics [22] as well as census statistics on the geographical distribution of workers [23]. Rather than representing entire workplace institutions, we formed workgroups of size 20 to represent the typical number of co-workers an individual is likely to have close contact with during the day. Average playgroup, daycare, and lower and upper secondary school (i.e., middle and high school) contact group sizes were chosen for similar reasons; see the Appendix.

Figure 1
figure 1

Age distribution for simulated population

Figure 2
figure 2

Household size distribution for simulated population

Susceptible people are assumed to have daily contacts with other individuals in their contact groups, i.e., their household and school or workgroups, as well as with people in their neighborhood and community. Infection of susceptibles depends on the number of infected persons in their contact groups, on the vaccine and antiviral-use status of susceptibles and their infectious contacts, and on age- and contact-group-specific per-contact transmission probabilities (Table 1). This disease transmission model is based on previously described models [7, 16], and is detailed in the Appendix. People infected with influenza first pass through a latent / incubation period, during which they do not have influenza symptoms. They are not infectious until the last day of the period; at that point, they become half as infectious as if they were to develop symptoms in the subsequent period. During that subsequent infectious period, 67% will develop influenza symptoms and 33% will be asymptomatic (and will be half as infectious as those who are symptomatic) [7]. The model allows for people to withdraw from all of their mixing groups, except the household, if they become infected or have an infected child.

Table 1 Per-contact influenza infection transmission probabilities within contact groups

The simulator is calibrated to match documented illness attack rates and basic reproductive numbers (R0). Baseline (no-intervention) scenario age-group-specific attack rates were derived using 2009 estimates for the H1N1 basic reproductive number in Ontario [14, 24, 25] (see Table 2). These rates take into account reduced risk in adults born prior to 1957 [24]. A compartmental model parameterized in this way was well-calibrated to observed attack rates during the Fall pandemic wave in Ontario [25]. The simulator’s R0 value of 1.4 is also consistent with other published reports [4, 26, 27].

Table 2 Age-group-specific H1N1 influenza illness attack rates in Ontario, Canada, 2009, and calibrated attack rates

Intervention strategies

We modeled a baseline case where no intervention takes place, along with strategies representing various combinations of vaccination, antiviral treatment and household prophylaxis, school closure, and general social distancing (see the results in Tables 3 and 4 and Supplementary Data Table S1 provided in “Additional File 1”). Each component of the strategies is described in detail below. Interventions are triggered in a particular simulation run when the overall illness attack rate reaches 0.01%. Twenty runs of the simulator were performed for each intervention strategy, from which average illness attack rates were calculated. We briefly describe the interventions under consideration.

Vaccination

We model both pre-vaccination as well as reactive strategies, with reactive vaccination programs beginning immediately, 30 days, or 60 days after the trigger. The delays model disruptions in vaccine production and supply chains. We allow enough doses to cover either 35% or 70% of the population. In reactive strategies, we consider cases where (i) all vaccines become available at the same time, and (ii) the doses become available in three equal-sized batches, two weeks apart, due to additional production and supply-chain disruptions. We study a low-efficacy single-dose vaccine (efficacy against susceptibility to infection, VEs = 0.3, and efficacy against infectiousness, VEi = 0.2) as well as a moderate-efficacy vaccine (VEs = 0.4, VEi = 0.5) [28]. Vaccine efficacy refers to the reduction, after vaccination, in the probability of becoming infected due to contact with an infected person (VEs), or to the reduction, after vaccination, in the probability of infecting a susceptible contact (VEi). Vaccine efficacy does not refer to the fraction of individuals having an immunogenic response to the vaccine (which is typically much larger than our measures).

Each day, our model randomly vaccinates any remaining unvaccinated individuals who are either uninfected or in the latent or asymptomatic phases of infection, all with equal probability based on the number of available doses. Moreover, protection from the vaccine builds over time, with 50% of the vaccine’s efficacy realized upon vaccination, and full protection after two weeks.

Antiviral treatment and household prophylaxis

We investigate strategies involving treatment of infected individuals with a five-day course of antivirals, as well as strategies that also allow for ten-day prophylaxis of the infected individuals’ household members. We assume that 1% of individuals do not complete their course. We use an antiviral efficacy against susceptibility (AVEs) of 0.3 and against infectiousness (AVEi) of 0.7 [16]. Individuals receive direct benefit from antivirals only while they are taking them. Antiviral use is considered alone and in combination with other intervention strategies. It is assumed that antiviral courses are available for 10% of the population and that they are distributed to infected individuals and their household members until the supply is exhausted.

School closure and social distancing

We implement a rolling school closure model, where a daycare or school closes for five days if five or more cases are identified in that group. Given that infected individuals are on average infectious for 4.1 days (see Figure 3), closing schools for fewer than 5 days is unlikely to be very effective. It is possible for these groups to close more than once during the simulation. We also model a reduction in workplace and general community contacts of 20% (i.e., 20% of infected individuals in each contact group will not infect other members of the group). This represents the exercise of a general level of caution, including a modest limitation of contacts within workgroups (e.g., by invoking occasional telecommuting and other self-limiting behaviors, holding fewer large meetings, etc.) and also within the general community (e.g., reduction in attendance in social groups and larger community events, etc.).

Figure 3
figure 3

Simulation flowchart and modeled influenza natural history

Economic cost estimation

We determine economic costs associated with the influenza outbreaks and modeled intervention strategies using methods described by Meltzer et al. [29]. We include medical spending due to illness, costs of antivirals and vaccines, and costs associated with teachers and other working adults staying home due to their own illness, illness of dependent children, or due to school closure. Medical spending includes co-payments and net payments for outpatient visits and hospitalization, as well as prescription and over-the-counter medications for influenza and complications or secondary infections. Costs are stratified by age-group and by low- or high-risk status of individuals with respect to complications of influenza. We also include the present value of earnings lost due to premature mortality.

Cost estimates and probabilities of risk status and of complications and death were taken from Meltzer et al. [29], with costs inflated using 2008 consumer price index and medical price index estimates [3033]. These costs are combined with the data on age-specific attack rates, utilized vaccination doses, and days of school closure obtained from our simulation model. Details of the cost calculations are given in the Appendix.

Results

With no intervention, the average overall illness attack rate is 34.1%, with an estimated total cost of $81.1 million (Table 3). Pre-vaccination of 35% of the population with a low-efficacy vaccine reduces the average overall illness attack rate to 26.1% (total cost $71.1 million), and with a moderate-efficacy vaccine to 18.8% (total cost $53.7 million). Not surprisingly, pre-vaccination of 70% of the population is more effective (overall average illness attack rate 12.0%, total cost $47.0 million for a low-efficacy vaccine; and 0.2% and $19.3 million with a moderate-efficacy vaccine; see Table 4).

Table 3 Average overall illness attack rates and total costs of interventions with 35% vaccination coverage
Table 4 Average overall illness attack rates and total costs of interventions with 70% vaccination coverage

Reactive vaccination alone, of 35% of the population with a low-efficacy vaccine delivered in three batches, reduces the overall average illness attack rate to 28.8% (or 22.8% with a moderate-efficacy vaccine), with a total cost of $77.7 million ($63.1 million with a moderate-efficacy vaccine). Thirty- and 60-day delays in initiation of reactive vaccination, with vaccines delivered in three batches, result in attack rates of 29.5% (total cost $79.3 million) and 32.2% (total cost $86.0 million), respectively, for a low-efficacy vaccine, and 24.6% (total cost $67.5 million) and 30.8% (total cost $82.5 million), respectively, for a moderate-efficacy vaccine. Clearly, with a 60-day delay, interventions occur too late in the epidemic to have any meaningful effect (see Figure 4).

Figure 4
figure 4

Daily attack rates for (i) the case of 70% coverage of low-efficacy vaccine with 60-day initial delay, and (ii) the baseline case. For case (i), the vaccine is given on the 60th day followed by receipt of vaccine after two additional two-week delays (see arrows). Note that vaccine given on the 60th day decreases the attack rate compared to the baseline; but the two subsequent receipts of vaccine do not result in additional benefits.

Antiviral use at low (10%) coverage alone results in an overall attack rate of 31.3% (total cost $75.9 million). School closure and social distancing alone result in an attack rate of 24.0%, with a total cost of $125.0 million.

Suppose we combine reactive low-efficacy vaccination of 35% of the population delivered in three batches, antivirals (10% coverage), and school closure and social distancing. Then the overall average illness attack rate is 4.5% (total cost $32.2 million) if no delays occur in the initiation of vaccination, and 5.4% (total cost $36.8 million) if a 60-day delay occurs. With a moderate-efficacy vaccine, the attack rate for this last scenario reduces to 2.4% (total cost $22.0 million). Similar relationships between interventions are apparent for interventions with 70% vaccination coverage, shown in Table 4. Vaccination coverage of 70% with a moderate-efficacy vaccine, combined with antiviral treatment and school closure, is highly effective, even with an initial 60-day delay and additional supply-chain disruptions (average illness attack rate 1.4%, total cost $27.4 million).

We note that the results when all vaccines are available at the same time are better than those involving delivery in batches, and sometimes significantly so, especially for a moderate-efficacy vaccine (Tables 3 and 4). Figures 5A through 5D illustrate the comparative illness attack rates of the various intervention strategies discussed above for all combinations of low/moderate-efficacy vaccine delivered in three batches and at 35% / 70% coverage as a function of the initial delay in vaccination implementation due to supply-chain disruptions. The impact of vaccinating 70% of the population, rather than 35%, ranges from moderate to substantial, with the increased coverage being most beneficial when the vaccine is delivered in a timely manner, and the vaccine is either of moderate efficacy or of low efficacy applied in combination with other intervention strategies.

Figure 5
figure 5

Average overall illness attack rates (%) for modeled interventions. Average overall illness attack rates for the following scenarios: no intervention; pre-vaccination; reactive vaccination with delays in initiation of 0, 30, and 60 days after the intervention trigger of a 0.01% overall illness attack rate; antiviral treatment or household prophylaxis with 10% population coverage (intervention “A”); rolling, as-needed five-day individual school closures and social distancing (20% reduction in workgroup and general community contacts—intervention “S”); antiviral use plus vaccination; school closure, and social distancing plus vaccination; antiviral use, school closure and social distancing, plus vaccination (“A+S”). Vaccination coverage is 35% of the population in Figures 5A and 5B; it is 70% of the population in Figures 5C and 5D. In reactive vaccination scenarios, additional supply-chain disruptions are assumed, such that vaccines are available in three equal batches, spaced two weeks apart, after initiation of vaccination programs. In Figures 5A and 5C, a low-efficacy vaccine is assumed (efficacy against susceptibility, VEs, 0.3; efficacy against infectiousness, VEi, 0.2). In Figures 5B and 5D, a moderate-efficacy vaccine is assumed (VEs, 0.5; VEi, 0.5).

Complete (age-stratified and overall) average illness attack results for all modeled interventions are given in Supplementary Data Table S1. The comparative effectiveness of interventions is similar when age-group-specific results are studied.

Figure 6A illustrates attack rate and total cost combinations for interventions that result in at least a 75% reduction in cost compared to no intervention. The closer to the origin, the more desirable an intervention is in terms of total cost and average illness attack rate. Aside from pre-vaccination strategies, we see that 70% reactive vaccination with a moderate-efficacy vaccine and school closure and social distancing, or even 35% reactive vaccination with a moderate-efficacy vaccine, antiviral use, and school closure, also result in substantial reductions in cost and attack rates. Figure 6B illustrates attack rate and cost results for interventions that result in more-modest 50%–75% reductions in cost compared to no intervention. Once again, several strategies combining vaccination, antiviral use, and school closure/social distancing are competitive with pre-vaccination.

Figure 6
figure 6

Total cost of modeled intervention strategies (US$m) vs. average illness attack rate (%) Figure 6A shows results for interventions with cost reductions of more than 75% compared with no intervention, and Figure 6B shows results for interventions with cost reductions of 50%−75% compared with no intervention. Abbreviations for modeled interventions: PV (pre-vaccination), V (vaccination), L (low-efficacy), M (moderate efficacy), 35 (35% coverage of population), 70 (70% coverage), A (antiviral treatment and household prophylaxis of up to 10% of the population), S (school closure and social distancing). Multiple occurrences of each plotting symbol may occur; occurrences at higher costs and illness attack rates represent interventions with longer supply-chain delays.

Discussion

Previously published research has shown that pre-vaccination of 60%–70% of the population can contain seasonal as well as pandemic influenza, but that delays in vaccination can greatly reduce the effectiveness of the vaccination programs [47]. Our model confirms these results for moderate-efficacy vaccines (Tables 3, 4, and S1). However, vaccination efforts in countries such as the US, Canada, and others began well after the first waves of H1N1 activity, and it is reasonable to believe that the same will be true in future outbreaks of pandemic influenza. In particular, in the event of an outbreak, it will likely take time to achieve high levels of vaccination coverage, and, if past experience with seasonal influenza vaccination campaigns is an indication, it is plausible that only low or moderate coverage will eventually be achieved. The results of our simulation model show that delayed and low-coverage reactive vaccination strategies (with a low-efficacy vaccine, plus limited use of antivirals) will not be enough to mitigate the pandemic or to significantly reduce total costs associated with influenza morbidity and mortality (based on results from Table 3, average illness attack rates are only reduced by 26% and total costs by 13%, compared to no intervention).

According to our model, combining rolling, limited-duration, as-needed closures of individual schools and a practical social distancing policy with 35% reactive low-efficacy vaccination coverage and low-level (10%) antiviral use can reduce illness attack rates by 89% compared to no intervention, as well as total costs by 64%. Similarly, combining interventions in this manner reduces overall attack rates by 99% and costs by 84% when a moderate-efficacy vaccine is available. This strategy remains highly effective even when delays in implementing vaccination of up to 60 days occur. Previously published results have left open the question of how costly interventions involving school closure might be [5]. Our results show that reactive combination strategies that include practical school closure measures, when diligently implemented, can reduce total costs associated with influenza morbidity and mortality substantially.

Our model has several limitations. We do not consider vaccination strategies targeted to high-risk groups, which could reduce costs associated with complications from influenza. We have not modeled co-circulating strains of seasonal and pandemic influenza or possible resistance to antiviral drugs (although, to mitigate this limitation, our model assumes only low coverage with antivirals, as well as interventions without antivirals). As is always the case with simulation models, continuing follow-up analyses are needed, including: (i) sensitivity to model parameters; (ii) sensitivity to model intervention triggers (e.g., overall illness attack rate, numbers of cases detected in schools, etc.); (iii) sensitivity to R0, which can be heterogeneous across cities and countries; and (iv) results for new H1N1 natural history and transmission parameters, and new cost estimates for complications resulting from H1N1 illness, as they become known.

Our model has several strengths. We model a large, realistic, heterogeneous population, base the simulation model on well-studied and documented stochastic simulators, calibrate to actual H1N1 attack rates and most-likely R0 values, and have the ability to model large numbers of scenarios in a relatively short amount of time on a desktop platform. The model also provides cost estimates that are useful for making policy decisions about potentially expensive interventions. In particular, we model and analyze a variety of interventions and combinations of interventions in terms of costs and efficacy. We also take into consideration reactive strategies incorporating supply-chain delays, and we identify strategies that effectively contain outbreaks and costs even in the presence of supply-chain delays, low vaccine efficacy, and low vaccine coverage.

Conclusions

Our model illustrates the epidemiological effectiveness of a combination strategy involving short-term closures of individual schools on an as-needed basis, other practical social distancing activities, reactive vaccination of 35% or more of the population, and limited use of antivirals for treatment and prophylaxis. The model also quantifies the cost savings for this and alternative reactive strategies. Public health authorities should consider placing renewed emphasis on such combination strategies when reacting to possible additional waves of the current pandemic, or to new waves of future pandemics.

Appendix

In this Appendix, we provide details on the simulation model as well as economic cost considerations.

Simulation model

Our simulator is similar to those developed by Longini et al. for high-end computing platforms [7, 16]; our simulator is programmed in C++ and runs on desktop platforms. Population structure and influenza transmission model details are given below.

Population structure

As discussed in the main text, the stochastically generated attributes for each person in our population of 649,565 included: age, household, playgroup or daycare attended (for pre-school children), school attended (for children 5–18 years of age), workgroup (for working adults and working 16–18 year old children), household census tract and workplace census subdivision, community (approximately 2000 people), and neighborhood (approximately 500 people). Thus individuals belong to three or four contact groups. In particular, each individual belongs to a household, neighborhood, and community. In addition, children younger than 16 belong to either a playgroup, daycare, or school, depending on age; most children in age range 16–18 belong to a school or workgroup; and most adults in age range 19–59 belong to a workgroup. Preschool children were categorized as belonging to a playgroup / daycare, each with 50% probability. We separated secondary schools into middle schools and high schools based on grade to allow different contact group sizes and to make our model more representative of mid-sized US cities. The numbers of playgroups, daycares, elementary, middle, and high schools in each community were based on Longini et al. [16], and were combined with the number of individuals in each category in our simulation population to obtain the contact group sizes. The number of working adults (19–59 years old) was based on census data [23]; and the number of working children (16–18 years old) was based on Ontario data on drop-out rates [34] and the employment rate for ages 15–24 [23].

Influenza transmission model

The simulator models influenza transmission over a 180-day period, within the contact groups previously defined. Figure 3 depicts a flowchart of the model. The modeled natural history and simulator dynamics parameters, described below and shown in Figure 3, were based on Longini et al. [7, 19].

To initiate influenza outbreaks, simulations are seeded with approximately 100 randomly selected initial infectives, with all other individuals considered susceptible (state 0). Susceptible people have the opportunity, each day, to become infected in their contact groups. As discussed in the main text, the daily probability of infection for each susceptible person is determined by the number of infectious contacts in his contact groups, and on the per-contact probability of transmission for each type of contact. For example, the probability of a susceptible child who attends daycare being infected on a particular day is:

1 – [Pr(child is not infected in the household)

× Pr(child is not infected in the neighborhood)

× Pr(child is not infected in the community)

× Pr(child is not infected at the daycare center)].

Within each contact group, the probability of infection of a susceptible individual depends on the number of infectious individuals in the group. For example, suppose that k1 children and k2 adults in a household are infectious on a particular day. Then the probability of a susceptible household member being infected in that household on that day is:

1 – [Pr(not infected by a particular infected child in the household)k1

× Pr(not infected by a particular infected adult in the household)k2].

The number of infectious people in the contact groups (e.g., k1 and k2), are random variables that are updated at the beginning of each day.

Age- and contact-group-specific per-contact probabilities of transmission of infection are given in Table 1. The probability that infection is transmitted from an infected person to a susceptible person also depends on whether the infectious person is symptomatic or asymptomatic. Table 1 shows the rates for symptomatic individuals. The transmission rates for asymptomatic individuals are half of those shown in Table 1. These probabilities are based on Longini et al. [7, 16], with adjustments made to calibrate baseline (no intervention) results to age-group-specific illness attack rates and R0 estimates for novel A (H1N1) in Ontario [14, 24, 25]; see Table 2.

Once infected, people enter a 1–3 day latent period (state 1; average length 1.9 days). They are assumed to become infectious on the last day of the latent period, and are half as infectious as they will be after the latent period ends. After the latent period, 67% of infectives become symptomatic (state 2), and 33% are asymptomatic (state 3). These infectious states last between 3 and 6 days. Symptomatic infectives are assumed to be twice as infectious as asymptomatics, and have a chance of withdrawing home during each day of illness (see Figure 3); upon withdrawal, they only make contacts within their household and neighborhood, with transmission probabilities doubled in the household contact group, until they recover. If a school child withdraws home due to illness, one adult in the household also stays home. Each day in states 2 and 3, an infectious person has a chance to exit the state and be removed from the simulation (i.e., to recover or die — state 4). Probabilities for transition into and out of states are given in Figure 3 and are based on Longini et al. [7, 16].

Economic cost calculations

The total cost of each intervention scenario includes the cost of vaccine doses and antiviral courses used, if any; costs associated with parents staying at home with sick children and school teachers, parents, and children staying home due to school closure; costs due to illness-related absence from work; medical costs associated with illness, including outpatient visits, prescription and over-the-counter drugs, and hospitalization; and lost earnings due to death.

We use methods described by Meltzer et al. [29] to quantify most medical and work-loss costs (see also [33]). Table 5 shows the proportions of illnesses assumed to be at high risk for complications among children (0–18 years old), younger adults (19–59 years old) and seniors (over 60). Table 6 shows estimated rates of outpatient visits, hospitalizations, and death used in our calculations for children, adults, and seniors at high risk and not at high risk of complications. We chose the ‘low’ rate estimates presented in Meltzer et al. [29], which we believe to be most consistent with the relatively low R0 (1.4) for our model. Outpatient visit, hospitalization, and death costs are shown in Table 7; cost figures from Meltzer et al. [29] have been inflated using 2008 consumer price and medical price indexes [3032]. All the above costs were combined with age-specific attack rates obtained from our simulation model. In addition, we assume average costs of $25 per vaccine dose or antiviral course used, consistent with previous reports [35]. Table 8 shows other costs associated with vaccination (i.e., the cost of lost time, travel, and side effects). These costs are based on [34], inflated as described above. The vaccination costs are combined with the number of used vaccination doses obtained from our simulation model. We assume that 1% of antiviral users discontinue use due to side effects; medical and other costs associated with these side effects are not included in our model.

Table 5 Proportions of influenza cases at high risk for complications1
Table 6 Outpatient visit, hospitalization, and death rates, by age group and risk status for complications1
Table 7 Frequency and costs (in US$) associated with influenza-related outpatient visits, hospitalizations, and deaths1
Table 8 Costs and impacts of vaccination1

To estimate costs of ill individuals staying home and work-loss associated with parents staying at home with sick children, we multiplied the number of days (obtained from our simulation model) with the inflation-adjusted average value of lost days from Table 7. Similarly, we estimated the average number of teachers at schools and daycares by dividing the total number of such teachers in Hamilton [36] among the schools and daycares in our model. To estimate the cost of lost teacher productivity due to school closures, we multiplied the number of days schools and daycares are closed in our simulation model by the average number of teachers at Hamilton schools and daycares and by the average value of a day of lost work obtained from Table 7.

Table S1 shows age-stratified and overall illness attack rates for all modeled scenarios, along with total cost estimates. Figure 7 depicts the total cost (US$) plotted vs. average overall illness attack rate for each intervention.

Figure 7
figure 7

Total cost of modeled intervention strategies versus the average illness attack rate

References

  1. World Health Organization: World now at the start of 2009 influenza pandemic. (Accessed October 26, 2009.), [http://www.who.int/mediacentre/news/statements/2009/h1n1_pandemic_phase6_20090611/en/index.html]

  2. World Health Organization: Pandemic (H1N1) 2009 - update 71. (Accessed October 26, 2009.), [http://www.who.int/csr/don/2009_10_23/en/index.html]

  3. The Centers for Disease Control and Prevention: 2009 H1N1 Flu Situation Update. (Accessed October 26, 2009.), [http://www.cdc.gov/H1n1flu/update.htm]

  4. Yang Y, Sugimoto JD, Halloran ME, Basta NE, Chao DL, Matrajt L, Potter G, Kenah E, Longini IM: The transmissibility and control of pandemic influenza A (H1N1) virus. Science. 2009, 326: 729-733. 10.1126/science.1177373.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Gojovic MZ, Sander B, Fisman D, Krahn MD, Bauch CT: Modelling mitigation strategies for pandemic (H1N1) 2009. Canadian Medical Association Journal. 2009, 181 (10): 673-680. 10.1503/cmaj.091641.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, Iamsirithaworn S, Burke SD: Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005, 437: 209-214. 10.1038/nature04017.

    Article  CAS  PubMed  Google Scholar 

  7. Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings DAT, Halloran ME: Containing pandemic influenza at the source. Science. 2005, 309: 1083-1087. 10.1126/science.1115717.

    Article  CAS  PubMed  Google Scholar 

  8. The Centers for Disease Control and Prevention: 2009 H1N1 vaccine doses allocated, ordered, and shipped by project area. (Accessed October 26, 2009.), [http://www.cdc.gov/h1n1flu/vaccination/updates/102309.htm]

  9. The Public Health Agency of Canada: H1N1 flu virus in Canada. (Accessed October 26, 2009.), [http://www.phac-aspc.gc.ca/alert-alerte/h1n1/faq_rg_h1n1-eng.php]

  10. The New York Times: A nation battling flu and short vaccine supplies. 2009, (Accessed October 26, 2009.), [http://www.nytimes.com/2009/10/26/health/26flu.html?hpw]

    Google Scholar 

  11. The New York Times: City parents opting out of swine flu vaccine. 2009, (Accessed November 6, 2009.), [http://www.nytimes.com/2009/10/29/nyregion/29vaccine.html?_r=1]

    Google Scholar 

  12. The New York Times: Swine flu shots revive a debate about vaccines. 2009, (Accessed October 21, 2010.), [http://query.nytimes.com/gst/fullpage.html?res=9903EFDA1739F935A25753C1A96F9C8B63]

    Google Scholar 

  13. The Centers for Disease Control and Prevention: 2009 H1N1 vaccination recommendations. (Accessed December 9, 2009.), [http://www.cdc.gov/h1n1flu/vaccination/acip.htm]

  14. Tuite AR, Greer AL, Whelan M, Winter A, Lee B, Yan P, Wu J, Moghadas S, Buckeridge D, Pourbohloul B, Fisman DN: Estimated epidemiologic parameters and morbidity associated with pandemic H1N1 influenza. Canadian Medical Association Journal. DOI:10.1503/cmaj.091807. (Advance access December 3, 2009.)

  15. Halloran ME, Ferguson NM, Eubank S, Longini IM, Cummings DAT, Lewis B, Xu S, Fraser C, Vullikanti A, Germann TC, Wagener D, Beckman R, Kadau K, Barrett C, Macken CA, Burke DS, Cooley P: Modeling targeted layered containment of an influenza pandemic in the United States. Proceedings of the National Academy of Sciences of the USA. 2008, 105 (12): 4639-4644. 10.1073/pnas.0706849105. (Accessed October 22, 2010.), [http://www.pnas.org/cgi/doi/10.1073/pnas.0706849105]

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Longini IM, Halloran ME, Nizam A, Yang Y: Containing pandemic influenza with antiviral agents. American Journal of Epidemiology. 2004, 159: 623-633. 10.1093/aje/kwh092.

    Article  PubMed  Google Scholar 

  17. Aleman DM, Wibisono TG, Schwartz B: Accounting for individual behaviors in a pandemic disease spread model. Proceedings of the 2009 Winter Simulation Conference. Edited by: Rossetti MD, Hill RR, Johansson B, Dunkin A, Ingalls RG. 2009, Piscataway, NJ: Institute for Electrical and Electronics Engineers, 1977-1985.

    Chapter  Google Scholar 

  18. Stroud P, Del Valle S, Sydoriak S, Riese J, Mniszewski S: Spatial dynamics of pandemic influenza in a massive artificial society. Journal of Artificial Societies and Social Simulation. 2007, 10 (4): 9-[http://jasss.soc.surrey.ac.uk/10/4/9.html]

    Google Scholar 

  19. Statistics Canada 2006 Census of Population, CANSIM using CHASS: Table 510036 - Estimates of population, by sex and age group, census divisions and census metropolitan areas, 2006 Census boundaries, annually (Persons unless specified). (Accessed May 6, 2009.), [http://dc1.chass.utoronto.ca/cansimdim/English/]

  20. Statistics Canada 2006 Census of Population, CANSIM using CHASS: Profile Tables/Marital status, families & household living arrangements, and housing. (Accessed May 6, 2009.), [http://dc1.chass.utoronto.ca/census/index.html and http://www.40.statcan.gc.ca/l01/cst01/famil54b-eng.htm]

  21. Statistics Canada 2006 Census of Population, CANSIM using CHASS: 2006 Profile of Census Tracts / Age and sex. (Accessed May 6, 2009.), [http://dc1.chass.utoronto.ca/cansimdim/English/]

  22. Statistics Canada 2006 Census of Population, CANSIM using CHASS: Labour force survey estimates (LFS), by census metropolitan area, sex and age group, annually (Persons unless specified). (Accessed May 6, 2009.), [http://dc1.chass.utoronto.ca/cansimdim/English/]

  23. Statistics Canada 2006 Census of Population, CANSIM using CHASS: 2006 Profile of Census Tracts/Labour market activity. (Accessed May 6, 2009.), [http://dc1.chass.utoronto.ca/cansimdim/English/]

  24. Hancock K, Veguilla V, Lu X, Zhong W, Butler EN, Sun H, Liu F, Dong L, DeVos JR, Gargiullo PM, Brammer TL, Cox NJ, Tumpey TM, Katz M: Cross-reactive antibody responses to the 2009 pandemic H1N1 Influenza Virus. New England Journal of Medicine. 2009, 361: 1945-1952. 10.1056/NEJMoa0906453.

    Article  CAS  PubMed  Google Scholar 

  25. Tuite AR, Greer A, Kwong JC, Fisman DN: Optimal pandemic influenza vaccine allocation strategies for the Canadian population. PLoS ONE. 5 (5): e10520-10.1371/journal.pone.0010520. (Accessed Sept. 12, 2010.), [http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2865540/pdf/pone.0010520.pdf]

  26. Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, Jombart T, Hinsley WR, Grassly NC, Balloux F, Ghani AC, Ferguson NM, Rambaut A, Pybus OG, Lopez-Gatell H, Alpuche-Aranda CM, Chapela IB, Zavala EP, Guevara DME, Checchi F, Garcia E, Hugonnet S, Roth C: The WHO Rapid Pandemic Assessment Collaboration: Pandemic potential of a strain of influenza A (H1N1): Early findings. Science. 2009, 324: 1557-1561. 10.1126/science.1176062.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. New Scientist: Warm weather may not halt swine flu. 2009, (Accessed November 13, 2009.), [http://www.newscientist.com/article/dn17100-warm-weather-may-not-halt-swine-flu.html]

    Google Scholar 

  28. Basta NE, Chao DL, Halloran ME, Matrajt L, Longini IM: Strategies for pandemic and seasonal influenza vaccination of schoolchildren in the United States. American Journal of Epidemiology. 2009, 170: 679-686. 10.1093/aje/kwp237.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Meltzer MI, Cox NJ, Fukuda K: The economic impact of pandemic influenza in the United States: Priorities for intervention. Emerging Infectious Diseases. 1999, 5: 659-671. 10.3201/eid0505.990507.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. InflationData.com: Historical CPI-U data from 1913 to the present. 2009, (Accessed November 29, 2009.), [http://www.inflationdata.com/Inflation/Consumer_Price_Index/HistoricalCPI.aspx]

    Google Scholar 

  31. Bureau of Labor Statistics: Consumer Price Index. (Accessed November 29, 2009.), [http://www.bls.gov/cpi/cpi_dr.htm]

  32. Swivel Preview: Consumer Price Index for Medical Care. (Accessed November 29, 2009.), [http://www.swivel.com/data_columns/spreadsheet/1875739]

  33. Medlock J, Galvani AP: Optimizing influenza vaccine distribution. Science. 2009, 325: 1705-1708. 10.1126/science.1175570.

    Article  CAS  PubMed  Google Scholar 

  34. Bowlby G: Provincial drop-out rates – Trends and consequences. Statistics Canada, (Accessed May 6, 2009.), [http://www.statcan.gc.ca/pub/81-004-x/2005004/8984-eng.htm#b]

  35. Sander B, Nizam A, Garrison L, Postma MJ, Halloran ME, Longini IM: Economic evaluation of influenza pandemic mitigation strategies in the US using a stochastic microsimulation transmission model. Value in Health. 2008, 12: 226-233. 10.1111/j.1524-4733.2008.00437.x.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Statistics Canada: Industry – North American Industry Classification System 2002 (433), Class of Worker (6) and Sex (3) for the Labour Force 15 Years and Over of Canada, Provinces, Territories, Census Metropolitan Areas and Census Agglomerations, 2006 Census – 20% Sample Data. (Accessed 11/29/2009.), [http://www.statcan.gc.ca/bsolc/olc-cel/olc-cel?lang=eng&catno=97-559-X2006009]

Download references

Acknowledgement

We thank the referees and Robert Smith? for their constructive comments.

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/1471-2458/11?issue=S1.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to David Goldsman.

Additional information

Authors' contributions

Study conception and design: AN, SA, DG, KLT

Simulation model development: AN, WC, KLT, SA, MLL, DG, BS, DNF

Analysis and interpretation of simulation results: AN, SA, DG, WC

Drafting of manuscript: SA, AN, DG

All authors read and approved the final manuscript.

Competing interests

DNF has received grant matching funds from Sanofi Pasteur, which manufactures a vaccine for use against influenza A (H1N1)-2009 outside Canada.

Electronic supplementary material

12889_2011_2823_MOESM1_ESM.doc

Additional file 1: Supplementary Data for Reactive Strategies for Containing Developing Outbreaks of Pandemic Influenza(DOC 283 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Andradóttir, S., Chiu, W., Goldsman, D. et al. Reactive strategies for containing developing outbreaks of pandemic influenza. BMC Public Health 11 (Suppl 1), S1 (2011). https://doi.org/10.1186/1471-2458-11-S1-S1

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2458-11-S1-S1

Keywords