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

bug in print_er_info arguments

Browse source
This commit is contained in:
Luis Aleixo 2021-09-07 09:56:42 +02:00
parent 5f9091f18b
commit b820b3d889
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
cara/plot_output.py
View file

@ -23,7 +23,7 @@ print('\n<<<<<<<<<<< Vlout for Talking, seated with chosen Cn,L >>>>>>>>>>>')
# Exhaled virions while breathing according to BLO model, seated #
print('\n<<<<<<<<<<< Vlout for Breathing, seated with chosen Cn,B >>>>>>>>>>>')
exposure_model_from_vl_breathing_cn()
#exposure_model_from_vl_breathing_cn()
print('\n')
############ Statistical Data ############
@ -34,7 +34,7 @@ exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE)
emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present(cn_B=0.06, cn_L=0.2)
breathing_er = [np.log10(er) for er in emission_rate]
print('\n<<<<<<<<<<< Breathing model statistics >>>>>>>>>>>')
print_er_info(breathing_er)
print_er_info(emission_rate, breathing_er)
############ Speaking model ############
exposure_mc = model_scenarios_paper.speaking_exposure()
@ -42,7 +42,7 @@ exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE)
emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present(cn_B=0.06, cn_L=0.2)
speaking_er = [np.log10(er) for er in emission_rate]
print('\n<<<<<<<<<<< Speaking model statistics >>>>>>>>>>>')
print_er_info(speaking_er)
print_er_info(emission_rate, speaking_er)
############ Shouting model ############
exposure_mc = model_scenarios_paper.shouting_exposure()
@ -50,7 +50,7 @@ exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE)
emission_rate = exposure_model.concentration_model.infected.emission_rate_when_present(cn_B=0.06, cn_L=0.2)
shouting_er = [np.log10(er) for er in emission_rate]
print('\n<<<<<<<<<<< Shouting model statistics >>>>>>>>>>>')
print_er_info(shouting_er)
print_er_info(emission_rate, shouting_er)
############ Plots with viral loads and emission rates