From f6d0f4ce244b7e6c2f9030249ff65d8b885bdcb0 Mon Sep 17 00:00:00 2001 From: markus Date: Tue, 16 Feb 2021 16:12:56 +0100 Subject: [PATCH] update usage example --- cara/montecarlo.py | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index b993ea02..2d4aaef3 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -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)