test_statis.py 文件源码

def test_RvPCA():
    data = get_wine_data()
    for d in data:
        d.cross_product()

    n_datasets = len(data)

    expected_output = np.array([[0.10360263],
                                [0.10363524],
                                [0.09208477],
                                [0.10370834],
                                [0.08063234],
                                [0.09907428],
                                [0.09353886],
                                [0.08881811],
                                [0.1110871],
                                [0.12381833]])
    output, _, _ = rv_pca(data, n_datasets)

    np.testing.assert_array_almost_equal(np.atleast_2d(output).T, expected_output)