Factors in the Control of Measles Epidemics: A Simulation
Factors in the Control of Measles Epidemics: A Simulation
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) 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, 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.
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).
(Enlarge Image)
Figure 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.
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.
(Enlarge Image)
Figure 2.
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.
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.
(Enlarge Image)
Figure 3.
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.
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.
(Enlarge Image)
Figure 4.
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.
Results
Sensitivity Analysis
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) 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.
Scenario Analysis
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, 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.
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).
(Enlarge Image)
Figure 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.
(Enlarge Image)
Figure 2.
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.
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.
(Enlarge Image)
Figure 3.
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.
(Enlarge Image)
Figure 4.
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.