docs: def resample_1d(arr, n_out=None, random_state=None):
      """Resample an array, with replacement.
  
      Parameters
      ==========
      arr: np.ndarray
          The array is resampled along the first axis.
      n_out: int, optional
          Number of samples to return. If not specified,
          return ``len(arr)`` samples.
      """
      if random_state is None:
          random_state = np.random.RandomState()
      arr = np.atleast_1d(arr)
      n = len(arr)
      if n_out is None:
          n_out = n
      idx = random_state.randint(0, n, size=n)
      return arr[idx]