03-22-2019_Number_of_Unique_Family.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "suppressPackageStartupMessages({\n",
    "    library(tidyverse)\n",
    "    library(reshape2)\n",
    "    library(ggpubr)\n",
    "    library(matrixStats)\n",
    "})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "df = read.csv(\"all_samples_fam.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<thead><tr><th scope=col>taxa</th><th scope=col>V1V2_HNO17_S8_R1</th><th scope=col>V1V2_HNO19_S9_R1</th><th scope=col>V1V2_HNO209_S14_R1</th><th scope=col>V1V2_HNO223_S15_R1</th></tr></thead>\n",
       "<tbody>\n",
       "\t<tr><td>Acholeplasmataceae</td><td> 0                </td><td>0                 </td><td>0                 </td><td>0                 </td></tr>\n",
       "\t<tr><td>Actinomycetaceae  </td><td>11                </td><td>9                 </td><td>5                 </td><td>5                 </td></tr>\n",
       "\t<tr><td>Aerococcaceae     </td><td> 1                </td><td>1                 </td><td>0                 </td><td>1                 </td></tr>\n",
       "\t<tr><td>Akkermansiaceae   </td><td> 1                </td><td>0                 </td><td>0                 </td><td>0                 </td></tr>\n",
       "\t<tr><td>Alteromonadaceae  </td><td> 0                </td><td>0                 </td><td>0                 </td><td>0                 </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "\\begin{tabular}{r|lllll}\n",
       " taxa & V1V2\\_HNO17\\_S8\\_R1 & V1V2\\_HNO19\\_S9\\_R1 & V1V2\\_HNO209\\_S14\\_R1 & V1V2\\_HNO223\\_S15\\_R1\\\\\n",
       "\\hline\n",
       "\t Acholeplasmataceae &  0                 & 0                  & 0                  & 0                 \\\\\n",
       "\t Actinomycetaceae   & 11                 & 9                  & 5                  & 5                 \\\\\n",
       "\t Aerococcaceae      &  1                 & 1                  & 0                  & 1                 \\\\\n",
       "\t Akkermansiaceae    &  1                 & 0                  & 0                  & 0                 \\\\\n",
       "\t Alteromonadaceae   &  0                 & 0                  & 0                  & 0                 \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "| taxa | V1V2_HNO17_S8_R1 | V1V2_HNO19_S9_R1 | V1V2_HNO209_S14_R1 | V1V2_HNO223_S15_R1 |\n",
       "|---|---|---|---|---|\n",
       "| Acholeplasmataceae |  0                 | 0                  | 0                  | 0                  |\n",
       "| Actinomycetaceae   | 11                 | 9                  | 5                  | 5                  |\n",
       "| Aerococcaceae      |  1                 | 1                  | 0                  | 1                  |\n",
       "| Akkermansiaceae    |  1                 | 0                  | 0                  | 0                  |\n",
       "| Alteromonadaceae   |  0                 | 0                  | 0                  | 0                  |\n",
       "\n"
      ],
      "text/plain": [
       "  taxa               V1V2_HNO17_S8_R1 V1V2_HNO19_S9_R1 V1V2_HNO209_S14_R1\n",
       "1 Acholeplasmataceae  0               0                0                 \n",
       "2 Actinomycetaceae   11               9                5                 \n",
       "3 Aerococcaceae       1               1                0                 \n",
       "4 Akkermansiaceae     1               0                0                 \n",
       "5 Alteromonadaceae    0               0                0                 \n",
       "  V1V2_HNO223_S15_R1\n",
       "1 0                 \n",
       "2 5                 \n",
       "3 1                 \n",
       "4 0                 \n",
       "5 0                 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df[1:5,1:5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "df2_mat = as.matrix( df[-1])\n",
    "rownames(df2_mat) = df$taxa"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<ol class=list-inline>\n",
       "\t<li>'V1V2'</li>\n",
       "\t<li>'V2V3'</li>\n",
       "\t<li>'V3V4'</li>\n",
       "\t<li>'V4V5'</li>\n",
       "\t<li>'V5V7'</li>\n",
       "\t<li>'V7V9'</li>\n",
       "</ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 'V1V2'\n",
       "\\item 'V2V3'\n",
       "\\item 'V3V4'\n",
       "\\item 'V4V5'\n",
       "\\item 'V5V7'\n",
       "\\item 'V7V9'\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 'V1V2'\n",
       "2. 'V2V3'\n",
       "3. 'V3V4'\n",
       "4. 'V4V5'\n",
       "5. 'V5V7'\n",
       "6. 'V7V9'\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] \"V1V2\" \"V2V3\" \"V3V4\" \"V4V5\" \"V5V7\" \"V7V9\""
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# sets\n",
    "primer_sets = unique( sapply(strsplit(colnames(df2_mat),'_'),'[[',1) )\n",
    "primer_sets"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[1] \"V1V2\"\n",
      "[1] \"V2V3\"\n",
      "[1] \"V3V4\"\n",
      "[1] \"V4V5\"\n",
      "[1] \"V5V7\"\n",
      "[1] \"V7V9\"\n"
     ]
    }
   ],
   "source": [
    "new_df = list();\n",
    "for( i in 1:length(primer_sets) ){\n",
    "    primer= primer_sets[i]\n",
    "    print(primer)\n",
    "    new_df[[primer]] = rowMeans(df2_mat[ , grep( primer, colnames(df2_mat))])\n",
    "    }"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "new_df = as.data.frame( new_df )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<thead><tr><th></th><th scope=col>V1V2</th><th scope=col>V2V3</th><th scope=col>V3V4</th><th scope=col>V4V5</th><th scope=col>V5V7</th><th scope=col>V7V9</th></tr></thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>Acholeplasmataceae</th><td>0.04347826</td><td>0.04347826</td><td>0.04347826</td><td>0.04347826</td><td>0.04347826</td><td>0.0000000 </td></tr>\n",
       "\t<tr><th scope=row>Actinomycetaceae</th><td>6.73913043</td><td>7.39130435</td><td>6.39130435</td><td>6.13043478</td><td>8.47826087</td><td>2.9130435 </td></tr>\n",
       "\t<tr><th scope=row>Aerococcaceae</th><td>0.52173913</td><td>0.08695652</td><td>0.34782609</td><td>0.34782609</td><td>1.17391304</td><td>0.4782609 </td></tr>\n",
       "\t<tr><th scope=row>Akkermansiaceae</th><td>0.39130435</td><td>0.08695652</td><td>0.08695652</td><td>0.08695652</td><td>0.17391304</td><td>0.0000000 </td></tr>\n",
       "\t<tr><th scope=row>Alteromonadaceae</th><td>0.00000000</td><td>0.00000000</td><td>0.00000000</td><td>0.04347826</td><td>0.00000000</td><td>0.0000000 </td></tr>\n",
       "\t<tr><th scope=row>Anaerolineaceae</th><td>0.00000000</td><td>0.04347826</td><td>0.00000000</td><td>0.00000000</td><td>0.00000000</td><td>0.0000000 </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "\\begin{tabular}{r|llllll}\n",
       "  & V1V2 & V2V3 & V3V4 & V4V5 & V5V7 & V7V9\\\\\n",
       "\\hline\n",
       "\tAcholeplasmataceae & 0.04347826 & 0.04347826 & 0.04347826 & 0.04347826 & 0.04347826 & 0.0000000 \\\\\n",
       "\tActinomycetaceae & 6.73913043 & 7.39130435 & 6.39130435 & 6.13043478 & 8.47826087 & 2.9130435 \\\\\n",
       "\tAerococcaceae & 0.52173913 & 0.08695652 & 0.34782609 & 0.34782609 & 1.17391304 & 0.4782609 \\\\\n",
       "\tAkkermansiaceae & 0.39130435 & 0.08695652 & 0.08695652 & 0.08695652 & 0.17391304 & 0.0000000 \\\\\n",
       "\tAlteromonadaceae & 0.00000000 & 0.00000000 & 0.00000000 & 0.04347826 & 0.00000000 & 0.0000000 \\\\\n",
       "\tAnaerolineaceae & 0.00000000 & 0.04347826 & 0.00000000 & 0.00000000 & 0.00000000 & 0.0000000 \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "| <!--/--> | V1V2 | V2V3 | V3V4 | V4V5 | V5V7 | V7V9 |\n",
       "|---|---|---|---|---|---|---|\n",
       "| Acholeplasmataceae | 0.04347826 | 0.04347826 | 0.04347826 | 0.04347826 | 0.04347826 | 0.0000000  |\n",
       "| Actinomycetaceae | 6.73913043 | 7.39130435 | 6.39130435 | 6.13043478 | 8.47826087 | 2.9130435  |\n",
       "| Aerococcaceae | 0.52173913 | 0.08695652 | 0.34782609 | 0.34782609 | 1.17391304 | 0.4782609  |\n",
       "| Akkermansiaceae | 0.39130435 | 0.08695652 | 0.08695652 | 0.08695652 | 0.17391304 | 0.0000000  |\n",
       "| Alteromonadaceae | 0.00000000 | 0.00000000 | 0.00000000 | 0.04347826 | 0.00000000 | 0.0000000  |\n",
       "| Anaerolineaceae | 0.00000000 | 0.04347826 | 0.00000000 | 0.00000000 | 0.00000000 | 0.0000000  |\n",
       "\n"
      ],
      "text/plain": [
       "                   V1V2       V2V3       V3V4       V4V5       V5V7      \n",
       "Acholeplasmataceae 0.04347826 0.04347826 0.04347826 0.04347826 0.04347826\n",
       "Actinomycetaceae   6.73913043 7.39130435 6.39130435 6.13043478 8.47826087\n",
       "Aerococcaceae      0.52173913 0.08695652 0.34782609 0.34782609 1.17391304\n",
       "Akkermansiaceae    0.39130435 0.08695652 0.08695652 0.08695652 0.17391304\n",
       "Alteromonadaceae   0.00000000 0.00000000 0.00000000 0.04347826 0.00000000\n",
       "Anaerolineaceae    0.00000000 0.04347826 0.00000000 0.00000000 0.00000000\n",
       "                   V7V9     \n",
       "Acholeplasmataceae 0.0000000\n",
       "Actinomycetaceae   2.9130435\n",
       "Aerococcaceae      0.4782609\n",
       "Akkermansiaceae    0.0000000\n",
       "Alteromonadaceae   0.0000000\n",
       "Anaerolineaceae    0.0000000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "new_df %>% head"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "67"
      ],
      "text/latex": [
       "67"
      ],
      "text/markdown": [
       "67"
      ],
      "text/plain": [
       "[1] 67"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "length( unique((rownames(new_df[ new_df$V1V2 !=0,][\"V1V2\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "63"
      ],
      "text/latex": [
       "63"
      ],
      "text/markdown": [
       "63"
      ],
      "text/plain": [
       "[1] 63"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "length( unique((rownames(new_df[ new_df$V2V3 !=0,][\"V1V2\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "64"
      ],
      "text/latex": [
       "64"
      ],
      "text/markdown": [
       "64"
      ],
      "text/plain": [
       "[1] 64"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "length( unique((rownames(new_df[ new_df$V4V5 !=0,][\"V1V2\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "67"
      ],
      "text/latex": [
       "67"
      ],
      "text/markdown": [
       "67"
      ],
      "text/plain": [
       "[1] 67"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "length( unique((rownames(new_df[ new_df$V3V4 !=0,][\"V1V2\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "67"
      ],
      "text/latex": [
       "67"
      ],
      "text/markdown": [
       "67"
      ],
      "text/plain": [
       "[1] 67"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "length( unique((rownames(new_df[ new_df$V3V4 !=0,][\"V1V2\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "v1v2 = length( unique((rownames(new_df[ new_df$V1V2 !=0,][\"V1V2\"]))) )\n",
    "v2v3 = length( unique((rownames(new_df[ new_df$V2V3 !=0,][\"V2V3\"]))) )\n",
    "v3v4 = length( unique((rownames(new_df[ new_df$V3V4 !=0,][\"V3V4\"]))) )\n",
    "v4v5 = length( unique((rownames(new_df[ new_df$V4V5 !=0,][\"V4V5\"]))) )\n",
    "v5v7 = length( unique((rownames(new_df[ new_df$V5V7 !=0,][\"V5V7\"]))) )\n",
    "v7v9 = length( unique((rownames(new_df[ new_df$V5V7 !=0,][\"V7V9\"]))) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAHgCAIAAAD2dYQOAAAABmJLR0QA/wD/AP+gvaeTAAAg\nAElEQVR4nO3df1RUdf748dcwyIiDkKJoxgBxiDItsXLdpD3nM2lHah3oVLtulmXksi2nZfux\n4WptRzez2ErqbLkr2ZKbpz1lv/ixoKUznUIp8BRZTnJcBVxYKTdAWMcBgvn+Md/lsIgzIAxz\nL+/n46+4c9/D6z24PLnMMGvweDwCAICqQoI9AAAAwUQIAQBKI4QAAKURQgCA0gghAEBphBAA\noDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIA\nlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA\n0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQ\nGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBK\nI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBp\nhBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiN\nEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKUR\nQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRC\nCABQGiEEACgtNNgDjE89PT1lZWVutzvYgwzDxIkTb7rpJqPRGOxBxpoev1ii8NcLGHWEMCDK\ny8vT09ODPcWwlZSULFu2LNhTjDWdfrFE1a8XMOoIYUCcOXNGRBITE6dMmRLsWYaktbX12LFj\n3rFVo7svlqj99QJGHSEMoClTpsyaNSvYU2BI+GIByuLFMgAApRFCAIDSCCEAQGmEEACgNF4s\nA0A3xusffepxX+PpL1kJIQDdGK9/9KnTfY2bv2QlhAB0Y7z+0afu9jXO/pKVEALjkB5/1SZD\n/m3beP2jz/G6L+0jhMA4pNNftck4+m0bdIQQAuOQ7n7VJuPut23QEUIIjFv8qg0YCkKIoRrf\nTzsBUBYhxFDxtBOAcYkQYqh42gnAuEQIMTw87QRgnOG9RgEASiOEAAClEUIAgNL0EcLS0tL/\n+7//i46OnjRp0vz587du3dr/VpfLtX79+uTk5PDw8NjY2MzMzMbGxmCNCgDQFx28WOb5559/\n8MEH+z6sqam57777uru777//fhHp7u622Wx2u917a1NTU2FhYWlpaXV1dXx8fHAmBgDoh9av\nCI8cOfKb3/zGYDA8+eST//73v7/99ttnn302JCRk/fr13tfEFxQU2O12i8Vit9vdbvfRo0cz\nMjJOnjyZnZ0d7NkBADqg9RBu3ry5p6dnzZo169ati46Onj59+sMPP7x8+fLvvvvuq6++EpEt\nW7aIyPbt261Wq8lkSkxM3LFjh9lsLi8vP3HiRLDHBwBondZD+P7774eGhj788MP9D77++use\nj2fBggXNzc1Op9NisVit1r5bIyIilixZ4vF4KioqxnxeAIDOaPo5wpaWlmPHjs2fP7+np+fO\nO+8sLy93uVxz587Nzs5etWqVwWCora0VkZSUlAEL586dW1RUVFdX5/dTDOX9M3t7e2tra3/7\n29+GhYWNZDsAAA3SdAhPnjwpItHR0YsWLTp27Jj34IEDBzIzMz/++OO//OUvLS0tIjJ9+vQB\nC6dNmyYiHR0dfj+Fw+EY4vtnhoSEPPbYY8OaHwCgfZoOYVtbm4js2bNn3rx5e/fuXbhwYXd3\n9zvvvPPLX/6ysLBw1apVXV1dgy40GAwiYjab/X4Kq9VaXFzs+4qwsrIyPz//kksuOa9NAAA0\nTdMhNJlMImIwGMrKyvre3zIzM9PpdD733HN79uxZtGiRiHivC/trbW0VkZiYGL+fwmg02mw2\nv6fl5+eHhGj9+VQAwHnQ9Df3Cy+8UERiYmIGvMvzNddcIyKnTp1KSkoSkZqamgELnU6niHhv\nBQDAB02HcMaMGXFxcd988433RTF9PvnkExG55JJLkpKSYmNj6+vrq6qq+m5taWnZvXt3dHR0\namrqWE8MANAbTYdQRO69914RSU9Pf//999vb25ubm1944YWXXnrJbDb/9Kc/FZGsrCwRWbFi\nxf79+zs7Ow8ePGiz2drb21euXMn/KTkAwC9NP0coImvWrNm7d+9HH320dOnSvoMhISEvv/yy\n9ynA3NzckpKS6urq/td/s2fP3rBhQxDGBQDojdavCE0m0wcffPD0009ffvnlYWFhF1xwQVpa\nmsPhuP322/tOcDgcubm5CQkJYWFhcXFxOTk5+/bti4yMDO7kAABd0PoVoYiEhYWtWbNmzZo1\n5zrBbDbn5eXl5eWN5VQAgPFB61eEAAAEFCEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBp\nhBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiN\nEAIAlEYIAQBKI4QAAKURQgCA0gghAEBpOghhSkqK4SxpaWl9J7hcrvXr1ycnJ4eHh8fGxmZm\nZjY2NgZxYACAjoQGewA/PB7PkSNHfJzQ3d1ts9nsdrv3w6ampsLCwtLS0urq6vj4+DGZEQCg\nY1q/Ijx+/LjL5Vq7dq3nf+3atct7QkFBgd1ut1gsdrvd7XYfPXo0IyPj5MmT2dnZwZ0cAKAL\nWg9hbW2tiFx11VXnOmHLli0isn37dqvVajKZEhMTd+zYYTaby8vLT5w4MXaDAgD0SeshPHz4\nsIi0tbUtWrRo0qRJU6dOvfnmmz/77DPvrc3NzU6n02KxWK3WviURERFLlizxeDwVFRXBGRoA\noB9af47Qe0WYlZXl8XhE5MyZM0VFReXl5Tt37kxPT/fempKSMmDV3Llzi4qK6urq/N5/T09P\nWVmZ2+32cU5lZaWI9Pb2nvcuAACapYMQhoaGPvLII/fdd19MTExtbe3jjz9eXFx81113HTt2\nrKWlRUSmT58+YNW0adNEpKOjw+/9OxyO9PT0oUxy6NCh4Y8PANA6rYdwz549/T+cN2/e22+/\nvWjRourq6uLi4vDw8EFXGQwGETGbzX7v32q1FhcX+70izM/PnzNnznAGBwDog9ZDeLbQ0FCb\nzVZdXV1XV3fttdeKiPe6sL/W1lYRiYmJ8XtvRqPRZrP5PS0/Pz8kROvPpwIAzoMuv7l3d3eL\nyOTJk5OSkkSkpqZmwAlOp1NEvLcCAOCDpkO4d+9eg8GQkZHR/2B3d/ebb74pIgsXLkxKSoqN\nja2vr6+qquo7oaWlZffu3dHR0ampqWM9MQBAbzQdwmuvvTYyMrK4uHjdunWHDx92uVxffPGF\n98Wi8+bNu+6660QkKytLRFasWLF///7Ozs6DBw/abLb29vaVK1cajcZg7wAAoHWafo5w0qRJ\nL7744qpVq5566qmnnnqq73hUVNSrr77qfUVMbm5uSUlJdXV1/+u/2bNnb9iwIQgTAwD0RtNX\nhCKycuXKAwcO3HbbbbGxsWFhYRaL5Z577jlw4EDf3w6aTCaHw5Gbm5uQkBAWFhYXF5eTk7Nv\n377IyMjgTg4A0AVNXxF6zZ8/f+fOnT5OMJvNeXl5eXl5YzYSAGDc0PoVIQAAAUUIAQBKI4QA\nAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAA\noDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQms5C+NZbb4WE\nhKxevbr/QZfLtX79+uTk5PDw8NjY2MzMzMbGxmBNCADQl9BgDzAMBw4cuOuuuzweT/+D3d3d\nNpvNbrd7P2xqaiosLCwtLa2uro6Pjw/GmAAAPdHNFWFTU1NGRkZkZOSA4wUFBXa73WKx2O12\nt9t99OjRjIyMkydPZmdnB2VOAIC+6COELpcrPT19woQJW7duHXDTli1bRGT79u1Wq9VkMiUm\nJu7YscNsNpeXl584cSIYwwIA9EQHIfR4PHfdddeRI0dKSkpmzJjR/6bm5man02mxWKxWa9/B\niIiIJUuWeDyeioqKMR8WAKAzOniO8LHHHnv33XeLioquuOKKTz75pP9NtbW1IpKSkjJgydy5\nc4uKiurq6vzeeU9PT1lZmdvt9nFOZWWliPT29g57dACA5mk9hDt27Ni0adNzzz23bNmys29t\naWkRkenTpw84Pm3aNBHp6Ojwe/8OhyM9PX0okxw6dGgopwEA9EXTIayqqlq9evXq1asfeuih\nQU/o6uoa9LjBYBARs9ns91NYrdbi4mK/V4T5+flz5swZwsgAAJ3RdAjff//9zs7Obdu2bdu2\nrf/xV1555ZVXXlm6dOkDDzwg/70u7K+1tVVEYmJi/H4Ko9Fos9n8npafnx8SooPnUwEAw6Xv\nb+5JSUkiUlNTM+C40+nsuxUAAB80HcLHHnvM878+//xzEbn33ns9Hs+uXbuSkpJiY2Pr6+ur\nqqr6VrW0tOzevTs6Ojo1NTV4swMA9EHTIRyKrKwsEVmxYsX+/fs7OzsPHjxos9na29tXrlxp\nNBqDPR0AQOs0/RzhUOTm5paUlFRXV/e//ps9e/aGDRuCOBUAQC90f0VoMpkcDkdubm5CQkJY\nWFhcXFxOTs6+ffvOfjM2AADOprMrwpSUlAFvui0iZrM5Ly8vLy8vKCMBAHRN91eEAACMBCEE\nACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QA\nAKURQgCA0gghAEBphBAAoDRCCABQGiEEACiNEAIAlEYIAQBKI4QAAKURQgCA0gghAEBphBAA\noDQdhPCDDz5YunRpdHT0xIkTExISMjIy9u7d2/8El8u1fv365OTk8PDw2NjYzMzMxsbGYE0L\nANCX0GAP4MfOnTuXL1/u8Xi8HzY0NDQ0NBQXFxcUFPz85z8Xke7ubpvNZrfbvSc0NTUVFhaW\nlpZWV1fHx8cHbW4AgE5o/YowNzfX4/Fs2rSpqampq6vr2LFjDzzwgIg8++yz3hMKCgrsdrvF\nYrHb7W63++jRoxkZGSdPnszOzg7q4AAAfdB0CJubm+vr6xcvXrx27dpZs2ZNmDDh4osv3rx5\n86RJk/7zn/94z9myZYuIbN++3Wq1mkymxMTEHTt2mM3m8vLyEydOBHV8AIAOaDqEM2fO9Hg8\ne/bs8X7Y09PT0NCwbt06l8uVnp4uIs3NzU6n02KxWK3WvlURERFLlizxeDwVFRXBmRsAoB9a\nf46wz3333bd161bvf99xxx2bN28WkdraWhFJSUkZcPLcuXOLiorq6ur83m1PT09ZWZnb7fZx\nTmVlpYj09vae3+QAAC3TTQj7e/PNNxctWpSdnd3S0iIi06dPH3DCtGnTRKSjo8PvXTkcDu/F\npV+HDh0a/qQAAK3TTQj//Oc/5+fnHz9+vLCwMC8vLycnZ+nSpV1dXYOebDAYRMRsNvu9W6vV\nWlxc7PeKMD8/f86cOec3OQBAy3QTQhEJDw+/9NJLn3766fr6+jfeeKO8vDwpKUlEvNeF/bW2\ntopITEyM3/s0Go02m83vafn5+SEhmn4+FQBwfjT9zX3jxo0Gg6HvLyX6eC/O2travCGsqakZ\ncILT6RQR760AAPig6RBedtllIjLgfWRE5OOPPxaRpKSkpKSk2NjY+vr6qqqqvltbWlp2794d\nHR2dmpo6ltMCAPRI0yFMS0uLioratWvXo48+2tTU1NnZ6XQ6V61a9cEHH0ydOvWmm24Skays\nLBFZsWLF/v37Ozs7Dx48aLPZ2tvbV65caTQag70DAIDWafo5woiIiBdeeOGee+7ZtGnTpk2b\n+o6HhYW9+uqrkZGRIpKbm1tSUlJdXd3/+m/27NkbNmwIwsQAAL3R9BWhiNx9991///vfFy9e\nfMEFF0yYMGHWrFm33357VVVV3ytcTCaTw+HIzc1NSEgICwuLi4vLycnZt2+fN5MAAPim6StC\nrxtvvPHGG2/0cYLZbM7Ly8vLyxuzkQAA44bWrwgBAAgoQggAUBohBAAojRACAJRGCAEASiOE\nAAClEUIAgNIIIQBAaYQQAKA0QggAUBohBAAojRACAJRGCAEASiOEAAClEUIAgNIIIQBAaYQQ\nAKA0QggAUBohBAAojRACAJRGCAEASiOEAAClEUIAgNJ0EMJTp0498sgjycnJ4eHhkydPXrhw\n4bZt2/qf4HK51q9f7z0hNjY2MzOzsbExWNMCAPQlNNgD+HHq1KmFCxfW1tb2HamqqqqqqnI6\nnZs3bxaR7u5um81mt9u9tzY1NRUWFpaWllZXV8fHxwdnaACAfmj9inDz5s21tbVXXnllRUXF\nmTNnmpqaNmzYICLPP//8F198ISIFBQV2u91isdjtdrfbffTo0YyMjJMnT2ZnZwd7dgCADmg9\nhG+//baI7Ny5MzU1deLEibNmzXr88ceXL1/u8XjKyspEZMuWLSKyfft2q9VqMpkSExN37Nhh\nNpvLy8tPnDgR5OkBAJqn9RDW1dVddNFFycnJ/Q9efvnlIuJ2u5ubm51Op8VisVqtfbdGREQs\nWbLE4/FUVFSM9bgAAL3R+nOEp0+fHnCkq6vrvffeE5HrrrvO+9xhSkrKgHPmzp1bVFRUV1fn\n9/57enrKysrcbrePcyorK0Wkt7d3WJMDAHRB6yEcoKGhITMz8/PPP8/IyLjhhhveffddEZk+\nffqA06ZNmyYiHR0dfu/Q4XCkp6cP5VMfOnRo+PMCALRONyFsa2v7wx/+8Pzzz3d2dubk5Dzz\nzDMi0tXVNejJBoNBRMxms9+7tVqtxcXFfq8I8/Pz58yZc16DAwA0TQch9Hg8W7Zs+d3vftfa\n2pqWlvbUU0/1/S40KipKRFpaWgYsaW1tFZGYmBi/d240Gm02m9/T8vPzQ0K0/nwqAOA8aD2E\n3d3dy5cvf/fddxcsWJCfn5+amtr/1qSkJBGpqakZsMrpdPbdCgCAD1q/ytm4ceO77767atWq\nysrKARUUkaSkpNjY2Pr6+qqqqr6DLS0tu3fvjo6OPvt8AAAG0HQIe3p6CgoK4uPjX375ZaPR\nOOg5WVlZIrJixYr9+/d3dnYePHjQZrO1t7evXLnyXEsAAOij6RB+/fXXzc3NDQ0NEyZMMJxl\n48aNIpKbm7tgwYKjR496/+J+3rx5+/fvnz17tvcNaAAA8E3TIayvr/d7jslkcjgcubm5CQkJ\nYWFhcXFxOTk5+/bti4yMDPyAAADd0/SLZZYtW+bxePyeZjab8/Ly8vLyxmAkAMA4o+krQgAA\nAo0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAA\npRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACg\nND2FsLi42GAwfPjhhwOOu1yu9evXJycnh4eHx8bGZmZmNjY2BmNAAID+hAZ7gGEoLCw8+2B3\nd7fNZrPb7d4Pm5qaCgsLS0tLq6ur4+Pjx3ZAAID+6OCKsLW11eFwrFq16r333jv71oKCArvd\nbrFY7Ha72+0+evRoRkbGyZMns7Ozx35UAIDuaP2K0O12T5061ccJW7ZsEZHt27dbrVYRSUxM\n3LFjx8yZM8vLy0+cOHHhhReO0aAAAH3S+hXhxIkTPf916623Dri1ubnZ6XRaLBZvBb0iIiKW\nLFni8XgqKirGdlgAgP5o/YrQt9raWhFJSUkZcHzu3LlFRUV1dXV+76Gnp6esrMztdvs4p7Ky\nUkR6e3tHMCkAQKP0HcKWlhYRmT59+oDj06ZNE5GOjg6/9+BwONLT04fyuQ4dOjT8AQEAWqfv\nEHZ1dQ163GAwiIjZbPZ7D1artbi42O8VYX5+/pw5c85vSACAluk7hFFRUfLf68L+WltbRSQm\nJsbvPRiNRpvN5ve0/Pz8kBCtP58KADgP+v7mnpSUJCI1NTUDjjudzr5bAQDwQfchjI2Nra+v\nr6qq6jvY0tKye/fu6Ojo1NTUIM4GANAFfYdQRLKyskRkxYoV+/fv7+zsPHjwoM1ma29vX7ly\npdFoDPZ0AACt0/dzhCKSm5tbUlJSXV3d//pv9uzZGzZsCOJUAAC90P0Voclkcjgcubm5CQkJ\nYWFhcXFxOTk5+/bti4yMDPZoAAAd0NMV4VtvvTXocbPZnJeXl5eXN8bzAADGAd1fEQIAMBKE\nEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0Q\nAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFC\nAIDSxkMIXS7X+vXrk5OTw8PDY2NjMzMzGxsbgz0UAEAfQoM9wEh1d3fbbDa73e79sKmpqbCw\nsLS0tLq6Oj4+PrizAQC0T/dXhAUFBXa73WKx2O12t9t99OjRjIyMkydPZmdnB3s0AIAO6D6E\nW7ZsEZHt27dbrVaTyZSYmLhjxw6z2VxeXn7ixIlgTwcA0Dp9h7C5udnpdFosFqvV2ncwIiJi\nyZIlHo+noqIiiLMBAHRB388R1tbWikhKSsqA43Pnzi0qKqqrq/N7Dz09PWVlZW6328c5lZWV\nItLb2zvc8VpbW4e7JFiGPqqONiXK72tcbmpYZ2rBuNyXjkYdCn2HsKWlRUSmT58+4Pi0adNE\npKOjw+89OByO9PT0oXyuI0eODH2w8PBwETl27NjQl2iBd2zft+puU6LkvsblpoR9aYzvTemI\nvkPY1dU16HGDwSAiZrPZ7z1Yrdbi4mLfV4S9vb21tbW5ublDH+zGG2/0e7daM3HixBtvvNHH\nCXrclKi6r3G5KWFfWuJ3Uzqi7xBGRUXJf68L+/NetsfExPi9B6PRaLPZRn2wAN1tcI3LTck4\n3de43JSwLwSGvl8sk5SUJCI1NTUDjjudzr5bAQDwweDxeII9w4hYLJbGxsZPP/30Bz/4gfdI\nS0vLxRdfPGHChG+++cZoNAZ3PACAxun7ilBEsrKyRGTFihX79+/v7Ow8ePCgzWZrb29fuXIl\nFQQA+KX7K8LOzs4f/ehH1dXV/Q/Onj37k08+iYyMDNZUAAC90P0Voclkcjgcubm5CQkJYWFh\ncXFxOTk5+/bto4IAgKHQ/RUhAAAjofsrQgAARoIQAgCURggBAEojhAAApRFCAIDSCCEAQGmE\nEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0Q\nAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0Qal1xcbHBYPjwww+HterU\nqVOPPPJIcnJyeHj45MmTFy5cuG3btsAMODwjHOyDDz5YunRpdHT0xIkTExISMjIy9u7dG7hp\nx36qt956KyQkZPXq1aM74cidx2ApKSmGs6SlpQVuyIBO9eGHH569UAv7GpXBamtrb7/99pkz\nZ5pMposvvvhXv/pVc3NzoCfXEA+07eabbxYRh8Mx9CVtbW2XXnrp2V/rBx98MGBjjsVgb775\npsFgOHt5QUFBoCcfm6mqq6vDw8NF5N577w3EqOftPAbr7e2dNGnS2Q/L0qVLAzpq4KZyOBw+\nvpHecsstYzB/gAb78ssvo6KiBiyMj4//17/+NQbzawEh1KiWlha73X733Xd7/1EOK4SPP/64\niFx55ZUVFRVnzpxpamrasGGDiBgMhpqamoCNHPDBEhISRGTTpk1NTU1dXV3Hjh174IEHRCQ5\nOXkMhg/0VI2NjbNmzZoxY4bWQnh+g9XX14vI2rVrAzrbcI36VD09PWlpaYmJid98881o3eeo\nGNZgixYtEpFbb7310KFDLpfryy+/TE9PF5E77rhjDEbVAkIYBN5fVrz22mv9D/7zn/80GAyX\nXnqpx+M5c+bMgJ/O+ofQ7/I5c+aISG1tbf8Tli9f7v1+HcR9+R7M9/ITJ06IyOLFi/vf6v0B\nf9asWQHbk59N+Z3K72Pidfr06auuuio+Pv69994bmxCOcDC/y3fv3i0iO3fuDPRG+hvhVEN8\nTPpbu3atyWT6/PPPR2kHgxvhYL6XNzQ0iEhSUlJXV1ffrW63+8ILL5w4ceLp06cDsCHN4TnC\nILjzzjtF5J133ul/0Ps/zrvuuktEJk6c2PcVuvXWW4e7vK6u7qKLLkpOTu5/wuWXXy4ibrd7\n9PczSoP5Xj5z5kyPx7Nnzx7v8Z6enoaGhnXr1rlcLu9Pr0HZlN+p/D4mIuL97yNHjpSUlHgv\nvMbACAfzu/zw4cMi0tbWtmjRokmTJk2dOvXmm2/+7LPPArmnkU41lMekv3379j399NOPPvpo\nSkpKYDYkozKY7+Xeq+RrrrlmwoQJfbeaTKbk5GS32/3pp5+O/n40aAyji//v9OnTERERkyZN\n6v/T1g9/+EODwdDQ0DDgZG8I+18RDmu5V2dn5/z580Xk/fffH82d/K8RDjb05b/4xS/6/gHf\ncccdLpcrQDvyDGdTg041lOXr1q0LCQkpKSnxeDyVlZUyJleEIxzM7/Ls7GwRGfDsaVhYWFFR\nURA35XuqYf0Ddrvdl1122aWXXtrZ2Rm4HQ1xX74H87384MGDIjJ//vz+d9LZ2Tlz5kwRef31\n1wO2LQ0hhMHh/Rnt7bff9n54/Phxg8FgtVrPPvPsEA5rucfjqa+vv/7660UkIyNjNPcwmBEO\nNsTl/ZMzYcKEl156KTC78YzKVL6Xv/baayLy3HPPeT8csxCOfDDfyxcvXlpXgxYAAAYESURB\nVBwaGrp27dqGhoYzZ87U1NR4L5GjoqK+++67YG3K71RD/wecl5cnIm+88Ubg9tLfCAfzsbyn\npycuLk5EcnJy6urqvA/LsmXLvP+St27dGuCdaQIhDI5du3aJyJ133un98LnnnhORwsLCs88c\nNIRDXN7a2rp27drw8PCQkJCcnJwx+NF1hIMN/WFxuVyHDx9es2aNiBiNxn/84x8B2tHIp/Kx\n/NNPPzWZTKtXr+67h7EM4QgHG/rD4tXd3b1gwQLf54zcCKca+j/gKVOmXH755b29vQHYxCBG\nOJjv5aWlpaGhof2vkkNDQ70tfOuttwK6L40ghMHx/fffz5w584ILLvA+Qb1w4cJJkyZ1dHSc\nfeagIfS7vLe398UXX5wyZYqIpKWlBfrJ/NEabOgPSx/va23++Mc/BmZDozCVj+VPPPGEnFug\n/9JghIOdx8Py+9//XkQef/zxoGxqKFMNcfmTTz4pY/t3OyMczO/yioqKpUuXRkVFhYeHp6am\nfvjhhz/5yU9EpKKiIsA70wReLBMcRqPxZz/7WVtbm91uP378eFVV1S233BIRETEqy7u7u2+9\n9db7778/KSmpoqKivLw80E/mj9ZgPpZv3LjRYDA8++yzAz6j95WobW1tQdnUUKYa4dc6cAL6\nj3BQ3d3dIjJ58uRRmD4wUw1l+ffff/+nP/0pIiLijjvuCNxGBhjhYH6Xp6am7tq1q62tzeVy\nVVRUeFsYHh5+9dVXj8X2gi7YJVbXgQMHRCQrK+uZZ56Rc7+MZdArQt/LvX+ut2rVqu+//z5w\n85/LCAc71/KdO3eKSFpa2oDzb7jhBhH529/+Foi9jNZUQ/xaezyezz//XMbw7whHONi5lntf\nRpuent7/5K6uLu/bKXz00UeB2MtoTeX3MSkvLxeRFStWBHQXZxvhYOda3t7eLiIXX3xx/5O3\nb98uIj/+8Y8DsRENIoTBNHv27BkzZlxzzTUXXXRRT0/PoOecK4TnWu79HUh8fHx3d3fABvdj\nhIMNuryjo8P75hfr1q1rbGx0u92HDh3yvuHA1KlTT506FdAdjXyqoXytPWMewpEPNujy06dP\nR0ZGisjatWu//vrr06dP19TUeP+abd68eWPwvNoIp/L9mHjfZy4oL6cc4WDnWn7FFVeIyK9/\n/ev6+vpvv/32xRdf9L6RUKB/ZNEOQhhMGzdu9F6X5+bmnuscHyEcdPmXX37p4xcATzzxRCA2\nMrqDnethefXVV89+M7OwsLDi4uJgbWroUw3la+0JRghHONi5lv/1r38NCRn4zEtUVNTYPF09\nwql8Pybe11g2NjYGdg+DGeFg51peXl5uNBoHPCwPPfTQ6G9AqwhhMNXX13u/h3711VfnOsdH\nCAddXlJSMnhqRGSsQjjCwXw8LGVlZYsXL77gggsmTJgwa9as22+/fczeNG6EUw3la+0JRghH\nOJiP5Z999tltt90WGxsbFhZmsVjuueeeI0eOjPL05zDCqXws976XUGxsbACnP7cRDuZjud1u\nv/766ydPnmw2m6+++uqAvrJXgwwej8fHtycAAMY3XjUKAFAaIQQAKI0QAgCURggBAEojhAAA\npRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACg\nNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCU\nRggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDS\nCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAa\nIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEoj\nhAAApRFCAIDSCCEAQGmEEACgNEIIAFAaIQQAKI0QAgCURggBAEojhAAApRFCAIDSCCEAQGmE\nEACgNEIIAFAaIQQAKI0QAgCURggBAEr7f5nIE5cwCO9zAAAAAElFTkSuQmCC",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "options(repr.plot.height=4, repr.plot.width=5)\n",
    "barplot( cbind( v1v2, v2v3, v3v4, v4v5, v5v7, v7v9 ) , las=1)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R 3.5",
   "language": "R",
   "name": "ir35"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}