ecdf = np.arange(1, n + 1) / (n+1)