Measles

public health

How International Travel and Vaccine Resistance have led to a resurgence of Measles in the U.S.

This resurgence of measles cases is due to two main risk factors. The first, as previously noted by Olive et. al. and reiterated by Hotez, is a fall in vaccination rates fueled by anti-vaccination movements in both the U.S. and Europe. In the U.S. in particular there has been an upsurge of non-medical exemptions (NMEs) in the 18 states which allow it. Compounding the decrease in vaccination rates is the increase in international travel into these localities (where measles is supposed to have been eliminated through vaccination) from regions where measles remains endemic, and outbreaks are ongoing. The fact that measles is a highly contagious virus capable of airborne spread further increases the seriousness of this threat.

Olive et al. previously presented a “heat map” of measles risk based on NMEs for vaccination in U.S. counties. In collaboration with Sahotra Sarkar at UT Austin and others (Aleksa Zlojutro at UNSW Sydney and Kamran Khan at University of Toronto), we further investigated how this “heat map” of measles risk in the US might shift geographically when county level NME rates (that reveal pockets of low vaccination) were compounded with incoming international air travel volume weighted by the size of measles outbreaks abroad.

To do so, we performed an analysis that compounds four factors: 1) international air travel volume arriving from measles affected countries into each U.S. county, 2) county NME rates, 3) county population, and 4) the size (incidence rate) of the measles outbreak at travel origin country.  The analysis was conducted for each year between 2011 and 2019. The detailed methods and results of the analysis are soon to be published in Lancet Infectious Diseases. An interactive map illustrating the results of our analysis and corresponding data is available here. Additionally, the figures below illustrate the U.S. measles risk for 2019 computed by the model, and the corresponding global level measles outbreak data from WHO. (NOTE: the 2019 analysis uses WHO reported cases from each country up to April 19.)

Results from the 2019 analysis reveal the U.S. counties at highest risk of a measles outbreak at present, and are spatially consistent with the U.S. measles cases reported to date (April 19): either the counties we identify or those immediately adjacent to them are the ones that have reported measles cases.

Critically, the results correctly predict the areas in Washington, Oregon, and New York that have seen major measles outbreaks. The risk analysis also reveals U.S. localities that have not experienced a measles outbreak in 2019 but are at risk of imported measles cases resulting in a local outbreak because they lie adjacent to a county that has and/or is served by a major international airport. These include Travis (TX), Maricopa (AZ), Clayton (GA), Honolulu County (HI), Wayne (MI), Salt Lake (UT), Hennepin (MN), Suffolk (MA), Loudon (VA), San Diego (CA) and multiple counties in Florida. Additional surveillance should also target Cook (IL) and Los Angeles (CA), which have seen only one case so far, but due to the presence of major international airports, they may serve as the fulcrum of continuous importation of measles into the United States.

Lastly, this analysis also reveals the set of countries that contribute most to measles risk across the U.S. in 2019. The top countries ordered by risk posed are: Ukraine, Mexico, Cuba, Israel, Japan, Thailand and Philippines. Therefore, we recommend that surveillance should also be directed towards those U.S. counties with the high incoming passenger volume from these countries, which should themselves be targeted for vaccination efforts. Prof. Sarkar suggests additional actions for prevention in his blog.

For further information please see our publication in Lancet Infectious Diseases (Published Online May 9, 2019).

Figure 1: Top 25 U.S. Counties Predicted to be at Highest Risk for Measles in 2019.

This results in this figure correspond to our publication in Lancet Infectious Diseases (Published Online May 9, 2019).  Risk is measured by the expected relative size of a measles outbreak in a county.

Figure 2: Annual Measles Cases by Country as reported by WHO
public health

Media Response

WLRN Podcast

List of reports

sElected links

Inverse optimization is an area of study where the purpose is to infer the unknown parameters of an optimization problem when a set of observations is available on the previous decisions made in the settings of the problem. We develop a framework to effectively and efficiently infer the cost vector of a linear optimization problem based on multiple observations on the decisions made previously. 

We then test our models in the setting of a diet problem on a data-set obtained from NHANES; The data-set is accessible via the link bellow:

https://github.com/CSSEHealthcare/Dietary-Behavior-Dataset

A set of female individuals with the above criteria were considered. Further demographic and diet considerations (in order to select similar patients) led to selecting 11 different individuals’ one day of intake as the initial dataset for the model. In another setting, we only considered people that have consumed a reasonable amount of sodium and water. We consider these two nutrients as the main constraints in the DASH diet. 



In order to compare different potential data and their performance with the model, we used different data groups from the NHANES database. A group of middle-aged women with certain similar characteristics and a group of people with certain attributes in their diets. In the first group, we did not consider how the individual’s daily diet is reflecting on the constraints that the forward problem had and we relied on their own personal answer to questions regarding hypertension and also how prone they thought they were to type-2 diabetes. The result was a sparse set of variables and an inconclusive optimal solution in regards to the preferences. In the second group, we tried to obtain sub-optimal data. We prioritized the maximum sodium intake constraint and the water intake constraints as our main and most important constraints. 

We introduce a new approach that combines inverse optimization with conventional data analytics to recover the utility function of a human operator. In this approach, a set of final decisions of the operator is observed. For instance, the final treatment plans that a clinician chose for a patient or the dietary choices that a patient made to control their disease while also considering her own personal preferences. Based on these observations, we develop a new framework that uses inverse optimization to infer how the operator prioritized different trade-offs to arrive at her decision. 

We develop a new inverse optimization framework to infer the constraint parameters of a linear (forward) optimization based on multiple observations of the system. The goal is to find a feasible region for the forward problem such that all given observations become feasible and the preferred observations become optimal. We explore the theoretical properties of the model and develop computationally efficient equivalent models. We consider an array of functions to capture various desirable properties of the inferred feasible region. We apply our method to radiation therapy treatment planning—a complex optimization problem in itself—to understand the clinical guidelines that in practice are used by oncologists. These guidelines (constraints) will standardize the practice, increase planning efficiency and automation, and make high-quality personalized treatment plans for cancer patients possible.

Assume that a decision-maker’s uncertain behavior is observed. We develop a an inverse optimization framework to impute an objective function that is robust against misspecifications of the behavior. In our model, instead of considering multiple data points, we consider an uncertainty set that encapsulates all possible realizations of the input data. We adopt this idea from robust optimization, which has been widely used for solving optimization problems with uncertain parameters. By bringing robust and inverse optimization together, we propose a robust inverse linear optimization model for uncertain input observations. We aim to find a cost vector for the underlying forward problem such that the associated error is minimized for the worst-case realization of the uncertainty in the observed solutions. That is, such a cost vector is robust in the sense that it protects against the worst misspecification of a decision-maker’s behavior. 

As an example, we consider a diet recommendation problem. Suppose we want to learn the diet patterns and preferences of a specific person and make personalized recommendations in the future. The person’s choice, even if restricted by nutritional and budgetary constraints, may be inconsistent and vary over time. Assuming the person’s behavior can be represented by an uncertainty set, it is important to find a cost vector that renders the worst-case behavior within the uncertainty set as close to optimal as possible. Note that the cost vector can have a general meaning and may be interpreted differently depending on the application (e.g., monetary cost, utility function, or preferences). Under such a cost vector, any non-worst-case diet will thus have a smaller deviation from optimality.  

Radiation therapy is frequently used in diagnosing patients with cancer. Currently, the planning of such treatments is typically done manually which is time-consuming and prone to human error. The new advancements in computational powers and treating units now allow for designing treatment plans automatically.

To design a high-quality treatment, we select the beams sizes, positions, and shapes using optimization models and approximation algorithms. The optimization models are designed to deliver an appropriate amount of dose to the tumor volume while simultaneously avoiding sensitive healthy tissues. In this project, we work on finding the best beam positions for the radiation focal points for Gamma Knife® Perfexion™, using quadratic programming and algorithms such as grassfire and sphere-packing.

In radiation therapy with continuous dose delivery for Gamma Knife® Perfexion™, the dose is delivered while the radiation machine is in movement, as oppose to the conventional step-and-shoot approach which requires the unit to stop before any radiation is delivered. Continuous delivery can increase dose homogeneity and decrease treatment time. To design inverse plans, we first find a path inside the tumor volume, along which the radiation is delivered, and then find the beam durations and shapes using a mixed-integer programming optimization (MIP) model. The MIP model considers various machine-constraints as well as clinical guidelines and constraints.

Perioperative services are one of the vital components of hospitals and any disruption in their operations can leave a downstream effect in the rest of the hospital. A large body of evidence links inefficiencies in perioperative throughput with adverse clinical outcomes. A regular delay in the operating room (OR), may lead to overcrowding in post-surgical units, and consequently, more overnight patients in the hospital. Conversely, an underutilization of OR is not only a waste of an expensive and high-demand resource, but it also means that other services who have a demand are not able to utilize OR. This mismatch in demand and utilization may, in turn, lead to hold-ups in the OR and cause further downstream utilization. We investigate the utilization of operating rooms by each service. The null hypothesis of this work is that the predicted utilization of the OR, i.e., the current block schedule, matches completely with the actual utilization of the service. We test this hypothesis for different utilization definitions, including physical and operational utilization and reject the null hypothesis. We further analyze why a mismatch may exist and how to optimize the schedule to improve patient flow in the hospital.

Primary care is an important piece in the healthcare system that affects the downstream medical care of patients heavily. There are specific challenges in primary care as healthcare shifts from fee-for-service to population health management and medical home, focuses on cost savings and integrates quality measures. We consider the primary care unit at a large academic center that is facing similar challenges. In this work we focus on the imbalance in workload, which is a growing regulatory burden and directly concerns any staff in primary care. It can result in missed opportunities to deliver better patient care or providing a good work-environment for the physicians and the staff. We consider the primary care unit at the large academic center and focus on their challenge in balancing staff time with quality of care through a redesign of their system. We employ optimization models to reschedule providers’ sessions to improve the patient flow, and through that, a more balanced work-level for the support staff. 

This work was performed with the MIT/MGH Collaboration.

In many healthcare services, care is provided continuously, however, the care providers, e.g., doctors and nurses, work in shifts that are discrete. Hence, hand-offs between care providers is inevitable. Hand-offs are generally thought to effect patient care, although it is often hard to quantify the effects due to reverse causal effects between patients’ duration of stay and the number of hand-off events. We use a natural randomized control experiment, induced by physicians’ schedules, in teaching general medicine teams. We employ statistical tools to show that between the two randomly assigned groups of patients, a subset who experiences hand-off experience a different length of stay compared to the other group.

This work was performed with the MIT/MGH Collaboration.

Many outpatient facilities with expensive resources, such as infusion and imaging centers, experience surge in their patient arrival at times and are under-utilization at other times. This pattern results in patient safety concerns, patient and staff dissatisfaction, and limitation in growth, among others. Scheduling practices is found to be one of the main contributors to this problem.

We developed a real-time scheduling framework to address the problem, specifically for infusion clinics. The algorithm assumes no knowledge of future appointments and does not change past appointments. Operational constraints are taken into account, and the algorithm can offer multiple choices to patients.

We generalize this framework to a new scheduling model and analyze its performance through competitive ratio. The resource utilization of the real-time algorithm is compared with an optimal algorithm, which knows the entire future. It can be proved that the competitive ratio of the scheduling algorithm is between 3/2 and 5/3 of an optimal algorithm.

This work was performed with the MIT/MGH Collaboration.

Tracking COVID-19

We are tracking the COVID-19 spread in real-time on our interactive dashboard with data available for download. We are also modeling the spread of the virus. Preliminary study results are discussed on our blog.