*** Wartungsfenster jeden ersten Mittwoch vormittag im Monat ***

Skip to content
Snippets Groups Projects
Commit 20013df3 authored by Patrick Kappl's avatar Patrick Kappl
Browse files

Change HDF5 file and dataset names

parent 7ce8fc70
No related branches found
No related tags found
No related merge requests found
import numpy as np
import h5py
import time
# TODO: Exceptions!
......@@ -8,7 +9,7 @@ import h5py
# TODO: Rewrite docstring with argument descriptions
def do_jackknife_estimation(input_sample_generator, n_samples,
transformation_function,
output_file_name="jackknife.hdf5"):
output_file_name=""):
"""Do a jackknife estimation for the given samples and function.
Estimate and return the sample mean, variance and standard deviation
......@@ -25,6 +26,7 @@ def do_jackknife_estimation(input_sample_generator, n_samples,
print("At least two samples are needed for the jackknife estimation.")
exit()
date_and_time = time.strftime("%Y-%m-%d_%H-%M-%S")
gen = input_sample_generator
f = transformation_function
samples = gen()
......@@ -55,10 +57,12 @@ def do_jackknife_estimation(input_sample_generator, n_samples,
output_variance /= (n_samples - 1)
output_standard_deviation = np.lib.scimath.sqrt(output_variance)
if output_file_name == "":
output_file_name = "jackknife_" + date_and_time + ".hdf5"
out_file = h5py.File(output_file_name, "w")
out_file.create_dataset("estimated_mean", data=output_mean)
out_file.create_dataset("estimated_variance", data=output_variance)
out_file.create_dataset("estimated_standard_deviation",
out_file.create_dataset("mean", data=output_mean)
out_file.create_dataset("variance", data=output_variance)
out_file.create_dataset("standard_deviation",
data=output_standard_deviation)
out_file.create_dataset("transformed_input_mean", data=transformed_mean)
out_file.close()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment