docs: def bootstrap_params(rv_cont, data, n_iter=5, **kwargs):
      """Bootstrap the fit params of a distribution.
  
      Parameters
      ==========
      rv_cont: scipy.stats.rv_continuous instance
          The distribution which to fit.
      data: array-like, 1d
          The data on which to fit.
      n_iter: int [default=10]
          Number of bootstrap iterations.
      """
      fit_res = []
      for _ in range(n_iter):
          params = rv_cont.fit(resample_1d(data, **kwargs))
          fit_res.append(params)
      fit_res = np.array(fit_res)
      return fit_res