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

incorporate print_qr_info in present_model

Browse source
This commit is contained in:
markus 2021-02-08 15:52:59 +01:00
parent 18c24ff063
commit f5230199cb
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/montecarlo.py
View file

@ -388,8 +388,8 @@ def logscale_hist(x: typing.Iterable, bins: int) -> None:
plt.xscale('log')
def print_qr_info(qr_values: np.ndarray) -> None:
log_qr = np.log10(qr_values)
def print_qr_info(log_qr: np.ndarray) -> None:
qr_values = 10 ** log_qr
print(f"MEAN of log_10(qR) = {np.mean(log_qr)}\n"
f"MEAN of qR = {np.mean(qr_values)}")
@ -464,7 +464,7 @@ def present_model(model: MCConcentrationModel, bins: int = 200) -> None:
std_patch = patches.Patch(color='lightgrey', linestyle='dashed', label='Standard deviations')
fig.legend(handles=(mean_patch, std_patch, median_patch))
# TODO: call print_qr_info
print_qr_info(qRs)
plt.show()
@ -519,7 +519,7 @@ models = [MCExposureModel(
present_model(models[0].concentration_model)
original_pi, english_pi = [model.infection_probability() for model in models]
print(f"Median(P_i) - Original: {'{:.2f}'.format(np.median(original_pi))}%\n"
print(f"\nMedian(P_i) - Original: {'{:.2f}'.format(np.median(original_pi))}%\n"
f"Mean(P_i) - Original: {'{:.2f}'.format(np.mean(original_pi))}%\n\n"
f"Median(P_i) - English: {'{:.2f}'.format(np.median(english_pi))}%\n"
f"Mean(P_i) - English: {'{:.2f}'.format(np.mean(english_pi))}%\n")