def METHOD_NAME(self):
    x = np.array(range(1, 21), dtype=float)
    y = x + np.log(x)
    stats = ddm.distance_statistics(x, y)
    # Values were obtained using the R `energy` package
    assert_almost_equal(stats.test_statistic, 398.94623, 5)
    assert_almost_equal(stats.distance_correlation, 0.9996107, 4)
    assert_almost_equal(stats.distance_covariance, 4.4662414, 4)
    assert_almost_equal(stats.dvar_x, 4.2294799, 4)
    assert_almost_equal(stats.dvar_y, 4.7199304, 4)
    assert_almost_equal(stats.S, 49.880200, 4)