docs: def empirical_rate(data, sigma=3.0):
          """
          Smooth count data to get an empirical rate
          """
          from scipy.ndimage.filters import gaussian_filter1d
          return 0.001 + gaussian_filter1d(data.astype(np.float), sigma, axis=0)