docs: def percentile_index(a, q):
      """
      Returns the index of the value at the Qth percentile in array a.
      """
      return np.where(
          a==np.percentile(a, q, interpolation='nearest')
      )[0][0]