code: def detect_mean_shift(self, ts, B=1000):
          
          x = np.arange(0, len(ts))
          stat_ts_func = self.compute_balance_mean_ts
          null_ts_func = self.shuffle_timeseries
          stats_ts, pvals, nums = self.get_ts_stats_significance(x, ts, stat_ts_func, null_ts_func, B=B, permute_fast=True)
          return stats_ts, pvals, nums
docs: Detect mean shift in a time series. B is number of bootstrapped
              samples to draw.