Extract the rotation matrix from dimRedResult objects derived from PCA and FastICA
Arguments
- x
of type
dimRedResult
Details
The data has to be pre-processed the same way as the method does, e.g. centering and/or scaling.
Examples
dat <- loadDataSet("Iris")
pca <- embed(dat, "PCA")
rot_pca <- getRotationMatrix(pca)
scale(getData(dat), TRUE, FALSE) %*% rot_pca - getData(getDimRedData(pca))
#> PC1 PC2
#> [1,] 0 0
#> [2,] 0 0
#> [3,] 0 0
#> [4,] 0 0
#> [5,] 0 0
#> [6,] 0 0
#> [7,] 0 0
#> [8,] 0 0
#> [9,] 0 0
#> [10,] 0 0
#> [11,] 0 0
#> [12,] 0 0
#> [13,] 0 0
#> [14,] 0 0
#> [15,] 0 0
#> [16,] 0 0
#> [17,] 0 0
#> [18,] 0 0
#> [19,] 0 0
#> [20,] 0 0
#> [21,] 0 0
#> [22,] 0 0
#> [23,] 0 0
#> [24,] 0 0
#> [25,] 0 0
#> [26,] 0 0
#> [27,] 0 0
#> [28,] 0 0
#> [29,] 0 0
#> [30,] 0 0
#> [31,] 0 0
#> [32,] 0 0
#> [33,] 0 0
#> [34,] 0 0
#> [35,] 0 0
#> [36,] 0 0
#> [37,] 0 0
#> [38,] 0 0
#> [39,] 0 0
#> [40,] 0 0
#> [41,] 0 0
#> [42,] 0 0
#> [43,] 0 0
#> [44,] 0 0
#> [45,] 0 0
#> [46,] 0 0
#> [47,] 0 0
#> [48,] 0 0
#> [49,] 0 0
#> [50,] 0 0
#> [51,] 0 0
#> [52,] 0 0
#> [53,] 0 0
#> [54,] 0 0
#> [55,] 0 0
#> [56,] 0 0
#> [57,] 0 0
#> [58,] 0 0
#> [59,] 0 0
#> [60,] 0 0
#> [61,] 0 0
#> [62,] 0 0
#> [63,] 0 0
#> [64,] 0 0
#> [65,] 0 0
#> [66,] 0 0
#> [67,] 0 0
#> [68,] 0 0
#> [69,] 0 0
#> [70,] 0 0
#> [71,] 0 0
#> [72,] 0 0
#> [73,] 0 0
#> [74,] 0 0
#> [75,] 0 0
#> [76,] 0 0
#> [77,] 0 0
#> [78,] 0 0
#> [79,] 0 0
#> [80,] 0 0
#> [81,] 0 0
#> [82,] 0 0
#> [83,] 0 0
#> [84,] 0 0
#> [85,] 0 0
#> [86,] 0 0
#> [87,] 0 0
#> [88,] 0 0
#> [89,] 0 0
#> [90,] 0 0
#> [91,] 0 0
#> [92,] 0 0
#> [93,] 0 0
#> [94,] 0 0
#> [95,] 0 0
#> [96,] 0 0
#> [97,] 0 0
#> [98,] 0 0
#> [99,] 0 0
#> [100,] 0 0
#> [101,] 0 0
#> [102,] 0 0
#> [103,] 0 0
#> [104,] 0 0
#> [105,] 0 0
#> [106,] 0 0
#> [107,] 0 0
#> [108,] 0 0
#> [109,] 0 0
#> [110,] 0 0
#> [111,] 0 0
#> [112,] 0 0
#> [113,] 0 0
#> [114,] 0 0
#> [115,] 0 0
#> [116,] 0 0
#> [117,] 0 0
#> [118,] 0 0
#> [119,] 0 0
#> [120,] 0 0
#> [121,] 0 0
#> [122,] 0 0
#> [123,] 0 0
#> [124,] 0 0
#> [125,] 0 0
#> [126,] 0 0
#> [127,] 0 0
#> [128,] 0 0
#> [129,] 0 0
#> [130,] 0 0
#> [131,] 0 0
#> [132,] 0 0
#> [133,] 0 0
#> [134,] 0 0
#> [135,] 0 0
#> [136,] 0 0
#> [137,] 0 0
#> [138,] 0 0
#> [139,] 0 0
#> [140,] 0 0
#> [141,] 0 0
#> [142,] 0 0
#> [143,] 0 0
#> [144,] 0 0
#> [145,] 0 0
#> [146,] 0 0
#> [147,] 0 0
#> [148,] 0 0
#> [149,] 0 0
#> [150,] 0 0
if(requireNamespace("fastICA", quietly = TRUE)) {
ica <- embed(dat, "FastICA")
rot_ica <- getRotationMatrix(ica)
scale(getData(dat), TRUE, FALSE) %*% rot_ica - getData(getDimRedData(ica))
}
#> ICA1 ICA2
#> [1,] -6.661338e-16 1.110223e-16
#> [2,] -8.881784e-16 1.110223e-16
#> [3,] -8.881784e-16 0.000000e+00
#> [4,] -8.881784e-16 1.110223e-16
#> [5,] -6.661338e-16 1.110223e-16
#> [6,] -4.440892e-16 4.440892e-16
#> [7,] -6.661338e-16 1.665335e-16
#> [8,] -6.661338e-16 1.387779e-16
#> [9,] -6.661338e-16 6.661338e-16
#> [10,] -4.440892e-16 5.551115e-17
#> [11,] -8.881784e-16 2.220446e-16
#> [12,] -6.661338e-16 1.249001e-16
#> [13,] -4.440892e-16 1.110223e-16
#> [14,] -4.440892e-16 4.440892e-16
#> [15,] -4.440892e-16 0.000000e+00
#> [16,] -4.440892e-16 0.000000e+00
#> [17,] -6.661338e-16 2.220446e-16
#> [18,] -4.440892e-16 1.665335e-16
#> [19,] -2.220446e-16 4.440892e-16
#> [20,] -4.440892e-16 1.110223e-16
#> [21,] -4.440892e-16 4.440892e-16
#> [22,] -4.440892e-16 1.110223e-16
#> [23,] -8.881784e-16 1.249001e-16
#> [24,] -4.440892e-16 2.636780e-16
#> [25,] -6.661338e-16 1.387779e-16
#> [26,] -4.440892e-16 2.220446e-16
#> [27,] -4.440892e-16 1.110223e-16
#> [28,] -6.661338e-16 1.110223e-16
#> [29,] -6.661338e-16 1.665335e-16
#> [30,] -6.661338e-16 2.220446e-16
#> [31,] -6.661338e-16 2.220446e-16
#> [32,] -4.440892e-16 2.220446e-16
#> [33,] -4.440892e-16 0.000000e+00
#> [34,] -6.661338e-16 0.000000e+00
#> [35,] -6.661338e-16 1.110223e-16
#> [36,] -6.661338e-16 2.550044e-16
#> [37,] -8.881784e-16 0.000000e+00
#> [38,] -4.440892e-16 1.110223e-16
#> [39,] -8.881784e-16 2.220446e-16
#> [40,] -6.661338e-16 1.665335e-16
#> [41,] -6.661338e-16 0.000000e+00
#> [42,] -6.661338e-16 0.000000e+00
#> [43,] -8.881784e-16 0.000000e+00
#> [44,] -4.440892e-16 1.110223e-16
#> [45,] -6.661338e-16 1.110223e-16
#> [46,] -4.440892e-16 1.110223e-16
#> [47,] -6.661338e-16 2.220446e-16
#> [48,] -6.661338e-16 1.110223e-16
#> [49,] -6.661338e-16 2.220446e-16
#> [50,] -8.881784e-16 1.249001e-16
#> [51,] -6.106227e-16 2.220446e-16
#> [52,] -6.106227e-16 2.220446e-16
#> [53,] -6.661338e-16 0.000000e+00
#> [54,] -5.551115e-16 0.000000e+00
#> [55,] -5.551115e-16 1.110223e-16
#> [56,] -5.551115e-16 2.220446e-16
#> [57,] -6.106227e-16 2.220446e-16
#> [58,] -6.106227e-16 0.000000e+00
#> [59,] -5.551115e-16 1.110223e-16
#> [60,] -6.106227e-16 2.220446e-16
#> [61,] -5.967449e-16 4.440892e-16
#> [62,] -6.106227e-16 2.220446e-16
#> [63,] -6.106227e-16 2.220446e-16
#> [64,] -5.551115e-16 2.775558e-16
#> [65,] -5.863365e-16 3.330669e-16
#> [66,] -6.106227e-16 0.000000e+00
#> [67,] -6.106227e-16 2.220446e-16
#> [68,] -5.828671e-16 2.220446e-16
#> [69,] -5.551115e-16 0.000000e+00
#> [70,] -5.828671e-16 6.661338e-16
#> [71,] -5.551115e-16 2.359224e-16
#> [72,] -6.383782e-16 2.498002e-16
#> [73,] -5.551115e-16 2.220446e-16
#> [74,] -5.551115e-16 2.775558e-16
#> [75,] -5.551115e-16 2.220446e-16
#> [76,] -6.106227e-16 1.110223e-16
#> [77,] -5.551115e-16 0.000000e+00
#> [78,] -4.440892e-16 1.110223e-16
#> [79,] -6.106227e-16 2.220446e-16
#> [80,] -5.828671e-16 0.000000e+00
#> [81,] -5.967449e-16 2.220446e-16
#> [82,] -5.967449e-16 2.220446e-16
#> [83,] -6.383782e-16 2.220446e-16
#> [84,] -6.661338e-16 2.220446e-16
#> [85,] -6.106227e-16 2.220446e-16
#> [86,] -6.106227e-16 2.220446e-16
#> [87,] -6.661338e-16 1.110223e-16
#> [88,] -5.551115e-16 3.330669e-16
#> [89,] -6.106227e-16 3.330669e-16
#> [90,] -6.106227e-16 0.000000e+00
#> [91,] -5.551115e-16 0.000000e+00
#> [92,] -6.106227e-16 2.393918e-16
#> [93,] -6.106227e-16 2.220446e-16
#> [94,] -6.245005e-16 0.000000e+00
#> [95,] -5.551115e-16 2.220446e-16
#> [96,] -6.106227e-16 2.775558e-16
#> [97,] -6.106227e-16 2.220446e-16
#> [98,] -5.551115e-16 2.359224e-16
#> [99,] -6.106227e-16 2.220446e-16
#> [100,] -6.106227e-16 3.330669e-16
#> [101,] -6.661338e-16 2.081668e-16
#> [102,] -6.661338e-16 0.000000e+00
#> [103,] -6.661338e-16 1.110223e-16
#> [104,] -6.661338e-16 2.775558e-16
#> [105,] -4.440892e-16 1.179612e-16
#> [106,] -6.661338e-16 4.440892e-16
#> [107,] -5.551115e-16 0.000000e+00
#> [108,] -4.440892e-16 2.220446e-16
#> [109,] -4.440892e-16 1.665335e-16
#> [110,] -6.661338e-16 0.000000e+00
#> [111,] -6.661338e-16 0.000000e+00
#> [112,] -5.551115e-16 2.220446e-16
#> [113,] -6.661338e-16 1.110223e-16
#> [114,] -5.551115e-16 0.000000e+00
#> [115,] -6.661338e-16 2.220446e-16
#> [116,] -5.551115e-16 2.220446e-16
#> [117,] -6.661338e-16 1.387779e-16
#> [118,] -4.440892e-16 8.881784e-16
#> [119,] -6.661338e-16 3.330669e-16
#> [120,] -5.551115e-16 0.000000e+00
#> [121,] -4.440892e-16 1.110223e-16
#> [122,] -6.661338e-16 2.220446e-16
#> [123,] -6.661338e-16 4.440892e-16
#> [124,] -6.661338e-16 2.775558e-16
#> [125,] -6.661338e-16 2.220446e-16
#> [126,] -4.440892e-16 4.440892e-16
#> [127,] -5.551115e-16 2.220446e-16
#> [128,] -6.661338e-16 2.498002e-16
#> [129,] -4.440892e-16 2.220446e-16
#> [130,] -6.661338e-16 2.220446e-16
#> [131,] -4.440892e-16 3.330669e-16
#> [132,] -8.881784e-16 8.881784e-16
#> [133,] -6.661338e-16 2.220446e-16
#> [134,] -6.661338e-16 2.220446e-16
#> [135,] -5.551115e-16 2.220446e-16
#> [136,] -6.661338e-16 4.440892e-16
#> [137,] -4.440892e-16 4.440892e-16
#> [138,] -5.551115e-16 2.498002e-16
#> [139,] -5.551115e-16 2.220446e-16
#> [140,] -6.661338e-16 1.110223e-16
#> [141,] -6.661338e-16 2.220446e-16
#> [142,] -5.551115e-16 1.110223e-16
#> [143,] -6.661338e-16 0.000000e+00
#> [144,] -4.440892e-16 1.110223e-16
#> [145,] -4.440892e-16 2.220446e-16
#> [146,] -6.661338e-16 1.110223e-16
#> [147,] -5.551115e-16 3.330669e-16
#> [148,] -6.661338e-16 1.665335e-16
#> [149,] -6.661338e-16 3.330669e-16
#> [150,] -6.661338e-16 2.220446e-16