cppEDM/tests/MultiviewTest.cc at master · SugiharaLab/cppEDM · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
// Multiview Test

#include "TestCommon.h"

int main( int argc, char *argv[] ) {

    //---------------------------------------------------------
    // Load valid cppEDM output
    //---------------------------------------------------------
    DataFrame < double > validCppPredict( "./data/",
                                          "Multiview_pred_valid.csv" );

    DataFrame < double > validCppCombos( "./data/",
                                         "Multiview_combos_valid.csv",
                                         true ); // noTime = true

    // Generate new cppEDM output
    MultiviewValues MV = Multiview( "../data/",
                                    "block_3sp.csv",
                                    "./",
                                    "MultiviewBlock3sp.csv",
                                    "1 100",
                                    "101 198",
                                    0,             // D
                                    3,             // E
                                    1,             // Tp
                                    0,             // knn
                                    -1,            // tau
                                    "x_t y_t z_t", // columns
                                    "x_t",         // target,
                                    0,             // multiview
                                    0,             // exclusionRadius
                                    true,          // trainLib
                                    false,         // excludeTarget
                                    false,         // parameterList
                                    false,         // verbose,
                                    1 );           // nThreads

    DataFrame< double > combos = MV.ComboRho;
    DataFrame< double > output = MV.Predictions;

    combos.WriteData( "./", "Multiview_combos.csv" );
    output.WriteData( "./", "Multiview_pred.csv"   );

    //VectorError ve = ComputeError( output.VectorColumnName("Observations"),
    //                               output.VectorColumnName("Predictions"));

    //std::cout << "Multiview(): rho " << ve.rho
    //          << "  MAE " << ve.MAE << "  RMSE " << ve.RMSE << std::endl;

    // comparison
    MakeTest ("Multiview: combos test",     validCppCombos,  combos );
    MakeTest ("Multiview: prediction test", validCppPredict, output );

    return 0;
}