- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Influenza Research and Treatment
Volume 2011 (2011), Article ID 579597, 14 pages
Predictive and Reactive Distribution of Vaccines and Antivirals during Cross-Regional Pandemic Outbreaks
Department of Industrial and Management Systems Engineering, University of South Florida, Tampa, FL 33620, USA
Received 14 July 2010; Revised 30 January 2011; Accepted 21 March 2011
Academic Editor: Lance Jennings
Copyright © 2011 Andrés Uribe-Sánchez and Alex Savachkin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
As recently pointed out by the Institute of Medicine, the existing pandemic mitigation models lack the dynamic decision support capability. We develop a large-scale simulation-driven optimization model for generating dynamic predictive distribution of vaccines and antivirals over a network of regional pandemic outbreaks. The model incorporates measures of morbidity, mortality, and social distancing, translated into the cost of lost productivity and medical expenses. The performance of the strategy is compared to that of the reactive myopic policy, using a sample outbreak in Fla, USA, with an affected population of over four millions. The comparison is implemented at different levels of vaccine and antiviral availability and administration capacity. Sensitivity analysis is performed to assess the impact of variability of some critical factors on policy performance. The model is intended to support public health policy making for effective distribution of limited mitigation resources.
As of July 2010, WHO has reported 501 confirmed human cases of avian influenza A/(H5N1) which resulted in 287 deaths worldwide . At the same time, the statistics for the H1N1 2009 outbreak has so far included 214 countries with a total reported number of infections and deaths of 419,289 and 18,239, respectively . Today, an ominous expectation exists that the next pandemic will be triggered by a highly pathogenic virus, to which there is little or no pre-existing immunity in humans .
The nation's ability to mitigate a pandemic influenza depends on the available emergency response resources and infrastructure, and, at present, challenges abound. Predicting the exact virus subtype remains a difficult task, and even when identified, reaching an adequate vaccine supply can currently take up to nine months [4, 5]. Even if the existing vaccines prove to be potent, their availability will be limited by high production and inventory costs [6, 7] and also will be constrained by the supply of antiviral drugs, healthcare providers, hospital beds, medical supplies, and logistics. Hence, pandemic mitigation will have to be done amidst limited availability of resources and supporting infrastructure. This challenge has been acknowledged by WHO  and echoed by the HHS and CDC [8, 9].
The existing models on pandemic influenza (PI) containment and mitigation aims to address various complex aspects of the pandemic evolution process including: (i) the mechanism of disease progression, from the initial contact and infection transmission to the asymptomatic phase, manifestation of symptoms, and the final health outcome [10–12], (ii) the population dynamics, including individual susceptibility [13, 14] and transmissibility [10, 15–17], and behavioral factors affecting infection generation and effectiveness of interventions [18–20], (iii) the impact of pharmaceutical and nonpharmaceutical measures, including vaccination [21–23], antiviral therapy [24–26], social distancing [27–31] and travel restrictions, and the use of low-cost measures, such as face masks and hand washing [26, 32–34].
Recently, the modeling efforts have focused on combining pharmaceutical and nonpharmaceutical interventions in search for synergistic strategies, aimed at better resource utilization. Most of such approaches attempt implementing a form of social distancing followed by application of pharmaceutical measures. For significant contributions in this area see [33, 35–41]. One of the most notable among these efforts is a 2006-07 initiative by MIDAS , which cross-examined independent simulation models of PI spread in rural areas of Asia [43, 44], USA and UK [45, 46], and the city of Chicago , respectively. MIDAS cross-validated the models by simulating the city of Chicago, with 8.6M inhabitants and implementing a targeted layered containment [48, 49]. The research findings of MIDAS and some other groups [12, 33] were used in a recent “Modeling Community Containment for Pandemic Influenza” report by IOM, to formulate a set of recommendations for PI mitigation . These findings were also used in a pandemic preparedness guidance developed by CDC .
At the same time, The IOM report  points out several limitations of the MIDAS models, observing that “because of the significant constraints placed on the models … the scope of models should be expanded.” The IOM recommends “to adapt or develop decision-aid models that can … provide real-time feedback … and include the costs and benefits of intervention strategies.” Our literature review yields a similar observation that most existing approaches focus on assessment of a priori defined strategies, and virtually none of the models are capable of “learning,” that is, adapting to changes in the pandemic progress, or even predicting them, to generate dynamic strategies. Such a strategy has the advantage of being developed dynamically, as the pandemic spreads, by selecting a mix of available mitigation options at each decision epoch, based on both the present state of the pandemic and its predicted evolution.
In an attempt to address the IOM recommendations, we present a simulation optimization model for developing predictive resource distribution over a network of regional outbreaks. The underlying simulation model mimics the disease and population dynamics of each of the affected regions (Sections 2.1 and 2.2). As the pandemic spreads from region to region, the optimization model distributes mitigation resources, including stockpiles of vaccines and antiviral and administration capacities (Section 2.3). The model seeks to minimize the impact of ongoing outbreaks and the expected impact of potential outbreaks, using measures of morbidity, mortality, and social distancing, translated into the cost of lost productivity and medical expenses. The methodology is calibrated and implemented on a sample outbreak in Fla, USA with over 4M inhabitants (Section 3). The strategy is compared to the reactive myopic policy, which allocates resources from one actual outbreak region to the next, each time trying to cover the entire regional population at risk, regardless of the resource availability. The comparison is done at different levels of vaccine and antiviral availability and administration capacity. We also present a sensitivity analysis for assessing the impact of variability of some critical factors, including: (i) antiviral efficacy, (ii) social distancing conformance, and (iii) CDC response delay.
The objective of our methodology is to generate a progressive allocation of the total resource availability over a network of regional outbreaks. The methodology incorporates (i) a cross-regional simulation model, (ii) a set of single-region simulation models, and (iii) an embedded optimization model.
We consider a network of regions with each of which classified as either unaffected, ongoing outbreak, or contained outbreak (Figure 1). The cross-regional simulation model connects the regions by air and land travel. The single-region simulation models mimic the population and disease dynamics of each ongoing region, impacted by intervention measures. The pandemic can spread from ongoing to unaffected regions by infectious travelers who pass through regional border control. At every new regional outbreak epoch, the optimization model allocates available resources to the new outbreak region (actual distribution) and unaffected regions (virtual distribution). Daily statistics is collected for each ongoing region, including the number of infected, deceased, and quarantined cases, for different age groups. As a regional outbreak is contained, its societal and economic costs are calculated.
2.1. Cross-Regional Simulation Model
A schematic of the cross-regional simulation model is shown in Figure 2. The model is initialized by creating population entities and mixing groups, for each region. A pandemic is started by an infectious case injected into a randomly chosen region. The details of the resulting regional contact dynamics and infection transmission are given in Section 2.2. As the infected cases start seeking medical help, a new regional outbreak is detected. A resource distribution is then determined and returned to the single-region model. The outbreak can spread to unaffected regions as some infectious travelers pass undetected through the border control. By tracing these travelers, the model determines which of the unaffected regions, if any, become new outbreaks. The model also determines if any ongoing outbreaks have been contained. The simulation stops when all outbreaks are contained.
2.2. Single-Region Simulation Model
The single-region model subsumes the following components (see Figure 3): (i) population dynamics (mixing groups and schedules), (ii) contact and infection process, (iii) disease natural history, and (iv) mitigation strategies, including social distancing, vaccination, and antiviral application. The model collects detailed statistics, including number of infected, recovered, deceased, and quarantined cases, for different age groups. For a contained outbreak, its societal and economic costs are calculated. The societal cost includes the cost of lost lifetime productivity of the deceased; the economic cost includes the cost of medical expenses of the recovered and deceased and the cost of lost productivity of the quarantined .
2.2.1. Mixing Groups and Schedules
Each region is modeled as a set of population centers formed by mixing groups or places where individuals come into contact with each other during the course of their social interaction. Examples of mixing groups include households, offices, schools, universities, shopping centers, entertainment centers, and so forth, . Each individual is assigned a set of attributes such as age, gender, parenthood, workplace, infection susceptibility, and probability of travel, among others. Each person is also assigned time-discrete (e.g., = 1 hour) weekday and weekend schedules, which depend on: (i) person's age, parenthood, and employment status, (ii) disease status, (iii) travel status, and (iv) person's compliance to social distancing decrees . As their schedules advance, the individuals circulate throughout the mixing groups and come into contact with each other (see Section 2.2.2).
It is assumed that at any point of time, an individual belongs to one of the following compartments (see Figure 4): susceptible, contacted (by an infectious individual), infected (asymptomatic or symptomatic), and recovered/deceased. In what follows, we present the infection transmission and disease natural history model, which delineates the transitions between the above compartments.
2.2.2. Contact and Infection Process
Infection transmission occurs during contact events between susceptible and infectious cases, which take place in the mixing groups. At the beginning of every period (e.g., one hour), for each mixing group , the simulation tracks the total number of infectious cases, , present in the group. It is assumed that each infectious case generates per unit of time new contacts , chosen randomly (uniformly) from the pool of susceptibles present in the group. We also assume the following: (i) during period, a susceptible may come into contact with at most one infectious case and (ii) each contact exposure lasts units of time. Once a susceptible has started her contact exposure at time , she will develop infection at time with a certain probability that is calculated as shown below.
Let be a nonnegative continuous random variable that represents the duration of contact exposure, starting at time , required for susceptible to become infected. We assume that is distributed exponentially with mean , where represents the instantaneous force of infection applied to susceptible at time [55–57]. The probability that susceptible , whose contact exposure has started at time , will develop infection at time is then given as
2.2.3. Disease Natural History
A schematic of the disease natural history is shown in Figure 5. During the incubation phase, the infected case stays asymptomatic. At the end of the latency phase, she enters the infectious phase [44, 46, 48]. She becomes symptomatic at the end of the incubation period. At the end of the infectious phase, she enters the period leading to a health outcome, which culminates in her recovery or death.
Mortality for influenza-like diseases is a complex process affected by many factors and variables, most of which have limited accurate data support available from past pandemics. Furthermore, the time of death can sometimes be weeks following the disease episode (which is often attributable to pneumonia-related complications ). Because of the uncertainty underlying the mortality process, we adopted an age-based form of the mortality probability of infected , as follows: where is the age-dependent base mortality probability of infected , is her status of antiviral therapy (0 or 1), and is the antiviral efficacy measured as the relative decrease in the base probability . We assume that a recovered case develops full immunity but continues circulating in the region.
2.2.4. Mitigation Strategies
Mitigation options include pharmaceutical and nonpharmaceutical interventions. Mitigation is initiated upon detection of a critical number of confirmed infected cases , which triggers resource distribution and deployment. The model incorporates a certain delay for deploying field responders.
Pharmaceutical intervention (PHI) includes vaccination and antiviral application. Vaccination is targeted at individuals at risk  to reduce their infection susceptibility. The vaccine takes a certain period to become effective . Vaccination is constrained by the allocated stockpile and administration capacity, measured in terms of the immunizer-hours. We assume that as some symptomatic cases seek medical help [62, 63], those at risk of them will receive an antiviral. The process is constrained by the allocated stockpile and administration capacity, measured in terms of the number of certified providers.
Both vaccination and antiviral application are affected by a number of sociobehavioral factors, including conformance of the target population, degree of risk perception, and compliance of healthcare personnel [64–66]. The conformance level of the population at risk can be affected, among other factors, by the demographics and income level [67–71] as well as by the quality of public information available . The degree of risk perception can be influenced by the negative experience developed during previous pharmaceutical campaigns [72, 73], as well as by public fear and rumors [74, 75].
Nonpharmaceutical intervention (NPI) includes social distancing and travel restrictions. We adopted a CDC guidance , which establishes five categories of pandemic severity and recommends quarantine and closure options according to the category. The categories are determined based on the value of the case fatality ratio (CFR), the proportion of fatalities in the total infected population. For CFR values lower than (Category 1), voluntary at-home isolation of infected cases is implemented. For CFR values between and (Categories 2 and 3), in addition to at-home isolation, the following measures are recommended: (i) voluntary quarantine of household members of infected cases and (ii) child and adult social distancing. For CFR values exceeding (Categories 4 and 5), all the above measures are implemented. As the effectiveness of social distancing is affected by some of the behavioral factors listed above , we assume a certain social distancing conformance level. Travel restrictions considered in the model included regional air and land border control for infected travelers.
2.3. Optimization Model
As presented in Figure 2, the optimization model is invoked at the beginning of every th new regional outbreak epoch (), starting from the initial outbreak region (). The objective of the model is to allocate some of the available mitigation resources to the new outbreak region (actual distribution) while reserving the rest of the quantities for potential outbreak regions (virtual distribution). By doing so, the model seeks to progressively minimize the impact of ongoing outbreaks and the expected impact of potential outbreaks, spreading from the ongoing locations. Mitigation resources can include stockpiles of vaccines and antivirals, administration capacity, hospital beds, medical supplies, and social distancing enforcement resources, among others. The predictive mechanism of the optimization model is based on a set of regression equations obtained using single-region simulation models. In what follows, we present the construction of the optimization model and explain the solution algorithm for the overall simulation-based optimization methodology.
We introduce the following general terminology and notation:: Set of all network regions: Set of regions in which pandemic is contained at the th outbreak epoch (): Set of ongoing regions at the th outbreak epoch: Set of unaffected regions at the th outbreak epoch: Set of available types of mitigation resources (: Amount of resource allocated to region : Available amount of resource at the th outbreak epoch: Set of age groups.
The optimization criterion (objective function) of the model incorporates measures of expected societal and economic costs of the pandemic: the societal cost includes the cost of lost lifetime productivity of the deceased; the economic cost includes the cost of medical expenses of the recovered and deceased and the cost of lost productivity of the quarantined. To compute these costs, the following impact measures of morbidity, mortality, and quarantine are used, for each region :: Total number of infected cases in age group who seek medical assistance: Total number of infected cases in age group who do not seek medical assistance: Total number of deceased cases in age group : Total number of person-days of cases in age group who comply with quarantine.
To estimate these measures, we use the following regression models obtained using a single-region simulation of each region : where denotes the regression coefficient associated with resource and is the regression coefficient for the interaction between resources and . Similar models are used for , , and .
The above relationships between the impact measures and the resource distributions ought to be determined a priori of implementing a cross-regional scenario (see Section 3). Here, we consider each region as the initial outbreak region. We assume, however, that as the pandemic evolves, the disease infectivity will naturally subside. Hence, the regression equations need to be re-estimated at every new outbreak epoch, for each region , using the single-region simulation models, where each simulation must be initialized to the current outbreak status in region in the cross-regional simulation. As an alternative to using a computationally burdensome approach of re-estimating the regression equations, a modeler may choose to use a certain decay factor  to adjust the estimates of the regional impact measures at every th outbreak epoch, in the following way:
In addition, we use the following regression model to estimate the probability of pandemic spread from affected region to unaffected region , as a function of resources allocated to region , which, in turn, impact the number of outgoing infectious travelers from the region: where denotes the regression coefficient associated with resource , is the regression coefficient associated with interaction between resources and , and represents the intercept. Consequently, the total outbreak probability for unaffected region can be found as . As in the case of the impact measures, the estimates of the regional outbreak probabilities need to be progressively re-estimated or adjusted using a scheme similar to (4), as follows:
Finally, we calculate the total cost of an outbreak in region at the th decision epoch as follows: where is total medical cost of an infected case in age group over his/her disease period, is total cost of lost wages of an infected case in age group over his/her disease period, is cost of lost lifetime wages of a deceased case in age group , and is daily cost of lost wages of a non-infected case in age group who complies with quarantine.
The optimization model has the following form. The first term of the objective function represents the total cost of the new outbreak , estimated at the th outbreak epoch, based on the actual resource distribution (see (7)). The second term represents the total expected cost of outbreaks in currently unaffected regions, based on the virtual distributions (7) and the regional outbreak probabilities (6). The set of constraints assures that for each resource , the total quantity allocated (current and virtual, both nonnegative) does not exceed the total resource availability at the th decision epoch. Note that both the objective function and the availability constraints are nonlinear in the decision variables.
2.4. Solution Algorithm
The solution algorithm for our dynamic predictive simulation optimization (DPO) model is given below.(1) Estimate regression equations for each region using the single-region simulation model.(2) Begin the cross-regional simulation model. (3) Initialize the sets of regions: , , . (4) Select randomly the initial outbreak region . Set . (5) Update sets of regions: and . (6) Solve the resource distribution model for region . Update the total resource availabilities. (7) If , do step 8. Else, do step 10.(8)(a) For each ongoing region, implement a next day run of its single-region simulation. (b) Check the containment status of each ongoing region. Update sets and , if needed. (c) For each unaffected region, calculate its outbreak probability. (d) Based on the outbreak probability values, determine if there is a new outbreak region(s) . If there is no new outbreak(s), go to step 7. Otherwise, go to step 9. (9) For each new outbreak region , (a) Increment . (b) Update sets and . (c) Re-estimate regression equations for each region using the single-region simulations, where each simulation is initialized to the current outbreak status in the region (alternatively, use (4) and (6)). (d) Solve the resource distribution model for region . (e) Update the total resource availabilities. (10) Calculate the total cost for each contained region and update the overall pandemic cost.
3. Testbed Illustration
To illustrate the use of our methodology, we present a sample H5N1 outbreak scenario including four counties in Fla, USA: Hillsborough, Miami Dade, Duval, and Leon, with populations of 1.0, 2.2, 0.8, and 0.25 million people, respectively. A basic unit of time for population and disease dynamics models was taken to be hour. Regional simulations were run for a period (up to 180 days) until the daily infection rate approached near zero (see Section 3.3). Below, we present the details on selecting model parameter values. Most of the testbed data can be found in the supplement .
3.1. Parameter Values for Population and Disease Dynamics Models
Each infected person was assigned a daily travel probability of 0.24% , of which was by air and by land. The probabilities of travel among the four regions were calculated using traffic volume data [80–83], see Table 1. Infection detection probabilities for border control for symptomatic cases were assumed to be 95% and 90%, for air and land, respectively .
The instantaneous force of infection applied to contact at time ((1), ) was modeled as where is the age-dependent base instantaneous infection probability of contact , is her status of vaccination at time (0 or 1), and is the vaccine efficacy, measured as the reduction in the base instantaneous infection probability (achieved after 10 days ).
The values of age-dependent base instantaneous infection probabilities were adopted from  (see Table 2). The disease natural history included a latent period of 29 hours (1.21 days), an incubation period of 46 hours (1.92 days), an infectiousness period from 29 to 127 hours (1.21 to 5.29 days), and a period leading to health outcome from 127 to 240 hours (5.29 to 10 days) .
Base mortality probabilities ( in (2)) were found using the statistics recommended by the Working Group on Pandemic Preparedness and Influenza Response . This data shows the percentage of mortality for age-based high-risk cases (HRC) (Table 3, columns 1–3). Mortality probabilities (column 4) were estimated under the assumption that high-risk cases are expected to account for of the total number of fatalities, for each age group .
3.2. Calibration of the Single-Region Models
Single-region simulation models were calibrated using two common measures of pandemic severity [35, 45, 46]: the basic reproduction number () and the infection attack rate (IAR). is defined as the average number of secondary infections produced by a typical infected case in a totally susceptible population. IAR is defined as the ratio of the total number of infections over the pandemic period to the size of the initial susceptible population. To determine , all infected cases inside the simulation were classified by generation of infection, as in [33, 43]. The value of was calculated as the average reproduction number of a typical generation in the early stage of the pandemic, with no interventions implemented (the baseline scenario) . Historically, values for PI ranged between 1.4 and 3.9 [37, 43]. To attain similar values, we calibrated the hourly contact rates of mixing groups  (original rates were adopted from ). For the four regions, the average baseline value of was 2.54, which represented a high transmissibility scenario. The values of regional baseline IAR averaged 0.538.
3.3. Parameters of Mitigation
The vaccination risk group included healthcare providers , and individuals younger than 5 years (excluding younger than 12 months old) and older than 65 years . The risk group for antiviral included symptomatic individuals below 15 years and above 55 years [60, 86]. The efficacy levels for the vaccine ( in (9)) and antiviral ( in (2)) were assumed to be [44, 87] and , respectively. We did not consider the use of antiviral for a mass prophylactic reduction of infection susceptibility due to the limited antiviral availability  and the risk of emergence of antiviral resistant transmissible virus strains . We assumed a target population conformance for both vaccination and antiviral treatment . The immunity development period for the vaccine was taken as 10 days .
A version of the CDC guidance for quarantine and isolation for Category 5 was implemented (Section 2.2.4, ). Once the reported CFR value had reached , the following policy was declared and remained in effect for 14 days : (i) individuals below a certain age (22 years) stayed at home during the entire policy duration, (ii) of the remaining population, a certain proportion  stayed at home and was allowed a one-hour leave, every three days, to buy essential supplies, and (iii) the remaining noncompliant proportion followed a regular schedule. All testbed scenarios considered the quarantine conformance level equal to .
An outbreak was considered contained, if the daily infection rate did not exceed five cases, for seven consecutive days. Once contained, a region was simulated for an additional 10 days for accurate estimation of the pandemic statistics. A 25 statistical design of experiment  was used to estimate the regression coefficient values of the significant decision factors and their interactions (see Section 2.3; the values of adjusted ranged from to ).
The simulation code was developed using C++. The running time for a cross-regional simulation replicate involving over four million inhabitants was between 17 and 26 minutes (depending on the initial outbreak region, with a total of 150 replicates) on a Pentium 3.40 GHz with 4.0 GB of RAM.
3.4. Comparison of DPO and Myopic Strategies
The performance of the DPO and myopic policies is compared at different levels of resource availability.
Table 4 summarizes the total vaccine and antiviral requirements for each region, based on the composition of the regional risk groups (see Section 3.3). Table 5 shows the per capita costs of lost productivity and medical expenses, which were adopted from  and adjusted for inflation for the year of 2010 .
Comparison of the two strategies is done at the levels of 20%, 50%, and 80% of the total resource requirement shown in Table 4. Figures 6(a) and 6(b) show the policy comparison in the form of the confidence intervals (CI) for the average number of infected and deceased, respectively.
Figure 7 also shows the policy comparison using the CI for the average total pandemic cost, calculated using the pandemic statistics, and the per capita costs from Table 5. For illustrative purposes, we also show the average number of regional outbreaks, for each policy, at different levels of resource availability, in the testbed scenario involving four regions, with the Hillsborough as the initial outbreak region (Table 6).
It can be observed that the values of all impact measures exhibit a downward trend, for both DPO and myopic policies, as the total resource availability increases from to .
An increased total resource availability not only helps alleviating the pandemic impact inside the ongoing regions but also reduces the probability of spread to the unaffected regions. For both policies, as the total resource availability approaches the total resource requirement (starting from approximately 60%), the impact numbers show a converging behavior, whereby the marginal utility of additional resource availability diminishes. This behavior can be explained by noting that the total resource requirements were determined assuming the worst case scenario when all (four) regions would be affected and ought to provided with enough resources to cover their respective regional populations at risk. It can also be seen that on average, the DPO policy outperforms the myopic approach at all levels, which can attest to a more efficient resource utilization achieved by the DPO policy (see also Table 6). The difference in the policy performance is particularly noticeable at the lower levels of resource availability, and it gradually diminishes, as the resource availability increases and becomes closer to be sufficient to cover the entire populations at risk in all regions. It can also be noted that the variability in the performance of the DPO strategy is generally smaller than that of the myopic policy. In general, for both strategies, the performance variability decreases with higher availability of resources.
3.5. Sensitivity Analysis
In this section, we assess the marginal impact of variability of some of the critical factors. The impact was measured separately by the change in the total pandemic cost and the number of deaths (averaged over multiple replicates), resulting from a unit change in a decision factor value, one factor at a time. Factors under consideration included: (i) antiviral efficacy, (ii) social distancing conformance, and (iii) CDC response delay. We have used all four regions, separately, as initial outbreak regions for each type of sensitivity analysis. The results (patterns) were rather similar. Due to limited space, we have opted to show the results for only one initial region, chosen arbitrarily, for each of the three types of sensitivity studies. While Duval County was selected as the initial outbreak region to show the sensitivity results on antiviral efficacy, Hillsborough and Miami Dade were used as the initial regions to show the results on, respectively, social distancing conformance and CDC response delay.
3.5.1. Antiviral Efficacy
Figure 8 depicts the sensitivity of the average total cost and average total deaths to antiviral efficacy values between and . As expected, for both policies, the curves for the average number of deaths exhibit a decreasing trend which is almost linear for the values of between 0% and 40%. As the value of approaches 70%, the curves start exhibit a converging behavior. The curves for the average total pandemic cost exhibit a similar pattern for both policies.
It can be noted that the performance of both policies is somewhat identical for low antiviral efficacy (between and ). However, the performance of the DPO policy improves consistently as increases which can be attributed to a more discretionary allocation of the antiviral stockpile by the DPO policy.
3.5.2. Social Distancing Conformance
Reduction of the contact intensity through quarantine and social distancing has proven to be one of the most effective containment measures, especially in the early stages of the pandemic [27, 30, 31, 41].
Figure 9 shows the sensitivity of the average total cost and average total deaths to the social distancing conformance ranging between and . We observed that for both impact measures, the DPO policy demonstrated a better performance with the difference ranging from $3B to $26B in the total cost and from 1,400 to 20,000 in the number of fatalities. The biggest difference in performance was achieved at the lower-to-medium levels of conformance (between and ). As the conformance level approached , the dominating impact of social distancing masked the effect of better utilization of vaccines and antivirals achieved by the DPO strategy.
3.5.3. CDC Response Delay
The CDC response delay corresponds to the interval of time from the moment an outbreak is detected to a complete deployment of mitigation resources. Depending on the disease infectivity, CDC response delay may represent one of the most critical factors in the mitigation process.
Figure 10 shows how the performance of both policies was significantly impacted by this factor. The DPO policy showed a uniformly better performance with the difference ranging between $3B to $4B in the average total cost, and between 800 to 1,800 in the average number of mortalities, over the range (24–72 hours) of the response delay that we examined. For both policies, there was no significant difference when the delay was between 24 and 48 hours. However, for the delay values exceeding 48 hours, the average number of deaths and total cost increased at a high rate.
As recently pointed by the IOM, the existing models for PI mitigation fall short of providing dynamic decision support which would incorporate “the costs and benefits of intervention” . In this paper, we present a large-scale simulation optimization model which is attempted at filling this gap.
The model supports dynamic predictive resource distribution over a network of regions exposed to the pandemic. The model aims to balance both the ongoing and potential outbreak impact, which is measured in terms of morbidity, mortality, and social distancing, translated into the cost of lost productivity and medical expenses. The model was calibrated using historic pandemic data and compared to the myopic strategy, using a sample outbreak in Fla, USA, with over 4 million inhabitants.
Summary of the main results. In the testbed scenario, for both strategies, the marginal utility of additional resource availability was found to be diminishing, as the total resource availability approached the total requirement.
In the testbed scenario, the DPO strategy on average outperformed the myopic policy. As opposed to the DPO strategy, the myopic policy is reactive, rather than predictive, as it allocates resources regardless of the remaining availability and the overall cross-regional pandemic status. In contrast, the DPO model distributes resources trying to balance the impact of actual outbreaks and the expected impact of potential outbreaks. It does so by exploiting region-specific effectiveness of mitigation resources and dynamic reassessment of pandemic spread probabilities, using a set of regression submodels. Hence, we believe that in scenarios involving regions with a more heterogeneous demographics, the DPO policy will likely to perform even better and with less variability than the myopic strategy. We also note that the difference in the model performance was particularly noticeable at lower levels of resource availability, which is in accordance with a higher marginal utility of additional availability at that levels. We thus believe that the DPO model can be particularly useful in scenarios with very limited resources.
Contributions of the paper. The simulation optimization methodology presented in this paper is one of the first attempts to offer dynamic predictive decision support for pandemic mitigation, which incorporates measures of societal and economic costs. Our comparison study of the DPO versus myopic cross-regional resource distribution is also novel. Additionally, our simulation model represents one of the first of its kind in considering a broader range of social behavioral aspects, including vaccination and antiviral treatment conformance. The simulation features a flexible design which can be particularized to a broader range of PHI and NPI and even more granular mixing groups.
We also developed a decision-aid simulator which is made available to the general public through our web site at http://imse.eng.usf.edu/pandemics.aspx. The tool is intended to assist public health decision makers in implementing what-if analysis for assessment of mitigation options and development of policy guidelines. Examples of such guidelines include vaccine and antiviral risk groups, social distancing policies (e.g., thresholds for declaration/lifting and closure options), and travel restrictions.
Limitations of the model. Lack of reliable data prevented us from considering geo-spatial aspects of mixing group formation. We also did not consider the impact of public education and the use of personal protective measures (e.g., face masks) on transmission, again due to a lack of effectiveness data . We did not study the marginal effectiveness of individual resources due to a considerable uncertainty about the transmissibility of an emerging pandemic virus and efficacy of vaccine and antiviral. For the same reason, the vaccine and antiviral risk groups considered in the testbed can be adjusted, as different prioritization schemes have been suggested. The form of social distancing implemented in the testbed can also be modified as a variety of schemes can be found in the literature, including those based on geographical and social targeting. Effectiveness of these approaches is substantially influenced by the compliance factor, for which limited accurate data support exists. It will thus be vital to gather the most detailed data on the epidemiology of a new virus and the population dynamics early in the evolution of a pandemic, and expeditiously analyze the data to adjust the interventions accordingly.
The authors would like to acknowledge with thanks the many helpful suggestions made by Professor Yiliang Zhu, Department of Epidemiology and Biostatistics at the University of South Florida, Tampa, Fla, USA.
- World Health Organization (WHO), “Cumulative number of confirmed human cases of Avian Influenza A(H5N1) reported to WHO,” 2010, http://www.who.int/csr/disease/avian_Influenza/country/cases_table_2010_07_05/en/index.html.
- World Health Organization, “Pandemic (H1N1),” 2009, http://www.who.int/csr/don/2010_07_02/en/index.html.
- A. Schoenstadt, “Spanish Flu,” 2009, http://flu.emedtv.com/spanish-flu/spanish-flu.html.
- D. S. Fedson, “Pandemic Influenza and the global vaccine supply,” Clinical Infectious Diseases, vol. 36, no. 12, pp. 1552–1561, 2003.
- J. Aunins, A. Lee, and D. Volkin, Vaccine Production, CRC Press, Boca Raton, Fla, USA, 2nd edition, 2000.
- World Health Organization, “Pandemic (h1n1) vaccine deployment,” 2009, http://www.who.int/csr/disease/swineflu/vaccines/h1n1_vaccination_deployment_update_20091217.pdf.
- WHO Global Influenza Programme, “Pandemic Influenza preparedness and response,” Tech. Rep., World Health Organization, Lyon, France, 2009.
- Centers for Disease Control and Prevention (CDC), “Preparing for pandemic Influenza,” 2007, http://www.cdc.gov/flu/pandemic/preparedness.htm.
- U.S. Department of Health & Human Services, “HHS pandemic influenza plan,” 2007, http://www.hhs.gov/pandemicflu/plan/.
- A. Handel, I. Longini, and R. Antia, “Towards a quantitative understanding of the within-host dynamics of Influenza A infections,” Epidemics, vol. 1, no. 3, pp. 185–195, 2009.
- B. Pourbohloul, A. Ahued, B. Davoudi et al., “Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America,” Influenza and other Respiratory Viruses, vol. 3, no. 5, pp. 215–222, 2009.
- M. P. Atkinson and L. M. Wein, “Quantifying the routes of transmission for pandemic Influenza,” Bulletin of Mathematical Biology, vol. 70, no. 3, pp. 820–867, 2008.
- V. E. Pitzer, S. J. Olsen, C. T. Bergstrom, S. F. Dowell, and M. Lipsitch, “Little evidence for genetic susceptibility to Influenza A (H5N1) from family clustering data,” Emerging Infectious Diseases, vol. 13, no. 7, pp. 1074–1076, 2007.
- V. E. Pitzer, G. M. Leung, and M. Lipsitch, “Estimating variability in the transmission of severe acute respiratory syndrome to household contacts in Hong Kong, China,” American Journal of Epidemiology, vol. 166, no. 3, pp. 355–363, 2007.
- Y. Yang, J. D. Sugimoto, M. Elizabeth Halloran et al., “The transmissibility and control of pandemic Influenza a (H1N1) virus,” Science, vol. 326, no. 5953, pp. 729–733, 2009.
- Y. Yang, M. E. Halloran, J. D. Sugimoto, and I. M. Longini, “Detecting human-to-human transmission of avian Influenza A (H5N1),” Emerging Infectious Diseases, vol. 13, no. 9, pp. 1348–1353, 2007.
- S. Cauchemez, F. Carrat, C. Viboud, A. J. Valleron, and P. Y. Boëlle, “A Bayesian MCMC approach to study transmission of Influenza: application to household longitudinal data,” Statistics in Medicine, vol. 23, no. 22, pp. 3469–3487, 2004.
- V. Colizza, A. Barrat, M. Barthélemy, and A. Vespignani, “The role of the airline transportation network in the prediction and predictability of global epidemics,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 7, pp. 2015–2020, 2006.
- J. M. Epstein, J. Parker, D. Cummings, and R. A. Hammond, “Coupled contagion dynamics of fear and disease: mathematical and computational explorations,” PLoS ONE, vol. 3, no. 12, Article ID e3955, 2008.
- M. E. Halloran, “Invited commentary: challenges of using contact data to understand acute respiratory disease transmission,” American Journal of Epidemiology, vol. 164, no. 10, pp. 945–946, 2006.
- F. G. Ball and O. D. Lyne, “Optimal vaccination policies for stochastic epidemics among a population of households,” Mathematical Biosciences, vol. 177-178, pp. 333–354, 2002.
- N. G. Becker and D. N. Starczak, “Optimal vaccination strategies for a community of households,” Mathematical Biosciences, vol. 139, no. 2, pp. 117–132, 1997.
- F. Carrat, A. Lavenu, S. Cauchemez, and S. Deleguer, “Repeated Influenza vaccination of healthy children and adults:borrow now, pay later?” Epidemiology and Infection, vol. 134, no. 1, pp. 63–70, 2005.
- N. M. Ferguson, S. Mallett, H. Jackson, N. Roberts, and P. Ward, “A population-dynamic model for evaluating the potential spread of drug-resistant Influenza virus infections during community-based use of antivirals,” Journal of Antimicrobial Chemotherapy, vol. 51, no. 4, pp. 977–990, 2003.
- P. Y. Lee, D. B. Matchar, D. A. Clements, J. Huber, J. D. Hamilton, and E. D. Peterson, “Economic analysis of Influenza vaccination and antiviral treatment for healthy working adults,” Annals of Internal Medicine, vol. 137, no. 4, pp. 225–231, 2002.
- M. Lipsitch, T. Cohen, M. Murray, and B. R. Levin, “Antiviral resistance and the control of pandemic Influenza,” PLoS Medicine, vol. 4, no. 1, pp. 111–121, 2007.
- G. J. Milne, J. K. Kelso, H. A. Kelly, S. T. Huband, and J. McVernon, “A small community model for the transmission of infectious diseases: comparison of School closure as an intervention in individual-based models of an Influenza pandemic,” PLoS ONE, vol. 3, no. 12, Article ID e4005, 2008.
- H. Yasuda and K. Suzuki, “Measures against transmission of pandemic H1N1 Influenza in Japan in 2009: simulation model,” Euro Surveillance : Bulletin Européen sur les Maladies Transmissibles = European Communicable Disease Bulletin, vol. 14, no. 44, 2009.
- G. Miller, S. Randolph, and J. E. Patterson, “Responding to simulated pandemic Influenza in San Antonio, Texas,” Infection Control and Hospital Epidemiology, vol. 29, no. 4, pp. 320–326, 2008.
- J. K. Kelso, G. J. Milne, and H. Kelly, “Simulation suggests that rapid activation of social distancing can arrest epidemic development due to a novel strain of Influenza,” BMC Public Health, vol. 9, article 117, 2009.
- N. Halder, J. Kelso, and G. Milne, “Analysis of the effectiveness of interventions used during the 2009 A/H1N1 Influenza pandemic,” BioMed Central, vol. 10, article 168, 2010.
- D. O. Scharfstein, M. E. Halloran, H. Chu, and M. J. Daniels, “On estimation of vaccine efficacy using validation samples with selection bias,” Biostatistics, vol. 7, no. 4, pp. 615–629, 2006.
- R. J. Glass, L. M. Glass, W. E. Beyeler, and H. J. Min, “Targeted social distancing design for pandemic Influenza,” Emerging Infectious Diseases, vol. 12, no. 11, pp. 1671–1681, 2006.
- K. R. Nigmatulina and R. C. Larson, “Living with Influenza: impacts of government imposed and voluntarily selected interventions,” European Journal of Operational Research, vol. 195, no. 2, pp. 613–627, 2009.
- I. M. Longini, M. E. Halloran, A. Nizam, and Y. Yang, “Containing pandemic Influenza with antiviral agents,” American Journal of Epidemiology, vol. 159, no. 7, pp. 623–633, 2004.
- J. T. Wu, S. Riley, C. Fraser, and G. M. Leung, “Reducing the impact of the next Influenza pandemic using household-based public health interventions,” PLoS Medicine, vol. 3, no. 9, pp. 1532–1540, 2006.
- C. E. Mills, J. M. Robins, and M. Lipsitch, “Transmissibility of 1918 pandemic Influenza,” Nature, vol. 432, no. 7019, pp. 904–906, 2004.
- R. Patel, I. M. Longini, and M. E. Halloran, “Finding optimal vaccination strategies for pandemic Influenza using genetic algorithms,” Journal of Theoretical Biology, vol. 234, no. 2, pp. 201–212, 2005.
- V. Colizza, A. Barrat, M. Barthelemy, A. J. Valleron, and A. Vespignani, “Modeling the worldwide spread of pandemic Influenza: baseline case and containment interventions,” PLoS Medicine, vol. 4, no. 1, pp. 95–110, 2007.
- P. Cooley, L. Ganapathi, G. Ghneim, S. Holmberg, W. Wheaton, and C. R. Hollingsworth, “Using Influenza-like illness data to reconstruct an Influenza outbreak,” Mathematical and Computer Modelling, vol. 48, no. 5-6, pp. 929–939, 2008.
- D. L. Chao, M. E. Halloran, V. J. Obenchain, and I. M. Longini, “FluTE, a publicly available stochastic Influenza epidemic simulation model,” PLoS Computational Biology, vol. 6, no. 1, Article ID e1000656, 2010.
- National Institute of General Medical Sciences, “Models of Infectious Disease Agent Study,” 2008, http://www.nigms.nih.gov/Initiatives/MIDAS/.
- N. M. Ferguson, D. A. T. Cummings, S. Cauchemez et al., “Strategies for containing an emerging Influenza pandemic in Southeast Asia,” Nature, vol. 437, no. 7056, pp. 209–214, 2005.
- I. M. Longini, A. Nizam, S. Xu et al., “Containing pandemic Influenza at the source,” Science, vol. 309, no. 5737, pp. 1083–1087, 2005.
- N. M. Ferguson, D. A. T. Cummings, C. Fraser, J. C. Cajka, P. C. Cooley, and D. S. Burke, “Strategies for mitigating an Influenza pandemic,” Nature, vol. 442, no. 7101, pp. 448–452, 2006.
- T. C. Germann, K. Kadau, I. M. Longini, and C. A. Macken, “Mitigation strategies for pandemic Influenza in the United States,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 15, pp. 5935–5940, 2006.
- S. Eubank, H. Guclu, V. S. A. Kumar et al., “Modelling disease outbreaks in realistic urban social networks,” Nature, vol. 429, no. 6988, pp. 180–184, 2004.
- M. E. Halloran, N. M. Ferguson, S. Eubank et al., “Modeling targeted layered containment of an Influenza pandemic in the United States,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 12, pp. 4639–4644, 2008.
- MIDAS, “Report from the Models of Infectious Disease Agent Study (MIDAS) Steering Committee,” 2004, http://www.nigms.nih.gov/News/Reports/midas_steering_050404.htm.
- Committee on Modeling Community Containment for Pandemic Influenza, “Modeling community containment for pandemic Influenza: a letter report,” 2006, http://www.nap.edu/catalog/11800.html.
- Centers for Disease Control and Prevention (CDC), “Interim pre-pandemic planning guidance: community strategy for pandemic Influenza mitigation in the United States,” 2007, http://www.pandemicflu.gov/professional/community/community_mitigation.pdf.
- M. I. Meltzer, N. J. Cox, and K. Fukuda, “The economic impact of pandemic Influenza in the United States: priorities for intervention,” Emerging Infectious Diseases, vol. 5, no. 5, pp. 659–671, 1999.
- T. Das and A. Savachkin, “A large scale simulation model for assessment of societal risk and development of dynamic mitigation strategies,” IIE Transactions, vol. 40, no. 9, pp. 893–905, 2008.
- Colardo Department of Human Services Division of Mental Health, “Pandemic Influenza: quarantine, isolation and social distancing,” 2009, http://www.flu.gov/news/colorado_toolbox.pdf.
- J. T. Wu, S. Riley, and G. M. Leung, “Spatial considerations for the allocation of pre-pandemic Influenza vaccination in the United States,” Proceedings of the Royal Society B: Biological Sciences, vol. 274, no. 1627, pp. 2811–2817, 2007.
- J. Lawless and J. Lawless, Statistical Models and Methods for Lifetime Data, John Wiley & Sons, New York, NY, USA, 1982.
- O. Diekmann and J. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation, John Wiley & Sons, New York, NY, USA, 2000.
- J. F. Brundage and G. D. Shanks, “Deaths from bacterial pneumonia during 1918-19 Influenza pandemic,” Emerging Infectious Diseases, vol. 14, no. 8, pp. 1193–1199, 2008.
- Centers for Disease Control and Prevention (CDC), “CDC Influenza operational plan,” 2006, http://www.cdc.gov/flu/pandemic/cdcplan.htm.
- World Health Organization (WHO), “WHO guidelines on the use of vaccine and antivirals during Influenza pandemics,” 2004, http://www.who.int/csr/resources/publications/Influenza/11_29_01_A.pdf.
- S. Pasteur, “Influenza A(H1N1) Monovalent Vaccine,” 2009, http://www.fda.gov/downloads/biologicsbloodvaccines/vaccines/approvedproducts/ucm182404.pfd.
- R. J. Blendon, L. M. Koonin, J. M. Benson et al., “Public response to community mitigation measures for pandemic Influenza,” Emerging Infectious Diseases, vol. 14, no. 5, pp. 778–786, 2008.
- M. Z. Sadique, W. J. Edmunds, R. D. Smith et al., “Precautionary behavior in response to perceived threat of pandemic Influenza,” Emerging Infectious Diseases, vol. 13, no. 9, pp. 1307–1313, 2007.
- R. Maunder, J. Hunter, L. Vincent et al., “The immediate psychological and occupational impact of the 2003 SARS outbreak in a teaching hospital,” Canadian Medical Association Journal, vol. 168, no. 10, pp. 1245–1251, 2003.
- E. Robertson, K. Hershenfield, S. L. Grace, and D. E. Stewart, “The psychosocial effects of being quarantined following exposure to SARS: a qualitative study of Toronto health care workers,” Canadian Journal of Psychiatry, vol. 49, no. 6, pp. 403–407, 2004.
- M. Pearson, C. Bridges, and S. Harper, “Influenza vaccination of health-care personnel,” MMWR, vol. 55, no. 2, pp. 1–16, 2006.
- M. T. Keane, M. V. Walter, B. I. Patel et al., “Confidence in vaccination: a parent model,” Vaccine, vol. 23, no. 19, pp. 2486–2493, 2005.
- V. P. Niederhauser, G. Baruffi, and R. Heck, “Parental decision-making for the varicella vaccine,” Journal of Pediatric Health Care, vol. 15, no. 5, pp. 236–243, 2001.
- S. D. Rhodes and K. C. Hergenrather, “Exploring hepatitis B vaccination acceptance among young men who have sex with men: facilitators and barriers,” Preventive Medicine, vol. 35, no. 2, pp. 128–134, 2002.
- S. L. Rosenthal, R. K. Kottenhahn, F. M. Biro, and P. A. Succop, “Hepatitis B vaccine acceptance among adolescents and their parents,” Journal of Adolescent Health, vol. 17, no. 4, pp. 248–254, 1995.
- M. S. Smailbegovic, G. J. Laing, and H. Bedford, “Why do parents decide against immunization? The effect of health beliefs and health professionals,” Child: Care, Health and Development, vol. 29, no. 4, pp. 303–311, 2003.
- T. Safranek, D. Lawrence, L. Kuriand, et al., “Reassessment of the association between Guillain-Barré syndrome and receipt of swine Influenza vaccine in 1976-1977: results of a two-state study,” American Journal of Epidemiology, vol. 133, no. 9, p. 940, 1991.
- K. M. Cummings, A. M. Jette, B. M. Brock, and D. P. Haefner, “Psychosocial determinants of immunization behavior in a swine Influenza campaign,” Medical Care, vol. 17, no. 6, pp. 639–649, 1979.
- The New Yorker, “The fear factor,” 2009, http://www.newyorker.com/talk/comment/2009/10/12/091012taco_talk_specter.
- The New York Times, “Doctors swamped by swine flu vaccine fears,” 2009, http://www.msnbc.msn.com/id/33179695/ns/health-swine_flu/.
- A. Gosavi, T. K. Das, and S. Sarkar, “A simulation based learning automata framework for solving semi-Markov decision problems under long run average reward,” IIE Transactions in Operations Engineering, vol. 36, no. 6, pp. 557–567, 2001.
- A. Savachkin, A. Uribe-Sánchez, T. Das, D. Prieto, A. Santana, and D. Martinez, “Supplemental data and model parameter values for cross-regional simulation-based optimization testbed,” 2010, http://imse.eng.usf.edu/pandemic/supplement.pdf.
- U.S Census Bureau, “2001 American community survey,” 2000, http://www.census.gov/prod/2001pubs/statab/sec01.pdf.
- Bureau of transportation statistics, “2001 National household travel survey (NTHS),” 2002, http://www.bts.gov/programs/national_household_travel_survey/.
- “Tampa International Airport: daily traffic volume data,” 2010, http://www.tampaairport.com.
- “Miami International Airport: daily traffic volume data,” 2010, http://www.miami-airport.com.
- “Miami International Airport: daily traffic volume data,” 2010, http://www.jaa.aero/General/Default.
- “Tallahassee Regional Airport: daily traffic volume data,” 2010, http://www.talgov.com/airport/index.cfm.
- “How thermal-imaging cameras spot flu fevers,” 2010, http://www.msnbc.msn.com/id/30523865/.
- Writing committee of the World Health Organization (WHO), “Consultation on human Influenza A/H5.Avian Influenza A(H5N1) infection in humans,” The New England Journal of Medicine, vol. 353, no. 13, pp. 1374–1385, 2005.
- Institute of Medicine (IOM), Antivirals for Pandemic Influenza: Guidance on Developing a Distribution and Dispensing Program, The National Academies Press, Washington, DC, USA, 2008.
- J. J. Treanor, J. D. Campbell, K. M. Zangwill, T. Rowe, and M. Wolff, “Safety and immunogenicity of an inactivated subvirion Influenza A (H5N1) vaccine,” The New England Journal of Medicine, vol. 354, no. 13, pp. 1343–1351, 2006.
- R. J. Blendon, C. M. DesRoches, M. S. Cetron, J. M. Benson, T. Meinhardt, and W. Pollard, “Attitudes toward the use of quarantine in a public health emergency in four countries,” Health Affairs, vol. 25, no. 2, pp. 15–25, 2006.
- D. Montgomery, Design and Analysis of Experiments, John Wiley & Sons, New York, NY, USA, 2007.
- T. Halfhill, “Inflation calculator,” 2009, http://www.halfhill.com/inflation.html.
- D. Bell, A. Nicoll, K. Fukuda et al., “Nonpharmaceutical interventions for pandemic Influenza, national and community measures,” Emerging Infectious Diseases, vol. 12, no. 1, pp. 88–94, 2006.