text:
The role of vaccination coverage, individual behaviors, and the public health response in the control of measles epidemics: an agent-based simulation for California
Fengchen Liu1,
Wayne T A Enanoria2,
Jennifer Zipprich3,
Seth Blumberg1,
Kathleen Harriman3,
Sarah F Ackley1,
William D Wheaton4,
Justine L Allpress4 &
Travis C Porco1,2,5
Measles cases continue to occur among susceptible individuals despite the elimination of endemic measles transmission in the United States. Clustering of disease susceptibility can threaten herd immunity and impact the likelihood of disease outbreaks in a highly vaccinated population. Previous studies have examined the role of contact tracing to control infectious diseases among clustered populations, but have not explicitly modeled the public health response using an agent-based model.
We developed an agent-based simulation model of measles transmission using the Framework for Reconstructing Epidemiological Dynamics (FRED) and the Synthetic Population Database maintained by RTI International. The simulation of measles transmission was based on interactions among individuals in different places: households, schools, daycares, workplaces, and neighborhoods. The model simulated different levels of immunity clustering, vaccination coverage, and contact investigations with delays caused by individuals' behaviors and/or the delay in a health department's response. We examined the effects of these characteristics on the probability of uncontrolled measles outbreaks and the outbreak size in 365 days after the introduction of one index case into a synthetic population.
We found that large measles outbreaks can be prevented with contact investigations and moderate contact rates by having (1) a very high vaccination coverage (≥ 95%) with a moderate to low level of immunity clustering (≤ 0.5) for individuals aged less than or equal to 18 years, or (2) a moderate vaccination coverage (85% or 90%) with no immunity clustering for individuals (≤18 years of age), a short intervention delay, and a high probability that a contact can be traced. Without contact investigations, measles outbreaks may be prevented by the highest vaccination coverage with no immunity clustering for individuals (≤18 years of age) with moderate contact rates; but for the highest contact rates, even the highest coverage with no immunity clustering for individuals (≤18 years of age) cannot completely prevent measles outbreaks.
The simulation results demonstrated the importance of vaccination coverage, clustering of immunity, and contact investigations in preventing uncontrolled measles outbreaks.
Measles is a highly infectious, acute viral illness that can cause severe pneumonia, diarrhea, encephalitis, as well as death in rare cases [1]. Measles cases continue to occur among susceptible individuals despite the elimination of endemic measles transmission in the United States as a result of importation [2,3]. The introduction of one imported measles case into a group or community with a large number of susceptible individuals can yield large outbreaks [4]. In California, the number of measles cases associated with visits to one of two Disney theme parks was 68 on January 23, 2015 (www.cnn.com/2015/01/23/health/california-measles-outbreak/index.html). In 2013, an intentionally unvaccinated adolescent returned to New York from London, United Kingdom while infectious with measles. As a result, the imported case led to the identification of 58 cases and six generations of measles infection among members of an orthodox Jewish community living in two neighborhoods in Brooklyn [5], resulting in the largest outbreak of measles in the United States since 1996 [5,6]. The outbreak was propagated by a few extended families that either refused the measles, mumps, and rubella (MMR) vaccine or delayed receipt of the vaccine. High population-level vaccination coverage outside of this community coupled with the insular nature of the community itself likely prevented further measles transmission [5]. Despite high vaccination coverage of a population, clustering of disease susceptibility due to personal beliefs and other exemptions to mandatory childhood vaccination can threaten herd immunity and impact the likelihood of disease outbreaks [7].
Previous researchers have used mathematical models to explore these empirical observations. The measles vaccination threshold to attain herd immunity in a population is believed to be approximately 90%--95%, but these estimates were based on models with relatively simple assumptions about contacts among individuals in the population, e.g., random mixing [8]. Mathematical models have examined infectious disease dynamics over specific contact networks [9] and household structure [10,11], particularly preferential mixing among small closely connected groups to show its effects on the reproductive numbers of infectious diseases [11]. In general, population-level disease dynamics are influenced by heterogeneous contact patterns [12]. Furthermore, clustering may decrease the size of epidemics, but also may decrease the epidemic threshold, making it easier for a disease to spread [13]. Furthermore, outbreaks can occur even with very high vaccination rates as the clustering of vaccination increases [8].
With each identified measles case, public health departments investigate the close contacts of the case in order to identify anyone who may have been exposed while the case was infectious. If susceptible close contacts are identified, contact tracing efforts can lead to public health interventions (i.e., giving exposed contacts post-exposure vaccine or post-exposure immune globulin) with the hope of minimizing the chance of disease development. In addition, exposed individuals are asked to monitor their symptoms for at least one incubation period, and report any symptoms to public health authorities, should they develop symptoms consistent with measles. However, people do not always comply with public health interventions or recommendations. Also, public health authorities cannot always find identified contacts, or do not always reach susceptible contacts soon enough in order for post-exposure prophylaxis to be effective; contacts would be placed in quarantine in these instances. With vaccination coverage in the community, clustering of susceptibility, individual behaviors, and contacts occurring over networks of individuals, one question that remains is under which circumstances do contact tracing activities and the associated interventions given by a local or state health department help to control measles epidemics? The purpose of our study was to identify the characteristics of population vaccination coverage, clustering of susceptibility, individuals' behaviors, and the public health response that impact the occurrence of measles epidemics.
The project received a determination of not being research involving human subjects when reviewed by the Committee for the Protection of Human Subjects at the California Department of Public Health.
Model structure
We developed an agent-based simulation model of measles transmission after the introduction of one index case into a population in order to study the effects of contact tracing and public health interventions on the occurrence of measles epidemics. The agent-based simulation model has two main components: (1) the agent-based simulation model of measles transmission using the Framework for Reconstructing Epidemiological Dynamics (FRED) [14,15], and (2) a synthetic population of individuals adapted from the Synthetic Population Database maintained by RTI International [16]. In this section, we describe the various features of our model: the natural history of measles, the synthetic population, the contact network and transmission assumptions, clusters of immunity, contact investigations, public health interventions, and the analysis of the results. Briefly, all individuals in the synthetic population were in a social network and were given characteristics consistent with the U.S. Census and the California Department of Finance population estimates (www.dof.ca.gov), and measles could be transmitted through interactions among individuals in this network. The model simulated different levels of vaccination coverage, different levels of immunity clustering, and delays in contact investigations. We examined the effect of changes in the various features of the model on changes in the number of secondary cases, the total outbreak size (i.e., the total number of diagnosed cases), and the probability of uncontrolled outbreak scenarios over one year (i.e., chains of transmission continuing longer than one year).
Natural history of measles
The natural history of measles was modeled as follows: an infected individual was assumed to be infectious for approximately nine days (four days before rash onset through four days after the onset of rash [1,17]) after a latent period of approximately eight days [18]. A ten-day symptomatic period [19-24] began 8 to 14 days after exposure [25,26] with rash onset beginning 7 to 19 after exposure [27]. An individual who recovered from an infection was no longer susceptible to subsequent infection. In the agent-based model, we used a standard SEIR model (susceptible-exposed/latent-infectious-recovered) [28] to represent an individual's health states.
Synthetic population
Our model incorporated a synthetic population generated from a Synthetic Population Database developed by RTI International [16] that represented every individual in a specific geographic region. The synthetic population was generated based on the U.S. Census and the California Department of Finance population estimates, and each individual had associated demographic information (e.g., age, sex) and locations for social activity (e.g., household, neighborhood, and possibly school or workplace). To reflect measles transmission among preschool-aged children [29], we added daycares to the synthetic population (using daycare data from the California Department of Social Services' Community Care Licensing Division: www.ccld.ca.gov) according to location information for daycare-aged children. Further daycare details in the synthetic population are given in the Additional file 1.
Contact network and transmission
The simulation of measles transmission was based on interactions among individuals in different places: households, schools, daycares, workplaces, and neighborhoods. Depending on the assigned locations of social activity based on the synthetic population data, each individual visited his or her household and neighborhood on each simulated day. The individual could have visited his/her school, daycare, or workplace if the individual was a student, a daycare-aged child or employed, respectively. In addition, the individual could interact with other individuals who shared the same activity locations for the same simulated day.
For an infectious individual, the number of contacts on a given day depends on the contact rate per day (the time unit of simulation is one day) in each of the locations the infectious individual has visited, as well as the probability of transmission per contact. The contact rate per day depends on the type of place of the contact, as shown in Table 1. We assume that, on a given day, a child attending a daycare may have 3 to 20 random contacts with other children who are in the same daycare; similarly, a student attending a school may randomly contact 3 to 20 schoolmates (base on the child-staff ratio for daycares in California: www.daycare.com/california, and the average class size in California public schools: www.cde.ca.gov/ds/sd/dr/cefteachavgclssize.asp); a worker may contact 0 to 7 co-workers (according to the RTI's synthetic population used by the model, 61.4% workplaces have the number of workers ≤7); each individual may have 0.5 to 7 neighborhood contacts; each individual contacts each of his/her household members with a probability ranging from 0.01 to 1. We assume that the transmission probability per contact is very high, as measles is a highly contagious disease [30].
Table 1 Parameters
Whenever an infective individual contacts another individual, the chance that disease transmission will actually occur depends on the immunity level of the contact. The immune status of each simulated individual is a function, in part, of the age-specific vaccination coverage and the efficacy of vaccine (see Section: Clusters of immunity). Once the susceptible individual is infected, the individual is infectious for approximately nine days [1,17] starting the day after the latent period, which is approximately eight days long [18]. The probability distributions of the numbers of days spent in the latent period and infectious period are given in the Additional file 1. For each transmission episode, data about the episode (i.e., date, infector, infectee, and location) are recorded by the simulation program, and this transmission record is used in the simulations involving contact tracing and public health interventions (see Section: Contact investigation and public health intervention).
Measles transmission is initiated by introducing one randomly selected index case, whose age is consistent with the observed age distribution of measles cases in California, into a highly vaccinated population; the initialization of vaccination will be discussed in Section Clusters of immunity. After the latent period, the index case may transmit measles by contacting susceptible individuals in different locations (household, neighborhood, school, daycare, or workplace settings) on each day of the infectious period. Each new infectious individual caused by the index case may then transmit measles via the individual's contact network during his or her infectious period. To model an infected individual's daily behavior during the symptomatic period, we assume that the individual may decide to stay home with a specified probability (ranging from 0 to 1, Table 1) on each day of the symptomatic period until recovery.
To reflect the difference of social activities between weekdays and weekends during the simulation, schools and daycares are assumed to be closed on weekends. Approximately 34% of workers are randomly chosen to continue visiting their workplaces on weekends (based on data from the Bureau of Labor Statistics: www.bls.gov/news.release/atus.nr0.htm) and the number of neighborhood contacts per day on weekends is assumed to increase by 50% [15].
Clusters of immunity
Despite high vaccination coverage, clustering of disease susceptibility due to personal beliefs and other exemptions to mandatory childhood vaccination can threaten herd immunity and impact the likelihood of disease outbreaks [7]. Salathe et al. [31] showed that the effect of clustered susceptible individuals on disease outbreak probabilities is strong when the vaccination coverage is close to the level required to provide herd immunity under the assumption of random mixing.
To model the clusters of susceptible children in a highly vaccinated population (i.e., vaccination coverage of two-dose of MMR was assumed to range between 95% to 99% (V) for children less than or equal to 18 years of age [32-34] and 92.4% to 97.9% for individuals greater than 18 years of age [34]), we proceed as follows. On one extreme, we simply assume random vaccination of children; under this assumption, the vaccination status of each child is independent of the status of all other children in the household; we denote the coverage level by V. At the other extreme, we assume that while the overall coverage is still V, all children in a given household have the same vaccine status. In the model, we interpolate between these two extremes, with Ω = 0 representing random vaccination, and Ω = 1 denoting completely positive assortative vaccine status. The interpolation formula is given in the Additional file 1; values of Ω between 0 and 1 represent intermediate configurations between random vaccination and completely positive assortative vaccination of children. For each synthetic household, we assign vaccination status to each household member in this way to initialize the clusters of immunity in the synthetic population with a given level of clustering, Ω, (assuming that the efficacy of two doses of MMR vaccine is 99% [35,36], Table 1). The simulation of measles transmission discussed in Section Contact network and transmission is in a highly vaccinated and susceptibility-cluster-specific population. Note that detailed age-specific vaccine coverage data by household are not available.
Contact investigation and public health intervention
For an identified case, the public health department cannot know all the contacts of the case, and may not be able to find every contact of the case. To simulate a contact investigation, which may depend on the type of place and the probability that the contacts of an infectious case can be traced, we proceed as follows. First, we assume that neighborhood contacts are not traceable. However, each of the contacts from the household, school, daycare, and workplace settings can be identified by the infectious individual with some specified probability. In addition, the identified contact can be found by the public health department with a probability ranging from 0.7 to 1 (in Table 1).
For an individual who is traced from an infectious case by a contact investigation and receives a public health intervention, several days may have elapsed since exposure. For an infectious individual (either an index case or a case not prevented by the contact tracing and public health interventions) who visits the doctor and is confirmed as a measles case, there may have been several days since the first day of symptoms. To simulate a contact investigation, we assume there are three delays: (i) the self report delay for an infectious individual, i.e., the time between symptoms and diagnosis, (ii) the intervention delay for contacts of an index case, and (iii) the delay for finding a contact of an infectious case.
The self report delay for an infectious individual, i.e., the time between symptoms and diagnosis, is assumed to be between one and six days (J. Zipprich, personal communication, Table 1). The intervention delay for contacts of an index case occurs if the infectious individual is an index case, there may be an additional one to three days (intervention delay in Table 1), which is the time between diagnosis and initiation of contact investigation. If the infectious individual is not the index case but instead is a subsequent case, we assume a further round of contact investigation is inaugurated without delay. The delay for finding a contact of an infectious case occurs when the health department may need an additional one to three days (contact tracing delay in Table 1) to obtain the names of the infectious individual's contacts, trace them, test immunity and provide interventions.
For an individual who is traced by a contact investigation and receives a public health intervention, to simulate the delay from the first day of symptoms to the day of receiving a public health intervention, we assume that there are four scenarios:
If the individual is an index case, the individual receives the public health intervention at a time ranging from 2 to 9 days after the first day of symptoms, i.e., the sum of the self report delay and the intervention delay;
If the individual is the infector and an infectious case but not an index case, the individual receives the public health intervention at a time ranging from 1 to 6 days (the self report delay) after the first day of symptoms;
If the individual is the infectee and the infector is an index case, the individual receives the public health intervention at a time ranging from 3 to 12 days (the sum of three delays) after exposure;
If the individual is the infectee and the infector is an infectious but not an index case, the individual receives the public health intervention at a time ranging from 2 to 9 days (the sum of the self report delay and the contact tracing delay) after exposure.
According to California Department of Public Health [37], post-exposure prophylactic (PEP) MMR vaccine may be given within 72 hours of exposure to persons ≥ 6 months of age. Post-exposure prophylactic immune globulin (IG) may be given within 6 days to exposed persons of any age, though persons weighing less than 30 kg are prioritized. Home quarantine may be recommended for exposed persons contacted seven days or more from exposure. To simulate these three public health interventions (PEP MMR, PEP IG and home quarantine), we assume that one dose of PEP MMR is given to an individual if it has not been more than three days since the first day of exposure, the individual is willing to accept PEP MMR, and the individual's age is greater than one year of age (using one year as the cutoff is consistent with the public health department's practice). We assume that PEP IG is given to an individual if the individual's age is ≤ 1 year and it has not been more than six days since the first day of exposure, or if the individual's age is >1 year, the individual is cooperative, and it has been more than three days but not more than six days since exposure. For an individual who is not willing to accept PEP MMR or PEP IG, or if it has been more than six days since the first day of exposure, we assume that home quarantine is recommended.
Briefly, we simulated contact tracing activities and public health interventions in four steps:
Record all contacts. On each simulated day, the simulation of measles transmission (Section Contact network and transmission) records the information for every contact, including the day when a transmissible contact occurs, the infector (the identity of the infectious individual), the infectee (the identity of the individual contacted by the infector), and the place (location/setting where transmissible contact occurred).
Trace contacts. When a case is diagnosed, the case is interviewed and contacts elicited, which is modeled by choosing a random sample of the true contacts (excluding neighborhood contacts).
Model public health interventions for elicited contacts (including simulated time delays and loss to follow-up). For each elicited contact, a random Bernoulli trial is simulated with the contact finding probability, and if the value is 1 ("success"), then a random waiting time is added. The individual will then be contacted for any interventions at that time in the future (a random variable with an exponential distribution). Only contacts within the last seven days are elicited (by assumption).For an individual who has been given any intervention in a previous contact investigation, we do not schedule another public health intervention within seven days since the last intervention.
Give public health interventions. First, not all individuals cooperate with investigators; we, therefore, generate a random Bernoulli trial with the cooperation probability. If the result is 1 ("success"), then we determine how much time has elapsed since transmission for each individual with a scheduled public health intervention. Either PEP MMR, PEP IG, or a home quarantine recommendation is provided, depending on the elapsed time.
We analyzed the model in two steps. First, we selected a Latin Hypercube sample [38,39] of parameters chosen uniformly from the parameter ranges given in Table 1, and simulated one year of transmission with contact investigation and public health intervention. For each parameter set, we introduced a single index case into the synthetic population, and determined the total number of new cases diagnosed during the year, and whether or not active transmission was still occurring at the end of the year. In many cases, transmission chains die out rapidly, while for other simulated scenarios, transmission continues to occur. For each parameter set, we replicated the epidemic 256 times (28, chosen arbitrarily, but sufficient to ensure a 50% proportion has an estimated standard deviation less than 5%). Thus, for each scenario (parameter set), we determined (1) the number of secondary cases, (2) the outbreak size, and (3) the probability that there were over 500 cases within the year (uncontrolled transmission or "escape" probability). For each replication, we also calculated the number of symptomatic cases directly caused by the index case (allowing the expected number of secondary cases per introduction to be estimated). For a given parameter set, the simulation also monitored the outbreak size for each of 256 iterations; if the outbreak size was more than 500, the iteration was classed as uncontrolled transmission. The probability of uncontrolled transmission for each scenario was estimated as the number of introductions leading to uncontrolled transmission, divided by the total number of introductions (256 in each case). Second, we plotted the results after one year of simulated measles transmission with contact investigations and public health interventions for selected parameter sets.
To determine which parameters were the most important in determining the number of secondary cases, the outbreak size, and the probability of uncontrolled outbreak in a one year of measles transmission, we selected a Latin Hypercube sample of 1024 different parameter sets from the parameter ranges indicated in Table 1. From each, the number of secondary cases, the number of outbreak size, and the probability of uncontrolled measles outbreaks (over one year) were computed from 256 replications. For definiteness, we chose Humboldt County in California (population size 118,261). The simulation was started by introducing one index case who was randomly selected according to the observed age distribution of measles cases in California. We then computed the partial rank correlation coefficient (PRCC) [40,41] between each input parameter and each of three simulation outcomes (the number of secondary cases, the total outbreak size and the probability of uncontrolled measles outbreaks in 365 days). When the PRCC is close to zero, the value of the parameter has little relation to the simulation output; when the PRCC is close to +1 or −1, the value of the parameter is highly important in determining the simulation output.
We chose 4 parameters (PRCC ≥ 0.1, in Table 2), which are (1) vaccination coverage, (2) level of immunity clustering (Ω), (3) intervention delay, and (4) contact finding probability to explore factors which contribute to the successful control of measles (i.e., the probability of uncontrolled measles outbreak). In this section, we examine (1) the combinations of vaccination coverage and the level of immunity clustering, (2) the combinations of intervention delay and contact finding probability based on different levels of immunity clustering and vaccination coverage, (3) the combinations of vaccination coverage and contact rates in the neighborhood, and (4) the combinations of vaccination coverage and contact rates in all places. We report the probability of escape, or the mean outbreak size (when no scenarios showed a total number of cases exceeding 500 by one year after measles introduction). Due to the critical threshold level for the fraction of susceptible individuals below which introduction of infections can only lead to minor outbreaks and above which minor and major outbreaks can occur [42], we present the probability of uncontrolled outbreaks for scenarios in which minor and major outbreaks occurred and total outbreak size for scenarios that only led to minor outbreaks [43].
Table 2 PRCC of parameters with secondary cases, outbreak size and probability of uncontrolled measles outbreak in 365 days
Effects of vaccination coverage and the level of immunity clustering
For combinations of vaccination coverage ranging from 0.8 to 1 and the level of immunity clustering ranging from 0 to 1 (all other parameter values are shown in Table 1), the results are shown in Figure 1. When at least one escape scenario occurred (out of 256 trials), the result is shown in red; when there were no escape scenarios, the mean outbreak size is shown in blue; the border between an escape scenario and a no escape scenario is shown in black. Figure 1 shows that the vaccination coverage required to prevent measles outbreaks was approximately 95% for all levels of immunity clustering. For the scenario with the lowest level of immunity clustering, the threshold vaccination coverage was 80%. For the combinations with which there were no measles outbreaks, Figure 1 also shows that the outbreak size dramatically decreased with an increase in vaccination coverage for a given level of immunity clustering. These simulations simply illustrate the importance of high vaccination coverage and that clustering of unvaccinated individuals is unfavorable for disease control. The results assumed the contact rates, self-report delays, intervention delays for contacts of the index case, and the probability that a contact can be traced were constant (Table 1).
Effects of vaccination coverage and clustering of immunity on the control of measles epidemics. For each combination of vaccination coverage V and the level of immunity clustering Ω (all other parameters' values are shown in Table 1), we ran 256 iterations to obtain the outbreak size and the uncontrolled outbreak probability in 365 simulated days. The combinations with which the simulated uncontrolled outbreak probabilities >0 are represented by red cells and scaled from light red (lower uncontrolled outbreak probability) to dark red (higher uncontrolled outbreak probability); the values of simulated uncontrolled outbreak probabilities are shown in red cells. The combinations without uncontrolled outbreaks (the simulated uncontrolled outbreak probabilities =0) are shown by blue cells and scaled from light blue (higher outbreak size) to dark blue (lower outbreak size); the values of simulated outbreak sizes are shown in blue cells. The frontiers between adjacent combinations with and without uncontrolled outbreaks are shown by the black lines. These simulations suggest that the vaccination coverage is important in the control of measles epidemics (the higher vaccination coverage, the lower the probability of uncontrolled measles outbreaks and the smaller outbreak size); for a given vaccination coverage, a lower level of immunity clustering (i.e., the lower the chance of unvaccinated individuals clustered together in a household) may have better control of measles epidemics.
Effects of intervention delay and contact finding probability based on different levels of immunity clustering and vaccination coverage
To expand further on the relationship between vaccination coverage and immunity clustering and their role in the control of measles outbreaks, we examined the influence of two parameters related to contact investigations, intervention delay for contacts of an index case and the contact finding probability, on the occurrence of measles outbreaks. In Figure 2, we present the results of our findings for three levels of vaccine coverage (V = 0.85, 0.90, and 0.95 in the three columns) and three levels of immunity clustering (Ω=0, 0.5, and 1.0 in the three rows). For each coverage-clustering combination, we show the results of our simulations across values of the contact finding probability (x-axis ranging from 0.50 to 1.0) and the intervention delay for contacts of the index case (y-axis ranging from 1 to 7 days). In the setting with the lowest vaccination coverage and the highest level of immunity clustering (V=0.85 and Ω=1.0, Figure 2(A)), we found that escape scenarios could occur for all values of the intervention delay and the contact finding probability, and thus report the probability of escape scenarios denoted in red. In the setting with the highest vaccination coverage and the lowest level of immunity clustering (V=0.95 and Ω=0.0, Figure 2(I)), we obtained no uncontrolled outbreaks, and report the total outbreak size for all values of the intervention delay and the contact finding probability (shown in blue). The settings where the levels of the intervention delay and the contact finding probability become important were in Figure 2(C), (E), and (G), when the coverage-clustering combinations were either both high (Figure 2(C)), both moderate (Figure 2(E)), or both low (Figure 2(G)). For these three settings, the level of intervention delay for contacts of the index case and the probability a contact can be found determines whether we have controlled or uncontrolled measles outbreaks. Even with stochastic simulation error, higher contact finding probabilities with lower intervention delays lead to controlled measles outbreaks.
Effects of intervention delays and the contact finding probability on the control of measles epidemics. We simulated the main outcomes for each combination of the intervention delay for contacts of the index case and the contact finding probability under 9 combinations (A, B, C, D, E, F, G, H and I) of vaccination coverage (V= 0.85, 0.9 and 0.95) and the level of immunity clustering (Ω=0, 0.5 and 1). Cells in red indicate the combinations with which the simulations had uncontrolled outbreaks with the uncontrolled outbreak probabilities ranging from low to high. Cells in light blue show the combinations without uncontrolled outbreaks but with higher outbreak sizes than the combinations without uncontrolled outbreaks represented by dark blue cells. The frontiers between adjacent combinations with and without uncontrolled outbreaks are shown by the black lines. The results in (A), (B) and (C) suggest that increasing vaccination coverage levels and the contact finding probability while reducing intervention delays may reduce the uncontrolled outbreak probability. However, measles outbreaks may not be prevented with the highest level of immunity clustering. The results in (D), (E) and (F) suggest that scenarios with a lower level of immunity clustering, increasing vaccination coverage and contact finding probability, and reducing intervention delays have smaller uncontrolled outbreak probabilities than the highest level of immunity clustering (shown in Figures 2(A), (B) and (C)), and 95% of vaccination may be enough to prevent measles outbreak which is consistent with the result shown in Figure 1 (the cell with V=95% and Ω=0.5); when V=90%, measles outbreaks may be prevented by the combinations of low intervention delays and a high contact finding probability. The results in (G), (H) and (I) suggest that with the lowest level of immunity clustering, increasing vaccination coverage may dramatically reduce measles outbreaks.
However, Figure 2 does not show the effects of the speed of contact investigation on the uncontrolled probability and the outbreak size. To study the effects of the speed of contact investigations, we varied the self report delay, the intervention delay and the contact tracing delay and found that the uncontrolled outbreak probabilities and the outbreak sizes were reduced by increasing the speed of contact investigation (as shown in Table 3). In general, the speed of contact investigations does affect how well it contributes to disease control.
Table 3 Effects of the speed of contact investigation on outbreak sizes and uncontrolled outbreak probabilities
Effects of vaccination coverage and contact rate in neighborhood
The PRCC value for the neighborhood contact rate (in Table 2) indicates that the transmission in the neighborhood setting is highly important in determining the probability of uncontrolled measles outbreaks and the total outbreak size because of the assumption that contact tracing does not happen in the neighborhood (i.e., people are not likely to identify the people they encounter in the neighborhood setting). We used the parameters' values in Table 1 and combinations of vaccination coverage (ranging from 0 to 1) and the neighborhood contact rate (ranging from 1 to 7) to study their effects on the uncontrolled measles outbreak probability and total outbreak size. The simulation for each combination was run under three scenarios: (1) without contact investigation, (2) with contact investigation and little intervention delay for contacts of the index case, and (3) with contact investigation and more intervention delay for contacts of the index case. For the scenario without contact investigation, Figure 3(A) shows that more neighborhood contacts caused more outbreaks; increasing vaccination coverage could reduce the number of measles cases, but even with the highest vaccination coverage (100%), it still could not completely prevent the measles outbreak for high contact rates in neighborhood. For the scenario with contact investigation and less intervention delay for contacts of the index case, as shown Figure 3(B), with lower contact rates in the neighborhood, measles could be controlled even with 80% vaccination coverage; and 100% of vaccination coverage could prevent measles outbreaks for all contact rates in the neighborhood. For the scenario with contact investigations but with more intervention delays for contacts of the index case, Figure 3(C) shows that measles outbreaks were prevented by the combinations of low contact rates in the neighborhood with all vaccination coverage levels, and the combinations of the highest vaccination coverage level with all contact rates in neighborhood. Most of the uncontrolled outbreak probabilities and outbreak sizes were higher than those of the scenario with contact investigation and less intervention delays for contacts of the index case, i.e., there is less of a chance for a traced contact to receive PEP MMR or PEP IG with more intervention delays. By comparing the same combinations of the scenario without contact investigation in Figure 3(A), the probabilities of uncontrolled outbreaks and the outbreak sizes among the other two scenarios in Figures 3(B) and (C) were dramatically reduced by contact investigations. These simulations suggest: (1) contact investigation plays an important role in preventing measles outbreaks and reducing the outbreak size; (2) with contact investigations, reducing the contact rate in the neighborhood may lower the vaccination coverage required to prevent measles outbreaks; (3) with contact investigation and the highest vaccination coverage, measles outbreaks may be prevented even with very high contact rates in neighborhood; and (4) less intervention delays for contacts of the index case may help contact investigations reduce the probability of uncontrolled measles outbreaks and the outbreak size.
Effects of vaccination coverage and the contact rate in neighborhoods on the control of measles epidemics. We used the parameters' values in Table 1 and combinations of vaccination coverage and contact rates in the neighborhood, and ran 256 iterations for each combination to simulate the uncontrolled outbreak probability and the outbreak size in 365 days under three scenarios: (A) without contact investigations, (B) with contact investigations and less intervention delay for contacts of the index case, and (C) with contact investigations and more intervention delay for contacts of the index case. For each of the scenarios A, B, and C, red cells show the combinations with uncontrolled outbreaks (the simulated uncontrolled outbreak probabilities >0; light red cells indicate combinations with lower uncontrolled outbreak probabilities; dark red cells indicate combinations with higher uncontrolled outbreak probabilities); and blue cells represent the simulated outbreak sizes of the combinations without uncontrolled outbreaks (the simulated uncontrolled outbreak probabilities =0), and scale the outbreak sizes from low (dark blue) to high (light blue). The frontiers between adjacent combinations with and without uncontrolled outbreaks are shown by the black lines. These simulations suggest: contact investigation plays an important role in preventing measles uncontrolled outbreaks and reducing the total outbreak size; with contact investigations, reducing the contact rates in the neighborhood may lower the vaccination coverage required to prevent uncontrolled measles outbreaks; with contact investigations and the highest vaccination coverage, measles outbreaks may be prevented even with very high contact rates in the neighborhood; less intervention delays for contacts of the index case may help contact investigations reduce the probability of uncontrolled measles outbreaks and the total outbreak size.
Effects of vaccination coverage and contact rates in all places
To study further the effects of contact rates in all places and vaccination coverage on the uncontrolled measles outbreak probability and outbreak size, we used a parameter ranging from 0 and 1 to scale contact rates in school, daycare, workplace, household and neighborhood settings (i.e., when the scale parameter is 0, each place has the lowest contact rate; when the scale parameter is 1, each place has the highest contact rate; when the scale parameter is between 0 and 1, the contact rate in each place is between its lower and upper bounds shown in Table 1). Similarly, for each combination of vaccination coverage (ranging from 0 to 1) and the scale parameter (ranging from 0 to 1), we ran the simulations under three scenarios (all other the parameters' values are shown in Table 1): (1) without contact investigation, (2) with contact investigation and less intervention delay for contacts of the index case, and (3) with contact investigation and more intervention delay for contacts of the index case. For the scenario without contact investigation, Figure 4(A) shows that the more contact rates in all places the more outbreaks and the higher number of outbreak sizes; increasing vaccination coverage could reduce the number of measles cases; when scale of contact rates in all places was very high, even the highest vaccination coverage (100%) could not completely prevent the measles outbreak; however, for the lowest scale of contact rates, even 80% of vaccination could prevent measles outbreak. For the scenario with contact investigation and less intervention delay for contacts of the index case, Figure 4(B) shows that with lower scale of contact rates in all places, measles could be controlled even with 80% of vaccination coverage; and 100% of vaccination coverage could prevent measles outbreak for all scales of contact rates in all places. For the scenario with contact investigations but with more intervention delays among contacts of the index case, Figure 4(C) shows that measles outbreaks were prevented by the combinations of low contact rates in all places and all vaccination coverage levels, and the combinations of the highest vaccination coverage and all contact rates in all places; most of the uncontrolled outbreak probabilities and the outbreak sizes were higher than those of the scenario with contact investigations and less intervention delays for contacts of the index case (in Figure 4(B)) because the effects of intervention delays. By comparing this to the same combinations of the scenario without contact investigations in Figure 4(A), the probabilities of an uncontrolled outbreak and the outbreak sizes in the other two scenarios in Figures 4(B) and (C) were dramatically reduced by the contact investigations. These simulations suggest: (1) contact investigations play an important role in preventing measles outbreaks and reducing the total outbreak sizes; (2) without contact investigations but with the lowest scale of contact rates in all places, even an 80% vaccination coverage may prevent measles outbreaks; (3) with contact investigations, reducing the contact rates in all places may lower the vaccination coverage required to prevent measles outbreaks; (4) with contact investigations and the highest vaccination coverage level, measles outbreaks may be prevented even with very high contact rates in all places; (5) less intervention delays for contacts of the index case may help contact investigations reduce the number of measles cases.
Effects of vaccination coverage and contact rates in all places on the control of measles epidemics. We used a scale parameter ranging from 0 to 1 to scale contact rates in school, daycare, workplace, household and neighborhood. For each combination of vaccination coverage V and the scale of all contact rates, the main outcomes were obtained by running the simulation for each of three scenarios: (A) without contact investigations, (B) with contact investigations and less intervention delays for contacts of the index case, and (C) with contact investigations and more intervention delays for contacts of the index case. For scenarios A, B, and C, red cells show the combinations with uncontrolled outbreaks and blue cells represent the simulated outbreak sizes of the combinations without uncontrolled outbreaks. The values of simulated uncontrolled outbreak probabilities are shown with blue numbers in red cells; the values of simulated outbreak sizes are shown with red numbers in blue cells. The frontiers between adjacent combinations with and without uncontrolled outbreaks are shown by the black lines. These simulations suggest: (1) contact investigations play an important role in preventing uncontrolled measles outbreaks and reducing the total outbreak size; (2) without contact investigations but with the lowest scale of contact rates in all places, even an 80% vaccination coverage may prevent uncontrolled measles outbreaks; (3) with contact investigations, reducing the contact rates in all settings may lower the vaccination coverage required to prevent uncontrolled measles outbreaks; (4) with contact investigations and the highest vaccination coverage level, uncontrolled measles outbreaks may be prevented even with very high contact rates in all settings; (5) less intervention delay for contacts of the index case may help contact investigations reduce the number of measles cases.
We constructed a stochastic agent-based model to explore the impact of various characteristics (population vaccination coverage, clustering of immunity, individual behavior, and the public health response) on the occurrence of measles epidemics. Our analysis focused on contact investigations and their associated interventions (post-exposure prophylaxis) to control and prevent measles transmission. Our results confirm previous observations that vaccination coverage and clustering of immunity are important characteristics that influence the ability to control a measles epidemic. We were able to expand on our understanding of vaccination coverage and clustering by demonstrating that the prevention of uncontrolled measles outbreaks may depend on the level of intervention delays and the probability of finding a contact when vaccination coverage is high with a low to moderate level of clustering. In a highly immunity-clustered population, the prevention of uncontrolled measles outbreaks may only be possible with high vaccination coverage levels (e.g., 95% vaccination coverage). Furthermore, we showed that contact rates are important drivers of transmission, especially neighborhood contact rates given our assumptions about our ability to identify contacted individuals in the neighborhood setting. Reducing contact rates in all settings may lower the vaccination coverage level required to prevent measles outbreaks in the presence of contact investigations and their associated interventions. Overall, the results demonstrated the importance of contact tracing activities to prevent uncontrolled measles outbreaks or reduce the total outbreak size.
Previous mathematical epidemiology studies have examined the factors that are needed to sustain measles transmission in communities [44-46], including critical community size [47-49] and herd immunity [50,51]. In addition, the effects of contact patterns across network structures, including clustered and complex networks, have also been elucidated using mathematical models in order to design effective control strategies [8,9,13,52-55]. While there have been studies that have examined the role of contact tracing to control infectious diseases among clustered populations [56-58], none have modeled the public health response using an agent-based model as we have done in this study. The agent-based model created using an open-source software system (FRED) allows refinements to the model to be made for further explorations of the effects of the public health response, individual behaviors, population vaccination coverage, and network structures on the ability to control and prevent measles outbreaks.
However, the study is subject to several limitations. First, we do not have a lot of empirical evidence for many of the behavior-related parameters in our model. We obtained, where possible, expert opinion and/or published estimates from the literature to inform our parameter values. We incorporated a broad range of values for the behavioral-related parameters in order to explore the wide range of outcomes and conducted a sensitivity analysis of the simulation results. Second, we maintained a balance between core individual behaviors and public health response behaviors in our model with specific nuanced scenarios or possibilities. For example, we understand that children ages 6 months to 11 months can be given one dose of MMR vaccine instead of intramuscular IG if it is given within 72 hours of exposure so long as the child is vaccinated again after 12 months of age; our stochastic agent-based model does not allow for this possibility with this age group. Third, we only considered vaccination clustering status among children less than 18 years of age. Clustering among children occur because parents make decisions for their children regarding vaccination. However, there are little data on vaccination clustering among adults. Finally, the contact network of this model did not include the healthcare settings where measles cases interact with healthcare workers and healthcare workers who become cases and work while infectious may have a large number of contacts.
Our results support the use of contact investigations in the control and prevention of measles epidemics; the results elucidate the circumstances in which control is possible for various levels of vaccination coverage and clustering of susceptibility among households, and the simulation results are consistent with the findings in [59] in which the authors found that minimal changes in the level of aggregation of unvaccinated individuals can lead to sustained transmission in highly vaccinated populations.
We used an agent-based model to simulate measles transmission after the introduction of one index case into a population and study the effects of the level of immunity clustering, vaccination coverage, contact rates, contact tracing and public health interventions on the uncontrolled measles outbreak probability and the total outbreak size in 365 days. The simulated results show that the level of immunity clustering, vaccination coverage and contact investigation are important in the prevention of measles outbreaks. With contact investigations and moderate contact rates in all places, measles outbreaks can be prevented by (1) a very high vaccination coverage (≥95%) with a moderate to low level of immunity clustering (≤0.5), or (2) a moderate vaccination coverage (85% or 90%) with the lowest level of immunity clustering and a short intervention delay for contacts of the index case and a high probability that a contact can be traced. Without contact investigations, measles outbreaks may be prevented by the highest vaccination coverage with the lowest level of immunity clustering and moderate contact rates. For the highest contact rates in all places, even the lowest level of immunity clustering with the highest vaccination coverage cannot completely prevent measles outbreaks. Two directions for future work include using this model (with calibration) to study: (1) the effects of public health intervention strategies on the control and prevention of measles epidemics, and (2) the effects of clustering of susceptibility in school and daycare settings on vaccination coverage required to prevent measles outbreaks.
Strebel PM, Papania MJ, Halsey NA. In: Plotkin SA, Orenstein WA, Offit PA, editors. Vaccines. Philadelphia: Saunders; 2004. p. 389–440.
Katz SL, Hinman AR. Summary and conclusions: measles elimination meeting, 16–17 March 2000. J Infect Dis. 2004;189 Suppl 1:S43–7.
Orenstein WA, Papania MJ, Wharton ME. Measles elimination in the United States. J Infect Dis. 2004;189 Suppl 1:S1–3.
Parker AA, Staggs W, Dayan GH, Ortega-Sánchez IR, Rota PA, Lowe L, et al. Implications of a 2005 measles outbreak in Indiana for sustained elimination of measles in the United States. N Engl J Med. 2006;355(5):447–55.
Centers for Disease Control and Prevention. Notes from the field: measles outbreak among members of a religious community - Brooklyn, New York, March-June 2013. MMWR Morb Mortal Wkly Rep. 2013;62(36):752–3.
Centers for Disease Control and Prevention. Measles outbreak–Southwestern Utah, 1996. MMWR Morb Mortal Wkly Rep. 1997;46(33):766–9.
May T, Silverman RD. 'Clustering of exemptions' as a collective action threat to herd immunity. Vaccine. 2003;21(11–12):1048–51.
McBryde ES. Network structure can play a role in vaccination thresholds and herd immunity: a simulation using a network mathematical model. Clin Infect Dis. 2009;48(5):685–6. author reply 686.
Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2(4):295–307.
Ball F, Sirl D, Trapman P. Analysis of a stochastic SIR epidemic on a random network incorporating household structure. Math Biosci. 2010;224(2):53–73.
Goldstein E, Paur K, Fraser C, Kenah E, Wallinga J, Lipsitch M. Reproductive numbers, epidemic spread and control in a community of households. Math Biosci. 2009;221(1):11–25.
Bansal S, Grenfell BT, Meyers LA. When individual behaviour matters: homogeneous and network models in epidemiology. J R Soc Interface. 2007;4(16):879–91.
Newman MEJ. Properties of highly clustered networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2003;68(2 Pt 2):026121.
Grefenstette J, Brown S, DePasse J, Galloway D, Lee B, Weng Y-T, et al. FRED user's guide. Version 24 January 2013. University of Pittsburgh, Carnegie-Mellon University, Pittsburgh Supercomputing Center and RTI International. 2013. http://fred.publichealth.pitt.edu/downloads/documentation/pdf/FRED_UserGuide.pdf.
Grefenstette JJ, Brown ST, Rosenfeld R, Depasse J, Stone NT, Cooley PC, et al. FRED (A Framework for Reconstructing Epidemic Dynamics): an open-source software system for modeling infectious diseases and control strategies using census-based populations. BMC Public Health. 2013;13:940.
Wheaton W. U.S. synthetic population 2010 version 1.0: quick start guide. RTI International. 2014. http://www.epimodels.org/10_Midas_Docs/SynthPop/2010_synth_pop_ver1_quickstart.pdf.
Kutty P, Rota J, Bellini W, Redd SB. In: Roush SW, McIntyre L, Baldy LM, editors. Manual for the surveillance of vaccine-preventable diseases. 6th ed. Atlanta, GA: Centers for Disease Control and Prevention; 2013. p. 1–21. http://www.cdc.gov/vaccines/pubs/surv-manual/chpt07-measles.pdf.
Keeling MJ, Rohani P. Modeling infectious diseases in humans and animals. New Jersey: Princeton University Press; 2008.
Sartwell PE. The distribution of incubation periods of infectious disease. Am J Epidemiol 1995. 1949;141(5):386–94. discussion 385.
Sartwell PE. The distribution of incubation periods of infectious disease. Am J Hyg. 1950;51(3):310–8.
Simpson REH. Infectiousness of communicable diseases in the household (measles, chickenpox, and mumps). Lancet. 1952;260(6734):549–54.
Simpson REH. The period of transmission in certain epidemic diseases; an observational method for its discovery. Lancet. 1948;2(6533):755–60.
Nishiura H, Yan P, Sleeman CK, Mode CJ. Estimating the transmission potential of supercritical processes based on the final size distribution of minor outbreaks. J Theor Biol. 2012;294:48–55.
Goodall EW. Incubation period of measles. Brit Med J. 1931;1:73–4.
Lessler J, Reich NG, Brookmeyer R, Perl TM, Nelson KE, Cummings DAT. Incubation periods of acute respiratory viral infections: a systematic review. Lancet Infect Dis. 2009;9:291–300.
Anderson RM, May RM. Directly transmitted infections diseases: control by vaccination. Science. 1982;215(4536):1053–60.
Katz SL. In: Gershon AA, Hotez PJ, Katz SL, editors. Infectious diseases of children. 11th ed. Philadelphia: Mosby; 2004. p. 353–71.
Schenzle D. An age-structured model of pre- and post-vaccination measles transmission. IMA J Math App Med Biol. 1984;1:169–91.
Yip FY, Papania MJ, Redd SB. Measles outbreak epidemiology in the United States, 1993–2001. J Infect Dis. 2004;189 Suppl 1:S54–60.
Centers for Disease Control and Prevention. In: Atkinson W, Wolfe S, Hamborsky J, editors. Epidemiology and prevention of vaccine-preventable diseases.12th edition, second printing. Washington DC: Public Health Foundation; 2012. p. 173–92.
Salathé M, Bonhoeffer S. The effect of opinion clustering on disease outbreaks. J R Soc Interface. 2008;5(29):1505–8.
Centers for Disease Control and Prevention. National, state, and local area vaccination coverage among children aged 19–35 months-United States, 2012. MMWR Morb Mortal Wkly Rep. 2013;62(36):733.
California Department of Public Health. Vaccine-preventable diseases and immunization coverage in California, 2001–2006. Center for Infectious Diseases, Division of Communicable Disease Control, Immunization Branch. 2006. http://www.cdph.ca.gov/programs/immunize/Pages/VaccinePreventableDisease.aspx.
McQuillan GM, Kruszon-Moran D, Hyde TB, Forghani B, Bellini W, Dayan GH. Seroprevalence of measles antibody in the US population, 1999–2004. J Infect Dis. 2007;196(10):1459–64.
Uzicanin A, Zimmerman L. Field effectiveness of live attenuated measles-containing vaccines: a review of published literature. J Infect Dis. 2011;204 Suppl 1:S133–48.
Demicheli V, Jefferson T, Rivetti A, Price D. Vaccines for measles, mumps and rubella in children. Cochrane Database Syst Rev. 2005;4:CD004407.
California Department of Public Health. Measles investigation quicksheet. California Department of Public Health. 2015. http://www.cdph.ca.gov/programs/immunize/Documents/CDPHMeaslesInvestigationQuicksheet.pdf.
Inman RL, Helson JC, Campbell JE. An approach to sensitivity analysis of computer models: part I-ranking of input variables, response surface validation, distribution effect and technique synopsis. J Qual Technol. 1981;13(4):174–83.
Inman RL, Helson JC, Campbell JE. An approach to sensitivity analysis of computer models: part II-ranking of input variables, response surface validation, distribution effect and technique synopsis. J Qual Technol. 1981;13(4):232–40.
Blower SM, Dowlatabadi H. Sensitivity and uncertainty analysis of complex models of disease transmission: an HIV model, as an example. Int Stat Rev. 1994;62(2):229–43.
Porco TC, Holbrook KA, Fernyak SE, Portnoy DL, Reiter R, Aragón TJ. Logistics of community smallpox control through contact tracing and ring vaccination: a stochastic network model. BMC Public Health. 2004;4(1):34.
Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proc R Soc London A. 1927;115:700–21.
Wallinga J, Heijne JCM, Kretzschmar M. A measles epidemic threshold in a highly vaccinated population. PLoS Med. 2005;2(11), e316.
Keeling MJ, Grenfell BT. Understanding the persistence of measles: reconciling theory, simulation and observation. Proc Biol Sci. 2002;269(1489):335–43.
Bolker B, Grenfell B. Space, persistence and dynamics of measles epidemics. Philos Trans R Soc Lond B Biol Sci. 1995;348(1325):309–20.
Keeling MJ. Modelling the persistence of measles. Trends Microbiol. 1997;5(12):513–8.
Bartlett MS. The critical community size for measles in the United States. J R Stat Soc. 1960;123(1):37–44.
Black FL. Measles endemicity in insular populations: critical community size and its evolutionary implication. J Theor Biol. 1966;11(2):207–11.
Keeling MJ, Grenfell BT. Disease extinction and community size: modeling the persistence of measles. Science. 1997;275(5296):65–7.
Black FL. The role of herd immunity in control of measles. Yale J Biol Med. 1982;55(3–4):351–60.
Fine PE. Herd immunity: history, theory, practice. Epidemiol Rev. 1993;15(2):265–302.
Kiss IZ, Green DM. Comment on "properties of highly clustered networks". Phys Rev E Stat Nonlin Soft Matter Phys. 2008;78(4 Pt 2):048101.
Miller JC. Spread of infectious disease through clustered populations. J R Soc Interface. 2009;6(41):1121–34.
House T. Exact epidemic dynamics for generally clustered, complex networks. 2010. http://arxiv.org/abs/1006.3483.
Salathé M, Jones JH. Dynamics and control of diseases in networks with community structure. PLoS Comput Biol. 2010;6(4), e1000736.
Eames KTD, Keeling MJ. Contact tracing and disease control. Proc Biol Sci. 2003;270(1533):2565–71.
Eames KTD. Contact tracing strategies in heterogeneous populations. Epidemiol Infect. 2007;135(3):443–54.
House T, Keeling MJ. The impact of contact tracing in clustered populations. PLoS Comput Biol. 2010;6(3), e1000721.
Dallaire F, Serres GD, Tremblay F-W, Markowski F, Tipples G. Long-lasting measles outbreak affecting several unrelated networks of unvaccinated persons. J Infect Dis. 2009;200(10):1602–5.
Ordman CW, Jennings CG, Janeway CA. Chemical, clinical, and immunological studies on the products of human plasma fractionation. xii. The use of concentrated normal human serum gamma globulin (human immune serum globulin) in the prevention and attenuation of measles. J Clin Invest. 1944;23(4):541–9.
Stokes J, Maris EP, Gellis SS. Chemical, clinical, and immunological studies on the products of human plasma fractionation. xi. The use of concentrated normal human serum gamma globulin (human immune serum globulin) in the prophylaxis and treatment of measles. J Clin Invest. 1944;23(4):531–40.
Young MK, Nimmo GR, Cripps AW, Jones MA. Post-exposure passive immunisation for preventing measles. Cochrane Database of Systematic Reviews 2014(4):CD010056.
This study was supported by a Models of Infectious Disease Agent Study (MIDAS) grant from the National Institute of General Medical Sciences at the National Institutes of Health to the University of California, San Francisco (Award Number U01GM087728). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of General Medical Sciences, the University of California at San Francisco, or RTI.
Francis I. Proctor Foundation, University of California, San Francisco, CA, USA
Fengchen Liu
, Seth Blumberg
, Sarah F Ackley
& Travis C Porco
Department of Epidemiology and Biostatistics, University of California, San Francisco, San Francisco, CA, USA
Wayne T A Enanoria
California Department of Public Health, Immunization Branch, Richmond, CA, USA
Jennifer Zipprich
& Kathleen Harriman
RTI Research Triangle Institute International, San Francisco, CA, USA
William D Wheaton
& Justine L Allpress
Department of Ophthalmology, University of California, San Francisco, San Francisco, CA, USA
Travis C Porco
Search for Fengchen Liu in:
Search for Wayne T A Enanoria in:
Search for Jennifer Zipprich in:
Search for Seth Blumberg in:
Search for Kathleen Harriman in:
Search for Sarah F Ackley in:
Search for William D Wheaton in:
Search for Justine L Allpress in:
Search for Travis C Porco in:
Correspondence to Travis C Porco.
TCP, WTAE, FL, JZ, KH conceived and designed the experiments. WDW and JLA provided the synthetic population database. FL, WTAE, TCP performed the experiments. TCP, WTAE, FL, SB, SFA, JZ and KH analyzed the data. FL, WTAE, TCP, JZ and SB drafted the manuscript. All authors read and approved the final manuscript.
Fengchen Liu and Wayne T A Enanoria contributed equally to this work.
Appendix for "The role of vaccination coverage, individual behaviors, and the public health response in the control of measles epidemics: an agent-based simulation".
Liu, F., Enanoria, W.T.A., Zipprich, J. et al. The role of vaccination coverage, individual behaviors, and the public health response in the control of measles epidemics: an agent-based simulation for California. BMC Public Health 15, 447 (2015) doi:10.1186/s12889-015-1766-6
Measle
Public Health Intervention
Contact Rate
Contact Investigation
Infectious Disease epidemiology


We're going retro today with a blast from the past! Today we are visiting a few of our vintage lines! When we threw this idea out to the design team we didn't put any limits on it...just pick ANY collection older than a year and create something fun! You know our collections never go out of style. They are truly timeless like your favorite sweater or a great classic book. Now come see which Prima treasures they dug up!
I so much love the effect of the Plaster Paste that I used on my background with an amazing floral patterns stencil.
"The layout is done about my younger daughter who is a cunning little fox. She's bright and jumps to conclusions fast – in good and in bad. She learns fast, but also gets bored easily and that's usually when the trouble starts… lol!"
Anita chose one of our most popular collections to frame her gorgeous kitty cat! Rose Quartz is one of those quintessential lines with beautiful imagery, a touch of distressing, and overall a very versatile color palette...which makes it PURRRRFECT for a pet project, right?


text:
How to lose a guy (or gal) in 10 years: Common problems in couples' communication and how to solve them
I'm not a couples' psychologist by training. I don't specialize in relationship communication. But I've worked in a fertility clinic, where even the most loving and glowing relationships are put through the grinder. This is not surprising, given the extremes of emotional, physical, social, and financial pressures couples must deal with as a team in an often unforgiving landscape of loss and uncertainty. In this land of fertility treatment, the most common frictions between couples are highlighted and blown up like textbook examples of what can go wrong in couples' communication.
This week, I will share some of the most universal themes I've observed, along with advice that applies to every couple.
Photo by freestocks.org from Pexels
1. The continental drift.
Sam feels upset about something. Raj can see that she's upset, and he doesn't want to touch the wrong chord, so he carefully avoids acknowledging her bad mood. She sees that he clearly notices, but for some reason he doesn't even care enough to ask. He sees that she's fuming harder, so he retreats to another room, not wanting to set off a bomb. She feels even more pissed off, and there's no place to vent. Eventually, someone says something "wrong" and a fight erupts.
Understandably, people are afraid of "setting of" unpleasantness. Sometimes, we think that if we just give bad emotions a wide berth, they will eventually dissipate and everything will be fine again. But avoiding unpleasantness has a tendency to amplify it, because your avoidance makes it look like you don't care. Remember, you don't necessarily have to solve the problem. Most of the time, partners who want to vent don't even want a solution. Just be a good listener.
Don't think that you've been a good listener just because you remained in the room and allowed sounds to enter your ears. That is actually one of the most frustrating things you can do. You'll be a good listener—and help your partner to feel cared for—when you actively engage in the conversation by reflecting back what you heard ("it sounds like you've had a frustrating day"), asking follow-up questions ("what do you make of that?"), and expressing empathy ("wow, that really sucks!").
2. The faucet and the plumber.
John is sad that his mother's illness is making a turn for the worse. He's looking for a shoulder to cry on and tells Jack about his worries. Jack immediately puts on his problem-solving hat. "Has she gotten a second opinion yet?" "Maybe she can try alternative medicine." "You shouldn't have to shoulder all this by yourself. Ask your sister to help more." Meanwhile, John feels no better (in part because he has already considered all of these solutions), and eventually gets frustrated and walks away. Jack feels indignant and hurt because he's just trying to help.
Are you the "fixer" in the relationship? When there's a leaky faucet, it's very helpful to call a plumber. But many stressors in life and in relationships are not fixable—heartache, loss, frustration, deterioration. When you call a plumber to fix a heartache, how far could he get? In the process of wielding his awkward wrenches and screwdrivers, he may even make the heartache worse, by making his partner feel unheard and invalidated.
Instead of putting on your problem-solving, fixing, plumber hat… just slow down and listen. In fact, listen a lot, and repeat back what you're hearing so your partner knows you've been paying attention and understanding. And if the problem isn't fixable, simply offer a hug and say, "I'm so sorry. Tell me how I can help."
3. The unrelenting optimist.
Clara is feeling hopeless about the future because she and her partner Jamie have tried several rounds of fertility treatment with no success. She can't help crying each time she gets a negative pregnancy test. Jamie thinks he needs to be the emotional rock of the relationship, to be optimistic enough for both of them. Each time Clara voices a fear or doubt, he immediately counters: "I know for a fact that we'll get pregnant." "Third time is the charm, remember?" "I'm not worried at all, because we're meant to be parents." Clara wonders if she's crazy for feeling so hopeless, or if Jamie is not taking this seriously. Eventually she stops talking to him about her feelings, and feels even more alone.
This one is similar to the "plumber" problem above. Jamie is trying to solve the problem of hopelessness by countering it with unrelenting optimism. I understand—you're trying to pull up the average positivity. You're trying to lead by example. But what's really happening is that you're driving your partner away and leaving them feeling isolated. The two of you are farther and farther from being on the same page. Often, situations are truly bad—perhaps not totally hopeless, but it would be a lie to say "I'm sure everything will turn out well." In these cases, your unrealistic optimism makes your partner wonder if he/she is crazy, or feel guilty about dragging you both down. Or they may simply feel like you don't understand the situation or don't care.
Instead of being blindly optimistic… acknowledge the difficulties you face together, tell your partner that you understand why they feel so down, and simply offer, "We're in this together. I've got you."
4. Apples vs. oranges.
Paul asks Mimi if they can review their finances—he's worried that they haven't saved enough for a rainy day. As they come across a big purchase from last year, Mimi is reminded of how embarrassed she was by Paul's scolding on that occasion. Paul remembers this too, but he doesn't want to go there ("what's past is past"). Mimi says she's upset that he hasn't apologized. He's frustrated that they can't just figure out the current budget. They end up putting down the spreadsheets and accusing each other of having bad intentions. Paul never brings up finances again for fear of triggering another meandering fight. Mimi thinks he's incapable of talking about emotions.
Often, couples will say that one of them is the more "emotional" of the pair, while the other is more "rational." This is a false dichotomy. Almost everyone has emotional needs. Almost everyone is capable of rational thinking. Often, the problem is that during one particular moment, one person is trying to prioritize emotional processing, while the other is trying to prioritize problem-solving. Both can be achieved! But you have to take turns. If you agree ahead of time on the rules of a particular engagement, both parties will rest assured that their needs will be met, and will be able to focus on the task at hand.
Separate emotion-processing conversations from problem-solving ones. Instead of sitting down to a serious conversation that is totally open-ended ("I think we should talk about the future")… propose a specific goal for the conversation. For example, "I feel that we ended our conversation last night with bad feelings. Can we set aside problem-solving today and just talk through our feelings about that?" Be generous with your partner's goals too, for example, "I know we need to talk about our budget—how about we set aside 30 minutes tomorrow after dinner to look over spreadsheets?"
Of course, all of these good communication strategies are based on a foundation of good faith and mutual respect. If you care for one another, then open ears, generous empathy, and a little bit of pragmatism will be all the glue you need.
Photo by Jasmine Wallace Carter from Pexels
Jade Wu October 2, 2018
How much sleep do you need?
Jade Wu May 30, 2019
Insomnia 101 - How the brain regulates sleep and how you can help
Jade Wu September 25, 2018

redpajama_set_name: RedPajamaCommonCrawl


text:
Home Celebrities Is Elliott Page Transgender? How His Loved Ones Reacted to News of...
Is Elliott Page Transgender? How His Loved Ones Reacted to News of His Transition?
Fannie Jones
Thursday, 8 December 2022, 11:33 GMT+0530
Canada-born actor Elliot Page. He has garnered numerous honors, such as an Academy Award nomination, two BAFTA Award nominations, and a Primetime Emmy Award nomination, as well as a Satellite Award nomination. In December 2020, Elliot Page came out publicly as transgender. In March 2021, he was the first openly transgender guy to feature on Time's cover.
Elliot Page, the Oscar-nominated actor from "Juno" and "The Umbrella Academy" on Netflix, has revealed he is transgender.
Formerly known as Ellen Page, Elliot addressed his social media fans with the following message:
Hi buddies, I would like to inform you that I am transgender, that my pronouns are he/they, and that my name is Elliot. I am fortunate to be writing this. Being here. To have reached this point in life.
I am indebted to the extraordinary people who have supported me during this trip. I cannot begin to describe how extraordinary it feels to finally love myself enough to follow myself.
I've been incessantly motivated by so many members of the trans community. Thank you for your bravery, generosity, and tireless efforts to make this world more inclusive and compassionate.
I will provide whatever assistance I can and continue to work towards a society that is more caring and equitable,' he added.
A post shared by @elliotpage
I adore being trans. And I enjoy being homosexual. And the more I cherish myself and truly accept who I am, the more I dream, the more my heart expands, and the more I flourish. Page said, "To all the trans people who face daily harassment, self-loathing, abuse, and the danger of violence: I see you, I love you, and I will do everything I can to improve this world."
Page utilizes both he/him and they/them pronouns and identifies as transgender and non-binary, which means his gender identification is neither male nor female.
Director of Transgender Media at GLAAD, Nick Adams, stated, "Elliot Page has given us amazing characters on-screen and has been an ardent advocate for all LGBTQ individuals.
He will now serve as an example for innumerable trans and non-binary individuals. Everyone who identifies as transgender deserves the opportunity to be themselves and be accepted for who they are. Today, we honor the amazing Elliot Page."
In 2005, Page made a name for himself outside of his own Canada with the release of his vengeance thriller "Hard Candy." Two years later, he starred in the classic independent film "Juno" by Jason Reitman and Diablo Cody, for which he was nominated for an Academy Award.
His noteworthy performances include Kitty Pryde in the "X-Men" franchise, Christopher Nolan's "Inception," "Whip It!" and Sony's "Flatliners" revival. He has produced and acted in films such as "Tallulah" and "Freeheld," and made his directorial debut last year with a film on the environmental racism faced by people of color and First Nations communities in Canada.
Is Josh Groban Gay? What is The Sexuality of American Actor Josh Groban?
Is Dixie D'Amelio Gay? What is The Sexuality of American Singer Dixie D'Amelio?
Page was one of the most visible openly gay performers in Hollywood prior to his transformation. The declaration on Tuesday enhances his legacy and adds him to a tiny but rising group of out trans artists and Hollywood icons.
This includes both performers (Laverne Cox, Trace Lynette, and Brian Michael Smith) and creators (The Wachowski Sisters and Joey Soloway, creator of "Transparent")
Variety revealed the news of Netflix's breakthrough overall content production deal with Janet Mock, a prominent trans storyteller, and director whose credits include "Pose" and "Hollywood."
How His Loved Ones Reacted to News of His Transition?
Elliot Page beamed as he said that his life "changed drastically" after he decided to live his truth.
The actor, best known for his role in Netflix's smash hit "The Umbrella Academy," came out as transgender in December 2020, saying that his pronouns would be "he/they." In a statement on Instagram, Page said his new name would be Elliot.
Even though he got a lot of mixed reactions from strangers on social media, Page said that the people closest to him welcomed the news with open arms. "I feel really thankful for that," the 35-year-old said on "Late Night with Seth Meyers."
"Most people were not surprised. It's something I've been talking about with a lot of people in my life." Page says he has learned to ignore negative comments from online "keyboard warriors" who don't like him.
"What I want to focus on right now, and what has been so amazing, is the level of joy and presence I feel," he said.
"I feel in a way that I really never thought was possible for a very long time. So that's really what I'm focusing on and enjoying the most. But, of course, times can be overwhelming."
Page said it's "unfortunate" that he still has to deal with negative comments about his decision to change his gender.
"It's like, we're all on the same team here, ya know? "Because people are so obsessed with the binary and how we're all supposed to live our lives, we all have these expectations, these limits, and these constraints," he said. Page said that his life has gotten "dramatically" better.
"And I hope that people who don't like me can maybe try to hear that or accept it on some level," the actor said.
Previous articleHow Did Joan Crawford Die: The Cause Of Death Is Heart Attack!
Next articleFormer Theranos Executive Sunny Balwani Has Been Sentenced to Over 13 Years in Prison


text:
Closer To Home Series
To raise awareness of the diverse range of water problems faced across the world, our content series highlights recent water challenges in the home countries and cities of the current Manchester City team, Manchester City Women and New York City FC.
Nowhere on the planet is without water challenges – something we're dramatising with the help of the Manchester City men's and women's team, and New York City FC players. Discover water challenges Closer to Home.
Manchester City | Mens
Sergio Agüero kicks off the series by sharing some of the recent water issues in his home city of Buenos Aires.
David Silva talks about the challenges faced in his native Gran Canaria, and mainland Spain.
Ilkay Gündoğan discusses problems caused by recent droughts in his homeland, showing how water problems can affect everywhere and everyone.
Manchester City | Womens
Jill Scott, from South Shields in the northeast of England, introduces the serious water issues the region has faced.
Caroline Weir introduces some of the serious water problems faced north of the border, including flood, drought and contamination.
Steph Houghton shares how Britain's past has contributed to water contamination issues in the north of England.
From wildfires to snowmelt, Anton Tinnerholm delves into the wide-ranging effects of drought and rising temperatures in his home country of Sweden.
Ben Sweat brings to the fore some of the water challenges faced in his hometown of Florida, revealing the risk of coastal flooding and water pollution faced by inhabitants of the sunshine state.
Ismael Tajouri-Shradi dissects some of the water challenges facing his home country of Switzerland, drawing attention to increasing levels of water toxicity, and the likelihood of glacial melting.
Jesus Medina's home country of Paraguay has experienced water challenges at both ends of the scale in recent years, with river flooding and drought posing serious threats to humans and wildlife.
Sebastian Ibeagha unveils some of the water challenges faced in his home country of Nigeria, discussing the lack of access to clean drinking water, flooding, and crude oil pollution.
Closer to home World Water Facts
Learn about the serious water challenges faced in the home cities of a range of Manchester City players - with individual information and statistics.
Find out how you can help with our Top 12 water solving facts and tips.
Want to know about the Raining Champions Rainwater Brew?


text:
Ла́скер-Ринк () — ледовый каток на севере Центрального парка Нью-Йорка, южнее 110 улицы, приблизительно между 106 и 108 улицами Манхэттена.

Открыт в 1966 году.

Летом здесь размещается бассейн, зимой — каток. Зимний режим катка функционирует с октября по март.

Хоккей
На ледовом катке действует Central Park Ice Hockey, устраиваются массовые катания на коньках, в том числе школьные занятия. Central Park Ice Hockey занимается подготовкой и обучением хоккею любого уровня. Каток также является домашним стадионом мужской хоккейной команды Central Park Hawks и женской хоккейной команды The Lady Hawks.

Стадион принимает хоккейную команду Гарлема (), IHIH. За команду выступают любители игры в хоккей. Команда тренируется с 8 до 11 вечера несколько раз в неделю.

Фигурное катание
На стадионе проводятся занятия по фигурному катанию любого уровня, в том числе для тех, кто только начинает учиться и также для тех, кто готовится к соревнованиям разного уровня. Существуют программы обучения для детей, частные уроки для участников разного возраста и учёба танцам на льду для тех кому за…

Факты
Лауреат звания Лучший ледовый каток — 2011, присужденный газетой New York Daily News.

Билеты: взрослые — $6.50, пожилые (60 лет и старше) — $2.50, дети (12 лет и младше) — $3.50.

Стадион принимает соревнования хоккеистов на кубок Ласкера.

Примечания

Ссылки
 Официальный сайт

Здания и сооружения Манхэттена
Центральный парк

redpajama_set_name: RedPajamaWikipedia


text:
Abd-Al·lah ibn Issam al-Khawlaní () fou valí de Mayurqa (913-946). Succeí en el càrrec al seu pare Issam al-Khawlani i el va abandonar per peregrinar a la Meca.

Bibliografia
 Pablo Pi Ferrer, José Maria Quadrado. Las Islas Baleares.

Valís de Mallorca


text:
Zbytka je přírodní rezervace jihozápadně od obce Bohuslavice v okrese Náchod. Oblast spravuje Královéhradecký kraj. Důvodem ochrany je zachování ekosystému lužního lesa, slatin a luk na náplavech Zlatého potoka s výskytem zvláště chráněných druhů rostlin a živočichů, s důrazem na zachování jedinečných půdních profilů.

Galerie

Odkazy

Reference

Související články
 Seznam chráněných území v okrese Náchod

Externí odkazy

Přírodní rezervace v okrese Náchod
Přírodní rezervace v okrese Rychnov nad Kněžnou
Chráněná území v Orlické tabuli
Chráněná území vyhlášená roku 1994
Evropsky významné lokality v Královéhradeckém kraji
Lužní lesy v Královéhradeckém kraji
Lesy v okrese Náchod
Lesy v okrese Rychnov nad Kněžnou
Lesy v Orlické tabuli
Louky v okrese Náchod
Louky v okrese Rychnov nad Kněžnou
Louky v Orlické tabuli
Povodí Dědiny
Bohuslavice (okres Náchod)