sdata = data[data[:, (0)].argsort()]
unq_x, breaks = np.unique(sdata[:, (0)], return_index=True)
breaks = np.append(breaks, data.shape[0])
idx_range = np.arange(data.shape[0])
avg_y, _, _ = bstat(x=idx_range, values=sdata[:, (1)], statistic='mean',
    bins=breaks)
std_y, _, _ = bstat(x=idx_range, values=sdata[:, (1)], statistic='std',
    bins=breaks)
Выполнение группового среднего и стандартного отклонения с массивами NumPy