From a4b1e5adf2d0a46dfcc2c82254fb2caff5215a0f Mon Sep 17 00:00:00 2001 From: Andrejh Date: Mon, 6 Sep 2021 15:36:54 +0200 Subject: [PATCH 1/2] legend edits --- cara/montecarlo.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 0a12f297..207eeb78 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -121,10 +121,10 @@ def exposure_model_from_vl_breathing_cn(): ax.plot(viral_loads, er_means, color=cmap.to_rgba( cn, alpha=0.75), linewidth=1, ls='--') plt.text(viral_loads[int(len(viral_loads)*0.9)], 10**4.2, - r"$\mathbf{c_{n,B}=0.06}$", color=cmap.to_rgba(cn), size='small') + r"$\mathbf{c_{n,B}=0.06}$", color=cmap.to_rgba(cn), fontesize=12) fig.colorbar(cmap, ticks=[0.01, 0.1, 0.25, 0.5], - label="Particle emission concentration, ${c_{n,B}}$") + label="Particle emission concentration, ${c_{n,B}}$", fontesize=14) ax.set_yscale('log') ############# Coleman ############# @@ -186,7 +186,7 @@ def exposure_model_from_vl_talking(): build_talking_legend(fig) ############ Plot ############ - plt.title('Exhaled virions while talking for 15min', + plt.title('Exhaled virions while speaking (vocalization) for 15min', fontsize=16, fontweight="bold") plt.ylabel( 'Aerosol viral load, $\mathrm{vl_{out}}$\n(RNA copies)', fontsize=14) @@ -231,10 +231,10 @@ def exposure_model_from_vl_talking_cn(): ax.plot(viral_loads, er_means, color=cmap.to_rgba( cn, alpha=0.75), linewidth=1, ls='--') plt.text(viral_loads[int(len(viral_loads)*0.93)], 10**5.5, - r"$\mathbf{c_{n,L}=0.2}$", color=cmap.to_rgba(cn), size='small') + r"$\mathbf{c_{n,L}=0.2}$", color=cmap.to_rgba(cn), fontesize=12) fig.colorbar(cmap, ticks=[0.01, 0.5, 1.0, 2.0], - label="Particle emission concentration, ${c_{n,L}}$") + label="Particle emission concentration, ${c_{n,L}}$", fontesize=14) ax.set_yscale('log') ############# Coleman ############# @@ -244,7 +244,7 @@ def exposure_model_from_vl_talking_cn(): build_talking_legend(fig) ############ Plot ############ - plt.title('Exhaled virions while talking for 15min', + plt.title('Exhaled virions while speaking (vocalization)\nfor 15min', fontsize=16, fontweight="bold") plt.ylabel( 'Aerosol viral load, $\mathrm{vl_{out}}$\n(RNA copies)', fontsize=14) From 16e705f5bb0d84ddafcff67c9778a6199f7e5822 Mon Sep 17 00:00:00 2001 From: Andrejh Date: Mon, 6 Sep 2021 17:00:30 +0200 Subject: [PATCH 2/2] increase sample size + add means of log(er) --- cara/results_paper.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/cara/results_paper.py b/cara/results_paper.py index 7b13d28f..4c5846fe 100644 --- a/cara/results_paper.py +++ b/cara/results_paper.py @@ -14,7 +14,7 @@ import matplotlib as mpl ######### Plot material ######### -SAMPLE_SIZE = 50000 +SAMPLE_SIZE = 500000 viral_loads = np.linspace(2, 12, 600) ############# Markers (for legend) ############# @@ -394,14 +394,16 @@ def build_breathing_legend(fig): label.set_position((-3*width, 0)) -def print_er_info(er: np.array): +def print_er_info(er: np.array, log_er: np.array): """ Prints statistical parameters of a given distribution of ER-values :param er: A numpy-array of the ER-values :return: Nothing, parameters are printed """ print(f"MEAN of ER = {np.mean(er)}\n" + f"MEAN of log ER = {np.mean(log_er)}\n" f"SD of ER = {np.std(er)}\n" + f"SD of log ER = {np.std(log_er)}\n" f"Median of ER = {np.quantile(er, 0.5)}\n") print(f"Percentiles of ER:")