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

update usage example

Browse source
This commit is contained in:
markus 2021-02-16 16:12:56 +01:00
parent edcefab54e
commit f6d0f4ce24
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
cara/montecarlo.py
View file

@ -833,8 +833,13 @@ exposure_models = [MCExposureModel(
) for qid in (100, 60)]
rs = [model.expected_new_cases() for model in large_population_baselines]
compare_infection_probabilities_vs_viral_loads(*exposure_models)
print(f"R0 - original variant:\t{np.mean(rs[0])}")
print(f"R0 - english variant:\t{np.mean(rs[1])}")
print(f"Ratio between R0's:\t\t{np.mean(rs[1]) / np.mean(rs[0])}")
# compare_infection_probabilities_vs_viral_loads(*exposure_models)
# present_model(exposure_models[0].concentration_model)