def scott_binning(data, range=None, **kwargs) -> StaticBinning: """Binning schema based on Scott's rule (from astropy). Parameters ---------- data: arraylike range: Optional[tuple] See also -------- astropy.stats.histogram.scott_bin_width astropy.stats.histogram.histogram """ from astropy.stats.histogram import scott_bin_width if range is not None: data = data[(data >= range[0]) & (data <= range[1])] _, edges = scott_bin_width(data, True) return NumpyBinning(edges, **kwargs)