From 90ae11ebc4a13612bf493c28c336f9cd16eca6c2 Mon Sep 17 00:00:00 2001 From: Andrejh Date: Thu, 19 Aug 2021 11:49:27 +0200 Subject: [PATCH] normalize data from per hour to per 30min (breathing) and per 15min (talking) --- cara/model_scenarios.py | 16 ++++++++++------ 1 file changed, 10 insertions(+), 6 deletions(-) diff --git a/cara/model_scenarios.py b/cara/model_scenarios.py index c05a91cb..eb6af442 100644 --- a/cara/model_scenarios.py +++ b/cara/model_scenarios.py @@ -74,17 +74,19 @@ def exposure_model_from_vl_talking(viral_loads): ), ) exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE) - emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present() + # divide by 4 to have in 15min (quarter of an hour) + emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present()/4 er_means.append(np.mean(emission_rate)) er_medians.append(np.median(emission_rate)) lower_percentiles.append(np.quantile(emission_rate, 0.01)) upper_percentiles.append(np.quantile(emission_rate, 0.99)) + # divide by 4 to have in 15min (quarter of an hour) build_plot(viral_loads, er_means, - lower_percentiles, upper_percentiles, coleman_etal_vl_talking, coleman_etal_er_talking, milton_vl, milton_er, yann_vl, yann_er) + lower_percentiles, upper_percentiles, coleman_etal_vl_talking, [x/4 for x in coleman_etal_er_talking]) ############ Plot ############ - plt.title('Exhaled virions while talking for 1h', + plt.title('Exhaled virions while talking for 15min', fontsize=16, fontweight="bold") plt.ylabel( 'Aerosol viral load, $\mathrm{vl_{out}}$\n(RNA copies)', fontsize=14) @@ -126,17 +128,19 @@ def exposure_model_from_vl_breathing(viral_loads): ), ) exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE) - emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present() + # divide by 2 to have in 30min (half an hour) + emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present() /2 er_means.append(np.mean(emission_rate)) er_medians.append(np.median(emission_rate)) lower_percentiles.append(np.quantile(emission_rate, 0.01)) upper_percentiles.append(np.quantile(emission_rate, 0.99)) + # divide by 2 to have in 30min (half an hour) build_plot(viral_loads, er_means, - lower_percentiles, upper_percentiles, coleman_etal_vl_breathing, coleman_etal_er_breathing, milton_vl, milton_er, yann_vl, yann_er) + lower_percentiles, upper_percentiles, coleman_etal_vl_breathing, [x/2 for x in coleman_etal_er_breathing], milton_vl, [x/2 for x in milton_er], yann_vl, [x/2 for x in yann_er]) ############ Plot ############ - plt.title('Exhaled virions while breathing for 1h', + plt.title('Exhaled virions while breathing for 30 min', fontsize=16, fontweight="bold") plt.ylabel( 'Aerosol viral load, $\mathrm{vl_{out}}$\n(RNA copies)', fontsize=14)