Skip to content
Snippets Groups Projects
Commit 5674a740 authored by Kipp Cannon's avatar Kipp Cannon
Browse files

rate_estimation.py: remove overflow bins from rate binnings

parent 0106043a
No related branches found
No related tags found
No related merge requests found
......@@ -340,7 +340,7 @@ def binned_rates_from_samples(samples):
"""
lo, hi = math.floor(samples.min()), math.ceil(samples.max())
nbins = int(math.sqrt(len(samples)) / 40.)
binnedarray = rate.BinnedArray(rate.NDBins((rate.LogarithmicPlusOverflowBins(lo if lo !=0. else samples.min(), hi, nbins),)))
binnedarray = rate.BinnedArray(rate.NDBins((rate.LogarithmicBins(lo, hi, nbins),)))
for sample in samples:
binnedarray[sample,] += 1.
rate.filter_array(binnedarray.array, rate.gaussian_window(5), use_fft = False)
......
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