# 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.