vaccinated people with prior infected access natural mapping
@@ -493,6 +493,7 @@ def check_immunity(people, strain, sus=True, inds=None):
was_inf_same = cvu.true((people.recovered_strain == strain) & (people.t >= date_rec)) # Had a previous exposure to the same strain, now recovered
was_inf_diff = np.setdiff1d(was_inf, was_inf_same) # Had a previous exposure to a different strain, now recovered
is_sus_vacc = np.intersect1d(is_sus, is_vacc) # Susceptible and vaccinated
+ is_sus_vacc = np.setdiff1d(is_sus_vacc, was_inf) # Susceptible, vaccinated without prior infection
is_sus_was_inf_same = np.intersect1d(is_sus, was_inf_same) # Susceptible and being challenged by the same strain
is_sus_was_inf_diff = np.intersect1d(is_sus, was_inf_diff) # Susceptible and being challenged by a different strain
