docs: def bin_width_doane(a):
      """Compute accuracy (bin width) based on Doane's formula
  
      References
      ----------
      - `<https://en.wikipedia.org/wiki/Histogram#Number_of_bins_and_width>`_
      - `<https://stats.stackexchange.com/questions/55134/
        doanes-formula-for-histogram-binning>`_
      """
      bad = np.isnan(a) | np.isinf(a)
      data = a[~bad]
      n = data.size
      g1 = skew(data)
      sigma_g1 = np.sqrt(6 * (n - 2) / ((n + 1) * (n + 3)))
      k = 1 + np.log2(n) + np.log2(1 + np.abs(g1) / sigma_g1)
      acc = (data.max() - data.min()) / k
      return acc