view release on metacpan or  search on metacpan

perl/t/16_pca.t  view on Meta::CPAN

is (sprintf("%7.4f", $coordinates->[5][1]), "-0.2841");
is (sprintf("%7.4f", $coordinates->[6][0]), "-3.2019");
is (sprintf("%7.4f", $coordinates->[6][1]), " 0.0197");
is (sprintf("%7.4f", $coordinates->[7][0]), "-3.2019");
is (sprintf("%7.4f", $coordinates->[7][1]), " 0.0197");
is (sprintf("%7.4f", $coordinates->[8][0]), " 0.4698");
is (sprintf("%7.4f", $coordinates->[8][1]), "-0.1778");
is (sprintf("%7.4f", $coordinates->[9][0]), "-2.5550");
is (sprintf("%7.4f", $coordinates->[9][1]), " 0.1973");
is (sprintf("%7.4f", $coordinates->[10][0]), "-2.5034");
is (sprintf("%7.4f", $coordinates->[10][1]), "-0.1595");
is (sprintf("%7.4f", $coordinates->[11][0]), "-2.4366");
is (sprintf("%7.4f", $coordinates->[11][1]), "-0.2339");
is (sprintf("%7.4f", $coordinates->[12][0]), "-2.2802");
is (sprintf("%7.4f", $coordinates->[12][1]), " 0.0409");
is (sprintf("%7.4f", $pc->[0][0]), "-0.6681");
is (sprintf("%7.4f", $pc->[0][1]), "-0.7441");
is (sprintf("%7.4f", $pc->[1][0]), " 0.7441");
is (sprintf("%7.4f", $pc->[1][1]), "-0.6681");
is (sprintf("%7.4f", $eigenvalues->[0]), " 9.3110");
is (sprintf("%7.4f", $eigenvalues->[1]), " 1.4437");


($mean, $coordinates, $pc, $eigenvalues) =  Algorithm::Cluster::pca($data2);

is (sprintf("%7.4f", $mean->[0]), " 2.7500");
is (sprintf("%7.4f", $mean->[1]), " 5.8500");
is (sprintf("%7.4f", $mean->[2]), " 3.9500");
is (sprintf("%7.4f", $mean->[3]), " 6.0500");
is (sprintf("%7.4f", $mean->[4]), " 6.2750");
is (sprintf("%7.4f", $mean->[5]), " 8.0750");
is (sprintf("%7.4f", $coordinates->[0][0]), " 2.6461");
is (sprintf("%7.4f", $coordinates->[0][1]), "-2.1422");
is (sprintf("%7.4f", $coordinates->[0][2]), "-0.5662");
is (sprintf("%7.4f", abs($coordinates->[0][3])), " 0.0000");
is (sprintf("%7.4f", $coordinates->[1][0]), " 2.0644");
is (sprintf("%7.4f", $coordinates->[1][1]), " 0.5554");
is (sprintf("%7.4f", $coordinates->[1][2]), " 1.4819");
is (sprintf("%7.4f", abs($coordinates->[1][3])), " 0.0000");
is (sprintf("%7.4f", $coordinates->[2][0]), " 1.0687");
is (sprintf("%7.4f", $coordinates->[2][1]), " 1.9994");
is (sprintf("%7.4f", $coordinates->[2][2]), "-1.0007");
is (sprintf("%7.4f", abs($coordinates->[2][3])), " 0.0000");
is (sprintf("%7.4f", $coordinates->[3][0]), "-5.7792");
is (sprintf("%7.4f", $coordinates->[3][1]), "-0.4127");
is (sprintf("%7.4f", $coordinates->[3][2]), " 0.0851");
is (sprintf("%7.4f", abs($coordinates->[3][3])), " 0.0000");
is (sprintf("%7.4f", $pc->[0][0]), "-0.2638");
is (sprintf("%7.4f", $pc->[0][1]), " 0.0648");
is (sprintf("%7.4f", $pc->[0][2]), "-0.9176");
is (sprintf("%7.4f", $pc->[0][3]), " 0.2615");
is (sprintf("%7.4f", $pc->[1][0]), " 0.0507");
is (sprintf("%7.4f", $pc->[1][1]), " 0.6862");
is (sprintf("%7.4f", $pc->[1][2]), " 0.1382");
is (sprintf("%7.4f", $pc->[1][3]), " 0.1978");
is (sprintf("%7.4f", $pc->[2][0]), "-0.6300");
is (sprintf("%7.4f", $pc->[2][1]), " 0.0912");
is (sprintf("%7.4f", $pc->[2][2]), " 0.0456");
is (sprintf("%7.4f", $pc->[2][3]), "-0.6746");
# The last eigenvalue is zero. The corresponding eigenvector is strongly
# affected by roundoff error, so we don't test it. For PCA, it doesn't matter
# since the coordinates along this eigenvector are zero anyway.
is (sprintf("%7.4f", $eigenvalues->[0]), " 6.7679");
is (sprintf("%7.4f", $eigenvalues->[1]), " 3.0109");
is (sprintf("%7.4f", $eigenvalues->[2]), " 1.8776");
is (sprintf("%7.4f", abs($eigenvalues->[3])), " 0.0000");