StatsAnalysis_7. 主成分分析、典型相关分析R代码.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A matrix: 20 × 6 of type dbl</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>X1</th><th scope=col>X2</th><th scope=col>X3</th><th scope=col>Y1</th><th scope=col>Y2</th><th scope=col>Y3</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td> 0.5022173</td><td> 0.1873845</td><td>-0.84600343</td><td>-0.8418021</td><td> 0.2629199</td><td>-0.2008679</td></tr>\n",
       "\t<tr><td> 0.5832201</td><td> 0.8119993</td><td> 0.26350927</td><td> 0.4823810</td><td>-0.7120415</td><td> 0.5987035</td></tr>\n",
       "\t<tr><td> 0.4212145</td><td>-0.1249230</td><td>-1.40075978</td><td> 0.6715500</td><td> 0.1510391</td><td>-0.2398714</td></tr>\n",
       "\t<tr><td> 1.3122453</td><td> 0.8119993</td><td>-0.01386891</td><td>-0.2742951</td><td>-0.7120415</td><td>-0.6299063</td></tr>\n",
       "\t<tr><td>-0.1053036</td><td>-1.3741527</td><td> 2.48253466</td><td> 1.0498880</td><td> 0.8702730</td><td>-0.5909028</td></tr>\n",
       "\t<tr><td>-0.3888134</td><td>-0.4372304</td><td>-0.84600343</td><td> 1.4282260</td><td>-0.4083650</td><td>-0.6299063</td></tr>\n",
       "\t<tr><td>-0.9963344</td><td>-0.4372304</td><td> 1.09564379</td><td> 0.8607190</td><td> 1.1100176</td><td> 0.6767105</td></tr>\n",
       "\t<tr><td> 0.5832201</td><td> 0.1873845</td><td>-1.40075978</td><td>-0.6526331</td><td>-1.2075138</td><td>-0.7664185</td></tr>\n",
       "\t<tr><td>-0.1053036</td><td> 0.4996919</td><td>-0.29124708</td><td>-1.0309711</td><td>-1.3673435</td><td>-0.8834289</td></tr>\n",
       "\t<tr><td>-0.9153316</td><td>-0.7495378</td><td>-0.29124708</td><td> 1.0498880</td><td> 1.2698474</td><td> 0.0526547</td></tr>\n",
       "\t<tr><td> 0.4212145</td><td> 0.4996919</td><td>-0.56862526</td><td>-1.4093091</td><td>-0.5681948</td><td>-0.2008679</td></tr>\n",
       "\t<tr><td>-0.6723232</td><td>-0.1249230</td><td> 0.81826561</td><td> 0.4823810</td><td>-0.6481096</td><td>-0.6494080</td></tr>\n",
       "\t<tr><td> 0.1377048</td><td> 0.1873845</td><td>-0.01386891</td><td>-1.0309711</td><td>-0.7120415</td><td>-0.5518993</td></tr>\n",
       "\t<tr><td>-0.4698162</td><td>-0.4372304</td><td> 0.54088744</td><td>-0.6526331</td><td>-0.3284501</td><td>-0.5909028</td></tr>\n",
       "\t<tr><td>-0.9963344</td><td>-0.7495378</td><td>-0.01386891</td><td> 1.4282260</td><td> 1.6854047</td><td> 3.5044631</td></tr>\n",
       "\t<tr><td>-0.5103176</td><td>-0.7495378</td><td>-0.56862526</td><td> 0.6715500</td><td> 1.0301027</td><td> 0.8717279</td></tr>\n",
       "\t<tr><td> 2.7702957</td><td> 3.3104588</td><td>-0.84600343</td><td>-1.5984781</td><td>-1.5271733</td><td>-0.3958854</td></tr>\n",
       "\t<tr><td> 0.9477327</td><td> 0.4996919</td><td> 0.81826561</td><td> 0.4823810</td><td> 1.0301027</td><td> 0.9692366</td></tr>\n",
       "\t<tr><td>-0.8748302</td><td>-1.0618453</td><td>-0.56862526</td><td> 0.2932120</td><td> 1.3497622</td><td> 0.1891669</td></tr>\n",
       "\t<tr><td>-1.6443568</td><td>-0.7495378</td><td> 1.65040014</td><td>-1.4093091</td><td>-0.5681948</td><td>-0.5323976</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A matrix: 20 × 6 of type dbl\n",
       "\\begin{tabular}{llllll}\n",
       " X1 & X2 & X3 & Y1 & Y2 & Y3\\\\\n",
       "\\hline\n",
       "\t  0.5022173 &  0.1873845 & -0.84600343 & -0.8418021 &  0.2629199 & -0.2008679\\\\\n",
       "\t  0.5832201 &  0.8119993 &  0.26350927 &  0.4823810 & -0.7120415 &  0.5987035\\\\\n",
       "\t  0.4212145 & -0.1249230 & -1.40075978 &  0.6715500 &  0.1510391 & -0.2398714\\\\\n",
       "\t  1.3122453 &  0.8119993 & -0.01386891 & -0.2742951 & -0.7120415 & -0.6299063\\\\\n",
       "\t -0.1053036 & -1.3741527 &  2.48253466 &  1.0498880 &  0.8702730 & -0.5909028\\\\\n",
       "\t -0.3888134 & -0.4372304 & -0.84600343 &  1.4282260 & -0.4083650 & -0.6299063\\\\\n",
       "\t -0.9963344 & -0.4372304 &  1.09564379 &  0.8607190 &  1.1100176 &  0.6767105\\\\\n",
       "\t  0.5832201 &  0.1873845 & -1.40075978 & -0.6526331 & -1.2075138 & -0.7664185\\\\\n",
       "\t -0.1053036 &  0.4996919 & -0.29124708 & -1.0309711 & -1.3673435 & -0.8834289\\\\\n",
       "\t -0.9153316 & -0.7495378 & -0.29124708 &  1.0498880 &  1.2698474 &  0.0526547\\\\\n",
       "\t  0.4212145 &  0.4996919 & -0.56862526 & -1.4093091 & -0.5681948 & -0.2008679\\\\\n",
       "\t -0.6723232 & -0.1249230 &  0.81826561 &  0.4823810 & -0.6481096 & -0.6494080\\\\\n",
       "\t  0.1377048 &  0.1873845 & -0.01386891 & -1.0309711 & -0.7120415 & -0.5518993\\\\\n",
       "\t -0.4698162 & -0.4372304 &  0.54088744 & -0.6526331 & -0.3284501 & -0.5909028\\\\\n",
       "\t -0.9963344 & -0.7495378 & -0.01386891 &  1.4282260 &  1.6854047 &  3.5044631\\\\\n",
       "\t -0.5103176 & -0.7495378 & -0.56862526 &  0.6715500 &  1.0301027 &  0.8717279\\\\\n",
       "\t  2.7702957 &  3.3104588 & -0.84600343 & -1.5984781 & -1.5271733 & -0.3958854\\\\\n",
       "\t  0.9477327 &  0.4996919 &  0.81826561 &  0.4823810 &  1.0301027 &  0.9692366\\\\\n",
       "\t -0.8748302 & -1.0618453 & -0.56862526 &  0.2932120 &  1.3497622 &  0.1891669\\\\\n",
       "\t -1.6443568 & -0.7495378 &  1.65040014 & -1.4093091 & -0.5681948 & -0.5323976\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A matrix: 20 × 6 of type dbl\n",
       "\n",
       "| X1 | X2 | X3 | Y1 | Y2 | Y3 |\n",
       "|---|---|---|---|---|---|\n",
       "|  0.5022173 |  0.1873845 | -0.84600343 | -0.8418021 |  0.2629199 | -0.2008679 |\n",
       "|  0.5832201 |  0.8119993 |  0.26350927 |  0.4823810 | -0.7120415 |  0.5987035 |\n",
       "|  0.4212145 | -0.1249230 | -1.40075978 |  0.6715500 |  0.1510391 | -0.2398714 |\n",
       "|  1.3122453 |  0.8119993 | -0.01386891 | -0.2742951 | -0.7120415 | -0.6299063 |\n",
       "| -0.1053036 | -1.3741527 |  2.48253466 |  1.0498880 |  0.8702730 | -0.5909028 |\n",
       "| -0.3888134 | -0.4372304 | -0.84600343 |  1.4282260 | -0.4083650 | -0.6299063 |\n",
       "| -0.9963344 | -0.4372304 |  1.09564379 |  0.8607190 |  1.1100176 |  0.6767105 |\n",
       "|  0.5832201 |  0.1873845 | -1.40075978 | -0.6526331 | -1.2075138 | -0.7664185 |\n",
       "| -0.1053036 |  0.4996919 | -0.29124708 | -1.0309711 | -1.3673435 | -0.8834289 |\n",
       "| -0.9153316 | -0.7495378 | -0.29124708 |  1.0498880 |  1.2698474 |  0.0526547 |\n",
       "|  0.4212145 |  0.4996919 | -0.56862526 | -1.4093091 | -0.5681948 | -0.2008679 |\n",
       "| -0.6723232 | -0.1249230 |  0.81826561 |  0.4823810 | -0.6481096 | -0.6494080 |\n",
       "|  0.1377048 |  0.1873845 | -0.01386891 | -1.0309711 | -0.7120415 | -0.5518993 |\n",
       "| -0.4698162 | -0.4372304 |  0.54088744 | -0.6526331 | -0.3284501 | -0.5909028 |\n",
       "| -0.9963344 | -0.7495378 | -0.01386891 |  1.4282260 |  1.6854047 |  3.5044631 |\n",
       "| -0.5103176 | -0.7495378 | -0.56862526 |  0.6715500 |  1.0301027 |  0.8717279 |\n",
       "|  2.7702957 |  3.3104588 | -0.84600343 | -1.5984781 | -1.5271733 | -0.3958854 |\n",
       "|  0.9477327 |  0.4996919 |  0.81826561 |  0.4823810 |  1.0301027 |  0.9692366 |\n",
       "| -0.8748302 | -1.0618453 | -0.56862526 |  0.2932120 |  1.3497622 |  0.1891669 |\n",
       "| -1.6443568 | -0.7495378 |  1.65040014 | -1.4093091 | -0.5681948 | -0.5323976 |\n",
       "\n"
      ],
      "text/plain": [
       "      X1         X2         X3          Y1         Y2         Y3        \n",
       " [1,]  0.5022173  0.1873845 -0.84600343 -0.8418021  0.2629199 -0.2008679\n",
       " [2,]  0.5832201  0.8119993  0.26350927  0.4823810 -0.7120415  0.5987035\n",
       " [3,]  0.4212145 -0.1249230 -1.40075978  0.6715500  0.1510391 -0.2398714\n",
       " [4,]  1.3122453  0.8119993 -0.01386891 -0.2742951 -0.7120415 -0.6299063\n",
       " [5,] -0.1053036 -1.3741527  2.48253466  1.0498880  0.8702730 -0.5909028\n",
       " [6,] -0.3888134 -0.4372304 -0.84600343  1.4282260 -0.4083650 -0.6299063\n",
       " [7,] -0.9963344 -0.4372304  1.09564379  0.8607190  1.1100176  0.6767105\n",
       " [8,]  0.5832201  0.1873845 -1.40075978 -0.6526331 -1.2075138 -0.7664185\n",
       " [9,] -0.1053036  0.4996919 -0.29124708 -1.0309711 -1.3673435 -0.8834289\n",
       "[10,] -0.9153316 -0.7495378 -0.29124708  1.0498880  1.2698474  0.0526547\n",
       "[11,]  0.4212145  0.4996919 -0.56862526 -1.4093091 -0.5681948 -0.2008679\n",
       "[12,] -0.6723232 -0.1249230  0.81826561  0.4823810 -0.6481096 -0.6494080\n",
       "[13,]  0.1377048  0.1873845 -0.01386891 -1.0309711 -0.7120415 -0.5518993\n",
       "[14,] -0.4698162 -0.4372304  0.54088744 -0.6526331 -0.3284501 -0.5909028\n",
       "[15,] -0.9963344 -0.7495378 -0.01386891  1.4282260  1.6854047  3.5044631\n",
       "[16,] -0.5103176 -0.7495378 -0.56862526  0.6715500  1.0301027  0.8717279\n",
       "[17,]  2.7702957  3.3104588 -0.84600343 -1.5984781 -1.5271733 -0.3958854\n",
       "[18,]  0.9477327  0.4996919  0.81826561  0.4823810  1.0301027  0.9692366\n",
       "[19,] -0.8748302 -1.0618453 -0.56862526  0.2932120  1.3497622  0.1891669\n",
       "[20,] -1.6443568 -0.7495378  1.65040014 -1.4093091 -0.5681948 -0.5323976"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<dl>\n",
       "\t<dt>$cor</dt>\n",
       "\t\t<dd><style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>0.795608154419992</li><li>0.200556041107123</li><li>0.0725702862103672</li></ol>\n",
       "</dd>\n",
       "\t<dt>$xcoef</dt>\n",
       "\t\t<dd><table>\n",
       "<caption>A matrix: 3 × 3 of type dbl</caption>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>X1</th><td>-0.17788841</td><td>-0.43230348</td><td>-0.04381432</td></tr>\n",
       "\t<tr><th scope=row>X2</th><td> 0.36232695</td><td> 0.27085764</td><td> 0.11608883</td></tr>\n",
       "\t<tr><th scope=row>X3</th><td>-0.01356309</td><td>-0.05301954</td><td> 0.24106633</td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "</dd>\n",
       "\t<dt>$ycoef</dt>\n",
       "\t\t<dd><table>\n",
       "<caption>A matrix: 3 × 3 of type dbl</caption>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>Y1</th><td>-0.08018009</td><td>-0.08615561</td><td>-0.29745900</td></tr>\n",
       "\t<tr><th scope=row>Y2</th><td>-0.24180670</td><td> 0.02833066</td><td> 0.28373986</td></tr>\n",
       "\t<tr><th scope=row>Y3</th><td> 0.16435956</td><td> 0.24367781</td><td>-0.09608099</td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "</dd>\n",
       "\t<dt>$xcenter</dt>\n",
       "\t\t<dd><style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>X1</dt><dd>2.28983498828939e-16</dd><dt>X2</dt><dd>4.3159920082303e-16</dd><dt>X3</dt><dd>-1.77895892461422e-16</dd></dl>\n",
       "</dd>\n",
       "\t<dt>$ycenter</dt>\n",
       "\t\t<dd><style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Y1</dt><dd>1.47104550762833e-16</dd><dt>Y2</dt><dd>-1.77635683940025e-16</dd><dt>Y3</dt><dd>4.9960036108132e-17</dd></dl>\n",
       "</dd>\n",
       "</dl>\n"
      ],
      "text/latex": [
       "\\begin{description}\n",
       "\\item[\\$cor] \\begin{enumerate*}\n",
       "\\item 0.795608154419992\n",
       "\\item 0.200556041107123\n",
       "\\item 0.0725702862103672\n",
       "\\end{enumerate*}\n",
       "\n",
       "\\item[\\$xcoef] A matrix: 3 × 3 of type dbl\n",
       "\\begin{tabular}{r|lll}\n",
       "\tX1 & -0.17788841 & -0.43230348 & -0.04381432\\\\\n",
       "\tX2 &  0.36232695 &  0.27085764 &  0.11608883\\\\\n",
       "\tX3 & -0.01356309 & -0.05301954 &  0.24106633\\\\\n",
       "\\end{tabular}\n",
       "\n",
       "\\item[\\$ycoef] A matrix: 3 × 3 of type dbl\n",
       "\\begin{tabular}{r|lll}\n",
       "\tY1 & -0.08018009 & -0.08615561 & -0.29745900\\\\\n",
       "\tY2 & -0.24180670 &  0.02833066 &  0.28373986\\\\\n",
       "\tY3 &  0.16435956 &  0.24367781 & -0.09608099\\\\\n",
       "\\end{tabular}\n",
       "\n",
       "\\item[\\$xcenter] \\begin{description*}\n",
       "\\item[X1] 2.28983498828939e-16\n",
       "\\item[X2] 4.3159920082303e-16\n",
       "\\item[X3] -1.77895892461422e-16\n",
       "\\end{description*}\n",
       "\n",
       "\\item[\\$ycenter] \\begin{description*}\n",
       "\\item[Y1] 1.47104550762833e-16\n",
       "\\item[Y2] -1.77635683940025e-16\n",
       "\\item[Y3] 4.9960036108132e-17\n",
       "\\end{description*}\n",
       "\n",
       "\\end{description}\n"
      ],
      "text/markdown": [
       "$cor\n",
       ":   1. 0.795608154419992\n",
       "2. 0.200556041107123\n",
       "3. 0.0725702862103672\n",
       "\n",
       "\n",
       "\n",
       "$xcoef\n",
       ":   \n",
       "A matrix: 3 × 3 of type dbl\n",
       "\n",
       "| X1 | -0.17788841 | -0.43230348 | -0.04381432 |\n",
       "| X2 |  0.36232695 |  0.27085764 |  0.11608883 |\n",
       "| X3 | -0.01356309 | -0.05301954 |  0.24106633 |\n",
       "\n",
       "\n",
       "$ycoef\n",
       ":   \n",
       "A matrix: 3 × 3 of type dbl\n",
       "\n",
       "| Y1 | -0.08018009 | -0.08615561 | -0.29745900 |\n",
       "| Y2 | -0.24180670 |  0.02833066 |  0.28373986 |\n",
       "| Y3 |  0.16435956 |  0.24367781 | -0.09608099 |\n",
       "\n",
       "\n",
       "$xcenter\n",
       ":   X1\n",
       ":   2.28983498828939e-16X2\n",
       ":   4.3159920082303e-16X3\n",
       ":   -1.77895892461422e-16\n",
       "\n",
       "\n",
       "$ycenter\n",
       ":   Y1\n",
       ":   1.47104550762833e-16Y2\n",
       ":   -1.77635683940025e-16Y3\n",
       ":   4.9960036108132e-17\n",
       "\n",
       "\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "$cor\n",
       "[1] 0.79560815 0.20055604 0.07257029\n",
       "\n",
       "$xcoef\n",
       "          [,1]        [,2]        [,3]\n",
       "X1 -0.17788841 -0.43230348 -0.04381432\n",
       "X2  0.36232695  0.27085764  0.11608883\n",
       "X3 -0.01356309 -0.05301954  0.24106633\n",
       "\n",
       "$ycoef\n",
       "          [,1]        [,2]        [,3]\n",
       "Y1 -0.08018009 -0.08615561 -0.29745900\n",
       "Y2 -0.24180670  0.02833066  0.28373986\n",
       "Y3  0.16435956  0.24367781 -0.09608099\n",
       "\n",
       "$xcenter\n",
       "           X1            X2            X3 \n",
       " 2.289835e-16  4.315992e-16 -1.778959e-16 \n",
       "\n",
       "$ycenter\n",
       "           Y1            Y2            Y3 \n",
       " 1.471046e-16 -1.776357e-16  4.996004e-17 \n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A matrix: 20 × 3 of type dbl</caption>\n",
       "<tbody>\n",
       "\t<tr><td>-0.009969788</td><td>-0.121501078</td><td>-0.20419401</td></tr>\n",
       "\t<tr><td> 0.186887139</td><td>-0.046163013</td><td> 0.13223387</td></tr>\n",
       "\t<tr><td>-0.101193522</td><td>-0.141661215</td><td>-0.37063341</td></tr>\n",
       "\t<tr><td> 0.060964112</td><td>-0.346616669</td><td> 0.03342558</td></tr>\n",
       "\t<tr><td>-0.512831098</td><td>-0.458299483</td><td> 0.44354554</td></tr>\n",
       "\t<tr><td>-0.077780541</td><td> 0.094512914</td><td>-0.23766491</td></tr>\n",
       "\t<tr><td> 0.003955674</td><td> 0.254201102</td><td> 0.25701898</td></tr>\n",
       "\t<tr><td>-0.016855040</td><td>-0.127105942</td><td>-0.34147617</td></tr>\n",
       "\t<tr><td> 0.203734347</td><td> 0.196310283</td><td>-0.00758741</td></tr>\n",
       "\t<tr><td>-0.104800666</td><td> 0.208124774</td><td>-0.11711820</td></tr>\n",
       "\t<tr><td> 0.113834968</td><td>-0.016598895</td><td>-0.09752299</td></tr>\n",
       "\t<tr><td> 0.063237343</td><td> 0.213427257</td><td> 0.21221151</td></tr>\n",
       "\t<tr><td> 0.043586465</td><td>-0.008040409</td><td> 0.01237648</td></tr>\n",
       "\t<tr><td>-0.082181602</td><td> 0.055998387</td><td> 0.10021686</td></tr>\n",
       "\t<tr><td>-0.094153311</td><td> 0.228436101</td><td>-0.04670258</td></tr>\n",
       "\t<tr><td>-0.173085857</td><td> 0.047742282</td><td>-0.20173015</td></tr>\n",
       "\t<tr><td> 0.718139369</td><td>-0.256090676</td><td> 0.05898572</td></tr>\n",
       "\t<tr><td> 0.001362964</td><td>-0.317746855</td><td> 0.21374067</td></tr>\n",
       "\t<tr><td>-0.221400693</td><td> 0.120731486</td><td>-0.22201469</td></tr>\n",
       "\t<tr><td>-0.001450263</td><td> 0.420339649</td><td> 0.38288931</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A matrix: 20 × 3 of type dbl\n",
       "\\begin{tabular}{lll}\n",
       "\t -0.009969788 & -0.121501078 & -0.20419401\\\\\n",
       "\t  0.186887139 & -0.046163013 &  0.13223387\\\\\n",
       "\t -0.101193522 & -0.141661215 & -0.37063341\\\\\n",
       "\t  0.060964112 & -0.346616669 &  0.03342558\\\\\n",
       "\t -0.512831098 & -0.458299483 &  0.44354554\\\\\n",
       "\t -0.077780541 &  0.094512914 & -0.23766491\\\\\n",
       "\t  0.003955674 &  0.254201102 &  0.25701898\\\\\n",
       "\t -0.016855040 & -0.127105942 & -0.34147617\\\\\n",
       "\t  0.203734347 &  0.196310283 & -0.00758741\\\\\n",
       "\t -0.104800666 &  0.208124774 & -0.11711820\\\\\n",
       "\t  0.113834968 & -0.016598895 & -0.09752299\\\\\n",
       "\t  0.063237343 &  0.213427257 &  0.21221151\\\\\n",
       "\t  0.043586465 & -0.008040409 &  0.01237648\\\\\n",
       "\t -0.082181602 &  0.055998387 &  0.10021686\\\\\n",
       "\t -0.094153311 &  0.228436101 & -0.04670258\\\\\n",
       "\t -0.173085857 &  0.047742282 & -0.20173015\\\\\n",
       "\t  0.718139369 & -0.256090676 &  0.05898572\\\\\n",
       "\t  0.001362964 & -0.317746855 &  0.21374067\\\\\n",
       "\t -0.221400693 &  0.120731486 & -0.22201469\\\\\n",
       "\t -0.001450263 &  0.420339649 &  0.38288931\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A matrix: 20 × 3 of type dbl\n",
       "\n",
       "| -0.009969788 | -0.121501078 | -0.20419401 |\n",
       "|  0.186887139 | -0.046163013 |  0.13223387 |\n",
       "| -0.101193522 | -0.141661215 | -0.37063341 |\n",
       "|  0.060964112 | -0.346616669 |  0.03342558 |\n",
       "| -0.512831098 | -0.458299483 |  0.44354554 |\n",
       "| -0.077780541 |  0.094512914 | -0.23766491 |\n",
       "|  0.003955674 |  0.254201102 |  0.25701898 |\n",
       "| -0.016855040 | -0.127105942 | -0.34147617 |\n",
       "|  0.203734347 |  0.196310283 | -0.00758741 |\n",
       "| -0.104800666 |  0.208124774 | -0.11711820 |\n",
       "|  0.113834968 | -0.016598895 | -0.09752299 |\n",
       "|  0.063237343 |  0.213427257 |  0.21221151 |\n",
       "|  0.043586465 | -0.008040409 |  0.01237648 |\n",
       "| -0.082181602 |  0.055998387 |  0.10021686 |\n",
       "| -0.094153311 |  0.228436101 | -0.04670258 |\n",
       "| -0.173085857 |  0.047742282 | -0.20173015 |\n",
       "|  0.718139369 | -0.256090676 |  0.05898572 |\n",
       "|  0.001362964 | -0.317746855 |  0.21374067 |\n",
       "| -0.221400693 |  0.120731486 | -0.22201469 |\n",
       "| -0.001450263 |  0.420339649 |  0.38288931 |\n",
       "\n"
      ],
      "text/plain": [
       "      [,1]         [,2]         [,3]       \n",
       " [1,] -0.009969788 -0.121501078 -0.20419401\n",
       " [2,]  0.186887139 -0.046163013  0.13223387\n",
       " [3,] -0.101193522 -0.141661215 -0.37063341\n",
       " [4,]  0.060964112 -0.346616669  0.03342558\n",
       " [5,] -0.512831098 -0.458299483  0.44354554\n",
       " [6,] -0.077780541  0.094512914 -0.23766491\n",
       " [7,]  0.003955674  0.254201102  0.25701898\n",
       " [8,] -0.016855040 -0.127105942 -0.34147617\n",
       " [9,]  0.203734347  0.196310283 -0.00758741\n",
       "[10,] -0.104800666  0.208124774 -0.11711820\n",
       "[11,]  0.113834968 -0.016598895 -0.09752299\n",
       "[12,]  0.063237343  0.213427257  0.21221151\n",
       "[13,]  0.043586465 -0.008040409  0.01237648\n",
       "[14,] -0.082181602  0.055998387  0.10021686\n",
       "[15,] -0.094153311  0.228436101 -0.04670258\n",
       "[16,] -0.173085857  0.047742282 -0.20173015\n",
       "[17,]  0.718139369 -0.256090676  0.05898572\n",
       "[18,]  0.001362964 -0.317746855  0.21374067\n",
       "[19,] -0.221400693  0.120731486 -0.22201469\n",
       "[20,] -0.001450263  0.420339649  0.38288931"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A matrix: 20 × 3 of type dbl</caption>\n",
       "<tbody>\n",
       "\t<tr><td>-0.02909460</td><td> 0.031027608</td><td> 0.344302062</td></tr>\n",
       "\t<tr><td> 0.23190170</td><td> 0.084158321</td><td>-0.403047146</td></tr>\n",
       "\t<tr><td>-0.12979237</td><td>-0.112030106</td><td>-0.133855684</td></tr>\n",
       "\t<tr><td> 0.09063830</td><td>-0.150034732</td><td>-0.059921010</td></tr>\n",
       "\t<tr><td>-0.39173848</td><td>-0.209788233</td><td>-0.008592976</td></tr>\n",
       "\t<tr><td>-0.11930102</td><td>-0.288113119</td><td>-0.480186091</td></tr>\n",
       "\t<tr><td>-0.22619839</td><td> 0.122191086</td><td>-0.006091381</td></tr>\n",
       "\t<tr><td> 0.21834490</td><td>-0.164740837</td><td>-0.074849953</td></tr>\n",
       "\t<tr><td> 0.26809619</td><td>-0.165185828</td><td> 0.003582504</td></tr>\n",
       "\t<tr><td>-0.38258341</td><td>-0.041647344</td><td> 0.042948559</td></tr>\n",
       "\t<tr><td> 0.21737727</td><td> 0.056375494</td><td> 0.277291769</td></tr>\n",
       "\t<tr><td> 0.01130349</td><td>-0.218167527</td><td>-0.264987331</td></tr>\n",
       "\t<tr><td> 0.16412985</td><td>-0.065834278</td><td> 0.157664098</td></tr>\n",
       "\t<tr><td> 0.03462910</td><td>-0.097067107</td><td> 0.157711712</td></tr>\n",
       "\t<tr><td> 0.05393456</td><td> 0.778658842</td><td>-0.283334468</td></tr>\n",
       "\t<tr><td>-0.15965387</td><td> 0.183746435</td><td> 0.008766141</td></tr>\n",
       "\t<tr><td> 0.43237930</td><td>-0.002016448</td><td> 0.080198839</td></tr>\n",
       "\t<tr><td>-0.12845980</td><td> 0.223805116</td><td> 0.055667431</td></tr>\n",
       "\t<tr><td>-0.31879992</td><td> 0.059073577</td><td> 0.277587462</td></tr>\n",
       "\t<tr><td> 0.16288721</td><td>-0.024410919</td><td> 0.309145462</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A matrix: 20 × 3 of type dbl\n",
       "\\begin{tabular}{lll}\n",
       "\t -0.02909460 &  0.031027608 &  0.344302062\\\\\n",
       "\t  0.23190170 &  0.084158321 & -0.403047146\\\\\n",
       "\t -0.12979237 & -0.112030106 & -0.133855684\\\\\n",
       "\t  0.09063830 & -0.150034732 & -0.059921010\\\\\n",
       "\t -0.39173848 & -0.209788233 & -0.008592976\\\\\n",
       "\t -0.11930102 & -0.288113119 & -0.480186091\\\\\n",
       "\t -0.22619839 &  0.122191086 & -0.006091381\\\\\n",
       "\t  0.21834490 & -0.164740837 & -0.074849953\\\\\n",
       "\t  0.26809619 & -0.165185828 &  0.003582504\\\\\n",
       "\t -0.38258341 & -0.041647344 &  0.042948559\\\\\n",
       "\t  0.21737727 &  0.056375494 &  0.277291769\\\\\n",
       "\t  0.01130349 & -0.218167527 & -0.264987331\\\\\n",
       "\t  0.16412985 & -0.065834278 &  0.157664098\\\\\n",
       "\t  0.03462910 & -0.097067107 &  0.157711712\\\\\n",
       "\t  0.05393456 &  0.778658842 & -0.283334468\\\\\n",
       "\t -0.15965387 &  0.183746435 &  0.008766141\\\\\n",
       "\t  0.43237930 & -0.002016448 &  0.080198839\\\\\n",
       "\t -0.12845980 &  0.223805116 &  0.055667431\\\\\n",
       "\t -0.31879992 &  0.059073577 &  0.277587462\\\\\n",
       "\t  0.16288721 & -0.024410919 &  0.309145462\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A matrix: 20 × 3 of type dbl\n",
       "\n",
       "| -0.02909460 |  0.031027608 |  0.344302062 |\n",
       "|  0.23190170 |  0.084158321 | -0.403047146 |\n",
       "| -0.12979237 | -0.112030106 | -0.133855684 |\n",
       "|  0.09063830 | -0.150034732 | -0.059921010 |\n",
       "| -0.39173848 | -0.209788233 | -0.008592976 |\n",
       "| -0.11930102 | -0.288113119 | -0.480186091 |\n",
       "| -0.22619839 |  0.122191086 | -0.006091381 |\n",
       "|  0.21834490 | -0.164740837 | -0.074849953 |\n",
       "|  0.26809619 | -0.165185828 |  0.003582504 |\n",
       "| -0.38258341 | -0.041647344 |  0.042948559 |\n",
       "|  0.21737727 |  0.056375494 |  0.277291769 |\n",
       "|  0.01130349 | -0.218167527 | -0.264987331 |\n",
       "|  0.16412985 | -0.065834278 |  0.157664098 |\n",
       "|  0.03462910 | -0.097067107 |  0.157711712 |\n",
       "|  0.05393456 |  0.778658842 | -0.283334468 |\n",
       "| -0.15965387 |  0.183746435 |  0.008766141 |\n",
       "|  0.43237930 | -0.002016448 |  0.080198839 |\n",
       "| -0.12845980 |  0.223805116 |  0.055667431 |\n",
       "| -0.31879992 |  0.059073577 |  0.277587462 |\n",
       "|  0.16288721 | -0.024410919 |  0.309145462 |\n",
       "\n"
      ],
      "text/plain": [
       "      [,1]        [,2]         [,3]        \n",
       " [1,] -0.02909460  0.031027608  0.344302062\n",
       " [2,]  0.23190170  0.084158321 -0.403047146\n",
       " [3,] -0.12979237 -0.112030106 -0.133855684\n",
       " [4,]  0.09063830 -0.150034732 -0.059921010\n",
       " [5,] -0.39173848 -0.209788233 -0.008592976\n",
       " [6,] -0.11930102 -0.288113119 -0.480186091\n",
       " [7,] -0.22619839  0.122191086 -0.006091381\n",
       " [8,]  0.21834490 -0.164740837 -0.074849953\n",
       " [9,]  0.26809619 -0.165185828  0.003582504\n",
       "[10,] -0.38258341 -0.041647344  0.042948559\n",
       "[11,]  0.21737727  0.056375494  0.277291769\n",
       "[12,]  0.01130349 -0.218167527 -0.264987331\n",
       "[13,]  0.16412985 -0.065834278  0.157664098\n",
       "[14,]  0.03462910 -0.097067107  0.157711712\n",
       "[15,]  0.05393456  0.778658842 -0.283334468\n",
       "[16,] -0.15965387  0.183746435  0.008766141\n",
       "[17,]  0.43237930 -0.002016448  0.080198839\n",
       "[18,] -0.12845980  0.223805116  0.055667431\n",
       "[19,] -0.31879992  0.059073577  0.277587462\n",
       "[20,]  0.16288721 -0.024410919  0.309145462"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAMFBMVEUAAABNTU1oaGh8fHyM\njIyampqnp6eysrK9vb3Hx8fQ0NDZ2dnh4eHp6enw8PD////QFLu4AAAACXBIWXMAABJ0AAAS\ndAHeZh94AAAS6klEQVR4nO3d4ULaShhF0QkgKgJ9/7e9NdZeFCiKZ5JMstYPRdv6NcBWMiFY\nfgE/Vsb+D8AcCAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBAS\nBAgJAoQEAUKCACFBgJAgQEgQUD2kAjMyXki1B8BwhAQBQoIAIUGAkCBASBAgJAgQEgQICQKE\nBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJDgq+55jRMhwQd9RddSEhJ8TTl5e+UP\nv/MHKUKiKeXT+8t/+vU/SBESTRESBAgJEuwjQYBVO4hwHAnqEhIECAkChAQBQoIAIUGAkCBA\nSBBQLaTrh65CA2BChAQB4ZDKR3dNhvaEQ3rphMQSpR/aHTdlfej/vYd2LEh+H+m5lOdfQmJZ\nKiw2HNZlcxQSi1Jl1e6xdDshsSR1lr/3qxsrDT8dANNS6zjSg5BYEk8RggAhQUDNkByQZTGG\nDenLT3uAtnhoBwFCggAhQUA8pJfHTb8HtNm+3DkZ2hMO6bg6WU1Y3zcZ2hMOaVu6531/6bDr\nyvauydCecEhd2f+9vC/dXZOhPfFTza99EBoAk+QnEgTk95F2/Znm9pFYlPTy9/pk1W51vGsy\ntCd/HGnbH0fqNo+OI7EcntkAAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAg\nJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkC\nhAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAk\nCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEheU4ub5\nHiFxpq9ISt8iJM6Uk7d8jZD4rHx6zxcIic+EdAch8ZmQ7iAkzthH+j4hccaq3fcJiQscR/ou\nIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICT+56Ua7iYk3nnxoB8QEu+8nN0PCIk/vMDqTwiJP4T0E0LiDyH9hJB4Zx/pB4TEO6t2PyAk\n/uc40t2EBAFCggAhTVvowZbHbLUJacpCu/9WEeoT0pSFFqSta9cnpAkLHSJ1pHUAQpowIbUj\nHdLxoZT17s+X+OfXcLPeJKR2hEM6duXV5u1LCOmH7CM1IxzStjz9rumpW/dfQkg/ZNWuGeGQ\nurd/duhWByElOI7UiHBI77fXcb0WEgsSDmlVju+X1kJiOcIhPZWHP5cOZS0kFiO9/L39W8/u\nxuNyITEj8QOy+837pcODkFgKz2yAACFBgJAgoGZIFhtYjGFDKqd+PgCmwkM7CBASBAgJAuIh\nvTxu3k5J2r7cORnakz6xb3WymrC+bzK0J35iX/e87y8ddl3Z3jUZ2hM/sW//9/K+dHdNhvZU\nOrHv/IPQAJgkP5EgIL+PtDv0l+wjsSTp5e/1yard6vivvykkZiR/HGnbH0fqNo+OI7EcntkA\nAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkBrnBQKnQUhtuNKL3w47FUJqwdVeyslb\nxiSkFlzrpXx6z2iE1ICrvQhpMoTUACFNn5AacL0X+0hTIaQWXO3Fqt1UCKkF/+jFcaRpEFIb\n9DJxQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAk\nCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQpqLUlyjIxLSPPQVSWk8QpqHcvKW\nEQhpFsqn9wxNSLMgpLEJaRaENDYhzYN9pJEJaR6s2o1MSHPhONKohAQBQoIAIUGAkCBASBAg\nJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJDG4PUVZkdIw/OKPzMk\npOF5DboZEtLgvCrqHAlpcEKaIyENTkhzJKTh2UeaISENz6rdDAlpDI4jzY6QIEBIECAkCBAS\nBAgJAoQEAUIaj0XwGRHSWByWnRUhjcUThWZFSCPx1NV5EdJIhDQvQhqJkOZFSGOxjzQrN0Iq\nHw0yeSGs2s3KjZCehFSP40gzcuuh3b5bDz0Z2nNzH2lftgNPhvbcXmx4KvthJ0N7rNpBwI2Q\ndiNMhvbcWv7utoehJ0N7boS0KqWs6/xYEhIzcmsf6bDtfre0rbDeICRm5AuLDS8Pv1NaPR2H\nmgzt+dqq3fP6d0sPX3qI9/K46Z8Esdm+3DkZ2vPV5e/j4+/dpe7mlzuuTp5Q9O/nRAgpzDOO\nxvSN40i7L9xQ29I9v+1PHXbdv58T4VaP8hzYcYV/InUnT4PY//vvu82jnJUxrvA+0odvif/+\n/ug2T3Ke4Mi+ENLuG6t2fiKNREgjuxXSy+txpO7Lx5F+7yPt3p4KYR9pUEIaWfqZDeuTVbvV\nP3+Iucmj7CON6+Zz7R6/eSD2ZdsfR+o2j44jDcmq3bhuhHQjhiqTuY/jSGNyPhIECAkCaobk\nOBKLMWxI1V7bC8bloR0ECAkChAQB8ZCc2McShUNyYh/LFA7JiX0sUzgkp1GwTOGQnNjHMvmJ\nBAH5fSQn9rFA6eVvJ/axSPnjSE7sY4E8swEChAQBQoIAIUGAkCBASBAgpOlzWn4DhDR1Xvmx\nCUKaOq9F3AQhTZxXx2+DkCZOSG0Q0sQJqQ1Cmjr7SE0Q0tRZtWuCkKbPcaQGCAkChAQBQoIA\nIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJ\nAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAh\nQYCQIEBIECCkUZXiapgHIY2or0hKsyCkEZWTt7RNSOMpn97TMCGNR0gzIqTxCGlGhDQi+0jz\nIaQRfWvVzkr5pAlpaB+C+HIdVsonTkjDujcIjwInTkjDujMI6xJTJ6RB3RuEkKZOSIMS0lwJ\naVB3B2EfaeKENKx7g7BqN3FCGtb9QTiONGlCGpogZklIECAkCBASBAgJAoQEAUKCACFBgJAg\nQEgQICQIEBIECAkChAQBQoIAIcU4P2LJhBTiFNZlE1KIF1VYNiFleJmfhRNShpAWTkgZQlo4\nIYXYR1o2IYVYtVs2IcU4jrRkQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQICAd0vGh\nlPXuz5f459cQEjMSDunYlVebty8hJJYiHNK2PP2u6alb919CSCxFOKTu7Z8dutVBSCxIOKT3\ndo7rtZBYkHBIq3J8v7QWEssRDumpPPy5dChrIbEY6eXv7d96djdOvRYSMxI/ILvfvF86PAiJ\npfDMBggQEgQICQJqhmSxgcUQEgR4aAcBQoKAYUMqp2oMgHHEQ3p53LydkrR9uXMytCd9Yt/q\n5EfO+r7J0J74iX3d876/dNh1ZXvXZGhP/MS+/d/L+9LdNRnaU+nEvvMPQgNgkvxEgoD8PtLu\n0F+yj8SSpJe/1yerdqvjv/6mkJiR/HGkbX8cqds8Oo7EcniKEAQICQJqhXT7qXRCYkaEBAFC\nggAhQYCQmuEUrikTUiP6K1RKk2X5uxHl5C3TI6Q2lE/vmRghtUFIEyekNghp4oTUCPtI0yak\nRli1mzYhNcNxpCkTEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJ\nAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAh\nQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIENI3ldLof5yqhPQtfUVS4oyQvqWc\nvIX/Cek7yqf38IeQvkNIXCGk7xASVwjpW+wjcZmQvsWqHZcJ6ZscR+ISIUGAkCBASBAgJAgQ\nEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQB\nQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJCG\nVcoyt3v2qoV08w6zxDtUf6VIaY6ENKRy8pZZCYdUPrpr8nyVT++Zj3BIL52QrhPSfKUf2h03\nZX3o/72HdmeENF/5faTnUp5/Ceki+0izVWGx4bAum+MPQprxCrFVu9mqsmr3WLrdvSHN/L42\n4+8Sy1Zn+Xu/un2HuRbSz6fD4GodR3q4MyT74zRpak8REhJNEhIE1AzprgOy9pFo0fRCmveq\nHTM1tYd2v6wQ06IJhgTtGTakLz+jFdoSD+nlcdNnstm+3DkZ2hMO6bg6+ZGzvm8ytCcc0rZ0\nz/v+0mHXle1dk6E94ZC6sv97eV+6uyZDe+Knml/7IDQAJslPJAjI7yPt+jPN7SOxKOnl7/XJ\nqt3qeNdkaE/+ONK2P47UbR4dR2I5PEUIAoQEAbVCuv1UOiExI0KCACFBgJAgQEgQICQIsPwN\nASOGBDMyWkiVDPr/HnLYXGfNd8NGGxnhPtDYrPlu2GgjI9wHGps13w0bbWSE+0Bjs+a7YaON\njHAfaGzWfDdstJER7gONzZrvho02MsJ9oLFZ892w0UZGuA80Nmu+GzbayAj3gcZmzXfDRhsZ\n4T7Q2Kz5bthoIyPcBxqbNd8NG21khPtAY7Pmu2GjjYT5ERIECAkChAQBQoIAIUGAkCBASBAg\nJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCGgtpG1Xuu3x4+deam3E2bCn1fn0KoMubWatWbU2\n6uKwX/VurrNZ+4dSHg51hp1pLKR1/6s1Vh8+d+wqbcTZsG3/iS5+pzsbdGkza82qtVEXh/2q\nd3OdzdpV3LBzbYX0Urr9r31XXk4/ubn9awMzw/bl4fet8lQeag+6uJmVZtXaqIvDXlW6uc5n\ndb8/cdyUbY1p59oKaVt2v98+l8eTzz3/47c/ZYdt3ubEx50NurSZtWbV2qiLw37Vu7nOZj33\nCR1LV2PaubZC2pTXh7z7svn/U4eyrhTShWG9+LizQdcm15j1R5XrcMib62zWQ9nXmHNNWyGV\n82+f63KoFNKFYa+OZV170JXJVWa9yW/UlWHVbq6zWavy67HrH7cOovWQHstzne+mV+9zT/1j\niKqDhg8pv1GXh9W7uS5ciZt+saHGsEvzB5qTcXZt9T/KBw3p0MUfcI0fUoWNujis4s114Up8\nXWx4qLGjeXH+MGNCzn9+v65uDhnSscs/Bho9pBobdXFYxZvrwpX4uo90qHMQ4cL8Ycb80Ptv\nZu8+XVsP/SOS8C1zbVhvXeF2ORt0cXKlWa9qbNSlYVVuriuzan43uqStkN6WZg5/l2bKXwMM\ne/1gta5wpPxs0IXJ1WbV2qhLw6rcXFdmVV3Xv6CNkN499t/Sdn8PstW8Zc6GvV6u8hDobND5\n5Hqzam3UpWE1b64rV+Kh3tZ91FZIlw/5D/XMhlo3yqjPbKh5Txvy5rqwYavj62LDc41p59oK\n6deq/37W3/Rnj4arD3uo9f30bKtOPlF7VrWNujTs18dLdWc91rsSL2gspGP/FN/+Yv2QPg+r\n9sDkbKtOPlF7Vs0Hx4PeXOezdutaV+IFjYUE0yQkCBASBAgJAoQEAUKCACFBgJAgQEgQICQI\nEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQE\nAUKCACFBgJAgQEgQIKSG/f3dd+8XntycY3HNN+xzSPtavwSUm1zzDfsU0r4T0mhc8w37GNJT\nWQtpNK75hn0MqWyr/X53bnLNN+xjSPtfQhqPa75hZ6t2QhqNa75hQpoO13zDhDQdrvmGCWk6\nXPMNW5VD//5QVm+fENJoXPMNeyhP/fun8vD2CSGNxjXfsF3pS3ou5eXtE0IajWu+ZdvyZvvn\nYyGNxjXftN3md0ab3fuHQhqNax4ChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQ\nIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUHAf4KwXe6xg3We\nAAAAAElFTkSuQmCC",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAMFBMVEUAAABNTU1oaGh8fHyM\njIyampqnp6eysrK9vb3Hx8fQ0NDZ2dnh4eHp6enw8PD////QFLu4AAAACXBIWXMAABJ0AAAS\ndAHeZh94AAASK0lEQVR4nO3d7ULiugKG0RQQFfm4/7s9I249KAKKb9ukXeuHojNjYqeP0JRK\nOQB/VsaeAEyBkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIEBIECAkChAQBvYdUYELGC6nvAWA4QoIAIUGAkCBASBAgJAgQEgQICQKEBAFC\nggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgXXfmFFvCFkC44ViQlfkhIF5STt3CL\nkL5XvryHq4T0PSHxK0L6npD4FSFd4BiJ3xDSBVbt+A0hXeQ8Ej8nJAgQEgQICQKEBAFCggAh\nQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChMRcRX8ph5CYp/CviRIS\n8xT+xYVCYpbSv0pXSMySkCBASJDgGAkCrNpBhPNIUBshQYCQIEBIECAkCBASBAgJAoQEAUKC\nACFBgJAgQEgQICQIEBIECAkChHQieoEKsyKkD+FLJq+MU923zp8J6UP4Iv5LowyUK8MS0rv0\nr5W5Okxt3zx/JaR3w4Q0UK4MTUjvhMQfCOnDIA+6hDRRQvowzDKAY6RpEtKJIRamrdpNk5AG\n5zzSFAkJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKE\nBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGA\nkCBASBAgJAgQEgQICQKEBAHpkPYPpSw3/32Jq19DSExIOKR9V16t3r6EkJiLcEjr8vSvpqdu\nefwSQmIuwiF1b/9s1y12QmJGwiG9t7NfLoXEjIRDWpT9+62lkJiPcEhP5eG/W7uyFBKzkV7+\nXn/UsylCYjbiJ2S3q/dbuwchMRee2QABQoIAIUFAnyFZbGA2hg2pnPr7AFALD+0gQEgQICQI\niIf08rh6uyRp/XLnyNCe9IV9i5PVhOV9I0N74hf2dc/b463dpivru0aG9sQv7Nt+3N6W7q6R\noT09Xdh3/kFoAKiSeyQIyB8jbXbHW46RmJP08vfyZNVusb/2N4XEhOTPI62P55G61aPzSMyH\nZzZAgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhTUjqlQm8\nwsHvCWkyjnt/IIHU15kXIU1GOXlbw9eZFyFNRfnyfuyvMzNCmgohjUpIUyGkUQlpMhwjjUlI\nk2HVbkxCmhDnkcYjJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAg\nQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBAS\nBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFC\nmp1SbPo8Ic3MsSIpxQlpZsrJW3KENC/ly3tChDQvQuqJkOZFSD0R0sw4RuqHkGbGql0/hDQ7\nziP1QUgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEsOa6DOUhMSQJvucWSExpMlexSEk\nBjTd6wqFxICENODITJeQBhyZCXOMNNzITJhVu+FGZtKcRxpqZGiPkHo30R/BfCKknk32oIBP\nhNSzyS5T8YmQ+jXdEyd8IqR+CWkmhNQvIc2EkHrmGGkehNQzq3bzIKTeOY80B0KCACFBgJAg\nQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQLiIb08rsqr1frlzpGhPeGQ9ovy\nf8v7Rob2hENal+55e7y123RlfdfI0J5wSF3Zftzelu6ukevhInF+KhzSpz3v+m5Y/T7q15bw\nc+6RLvKLtPi5/DHSZne81fwxkl/tyC+kl7+XJ6t2i/1dI1ei/pAcwlUkfx5pfTyP1K0eGz+P\nVHtIDuGq4pkNF1V+jFT59OZGSBfV/SO/9jvMuRHSFTUfhAipLn2G1Ph5pLoJqS7DhlRO/X2A\nWXOMVBUP7VpV9yHcn7T4U1ZI7Wpxf/uBNn9CCInKVPyY9cqPLhf2UZd6V1Gu3lW6sI+6VBzS\nydsLf/ibP7jKhX38UbUhXZ+YyyioTK3HSIOGNKUL+xhHrat27pFoTKXr+gMfI03lwj74bMhV\nuwld2AdfDXoeaSoX9sEveGYDBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGA\nkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBwA9DeurK4mmg\nkaE9t0Larkr3dHj8wauUx0YmptKXvpuiGyFtjwWty8P+sFuV6H2S/+K+1fpirJN0I6SH15ev\nXL+9GOy+LAYZmZCrr3lK1o2Q3n6eldXJB72PTEb58p4+/Sik57fHdNdfpTw2MhlCGtLNh3YP\n76+ovH+4/irlsZHJENKQboS07z4ez5XsHZL/3945RhrQzfNI6/d8uuj9kf/g/lm1G9CNkDYj\njEyM80iDubXY0K13Q48M7bkR0uL1CQ393C0JiQm5dYy0W3evT23YDjgytOcHT1p9efiX0uJp\nf+lvpke++A884KdaP3v29/Py3178kH2I98sqLEFRs59ej7R/XIRPJP02pHv+EQzkFxf2bcZ8\nrp3T9FStlXskIVG1Vo6RhETVfhDSpopVO8dI1OxWSC+v55G6Cs4jWbWjZg09s8F5JOp187l2\nj+mHdDdHhvbcCOllhJGhPX5BJAQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAk\nCBASBAgJAoQEARWG5NXLaU91IR0rkhKNqS+kYUaHqNpCKlf/FColJAgQEgTUFpJjJJpUX0hW\n7WhQdSE5j0SLKgwJ2iMkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGA\nkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBQDqk/UMpy81/\nX+Lq1xASExIOad+VV6u3LyEk5iIc0ro8/avpqVsev4SQmItwSN3bP9t1i52QmJFwSO/t7JdL\nITEj4ZAWZf9+aykk5iMc0lN5+O/WriyFxGykl7/XH/VsipCYjfgJ2e3q/dbuQUjMhWc2QICQ\nIEBIENBnSBYbmA0hQYCHdhAgJAgYNqRyqo8BYBzxkF4eV2+XJK1f7hwZ2pO+sG9xcpezvG9k\naE/8wr7ueXu8tdt0ZX3XyNCe+IV924/b29LdNTK0p6cL+84/CA0AVXKPBAH5Y6TN7njLMRJz\nkl7+Xp6s2i321/6mkJiQ/Hmk9fE8Urd6dB6J+fAUIQgQEgT0FdLtp9IJiQkREgQICQKEBAFC\nggAhQYDlbwgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQI\nEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhMaRy+2VK2iQkhnOsaJopCYnh\nlJO3EyMkBlO+vJ8SITEYIQ06MlMlpEFHZrIcIw05MpNl1W7IkZkw55GGGxnaIyQIEBIECAkC\nhAQBQqrXZFe4pkhItZrwOZcpElKtJvwsgCkSUqWm/Ly0KRJSpYTUFiFVSkhtEVKtHCM1RUi1\nsmrXFCHVy3mkhggJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQ\nEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQB\nQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAg\nQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIENBbSOXW\nFxASEyIkCAiHVD67a2RoTzikl05IzFH6od1+VZa747/30I4ZyR8jPZfyfBAS89LDYsNuWVZ7\nITErvazaPZZuIyTmpJ/l7+3ixkrDXweAuvR1HulBSMyJpwhBgJAgoM+QnJBlNoQEAR7a8Qu3\nF2PnSkj82LEiKX1r2JB+/IxWalRO3vJZPKSXx9Uxk9X65c6RqVT58p4T4ZD2i5O7nOV9I1Mp\nIV0RDmlduuft8dZu05X1XSNTKSFdEQ6pK9uP29vS3TUytXKMdFn8UvNLH4QGYERW7S5zj8Qv\nWGy9JH+MtDleae4YiVlJL38vT1btFvu7Rob25M8jrY/nkbrVo/NIzIenCEGAkCCgr5Bur+4I\niQkREgQICQKEBAFCggAhQYDlbwgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKC\nACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBI\nECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQImHtIpdQxDxo375COFUmJv5t5SCdv4Q9mHVL5\n8h7uJaRDFVOhcUI6VDEVGjfrkBwjkTLzkKzakTHvkJxHImTuIUGEkCBASBAgJAgQEgQICQKE\nBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIGDEkmJDRQmrjPskcIxqY\nYl9zFNKROUY0MEUh9cocIxqYopB6ZY4RDUxRSL0yx4gGpiikXpljRANTFFKvzDGigSkKqVfm\nGNHAFIXUK3OMaGCKQuqVOUY0MEUh9cocIxqYopB6ZY4RDUyx4ZBg+oQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQJ6DGndlW69//y5l8rCPZvj0+J80mM6\nm+B3W3Vc1W/DwxD7Yn879vL4KhiLT5/bd3WFdDbH9fETXTV7wdkEv9uq46p+Gx4G2Rd727Ff\nSrc9bLvycvrJ1eXXlxnD2Ry35eHf//9TeRhzVifOJvjtVh1V9dvwMMy+2NuOvS6bf2+fy+PJ\n556vvFDTGM7muHqbXjWzPJvgd1t1XNVvw8Mw+2Jv3+6q7A6vP55W///Urixr2rzfzvGomlme\nTfDSjMdT/TY8DLMv9vbtlvMfTMuyq2nzfjvHV/uyHGEy3zmb4IUZj6j6bXgYZl8cMKTH8lzT\nHnC4uBM8HR8J1KDdkOrZhodh9sXhQjres1a0Bxwu7QS7rpoHTs2GVNE2PAyzLw4X0uJ1QbSi\nPeBwYSfYd/U8KGk1pJq24WGYfTH+P/L+Iurdl8k/HO/r69gDLs3xaFnRSZqzCX4741FVvw0P\nw+yLvYX0tlKy+1gpKR/SI/7epTm+frBY7kab15mzCX4z45FVvw0Pw+yLve3Wj8fqN2X9PlBF\nIb37OsfX21U9Jjmb4PmMx1b9NjwMsy8O+8yGmh6THL6Z466yfaDFZzbUtg0Pw+yL/e3Yi2Py\nx416drBci69zfKjtXvNsI558ohLVb8PDIPtif9/u/viM27dBag3p6xyre/h5thFPPlGJ6rfh\nYZB9saZvF5olJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECCkhn285tzx\nxv6hlIftmPOZMyE17HNI3fEFJ5U0DiE17FNI6/Lw+mY16ozmS0gN+xRSV/ann2JYtnvDPj+0\ne7vVjTSXuRNSw85DWpensSYzc0Jq2NeQnktZjzebeRNSw76G9LTqyuN405k1ITXsm2OkB4/t\nxiGkhn0T0t5qwziE1LBF2R3f78ri43PWv8dhszfs/XHc0+u52LfzSKdNMSAhNWxTjiU9l/Ly\n3zMb9ivHSOMQUsvW5c1x0fvtuXbLsec0U0Jq2mb1r53V5u2DdVcW7o9GIiQIEBIECAkChAQB\nQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAg\nQEgQICQIEBIECAkChAQB/wOtx2zFtzyrgQAAAABJRU5ErkJggg==",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "test<-data.frame(\n",
    "  X1=c(191, 193, 189, 211, 176, 169, 154, 193, 176, 156, \n",
    "       189, 162, 182, 167, 154, 166, 247, 202, 157, 138), \n",
    "  X2=c(36, 38, 35, 38, 31, 34, 34, 36, 37, 33, \n",
    "       37, 35, 36, 34, 33, 33, 46, 37, 32, 33),\n",
    "  X3=c(50, 58, 46, 56, 74, 50, 64, 46, 54, 54,\n",
    "       52, 62, 56, 60, 56, 52, 50, 62, 52, 68), \n",
    "  Y1=c( 5, 12, 13,  8, 15, 17, 14,  6,  4, 15, \n",
    "        2, 12,  4,  6, 17, 13,  1, 12, 11,  2), \n",
    "  Y2=c(162, 101, 155, 101, 200, 120, 215,  70,  60, 225, \n",
    "       110, 105, 101, 125, 251, 210,  50, 210, 230, 110), \n",
    "  Y3=c(60, 101, 58, 38, 40, 38, 105, 31, 25, 73, \n",
    "       60, 37, 42, 40, 250, 115, 50, 120, 80, 43)\n",
    ")\n",
    "test<-scale(test)\n",
    "\n",
    "ca<-cancor(test[,1:3],test[,4:6]); ca\n",
    "\n",
    "U<-as.matrix(test[, 1:3])%*% ca$xcoef ; U\n",
    "V<-as.matrix(test[, 4:6])%*% ca$ycoef ; V\n",
    "\n",
    "opar <- par(mfrow = c(1, 1),mar = c(5,4,1,1))\n",
    "\n",
    "plot(U[,1], V[,1], xlab=\"U1\", ylab=\"V1\")\n",
    "plot(U[,3], V[,3], xlab=\"U3\", ylab=\"V3\")\n",
    "\n",
    "par(opar)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "4.0.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}