scikit-hep/hepstats

Bayesian blocking gamma prior doesn't make sense

Moelf opened this issue · 0 comments

Moelf commented

when gamma is supplied,

return -np.log(self.gamma)

A_R = fit_vec - prior.calc(R + 1)

this is a constant prior, this doesn't make sense because the prior should depend on N_b (https://arxiv.org/pdf/1708.00810.pdf):
image