docs: def cvm_stat(data):
      """
      Calculates the Cramer-von Mises statistic for sorted values from U(0, 1).
      """
      samples2 = 2 * len(data)
      minuends = arange(1, samples2, 2) / samples2
      return 1 / (6 * samples2) + ((minuends - data) ** 2).sum()