tests/testthat/_snaps/pca.md

can prepare pca

Code
  prepare_pca(df = test_count_matrix[3:6])
Output
  $data
  Standard deviations (1, .., p=4):
  [1] 218.22860  81.12065  33.51049  25.88166

  Rotation (n x k) = (4 x 4):
              PC1        PC2          PC3          PC4
  T0_- -0.8423279 -0.5387762  0.013514466 -0.004613692
  T3_A -0.3503926  0.5671455  0.641227641 -0.379997505
  T3_B -0.2739171  0.4206085 -0.009173189  0.864854762
  T3_C -0.3044309  0.4595157 -0.767176854 -0.328034833

  $variance_explained
     PC variance_explained
  1 PC1          85.046773
  2 PC2          11.751608
  3 PC3           2.005378
  4 PC4           1.196241

can prepare pca with scaling

Code
  prepare_pca(df = test_count_matrix[3:6], scaled = T)
Output
  $data
  Standard deviations (1, .., p=4):
  [1] 1.8290114 0.6218325 0.3945069 0.3352696

  Rotation (n x k) = (4 x 4):
              PC1        PC2         PC3          PC4
  T0_- -0.4553397 -0.8900893 -0.01886846 -0.007123416
  T3_A -0.5140102  0.2700273 -0.55776245  0.593110301
  T3_B -0.5171755  0.2754655 -0.21950875 -0.780041125
  T3_C -0.5108706  0.2427860  0.80022556  0.199261814

  $variance_explained
     PC variance_explained
  1 PC1          83.632072
  2 PC2           9.666893
  3 PC3           3.890893
  4 PC4           2.810142

can prepare pca with transform

Code
  prepare_pca(df = test_count_matrix[3:6], transform = T)
Output
  $data
  Standard deviations (1, .., p=4):
  [1] 2.732481e+03 4.775840e+02 3.689677e+02 3.439720e-13

  Rotation (n x k) = (498 x 4):
                   PC1           PC2           PC3           PC4
    [1,] -6.312612e-02  1.828884e-02  4.732036e-02  8.825290e-01
    [2,] -1.540520e-02 -7.086642e-02  4.699002e-02  3.068933e-01
    [3,] -3.957768e-02 -1.183371e-01  4.188281e-02  2.233352e-01
    [4,] -3.214614e-02 -4.537490e-02  7.243763e-02 -6.400365e-03
    [5,] -6.938214e-03 -6.560316e-02 -1.288658e-02  9.446512e-03
    [6,] -1.777424e-02  5.330689e-02  3.921640e-02 -1.110433e-02
    [7,] -8.951885e-02 -1.950371e-02 -8.877509e-02  9.225551e-03
    [8,] -4.482429e-02  7.258838e-02  4.431838e-02  1.093555e-02
    [9,] -3.309397e-02  2.029137e-02  1.392982e-02 -1.031638e-02
   [10,] -8.556100e-02  3.284823e-02 -1.565225e-02  3.534297e-03
   [11,] -1.748502e-02 -7.665229e-02 -8.047520e-04 -9.135440e-04
   [12,] -2.080069e-04 -7.202347e-02 -4.836932e-02  1.850834e-03
   [13,] -3.420846e-02 -4.093458e-03  4.581830e-02 -1.593698e-02
   [14,]  2.624416e-03 -3.390809e-02 -3.994736e-02  4.259657e-03
   [15,] -3.112586e-02 -7.814918e-03  2.415454e-02 -9.658904e-03
   [16,] -2.309372e-02 -3.034630e-02 -1.857048e-02  1.085345e-02
   [17,] -5.865665e-02 -5.299706e-02  2.744350e-02 -1.459890e-02
   [18,] -3.543211e-02 -7.708076e-02 -5.326167e-02 -1.347050e-03
   [19,] -4.038980e-02  6.248892e-02 -1.141207e-02 -1.193103e-02
   [20,] -4.891850e-02 -6.726534e-02  9.025462e-02 -2.297290e-02
   [21,]  5.664776e-05 -1.260434e-01  1.028370e-01 -2.422920e-02
   [22,] -3.089657e-02  1.676419e-02 -5.640369e-02  5.716813e-03
   [23,] -6.833143e-03  3.733575e-03  1.206621e-03 -2.488551e-03
   [24,] -4.784317e-02 -5.853481e-02  5.602155e-03  2.021375e-03
   [25,] -1.630722e-03 -7.585054e-02  3.283900e-02 -8.197578e-03
   [26,] -1.597661e-02 -1.618607e-02  6.710222e-06  1.341768e-03
   [27,] -2.560882e-02  1.743351e-02  2.400182e-02 -1.109032e-02
   [28,] -3.734785e-02 -5.266493e-02 -1.997084e-02 -2.313047e-02
   [29,] -7.059994e-03 -6.624292e-02  1.754160e-02 -1.186180e-03
   [30,] -1.104333e-01  1.068366e-01  7.377178e-02 -2.483492e-03
   [31,] -2.667432e-02  2.975190e-02 -2.176578e-02 -1.087562e-02
   [32,] -7.055934e-02 -1.391863e-02 -6.516461e-02 -1.846219e-02
   [33,] -1.643311e-02 -1.256770e-02 -1.443025e-02  4.987450e-03
   [34,] -1.494780e-02 -8.794791e-03  1.336786e-02 -1.012285e-03
   [35,] -5.693646e-02  5.807056e-02  1.839586e-02 -4.223290e-03
   [36,] -3.541844e-02 -2.393489e-02  4.072977e-02 -3.876576e-03
   [37,] -1.056079e-02 -2.670732e-02  4.021422e-02 -1.677647e-03
   [38,] -5.741585e-03 -2.366006e-02 -4.769311e-03 -1.127072e-03
   [39,] -2.590281e-02  4.046387e-02  6.492245e-02  3.254225e-03
   [40,] -2.075440e-02 -2.044729e-02 -7.897870e-04  6.520001e-03
   [41,] -5.214033e-02  1.209627e-02 -6.894807e-03 -1.462630e-02
   [42,] -2.749920e-02  7.114683e-03  9.845494e-03 -8.986598e-03
   [43,] -4.864495e-02  8.965257e-03 -2.516807e-02 -1.678366e-02
   [44,] -3.781367e-02  1.153702e-02  2.200027e-02 -1.038102e-02
   [45,] -4.890342e-02  2.666947e-02  5.695231e-03  1.741934e-02
   [46,] -6.040728e-02  5.688180e-02  1.072743e-01  3.131736e-03
   [47,] -3.047595e-02 -6.591744e-02  4.832229e-02  1.040586e-02
   [48,] -5.635674e-02 -7.901073e-02  8.121881e-03  5.289951e-03
   [49,] -9.763905e-02 -2.631455e-02 -4.303934e-02  3.776157e-02
   [50,] -3.179084e-02 -4.217708e-02 -8.541190e-02  2.706215e-03
   [51,] -1.543858e-02  1.109722e-03  3.686443e-02 -4.562293e-03
   [52,] -3.113861e-02  1.093086e-02  4.476826e-02 -7.199986e-03
   [53,] -5.696627e-03 -2.690888e-02  3.415335e-02  6.721439e-04
   [54,] -2.864619e-02 -2.725329e-02 -3.024180e-02 -2.252121e-03
   [55,] -6.502069e-03 -2.150499e-02 -1.674955e-02  4.296213e-03
   [56,] -1.641802e-02 -2.663657e-02  1.448872e-02  7.544101e-04
   [57,] -7.237063e-02  4.207408e-02  7.323077e-03 -2.029387e-02
   [58,] -2.297522e-02 -1.211669e-03  1.330126e-02 -5.757474e-03
   [59,] -4.792596e-02  2.612459e-02  1.178491e-02 -1.918985e-03
   [60,] -2.044904e-02  1.266306e-02 -4.545472e-03 -6.060721e-03
   [61,] -2.363222e-02 -1.359915e-02 -1.021883e-01  4.245837e-03
   [62,] -4.502862e-03 -4.387247e-03 -4.039799e-03  2.900829e-04
   [63,] -1.491981e-02 -3.250619e-03 -1.479585e-02  7.921967e-04
   [64,] -2.685595e-03 -1.078122e-01 -5.604804e-02  1.067741e-02
   [65,] -6.262641e-02  3.683284e-02  3.333741e-02  2.101961e-03
   [66,]  1.067225e-03  2.915837e-02  4.232044e-02 -2.839610e-03
   [67,] -4.844707e-02  2.580802e-02  1.391249e-02 -6.990143e-03
   [68,] -7.442851e-03 -7.225809e-03 -2.582764e-02  1.669868e-03
   [69,] -3.841762e-02 -2.599724e-02 -1.158152e-02  1.867376e-02
   [70,] -1.452113e-02 -4.900542e-02  4.590220e-02 -1.211060e-03
   [71,] -1.393558e-02  2.489095e-02 -4.979933e-02  2.863988e-03
   [72,] -1.222015e-02 -4.414236e-03 -2.255263e-02 -1.852338e-03
   [73,] -1.586680e-02  4.758438e-03  6.703153e-03 -4.447361e-04
   [74,] -1.386712e-02  6.092282e-02  1.580476e-01 -9.449121e-03
   [75,] -3.479929e-02 -9.428059e-03 -5.321405e-02  3.588147e-03
   [76,] -9.781339e-02  7.835485e-02  1.427098e-03  1.093746e-02
   [77,]  8.546662e-03 -4.598640e-02 -4.833527e-03 -4.959006e-04
   [78,] -2.966412e-02  1.442336e-02 -1.050950e-02 -4.341713e-03
   [79,] -3.118120e-02  2.652836e-02  4.808451e-02 -4.154429e-03
   [80,] -1.459305e-02  1.163231e-02 -6.696172e-03 -1.775277e-03
   [81,] -1.426710e-02  1.728451e-02 -7.813605e-03 -4.282879e-03
   [82,] -9.683183e-03  2.069149e-03  4.997030e-02 -1.233796e-03
   [83,] -1.489966e-02 -6.355378e-02  1.187196e-02 -9.106478e-04
   [84,] -4.880931e-02 -9.520976e-03  4.397121e-02 -8.838806e-03
   [85,] -1.874924e-02 -5.843354e-02  3.969319e-02 -1.008458e-02
   [86,] -8.898412e-03 -3.086572e-02 -1.444690e-02  4.370363e-03
   [87,] -2.041227e-02 -2.514705e-02  4.866825e-02 -7.203582e-03
   [88,] -5.566186e-02  8.757350e-03 -8.495471e-02  1.702433e-02
   [89,] -5.637458e-02 -2.170242e-02  3.828055e-02 -2.285379e-03
   [90,] -1.971857e-02  2.444846e-03 -8.788620e-03 -9.582356e-03
   [91,] -2.194073e-02 -2.294826e-03 -6.923952e-04  3.388832e-03
   [92,] -1.214884e-02 -6.169609e-02  1.246075e-02 -3.751012e-03
   [93,] -2.971130e-02  4.221220e-02  4.181805e-02  4.818582e-03
   [94,] -3.035320e-02 -1.768813e-02 -9.357736e-03  7.485960e-03
   [95,] -2.192053e-02 -1.381048e-02 -2.400702e-02  2.140439e-03
   [96,] -1.193357e-02 -2.700296e-02 -6.806616e-02  3.746697e-03
   [97,] -2.241590e-02  8.285377e-03  4.850086e-02  1.723588e-03
   [98,] -3.136043e-02  3.469761e-02  7.993695e-02 -9.034507e-03
   [99,] -1.947151e-03 -2.296614e-02  2.566020e-02 -1.907733e-03
  [100,] -1.111037e-01  4.966286e-02  3.210062e-02 -1.881561e-02
  [101,] -2.569543e-03  1.637901e-04  4.652086e-03 -1.526052e-03
  [102,] -2.323466e-02 -1.041866e-01 -2.045816e-02  6.218260e-04
  [103,] -2.307241e-02 -7.158902e-02  1.516625e-02  2.128100e-03
  [104,] -6.080119e-02  6.191613e-02  3.188570e-02 -1.136081e-02
  [105,] -3.023765e-02 -6.443423e-02  3.347982e-02  5.686602e-03
  [106,] -9.372036e-02  1.038309e-02 -2.268048e-02 -2.682249e-02
  [107,] -2.676843e-02 -4.206133e-02  5.343659e-02 -1.723435e-04
  [108,] -7.244495e-02  8.122125e-02  8.973766e-02 -7.497122e-03
  [109,] -1.305168e-03 -5.425584e-02  5.768695e-03  3.351603e-04
  [110,]  5.587640e-03 -1.824887e-02  3.762176e-03  2.237617e-03
  [111,]  1.154212e-03 -1.210426e-02  8.844292e-02 -6.014422e-03
  [112,] -1.851503e-02 -1.107758e-01 -1.394103e-02  6.828945e-03
  [113,] -5.997968e-02  1.694335e-02 -7.418846e-03 -8.179328e-03
  [114,] -5.293605e-02  1.406366e-02  3.621677e-03 -1.455444e-02
  [115,] -8.571378e-03  1.194735e-02 -2.930269e-02 -8.991941e-04
  [116,] -5.696341e-03 -3.242270e-02  4.659274e-02 -3.269346e-03
  [117,] -1.672725e-02 -6.134014e-02  9.835419e-02 -1.957479e-03
  [118,] -8.124886e-02  3.313667e-02 -6.583330e-02  3.922675e-02
  [119,] -3.558151e-02 -6.576321e-02 -4.111754e-03 -1.128614e-02
  [120,] -5.216732e-02 -1.192798e-01  8.327059e-03 -8.365020e-03
  [121,] -1.024293e-01 -9.672053e-03 -1.807651e-02  1.266986e-02
  [122,] -3.287871e-02 -1.957521e-02  1.148641e-02  8.209290e-03
  [123,] -8.437919e-02 -3.085894e-02 -7.320048e-02 -7.531646e-03
  [124,] -1.279990e-02 -2.265401e-02 -1.877593e-02 -1.000213e-03
  [125,] -4.901654e-03 -2.276348e-02 -1.678278e-04  4.808312e-04
  [126,] -4.248618e-02  1.463944e-02  1.050126e-01 -1.450264e-02
  [127,] -6.031026e-03 -2.049281e-02 -5.566195e-03 -2.615596e-04
  [128,] -2.595227e-02 -3.528600e-02  3.943811e-03 -6.154860e-03
  [129,] -4.759509e-02 -1.170059e-02  2.196627e-03  2.822998e-03
  [130,] -5.055345e-02  2.802459e-02  1.378974e-02 -5.456092e-04
  [131,] -1.866795e-02  1.678194e-02 -3.078759e-04 -1.378224e-03
  [132,] -3.261158e-02 -3.676107e-02  1.309220e-02 -4.771305e-04
  [133,] -5.439757e-02 -1.021215e-01 -2.656188e-02 -4.819435e-03
  [134,] -2.356935e-02 -1.916713e-02  1.925748e-02 -5.308268e-03
  [135,] -3.079070e-02  2.077196e-02  4.499012e-02 -6.710611e-03
  [136,] -5.708655e-02  4.972869e-02  7.184215e-02 -2.231945e-02
  [137,] -1.019425e-01  4.847508e-02 -3.847270e-02  4.258721e-03
  [138,] -1.832155e-02 -2.267113e-03 -5.006751e-03 -5.631770e-03
  [139,] -3.664292e-02  3.806276e-04 -1.108760e-02 -1.630553e-02
  [140,] -8.094367e-03 -1.028729e-03 -3.303475e-02 -5.069801e-04
  [141,] -6.594732e-02  3.907109e-02  6.984365e-02 -9.476470e-04
  [142,] -6.668858e-02  5.668065e-02 -1.628287e-02  1.208773e-02
  [143,] -4.206499e-03 -1.506893e-02 -2.721529e-03  1.077550e-04
  [144,] -8.558673e-02 -2.300817e-02 -5.063769e-02 -3.403871e-02
  [145,] -3.318418e-02 -4.225689e-02  1.728612e-03  7.763194e-03
  [146,] -9.203687e-02  6.660068e-02  1.394894e-02  1.927595e-02
  [147,] -1.070692e-01  2.789893e-02 -5.492786e-02  3.706702e-02
  [148,] -2.749150e-02  1.159219e-03  2.687778e-02 -1.870916e-02
  [149,] -5.926400e-02 -1.816379e-02  7.253285e-03  7.383675e-03
  [150,] -3.697872e-02 -3.252213e-02 -4.221438e-02 -5.228358e-03
  [151,] -1.227120e-02 -1.786422e-02 -1.738471e-02  7.018426e-04
  [152,] -3.894960e-02 -2.292149e-03  4.932718e-02  1.636273e-04
  [153,] -3.450812e-02 -4.562512e-03 -2.046947e-02  2.212597e-03
  [154,] -6.529797e-02  3.311195e-02  2.644397e-02 -8.001824e-03
  [155,] -5.712193e-03 -5.335677e-02  3.655773e-02 -4.777872e-03
  [156,] -8.283663e-02 -8.192563e-02 -3.063881e-03 -8.565499e-04
  [157,] -2.804421e-02 -1.801020e-02 -1.294603e-02  3.494853e-03
  [158,] -9.783153e-02  2.697166e-02 -3.003954e-02  2.813617e-02
  [159,] -3.583827e-02  2.686174e-02 -1.209046e-02  1.035043e-02
  [160,]  2.394258e-03 -8.501765e-03 -3.210216e-02  3.065908e-03
  [161,] -3.244973e-02  1.177902e-02  2.276373e-02 -1.923711e-02
  [162,] -6.455655e-02  5.386121e-02  7.610154e-02 -1.030458e-02
  [163,] -6.320660e-02  2.686380e-02  1.845503e-02 -8.238160e-04
  [164,] -3.899459e-02 -6.715099e-03  1.769832e-02  6.902570e-04
  [165,] -5.180516e-02 -6.585041e-02 -4.253894e-02 -2.525613e-04
  [166,]  9.186178e-03 -3.955050e-02  4.679200e-02 -2.076545e-03
  [167,] -4.558742e-03 -2.590427e-02  6.580111e-02 -5.199966e-03
  [168,] -2.286590e-02  1.266004e-02  1.592620e-02 -3.315910e-03
  [169,] -4.071682e-02 -1.376813e-02  3.883857e-02 -2.299513e-02
  [170,] -5.164300e-02  9.138117e-02  2.875791e-02 -5.398733e-03
  [171,] -6.014745e-03 -7.829477e-03  2.312790e-02 -3.919874e-03
  [172,] -1.027898e-02 -8.575575e-03  9.391235e-03  2.915216e-04
  [173,] -5.972204e-02 -5.110734e-02 -4.810614e-03 -1.357379e-03
  [174,] -2.361763e-03 -5.460467e-02  1.220337e-02 -1.824380e-03
  [175,] -1.109603e-02 -6.481199e-03  2.885963e-02 -9.391451e-03
  [176,] -2.085691e-02 -3.140476e-02 -1.640682e-02  8.133497e-03
  [177,] -3.198990e-02 -4.610643e-02  3.001831e-03  1.913660e-03
  [178,] -6.797427e-02 -6.986600e-02  1.315622e-02 -3.223414e-02
  [179,] -4.350937e-02 -1.509968e-02 -9.107386e-02 -6.658656e-03
  [180,] -1.027952e-01 -3.253752e-02 -8.623396e-02 -2.995532e-02
  [181,]  1.412810e-02  4.054313e-03 -4.857776e-02  5.333931e-03
  [182,] -1.527361e-02 -1.720390e-02  2.070499e-02  1.015362e-03
  [183,] -4.371945e-02  7.749547e-03  4.990130e-02 -1.593786e-02
  [184,]  1.001995e-02  3.863466e-03  3.198814e-02 -4.416777e-04
  [185,] -2.444499e-02  1.379155e-02 -9.607437e-03  3.918449e-03
  [186,] -2.194199e-02 -1.189303e-01  9.220356e-02 -9.458558e-03
  [187,] -3.934352e-02 -3.837355e-02  1.662724e-02 -8.024713e-03
  [188,] -2.049071e-02 -4.602697e-02 -7.037320e-02 -9.062653e-03
  [189,] -6.073740e-02  1.550460e-02  4.879270e-02 -1.864175e-02
  [190,] -2.236704e-02  1.430154e-02  1.993675e-05 -1.306738e-03
  [191,] -1.184641e-02 -6.335073e-02 -2.301777e-02  7.476655e-04
  [192,] -1.497688e-02 -5.149985e-02  5.526993e-02 -7.508763e-03
  [193,] -5.593903e-02 -4.733795e-02 -7.504198e-02 -1.586859e-02
  [194,] -2.152160e-02  1.808834e-03 -2.165929e-02  6.128455e-03
  [195,] -1.698890e-02  9.344372e-03 -3.220171e-04 -4.904276e-03
  [196,] -5.595937e-02  4.002420e-02  3.392752e-02  4.874193e-03
  [197,] -3.894996e-02  2.132212e-02  1.608052e-02 -2.168308e-02
  [198,] -1.930067e-02 -2.785546e-02 -2.223686e-02 -9.919391e-04
  [199,] -7.370431e-02  4.208921e-02 -3.360574e-02 -1.542901e-02
  [200,] -1.915113e-02 -4.192992e-02 -6.516708e-02  7.721800e-03
  [201,] -4.411583e-02  1.086388e-02 -8.123969e-02  2.549212e-02
  [202,] -1.096397e-01 -5.297081e-02  1.237093e-02 -1.941775e-02
  [203,] -5.379034e-02 -1.999258e-02  3.659456e-02 -6.623743e-03
  [204,] -3.519401e-02 -8.818417e-02 -5.973625e-02 -9.134005e-03
  [205,] -3.666824e-02  4.278704e-02  2.906576e-02 -5.248749e-03
  [206,] -8.511653e-03 -7.528271e-03  3.339333e-02 -1.538053e-03
  [207,] -6.657282e-02  2.252116e-02  1.818680e-02 -1.551829e-03
  [208,] -6.248050e-02  4.289367e-02 -1.151028e-01 -2.470859e-03
  [209,] -2.693295e-02  9.097293e-03  4.556648e-02 -1.195049e-03
  [210,] -1.081570e-02 -9.567836e-03 -1.417787e-02  1.162248e-03
  [211,] -2.964156e-02  1.525637e-02  8.414777e-03  6.335921e-03
  [212,] -1.300156e-02 -1.807468e-02  3.719239e-04 -1.384583e-03
  [213,] -6.995742e-02 -5.346305e-02  2.113596e-03 -7.439445e-03
  [214,] -3.969950e-02 -1.864492e-02 -3.387357e-02 -2.548330e-02
  [215,] -1.260958e-02 -2.612694e-03  9.492052e-03 -1.717066e-03
  [216,] -2.017536e-02 -6.298267e-02  4.241853e-02 -4.384336e-03
  [217,] -5.628708e-02 -1.477652e-02  4.578590e-02 -1.084562e-03
  [218,] -6.673904e-03 -5.489311e-02  6.238443e-02 -3.415006e-03
  [219,] -1.289090e-01  2.430175e-02  1.042268e-01 -1.643290e-02
  [220,] -9.192756e-03 -1.766507e-02  2.862192e-02 -5.489300e-03
  [221,] -8.016749e-04 -9.057790e-02 -3.104781e-02  7.705154e-04
  [222,] -7.490476e-02  4.734040e-02 -1.159565e-02 -2.318762e-02
  [223,] -7.225676e-02  3.608268e-02 -3.176951e-02  2.531870e-02
  [224,] -2.070960e-02  1.250478e-02 -3.481682e-03  1.283381e-02
  [225,] -2.213838e-02 -7.343891e-02  2.264886e-02 -1.960560e-02
  [226,] -4.535377e-02 -1.148536e-01  4.802006e-02 -1.679356e-02
  [227,] -1.903305e-02 -3.796853e-02 -1.655957e-02 -5.334384e-03
  [228,] -6.441090e-02 -1.452070e-03 -1.398831e-02  1.341950e-02
  [229,] -3.146234e-02 -1.926136e-02 -4.536452e-02  4.473772e-03
  [230,] -3.581034e-02 -1.305869e-01 -3.945765e-02  1.738876e-03
  [231,] -1.525276e-02 -5.784764e-02  4.307642e-02 -3.571577e-03
  [232,] -4.600217e-02  7.867258e-02  8.525275e-02  6.476166e-05
  [233,] -2.657513e-02 -8.611680e-02  4.965390e-02  1.566540e-03
  [234,] -6.163073e-03 -4.667851e-04  1.012281e-01 -2.538785e-04
  [235,] -5.175332e-02  2.224813e-02  6.067861e-02 -1.652000e-03
  [236,] -2.422966e-02 -1.073150e-02 -2.663842e-02 -9.041796e-03
  [237,] -5.829658e-02 -1.743012e-02 -5.322201e-02  4.107571e-03
  [238,] -1.948009e-02  8.842907e-03 -2.470518e-02  8.236023e-03
  [239,] -3.108283e-02  9.432590e-03 -3.191270e-03  2.845065e-03
  [240,] -6.768226e-02 -1.957472e-01  3.403306e-02 -2.385622e-02
  [241,] -3.678041e-02 -2.395011e-02  1.552527e-02 -2.092474e-03
  [242,] -5.563989e-04 -4.929175e-02  7.460674e-02  1.375401e-03
  [243,] -2.362569e-02 -1.131748e-01 -1.414148e-02 -3.312493e-03
  [244,] -5.554352e-03 -7.674257e-02  2.824478e-02  2.639362e-03
  [245,] -5.554406e-02  1.823257e-02 -4.878659e-02 -8.489058e-03
  [246,] -7.107421e-02 -3.520741e-02  1.044450e-01 -4.503281e-02
  [247,] -3.925650e-02  2.595291e-03 -2.262757e-02  1.010754e-02
  [248,] -1.863395e-02  1.600962e-02 -5.068079e-02 -9.248620e-03
  [249,] -3.092998e-02  4.762461e-02 -2.384063e-02  8.908142e-04
  [250,] -7.194577e-02  1.193110e-02 -1.289757e-02  1.189094e-03
  [251,] -2.899415e-02 -1.132218e-02 -5.856472e-02  4.178033e-03
  [252,] -6.197903e-02  7.567210e-02  6.552628e-02 -1.337652e-02
  [253,] -2.130341e-02 -3.292741e-03  7.620134e-03 -9.973180e-04
  [254,] -1.428580e-02  1.657458e-02  6.311040e-02  1.676783e-03
  [255,]  8.709728e-04 -1.503162e-01  4.327788e-02 -1.079782e-02
  [256,] -5.659593e-03 -3.344874e-03  2.901675e-02 -1.188000e-03
  [257,] -5.160152e-02  3.242115e-02 -1.234843e-01 -4.701174e-03
  [258,] -1.305137e-02 -3.056467e-02  3.341604e-02 -2.939021e-03
  [259,] -2.461906e-02 -1.258189e-03 -1.887528e-03 -3.405722e-03
  [260,] -3.854153e-02 -1.872729e-02 -3.905392e-02 -3.515800e-03
  [261,] -1.831501e-02 -2.728302e-02  4.956613e-03 -1.184916e-02
  [262,] -1.494347e-01 -1.843941e-04 -1.292694e-01 -4.207169e-02
  [263,] -2.951041e-02  2.629930e-03  4.155416e-02 -4.265589e-03
  [264,] -1.336545e-01  2.657894e-02  6.522149e-03 -4.227429e-02
  [265,] -4.934039e-03  2.991429e-02  3.324445e-02  3.397997e-04
  [266,] -2.898777e-02 -7.469691e-02 -1.213240e-02  3.184132e-03
  [267,] -5.398598e-02  1.938598e-02 -9.155443e-03 -3.081385e-02
  [268,] -7.459630e-03 -1.015216e-01  1.949026e-02 -3.454782e-03
  [269,] -3.017370e-02 -7.248694e-02  1.391787e-02 -1.564610e-02
  [270,] -9.620541e-02  8.082247e-04 -6.305316e-02 -7.755108e-05
  [271,] -4.972789e-02 -4.388410e-02 -7.130383e-02 -3.260810e-03
  [272,] -4.521245e-02 -6.315748e-02 -8.680367e-02  1.841245e-02
  [273,] -5.075743e-02 -5.516967e-03  1.879975e-02 -8.967307e-03
  [274,] -9.358747e-03 -2.116880e-02  4.393214e-02 -3.266554e-03
  [275,] -2.180919e-02  2.580200e-03  2.300501e-02 -1.403203e-02
  [276,] -6.136094e-02  6.068982e-02 -2.197298e-02 -2.560135e-02
  [277,] -6.261729e-02  3.822235e-02 -9.278627e-02  1.468150e-02
  [278,] -8.658815e-02  5.159019e-02 -1.787297e-02  1.462369e-02
  [279,] -3.096116e-02  9.359060e-02  4.474163e-02  3.625729e-03
  [280,] -7.734122e-02  8.894113e-02  1.438078e-02  4.376679e-02
  [281,] -5.413849e-02  2.446765e-02 -2.604912e-02 -2.484609e-02
  [282,] -4.754232e-02 -1.044231e-01 -2.020951e-02  2.096403e-02
  [283,] -7.072496e-02 -7.448060e-02  1.878713e-02 -1.985749e-02
  [284,] -8.138073e-02 -1.501203e-02 -5.198766e-02 -2.470063e-02
  [285,] -1.648745e-02  3.592116e-02  1.764274e-02 -1.691186e-03
  [286,] -4.399227e-02 -6.235776e-03 -5.161910e-02 -1.901005e-02
  [287,] -3.003215e-02 -2.002525e-03  3.339054e-02 -1.274447e-03
  [288,] -5.860430e-02  1.020043e-02  1.693348e-04 -4.408552e-03
  [289,] -1.573618e-02  1.083140e-02  7.669964e-03  1.013510e-04
  [290,] -4.301013e-02  2.573165e-03  5.353597e-02 -5.744956e-03
  [291,] -4.843326e-03 -6.998220e-03 -6.962551e-03 -4.795579e-04
  [292,] -6.438679e-02  3.594286e-02  2.562909e-03 -3.301589e-03
  [293,] -8.498680e-02 -5.739979e-02 -3.504037e-02 -1.790870e-02
  [294,] -1.529347e-02  3.758544e-02  6.476570e-03  5.345511e-03
  [295,] -4.527944e-02 -6.351516e-03 -2.060346e-02  2.498336e-02
  [296,] -3.780408e-02  1.232756e-02 -9.275812e-02 -1.200272e-02
  [297,] -3.145817e-02  1.646364e-02  1.423451e-02  1.920632e-03
  [298,] -2.594832e-02  2.129636e-02  3.651587e-02 -9.418202e-03
  [299,] -1.953740e-02 -5.966470e-02  6.102226e-02 -1.273311e-02
  [300,] -5.122941e-02  2.811277e-02 -5.050905e-02 -1.464628e-02
  [301,] -4.278718e-02 -9.804904e-03  3.272999e-02 -6.883365e-03
  [302,] -3.134992e-02  3.429021e-02 -1.209084e-02  7.528501e-03
  [303,] -8.254475e-02  1.715904e-02 -5.269917e-02 -3.417340e-02
  [304,] -6.668664e-03 -2.165284e-02 -1.902433e-02 -1.907093e-03
  [305,] -4.421635e-02 -5.246383e-03 -4.517584e-02  2.990775e-03
  [306,] -2.466520e-02  1.490403e-02  8.668411e-02 -1.483886e-02
  [307,] -3.683705e-02 -5.910366e-03  2.616663e-02  1.142692e-02
  [308,] -5.641232e-02 -3.148186e-02  1.231868e-02 -1.591494e-02
  [309,] -5.322292e-02  4.216621e-02  3.669581e-02 -2.422498e-02
  [310,] -3.691664e-02 -6.205112e-03  2.732569e-02 -1.208206e-02
  [311,] -1.758957e-02 -5.815489e-02  7.375433e-02 -1.239010e-02
  [312,] -8.385754e-02  2.384780e-02  4.871908e-02  1.956886e-02
  [313,] -8.465547e-02  5.182535e-02 -1.947229e-02  3.428696e-02
  [314,] -3.804232e-02 -8.726213e-03  1.128497e-01 -2.934666e-02
  [315,] -8.499840e-02  1.238028e-02  1.272055e-01 -3.419046e-02
  [316,] -2.183396e-02 -8.024635e-02  3.885123e-02  1.774968e-03
  [317,] -4.531747e-02 -1.061714e-02 -5.900473e-02 -9.302751e-03
  [318,] -6.048853e-02  3.045381e-02  9.729293e-02  1.337495e-02
  [319,] -2.356201e-02 -1.508327e-03  3.043095e-03 -1.564438e-03
  [320,] -4.264073e-02  5.064611e-02  8.336961e-03 -2.371501e-02
  [321,] -4.128828e-02 -7.875275e-03  4.183769e-02  5.006815e-03
  [322,] -4.431895e-02 -3.921915e-02 -3.891150e-02  1.195451e-02
  [323,] -7.791088e-02  1.167675e-02 -2.173969e-02  3.144526e-02
  [324,] -1.537851e-02 -8.536198e-03  9.741226e-02  2.604359e-03
  [325,] -2.099949e-02 -1.717298e-02  1.975099e-02 -7.067971e-03
  [326,] -6.635081e-02 -2.977470e-02  1.733887e-02 -2.873325e-02
  [327,] -4.576194e-02 -1.343011e-02  4.480417e-02  4.573212e-03
  [328,] -4.115753e-02 -7.764880e-02 -4.285037e-02 -6.487537e-03
  [329,] -6.459137e-02  5.757184e-03 -1.001207e-02 -9.713808e-03
  [330,] -6.290945e-03 -1.633522e-02  5.788799e-03 -4.239075e-04
  [331,] -2.267683e-02 -9.374380e-03  5.900663e-02 -1.073891e-03
  [332,] -3.248386e-02  8.986796e-02 -1.116655e-02 -4.770633e-03
  [333,] -6.257607e-02  6.406738e-02 -2.928987e-03 -2.078700e-02
  [334,] -2.615374e-02 -5.923579e-02  5.741242e-02 -1.393871e-02
  [335,] -7.412438e-02 -1.932241e-03  7.319213e-03  1.951964e-02
  [336,] -3.770569e-02  1.433436e-03  1.863038e-02 -3.561854e-03
  [337,] -5.071387e-02  2.647510e-02 -1.176834e-02 -1.736611e-02
  [338,] -1.272608e-02 -7.034409e-02  6.469605e-02 -6.360717e-03
  [339,] -2.549624e-02 -3.830539e-02  7.015472e-03 -1.043732e-02
  [340,] -4.929607e-02  7.230942e-02  8.545226e-02 -1.668970e-03
  [341,] -2.744733e-02  2.979540e-02  1.389468e-02 -2.283211e-03
  [342,] -1.429738e-02 -1.870878e-02 -1.081525e-03 -5.184195e-05
  [343,] -5.937606e-02 -2.076992e-02 -7.085722e-02 -1.830192e-02
  [344,] -4.809524e-02  1.147012e-02 -3.787437e-02  1.864370e-02
  [345,] -3.424896e-02 -4.372305e-02 -1.294051e-02 -1.007518e-02
  [346,] -3.898757e-02  8.785760e-03 -3.661667e-03  2.967948e-03
  [347,] -5.443720e-02 -4.408729e-02  1.183413e-03 -3.080599e-04
  [348,] -4.857582e-02  5.845490e-02  1.571823e-02 -4.305236e-03
  [349,] -6.538776e-02 -4.263069e-03 -2.493001e-03  9.325508e-04
  [350,] -1.192816e-02  9.682046e-03  1.940906e-02 -3.660457e-03
  [351,] -2.361120e-02  2.376154e-02  1.098356e-02  2.077725e-03
  [352,] -1.553103e-02  2.998943e-02  4.512373e-02  2.357566e-03
  [353,] -4.189652e-02  3.582770e-02 -8.376900e-03 -3.677179e-03
  [354,] -1.288688e-02 -1.483535e-02 -2.950355e-02  4.289471e-03
  [355,] -2.191034e-02 -5.128992e-02  8.417251e-02 -2.503558e-03
  [356,] -8.562904e-02 -1.292450e-02 -3.488205e-02  2.618776e-02
  [357,] -1.632307e-02 -1.248055e-02  1.929317e-02 -1.885260e-04
  [358,] -4.908731e-02 -4.140304e-02  9.271987e-03 -3.009009e-02
  [359,] -2.586947e-02 -5.740074e-03  4.694698e-02 -1.430592e-03
  [360,]  4.433342e-04 -1.004582e-02  5.178645e-02 -6.576769e-03
  [361,] -7.358653e-02  8.500840e-02 -3.283247e-02  4.042699e-02
  [362,] -1.475875e-02  2.614756e-03  2.105343e-02 -2.375928e-03
  [363,] -4.032746e-02  2.530770e-03 -1.401339e-02 -1.351233e-02
  [364,] -1.372702e-02 -2.831855e-02 -2.573715e-02 -1.690574e-03
  [365,] -4.083442e-02  2.278724e-02 -4.109229e-02  2.136057e-02
  [366,] -1.816157e-02  4.099709e-02 -3.350261e-02  5.668651e-03
  [367,]  1.644769e-03 -1.151421e-03  3.791938e-02 -3.174616e-03
  [368,] -4.016612e-02  2.882321e-03  3.427528e-02 -1.614957e-02
  [369,] -2.191312e-02 -1.425212e-02 -1.941413e-02  4.558717e-03
  [370,] -5.530318e-02  1.177407e-01  6.211652e-02  6.741826e-04
  [371,] -5.272918e-02  7.469846e-02  3.762829e-02 -5.046651e-03
  [372,] -5.258757e-02  1.287711e-02  2.872232e-02 -1.356235e-02
  [373,] -1.512539e-03 -1.219634e-01 -3.230942e-02 -9.603281e-03
  [374,] -3.713171e-02 -6.795724e-02 -2.778471e-02  1.453628e-02
  [375,] -9.725550e-03 -3.190712e-03  8.031352e-02 -4.596172e-04
  [376,] -3.617062e-02 -1.041808e-01 -2.850776e-02 -1.391916e-02
  [377,] -5.693887e-03 -3.670434e-02 -6.031915e-02  2.166295e-03
  [378,] -6.437123e-03 -4.885608e-02  7.226368e-03 -4.323586e-03
  [379,] -8.650511e-02 -4.625479e-02 -4.322547e-03  1.247603e-02
  [380,] -2.931102e-02  8.717223e-03 -4.171396e-02  8.894509e-03
  [381,] -1.086386e-02  7.716500e-02  1.218812e-01 -1.028487e-02
  [382,] -5.040648e-02  2.866045e-02  6.425793e-02 -6.567353e-03
  [383,] -4.465346e-02  5.071517e-02 -2.699601e-04  6.200191e-03
  [384,] -1.081784e-02 -7.621784e-02  6.005053e-03 -8.574521e-04
  [385,] -1.184057e-02 -3.493647e-02  8.311654e-02 -8.844748e-03
  [386,] -2.466242e-02 -2.066364e-02  2.031268e-02 -1.385682e-04
  [387,] -3.758868e-02  4.279366e-02  2.892744e-03 -3.640281e-03
  [388,] -7.478452e-02  1.141824e-02 -3.965078e-02 -1.182196e-03
  [389,] -8.345660e-02  2.340481e-03 -3.181552e-02 -3.221246e-02
  [390,] -5.003337e-02  3.996778e-02 -2.458192e-02  1.021108e-02
  [391,] -6.247898e-02  3.833985e-02 -7.478718e-02 -1.201176e-02
  [392,] -1.022100e-02 -2.640025e-03  4.744700e-03 -5.492798e-03
  [393,] -1.011139e-02 -2.772611e-02  5.520388e-02 -5.064201e-03
  [394,] -5.926944e-02  1.056903e-02  1.894642e-02 -3.000751e-04
  [395,] -7.440942e-03  5.721839e-03  5.045968e-03 -1.645970e-03
  [396,] -3.882907e-02  3.083158e-02 -4.437205e-02  1.129928e-02
  [397,] -1.064886e-02 -2.260557e-02  7.830090e-03 -6.921806e-03
  [398,] -7.134201e-02  4.904892e-03 -1.055111e-01  4.512171e-02
  [399,] -2.884536e-02  5.380575e-02 -3.342506e-02 -3.246811e-03
  [400,] -5.466651e-02 -1.778462e-03  1.095193e-02 -3.258524e-02
  [401,] -2.901218e-02 -4.445882e-05  1.535670e-02 -1.368544e-03
  [402,] -1.744394e-02  2.650936e-02  4.749336e-03 -8.914894e-03
  [403,] -2.129297e-02 -1.904367e-02 -6.982008e-02  3.445991e-03
  [404,] -4.625607e-02  4.453984e-02 -6.016468e-02 -9.374975e-03
  [405,] -2.178641e-02 -9.053372e-02 -2.291831e-02  1.076606e-02
  [406,] -7.720799e-02  2.108760e-02 -1.455489e-02 -6.505998e-03
  [407,] -6.273839e-02  1.612627e-02 -2.396582e-02 -1.391460e-02
  [408,] -2.232631e-02 -4.001581e-02 -6.068855e-03 -2.225132e-04
  [409,] -3.936296e-02  4.158902e-02 -1.433695e-02 -1.656674e-02
  [410,] -1.087366e-02 -8.239133e-02  6.125838e-02  1.817342e-04
  [411,] -7.762519e-02 -3.512382e-03  7.268047e-02  1.107036e-03
  [412,] -5.216486e-02 -5.047190e-02 -6.710261e-03  3.786038e-03
  [413,] -1.688645e-02  4.958313e-03  2.988260e-02 -6.990672e-03
  [414,] -2.849950e-02 -2.110360e-02  7.786997e-03 -1.178633e-02
  [415,] -1.063276e-02  1.858688e-02  2.203427e-03 -4.298729e-03
  [416,] -3.718529e-02  1.553457e-02 -1.459568e-02 -1.174422e-02
  [417,] -6.680531e-02 -6.475310e-02  8.997765e-02 -8.699687e-03
  [418,] -1.978828e-02 -3.601715e-02 -7.455370e-02 -7.606844e-04
  [419,] -1.713474e-02 -5.446252e-02  4.165601e-02 -1.197809e-02
  [420,] -2.179630e-02 -1.892249e-02  8.610978e-03 -2.063002e-03
  [421,] -4.544582e-02  6.087251e-03 -3.124612e-02 -7.968048e-03
  [422,] -6.083550e-02  6.053051e-02  7.711302e-02 -2.666618e-02
  [423,] -3.268555e-02  1.368593e-02 -7.659948e-02 -1.717925e-03
  [424,] -1.790135e-02 -2.118306e-03  2.976639e-03  2.606766e-03
  [425,] -7.895966e-02 -5.056219e-03  1.524076e-02 -7.847702e-03
  [426,] -7.702789e-03 -4.418395e-02  2.821600e-02 -3.741899e-03
  [427,] -3.651183e-02 -2.758935e-02  3.663937e-02  1.597299e-03
  [428,]  4.122708e-03  4.446700e-02  4.344992e-02 -5.686466e-03
  [429,] -5.400188e-02 -5.033559e-02  3.079198e-02  9.057186e-04
  [430,] -1.661446e-02 -2.993271e-02 -5.083537e-03  1.100905e-03
  [431,] -2.338550e-02 -1.064157e-01  9.183444e-03 -7.699808e-05
  [432,] -1.637521e-02  1.146829e-02 -3.988406e-02  8.082053e-04
  [433,]  2.093227e-03 -4.208800e-02  7.286710e-02  1.932019e-04
  [434,] -3.537604e-02 -1.100327e-02  3.092757e-03  2.387426e-03
  [435,] -3.661580e-02 -1.864947e-02 -3.388084e-02 -7.180598e-03
  [436,] -2.134794e-02 -8.314661e-03 -1.264344e-02  9.813677e-04
  [437,] -5.878444e-02  6.068336e-02 -3.339747e-02 -7.472261e-04
  [438,] -5.429675e-02  2.268021e-02 -1.000413e-03 -2.507479e-02
  [439,] -2.720793e-02  3.499552e-02  2.070871e-02  5.497973e-04
  [440,] -2.389663e-02 -4.879462e-02  1.438008e-02 -9.231296e-03
  [441,] -4.328022e-02 -2.267775e-02  2.750114e-02 -2.996724e-03
  [442,] -8.451624e-03 -6.449157e-02 -2.603447e-02 -3.857424e-03
  [443,] -7.659991e-07 -3.500294e-02  1.888396e-02 -3.505380e-03
  [444,] -2.690646e-02  2.539693e-02  1.397515e-01  6.187591e-03
  [445,] -3.369152e-02  3.358770e-02  7.596626e-02 -1.352358e-02
  [446,] -5.743780e-02  5.450906e-02 -2.145027e-02 -7.239489e-03
  [447,] -3.008201e-02  1.127970e-02  3.833359e-02 -3.363307e-03
  [448,] -2.418537e-02 -2.596796e-02  9.286824e-03 -3.190293e-03
  [449,] -4.159206e-02  3.248049e-03  3.592653e-02  1.059751e-02
  [450,] -8.844774e-02  5.089348e-02  7.049905e-04  2.929222e-03
  [451,] -3.765697e-02  1.020651e-02 -3.607024e-02  4.602829e-03
  [452,] -2.570730e-02 -2.813164e-02 -3.011045e-02 -1.494032e-03
  [453,] -1.112428e-02  1.160217e-03  3.729016e-02 -2.126338e-03
  [454,] -6.569410e-02  2.304070e-02 -8.496383e-02  1.019605e-02
  [455,] -2.061535e-03 -6.516327e-02  1.033699e-01 -1.238998e-02
  [456,] -6.885846e-02  2.523523e-02  5.721027e-02  8.301633e-03
  [457,] -4.100282e-02 -7.369351e-04 -6.072719e-02 -9.496816e-03
  [458,] -7.065030e-02 -7.035409e-02 -1.011705e-01  3.863118e-02
  [459,] -5.830450e-02  7.912572e-03  7.267680e-02 -2.355054e-02
  [460,] -2.456898e-02 -3.539776e-02  2.019639e-02 -6.569821e-03
  [461,] -3.035927e-02 -2.671516e-02  2.743905e-02 -1.582992e-02
  [462,] -3.027915e-02 -1.167586e-02  2.305771e-02 -1.117388e-02
  [463,] -6.521135e-02  1.546358e-02  3.931978e-02 -7.679907e-03
  [464,] -6.307924e-02 -3.986174e-03 -1.744659e-02  7.043012e-03
  [465,] -3.171831e-02  4.147859e-02 -1.437467e-03  9.806475e-03
  [466,] -8.282001e-02 -2.598861e-02 -1.191283e-02 -1.805805e-03
  [467,] -3.296934e-02 -7.631440e-04  7.250069e-02 -9.621228e-04
  [468,] -5.993873e-02  1.632836e-02 -6.456417e-02 -4.275273e-03
  [469,] -2.012540e-02 -5.324960e-02  1.559411e-02 -1.150920e-02
  [470,] -2.355414e-02  3.089503e-02 -1.639357e-02  6.620056e-03
  [471,] -4.206828e-02  5.889884e-02  1.348085e-02 -2.513953e-02
  [472,] -2.828940e-02 -3.775118e-02  2.947612e-02 -7.215256e-03
  [473,] -3.565157e-02 -7.521361e-03 -1.124893e-02 -1.759576e-02
  [474,] -9.532648e-03 -1.043933e-01 -4.129662e-02 -5.576737e-03
  [475,] -5.130801e-02 -3.259609e-02  3.687652e-02 -1.471237e-02
  [476,] -9.266476e-02 -3.413661e-03 -3.725817e-02 -3.018297e-02
  [477,] -4.727947e-02  7.439780e-04 -7.792518e-02 -1.667112e-02
  [478,] -1.350009e-02  2.355750e-02 -2.126486e-02  7.980522e-03
  [479,] -6.179868e-02 -6.875776e-02  3.424548e-02 -1.908315e-02
  [480,] -1.903872e-02 -1.374978e-01  1.242736e-01 -1.375474e-02
  [481,] -4.517951e-02 -2.032934e-02 -3.885745e-02 -9.513343e-03
  [482,] -3.696850e-02 -2.888583e-02  2.327650e-02 -3.132146e-03
  [483,] -4.836677e-02 -6.077138e-02 -4.802265e-02 -2.017598e-02
  [484,] -6.045000e-03 -3.615100e-02  2.256619e-02 -3.505473e-03
  [485,] -1.492059e-02 -4.809589e-03 -3.130675e-02 -1.101080e-03
  [486,] -1.606335e-02 -6.266873e-02  5.170091e-02 -7.687562e-03
  [487,] -3.698408e-02 -2.188976e-02 -9.713973e-03  8.683240e-03
  [488,] -2.629546e-02 -6.007535e-02 -3.906966e-02  7.632136e-03
  [489,] -3.830816e-02 -1.488244e-02 -2.736894e-03 -1.726482e-03
  [490,] -8.957217e-02  5.080250e-02 -5.585287e-02 -1.568902e-02
  [491,] -5.733182e-03 -9.956114e-03  7.966596e-03 -2.692559e-03
  [492,] -4.387702e-02 -1.402409e-02 -5.661580e-02  5.615252e-03
  [493,] -5.462598e-02  4.552289e-02  6.241694e-02 -4.248524e-03
  [494,] -3.193835e-02 -2.558365e-02  1.905417e-03  2.399014e-03
  [495,] -3.339690e-02  4.831719e-02  9.941972e-03 -1.270186e-02
  [496,] -7.985035e-03 -2.060099e-02  4.985037e-03 -2.938623e-03
  [497,] -8.988853e-03  3.397693e-02  2.804638e-03 -6.084407e-03
  [498,] -2.202186e-02 -3.915148e-02  2.839721e-03 -3.902339e-03

  $variance_explained
     PC variance_explained
  1 PC1       9.534876e+01
  2 PC2       2.912731e+00
  3 PC3       1.738511e+00
  4 PC4       1.510939e-30

can prepare pca with log transform

Code
  prepare_pca(df = test_count_matrix[3:6], log_transform = T)
Output
  $data
  Standard deviations (1, .., p=4):
  [1] 1.6429460 0.5557668 0.3476552 0.3018909

  Rotation (n x k) = (4 x 4):
              PC1        PC2         PC3           PC4
  T0_- -0.4557836 -0.8892568 -0.03850871 -0.0008873807
  T3_A -0.4798683  0.2626174 -0.40172521  0.7344218120
  T3_B -0.5236962  0.2875918 -0.42714101 -0.6786632635
  T3_C -0.5364048  0.2398863  0.80912576  0.0063237126

  $variance_explained
     PC variance_explained
  1 PC1          83.824392
  2 PC2           9.591995
  3 PC3           3.753369
  4 PC4           2.830244

can plot pca

Code
  p$data
Output
                         PC1        PC2        PC3        PC4         PC5
  T0_-          -33.62613183  -6.439070  0.7232680  0.4286996  -3.5860350
  T3_A           -0.21900736  -9.401945 -0.4013321 -1.0314620   0.4772895
  T3_B           14.57571048 -11.002677 -1.4252431  0.9706678   0.9382381
  T3_C            9.69563648 -10.394493 -0.9529347  0.3893537  -1.1364381
  T18_A          -0.85395841   7.365705  1.0010557 17.8694881   2.1329670
  T18_B           0.09028665   5.103602 -7.4847275 -2.9933557  13.2505468
  T18_C          -4.35851802   4.945270 -0.7934075 -6.3873665   5.8065618
  T18_A_Starved   8.08879125   6.198627 -9.2802264 -0.2555649 -12.0908182
  T18_B_Starved  -0.43986752  10.177639  0.5057274 -6.5955674  -5.5316478
  T18_C_Starved   7.04705828   3.447342 18.1078202 -2.3948928  -0.2606642
                        PC6        PC7        PC8         PC9         PC10
  T0_-            2.2636988 -1.6046747 -1.6375648  2.10273388 1.125142e-14
  T3_A           -0.8727215 -1.1261056  8.0739865 -5.48778331 8.735200e-15
  T3_B           -2.4099736 -0.1187338  1.5396611  7.89418460 7.202572e-15
  T3_C           -3.0574301  0.9170718 -7.7374868 -4.59650194 1.505393e-14
  T18_A          -3.2147517  1.4402947  0.4969632 -0.46250920 1.333113e-14
  T18_B           5.1641489 -7.7415508 -1.2311043 -0.15208017 1.286818e-14
  T18_C          -1.5394861 14.0452857  0.2371402  0.29894270 1.474515e-14
  T18_A_Starved   9.4109694  1.4096994  0.4248301 -0.07558311 1.113692e-14
  T18_B_Starved -12.7069309 -5.5116097  0.1887208  0.37825982 1.714601e-14
  T18_C_Starved   6.9624769 -1.7096770 -0.3551460  0.10033674 1.421259e-14
Code
  p$layers
Output
  [[1]]
  mapping: x = ~get(pc_x), y = ~get(pc_y) 
  geom_point: na.rm = FALSE
  stat_identity: na.rm = FALSE
  position_identity

can plot pca with color

Code
  p$data
Output
                         PC1        PC2        PC3        PC4         PC5
  T0_-          -33.62613183  -6.439070  0.7232680  0.4286996  -3.5860350
  T3_A           -0.21900736  -9.401945 -0.4013321 -1.0314620   0.4772895
  T3_B           14.57571048 -11.002677 -1.4252431  0.9706678   0.9382381
  T3_C            9.69563648 -10.394493 -0.9529347  0.3893537  -1.1364381
  T18_A          -0.85395841   7.365705  1.0010557 17.8694881   2.1329670
  T18_B           0.09028665   5.103602 -7.4847275 -2.9933557  13.2505468
  T18_C          -4.35851802   4.945270 -0.7934075 -6.3873665   5.8065618
  T18_A_Starved   8.08879125   6.198627 -9.2802264 -0.2555649 -12.0908182
  T18_B_Starved  -0.43986752  10.177639  0.5057274 -6.5955674  -5.5316478
  T18_C_Starved   7.04705828   3.447342 18.1078202 -2.3948928  -0.2606642
                        PC6        PC7        PC8         PC9         PC10
  T0_-            2.2636988 -1.6046747 -1.6375648  2.10273388 1.125142e-14
  T3_A           -0.8727215 -1.1261056  8.0739865 -5.48778331 8.735200e-15
  T3_B           -2.4099736 -0.1187338  1.5396611  7.89418460 7.202572e-15
  T3_C           -3.0574301  0.9170718 -7.7374868 -4.59650194 1.505393e-14
  T18_A          -3.2147517  1.4402947  0.4969632 -0.46250920 1.333113e-14
  T18_B           5.1641489 -7.7415508 -1.2311043 -0.15208017 1.286818e-14
  T18_C          -1.5394861 14.0452857  0.2371402  0.29894270 1.474515e-14
  T18_A_Starved   9.4109694  1.4096994  0.4248301 -0.07558311 1.113692e-14
  T18_B_Starved -12.7069309 -5.5116097  0.1887208  0.37825982 1.714601e-14
  T18_C_Starved   6.9624769 -1.7096770 -0.3551460  0.10033674 1.421259e-14
                      color
  T0_-                   T0
  T3_A                   T3
  T3_B                   T3
  T3_C                   T3
  T18_A                 T18
  T18_B                 T18
  T18_C                 T18
  T18_A_Starved T18_Starved
  T18_B_Starved T18_Starved
  T18_C_Starved T18_Starved
Code
  p$layers
Output
  [[1]]
  mapping: x = ~get(pc_x), y = ~get(pc_y), colour = ~color 
  geom_point: na.rm = FALSE
  stat_identity: na.rm = FALSE
  position_identity

can plot pca with shape

Code
  p$data
Output
                         PC1        PC2        PC3        PC4         PC5
  T0_-          -33.62613183  -6.439070  0.7232680  0.4286996  -3.5860350
  T3_A           -0.21900736  -9.401945 -0.4013321 -1.0314620   0.4772895
  T3_B           14.57571048 -11.002677 -1.4252431  0.9706678   0.9382381
  T3_C            9.69563648 -10.394493 -0.9529347  0.3893537  -1.1364381
  T18_A          -0.85395841   7.365705  1.0010557 17.8694881   2.1329670
  T18_B           0.09028665   5.103602 -7.4847275 -2.9933557  13.2505468
  T18_C          -4.35851802   4.945270 -0.7934075 -6.3873665   5.8065618
  T18_A_Starved   8.08879125   6.198627 -9.2802264 -0.2555649 -12.0908182
  T18_B_Starved  -0.43986752  10.177639  0.5057274 -6.5955674  -5.5316478
  T18_C_Starved   7.04705828   3.447342 18.1078202 -2.3948928  -0.2606642
                        PC6        PC7        PC8         PC9         PC10
  T0_-            2.2636988 -1.6046747 -1.6375648  2.10273388 1.125142e-14
  T3_A           -0.8727215 -1.1261056  8.0739865 -5.48778331 8.735200e-15
  T3_B           -2.4099736 -0.1187338  1.5396611  7.89418460 7.202572e-15
  T3_C           -3.0574301  0.9170718 -7.7374868 -4.59650194 1.505393e-14
  T18_A          -3.2147517  1.4402947  0.4969632 -0.46250920 1.333113e-14
  T18_B           5.1641489 -7.7415508 -1.2311043 -0.15208017 1.286818e-14
  T18_C          -1.5394861 14.0452857  0.2371402  0.29894270 1.474515e-14
  T18_A_Starved   9.4109694  1.4096994  0.4248301 -0.07558311 1.113692e-14
  T18_B_Starved -12.7069309 -5.5116097  0.1887208  0.37825982 1.714601e-14
  T18_C_Starved   6.9624769 -1.7096770 -0.3551460  0.10033674 1.421259e-14
                      shape
  T0_-                   T0
  T3_A                   T3
  T3_B                   T3
  T3_C                   T3
  T18_A                 T18
  T18_B                 T18
  T18_C                 T18
  T18_A_Starved T18_Starved
  T18_B_Starved T18_Starved
  T18_C_Starved T18_Starved
Code
  p$layers
Output
  [[1]]
  mapping: x = ~get(pc_x), y = ~get(pc_y), shape = ~shape 
  geom_point: na.rm = FALSE
  stat_identity: na.rm = FALSE
  position_identity

can plot pca with color and shape

Code
  p$data
Output
                         PC1        PC2        PC3        PC4         PC5
  T0_-          -33.62613183  -6.439070  0.7232680  0.4286996  -3.5860350
  T3_A           -0.21900736  -9.401945 -0.4013321 -1.0314620   0.4772895
  T3_B           14.57571048 -11.002677 -1.4252431  0.9706678   0.9382381
  T3_C            9.69563648 -10.394493 -0.9529347  0.3893537  -1.1364381
  T18_A          -0.85395841   7.365705  1.0010557 17.8694881   2.1329670
  T18_B           0.09028665   5.103602 -7.4847275 -2.9933557  13.2505468
  T18_C          -4.35851802   4.945270 -0.7934075 -6.3873665   5.8065618
  T18_A_Starved   8.08879125   6.198627 -9.2802264 -0.2555649 -12.0908182
  T18_B_Starved  -0.43986752  10.177639  0.5057274 -6.5955674  -5.5316478
  T18_C_Starved   7.04705828   3.447342 18.1078202 -2.3948928  -0.2606642
                        PC6        PC7        PC8         PC9         PC10
  T0_-            2.2636988 -1.6046747 -1.6375648  2.10273388 1.125142e-14
  T3_A           -0.8727215 -1.1261056  8.0739865 -5.48778331 8.735200e-15
  T3_B           -2.4099736 -0.1187338  1.5396611  7.89418460 7.202572e-15
  T3_C           -3.0574301  0.9170718 -7.7374868 -4.59650194 1.505393e-14
  T18_A          -3.2147517  1.4402947  0.4969632 -0.46250920 1.333113e-14
  T18_B           5.1641489 -7.7415508 -1.2311043 -0.15208017 1.286818e-14
  T18_C          -1.5394861 14.0452857  0.2371402  0.29894270 1.474515e-14
  T18_A_Starved   9.4109694  1.4096994  0.4248301 -0.07558311 1.113692e-14
  T18_B_Starved -12.7069309 -5.5116097  0.1887208  0.37825982 1.714601e-14
  T18_C_Starved   6.9624769 -1.7096770 -0.3551460  0.10033674 1.421259e-14
                      color shape
  T0_-                   T0    T0
  T3_A                   T3    T3
  T3_B                   T3    T3
  T3_C                   T3    T3
  T18_A                 T18   T18
  T18_B                 T18   T18
  T18_C                 T18   T18
  T18_A_Starved T18_Starved   T18
  T18_B_Starved T18_Starved   T18
  T18_C_Starved T18_Starved   T18
Code
  p$layers
Output
  [[1]]
  geom_point: na.rm = FALSE
  stat_identity: na.rm = FALSE
  position_identity


cancerit/RCRISPR documentation built on April 26, 2023, 10:12 p.m.