These are chat archives for astropy/astropy
astropy.stats.bayesian_blocksfunction (which would give me the binlimits I need) has an argument
sigmawhich should do what I need, I suppose. but I'm confused that I have to pass two arrays (t & x) as arguments (apart from the sigma/data errors of course), not just the one I want to bin (x). And in the documentation of the
stats.histrogramfunction I don't see any support for data uncertainties, not even through the (unimplemented) weights.