stat1, edges1, bc = binned_statistic(x, u, 'std', bins=10) stat2, edges2, bc = binned_statistic(x, u, np.std, bins=10) assert_allclose(stat1, stat2) def test_empty ...