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

usage example

Browse source
This commit is contained in:
markus 2021-02-11 17:15:59 +01:00
parent a05aad27f5
commit ef7080e4bd
1 changed files with 8 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
cara/montecarlo.py
View file

@ -669,10 +669,12 @@ exposure_models = [MCExposureModel(
)
) for e in (False, True)]
plot_pi_vs_viral_load(exposure_models, labels=['Without masks', 'With masks'])
plot_pi_vs_viral_load(exposure_models[0])
for model in exposure_models:
present_model(model.concentration_model)
plt.hist(model.infection_probability(), bins=200)
plt.xlabel('Percentage probability of infection')
plt.show()
# for model in exposure_models:
# present_model(model.concentration_model, title=f'Model summary - {"English" if model.concentration_model.infected.english_variant else "Original"} variant')
# plt.hist(model.infection_probability(), bins=200)
# plt.xlabel('Percentage probability of infection')
# plt.title(f'Probability of infection in baseline case - {"English" if model.concentration_model.infected.english_variant else "Original"} variant')
# plt.yticks([], [])
# plt.show()