mrestimator_triallength/README.md

Mr. Estimator Triallength Dependence

This is a collection of scripts that produce Fig.~4 of the paper on Mr. Estimator.

Dependencies

conda install h5py numpy matplotlib tqdm
pip3 install -U 'mrestimator[full]'

Outline

1. create a realization of a branching process of certain length, timescale and number of trials. This is done in run/triallength.py
2. calculate the correlation coefficients. This is done in run/triallength.py.
   • Because the first two steps take significant time, we attach a copy of the data in ./dat/ that contains precalculated correlation coefficients.
   • The run/triallength.py script takes an integer argument (think thread or job id) because we ran this on the cluster. Every id creates it's own hdf5 output file.
   • Merge the individual files with run/merge_hdf5.py to obtain data files similar to the ones provided here.
3. fit the exponentials. This is done with plt/plot_merged.py.

Note that the realizations for long trials (large targettau and targetlength in run/triallength.py) will take on the order of days computation and multiple gigabytes of ram. Adjust as needed.