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

delete 'in sputum'

Browse source
This commit is contained in:
Andrejh 2021-02-16 22:14:18 +01:00
parent 5d675d1d1a
commit 7fd295930b
1 changed files with 5 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
cara/montecarlo.py
View file

@ -542,7 +542,7 @@ def plot_pi_vs_viral_load(baselines: typing.Union[MCExposureModel, typing.List[M
plt.title(title)
plt.ylabel('Percentage probability of infection')
plt.xticks(ticks=[i for i in range(3, 13)], labels=['$10^{' + str(i) + '}$' for i in range(3, 13)])
plt.xlabel('Viral load in sputum')
plt.xlabel('Viral load')
if labels is not None:
plt.legend(labels)
@ -715,7 +715,7 @@ def compare_infection_probabilities_vs_viral_loads(baseline1: MCExposureModel, b
plt.plot(viral_loads, mean_ratios)
plt.ylabel(f"Ratio of mean P(i) values - P(i|qID = {baseline1.concentration_model.infected.qid}) / P(i|qID = {baseline2.concentration_model.infected.qid})")
plt.xlabel("Viral load in sputum")
plt.xlabel("Viral load")
plt.xticks(ticks=[i for i in range(3, 11)], labels=['$10^{' + str(i) + '}$' for i in range(3, 11)])
plt.show()
@ -723,7 +723,7 @@ def compare_infection_probabilities_vs_viral_loads(baseline1: MCExposureModel, b
plt.plot(viral_loads, p2_means)
plt.ylim(0, 100)
plt.ylabel("Percentage probability of infection")
plt.xlabel("Viral load in sputum")
plt.xlabel("Viral load")
plt.legend([f'qID = {baseline1.concentration_model.infected.qid}',
f'qID = {baseline2.concentration_model.infected.qid}'])
plt.xticks(ticks=[i for i in range(3, 11)], labels=['$10^{' + str(i) + '}$' for i in range(3, 11)])
@ -841,3 +841,5 @@ for model in exposure_models:
plt.title(f'Probability of infection in baseline case - {"English" if model.concentration_model.infected.qid == 60 else "Original"} variant')
plt.yticks([], [])
plt.show()