pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

change exposure_models to illustrate point

Browse source
This commit is contained in:
markus 2021-02-15 17:44:11 +01:00
parent fe00ce1b42
commit a3dffcf1e7
1 changed files with 11 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
cara/montecarlo.py
View file

@ -776,11 +776,11 @@ exposure_models = [MCExposureModel(
infected=MCInfectedPopulation(
number=1,
presence=models.SpecificInterval(((0, 4), (5, 9))),
masked=m,
masked=False,
virus=MCVirus(halflife=1.1),
expiratory_activity=1,
samples=2000000,
qid=300,
qid=qid,
breathing_category=1,
english_variant=False
)
@ -791,14 +791,18 @@ exposure_models = [MCExposureModel(
activity=models.Activity.types['Seated'],
mask=models.Mask.types['Type I']
)
) for m in (False, True)]
) for qid in (300, 100)]
compare_infection_probabilities_vs_viral_loads(*exposure_models)
# present_model(exposure_models[0].concentration_model)
plot_pi_vs_qid(fixed_vl_exposure_models, labels=['Viral load = $10^{' + str(i) + '}$' for i in range(6, 11)],
qid_min=5, qid_max=2000, qid_samples=200)
# plot_pi_vs_qid(fixed_vl_exposure_models, labels=['Viral load = $10^{' + str(i) + '}$' for i in range(6, 11)],
# qid_min=5, qid_max=2000, qid_samples=200)
plot_pi_vs_qid(fixed_vl_exposure_models, labels=['Viral load = $10^{' + str(i) + '}$' for i in range(6, 11)],
qid_min=20, qid_max=200, qid_samples=100)
# plot_pi_vs_qid(fixed_vl_exposure_models, labels=['Viral load = $10^{' + str(i) + '}$' for i in range(6, 11)],
# qid_min=100, qid_max=400, qid_samples=100)
# plot_pi_vs_viral_load(exposure_models, labels=['Without masks', 'With masks'])