diff --git a/basicrta/gibbs.py b/basicrta/gibbs.py index ebed117..ce0f907 100644 --- a/basicrta/gibbs.py +++ b/basicrta/gibbs.py @@ -844,8 +844,9 @@ def estimate_tau(self): taus = 1 / rp.rates[rp.labels == index] ci = confidence_interval(taus) citaus = taus[(taus > ci[0]) & (taus < ci[1])] - bins = np.exp(np.linspace(np.log(citaus.min()), np.log(citaus.max()), - 20)) + #bins = np.exp(np.linspace(np.log(citaus.min()), np.log(citaus.max()), + # 20)) + bins = 15 h = np.histogram(taus, bins=bins) indmax = h[0].argmax() val = 0.5 * (h[1][:-1][indmax] + h[1][1:][indmax])[0]