Abstract
Background Societal segregation of unvaccinated people from public spaces has been a novel and controversial COVID-era public health practice in many countries. Models exploring potential consequences of vaccination-status-based segregation have not considered how segregation influences the contact frequencies in the segregated groups. We systematically investigate implementing effects of segregation on population-specific contact frequencies and show this critically determines the predicted epidemiological outcomes, focusing on the attack rates in the vaccinated and unvaccinated populations and the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people.
Methods We describe a susceptible-infectious-recovered (SIR) two-population model for vaccinated and unvaccinated groups of individuals that transmit an infectious disease by person-to-person contact. The degree of segregation of the two groups, ranging from zero to complete segregation, is implemented using the like-to-like mixing approach developed for sexually-transmitted diseases, adapted for presumed SARS-CoV-2 transmission. We allow the contact frequencies for individuals in the two groups to be different and depend, with variable strength, on the degree of segregation.
Results Segregation can either increase or decrease the attack rate among the vaccinated, depending on the type of segregation (isolating or compounding), and the contagiousness of the disease. For diseases with low contagiousness, segregation can cause an attack rate in the vaccinated, which does not occur without segregation.
Interpretation There is no predicted blanket epidemiological advantage to segregation, either for the vaccinated or the unvaccinated. Negative epidemiological consequences can occur for both groups.
Introduction
Models can be used to investigate infectious disease dynamics under different hypotheses about the characteristics of a disease and the effects of health policy. In this endeavour, there are advantages to working with simplest-possible but sufficiently realistic models [1,2], where one should exclude simple models that are not sufficiently realistic for the intended application, either because of their structure or because of incorrect assumptions about the underlying mechanisms. Following this approach, researchers have extended the foundational simple susceptible-infectious-recovered (SIR)-type model to explore diseases with birth and death dynamics, maternal- or vaccine-derived immunity, latency of infection, patterns of contact mixing between different societal groups, and so on [3-7], and to study the effect of isolating vulnerable individuals from the general population during a pandemic, in the absence of vaccination [8].
Recently, SIR models of epidemic dynamics have been implemented with two interacting societal groups (vaccinated and unvaccinated) to examine epidemic outcomes for variable degrees of interaction between the two groups, including whether the unvaccinated put the vaccinated unduly or disproportionately at risk, using epidemiological parameters intended to be representative of SARS-CoV-2 [9-12]. These prior implementations regarding groups differentiated by vaccination status take the contact frequencies of the majority and socially-excluded groups to be equal and held constant, irrespective of the degree of segregation (or exclusion or “like-to-like mixing”), which is not realistic.
Here, we implement population-specific contact frequencies that can be different for the two groups and that can either increase or decrease with increasing segregation. This is necessary because, for example, in many actual regulatory policies the excluded unvaccinated group is barred from public venues or services where people gather and from public transport where people are in close proximity for various durations. In general, the contact frequency of the excluded group decreases with increasing segregation if isolation is in effect, and increases with increasing segregation if the excluded individuals are crowded together. Implementing this essential model feature gives rise to more complex behaviour of the attack rates in the vaccinated and unvaccinated populations (Av and Au, respectively), which can increase or decrease, or rise to a maximum before decreasing, as the two groups are increasingly segregated. This is also true for the share (Bv) of infections among vaccinated people that are due to contacts with infectious unvaccinated people.
This article was previously posted to the medRxiv preprint server on August 23, 2022, and revised versions were posted to medRxiv on November 28, 2022; July 6, 2023; and July 19, 2023.
Methods
Model Design
We adopt the standard SIR framework in a structure with two sub-populations. If a susceptible person (S) comes into contact with an infectious person (I), the susceptible person can become infectious, and infectious people eventually recover (R) and become permanently immune.
We divide the population into two groups: vaccinated and unvaccinated. Vaccination is “all or nothing”, such that a proportion VE of the vaccinated population is immune (is in the R state from the outset of the simulation), where the parameter VE represents vaccine efficacy. The model also includes a natural immunity parameter, NI, equal to the proportion of unvaccinated that are immune from the outset due to previous infection [9].
The model parameter η controls the degree of segregation of vaccinated and unvaccinated people. When η = 0, there is no segregation, and the two groups mix randomly. When η = 1, there is complete segregation, such that vaccinated only come into contact with other vaccinated, and unvaccinated only come into contact with other unvaccinated.
The parameter η follows from Garnett and Anderson [1], who modeled sexually-transmitted disease spread in a population divided into groups with different frequencies of sexual contacts. They take the contact frequency to be a constant characteristic of the individuals within a group. However, contact frequency is not generally and solely an intrinsic individual characteristic [13].
In our model, the population-specific contact frequencies of the vaccinated and unvaccinated individuals (cv and cu, respectively) can increase, decrease, or remain constant as the two groups are segregated. We implement a new approach to achieve this: we keep the first two terms in Taylor expansions of cv and cu versus η (
and 
; see Eqs. A3 in Appendix A1). Thus, mv and mu determine the degree of increase or decrease of the contact frequency in either group, as η is increased.
For example, when mu < 0, as segregation is increased, the contact frequency of unvaccinated people decreases. This corresponds to segregation policy that excludes unvaccinated people from public spaces such as restaurants, cinemas, workplaces, airplanes, trains, etc. [14-17]. Conversely, mu > 0 corresponds to segregation policy that increases contacts between unvaccinated people; for example, by requiring returning unvaccinated travelers to stay in designated facilities [18-20].
In principle, the vaccinated and unvaccinated contact frequencies may be different even when the two groups are completely unsegregated. The unsegregated (η = 0) contact frequencies are set by the parameters 
and 
. Similarly, the probability that contact between a susceptible and infectious person results in transmission is βv (βu) for a susceptible vaccinated (unvaccinated) person and the rates of recovery from infection for the vaccinated and unvaccinated individuals are γv and γu, respectively.
There are thus two “β parameters”, two “c parameters” and two “γ parameters” in our model. Since each β parameter always occurs as part of a product with its respective c parameter, the β parameters can freely be set equal to 1. We set βv = βu = 1 in this paper, without any loss of generality. This implies that, by definition (since β = 1), the contact frequencies “c” in our model are conceptually for contacts that are of sufficiently close proximity and long duration that an infection is guaranteed to occur when a susceptible and an infectious person meet. For a more contagious virus, more of an individual’s contacts are long and close enough that transmission would be guaranteed, corresponding to higher 
and 
.
The model of Fisman et al. [9] is the special case of our model with 
and γv = γu, in which case the equal contact frequencies of both vaccinated and unvaccinated remain constant regardless of the level of segregation. Such an implementation does not represent how segregation has been applied during the COVID era in Canada and many countries [14-17,21], since unvaccinated people were excluded from public spaces while vaccinated people were allowed access, thus changing venues and opportunities for contact as segregation is imposed.
Throughout this paper, “contact frequency” refers to frequency of infectious contacts, since the probability of infection per infectious-susceptible contact is set equal to 1 without loss of generality (see Appendix 1).
Model Parameterization
The parameters of our model are listed in Table 1; calculated quantities in Table 2. Technical details of the model are in Appendix 1.
Analysis
The attack rate among the vaccinated population is defined as the proportion of initially-susceptible vaccinated people who become infected during the epidemic:
where Sv(t0) is the number of susceptible vaccinated people at the beginning of the epidemic and Sv(tf) is the number of susceptible vaccinated people remaining once there are no longer any infectious people in the entire (vaccinated and unvaccinated) population. Au is defined equivalently, for the unvaccinated, replacing the v subscripts with u in Eq. 1. The overall attack rate for the full (vaccinated plus unvaccinated population) is:
We also define Bv as the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people (see Eq. A6 of Appendix 1).
We focus on segregation types that are targeted at the unvaccinated group. We assume, for simplicity, that segregation has no impact on the contact frequency of vaccinated people (mv = 0). We also assume that the contact frequencies in both groups are the same when there is no segregation 
. We use the same values as used by Fisman et al. [9] for the remaining parameters: Pv = 0.8, VE=0.8, NI=0.2, γ = 73 yr-1, and N = 107. These values were intended to be representative for COVID-19 and vaccination; in particular, the recovery rate of 73 yr-1 is equivalent to a recovery time of 5 days [22,23], and assumed to be the same for vaccinated and unvaccinated people.
Appendix 2 contains supplementary figures with results for different parameter combinations, including Pv ≠ 0.8, mv ≠ 0 and 
. In all results in this paper, simulations were initiated with a seed number of 100 infectious individuals distributed proportionately among the two sub-populations.
Results
Fig. 1 shows simulation results for a range of model parameters for different epidemiological conditions and degrees and types of societal segregation. Each row of panels is for a fixed value of 
. Moving from the top row (Figs. 1a.i to a.iv) to the bottom (Figs. 1e.i to e.iv), 
. The left column of panels shows how the attack rate among the vaccinated population, Av, changes with the degree of segregation, η. The second and third columns show Au and At as functions of η, respectively, and the right column shows how Bv, the share of vaccinated infections that were due to contacts with unvaccinated people, varies with η.
Attack rates Av (vaccinated population), AU (unvaccinated population), At (overall population) and share of vaccinated infections that were due to contacts with unvaccinated people, Bv, as functions of the degree of segregation, η, of the vaccinated and unvaccinated. Each row of panels shows Av, Au, At and Bv for a particular choice of 
. Values of fixed model parameters are indicated at the top of the figure. For reference, in a single-population (no vaccination) model, the corresponding R0 values for rows a-e of the figure are 8.2, 6.0, 4.1, 3.4 and 2.7, respectively.
Figs. 1c.i to c.iv show results for a moderate value of 
contacts/yr. For reference, in a single-population (no vaccination) model, c = 300 contacts/yr and γv = γu = 73 yr-1 corresponds to a basic reproduction number R0 = c/γ = 4.1.
In Fig. 1c.i, when mu =-1 and mu = -0.5 (reflecting large and moderate degrees of exclusion and isolation of unvaccinated people) the vaccinated attack rate, Av decreases with increasing segregation. However, when mu > 0 (compounding of unvaccinated people) or mu = 0 (segregation has no influence on contact frequency of unvaccinated people), there is a maximum in Av for moderate values of η. Therefore, with compounding segregation, large values of η are required for Av to be lower than its value for no segregation (η = 0). Fig. 1c.ii shows that the unvaccinated attack rate, Au, increases with segregation for anything other than strong isolating segregation (mu approaching -1). This produces a maximum in the overall attack rate At at moderate degrees of segregation, even for values of mu for which Av decreases monotonically (mu = -0.5). Fig. 1c.iv shows that Bv, the share of vaccinated infections that are due to unvaccinated people, has a shape similar to Av(η, mu). In all panels, 20% of the total population is unvaccinated (Pv = 0.8; Table 1).
Figs. 1c.i to c.iv therefore demonstrate that whether segregation increases or decreases the vaccinated-population attack rate depends on both the degree of segregation and how segregation affects contact frequency.
Figs. 1a.i to a.iv and 1b.i to b.iv show results for larger 
. Compared to Fig. 1c.i, Av in Figs. 1a.i and 1b.i does not increase much with η when mu > 0, and Av no longer has a maximum when mu = 0. It can also be seen that Av increases with increasing 
when there is no segregation (η = 0).
Reducing 
(Figs. 1d.i to d.iv and 1e.i to e.iv), decreases Av(η = 0), and larger η can dramatically increase Av. Even with an isolating segregation policy (mu = −0.5 in Fig. 1d.i), Av is increased for moderate values of η.
When 
are small enough (
contacts/yr in Fig. 1e.i, corresponding to R0 = 2.7 in a single population (no vaccination) model), there is no epidemic among the vaccinated in the absence of segregation (Av(η = 0) = 0). However, a non-zero vaccinated-population attack rate (Av > 0) occurs if η is sufficiently large, and emerges regardless of whether one isolates or compounds the unvaccinated. Therefore, for small enough values of 
, any segregation could increase infections among the vaccinated.
The main qualitative features of the above results for Pv = 0.8 hold for other values of Pv. Appendix 2 provides a detailed exploration of results for Pv = 0.1 through 0.99; and for two values of VE (0.4 and 0.8). When VE is decreased, Av is not strongly influenced by η, regardless of mu; therefore, any beneficial effect of segregation on Av is reduced as VE decreases.
Appendix 2 also explores 
. For example, when 
, the unvaccinated contact frequency is reduced even when there is no segregation; increasing η can then increase Av substantially compared to the case of 
, holding all other parameter values constant (see panels a.i and b.i in Figs. A2.28 and A2.31).
Discussion
Segregation can have substantially different and negative impacts on the outcome of an epidemic, depending on the type and degree of segregation, and depending on cultural and population-density factors, for example, that co-determine 
and 
.
Segregation that compounds the unvaccinated (mu > 0 and mv = 0) generally causes an increase in the vaccinated-population attack rate, Av, for small and intermediate degrees of segregation, η, while for large η, Av decreases below its value in an unsegregated society. Segregation that isolates and excludes the unvaccinated (mu < 0 and mv = 0) decreases Av for “more contagious viruses” (i.e. large 
, large R0); however, for “less contagious viruses” (smaller 
, smaller R0), both isolating and compounding types of segregation can increase Av beyond its value in an unsegregated society. For “viruses that are not very contagious” (small 
, small R0), applying segregation can cause a sizeable epidemic among the vaccinated even though virtually no vaccinated people would be infected in an unsegregated society. Segregation increases the unvaccinated attack rate, Au, for compounding and moderately isolating types of segregation, and Au is only decreased for strongly isolating segregation (mu approaching -1).
Except for large negative values of mu, and small unvaccinated population fractions, applying segregation has the effect of increasing the frequency of unvaccinated-to-unvaccinated contacts (see Appendix 1, Fig. A1.3). This increases the overall probability of a susceptible-infectious interaction, since the unvaccinated population has a higher fraction of susceptibles, and creates a form of core group dynamics [24-26]. At the same time, increasing segregation shields the vaccinated population from the increased prevalence of infection in the unvaccinated population. This trade-off causes the non-monotonic relationship between Av and η. The same dynamic causes the emergence of an epidemic for large η when 
(and thus R0) is small.
We find that Bv, the share of vaccinated infections that are due to contact with unvaccinated people, follows a similar trend to Av as a function of the degree of segregation, when segregation has no impact on the vaccinated contact frequency (mv = 0). For this type of segregation, Av and Bv either increase or decrease simultaneously with increasing η, depending on the value of mu, and Bv is minimized for complete segregation. When mv = 0, there is no type or degree of segregation that reduces the vaccinated attack rate while simultaneously increasing the risk to vaccinated people from unvaccinated people (Fig. 1). Therefore, there are no circumstances in which the unvaccinated cause a disproportionate risk to the vaccinated, contrary to conclusions in Fisman et al. [9].
In contrast, when mv ≠ 0, such that segregation affects the contact frequencies of vaccinated people, increasing segregation can cause Av to increase while Bv decreases and vice-versa (see Appendix 2, Figs. A2.25 and A2.26).
The impact of vaccination-status-based societal segregation on contact frequencies has not previously been considered to our knowledge, including in network-based models in which unvaccinated people cluster together in “cliques” or households [26-28].
Limitations
Our model assumes only two risk populations (vaccinated and unvaccinated), considers only the attack rates on epidemic completion (Av, Au and At), and takes the degree of segregation η to be time-independent, without variation due to public holidays and such. It does not consider other outcomes such as death or hospitalization. Our model assumes an all-or-nothing VE, without waning immunity or influence on infectiousness; and no possibility of reinfection. We do not consider the impact of segregation policies on vaccination rates. SIR models and their variations are based on the paradigm of transmission due to pairwise contact between a recently infected and a susceptible individual. However, this paradigm is unable to account for important features of viral respiratory disease incidence, in particular its rapid emergence and disappearance occurring at essentially the same time at widely dispersed locations [29]. Air-borne transmission via suspended aerosol particles is not directly compatible with pairwise transmission, since it occurs in built environments where many people may transit or be present [30]. A related and unavoidable limitation is the lack of reliable empirical evaluations of needed infectious contact frequencies, which is important because our calculated outcomes are sensitive to the chosen contact frequency values. Lastly, we do not consider the deleterious health impacts of the segregation policies themselves, which can be significant [32-38].
Conclusion
In the two-population mixing-model framework, vaccination-status-based societal segregation can lead to substantially different and counter-intuitive epidemic outcomes depending on the type and degree of segregation, and depending on complex cultural and physical factors that co-determine infectious contact frequencies (i.e., the products βc). Negative epidemiological consequences can occur for either segregated group, irrespective of the deleterious health impacts of the policies themselves.
Given the lack of reliable empirical evaluations of needed infectious contact frequency values, given the demonstrated outcome sensitivities to the infectious contact frequencies, and given the intrinsic limitations of SIR models in this application, we cannot recommend that SIR modelling be used to motivate or justify segregation policies regarding viral respiratory diseases, in the present state of knowledge.
Data Availability
This is a modeling paper, therefore all data is from simulations. Information for reproducing the simulations is contained in the manuscript and appendices.
Appendix
Appendix 1: Elaboration of the model
This is Appendix 1 of the pre-print “Viral respiratory epidemic modeling of societal segregation based on vaccination status” by J. Hickey & D.G. Rancourt, uploaded to https://www.medrxiv.org/ on 2023-10-31 (Version 5).
A1.1: Model differential equations and “mixing” rule
The model is a susceptible-infectious-recovered (SIR) model with two populations: vaccinated (subscript v) and unvaccinated (subscript u) people, consisting of the following six differential equations:
Su, Iu, and Ru represent the number of susceptible, infectious, and recovered unvaccinated people, at time t. Nu represents the total number of unvaccinated people. cu represents the contact frequency (number of contacts per unit time) of unvaccinated people. βu is the probability that a susceptible unvaccinated person becomes infected upon contact with an infectious person (regardless of whether the infectious person is vaccinated or unvaccinated). γu is the rate at which infected unvaccinated people recover from infection. The quantities Sv, Iv, Rv, Nv, cv, βv, and γv are defined equivalently, for vaccinated people.
There are thus two “β parameters”, two “c parameters” and two “γ parameters” in our model. Since each β parameter always occurs as part of a product with its respective c parameter, the β parameters can freely be set equal to 1: this imposes that the “contacts” considered in the model are, by definition, only those contacts that are of sufficiently close proximity and long duration that an infection is guaranteed to occur when a susceptible and an infectious person meet. We set βv = βu = 1 in the main text, without any loss of generality.
fij is the probability that a person of type i (either u or v) has a contact with a person of type j (either u or v), and is defined as follows:,
where 0 ≤ η ≤ 1 and δij is the Kronecker delta, such that δuu = δvv = 1 and δuv = δvu = 0. η is therefore a parameter that controls the degree of separation between the u and v sub-populations. For example, when η = 1, then fuu = fvv = 1 and fuv = fvu = 0, such that u people only ever have contacts with other u people and likewise for v people (complete separation). At the other extreme, when η = 0, then the probability that a u person has a contact with a v person is entirely determined by the relative proportions of u and v people, weighted by their respective contact frequencies (no separation, or “random mixing”).
A1.2: Variation of contact frequency with degree of separation
To allow contact frequency to vary with degree of separation, we define:
where 
and 
represent the contact frequencies for unvaccinated and vaccinated people when there is no separation (η = 0), and mu ≥ −1 and mv ≥ −1 are two parameters that control how non-zero separation impacts the contact frequencies of unvaccinated and vaccinated people.
Fig. A1.1 shows cu as a function of η, for different values of the parameter mu and for 
. As can be seen, when 
, such that the contact frequency is constant regardless of the degree of separation.
Variation of the unvaccinated contact frequency, cu, with degree of separation, η, for various values of the parameter mu. 
in the figure.
When mu < 0, the contact frequency of unvaccinated people decreases with increasing η. This represents a separation policy that excludes unvaccinated people from public spaces while also isolating them from themselves to some degree. For example, a pair of intermediate values of mu < 0 and 0 < η < 1 could represent a separation policy (such as with “vaccination passports”) that excludes unvaccinated people from recreational venues like restaurants and cinemas but not more essential services such as grocery stores and hospitals. On the other hand, in the extreme case of mu = −1 and η = 1, cu = 0 such that unvaccinated people are completely separated and isolated, having no contacts with anyone.
When mu > 0, the contact frequency of unvaccinated people increases with increasing η. This represents a separation policy that compounds unvaccinated people by placing them together in close quarters, for example in designated facilities for returning unvaccinated travelers.
Via the two parameters mu and η, the model therefore spans the full range of contact frequencies, from cu = 0 in the isolating extreme of mu = −1 and η = 1 to arbitrarily high contact frequency in the compounding extreme with η = 1 and a mu > 0.
The impact of separation on the vaccinated sub-population is independently controlled via the parameters 
and mv.
A1.3 Intra- and inter-population contact frequencies
The contact frequency ci represents the frequency of contacts for a person in group i (either u or v) irrespective of whether the contact partners belong to group u or v. We use the term “population-specific contact frequency” to refer to ci.
Given the probability fij that a person of type i has a contact with a person of type j (see Eq. A2), one can define the intra-population contact frequencies cii and the inter-population contact frequencies cij as follows:
Eqs. A4 also show that ci = cii + cij.
Fig. A1.2 shows how the population-specific contact frequencies cu and cv, the intra-population contact frequencies cuu and cvv, and the inter-population contact frequencies cuv and cvu vary with η, for the six values of mu explored in the main text and mv = 0, and for three values of the population fraction of vaccinated people, Pv.
Fig. 1 of the main text uses a value of Pv = 0.8, corresponding to the right column of panels in Fig. A1.2. Model results for Pv = 0.25 (left column of Fig. A1.2) and Pv = 0.5 (middle column of Fig. A1.2) are included in Appendix 2.
Normalized population-specific contact frequencies cu and cv, intra-population contact frequencies cuu and cvv, and inter-population contact frequencies cuv and cvu, as functions of η, for the six values of mu explored in the main text (each row of panels corresponds to a different value of mu), for mv = 0, and for Pv = 0.2 (left column of panels), Pv = 0.5 (middle column), and Pv = 0.8 (right column).
A1.4: Quantities calculated from simulation results
The attack rate among the vaccinated population, Av, is defined as the proportion of initially-susceptible vaccinated people who become infected during the epidemic (Eq. 1 of the main text, reproduced below):
where Sv(t0) is the number of susceptible vaccinated people at the beginning of the epidemic and Sv(tf) is the number of susceptible vaccinated people remaining once there are no longer any infectious people in the entire (vaccinated and unvaccinated) population. Au is defined equivalently, for the unvaccinated.
The quantity Bv is equal to the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people, i.e.:
Similarly, Bu represents the share of infections among unvaccinated people that were due to contacts with infectious vaccinated people, and is defined in the same way as Bv, (interchanging the v and u subscripts in Eq. A6).
Appendix 2: Supplementary figures
This is Appendix 2 of the pre-print “Viral respiratory epidemic modeling of societal segregation based on vaccination status” by J. Hickey & D.G. Rancourt, uploaded to https://www.medrxiv.org/ on 2023-10-31 (Version 5).
A2.1: Epidemic curves
Population fraction of infectious individuals as a function of time, for the vaccinated (solid lines) and unvaccinated (dashed lines) populations, for parameters Pv = 0.8, VE = 0.2, NI = 0.2, mv = 0, 
, γ = 73. Each panel (a-d) shows a different value of η, and each coloured line is for a different value of mu, as indicated in the legend.
Same as Fig. A2.1, except that 
.
Same as Fig. A2.1, except that 
.
Same as Fig. A2.1, except that 
.
Same as Fig. A2.1, except that 
.
A2.2: Epidemic outcomes for different values of Pv
A2.2.1: VE = 0.8
Av, Au, Bv, and Bu as functions of η. Each row of panels corresponds to a choice of 
, and each coloured line to a choice of mu as indicated in the legends. For reference, in a single-population (no vaccination) model, the corresponding R0 values for rows a-e of the figure are 8.2, 6.0, 4.1, 3.4 and 2.7, respectively. Parameter values Pv = 0.1, VE = 0.8, NI = 0.2, mv = 0, γ = 73.
Same as Fig. A2.6, except that Pv = 0.25.
Same as Fig. A2.6, except that Pv = 0.5.
Same as Fig. A2.6, except that Pv = 0.6.
Same as Fig. A2.6, except that Pv = 0.7.
Same as Fig. A2.6, except that Pv = 0.9.
Same as Fig. A2.6, except that Pv = 0.95.
Same as Fig. A2.6, except that Pv = 0.99.
A2.2.2: VE = 0.4
Av, Au, Bv, and Bu as functions of η. Each row of panels corresponds to a choice of 
, and each coloured line to a choice of mu as indicated in the legends. For reference, in a single-population (no vaccination) model, the corresponding R0 values for rows a-e of the figure are 8.2, 6.0, 4.1, 3.4 and 2.7, respectively. Parameter values Pv = 0.1, VE = 0.4, NI = 0.2, mv = 0, γ = 73.
Same as Fig. A2.15, except that Pv = 0.25.
Same as Fig. A2.15, except that Pv = 0.5.
Same as Fig. A2.15, except that Pv = 0.6.
Same as Fig. A2.15, except that Pv = 0.7.
Same as Fig. A2.15, except that Pv = 0.8.
Same as Fig. A2.15, except that Pv = 0.9.
Same as Fig. A2.15, except that Pv = 0.95.
Same as Fig. A2.15, except that Pv = 0.99.
A2.3: Epidemic outcomes for different values of mv
Av, Au, Bv, and Bu as functions of η. Each row of panels corresponds to a choice of 
, and each coloured line to a choice of mu as indicated in the legends. For reference, in a single-population (no vaccination) model, the corresponding R0 values for rows a-e of the figure are 8.2, 6.0, 4.1, 3.4 and 2.7, respectively. Parameter values Pv = 0.8, VE = 0.8, NI = 0.2, mv = 0, γ = 73. This figure is for the same parameters as Fig. 1 of the main text, such that the Av, Au, and Bv columns in this figure are reproductions of the same columns in Fig. 1 of the main text.
Same as Fig. A2.24, except that mv = 1.
Same as Fig. A2.24, except that mv = -1.
A2.4: Epidemic outcomes for different values of 
A2.4.1: Fixed , vary
, mu = mv = 0, γu = γv = 73, VE = 0.8, NI = 0.2, various choices of 
(see legend within left-column panels), showing Av, Au, Bv,, and Bu as functions of η. 
is fixed for each row in the figure and decreases moving down the rows.
Same as Fig. A2.27, except that mu = 1.
Same as Fig. A2.27, except that mu = -1.
A2.4.2: Fixed weighted sum
, mu = mv = 0, γu = γv = 73, VE = 0.8, NI = 0.2, various choices of 
(see legend within left-column panels), showing Av, Au, Bv,, and Bu as functions of η. The weighted sum 
is fixed for each row in the figure and decreases moving down the rows.
Same as Fig. A.30, except that mu = 1.
Same as Fig. A.30, except that mu = -1.
Footnotes
Formatting, new references, adjustments to the text, new title.
References
