def test_plot_with_cvpeople(create_pop, do_show=False, do_save=False):
"""
Test plotting method works on covasim.people.People object.
Notes:
With this pop type, you will need to supply more information to
tell the method where to look for expected data.
"""
sp.logger.info(
"Test that the age comparison plotting method works on cv.people.People and plotting styles can be easily updated."
)
pop = create_pop
popdict = pop.to_dict()
cvpopdict = cv.make_synthpop(population=popdict,
community_contacts=2) # array based
# Actually create the people
people_pars = dict(
pop_size=pars.n,
beta_layer={k: 1.0
for k in 'hswcl'
}, # Since this is used to define hat layers exist
beta=
1.0, # TODO: this is required for plotting (people.plot()), but shouldn't be (in covasim)
)
people = cv.People(people_pars,
strict=False,
uid=cvpopdict['uid'],
age=cvpopdict['age'],
sex=cvpopdict['sex'])
kwargs = sc.objdict(sc.mergedicts(pars, pop.loc_pars))
kwargs.datadir = sp.settings.datadir
kwargs.figname = f"test_ages_{kwargs.location}_cvpeople"
kwargs.do_show = do_show
kwargs.do_save = do_save
# modify some plotting styles
kwargs.color_1 = '#9966cc'
kwargs.color_2 = 'indigo'
kwargs.markersize = 4.5
fig, ax = sp.plot_ages(people, **kwargs)
# fig, ax = sp.plot_ages(people) # to plot without extra information
assert isinstance(fig, mplt.figure.Figure), 'Check failed.'
print('Check passed. Figure made.')
return fig, ax, people
def test_synthpops():
sim = cv.Sim(use_waning=True,
**sc.mergedicts(base_pars,
dict(pop_size=5000, pop_type='synthpops')))
sim.popdict = cv.make_synthpop(sim,
with_facilities=True,
layer_mapping={'LTCF': 'f'})
sim.reset_layer_pars()
# Vaccinate 75+, then 65+, then 50+, then 18+ on days 20, 40, 60, 80
sim.vxsubtarg = sc.objdict()
sim.vxsubtarg.age = [75, 65, 50, 18]
sim.vxsubtarg.prob = [.05, .05, .05, .05]
sim.vxsubtarg.days = subtarg_days = [20, 40, 60, 80]
pfizer = cv.vaccinate(days=subtarg_days,
vaccine='pfizer',
subtarget=vacc_subtarg)
sim['interventions'] += [pfizer]
sim.run()
return sim
def make_population(pop_size,
rand_seed=1,
max_pop_seeds=None,
do_save=True,
popfile=None,
cohorting=True,
n_brackets=20,
community_contacts=20,
**kwargs):
'''
Generate the synthpops population.
Args:
pop_size (int): number of people in the model
rand_seed (int): random seed to use for generating the population
max_pop_seeds (int): if supplied, take the random seed as modulus of this to limit number of populations generated
do_save (bool): whether to save the population
popfile (str): if so, where to save it to
cohorting (bool): whether to use cohorting
n_brackets (int): whether to use 16- or 20-bracket age bins
community_contacts (int): how many community contacts there are
kwargs (dict): passed to sp.make_population()
'''
sp.set_nbrackets(
n_brackets) # Essential for getting the age distribution right
pars = sc.objdict(
n=pop_size,
rand_seed=rand_seed,
with_facilities=True,
use_two_group_reduction=True,
average_LTCF_degree=20,
ltcf_staff_age_min=20,
ltcf_staff_age_max=60,
with_school_types=True,
average_class_size=20,
inter_grade_mixing=0.1,
average_student_teacher_ratio=20,
average_teacher_teacher_degree=3,
teacher_age_min=25,
teacher_age_max=75,
with_non_teaching_staff=True,
# if with_non_teaching_staff is False, but generate is True, then ,average_all_staff_ratio should be average_student_teacher_ratio or 0
average_student_all_staff_ratio=11,
average_additional_staff_degree=20,
staff_age_min=20,
staff_age_max=75,
)
pars.update(kwargs) # Update any parameters
# For reference re: school_types
# school_mixing_type = 'random' means that students in the school have edges randomly chosen from other students, teachers, and non teaching staff across the school. Students, teachers, and non teaching staff are treated the same in terms of edge generation.
# school_mixing_type = 'age_clustered' means that students in the school have edges mostly within their own age/grade, with teachers, and non teaching staff. Strict classrooms are not generated. Teachers have some additional edges with other teachers.
# school_mixing_type = 'age_and_class_clustered' means that students are cohorted into classes of students of the same age/grade with at least 1 teacher, and then some have contact with non teaching staff. Teachers have some additional edges with other teachers.
if cohorting:
strategy = 'clustered' # students in pre-k, elementary, and middle school are cohorted into strict classrooms
pars.school_mixing_type = {
'pk': 'age_and_class_clustered',
'es': 'age_and_class_clustered',
'ms': 'age_and_class_clustered',
'hs': 'random',
'uv': 'random'
}
else:
strategy = 'normal'
pars.school_mixing_type = {
'pk': 'age_clustered',
'es': 'age_clustered',
'ms': 'age_clustered',
'hs': 'random',
'uv': 'random'
}
if popfile is None:
popfile = os.path.join(
'inputs', f'kc_{strategy}_{int(pars.n)}_seed{pars.rand_seed}.ppl')
T = sc.tic()
print('Making population...')
# Make the population
population = sp.make_population(**pars)
# Convert to a popdict
popdict = cv.make_synthpop(population=sc.dcp(population),
community_contacts=community_contacts)
school_ids = [None] * int(pop_size)
teacher_flag = [False] * int(pop_size)
staff_flag = [False] * int(pop_size)
student_flag = [False] * int(pop_size)
school_types = {'pk': [], 'es': [], 'ms': [], 'hs': [], 'uv': []}
school_type_by_person = [None] * int(pop_size)
schools = dict()
for uid, person in population.items():
if person['scid'] is not None:
school_ids[uid] = person['scid']
school_type_by_person[uid] = person['sc_type']
if person['scid'] not in school_types[person['sc_type']]:
school_types[person['sc_type']].append(person['scid'])
if person['scid'] in schools:
schools[person['scid']].append(uid)
else:
schools[person['scid']] = [uid]
if person['sc_teacher'] is not None:
teacher_flag[uid] = True
elif person['sc_student'] is not None:
student_flag[uid] = True
elif person['sc_staff'] is not None:
staff_flag[uid] = True
popdict['school_id'] = np.array(school_ids)
popdict['schools'] = schools
popdict['teacher_flag'] = teacher_flag
popdict['student_flag'] = student_flag
popdict['staff_flag'] = staff_flag
popdict['school_types'] = school_types
popdict['school_type_by_person'] = school_type_by_person
assert sum(
popdict['teacher_flag']
), 'Uh-oh, no teachers were found: as a school analysis this is treated as an error'
assert sum(
popdict['student_flag']
), 'Uh-oh, no students were found: as a school analysis this is treated as an error'
# Actually create the people
people_pars = dict(
pop_size=pars.n,
beta_layer={k: 1.0
for k in 'hswcl'
}, # Since this is used to define hat layers exist
beta=
1.0, # TODO: this is required for plotting (people.plot()), but shouldn't be
)
people = cv.People(people_pars,
strict=False,
uid=popdict['uid'],
age=popdict['age'],
sex=popdict['sex'],
contacts=popdict['contacts'],
school_id=popdict['school_id'],
schools=popdict['schools'],
school_types=popdict['school_types'],
student_flag=popdict['student_flag'],
teacher_flag=popdict['teacher_flag'],
staff_flag=popdict['staff_flag'],
school_type_by_person=popdict['school_type_by_person'])
if do_save:
print(f'Saving to "{popfile}"...')
sc.saveobj(popfile, people)
sc.toc(T)
print('Done')
return people
sub2 = cv.Sim(pars)
sub2.update_pars(beta_layer={'a': 0.5})
sub2.run()
assert sub2['beta_layer']['a'] == 0.5
with pytest.raises(sc.KeyNotFoundError): # q is not ok
sub2.update_pars(beta_layer={'q': 0.5})
sub2.run()
sc.heading('SynthPops test')
ssp = cv.Sim(pars, pop_type='synthpops')
ssp.run()
sc.heading('SynthPops + LTCF test')
sf = cv.Sim(pop_size=2000, pop_type='synthpops')
sf.popdict = cv.make_synthpop(sf,
with_facilities=True,
layer_mapping={'LTCF': 'f'})
cv.save('synth.pop', sf.popdict)
with pytest.raises(ValueError): # Layer key mismatch, f is not initialized
sf.run()
sf.reset_layer_pars()
sf.run()
assert 'f' in sf.layer_keys()
assert sf['beta_layer']['f'] == 1.0
sc.heading('Population save test')
sps = cv.Sim(pars, pop_type='hybrid', popfile='test.pop', save_pop=True)
sps.run()
sc.heading('Population load test')
spl = cv.Sim(pars, pop_type='hybrid', popfile='test.pop', load_pop=True)
'''
Script for testing integration with new SynthPops features -- namely, long-term
care facilities.
'''
import covasim as cv
which = ['simple', 'complex'][1]
if which == 'simple':
sim = cv.Sim(pop_size=5000, pop_type='synthpops')
sim.popdict = cv.make_synthpop(sim,
with_facilities=True,
layer_mapping={'LTCF': 'f'})
sim.reset_layer_pars()
sim.run()
else:
pars = dict(
pop_size=20000,
pop_type='synthpops',
n_days=120,
)
sims = []
for ltcf in [False, True]:
print(f'Running LTCF {ltcf}')
sim = cv.Sim(pars)
sim.popdict = cv.make_synthpop(sim,
