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

add function for printing distribution information

Browse source
This commit is contained in:
markus 2021-01-20 15:00:44 +01:00
parent b67ffc5e2b
commit da718df7c5
1 changed files with 9 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
cara/montecarlo.py
View file

@ -67,3 +67,12 @@ def logscale_hist(x: typing.Iterable, bins: int) -> None:
plt.hist(x, bins=logscale_bins) plt.hist(x, bins=logscale_bins)
plt.xscale('log') plt.xscale('log')
plt.show() plt.show()
def print_qr_info(qr_values: np.ndarray) -> None:
log_qr = np.log10(qr_values)
print(f"MEAN of log_10(qR) = {np.mean(log_qr)}\n"
f"MEAN of qR = {np.mean(qr_values)}")
for quantile in (0.01, 0.05, 0.25, 0.50, 0.75, 0.95, 0.99):
print(f"qR_{quantile} = {np.quantile(qr_values, quantile)}")