3.visualize-coefficients.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Gregory Way 2018\n",
    "# Multiple Myeloma Classifier\n",
    "# 3.visualize_coefficients.ipynb\n",
    "#\n",
    "# Observes differences in coefficients across the multiclass classifier classes\n",
    "#\n",
    "# Usage: Run in command line\n",
    "#\n",
    "#     Rscript --vanilla visualize_coefficients.R\n",
    "#\n",
    "# Output:\n",
    "# Descriptive Coefficient Plot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message:\n",
      "“package ‘dplyr’ was built under R version 3.4.3”\n",
      "Attaching package: ‘dplyr’\n",
      "\n",
      "The following objects are masked from ‘package:stats’:\n",
      "\n",
      "    filter, lag\n",
      "\n",
      "The following objects are masked from ‘package:base’:\n",
      "\n",
      "    intersect, setdiff, setequal, union\n",
      "\n",
      "Warning message:\n",
      "“package ‘ggplot2’ was built under R version 3.4.3”Warning message:\n",
      "“package ‘ggrepel’ was built under R version 3.4.3”"
     ]
    }
   ],
   "source": [
    "library(dplyr)\n",
    "library(ggplot2)\n",
    "library(ggrepel)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<thead><tr><th scope=col>num</th><th scope=col>GENE_ID</th><th scope=col>co</th><th scope=col>SYMBOL</th><th scope=col>ALT</th><th scope=col>a</th><th scope=col>b</th></tr></thead>\n",
       "<tbody>\n",
       "\t<tr><td>1                                                                                             </td><td>ENSG00000000003                                                                               </td><td>1.1                                                                                           </td><td>TSPAN6                                                                                        </td><td>TSPAN6                                                                                        </td><td>tetraspanin 6 [Source:HGNC Symbol;Acc:HGNC:11858]                                             </td><td>ARRAYEXPRESS, ENSG                                                                            </td></tr>\n",
       "\t<tr><td>2                                                                                             </td><td>ENSG00000000005                                                                               </td><td>2.1                                                                                           </td><td>TNMD                                                                                          </td><td>TNMD                                                                                          </td><td>tenomodulin [Source:HGNC Symbol;Acc:HGNC:17757]                                               </td><td>ARRAYEXPRESS, ENSG                                                                            </td></tr>\n",
       "\t<tr><td>3                                                                                             </td><td>ENSG00000000419                                                                               </td><td>3.1                                                                                           </td><td>DPM1                                                                                          </td><td>DPM1                                                                                          </td><td>dolichyl-phosphate mannosyltransferase subunit 1, catalytic [Source:HGNC Symbol;Acc:HGNC:3005]</td><td>ARRAYEXPRESS, ENSG                                                                            </td></tr>\n",
       "\t<tr><td>4                                                                                             </td><td>ENSG00000000457                                                                               </td><td>4.1                                                                                           </td><td>SCYL3                                                                                         </td><td>SCYL3                                                                                         </td><td>SCY1 like pseudokinase 3 [Source:HGNC Symbol;Acc:HGNC:19285]                                  </td><td>ENSG, ARRAYEXPRESS                                                                            </td></tr>\n",
       "\t<tr><td>5                                                                                             </td><td>ENSG00000000460                                                                               </td><td>5.1                                                                                           </td><td>C1ORF112                                                                                      </td><td>C1ORF112                                                                                      </td><td>chromosome 1 open reading frame 112 [Source:HGNC Symbol;Acc:HGNC:25565]                       </td><td>ARRAYEXPRESS, ENSG                                                                            </td></tr>\n",
       "\t<tr><td>6                                                                                             </td><td>ENSG00000000938                                                                               </td><td>6.1                                                                                           </td><td>FGR                                                                                           </td><td>FGR                                                                                           </td><td>FGR proto-oncogene, Src family tyrosine kinase [Source:HGNC Symbol;Acc:HGNC:3697]             </td><td>ENSG, ARRAYEXPRESS                                                                            </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "\\begin{tabular}{r|lllllll}\n",
       " num & GENE\\_ID & co & SYMBOL & ALT & a & b\\\\\n",
       "\\hline\n",
       "\t 1                                                                                                  & ENSG00000000003                                                                                    & 1.1                                                                                                & TSPAN6                                                                                             & TSPAN6                                                                                             & tetraspanin 6 {[}Source:HGNC Symbol;Acc:HGNC:11858{]}                                              & ARRAYEXPRESS, ENSG                                                                                \\\\\n",
       "\t 2                                                                                                  & ENSG00000000005                                                                                    & 2.1                                                                                                & TNMD                                                                                               & TNMD                                                                                               & tenomodulin {[}Source:HGNC Symbol;Acc:HGNC:17757{]}                                                & ARRAYEXPRESS, ENSG                                                                                \\\\\n",
       "\t 3                                                                                                  & ENSG00000000419                                                                                    & 3.1                                                                                                & DPM1                                                                                               & DPM1                                                                                               & dolichyl-phosphate mannosyltransferase subunit 1, catalytic {[}Source:HGNC Symbol;Acc:HGNC:3005{]} & ARRAYEXPRESS, ENSG                                                                                \\\\\n",
       "\t 4                                                                                                  & ENSG00000000457                                                                                    & 4.1                                                                                                & SCYL3                                                                                              & SCYL3                                                                                              & SCY1 like pseudokinase 3 {[}Source:HGNC Symbol;Acc:HGNC:19285{]}                                   & ENSG, ARRAYEXPRESS                                                                                \\\\\n",
       "\t 5                                                                                                  & ENSG00000000460                                                                                    & 5.1                                                                                                & C1ORF112                                                                                           & C1ORF112                                                                                           & chromosome 1 open reading frame 112 {[}Source:HGNC Symbol;Acc:HGNC:25565{]}                        & ARRAYEXPRESS, ENSG                                                                                \\\\\n",
       "\t 6                                                                                                  & ENSG00000000938                                                                                    & 6.1                                                                                                & FGR                                                                                                & FGR                                                                                                & FGR proto-oncogene, Src family tyrosine kinase {[}Source:HGNC Symbol;Acc:HGNC:3697{]}              & ENSG, ARRAYEXPRESS                                                                                \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "num | GENE_ID | co | SYMBOL | ALT | a | b | \n",
       "|---|---|---|---|---|---|\n",
       "| 1                                                                                              | ENSG00000000003                                                                                | 1.1                                                                                            | TSPAN6                                                                                         | TSPAN6                                                                                         | tetraspanin 6 [Source:HGNC Symbol;Acc:HGNC:11858]                                              | ARRAYEXPRESS, ENSG                                                                             | \n",
       "| 2                                                                                              | ENSG00000000005                                                                                | 2.1                                                                                            | TNMD                                                                                           | TNMD                                                                                           | tenomodulin [Source:HGNC Symbol;Acc:HGNC:17757]                                                | ARRAYEXPRESS, ENSG                                                                             | \n",
       "| 3                                                                                              | ENSG00000000419                                                                                | 3.1                                                                                            | DPM1                                                                                           | DPM1                                                                                           | dolichyl-phosphate mannosyltransferase subunit 1, catalytic [Source:HGNC Symbol;Acc:HGNC:3005] | ARRAYEXPRESS, ENSG                                                                             | \n",
       "| 4                                                                                              | ENSG00000000457                                                                                | 4.1                                                                                            | SCYL3                                                                                          | SCYL3                                                                                          | SCY1 like pseudokinase 3 [Source:HGNC Symbol;Acc:HGNC:19285]                                   | ENSG, ARRAYEXPRESS                                                                             | \n",
       "| 5                                                                                              | ENSG00000000460                                                                                | 5.1                                                                                            | C1ORF112                                                                                       | C1ORF112                                                                                       | chromosome 1 open reading frame 112 [Source:HGNC Symbol;Acc:HGNC:25565]                        | ARRAYEXPRESS, ENSG                                                                             | \n",
       "| 6                                                                                              | ENSG00000000938                                                                                | 6.1                                                                                            | FGR                                                                                            | FGR                                                                                            | FGR proto-oncogene, Src family tyrosine kinase [Source:HGNC Symbol;Acc:HGNC:3697]              | ENSG, ARRAYEXPRESS                                                                             | \n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "  num GENE_ID         co  SYMBOL   ALT     \n",
       "1 1   ENSG00000000003 1.1 TSPAN6   TSPAN6  \n",
       "2 2   ENSG00000000005 2.1 TNMD     TNMD    \n",
       "3 3   ENSG00000000419 3.1 DPM1     DPM1    \n",
       "4 4   ENSG00000000457 4.1 SCYL3    SCYL3   \n",
       "5 5   ENSG00000000460 5.1 C1ORF112 C1ORF112\n",
       "6 6   ENSG00000000938 6.1 FGR      FGR     \n",
       "  a                                                                                             \n",
       "1 tetraspanin 6 [Source:HGNC Symbol;Acc:HGNC:11858]                                             \n",
       "2 tenomodulin [Source:HGNC Symbol;Acc:HGNC:17757]                                               \n",
       "3 dolichyl-phosphate mannosyltransferase subunit 1, catalytic [Source:HGNC Symbol;Acc:HGNC:3005]\n",
       "4 SCY1 like pseudokinase 3 [Source:HGNC Symbol;Acc:HGNC:19285]                                  \n",
       "5 chromosome 1 open reading frame 112 [Source:HGNC Symbol;Acc:HGNC:25565]                       \n",
       "6 FGR proto-oncogene, Src family tyrosine kinase [Source:HGNC Symbol;Acc:HGNC:3697]             \n",
       "  b                 \n",
       "1 ARRAYEXPRESS, ENSG\n",
       "2 ARRAYEXPRESS, ENSG\n",
       "3 ARRAYEXPRESS, ENSG\n",
       "4 ENSG, ARRAYEXPRESS\n",
       "5 ARRAYEXPRESS, ENSG\n",
       "6 ENSG, ARRAYEXPRESS"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Load ENSEMBL to Symbol dictionary\n",
    "gene_file <- file.path('data', 'raw', 'gprofiler_results_1002952837509.xlsx')\n",
    "gene_dict_df <- readxl::read_excel(gene_file, col_names = FALSE)\n",
    "colnames(gene_dict_df) <- c('num', 'GENE_ID', 'co', 'SYMBOL', 'ALT', 'a', 'b')\n",
    "gene_dict_df <- gene_dict_df[!duplicated(gene_dict_df$GENE_ID), ]\n",
    "head(gene_dict_df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Parsed with column specification:\n",
      "cols(\n",
      "  GENE_ID = col_character(),\n",
      "  wildtype = col_double(),\n",
      "  KRAS = col_double(),\n",
      "  NRAS = col_double()\n",
      ")\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<thead><tr><th scope=col>GENE_ID</th><th scope=col>wildtype</th><th scope=col>KRAS</th><th scope=col>NRAS</th><th scope=col>num</th><th scope=col>co</th><th scope=col>SYMBOL</th><th scope=col>ALT</th><th scope=col>a</th><th scope=col>b</th></tr></thead>\n",
       "<tbody>\n",
       "\t<tr><td>ENSG00000244405                                                                   </td><td>-1.25074291                                                                       </td><td>1.0092025                                                                         </td><td>0.3876048                                                                         </td><td>40217                                                                             </td><td>40217.1                                                                           </td><td>ETV5                                                                              </td><td>ETV5                                                                              </td><td>ETS variant 5 [Source:HGNC Symbol;Acc:HGNC:3494]                                  </td><td>ARRAYEXPRESS, ENSG                                                                </td></tr>\n",
       "\t<tr><td>ENSG00000153064                                                                   </td><td>-0.53002518                                                                       </td><td>0.8806176                                                                         </td><td>0.0000000                                                                         </td><td> 9726                                                                             </td><td> 9726.1                                                                           </td><td>BANK1                                                                             </td><td>BANK1                                                                             </td><td>B-cell scaffold protein with ankyrin repeats 1 [Source:HGNC Symbol;Acc:HGNC:18233]</td><td>ARRAYEXPRESS, ENSG                                                                </td></tr>\n",
       "\t<tr><td>ENSG00000131187                                                                   </td><td>-0.79863799                                                                       </td><td>0.6108468                                                                         </td><td>0.2938256                                                                         </td><td> 6476                                                                             </td><td> 6476.1                                                                           </td><td>F12                                                                               </td><td>F12                                                                               </td><td>coagulation factor XII [Source:HGNC Symbol;Acc:HGNC:3530]                         </td><td>ARRAYEXPRESS, ENSG                                                                </td></tr>\n",
       "\t<tr><td>ENSG00000184445                                                                   </td><td>-0.18836871                                                                       </td><td>0.4978895                                                                         </td><td>0.0000000                                                                         </td><td>16035                                                                             </td><td>16035.1                                                                           </td><td>KNTC1                                                                             </td><td>KNTC1                                                                             </td><td>kinetochore associated 1 [Source:HGNC Symbol;Acc:HGNC:17255]                      </td><td>ARRAYEXPRESS, ENSG                                                                </td></tr>\n",
       "\t<tr><td>ENSG00000143878                                                                   </td><td>-0.09074332                                                                       </td><td>0.4774349                                                                         </td><td>0.0000000                                                                         </td><td> 8573                                                                             </td><td> 8573.1                                                                           </td><td>RHOB                                                                              </td><td>RHOB                                                                              </td><td>ras homolog family member B [Source:HGNC Symbol;Acc:HGNC:668]                     </td><td>ENSG, ARRAYEXPRESS                                                                </td></tr>\n",
       "\t<tr><td>ENSG00000135046                                                                   </td><td>-0.22650054                                                                       </td><td>0.4408205                                                                         </td><td>0.0000000                                                                         </td><td> 7068                                                                             </td><td> 7068.1                                                                           </td><td>ANXA1                                                                             </td><td>ANXA1                                                                             </td><td>annexin A1 [Source:HGNC Symbol;Acc:HGNC:533]                                      </td><td>ARRAYEXPRESS, ENSG                                                                </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "\\begin{tabular}{r|llllllllll}\n",
       " GENE\\_ID & wildtype & KRAS & NRAS & num & co & SYMBOL & ALT & a & b\\\\\n",
       "\\hline\n",
       "\t ENSG00000244405                                                                        & -1.25074291                                                                            & 1.0092025                                                                              & 0.3876048                                                                              & 40217                                                                                  & 40217.1                                                                                & ETV5                                                                                   & ETV5                                                                                   & ETS variant 5 {[}Source:HGNC Symbol;Acc:HGNC:3494{]}                                   & ARRAYEXPRESS, ENSG                                                                    \\\\\n",
       "\t ENSG00000153064                                                                        & -0.53002518                                                                            & 0.8806176                                                                              & 0.0000000                                                                              &  9726                                                                                  &  9726.1                                                                                & BANK1                                                                                  & BANK1                                                                                  & B-cell scaffold protein with ankyrin repeats 1 {[}Source:HGNC Symbol;Acc:HGNC:18233{]} & ARRAYEXPRESS, ENSG                                                                    \\\\\n",
       "\t ENSG00000131187                                                                        & -0.79863799                                                                            & 0.6108468                                                                              & 0.2938256                                                                              &  6476                                                                                  &  6476.1                                                                                & F12                                                                                    & F12                                                                                    & coagulation factor XII {[}Source:HGNC Symbol;Acc:HGNC:3530{]}                          & ARRAYEXPRESS, ENSG                                                                    \\\\\n",
       "\t ENSG00000184445                                                                        & -0.18836871                                                                            & 0.4978895                                                                              & 0.0000000                                                                              & 16035                                                                                  & 16035.1                                                                                & KNTC1                                                                                  & KNTC1                                                                                  & kinetochore associated 1 {[}Source:HGNC Symbol;Acc:HGNC:17255{]}                       & ARRAYEXPRESS, ENSG                                                                    \\\\\n",
       "\t ENSG00000143878                                                                        & -0.09074332                                                                            & 0.4774349                                                                              & 0.0000000                                                                              &  8573                                                                                  &  8573.1                                                                                & RHOB                                                                                   & RHOB                                                                                   & ras homolog family member B {[}Source:HGNC Symbol;Acc:HGNC:668{]}                      & ENSG, ARRAYEXPRESS                                                                    \\\\\n",
       "\t ENSG00000135046                                                                        & -0.22650054                                                                            & 0.4408205                                                                              & 0.0000000                                                                              &  7068                                                                                  &  7068.1                                                                                & ANXA1                                                                                  & ANXA1                                                                                  & annexin A1 {[}Source:HGNC Symbol;Acc:HGNC:533{]}                                       & ARRAYEXPRESS, ENSG                                                                    \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "GENE_ID | wildtype | KRAS | NRAS | num | co | SYMBOL | ALT | a | b | \n",
       "|---|---|---|---|---|---|\n",
       "| ENSG00000244405                                                                    | -1.25074291                                                                        | 1.0092025                                                                          | 0.3876048                                                                          | 40217                                                                              | 40217.1                                                                            | ETV5                                                                               | ETV5                                                                               | ETS variant 5 [Source:HGNC Symbol;Acc:HGNC:3494]                                   | ARRAYEXPRESS, ENSG                                                                 | \n",
       "| ENSG00000153064                                                                    | -0.53002518                                                                        | 0.8806176                                                                          | 0.0000000                                                                          |  9726                                                                              |  9726.1                                                                            | BANK1                                                                              | BANK1                                                                              | B-cell scaffold protein with ankyrin repeats 1 [Source:HGNC Symbol;Acc:HGNC:18233] | ARRAYEXPRESS, ENSG                                                                 | \n",
       "| ENSG00000131187                                                                    | -0.79863799                                                                        | 0.6108468                                                                          | 0.2938256                                                                          |  6476                                                                              |  6476.1                                                                            | F12                                                                                | F12                                                                                | coagulation factor XII [Source:HGNC Symbol;Acc:HGNC:3530]                          | ARRAYEXPRESS, ENSG                                                                 | \n",
       "| ENSG00000184445                                                                    | -0.18836871                                                                        | 0.4978895                                                                          | 0.0000000                                                                          | 16035                                                                              | 16035.1                                                                            | KNTC1                                                                              | KNTC1                                                                              | kinetochore associated 1 [Source:HGNC Symbol;Acc:HGNC:17255]                       | ARRAYEXPRESS, ENSG                                                                 | \n",
       "| ENSG00000143878                                                                    | -0.09074332                                                                        | 0.4774349                                                                          | 0.0000000                                                                          |  8573                                                                              |  8573.1                                                                            | RHOB                                                                               | RHOB                                                                               | ras homolog family member B [Source:HGNC Symbol;Acc:HGNC:668]                      | ENSG, ARRAYEXPRESS                                                                 | \n",
       "| ENSG00000135046                                                                    | -0.22650054                                                                        | 0.4408205                                                                          | 0.0000000                                                                          |  7068                                                                              |  7068.1                                                                            | ANXA1                                                                              | ANXA1                                                                              | annexin A1 [Source:HGNC Symbol;Acc:HGNC:533]                                       | ARRAYEXPRESS, ENSG                                                                 | \n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "  GENE_ID         wildtype    KRAS      NRAS      num   co      SYMBOL ALT  \n",
       "1 ENSG00000244405 -1.25074291 1.0092025 0.3876048 40217 40217.1 ETV5   ETV5 \n",
       "2 ENSG00000153064 -0.53002518 0.8806176 0.0000000  9726  9726.1 BANK1  BANK1\n",
       "3 ENSG00000131187 -0.79863799 0.6108468 0.2938256  6476  6476.1 F12    F12  \n",
       "4 ENSG00000184445 -0.18836871 0.4978895 0.0000000 16035 16035.1 KNTC1  KNTC1\n",
       "5 ENSG00000143878 -0.09074332 0.4774349 0.0000000  8573  8573.1 RHOB   RHOB \n",
       "6 ENSG00000135046 -0.22650054 0.4408205 0.0000000  7068  7068.1 ANXA1  ANXA1\n",
       "  a                                                                                 \n",
       "1 ETS variant 5 [Source:HGNC Symbol;Acc:HGNC:3494]                                  \n",
       "2 B-cell scaffold protein with ankyrin repeats 1 [Source:HGNC Symbol;Acc:HGNC:18233]\n",
       "3 coagulation factor XII [Source:HGNC Symbol;Acc:HGNC:3530]                         \n",
       "4 kinetochore associated 1 [Source:HGNC Symbol;Acc:HGNC:17255]                      \n",
       "5 ras homolog family member B [Source:HGNC Symbol;Acc:HGNC:668]                     \n",
       "6 annexin A1 [Source:HGNC Symbol;Acc:HGNC:533]                                      \n",
       "  b                 \n",
       "1 ARRAYEXPRESS, ENSG\n",
       "2 ARRAYEXPRESS, ENSG\n",
       "3 ARRAYEXPRESS, ENSG\n",
       "4 ARRAYEXPRESS, ENSG\n",
       "5 ENSG, ARRAYEXPRESS\n",
       "6 ARRAYEXPRESS, ENSG"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Load Classifier Coefficients\n",
    "coef_file <- file.path('results', 'classifier', 'classifier_coefficients.tsv')\n",
    "coef_df <- readr::read_tsv(coef_file) %>%\n",
    "  dplyr::left_join(gene_dict_df, by = 'GENE_ID')\n",
    "head(coef_df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define cutoffs to label top 10 genes (by absolute value) for each classifier\n",
    "# Note that dplyr still has some serious issues with writing functions\n",
    "n <- 5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message:\n",
      "“package ‘bindrcpp’ was built under R version 3.4.3”"
     ]
    }
   ],
   "source": [
    "top_kras <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, KRAS) %>%\n",
    "  dplyr::select(KRAS) %>%\n",
    "  dplyr::arrange(desc(KRAS)) %>%\n",
    "  dplyr::filter(row_number() == n()))\n",
    "\n",
    "bot_kras <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, desc(KRAS)) %>%\n",
    "  dplyr::select(KRAS) %>%\n",
    "  dplyr::arrange(KRAS) %>%\n",
    "  dplyr::filter(row_number()==n()))\n",
    "\n",
    "top_wt <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, wildtype) %>%\n",
    "  dplyr::select(wildtype) %>%\n",
    "  dplyr::arrange(desc(wildtype)) %>%\n",
    "  dplyr::filter(row_number() == n()))\n",
    "\n",
    "bot_wt <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, desc(wildtype)) %>%\n",
    "  dplyr::select(wildtype) %>%\n",
    "  dplyr::arrange(wildtype) %>%\n",
    "  dplyr::filter(row_number()==n()))\n",
    "\n",
    "top_nras <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, NRAS) %>%\n",
    "  dplyr::select(NRAS) %>%\n",
    "  dplyr::arrange(desc(NRAS)) %>%\n",
    "  dplyr::filter(row_number() == n()))\n",
    "\n",
    "bot_nras <- as.numeric(coef_df %>%\n",
    "  dplyr::top_n(n = n, desc(NRAS)) %>%\n",
    "  dplyr::select(NRAS) %>%\n",
    "  dplyr::arrange(NRAS) %>%\n",
    "  dplyr::filter(row_number()==n()))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot gene coefficients scatter\n",
    "p <- ggplot2::ggplot(coef_df,\n",
    "                     aes(x = KRAS, y = NRAS, color = wildtype)) +\n",
    "  geom_point(alpha = 0.8, size = 0.1) +\n",
    "  scale_color_gradient2('Wildtype\\nGene Weight',\n",
    "                        low = \"blue\",\n",
    "                        mid = \"grey\",\n",
    "                        high = \"red\") +\n",
    "  xlab(\"KRAS - Gene Weight\") +\n",
    "  ylab(\"NRAS - Gene Weight\") +\n",
    "  geom_text_repel(data = subset(coef_df,\n",
    "                                (wildtype >= top_wt | wildtype <= bot_wt) |\n",
    "                                  (KRAS >= top_kras | KRAS <= bot_kras) |\n",
    "                                  (NRAS >= top_nras | NRAS <= bot_nras)\n",
    "  ),\n",
    "  arrow = arrow(length = unit(0.02, 'npc')),\n",
    "  segment.size = 0.3,\n",
    "  segment.alpha = 0.6,\n",
    "  box.padding = 0.17,\n",
    "  point.padding = 0.1,\n",
    "  size = 1.8,\n",
    "  fontface = 'italic',\n",
    "  aes(x = KRAS, y = NRAS, label = SYMBOL)) +\n",
    "  theme_bw() +\n",
    "  theme(axis.text = element_text(size = rel(0.5)),\n",
    "        axis.title = element_text(size = rel(0.6)),\n",
    "        axis.title.y = element_text(margin = \n",
    "                                      margin(t = 0, r = 0, b = 0, l = 0)),\n",
    "        axis.title.x = element_text(margin =\n",
    "                                      margin(t = 3, r = 0, b = 0, l = 0)),\n",
    "        legend.text = element_text(size = rel(0.4)),\n",
    "        legend.title = element_text(size = rel(0.5)),\n",
    "        legend.key = element_rect(size = 0.2),\n",
    "        legend.position = 'right',\n",
    "        legend.key.size = unit(0.4, 'lines'),\n",
    "        legend.margin = margin(l = -0.3, unit = 'cm'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {},
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAIAAAByhViMAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nOzdeXwc5Z3v+19V77tauyxZki3LK8bExOyLCYvZAkkmyZDcEJYkXHJNcjM5\nDJmZMDN5HSbJZMJMCDOBcxNOViBkTjIBD+uAASdhNQbbsvFuy5a1b93qvbuW+0d7FMfYkoGW\nSqr+vF+8eHU/Xd31q7JU+vZTVc+jmKYpAAAAmP1UqwsAAABAaRDsAAAAbIJgBwAAYBMEOwAA\nAJsg2AEAANgEwQ4AAMAmCHYAAAA2QbADAACwCafVBfzRv/zLv2zYsMHqKmY0wzBERFWJ45Mw\nDIO9NCnTNE3TVBRFURSra5nRiqO4s5cmZRgGP04ngwPUyZj+v3erVq268847p211U2cGBbuR\nkZG///u/X7lypdWFzFzDw8MOh6OiosLqQma6WCwWiUT4AzOxVCqVyWQqKiqczhl0HJiBstms\nYRh+v9/qQmY0XddHR0c9Hk8oFLK6lpluZGSksrLS6ipmung8XigUqqqqpuFIPjo6eumll7a1\ntU31iqYHXxoAAABsgmAHAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0A\nAIBNEOwAAABsgmAHAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBN\nEOwAAABsgmAHAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCacVhcA4P16/fl9K85udnud214/HAh7\n5i+p3fzSwflLasOVvvhIes/WvmWrmnwBd3Hh8Zbi09RYbuebnZEq/+IPzBGR/W8P9B+Om4a5\n5ION0epAMp7dvvFw8VVdN954Yb+ImKaZSuQu/tgpVm0vAOBECHbArJfLFBxOddOGA01tVfVz\nIyKSiGWCFV4RGepNzG2rig2nx4PdeEu4yiUiY6OZytpgMOIVEdMwdc04+7L2kYHk4X0j0erA\n0a86HOqZlywQkbc3dS/+QKNVGwsAmACnYoHZLZ/TTNN8+ZndC5bXF1OdYZgioqrK4f0jtU2R\nUNQbG0oVF35nS0NLRXIsG6nyiYiiKu2n1otI76FYQ0uFiMxpjY6/WtR9YNQfdFdU+6d3KwEA\nJ4UeO2B2GxvJOJ2OyrrgcF+iosovIslYNlTh0wp6585BX8BlmqLrhoi8o6Wi+AnpZN4f9BQf\n5zKFnW/1NLRUVNeH3vlqIa93Hxg540Nt072RAICTQ7ADZrf4SLplUfWc1spXn91TMyccjvri\nI+lIpW/31r7TzmsJhr0i8uqze0TknS0iYhimqirFx7HhdOfOgaUfbHJ7nO98VUQO7BiYt7hm\nOrcOAPCucCoWmN3iw+lIpV9RZMXZzVtePmjoZnwk7XCqmWSumOFExOVxDnSPHdOSSeW79sTf\neOFAMp5948X9IrLx+X25jLb5pYObXzooIru29L723N7xV0VksDdR3RC2YisBACeFHjtgdlt5\nwbziA3/Ic/5Vi0XklDPmikhze/X4MqdfME9EahvDR7ekUqm57ZGKigqn88hx4NJPLD/6kxet\naDhmXedevrD0GwAAKB167AAAAGyCYAcAAGATBDsAAACbINgBAADYBMEOKFO6ZqTieTGtrgMA\nUDrcFQuUqZ2b+gZ7ElpWbVtWb3UtAIDSoMcOKGOmmPTYAYCN0GMHlKmFH6irmetvaqmyuhAA\nQMnQYweUKZfbEYp6lKNmDAMAzHYEOwAAAJsg2AFlKpsy+zt1XeMiOwCwD66xA8rUq+uSh/fk\nc/HUyosrrK4FAFAa9NgBZcrpFjHF6eYaOwCwD3rsgDJ19jXBhWemm+YHrC4EAFAy9NgBZcrh\nUkJVqkKHHQDYCMEOAADAJgh2AAAANkGwAwAAsAmCHQAAgE0Q7AAbeuJ/D6YT+vjTJ//34Jvr\nx4qPe/fn7v8fh4rjEo/2Gc/+YiQZ08YX2/hMvPj4v34xlE0ZWzYkHru//9Ef9Pcfyk3vFgAA\n3guCHWBD6THdH3KMPzUNc7S/ICJiypYNYw1tHodTEVMGDuoLT/cPHMqPL9nbmTN0U0SyKcPt\nUwt549ov1p15ZcWuN1LTvhEAgHeNYAfYTTZteP3q0U89flXTTBHZ9Waqfp4nGHGIyP6tuYY2\nR2W9a6ArP/6uBSv8ezentbzpdCuqKh+8NCIi+zvS85f7LdoaAMC7wADFgN0MHc5XN7qPeTra\nX8gk9f1b08vOCopIPmfseDXjDYrLHdcLRxaranQvPD3w9E8Gw1XOqgaXiKTH9Neejs1f7m9q\n91q0NQCAd4FgB9jNYHe+5qhgN9idr53rdjiVDb8eOe3CcO+BXF2L+43/ip//8ZDLn6+oqHjy\ngZHxxZwupbLBteO1ZMtS30BXvuMPiXOviXoDdO0DwOxQXsfrjo6Ozs7O4uO33367UCh0dHRs\n3bp169atmzdv7urqKr6UTCa3b9+eyx25Wry/v7+jo2P79u1btmxJpbjSCDPd0OH8lt+NPf7D\ngcd/OJDLGEPd+epGd02TW0xpmO8Z7M473UpiRA9XHbkIz+tXx4a14mIisvy80I7XU9Vz3E/+\n74F0Ql//yPDzjwxbukEAgJNVdj128XjcNE1FUQqFgsvlMgxjxYoVIqJp2ltvvTV37lwRGR0d\nra+vHxsbq6mpEZGurq7TTz9dUZRsNlsoFCzeAGAyl15f/SdPP1MtIvWtnstvqhGRNZ+tFpE1\nN3jGv6Vc9tnq8cVEJBR1/j//3CwiN36jaRqrBgCUQBkFu2KSi0QiAwMD1dXVDoejUCik0+mO\njo58Pu/xeBYuXCgi/f39VVVVqqr29PQUg11NTc3OnTurq6urqqq8Xq40AgAAM1QZBbtkMhkM\nBuvq6rZt2+b1egOBQDKZnDt3blNT0+HDh51OZyQS0XW9u7vb4/GIiK4fGQaspaVF07ShoaHN\nmzefdtppqlpe568BAMBsUUYZpRjsVFUNBAJ9fX2hUCiZTPp8PhGpqakZGBgQkc7OzsWLFy9b\ntmzZsmWmaYrI4cOH8/m80+msr693OBz5fH6S1QAAAFikjHrsEolEQ0ODiDQ2Nr766qstLS3F\nc7Ii4vF4NE2Lx+PZbNbvPzJel8vlymaz4XB4x44dDofDMIxIJMKpWAAAMGOVUbBbunRp8YHX\n6129evXRLSKycuVKEYlEIscs7/V6i3dXAAAAzHBldCoWAADA3gh2AAAANkGwAwAAsAmCHQAA\ngE0Q7AAAAGyCYAcAAGATZRTsent7Ozo64vG41YUAAABMiTIKdkNDQ7FYbGRkxOpCAAAApkQZ\nDVDc1NQUCATq6+utLgQAAGBKlFGwi0aj0WjU6ioAAACmShmdigUAALA3gh0AAIBNEOyAMpXL\n5mOjCU3TrS4EAFAyZXSNHYCjHdjXFR9NKKK2tbdYXQsAoDTosQPKlNvtMk3T43FbXQgAoGTo\nsQPK1Pz25pr6yurqKqsLAQCUDD12QPlyOh1WlwAAKCWCHQAAgE0Q7AAAAGyCYAcAAGATBDsA\nAACb4K5YwDKZTHbvzs5CQTMMo7qmsrWtKZvJ7d3dqRU0wzCrqita5jeJyJZNb0cqwq1tTSIS\nGx3b9NrWiy47p+Otne9s3PrWTsMwVEVxOB2Ll7a53K79ew7FYwnDMOa2NNTWV1u8wQCAKUaw\nA6yh6/q2zbuWrVjo9/tE5MDernwu37F55ykrFvn8XhHZ+tbOsXgyHAmapplKpYvvOnSguyIa\nUVX1uI25bO7Mcz8gIoMDI7t3Hli6vD1cEZrf3jw6HO/vHyLYAYDtcSoWsMZA33BtfXUx1YnI\nvAVzh4djNXVVxVQnItHKSDKRKhQ0p8tpGIaI9PYMRKJhj9d9okaX21V8b01tZWIsqShKJBLc\n+MqWw4d62xe1WrCRAIDpRbADrJFOZ/wB35+0JDPBUGD8aTKRDIYCibFkKBx0OBz5fGGgbzgY\n9IfCwQkai+/N5/JOp1NEXG7XqrNXiEgmnZ3GjQMAWINgB1gjEPTHRseKjwf6h3t7BoIhf2w0\nXmyJxxLZTC4cCSbGkuFwMBQO7ty+t2VeY2IsFQ6fuDESFBHDMHa9vb+puWH3jv35XF5EVAe/\n6QBQFrjGDrBGXX1NbGTszY3bRCQQ8LUvnqco6ujI2Jsbt6mK4va4TzltsYgkxlKNTfWqqowO\nxyqi4a6DPU3N9d2H+47T2NU32D/Sc7hf1/U5jXX1c2r8AV/H5l1Ol9Pv9x7dFwgAsCuCHWAN\nRZHFy9qOaXxnyykrFolIJBo+NRoWkeWnLT5h42mLjnlvOBI8/czlU1I9AGBGmkHBzjCMXC6X\nyWSsLmRGMwyDXTQpwzCyWS4pm4SmaSKSy+UKhYLVtcxomqbxezep4q08uq6zoyZlmiZ7aVLF\nn6hMJqMoylSvy2Z/L2ZQsBMRRVGm4Z9wtmMXTYofpJPHvjoZ7KVJje8fdtTJYC+dpOn51bPZ\nP8cMCnaqqrrdbq/Xa3UhM1cqlVJVlV00qWw26/F4bPa7WnK6rhcKBbfbXbx/9rg63x7rfHvM\n4VTzWf2Dl9RGqj0bftNd3ehddlaViAwczjz38KE//x8Ls0lt8+8GCzlD18yGeYGlZ1aKSDJe\neOuFwXxW1zWzcUFg2VlVY8P5Tc8PVNZ5VlxQM33bWQqGYfB7NzFd19PptMPhYEdNKp1Os5cm\nlcvldF33er3TcCT3eDxTvYrpNIOCHYCZZsfrI5ff0KookooXchldREzTHBvOF1/duXGkdq7f\nNMyX1vWc8+E5wQqXiPz+0e7h3myk2v3Sup5zrm4IRd0isu3l4ULeGB3M1TT5wpVuC7cIAOyN\nQRAAnNDchaFXHu89tCvhDTor6735jO72OnTNFJED28dqGn3+kLNrd7KxPVhMdSJS1+yPDeYO\n7Uo0LwoVU52InHJOlcuttiwOjQ3nK+voqwCAqUKwA3BCp5xT9cFLaws5Y/3Dh3TNHOnPRms9\nTreaTetduxKRak+01jM2nI/W/DGrjfbnKmo8ydHCcXvmUvHCeAQEAJQcwQ7A8e3aNJpN626v\no+3UiMOlZtPaaH+ust5bWefZ+Ezf4jMqRweylfXeihrPwOEjs9YOdWdSiUJVgzdS7R44fOS+\nv67diQPbx0TE0E3VwYWPADCFuMYOwPFVNXhfWtfjdKmGYdbO9QXCrpGBbNuKiOpQ+g6ma5t8\nuzeNtp9W4fI4BrrSz/+qS1EVX8B53ofniEjz4vBAV//zv+oSkXCl+wMX1Xb8YWjgcGZsOP/7\nR7vP/0ij1RsHAPZEsANwfNVzfBdfN/folnM/PEdEahp9NR9tFJHzrp1TbF91Wd0x71WUYxuX\nn1c9hbUCAESEU7EAAAC2QbADypShm6m4bnUVAIBSItgBZWrzC7GX/iO2b+uY1YUAAEqGYAeU\nqUJOVxQp5AyrCwEAlAw3TwBl6tQLo0N9yXmLIlYXAgAoGXrsgDLlCzqqG90OJwPLAYB9EOwA\nAABsgmAHAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABs\ngmAHAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAH\nAABgEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAHAABg\nEwQ7AAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAHAABgEwQ7\nAAAAmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAHAABgEwQ7AAAA\nmyDYAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAHAABgEwQ7AAAAmyDY\nAQAA2ATBDgAAwCYIdgAAADZBsAMAALAJgh0AAIBNEOwAAABsgmAHAABgEwQ7AAAAmyDYAWXK\n0I1sumB1FQCAUnJaXQAAa+zd3jPYF9fy0rqgwepaAAClQY8dUKZ03RARXTOsLgQAUDL02AFl\nasGyOSND8aaWWqsLAQCUDD12QJlyuZ3hqF9VFasLAQCUDMEOAADAJgh2AAAANkGwAwAAsAmC\nHQAAgE0Q7AAAAGyCYAcAAGATBDsAAACbINgBAADYBMEOAADAJgh2AAAANkGwAwAAsAmCHQAA\ngE0Q7AAAAGyCYAeUL62gW10CAKCUCHZAmdq/u2v75n293YNWFwIAKBmCHVCmstm8mGYum7e6\nEABAyTitLgCANVI9wdigWlfpt7oQAJjItddKoSCqKiJyxhlyww2ydq2MjEh/vyxZIjU1UijI\ngw+KiJimXHSRPPSQrF0rmnZk+b/7OyuLn34EO6Bcma5cwq0qdNsDmNHicXnxxT9pefxx+dnP\nJJGQ224T05Szzz7S/rOfyeWXS0ODZLPy9NPTXujMQLADytRpH4rO7XM1t0esLgQATmhkRPbt\nk498RETE45Ff/epI+1tvySc+ISKiKBIKSSolpik/+Yk8+6xs3iw9PfKpT0llpXz/++Iss6TD\nl3WgTLk8SqTWqShW1wEAJ/bWW7J2rTz6qDz66B9TnYhs2SIrVhx5fPrpsnmzfOtb8ld/JW63\nNDfLunXyy19KICBPPmlJ1VYi2AEAgBlq82ZZvPjYRtOUVEqCwSNPzzhDnnpK9u6VK64QEXnm\nGZkzR0Qkn5fKymmsdWYoZQdlLpe75557RETX9dtuuy0cDh/96tatW7/1rW898sgjJVwjAACw\nsbfekt/+Vh54QERkzhz54Q9FRPbvl/nz/7jMmWfKTTfJm28eeVpZKVdfLeGwtLbKeedNe8VW\nK2Wwe/bZZ5cvX37llVe+8sorjz322PXXXz/+UjqdXrduXVNTUwlXBwAA7K14u+sx2trk6G6i\nxkaJx//49IorjnTdladSBrv9+/dffvnlItLW1vbMM88c/dJPfvKT66+//r777ju6cffu3X/5\nl385/rSmpiaRSIyOjpawJPvRdZ1dNCnDMGKxmNVVzHSGYYhIIpGwupCZzjRNEcnlclYXMgvk\n83kOUJMyTZO9NCld10Vkeo7kNvt7MVXX2Hk8nvHHr776an19fUtLyxStCwAAAFLaHrsFCxbs\n2bNn4cKFu3btamtrG29/5ZVXROSee+7p6el58MEHP/OZzxTbFy5c+Nhjj40vduedd4ZCoWg0\nWsKSbGZ4eNjhcFRUVFhdyEwXi8UikYjCDZ8TSqVSmUwmFAo5y20wgHcpm80ahuH3M5LzRIon\nE9xudygUsrqWmW5kZIS/dJOKx+OFQqGiooIj+btVygP6JZdc8v3vf7+jo0NEvvSlLx0+fPjH\nP/7x3/3d3/3FX/xFcYGvfe1r46kOAAAApVXKYOd2u4++Zi4QCPzdn07k8Z3vfKeEqwPwfgwM\nDIyOjqqqWlmG4wEAgE0xjh1QpkZHR7ldCQBshmAHlKm6urqqqqq6ujqrCwGAE9qxQy67TNau\nFV23upRZgmAHlKmKior6+nqv12t1IQBwQps2yWuvyU9/+icj1WEC3A0HAABmqA9/WAYGpL29\nHCcHe28IdgAAYIaKROSrX7W6iFmFU7EAAAA2QbADAACwCYIdAACATRDsAAAAbIJgBwAAYBME\nOwAAAJsg2AEAANgEwQ4AAMAmCHYAAAA2QbADAACwCYIdAACATRDsAAAAbIJgBwAAYBMEOwAA\nAJsg2AEAANgEwQ4AAMAmCHYAAAA2QbADgBnKNOX3/zH29E9jmaRhdS0AZgeCHQDMUKm4/ub6\n9O5N6e49eatrATA7OK0uAABwfMEKx4WfCCVjevMSt9W1AJgdCHYAMHOdtjpgdQkAZhOCHQC8\nL/EhfcOv49mUISIN893nfzQsIo/+YFjXRFFFK5hty72nXxp89AfDl10f9QXVF34Vr6h1rLw4\n2NeZf+XxRCFvnnZhYOHpPqu3A4AdEOyAMjUwMDAyMqIoSlVVldW1zGKFnPnEAyNX3hytqHWK\nyOtPJzJJwxdUtYL58a9Ui0g+Yz5y9+DplwZTY4bTrTz+o5GlZ/nbVnh1zfzDo4mPrK1UFBnq\n0azeDgA2QbADytTo6GgqlYrFYgS792Pn6+klZ/qLqU5Ezrg8JCLZlDHapz36g+FkzAhVOi75\ndEU2ZZiG/Ps/D136mYq6FpeIHNqZc3uUF34VVxT50HURK7cBgI1wVyxQpurr66uqqurr660u\nZHYbHdAq65wi0rU79+t7hn/2PwdScX3gUOH0S4MfWVu19Cxf2wrvnDb3wKGCy6O0LvN07c4V\n39h/sNDY7r70MxUev3pge87SjQBgHwQ7oExFIpG6ujqPx2N1IbNbZb3r4I6siMxd6PnobZVe\nnxqIOAa6CtE6p4i0r/Tt2pgRkYGuwooLA+deGz7QkR3qLoiIy6NUNThFJJ0wQlGHpRsBwD4I\ndgDw3i0726/l5f98b+i3/za87n+NnHFFUET6Dx0JdqGoI5cxMkmj/1Chdq5LUeSy6yv+6xcx\nXTOXnuV/64XUY/ePBCNq7VyX1dsBwCa4xg4A3jtFlQ996tgr5K76fHT88af/qubolkiNs9ji\nCyofvY2rGwGUGD12QJnqPTS8d1tvbDhpdSEAgJIh2AFlaqhvbGwkM9wft7oQAEDJcCoWKFNN\n86qHg2MNLdVWFwIAKBmCHVCmojUhb9DpD3BXLOzpl9+N65qoqojIkjM8cxe5/usXyfpW54f+\nPCAirz6Z2fNW3tDNSz4dbFzAn0LYBz/NAAAbig8at/7TH+9i2f5qbu4iV/Uch4gYuhRy5vVf\njxzaVdj6+2zjgqB1ZQIlxjV2ADBTFPKm1SXYRCZpegPK0S3LzvIM9+j1rU4RUR1y/kf9IrJz\nY37xGXRaw1bosQOAGaHjpey6/y950ScD51zts7qWWa+3szA6oD/4rbiI+ELKn30pLCKxQT1a\ne2Qs6GTMeOH/pBav8sxbxiCCsBWCHQDMCANduqJKX6dmdSF20Nepr/5E4AOrveMtuiYO55E+\nvJ792sZnMpd+JuAPcdoKdkOwA4AZ4Yw1vki1Ov8Ut9WF2EHfgcJ5H/WPP93w6/ShXYWhHu3/\n3DP28f83/Mjd8bpm57r/lfSHlWv+75CFdQIlR7ADgBkhFFU/eAknYUvjY18KH/30wo/7j376\n1fuY8wO2RS80UKb2bOve8tKhnq4hqwsBAJQMwQ6YEfL5/PDwsKZN3/VVfYcS6ZjZfWBk2tYI\nAJhqBDtgRjhw4MD27dsPHz48bWtsnF8RrHS0tNdO2xoBAFONa+xQYplMpqurKxKJ1NXVWV3L\nbOJyuUTE6Zy+X8n5i+trm4JVVRXTtkYAwFQj2KHEhoaG+vv7U6kUwe5dmT9//pw5c7xe7+SL\nlsihQ4cGBgYWLlzY0NAwbSsFAEwpgh1KrKqqKpVKVVTQD/SuTWeqE5FcLici2Wx2OlcKTKl0\nwmBoOpQ5gh1KzO/3L1682OoqMLmmpqZAIEB33fvU8VJu28tZl1vJJM3LbwjWNDnG557XNbnq\nc8HKescj3x3TNVNRRctL+0r32Vf5fvnduKGJooqINLa5Lvy4f/xdHp9yxU1Bf0hllvp369Wn\n0s/9Mn3xnwfOvopRY1C+OFgAZcrn88n0XtVnS688nv78N6OqKrFBPZ0wRSQ+ZNz6naiIHNhW\nePWpzJU3BQt58/qvR0QklzF/+o3Y2Vf5cmnzxr//k17txIhxy7ejIrL7zfx/PZi65pYQs9S/\nW/EhQ1VlpF+3uhDAShzTgTKl63o2mzVNZp1/Xxav8jx2f2LxKnf7B9wVNUomafr+e+75oR6t\nolrNJM3hXv2R744lRvVwlePKzwUzSXOkX3/ku2Mi4nDJJ74SziTN8ROI7R9wP/9Iilnq34Nz\nr/E3trlalzL3K8oawQ4oUwcPHhwZGTFNc+7cuVbXMotd8DF/JmXu2pj7+T/EP3tnpLezEB8y\nHvnumChSWad+6LrAoV2Fs670nXWl77WnMh6fMneha/+2/BmX+c77yB/nQujtLNTPO3I0TsUN\nj18RZql/94IR9ZRzCMEodwQ7oEwZhjH+f7w3rz+dWXa2JxBRT1vt7Xgpl4qbfZ36hR/3r7jg\nj/fB9HXqNU0OEVlyhuc/f5g4bbW3r1OvbnQc/Tl9nfqceU4R0Qrm0z9LrrrMxyz1AN4bgh3+\nqLu72zCMxsZGVeVvif21trbG4/H6+nqrC5nF5ixw/se/JdxeRdfM5sWuSLXad6Bw/scCRy/T\nd6CweJVbRMJVajZtphNG34HCzteNN5/Pikgoqn74llDfgcLuTcbmDdl81jztQu+ysz3fWzvM\nLPUA3gOCHY5IpVL79u0TkWAwGI1GrS4HU667c3hkcMztCNQ3MiH6e9S0wFW8K2LcMXPPH9Py\nubsqJl2miFnqAbw3dMzgCJ/PV1NTU11dHQxy/11ZGB0aSyeyI0NxqwsBAJQMPXY4QlXVJUuW\nWF0Fpo8roDs13eHRrC4EAFAy9NgBNhePx/fu3Ts2NnZMe1VNRSjqrqqutKQqAMBUoMcOsLme\nnp6hoSFN08LhP7mQK5eSfNJRyDKOHQDYBz12gM3V1NRUVlbW1NQc0z4yOpQzEoODg5ZUBQCY\nCvTYATZXXV1dXV39znav36GlxBdg8FsAsA967IAy5XdXOIxAKMg1dgBgH/TY2VlHR4dhGIqi\nmKa5cOFCn8+3c+fOpqamAwcOFNtFpKamJhKJ7N27NxQKzZs3z+qSMX3iw7l0zIwNpernMGQa\nANgEwc7OstnsqlWrRGR0dPTw4cPt7e3pdDoQCIy3Fw0MDEQiEb/ff+JPgg3Nbase6o83tpDq\nAMA+OBVrW4VCweU6cvlUOp32eDymaYqIpmlO558E+tra2nQ6HQoxbVF5iVYHm+ZX+gJMmg4A\n9kGPnW0lk8lsNtvR0SEiPp9v/vz5qVTK7/cX27ds2SIiLpdr6dKlIpLNZr1e7ySfCHsxDCOX\nyxXjPgDAHgh2tpVMJltbW4+e4j2ZTAaDwWJ7Q0PDeLthGKpK323Z2bt3bzwez7fADMUAACAA\nSURBVOfzXFsJALbBn3PbSiQSxwxIm0gkQqHQMe2dnZ0dHR3pdHr79u3TXiOslEqlTNOMx5kr\nFgDsgx472yqeYz1ae3u7iEQikaMbW1tbp60kzChVVVWjo6NH990CAGY7gh1Qppqbm6urq6PR\nqNWFAABKhlOxQPkqjmUIALANgh0AAIBNEOwAAABsgmAHAABgEwQ7AAAAmyDYAQAA2ATBDrCJ\nbDbb19dXKBROfvl4PK7r+pRWBQCYToxjB9jE/v37h4eHs9nsSQ46ffDgwUQiISJtbW1TWxkA\nYLrQYwfYhN/vNwzD5/Od5PK5XM4wjEwmM6VVAQCmEz12gE20trY2NTU5nSf7Sx0Oh3O5HDNP\nAICdEOwA+zj5VCciLS0tlZWVNTU1U1cPAGCacSoWKFOKorhcLqurAACUEsEOKFOmaRby3BIL\nALZCsAPK1O7NA2+s7+4+MGp1IQCAkiHYAWVqdGzQW50c6BuyuhAAQMlw8wRQplRvxqFK1hix\nuhAAQMnQYweUKdNQDF1UBwcBALAPjulAuTIdRkEVg4MAANgHx3SgTEUq/A6XNDTVWl0IAKBk\nuMYOKFPzF8wbHh6eO7fJ6kIAACVDsAPK1MGDB8fGxlwuV0tLi9W1AABKg1OxQJlSVdU0TVXl\nIAAA9kGPHVCmWltbq6qq6urqrC4EAFAyfFkHypTD4fB4PFZXAQAoJYIdys7o6OjmzZt7enqs\nLgQAgBIj2KHsxGKxRCIxPDxsdSEWi8VifX19mUzG6kIAACUzg66xMwwjn89ns1mrC5nRDMNg\nF03KNM1cLneiV6PRqK7rkUikzPdkb29vKpVyu90ul8vqWma0QqFgmmaZ/7RMyjAMEdF1nR01\nKX6cTkbxJyqbzSqKMtXrmuDvxWw0g4KdiJimaZqm1VXMdOyiSU38g+TxeBobG6Xs92Q0GjVN\nMxKJlPl+OBkcmiY1vn/YUSeDvXSSpmdH2eyfYwYFO1VVPR6Pz+ezupCZK51Oq6rKLppULpfz\ner3T8D1vVstkMul0ulAo8BM1sWw2axgGe2liuq5nMhmHw8GOmlQmk2EvTSqfz+u67vP5puFI\nbrMOVK6xA8pUKpUyDCORSFhdCACgZAh2QJlyOBwOh4ML7ADATmbQqVgA06mxsXFkZKS+vt7q\nQgAAJUOPHVCmhoeHh4aGxsbGrC4EAFAyBDugTMXjcdM0BwYGrC4EAFAyBDugTAWDQUVRKisr\nrS4EAFAyXGMHlKm2traxsbHa2lqrCwEAlAw9dkCZSo8Zo91mIW9YXQgAoGTosQPKVMfvRwcO\nZ/Wc+9TzaqyuBQBQGvTYAWXKEzTc4Yw3yEEAAOyDHjugTJnhw1GfljJMkSqrawEAlAZf1oEy\nZRiGoki+kLe6EABAyRDsgDLV3NwcDAYXL15sdSEAgJLhVCxQpmpra0OhUCAQsLoQAEDJ0GMH\nlKmDnYd27dwzMDBodSEAgJIh2AFlamhoKJfN9fX2WV0IAKBkCHZAueju7u7s7NR1vfi0qrrK\n43XPmdNgbVUAgBLiGjugLGSz2X379olIKBSqqqoSkVwuqxtaJpuxujQAQMnQYweUBY/HM2fO\nnOrq6nA4XGxJJBK6rg8Oco0dANgHPXZAWVAUZcGCBUe3qKpqGIbH47GqJABAydFjB5Qpj8ej\nKIrX67W6EABAydBjB5SpuXPnjo6ONjY2Wl0IAKBkCHZAmQqHww6Hkx47ALATTsUCs9vOnTu3\nbduWybzrm1v3bRvc/mp/X1dsKqoCAFiCYAfMYtlsdmBgYGRkJJlMvtv3JuPZTCqfjOemojAA\ngCU4FQvMYl6vd+HChYVCobKy8t2+t6HN7w7n6pqDU1EYAMASBDtgdquvr39vbxyND6XzyZHR\noWhlpLQlAQCsQrADylQ4HNZ1fXy8YgCADRDsgDLV0NBQUVFRUVFhdSEAgJLh5gkAAACbINgB\nZcowjEwmY5rmNK83Fc/LdK8TAMoFp2KBMrV///5YLGaaZnNz87StdO+W0c2/65u3tOL0ixum\nbaUAUD7osQMwfXJZTVGUfE63uhAAsCd67DA7aJrmcDgURbG6EPuYN29eLBZraJjWnrOFH6iM\n1voq65jHDACmBMEOs0Aymdy1a5fT6TzllFMcDofV5diEpmmapum67nK5pm2lLrdjzjyGRAaA\nqUKwwyyQz+fT6bRhGMV+O6vLsYlDhw7F43FFUebNm2d1LQCA0iDYYRaIRqNLly51u90ej8fq\nWuzD7XYburBLAcBOCHaYBRRFqaqqsroKu0n3h7Ldqtbos7oQAEDJcFcsUKaSo4V0whgbyVld\nCACgZOixA8pUw2Jx16Tr5tdZXQgAoGTosQPK1OBIT15P9vZ3W10IAKBkCHZAmdI0TURSqZTV\nhQAASoZgB5Qv0zQZ8xkA7IRgB5Sp5ubmQCDQ3t5udSEAgJLh5gmgTAUCgVwu5/UyuxcA2Ac9\ndkCZ6u7uHhgY6Ovrs7oQAEDJ0GMHlKlQKKRpWjgctroQAEDJTN5jd8MNNxz99Pzzz5+yYgBM\nH8MwTNM0DMPqQgAAJTNRj93DDz/8z//8z3v27NmyZUuxJZfLqSpnbwE7SCZTmXQ2mUzW1TFG\nMQDYxETB7tOf/vS111576623fv3rXy+2qKra3Nw8LYUBmFqm5tBzDsXkegwAsI9JjumBQOC+\n++578sknY7GYaZrFxltvvXXqCwMwtTRNF9MoFHSrCwEAlMzkX9avuuqqqqqqlpaWaagGwLRJ\np1OqRx/sH1q0ZL7VtQAASmPyYJfNZn/7299OQykAppMphiKKphesLgQAUDKT3wmxaNGigYGB\naSgFKB+maY6OjiaTSQtrCIZ9Ikpjc62FNQAASmuiHrtLLrlERHRdX7JkyfLly53OIws/99xz\n01EaYF+JRGLbtm1+v3/58uVut9uSGqqi1aY2XBmtsmTtAICpMFGw+8Y3vjFdZQDlxel0mqbp\ncDgcDodVNfT19+YL+UOHuqqqK62qAQBQWhMFu/POO09EfvKTnyiK8sc3OJ0ej2fVqlVTXhpg\nX36//6yzznr/wc40zaN/Pd/tuxVFURiYEgBsZPKbJzZs2PDMM89cdNFFqqq+8MIL559//g9+\n8IMPfehD3/zmN6ehPsCu3v8Z2K6urqGhodbW1mg0+h7eHgwF4/F4NFrxPssAAMwck39bHxsb\n27Fjx8MPP/zggw9u375dUZSXXnqJy+wAyxVvv4jH4+/t7aqqFk8Hl7YqAICFJg92PT09FRVH\nvtNXVFQcPHhQVdVCgSESAIu1tra2trY2NDS8t7cnEgkRGR4eLmlROI7/HtwdAKbc5KdiFy5c\nePHFF5933nmKorz22mv19fUPPfTQ8uXLp6E4ABMIh8PhcPg9vz0QCMTj8fGvbZgi/QezG58Z\nntPmW3kxN6kAmHKTB7sf//jHTzzxxNatW0XkpptuuvbaaxOJxCc/+cmprw3AFGpraxsbG6ur\nq7O6EJsb6s6N9OYVVXl/d7oAwEmZKNhdd911jzzyyEc+8pHxltdee+3nP//5448/PvWFAZha\nyXh2uC8dDuVD4cm/4OE9a1sRFJGaJg+pDrCHdevW5fP5ZDKZTCZvu+22YuOvf/3rnTt33nnn\nndbWJhMHuy9+8Ysi8pWvfGW6igEwffbtPpTLZxRVlq9st7oWO/MGHMvOiVhdBYB3Z9WqVY8/\n/nhdXd3LL7/853/+511dXSJy3333JRKJr33tayLy05/+9J3v2r1798MPP2ztMMATBbsLL7xQ\nRM4///yf//znhw4duuuuuzZt2nTKKadMV20AppCmpFSPkcqNWF0IAMw411xzzZNPPnnTTTc9\n+eSTy5Yt27Jly4oVK5588sm77777gQceyGazwWBQRHp6em644YbFixcnEokFCxZ885vf3L59\n+8UXX/y1r33tscceq6mpufvuu6urqzVNe/LJJxctWrR9+/Z//Md/bGlpufnmm0OhUCqVuvfe\ne2tqakpY+eR3xa5du/bNN9986qmnRGT9+vW33HJLCVcPwCqm6CJmJpu1uhAAmHGuvfbaJ554\nQkRef/31v/zLv3ziiSey2WxPT8/ixYuPXuynP/3pZz/72X/9139ta2sTkU984hNXX331+eef\nf/311z/44IMi8vjjj3/84x93Op3Nzc3f/va3b7/99nvvvfdHP/rRhRde+MADD3zqU5/6wQ9+\nUNrKJw92Bw4cuP/++4vJ9I477jh48GBpKwBgKYbiAIBjnXrqqXv37u3s7Kyvr7/gggt+97vf\nvfjii5deeukxi3V3d8+dO1dE5s+ff3T7pz/96UceeWT//v3Nzc3FBNXS0iIiDQ0Nvb29+/bt\nW7du3Y033vjLX/4yW+pv15NfNK1p2vjjTCaTSqVKWwEAC6kqV/QDwHFcdNFFd9111xVXXOFy\nuUKh0K9+9asvfOELxyxTV1dXvPxuz549TqdTURTDMEQkEoksWLDgzjvvvPnmm4tLHjhwQEQ6\nOzubmprmz5+/atWqz372syMjIyUfGHjyHrs1a9ZcddVVnZ2dt99++8qVK//sz/6stBUAeKdc\nLnf0d6opYpqmqjJZLAAcx7XXXvvQQw9ddtllInLZZZc999xzZ5111jHL3HjjjT/+8Y+/8IUv\n7Nu3zzTN9vb23/zmN8XBQ26++ebf//73H/rQh4pL9vX1ffnLX/7Wt771pS996XOf+9xTTz11\n8803f+Yzn+nv7y9t2RP12L3xxhunn3763/zN37z44osbNmzwer0PPPDAueeeW9oKABwjlUpt\n2rQpGAyecsop739K2RMJBALpdLq+vn6KPh8AZrXVq1ePnyf9whe+MN5d9/nPf/7oxV544YWj\nn27fvr34IJ/P33zzzeNfni+44ILxsVFE5Je//OUUlT1RsPv0pz8di8Uuu+yyNWvW3HrrrQxk\nitnIMIx0Ou33+2dR15Su6yKSTCZN5qICgFno7rvvfu655/793/99+lc9UbDbvXv3oUOH1q9f\n//TTT3/ta1+rr6+//PLL16xZUxwGBZgVOjs7Dx8+3NjYWLxlaVYIh8MrVqxwOp0ej2fq1pJM\npEWU/r6B9nbGsQOAUrr99ttvv/328ac33njjtK16kj6M5ubmm2666aGHHurp6bnrrrtefvnl\n1atXT0thQGkUO71mXddXJBIJBAJTuw5TURSjkDOmdi0AgJOzbt26X//618c0Pvjgg/fcc8/J\nf8gkwa6/v//hhx+++eab29vbv/Od76xevXrDhg3vulLAOq2trStWrDjmRnSIiNPhNE0lEmVS\nBAA4jt/85jeXXHLJ1VdfvXr16u985zvv+XNWrVpVvEPi5ZdfLo6NIiL33XffOz/zmmuu+fjH\nPz7xp+3evXvimS0mOhV76qmnJhKJT33qU9ddd92//du/+f3+SasHZhqHwxGJkF2OI58RLeNN\nKgxgBADH2rp16/e+971nnnkmEAgYhvEP//APAwMDgUDg6BkjnnjiiWeeeWbRokUbN278p3/6\np9bW1uPOJzHBJBapVOrotzz22GPZbPajH/3oTTfd1N7ens/nm5ub582b99JLL7399tsdHR0/\n/OEP77777uLkFueff/5xK3dMkPsKhUI6nd6wYUN3d/fY2FgkEqmqqpqK3Vf0/PPPL1mypKGh\nYepWMdtlMhlVVb1er9WFzHTZbNbr9TLn+sQO7D2kuDSX09Uyr8nqWmY0TdNM03S5XFYXMqOZ\nppnNZqf6wlB7yGQyPp/P6ipmulwuZxiG3++fhiN5Npv9xS9+0dLSsmbNmvHGn//852eeeebZ\nZ58tIoqiXHjhhYFA4L777mtra7vrrrvcbve6deuqq6tTqdS3v/1th8OxcePGLVu2HP3q+KVr\n0Wj0gQce+OQnP/mP//iPd9xxxwsvvLBq1ap77rnn61//+jEfGA6HNU177bXXVq9e/Vd/9Vev\nvvqq2+2uqKiIxWL333+/3+/fuHHjRRddFAwGb7jhhhNtzkQ9dl/+8pe//OUv67q+cePG9evX\nf/GLX+zt7T3jjDOOO/EtgNlF9WiqYqpu3epCAGDGURSleHH2zp07v/KVrySTyZtuumnfvn17\n9ux5/fXXc7lccSaJ4qlVr9ebyWQOHz58zKtFx0xi8d3vfnflypXFSSze+YEi0t3dXbxLtfhG\nEVmwYIGI+P3+XC43aeWTzzzhcDiWLFnS398/MDAwNDT0hz/84d3vHwAzjtPl0DTN56MDGACO\ndcUVV9x8883XXXfd4sWLn3766R/96Eejo6PHzBjx7LPPHv2WCeaTONEkFse85T//8z9FpLa2\ntru7W0S2bt36zqvgxie3OJGJbp544YUXvv71r5911lktLS0PPPDAwoUL/+M//qMYHgHMdoqi\nisjUDYCM96C3M93bmba6CgCybNmyO++882Mf+9jVV1999dVXb9y48cYbb5x4xogJXj3RJBbH\nfcsNN9zwwx/+cO3atb29ve88E3305BbHpUwwDMRpp522Zs2aNWvWnHfeedNw9C/uwZUrV071\nimav4eFhh8NRUVFhdSEzXSwWi0QiXGM3sRdeeFHEdCjOC1Yf/wpcFGWz2eK1PlO9ovhQ/omf\nHlREuez/aqpqmGU9qbquj46OejyeUChkdS0z3cjISGVlpdVVzHTxeLxQKFRVVU3DkXx0dPTS\nSy+94IIL/uVf/mWq13Uyenp6BgcHV6xYce+99/r9/mMmupjURKdiN2/e/P5qAzDTmTLLRviz\nMY9Praz1iChev8PqWgBYxuFwfPWrX21qakokEu/hrobJr7EDYEuKIoZhqk76NWcKb8B5yXVN\noojTNWumvwNQcnV1devXr3/PbyfYAeVLURSHg86hGcTpJtIBeF8mP4jkcrkf/ehHf/u3fysi\nmzZtOplbbQHMfMUrV7h5AgDsZPJgt3bt2jfffPOpp54SkfXr199yyy1TXxUAAADetcmD3YED\nB+6///5gMCgid9xxx8GDB6e+KgAlMME97yIipmKaYk40HBIAYJaZ/Bo7TdPGH2cymVSKmSWB\nWWDPnj3JZHLevHknGh/HMAwRM5vNTnNhADDDGYm4MTZyolcd9XMVx8y9RWHyytasWXPVVVd1\ndnbefvvtTzzxxATTkwHlQ9O0mXzngWEYyf92woEPFVMRheFOAOAYuY7X0ut/e6JXK277n45o\nzXTW865MHuz+5m/+5sUXX9ywYYPX633ggQfOPffcaSgLmMmy2ezbb78tIqeeeqrTORO/t6mq\nOn/+/GQyWVdXN8FijOEMAMdlKhNcq3acI+e2bdu+973vORyOz3/+82eccYaI3HvvvW+88UZl\nZeWHP/zhiy++eMoqPdbkf5OSyeTg4GB9fb1pmh0dHR0dHbfeeus0VAbMWIVCIZlMmqZZKBRm\nZrATkUgkEolEJl7GMIwZ2+kIAFYxVYfhcp3w5eN9I/7+97//ve99z+PxfP7zn//FL34hIkND\nQ+FwWNf11tbWqSr0eCb/m/ThD3+4oaGhtrZ2GqoBZoVQKLR8+XKHw+Hz+ayu5f1QFEV0Tbe6\nDACYWdRoTfCKT2V3b80f3HN0e+CMixyRSjleZ14ymQyHw3LUnQm33HJLVVVVNpv9i7/4i/cw\ngcR7NnmwmzNnzkMPPTQNpQCzSDQatbqEklBUlR47APgThURs7KX/EhFx/slIn2NvviQidTd+\n9Z1viUaj8Xjc6/WOzyu9Z8+ehoYGh8NRKBSmvOKjTB7sKioqMpnMLO+ZAHAsVVUNQw+GAlYX\nAgAzi6k6DOeJT8Ue71zsbbfd9uUvf9npdK5du/b555/P5XKqqt50002qqq5du3bqSn2nyYPd\nnj17Wltbly5dOn4tznPPPTfFVQGYcqqqmKbimuA6EgAoS6bq0J0nnpXneNfYLV269Gc/+9kx\njWvWrClpXSdl8mD313/91xz6AfspXggyMnLCsZoAoDyZqjpBj5153GQ3Y0w+88Q555yzY8eO\nZ5555rzzzvP5fKtWrZqGsgBMteJYJ5PMTgEA5afYY3ei/2RmDxQ1eY/d2rVrXS7Xxo0b77rr\nrvXr1997773v7GwsK0Z8JPXUL41M2kyNOarrnS0LC/veDn3iFsXjE5HEr+4PXPmpwsE9+e1v\niMttZlKByz7hqGkwYsOp5x81MynRNFf7Kb5zLtOH+9PP/tpRP9e/+hqrtwnlqBjp3lWwS8Qy\nDofDHzrxGQoAmP0MVZ3wVOwsD3YHDhxYv3796tWrReSOO+4oPjiuXC53zz33iIiu67fddlvx\nvt8J2mcpNVIZum5tbsurZj7rXbVaRPLbNxVTnYgYiZgaqsi+uj7yuTtEUY34iJFOmoVc4tGf\nBK+5wVFZIyKZPzxl5rN6f7dzbpujaqLxY2eRbCazf9/e4r0/4XBkXlubiGzv2BoMhlrmzROR\nHdu3L1i4UNe0YxYzDOPtjg4RMcXMpNNnnH2OpdtRRgxNEVEczpM9QiXj2d/95y4RufjPlnn9\nXJ4xo23euMMwTEURXTdr6qIt8+ds3rhj6akL3B6XiGx6dfupKxdpur5350GtoBuGUVVT0drW\nKCJvbdxhGqaiKLqmLzut3ef3vP9i3vj9nuWrWt0e59tvHgqEvK0L6zb9fk846m8/pVFE3np5\n3/wl9bu3dhfyWi5bCIZ9VbWhhubKnZu7CnldN4zahoq2pQ2psezbbx0KR/2LTm16/yUBEzOV\nia6xm+GnYks5V+yzzz67fPnyK6+88pVXXnnssceuv/76idtnNa2vy7N0pYgY8RG1orLYaGZS\nii8gIu7FK5Lrfu5etMLVvtwZqcxtfc2z+LRiqhMR33lXiIh76crkup95ln3wXa5Y04f6HFV1\nJf7GYBiF+IiroupEHzs6OOBwOMKVVf/dYObjMVcwVJwvT9f1HW9vX7x0WfHu6a6DBwuFQvHS\nzLGxeCqV8ng8hUJeTyd37Nm7YNHiQqFQWVk5vtgpK1aIyP59e1vnzy9+eiqV8vl8qnrCSwUy\nmUw+nx8fgFfX9Vwu53K5dF33er0iYhpGPpfN5XLhaOX4u0zTzKSSeiLuqap1u93pdNrj8Tgc\njnw+Xxxq2OPxFG9lKmTSBU13+3yGrhuGIYoSDAYVRRmLjTqdLlVMbzAkIsNDQ+l0ujoSMp3u\nfKGgqqrT6czlcpqmuVyucDisqmosFitW5ff7TdOMjY663W7TNA3THBsa8oZChUIh6Pen06mK\nqupYbNSl5RV/0O8PiMhAf38ukwmGQunREX9lpdvjHRoaKs5mFvD7o9GKxNiYIUpFRcXg4KAq\nIoaeLWgul6u6piY2OqqKxJOJfCYTCYWiNbUDvT2aKG6nU1TVNE2HKqZiOlRnMpnM5XJDg6Mi\nZjDk93g8uWx+dDitOPO1tTWVlZWxWGxsLJmKa+LRzJwMDQ0G/EGXy+H2urq7e51Oh6o6HU4l\nEPBnMplMsuD2uDw+p8/nczgcXV2Ho9Go1+vN53N+v1/TNE3TFEVNjWVU1ZlM5Ctr/D6/W1EU\nt9stItls1jTNVCrtdDqc/y0ej3s8nvF/ZdM0g8GgQ3VkM/lAyFds9/v96XTa6XRqmlb8d8xk\nMpqm19RUJ8ZSLqfLF/DEYjFDN0PhoGEYpqGqqpLNpTVNi0aj2WzW6XAVCloulw+Hg06XQ0Q0\nTSsUCuPnrDOZjMfjSSfzoYjvqB/FbCZViFYFjzuNh1bQx2LZQMjj8ToNw8hkMn6/P5PJFvJa\npCJ09JK5jK5rhtNjmqbp8fxJnDJNM53O+P2+Y1Yxwa+Jrhunn7VMRDRN3/hyRyRUVSjobo+r\nuG/TyWwymdmzo3P5B9p9fq9hmJtf3xGtDEeioUJeO+Pc5SKy++3OxFjqncEuncq6PS6n89hR\ncjKpnFYwPMfLgblsweFQ33p5X2NrdV3jkQnuRoeSpmEqqlLIaZFoYNWFCw8fGNIKeuvCOl0z\nXntx54oz5wdCXhHZu71HK+hjsXS0JhgMMT7DFMokdNM0/eEZOuj6dDJUhzbBXbEzOteVdK7Y\n/fv3X3755SLS1tb2zDPPTNq+b9++b3zjG+NPPR5PMpmMxWLvaUOmm9HdWVh5ocRiZuduidQU\nyzYP7TGjdbFYTE45W3KZwv63jZ/8s+OjnzN6u5S6puw7Ns0Y6tfEKe9mk/XnfxPf/7Zy0UeV\nRaeVcHNy2zbmd73lWbzSfbygOTLQ19112JVPz128LBCtEpFc94Hkltfc9U2hleeJyPDQUDgc\nyeVyuVxOREKRSCqV0jXNMEyHw7npjY1+f0DV8rufe8IVCHdkc3lNi0ajLS0t498TYiMjpim6\nbsRiscHBwd7e3qqqqsbGxuNXm8vt2LHDNM3Gxsbi0Nn79+9PJBJutzufzy9YsEAd6okf2JXK\nF/JuX2Dxipr6huIbe3t7B/r6wrEel8frW3Rad3d3MBhsamrasWOHYRiqqra2th48eNDjcroH\ne3RVSfmjqoiYmqiOOc2tDjH69uw0VYdH12ra2nOqq6urS0Q6tbwvl0wEKkVEUZTiyU1FUSKR\nSDAY7O7uNk1TVdXGxsbRkZFUKqWYhimKImKIuLS85vaKrvlyqT1un+ZwqWIGk7F0RY1hGIZh\niGk6Dd1UFM/BfdlgVNf14ioURZF9utvQdNWhq05FUQzDcBm6ahoFl6er65BhiqooDr2gOVxj\niWTvoYOay6Mauq46VDENUUQVRUTXs2+88UaxZj3ndLhM1WkahqEX1HzM29/X7/G68/m8iGLq\nisPt8lTnd+/daeqqlvA6faahmQ5fQVTdyDkdHkMUU0SyQz5PZU51mCKKaRoHDx4U3ak4dLfH\npWmaaZiGIVrGKYZqFBxaXvFV5p1efdGiRbqu79u3zzBMQxdRTEUxi4Evl8sdvWOL+1MRZz5j\ntrbOTWRGEolEbW3twMBA8VWPx5PNZg3DME3T6/blCzkt46yprxgdGzJN06m6dbOQj3k9Hq/m\niqkOMxAIZDKZQtKl51RnsOBWgstOb1UU2bdvXzqdnj9/vt/v7+zsHBoadXnzWgAAIABJREFU\ncqmBeI80L6yc0xoVkXw+v33722JKTbS5aV7lMT+opmlue607Ppz1+VzLz23q6+sZGYnV1lb3\nD/SbhjQ1NtfWH3mLVjA2PTucGM7XLc0YzlR7e/v4UFgi0tPdOzAwVF9fW9/wxw7+kZGRw4cP\nh8Phdw5qr2l6MpHa+PLWQl5ze1wOveKJnxyqXZJ87Q9bstlsPl/Qco7fP7G7cYE3l8/m8tk9\nWweH+pN6oWvekrp0KrPp1W25XMEf8NTOqTjmgJxMZHZtPRgM+xae0nx0ykwn81tfOeQLuBZ9\noEHX/2TIa62gawXtpWe3ty2r9QQkFotpBd0UI1Ll3bfrcGVd0DD14loG+0ar60OxWGywZyxS\n5fv/2bvzKDmq+/77t6p633u6Z5U0GmnEaEULwpLZEZvBQCDYhnjLZrzFODk8OTk+ySH87Pzi\nHOcP2xA/dhI/MTgx2ATvWLaxsVm8sQkEAoSEVkbS7NP7Xst9/mh5GLSMhNUz1VP9fh2OPX27\nuupbrZ6aT99bdUs3K5lMRQiRXBAoFPP+iHr4YC4Uc8+XvxEzsCyrCfeiUrQe/mpGKOKqv4wF\nIqc+/3621T9I2Wx2DrZ1/FbkjEOx877H7g+7V6z3hF/c3txeq9WOHDky9bCvr8+yrGMOCs1K\nCr1qqpowTWFJodeEaQohxUtPiY2Xmi/8Viw7W/iDYmC92P2Cmc+KaEIeOSAWnSWEEAdeFXpN\nDKwTpiFU7S3tr7QsRdelELJWEQ19oyy9JoRi1qonrMeo1YSQQlGN6tEFTF0XQliGYZqGEEql\nUg6FwqZpFouFidExXa/19S+rVMpen9fl9ohxaRhGUFN1tydg6kJaQshCNrt75yt9/ctcLpdp\nmul0qrdvSX3l9U5iXddP9uaYpln/G28YhmmaUx+b+p9/wzC0TMrIZ1XNLTx+yzCm1mOaphBS\nClWYZn042LKs+uamViiltCxLCqlIRVGkIoUUipTSMAxpmVIoQgpLCtMw9Om/2pZV///p+cM0\nzalS61ufqkQRQliW6nIJpX5fL3XqrA0phVTe+EWQQtRfLxVRf7ellJqmWZalqKq0lPrTlmUJ\nIUwpFSEUIaQQqqpKKZV6Mb8vs76qeomKokjLUhVFKIqUUkopFGGallCllELR6suov69fHj0n\nT1iqqkqpCCGEtKSiyPpBTlGkEKpQpJCKS0hLKC5FSCGFUISwLKFoov6eSykVoWqaYpqWIlRV\nqtISx/+LK79/P+uN088FrPdRWdJQXUpNN+r7Pn1soR7phBCKUCzT+v16rOkrUVSlXNQ9MUVK\nq74Gofx+D4U0DENRjq7HMIyjIVsIy7SEohi6Of2zWl/m+I+rlMI0pBCiXDIM3TDNaXUq0pz2\nEr1m1d//SrnmjSjHrO3otqZ9fqbWY765sa6QK3Z0xTt62saG05qmZo643b6i0EPLVnYfOnQo\nNVEzyoqiSc119PhjGqaiSUWohVyxvTPWtTAhhNiz81ClXHV73vQ3Qq/pQhGmaeq6oWnqtHZD\nSFEu6qZ57DE8lymrqhKO+9ITRX/IXW/xh9xtXcHXXhhxeVR/yFN/SSFbWdAfM02zVKiFIt7j\n96tSqrm96jz5G3EKTbgXes2qHztqFcMbtH/e8qPfNufkjbJ+fwB/o0XVDAefYyeEuPTSS2c4\ntW7KsmXL9uzZMzAwsHv37v7+/lO2r1y58tFHH516eMcdd0QikUQiIZqemR4vtXeH66XGNxce\n2mX9/FvCND2rNvpWnG2EQ6XHf6B4vMI0Xf0r/X39YvHS4s8eNH/2TSEULdERuPyPK7/7uX5o\nn5lLuR7/fvhdt57mdicnJ9UtN4T0smvh0voYaKNYF1yhD6xxd3Sr7hN8lOPxeCgQcGlqx+Kj\nQ6UyHqt0LfDE2jSfXwih12qlUjGRSCQSiUULF+144YXOzs7Dg4MdnZ3RaKxaKQsp2tvby9FY\n1ZLrehen0+nOjo6dL79cv0X96wcP9C1ZGm9rm9pcV1dXOBye4TasPp+vWq12dnbW/8aHw+FK\npVIfio1EIim3FmrvlMFwTXO39yycGquKxWLZjnYz3RHsWewLBru6ugKBgMfjiUQi5XLZ5/OF\nw+GOjg6XplWyqYpuBCIxyzJrtZqqaolkUlXVcCTidrncLi0UTwhF8bjdlXKppy1uuT0V3VAU\nxev1lsvlUqkUDAY7Ozs1TQuHw/VR41gsVqtVR0dGfP6AIkSlXC5kUm5/lyllxO8rFUuxjs7U\nxLi7VtUiSwbaElLKwddfr5XLwUCgVsz4ero6vIGxsbF6f5XP50vEY4Vc1lJd4XB4fHxclVK1\nzIphutyujq7uVColjVqpUhOG4XVr0Xh3JpWypNDc7nqSK+TzqhBC01auXFksFsdHJqVXDcW8\ngUCgUq6lxgrBBVZ718Lu7u7x8fGJ8ZRRVfwdrkJRxOPxoDfm9Xu8PvehwUNSmJrq8/ld/oC/\nUMiX89LT7vaH1FgspqrqgQMHgsFwLBrTjWooFDJNs1gsKopWzFYNXVQKRseCaDjmUVQlHA5L\nKdva2nRdL+SLmktzubRAIOByuUZHR8PhcP0WvS6XS1XVSCSSmkxpqqe3b4Gu6+VyORKJ5HI5\nVVUty6r32KXTGcMwFy1akEnlPG5fLBEaGR4xLautLW5ZljQ0zaVmcpOmaXZ3dxeLRVVx16q1\naqUWT8SCIV/9c1WtVr1er2VZK1asyOfzfl+gkKvGk6GpWOP3BYu5Ss+idlU7QSfH5i3h8ZFc\nPBEMx/ztnW2FQjEajWQy2Wq51tXTMb3TK/TOSK1q+kLSklY4/KZR2lgslsvlI5Hw9Bv7xuPx\n9vb2YDBYH8KeLp+ptHcmE4l4MBje+eLecy5Z4ovqic5QIpGIRCK7dx5QE17VJTSXqB9v5Rpt\n767aOecPDB4cbu9MJhJtQoh96lBbW1v9nLwpbW1t8ba4x+MOBH3T2xMJEQ6HSuViJBo8pvjs\n+Ej/yp7u3sQzj+/yaP5wLJAdH+laEE20hyc7yvnJWmdPIpGICyFUZaSjo10IUS0quXSxXtvI\n4bRpWAv6EpYlvd6JefEH4pRSqVRb27H9u7ZLJMQffSwqpUx0N+DEyjOXzWbrJ+2c8CSHxjr+\nfAZL0QztrU138vLLL3/xi1/UNO3WW2/dtGnTCVvmhjLzNXE//OEPb7jhBiHEjTfeWL8P2n//\n938f83s7pVar3X333fV8/clPfjKdTt9zzz133nnnMe3B4Ilnur/jjjtuuummc84550z3ybkm\nJyc1TYvFYnYXciwp5b49e0qlYv3Xo2fBwrZEYtfOV5YNLHe5XNVK5eknf7fpvPO9Xu/xiwkh\nXnz++XUbNjTwO1Amk4lGo3NwOJjXHnvsMSGEkHLLZZdNb5dSTk5OulyuJvyk2aI+qjt9bLSB\nyuXy4cOHY7FYe3t7o9b50vOv9S/vrWevZ367Y8PbVu56+UD/8kWBoF8I8eJzuwdW9vn83l0v\n7y+XKoqieLzugVV9brfrpedfq1ZqLrdmGFZHV1vvku7T36hpmul02uv1HvMH4oUn9y1bvSAU\n8ZUK1Ree3Pf2y1fueHr/6o19bo9WLtUe+9GLW65b6w96S4Xq7hcPbbhgmRBCSrHz+deL+YoQ\nIhj2rVy/aN+rw6nxQjFXjiVD51ywrFFvlF2aM9g1m3qwSyQSc3AkT6fTV1555cUXX/yFL3yh\n3nJk8PUjhwYNXZ8+FCCE8Hq9iqqefc5Gv//YA8KHP/zhz3/+816v99Zbb/3GN75xwpa5MVOv\nz5e//OX/+I//uOqqq/x+/3PPPff973//O9/5zuc+97nPfvazJ1ze4/H83d/93dTDYDB45513\nHt8O51EUZdnAwDGNK1atrv/g9fku3nI0Ohy/mBBiHWneDqqhWy6326wd057L5Xbu3CmEOPfc\nc2cpzWBK/YzSQqHQwGB39jlv/JZtumDtMS3rNi6v/7Dy7KUzvLBR1p93dJQmEPKef+Wq6S3+\ngOedt7xt6tkNv09siiJWb1w8fSX162eBOaN5vInOnnQ6U8nlprcn2zq8Xo+inKB7vt75Jaad\nFnJ8y9yYKdh95Stfefjhh+vXOWqadu65565YseK88847WbADWoQ0jerQoDue1ELzeO6eoF5W\nyxndFzqm3ev1BoNBTdO45cwcSCQSxWIxHo/bXQiAN1RK+oG99WsA3nS64dChcSFET0/P8S+J\nx+PZbLY+AcLJWubGTMHOMIxFixbVf/7Upz4lhAiFQrXasd/vgVZTeu2lzO9+4Wnvbr/+/U1+\nFu0M8r6Qy/Iff+WXz+dbv369qqqMZc+BYDC4cuVKu6sA8CbSUiz9pJcGH7167M1uu+22v/7r\nv3a5XJ/4xCceffTRarU6vWU2iz3WTMGuVCpNzUb28Y9/XAhRLBbL5fIclQY0K9Xrl1IqXt+p\nF21imjQVaWni2MvBhBDTT9IHgFYzc7A7oVWrVh1/Xy5b7tQ1U7C7/vrr77zzzs9+9rP1U91r\ntdonP/nJ9773vXNVG9Ck/EtXdCW7XKHw/O2uE/XJPRShzeddOBP1iUKOv6QUAKSpmLWTf79t\n7jtszxTsPve5z91www0rVqw499xzDcN46qmnNmzY8O///u9zVhzQtFyReX/FqKG6FCFP1GHn\nfJZl7dixI5/Pb9y48WTX6QPzyyP3jZimVFWhV+XF7+oIx11CiMmh6gtPZN7+zkQw6hJCPHL/\nSLLbu+GyuBDisQdHz78u+avvjV94Y7svoD31k4lIwr36vGh2Qn/qJxPJHu/GK1r30t0/YCi2\necwU7CKRyGOPPfbkk0/u3LlT07S///u/37Bhw5xVBmB2KWoDv3bWajW32z1fTsubmji6Pmdy\nYxUKhYmJiWQyGQode2EKMHuqJfO6jywQQjz908nUSDUcdwkphg+Uz9oQHj9SrQc7IcTooYpl\nSlVTKiXLG9DKBdPlVh7939GzNoR6VwSFEKmRatdiXzTZ0p3Z0nzLQ7HN49ST3J533nnnnXfe\nHJSCEzLy2cxTj1q1ijQMmezRVm/UM5OZJ3/hSXZF33aJEGL84W8LaQkhPB090Y0XGfls9tnH\nrWpVmoZv0dLIurcfXcYyhaJIw/AvXhZeu3n84W9L01RURQjhX7LC09aeff63Uq+FVp8TWLry\nmE3Akepz+dYvez9DmUxmx44diURi9erVZ762OeB2uwcGBnRdn7rdcAMdOXJkbGysXC5zVQTm\nTLVsFbLGow+MFnNGW5dn8cqgEGLfjsLCgYDmUnZvy/etClbLltevdvb6DrxcXLwq6HYr1bJl\nWfLHXxu68Ib25IKj8xIvWRN64rtjS9e29NcSaamm7sShWNhOGvrkL3+Q2PJHrmhcCDHym0ek\noevZCW/XIlesTQghpJSm0XHte6ct/8PE5X/kCseEEBOPfK82Puxp755axqpVxx66L7x2s5nP\ndr3n6B0vpGmOP/xg+zveI1RFT40LIfTJsTc2AYdSLNNlmrVq5cxXVb9YXtf1+i13z3yFc2D2\nutPi8Xi5XGYGWsylyaHqirdF1l8aF1I8/N/DpZzh9qo7n8rWO+oM3aovk+j29q8L//JbI+G4\nK97lmRyquj1qV59v+GB5KtgJIQppIxxv6dmOLEsxHTkUC9uV9u8KLF1RT3VCCPfKcxRNC/Sv\nnHzsR4Flq4QQtclRs1iY+MX3Va+/7YKrSgd2+/vOqqc6IYSvZ7GeGndFYkZmcvzhb5ulgisY\niV90tVUtK9Nu2lsdOqi6PenfPSIUJX7BVUKI6ZuAU3mqJc00LLUBZ5glk8lVq1YFAoH5kupm\nVUdHR0dHh91VoLVMDlfjnUcHT42apajKC4+nL35XRzTpFkI8/PXh+jLJBV6XW2nr9Lz2fH7R\nQGByuLpiU2Tp2aGH7x1asCzQ1ukRQlimVF1NHVzmgDTFDMHOOT12e/bsOeuss2avFBzPyKU9\n7Se4q4+Rz9TTmysUab/mZlc4mnnql+XBvUZm0tO5cGqx2uRoaOWG2sRoeO2m8Nmb8i89q3q8\n3s4FlSMHzVxm7MffEkJoPr8rlvR2LQqv3ZR5+rHKoX3+xWdN3wQcy+U2VMXnbcBQrKqqyWTy\nzNcD4A8zOVQd3FV67bm8XrOWnB2qVax82qinOiGE16/m0/rkcG35uREhxMrN0Qe/MLjukviB\nlwuLlgcVRVx4Y/tj3x679kM9O36VGXm9nJ3QH31g9LI/6bR1n+wkLcU0/pCvqcffH/bf/u3f\ntm3b1tbWdv31119++eWNrvQE3kKwu+aaa/bu3Tt7peB47niyOnK4nrTKB3YbmbS2ZIW0TEU7\n+g9XPrQ/0L9SCCFNS/X53W0d1ZFD/t5+IUR19IhZyHnau3MvPu2OJ4UQgaXLU0/8JLh8bW1i\nNHruxcEV6+orye94pr6AVSlpwYgQYvom4FS6UFWhFphyHJj/Lr352BA2PZZtuaVTCHHpe452\nJIdirr/8p6XTXxVuc//RRxcIITZcFheC+6AIy1LM2h8yFHv33Xd/8YtfnH5/2ImJiUgkYppm\nX1/fbJR6PP54N7XA0pXV4UNjP/6WEIo71uYaWF995dnx9LiRmZx45HvJK29Svf7xn/yv6vG6\nwlFv1yIhZXV4cOzHDyiqqgWCictvEELoEyOBJQNCCC0YsWpVq1LSJ0b8G86f2kpwYM3kY1sL\nO7e725KeZGf2uV9Xhw9NbcK2nccsU4QUQmEmYgA4hj/gXrU5OXKgNH7kTWch96+LBMIzDVQf\nf3/Yj3zkI4lEolKp3H777V//+tdnr+YpbyHYzfE9MSCEEIoSv/AdU48mJye9q98Wi70xQurv\n7a/3z51w+bp6vKvrvPHPjmkRQqi+QPs1N089jG68qDHFo7lZbq9lGOFweHh4uFwu9/b2ulx8\n0wMAUS5ae18sVkum8eYz7Q6/VlFdysA5J7ia/mtf+9r27dsVRTnm/rB79uzp7u7WNG02Jlc6\nobdwHL/99ttnrw4As0pKOb7nVUVRk/0DiqoKIUzLEqpaKldee+01VVVDoRCn/AOAEEKvicyo\nKYQQ4k3BLjdpCXHiiyc+9KEPCSF27tx5zB1jVVX9i7/4C1VV56x3jC/oQEsop1Oju3ZKRQQS\nyWDbGxc6WJbV29tbqVRmY0Y32x08eLBYLPb39/t88/vGvgDmkjRV4+Tz2M1wjt0J7xj7jncc\nO5I2qwh2QEvwRWPRBYtUTfWFjwY4RVGklFLKJUuW2FvbLDFNc3BwUFGURCLR1dVldzkA5g3L\nFMbJL55ocjPVncvlvvzlL9d/vu+++84///w/+7M/m5iYmJPCADSSqmm9GzcvXP82zX10BgRp\nWUIIB09XpWnawMDAggULEomE3bUAmE8sSzF09WT/yfk7j91HPvKR+uTp+/fvv/322++9997t\n27f/7d/+7fHdjADmISmEUr8fnVPRUQfgD2CZykw9dvP3zhOvvPLKjh07hBDf+9733v/+9193\n3XXXXXfdOeecM1e1AZhFiqJKIRzdZwcAf4h6j93Jnm3yHruZhmKDwaCiKEKIxx9//LLLLqs3\nyibfIQCnp35tbDActrsQAGgulqXounqy/5o8Bs3UY2dZ1sjIiGEYv/vd7771rW8JISYnJwl2\ngDNEIpF8Ps8UJwBwDNNUarWTjmY0eQ6aKdj9wz/8w5o1ayqVyv/5P/8nHA6n0+ktW7Z87GMf\nm7PiAMyebDYrhDhy5Ehvb6/dtQBAE6n32J306fl7jt2NN9549dVXF4vF+jVlkUjkc5/73KpV\nq+aqNgCzyOVy6bru9/vtLgQAmotpKrWZ7hV70heOj49/6lOfWrp06R133FFvefnll7/4xS9q\nmnbrrbdu2rSp4aUe7xTz2Pl8vvrEnqVS6Xvf+9699977yiuvjIyMzEFlwOmQpqloDbvfqaXr\no089LjSt6+2X1k9BczDLshRFMU3T7kIAoLmYppgp2J38hZZl/c3f/M2PfvSjqZa77777i1/8\notfrvfXWW7/xjW80tMwTO/UExU8++eS999774IMPGobx1a9+9d3vfvcclAWcjsrhA7lnn/Au\n6ItuurQhK6xmJgvDg0LK2qr13lhbQ9bZtEzTVBSlWq3aXQgANJe2Tu1Pbvdve8x45Zk3ffW9\n9k89ye76ZaUn1tnZmU6np7cUCoVIJCKEMAxjdoo91kx9Ep/73OdWrFhx8803t7e3P/fcc8uW\nLXvf+97n8XjmpjLglPTJMT01XhsbatS5rP5kZ3LtpvaNFzg+1QkhQqGQoig9PT12FwIAzWX4\nkPzP/6s/9xtZqanT//vufxn/+X910zxBsvva17522223ZTKZY9rj8Xg2m61Wq4FAYE5qn7HH\n7o477li9evVXvvKVLVu2zBRQgTNj5DPSMNzx5KkXfbPg8rWKy+3p7BGN+nwqSnzF2Y1ZVdOr\nVquWZaXT6b6+PrtrAYAmYlhK5S2eY/ehD31ICPHEE0/cc889e/bsaW9vP+uss6rV6m233fbX\nf/3XLpfrE5/4xOwVPN1Mwe7QoUNf//rXP/rRj5qm+ed//ue1Wm1uakJLMcul8V/+SE+Nd77z\nZm/Xwrf0WtXnD61mxuw/kK7rUspisWh3IadmGMbu3btVVR0YGNAad0olAJyQbijFykn7C6yT\nDxFdcskll1xyyTGNc3y/rpmGYru7u//+7/9+z54999577969ewcHB6+//vqHH354zopDK1AU\nRVGEFEI4/WKFZqMIoSiKOh864wuFQiqVGhsbK5fLdtcCwPlqhpIvqyf7bx7PYzelnkC/9KUv\nfetb3/rHf/zHq6++erbLQutQff7k5TdIw3C3wGltTSUajeby+YUL31ovqS2i0eiSJUtUVQ2F\nQnbXAsD5aoaSr5y0r8Gav/PYCSEsyxodHe3s7FRVNRqNfuxjH7viiivmpjK0DlcoYncJrahY\nLErTnEylehcvtruWU1AUZV4EUADO0Nkt1m48aXrzeueylrdspmD36KOP3nLLLYZhJJPJrVu3\nDgwMfOlLX/rsZz87Ojo6Z/UBmCWWZQkhDF23uxAAaC7XvUe57j1N3S03g5mC3d/+7d9u3bp1\n8+bNDz744Ec/+lEpZTgc3rZt25wVB2D2dPf0pCYm+pcts7sQAEDDzHS6erlc3rx5sxDi5ptv\n3rVr11/91V9t3bp10aJFc1UbcKZqh/dX9r/a7Hdstkk+ny9Xq8fMpQkAmNdO9zrErq6uW265\nZVZLARrLzGcnf/SN9M+/XTm01+5amlGpVJJS5vN5uwsBADTMaV0VC8xHqj/g7u4VhuEKx+yu\npRl5pCUtw8+0cLPp8OHDpVJp0aJFfr/f7loAtISZgt3evXuTyaM3A8hkMlM/T0xMzHpdwBlT\nXO7E9X8qpKW43HbX0oxKhiFV10Q6tdzuSpzKNM19+/ZJKYUQAwMDdpcDoCXMFOwOHz48Z3UA\ns0HRNCHokToxtxBCrwYjb/lObjhNmqYFg8FCoWCa5qmXBoBGmCnYMRco4GA+SzeNqsc07C7E\nyZYtWzY+Pt7Z2Wl3IQBaxUzBbs2aNdMfKopiGMbhw4cl1xgCb5b72YP+dee5uxZVdm03JkZD\nF16d+cG9wjSFokjD8PavDGy8+PgWIYQxdqT49KOhS6/XwjEhhJEaLzz+kKtjQejCWb+/S1Uq\niqoVq7X9+/d7PB5mAJ4NsVgsFuMUTwBzZ6Zgd/DgwamfpZT/8z//85nPfOaf/umfZr0oYL4J\nvv2Kwq9+HNh4ce31PZF3vEcIIQ09/u6PCCFkrZJ64N8DGy8+vkUIWTu0z7dqozFy+GiwGx9y\nL1jiamufg5prLo/U3C6hHDlyRErZ3t7ubfL51AEAp3JaV8Vu3779E5/4RHd39+OPP97b2zvb\nNQHzjhZtU9ye4tOPRv/oT4VQrErJTE9kfnCvVcip4VjkipuObxFCVF59wdO3XHG5yzue9p61\nRgjhW74u99MHfCvWz0HNXq+3UqnEolG3x+P1ekl1AOAApwh26XT6jjvueOKJJ+666y7uEguc\njJnPWMW8lFLRNCGEMXYksPGiwDkXlZ7/jeLxunsW1wb3HNMia9XS9t9qkZiQUhpv3NfLzKW1\naNsc1BwJhYxaLZlMdvf0zMHmAABzYKYJiv/rv/5r48aNfX1927dvJ9UBJyNr1fwj342842ZX\nPFk7tE8IYYwNafF2IYT3rLOru188YUvxqV9Er7klet0Hotd/UP7+qklpmmKuJpYbn5gwLGvP\na6/NzeYAAHNgph67D3/4w0uXLn3ggQfuv//+6e0vvPDCLFcFzB/Syv3swdCF16ihSGDTZfmf\nf9uzqF8fPeJdtkYIoYWjVrVilYvHtOhDB81cuh71hBCqL2BmU5Wdz9UOHzDTY9kffSN6/Qfn\noHZFURR1vt7oGgBwvJmC3fDw8JzVAcxXijoVwrRwNPauW4UQ0WvfN/V82/tuO3FLT99US/S6\n9wshguddGZyLio/SXC7DMELh6BxuEwAwu2YKdl1dXXNWB4A5ZpqGEKLIvWIBwEFmOscOgINJ\nSwohFMW2odhCoTA0NKTr+qkXBQCcntOa7gSA82hCSClV+4Ld/v37s9msYRhMogQAjUKPHdCi\nuhYs8Pl8y5Yvt6uAcDgcCAS4dSEANBA9dkAjGZMjped+7Vm0zLdyg921nEK+UKgaRi6Xa2+f\nixtdHG/JkiV9fX02jgUDgPPQYwc0UnXPK9XdL5ZffNLuQk6tVCqZppnNZm2sgVQHAI1Fjx3Q\nSN5lq83MpGdRv92FnJpHEaa0Ql633YUAABqGYAc0kivZFbn6ZrurOC2maahCVipVuwsBADQM\nQ7FAiwq4NXetHAr47C4EANAwBDugRRmGqQhRqzGNHAA4B0OxQItyR2IlU3ojMbsLAQA0DMEO\naFGJZLslRSKRsLsQAEDDMBQLtKixsbF8Pj8+Pm53IQCAhiHYAS3K4/FomubzcfEEADgHwQ5o\nUbVazTTNSqVidyEAgIYh2AEtKiiNUDkX0Lj3AwA4B8EOaFFWNqUC5xFeAAAgAElEQVSWi0Y2\nbXchAICGIdgBLSpnKVW3L12u2V0IAKBhCHZAE6lNjOR2PGsWcnOwLd0bKAUihot7xQKAcxDs\ngCaS3f5k9vnf5F7eNgfbCoXDiqIk2tvnYFsAgLnBBMVAE/F1L7LKJW/nwjnYVl9fXyaT6e7u\nnoNtAQDmBsEOaCLhNeeGVm9UlLm4UjWbzabT6Vgs5vF45mBzAIA5wFAs0FzmJtUJISYmJvL5\n/OTk5NxsDgAwBwh2QIvq7OyMxWIdHR12FwIAaBiCHdCi/H5/KBRiHBYAnIRz7IAWdejQoVwu\npyjKkiVL7K4FANAY9NgBLcrtdksp6bEDACehxw5oUX19fYlEIplM2l0IAKBh6LEDWpfLxVc7\nAHAUgh3QokzTLJfLdlcBAGgkvq8DLergwYPpdFpK2dvba3ctAIDGoMcOaFHValVKWalU7C4E\nANAwBDvAyUrDg5MvPGWUCsc/RbADAOdhKBZwMJl+aVt1ckx1u+OrNx7znGmaQgiCHQA4CcEO\ncDAluGip6nYHuk9wFp1LCF0Ij0a3PQA4B8EOcLLYyvWxletP+FR7LVcqlTqj/jkuCQAwewh2\nQIsKxhOiVg0kO+wuBADQMIzCAC0qk06Va7XUxLjdhQAAGoYeO6BF1VSX1Fy6otldCACgYQh2\nQIuyYu151RuLttldCACgYRiKBVqU2+1WVJXbxQKAk3BMB1rU4sWLY7FYd3e33YU0TLVaPXz4\ncDgc7ujgihAALYoeO6BFqaqqKIqiKHYX0jCpVGpoaOjQoUP1uZcBoAXRYwe0qAMHDqTTaSll\nb+8Jpi+ej2KxWFtbWyQS0TSuCAHQogh2ABzC7/evXr3a7ioAwE4EO6BF9fb2hkIhJ51jBwDg\nHDugRR05cuTQoUNjY2N2FwIAaBiCHdCidF1XFEXXdbsLAQA0DEOxQIvq7e3NZrMMxQKAk9Bj\nB7So4uRY8chgpZCzuxAAQMMQ7IAWNbZndy01NrjzJbsLAQA0TBMNxUopdV2vVqt2F9LUpJS8\nRadUf5ecNPXubDBcHkVK1eXhEzUzwzAsy+JdmpllWUII0zR5o06Jw/jpqH+i5uZIXqvVZnsT\nc6m5gp1pmoZh2F1IU5NS8hadUv2zZHcVzU5G4uVKJRiJ84mamWma/N6dkpRScIA6bbxLp1T/\nRM3Nkdxh/xxNFOxUVfX5fMFg0O5CmlelUlFVlbfolHRdDwQC9NjNrKenJ5VKLViwgE/UzCqV\nimVZgUDA7kKammmalUrF5XLxcTqlarXKu3RK9Z7yuTmSO6zHjnPsgBalaZrL5eLuWwDgJAQ7\noEUN791dHtw3dHC/3YUAABqGYAe0KLWQdtfK5fFhuwsBADQMwQ5oUdIXNtweXzRhdyEAgIZp\noosnAMyls87dnE2nunoW2F0IAKBh6LEDWpTmcvlDYUXlIAAAzsExHQAAwCEIdgAAAA5BsAMA\nAHAIgh0AAIBDEOwAvEmpVOJOuwAwTxHsgBaVGh8bHjxYyOWmN05OTm7btm3nzp31O3ADAOYX\ngh3QolLjY/lMOjM5Mb3Rsqz6/xLsAGA+YoJioEUFI1HdMEKx2PTGZDK5du1av9+vMr8dAMxD\nHLuBFjWeSheq+sjo2PRGRVFisZjX67WrKgDAmSDYAS1K13UpZalUsrsQAEDDEOyAFhWLxdxu\nd09Pj92FAAAahmAHtChfaiQ6ckDNpewuBADQMAQ7oEVZmQmXUa2MDdtdCACgYQh2QIuKrVjr\n6ujpWL3O7kIAAA3DdCdAiwovWuJKdvnePN0JAGBeo8cOAADAIQh2QIsaHxs99PrBXC5rdyEA\ngIYh2AEtanJ0pJRJT46O2l0IAKBhCHZAi1JNQ7UMlXvCAoCDEOyAFmVpLkuoknvCAoCDcEwH\nWpTLNDy1siotuwsBADQMwQ5oUVox69UrSjFndyEAgIZhHjugRRVUj+o2paXYXQgAoGHosQNa\nlCeeLAVjkWSH3YUAABqGHjugRfUuXjw5Odnd02N3IQCAhqHHDmhRhw8fHh0dHR4etrsQAEDD\nEOwAe1THh1O/+0Vl6HW7CvB4PFJKj8djVwEAgIZjKBawR/G1lwv7dpqVsq9nsS0FdHV1eTye\nZDJpy9YBALOBYAfYI7BkwKxVgktX2lXA4OBgLpdzuVx9fX121QAAaCyCHWAPX89iu/rqjqqU\npGmKSsnOGgAADUWwA1pUPD/hzaYDMb/dhQAAGoZgB7So0LKVysiR8JIBuwsBADQMwQ5oUb6O\nHhmOu0MRuwsBADQM050AAAA4BMEOaFGvv/76a6+9NjExYXchAICGIdgBLWpiYqJarQ4NDdld\nCACgYQh2QItKJpNer3fBggV2FwIAaBgungBa1OLFizs6OmKxmN2FAAAahh47AAAAhyDYAQAA\nOATBDgAAwCEIdkALk5bdFQAAGolgB7SozO6X0s//rjA0aHchAICGIdgBLaqaTevFfDWbsbsQ\nAEDDMN0J0KLiK9YWxkeiSwfsLgQA0DAEO6BFeSIxv9vr8vrsLgQA0DAMxQItqlIs5MbHjFrN\n7kIAAA1DsAPmG8syJ0eFdaYXtI4c2D9xaHDs8KGGFAUAaAYEO2CeKT71i9R9d+d//ZMzXI8/\nFPb4/f5gqCFVAQCaAefYAfOM1GtCUYShn+F6YrGoq1YKR8INqQoA0AwIdsA8E3j7Fd4lK1xd\ni85wPdnXXi6nJ7wul+/scxpSGADAdgQ7YJ5RvT61d9mZr8cTTxp6zRdPnPmqAABNgmAHtKjY\nWau8C/uCsbjdhQAAGoaLJ4BWpthdAACgkQh2AAAADkGwAwAAcAiCHQAAgENw8QQw6yq7tlde\nfUFxu61KKbzlBleiM/ODeyNXvUcNhIQQ6e98NXrtB8ov/FYfHpSGHlh/vndgbeYH90auerfq\nD+Yfe0iLJQPnXJj5/r3SMhVVlbVq5OpbtGjbGy2mEbniJi2WNFLjhccfcnUsCF14td07DQCw\nAcEOmHWlbb9qe/8nhaKaubRVLgohZK1aT3VCSquQU70+V9ei4HlX1g7tq+55yTuw1ioWFJcn\nu/V+36qN3v5VQgirUmp77yeEEPnHf2SMD2vRNjOfTvzp/yOEqA3uLT336/Dlf2yMD7kXLHG1\ntdu5twAA+xDsgFnnXbY6+9P/9Z21xrN0pRaJW5WSMTmW/s7/J4SQek2LJoSqursWpb71/2rh\nWOQdN1uVkpBW+sH/CF/5LnfnQiGEVSlZ+XR2631mPuNq7/EuW21VSqovUF+/kRrTInEhhG/5\nutxPH/CtWG/jzgIAbESwA2Zd8O1XWJVyde/L6Qf/M/6ejxpjRwIbzg++/QohROW1Hcb4sBBC\n9Qfb3ntbdut9ZjZllfKK2+NesEQ/tK8e7IyxI/515wU3Xy6ETH/nv6xCzkiNmrl05gf3CqFo\nsUToomvq2zJzaS3aZuPOAgBsRLADZldp+299y9epgZB/zdsqr263SgVjbMiV7Ko/a06MuDsX\n5B97KLhpixoMKy63kNIYG/KvO8+3fF36O1/19A24kt3G2JArcfQlUq8JVTHGhkLnX+VbtXH6\ntqRpCk2b6z0EADQNgh0wu9xdi7I//qbi8QjL8ixcokVi+ugRb//q+rPG5Khv9blqOJ798TcV\nr88VS7jau4vPPBZculIoSuSqd2d//M34LR/XR49U9+0sv/yMrNV8y9epgbA+eiT49sunb6j4\n5CO1wwfM9Fj2R9+IXv/BUxZWKBRSqZTf73e5OA4AgENwQAdml7u7N/6ej0xviV77vjd+vv6D\nQggtKuI3f/T4BbRoou19nzzmJcevpC543pXBt1LYgQMHKpWKEGLlypVv5XUAgObFPHZAi7Is\nS0ppWZbdhQAAGoZgB7So3t7etra2RYsW2V0IAKBhCHZAi9JU1aVpmspBAACcg2M60KJGhodT\nExMjIyN2FwIAaBgungBaVCwetywrHo/bXQgAoGEIdkCL6ujsDEcikWjU7kIAAA3DUCwAAIBD\nEOwAAAAcgmAHAADgEAQ7AAAAhyDYAQAAOATBDgAAwCEIdgAAAA5BsAMAAHAIgh0AAIBDcOcJ\noFWktt4fu+wGNRCaapHp8dwzj0QveqcWjtYXkKapKIo0DW/f8tCG81Nb75OmqaiqNM3opde7\nYgkjPZH79U/d7d3h864QQryxgF6LXXGTFuEGZQBgJ4Id0CqsUn56qhNCWqOHAwPr9LEj9WAn\nDT1x458LIWStOvHd/wptON/MZ9vf+wkhRPXw/uKLT0YvuU6fGPH09Lri7UfXWSkn3/1hIUTu\nNw/rEyMEOwCwF8EOaAlWpax4A9Nb9H071e7F7ra2yq7tvv5VVqVsZCZTW+83izktFI1eer1V\nKatef31hIzWuhWNCCP9ZazKPfNc/sLa+TjOfTT/8v2Y+5052+paunPv9AgBMR7ADWoI+PuRu\n7556KGvV6s7npD9YcHuEqdcXCK0/P7j+vOKLTyoer6e7t3pon5nPpLbeLxThiibqY69CCCOX\nqffM6eNDwTVvC73tEiHk5A//xyzmtWDYlr0DANQR7ICWoI8Pu9u7ph7mn33Cd9E7c5aa6OrK\n/eSbRxdIdAghfP2rM7/8QWDlOfr4cHjzZf4V66evR5qmomlT63QlOo626zVFUeZoZwAAJ0Gw\nA1qCPj5c2fdq+dUXhBChzVvMfHpwIpPJZKSUYa/fzKX18WFf/yohhBaKyFrFKpf08WHf2y6d\nvpL8M4/Vhl430hPpn/5v/Jpb9PHh6sHdpZ3PS73mP+vsN5/ABwCwAcEOaAnxd7xn+kPP1bcc\n2rHj6FNX33zMAsn3fOT4lwghwpu2zLBOAIDtCHZAi+rr64vFYt3d3adeFAAwTzBBMdCiNE0L\nBAKcGAcATkKwAwAAcAiCHQAAgEMQ7IAWJaWUtardVQAAGqmRF09Uq9W77rpLCGGa5m233RaJ\nRKY/u2PHjn/5l3954IEHGrhFAH+w/M+/Y44eyp9zYXzzllMvDQCYDxrZY/fII4+cffbZn/rU\np7Zs2fLDH/5w+lOlUumhhx5auHBhAzcH4EzI9JhiVGvDr9tdCACgYRrZY7d///6rr75aCNHf\n3/+zn/1s+lP33nvvBz/4wa985SvTG/ft2/fpT3966qHX6y0UCplMpoElOY9pmrxFp2SaZjab\ntbuKZudaf4F5+IB73Xl8omZmWZYQolar2V1IU5NSCiFqtRofp1OyLIt36ZRM0xRCzM2R3GF/\nLxoQ7LZu3bpt27a1a9dOb/R6vVM/P/XUU11dXYsXLz7mheVy+dVXX516uHr1atM0DcM485Kc\njbfodPAunZK6eLm6eLklhMV7dRrq8Q4zk1Lyq3c6eJdO09y8UfUQ6RgNCHbXXXfdddddJ4T4\nyU9+smfPnoGBgd27d/f3908t8OSTTwoh7rrrrqGhofvuu+8DH/hAvX3NmjXbtm2bWuyOO+6I\nRqPJZPLMS3KqyclJTdNisZjdhTS7TCYTjUaZoW1mxWKxXC7HYjGXi4nKZ1KpVCzLCgQCdhfS\n1EzTTKfTXq83HA7bXUuzS6VSbW1tdlfR7LLZrK7riURiDo7k2u/vf+0MjTygX3HFFXffffdL\nL70khPjkJz95+PDhe+65584777z99tvrC3zqU5+aSnUAAABorEYGO4/H83d/93dTD4PB4J13\n3jl9gX/9139t4OYAAAAwHfPYAQAAOATBDgAAwCEIdgAAAA5BsAMAAHAIgh0AAIBDEOwAAAAc\ngmAHAADgEAQ7AAAAhyDYAWg5lmVVKhW7qwCAxiPYAWg5u3fvfuaZZ8bGxuwuBAAajGAHoOUY\nhqEoimEYdhcCAA3WyHvFAmgGIy8+W81nhRClibGFmy9OH3jNF0t0rFonhDj89K+6N2w+/Oxv\npGVZlmWaprK4vzhyREpLCOFva68vJk1z/2M/bV+5NrKg9/iHDjAwMFAqlaLRqN2FAECDEewA\np+la9zYhRGr/a4FER7h7YfrAa6XJMWlZiqoa1Yrm8eql4rIr/6hYLJbL5Wg0mh8a7Lvoyulr\nGHlpm6IqvljbCR86gNfr9Xq9dlcBAI3HUCzgQNlDB2qFXPvKtWatqnm80YV92cMHLdNQXS6z\nVtXcnqklK9m0US4deuqJoeefqvfbTe7ZGepaKBTVEwwd/xAA0MwIdoDT5EeOFMeGu9aeK4Qo\nZ1K+aFu0d2l2cH8lk/ZF4+VMqlbMH/z1IyPP/jq7c7vbH+i94PJFb79EdbkKw0cKo0OWYYQ6\nu1VVFUIc8xAA0OQYigUcpTgxmn1934JNF9UfVjIpfzyhapovGk8f3BPuWljJpDpWr4/3nVUf\nii2MDbct7hdCSMvSPJ7x3S9rbs+h3z2ul0ulybHJva9OfxhIdNi6cwCAUyDYAc5h6rXB3/zS\n35Yc/O0vhRCRnt5KJhVfcpYQoq1/xWs//W77irW5wwfbV66deonm8b7+m1+obrcnEAokOxcn\nO4UQmYN7pRCBRMfiCy6f/tCm3QIAnC6CHeAcmtuz8sb3TW+JLx2o/+AOBFe/60+FEAs3Xzx9\ngVBnT+y4a11jfctmeAgAaFqcNwO0LtM07S4BANBIBDugRQ0ODu7atWt8fNzuQgAADUOwA1pU\ntVoVQnDLVABwEs6xA1rUwoULg8Fgd3e33YUAABqGHjugVRVz7vSorJbtrgMA0DAEO6BFZffv\nLh4+mHl9v92FAAAahqFYoEUFOhdIRQ11MhQLAM5BsANaVGRxv7ujxx+L2V0IAKBhGIoFAABw\nCIIdAACAQxDsAAAAHIJgBwAA4BAEOwAAAIcg2AEAADgEwQ4AAMAhCHYAAAAOQbADWtTI0NDB\nvXsz6bTdhQAAGoZgB7SofDZbKZfz2azdhQAAGoZbigEtqmfRotTkZEc394oFAOcg2AEtKhgK\nqZrm9XrtLgQA0DAMxQIAADgEwQ4AAMAhCHYAAAAOQbADAABwCC6eAOa3sdHRsZERVdMMXV82\nMBAIBiuVyv69ew1dtywrkUwuWrz4lR07pJRCiHAksnjJkpdffFFKaVmWlDIYDLpcHAcAwCE4\noAPz2+HBwQ3nnqsoSqVS0XXdNM1dr7yyfNUqv98vhBg8eNAwDMuyzl6/fuoltWr1nE2bisXi\n2OjowX37Vq5ZY1/5AIBGYigWmN+S7e27d+6cGBvzeDzhcHhifDzZ3l5PdUKI3r6+SrlcrVZf\nffnlPbt3SykNXXd7PPVnw5FIsVCwr3YAQIMR7ID5rbevb9ny5YZp7ti+3bKscqnkDwSmL+D1\n+dasW7dyzRpN01KTk/l8PhQO158yDENjHBYAHIRgB8xjRw4f1ms1l8vV1d2tqqpeqwWDwdzv\n7xI2MT4+NjKSnpz0ejxCCMuyXC5XsVCoBzvLsoaPHOnq6bFzBwAADcWXdWAeC0cir77yiqZp\nUspoLOb1+dq93mwm89ILLwgh/IHA0mXLMun0Sy++6HK5vD5fNBYbOnJkcmJidHi4VqtFYrH2\njg67dwIA0DAEO2Aei0QiazdseFOToixbvnx6Q1si0ZZITD1cuXp1/YdisVgul2e/RgDA3GEo\nFgAAwCEIdgAAAA5BsAMAAHAIgh0AAIBDEOwAAAAcgmAHtKhKZrJwaL9eKtpdCACgYQh2QItK\n7d1VGNyXHtxvdyEAgIZhHjugRfkT7aZhBNqSdhcCAGgYgh3QouJLBnxdi0KxmN2FAAAahqFY\nAAAAhyDYAQAAOATBDgAAwCEIdgAAAA5BsANwugaf/rWU0u4qAAAnRbADcAIHn3xi9NWX6j8P\nPvMbs1YtZ9K+aFxRlHIm9frTv9bLJXsrBAAcj2AH4MRKk+PSsoQQZq2qebypg3valiwTQhTH\nR+O9S0vpSbsLBAAci2AH4FhmrebyeKILe7NHBi3TUDWXqdekZbm8vsyhg+GuHl8kWibYAUDz\nYYJiAMcqZ1K+aDy2aMngU7/yBIK+aCx9cF/b4n7L0Cf27fb4A0IIyzTsLhMAcCx67AAcq5xN\n+6NxVdO80Vjq9X2+WFsxNRFItI/uennRuef1br6od/NF0uIqCgBoOgQ7AMeqZFK+WFwIkVw6\nkH59v1GphDu6q4WcXix4Q5H6MprHUysWbC0TAHAshmIBHGvR2y6o/+AOBM/+4/dNtfduvuiN\nnzddONdlAQBOhR47AAAAhyDYAQAAOATBDgAAwCEIdkCr0qvW5IiwLLvrAAA0DMEOaFGl3z5s\nPPrd8ivP2l0IAKBhCHZAS5PMRgcADsJ0J0CLClxwTems9YElZ9ldCACgYeixA1qV26O0dQhF\nsbsOAEDDEOwAAAAcgmAHAADgEAQ7AAAAhyDYAQAAOARXxQLzUurH97vbe8Kbtggh0j/7dvSS\na3O/+Vlw/Xn5p34hTVNRVWma0Uuvd8US+vhwYdsTkQuv0cJRIURq6/2xy25QA6H6eso7nq4N\n7hGWFb7gKk/nQjt3CQBwxgh2wHxVGx6UlqmomlUuqr6AkR53tbWb+Wz7ez8hhKge3l988cno\nJddWDx/wr9igjx2pBzurlJ9KdcKypF5L3PjnteHXy7teINgBwHzHUCww/1iVsuoN+Jetqex5\nWRq64vZIyxRCyFpN9frryxipcS0cK+/e4Vt8ljvRURs7Un+h4g28sSJVDWy8SAhR2feqr3+V\nDXsCAGgoeuyA+UcfH3K3d/mXr03/5AEtEncnOo3UuCvero8PmflMauv9QhGuaCK08eLU1vu0\ncFRIIY3a71/YPX1VVqmQ2faEv3+ld+FSm/YGANAwBDtg/tHHhz2dCxSX25XoLL36vG/xcn18\n2N3erY8Phzdf5l+xvr5Y7rc/j115kyuWFEJM/vC/6y90t3dNrUemxgoHd8YufqfqC5xwQwCA\n+YWhWGD+0SdGXMluIURw7abyzu3u9q56V5w+PuzuWFBfxshMmPl0PdUJIVSv38yl9fHh4gtP\npR76Ruqhb8haxfjVQ7JczPzi+5lfft+2nQEANA49dsD8E7/q3fUftHCs+7bPCCGiF18rhPD0\nLJ5axhVLxq++5Y2XXH2zECL+jvdMtRSLRfeNt0ZjMZeL4wAAOAQ9dgAAAA5BsAMAAHAIgh0A\nAIBDEOwAAAAcgmAHAADgEAQ7AAAAhyDYAQAAOATBDgAAwCGaaGJSKaVhGLqu211IU5NS8had\nUv1dUhTF7kKammVZQgjDMKSUdtfS1EzTtCyL37uZ1T9OvFGng8P46agfl+bmSO6wf44mCnb1\nI0K1WrW7kKYmpeQtOiUpZa1Ws7uKZmeaphBC13XDMOyupamZpsnv3SnV/wxblsUbdTp4l06p\n/lVhbo7kDvt70UTBTtM0v98fCoXsLqR5VatVVVV5i04pk8kEg0F67GZWLBYNw/D7/dxSbGaV\nSsWyrEAgYHchTc00zWq16nK5OECdUq1W4106pWw2a1nW3BzJHdZjxzl2AAAADkGwAwAAcAiC\nHQAAgEMQ7AAAAByCYAcAAOAQBDsAAACHINgBAAA4BMEOAADAIQh2AAAADkGwAwAAcAiCHQAA\ngEMQ7AAAAByCYAcAAOAQBDsAAACHINgBAAA4BMEOAADAIQh2AAAADkGwAwAAcAiCHQAAgEMQ\n7AAAAByCYAcAAOAQBDsAAACHINgBAAA4BMEOAADAIQh2AAAADkGwAwAAcAiCHQAAgEMQ7AAA\nAByCYAcAAOAQBDsAAACHINgBAAA4BMEOAADAIQh2AAAADkGwAwAAcAiCHQAAgEMQ7AAAAByC\nYAcAAOAQBDsAAACHINgBAAA4BMEOAADAIQh2AAAADkGwAwAAcAiCHQAAgEMQ7AAAAByCYAcA\nAOAQBDsAAACHINgBAAA4BMEOAADAIQh2QIsyJ8esfS9blZLdhQAAGsZldwEA7FF99jFzZLBi\nGZ7zrrC7FgBAYxDsgBalxtuVUkFr67C7EABAwxDsgBbl23yZLOS9iaTdhQAAGoZz7IAWpvHV\nDgAchWAHAADgEAQ7AAAAhyDYAQAAOATBDgAAwCE4dRqYr8q7X6zs3qG43Fa1HLnkWldbx/Et\n6Ye+IS1TKKrUa7Gr3qVF4lMtpl4TGy8VsVjh6UdrI4ekYQTXbfYtW2P3bgEA/nAEO2C+Kj7/\nm+SffFwoqpnLWJXiCVusSjlx80eEELlf/0QfH9YicTOfSb7/k0KI3J5XKru2iwW97s6Foc2X\n1Q7vr+x9hWAHAPMawQ6Yr3z9qzI//66vf5V3yXItEju+xaqUzXwm85MHzELWlezy9a+yKmXF\n56+/3MpMimBYqKq7c8Hkg/+phaLRK2+ydYcAAGeKc+yA+Sq0aUvk0utkrZr67j3SMI5v0ceH\nAmdvir3zTxI3f8TMpsxiTh8fMnOZ9I/uS//oPjOfVldsEEKo/mDi5o8KIcxc2uZdAgCcGYId\nMC8VX3zSKhVUr9+/6hzF5bLKheNbjLFhV+LoHcOkriuKaowNh99+efz6D8Sv/4Bv8xVC1Qq/\n+alVzAshFLdbSmnrPgEAzhRDscC85OlcmHn4QcXjlabpWdCnhWPHt+jjQ5UDu8o7n5e1qm/g\nbDUQ0seHQpsunb4e3/L1mYcfVLw+LdrmTnbZtDcAgMYg2AHzkrtrUdtNfzlzS+zqm4951fEt\nrvbutnd9aDYqBADMPYZiAQAAHIJgBwAA4BAEOwAAAIcg2AEAADgEwQ4AAMAhCHYAAAAOQbAD\nAABwCIIdAACAQxDsAAAAHIJgBwAA4BAEOwAAAIcg2AEAADgEwQ4AAMAhCHYAAAAOQbADAABw\nCIIdAACAQxDsAAAAHIJgBwAA4BAEOwAAAIcg2AEAADiEy+4C3mRwcNDv99tdRfPKZrOqqobD\nYbsLaXb5fD4UCimKYnchTa1cLler1XA4rGma3bU0tVqtZlmWz+ezu5CmZllWLpfzeDyBQMDu\nWppdLpeLRCJ2V9HsCoWCYRhjY2NzcCTP5/OzvYm51FzB7p//+Z/tLgEAALQcVXXIGGYTBbuL\nL764s7PT7iqa2n333ReNRq+//nq7C4ETPPvsszt37rzmmjs3Pv4AAAldSURBVGs6OjrsrgXz\nXi6X+/73v9/f33/hhRfaXQuc4Oc///nw8PD73vc+t9s9N1vcsmXL3GxotilSSrtrwOk6//zz\nlyxZcv/999tdCJzgC1/4wje/+c177rln7dq1dteCeW9wcPCmm2669tprP/OZz9hdC5zgr/7q\nr5555plf/epXDO6/VQ7peAQAAADBDgAAwCEYip1Pnn32Wb/fv2bNGrsLgRMMDg6OjIysWrUq\nFArZXQvmvUqlsmPHjmQyuXTpUrtrgRPs2rUrl8ude+65jrmmYc4Q7AAAAByCIAwAAOAQBDsA\nAACHaKJ57HBC1Wr1rrvuEkKYpnnbbbdNzVf+wgsvfPWrX12wYIEQ4o//+I9XrVplZ5WYD072\nWTpZOzADDk1oLCnlj3/8429/+9uf//znk8nkVDsHqLdK+/SnP213DZjJww8/3NXV9Rd/8Rcu\nl+vJJ59ct25dvX3v3r09PT0f/vCHL7744vb2dnuLxLxwss/SydqBGXBoQmPl83nDMFKp1Lnn\nnjt97joOUG8VQ7HNbv/+/cuWLRNC9Pf379u3b6q9UChs27btC1/4wl133ZXNZu0rEPPGyT5L\nJ2sHZsChCY0ViUQ2bNhwfDsHqLeKodgmtXXr1m3bth1zSwCv1zv186ZNm84555yOjo6nn376\nm9/85sc//vE5rxHz2PTP0um0AzPg0IS5wQHqdNBj16Suu+66T3/60zfddNOyZcv27NkjhNi9\ne3d/f//UAocPH65Wq0KIQCCg67pthWL+ONln6WTtwAw4NGFucIB6q5jHrtnVarW7777bNE0h\nxP/f3r2HNNXGcQD/HS9JiWYKqYSzjOZQ0zEt1GY703AWxazE0sTMBKcm0T8VdPsrKtEmhRcq\n0yRGChohOpDJZgkrLZiGornAMlIh1jSdhrf3jwN7Rd+8Ja81vp+/znnO8zvPc8Y4fHl2Obm5\nud+/f3/8+PH169f7+/sfPnzo7Ow8Pj6uUCi8vLzWe6bwp/vVe2leu7Oz83rPFP4CuDXB2jIa\njVVVVR8+fNixY4dQKAwLC8MNanUQ7AAAAABsBD6KBQAAALARCHYAAAAANgLBDgAAAMBGINgB\nAAAA2AgEOwBYmcHBQesDf4aHh0UiUU1NzeDgoIODg0AgEAgE/v7+oaGhOp3OWjIwMODg4HDn\nzh1rS11dXWRkZERExO7du0+dOjU6OrrkuD9+/MjNzeXxeP7+/kFBQaWlpWtyOSEhIVqtltu+\nd+/e3P+15/P57969m9e/ubn55MmTC8/T0tJy4MCBeY0qlWpNJgkAsEwIdgCwShMTE3K5PDs7\n+/jx40Tk5ubW3d3d3d3d09NTVFSUmJg4MzPD9ayoqAgICCgvL7cWpqen19TU6PX6jo4OFxeX\nJ0+eLDmcXC7fuHHjx48fe3p61Gp1UVHRixcvfv8q4uLiNBoNt93U1DQyMjI0NERE/f39ZrNZ\nJBLN6y+RSJ49e7acM1ssloKCgt+fIQDA8iHYAcBqTE9PJyUlHTp0KCMjY+HR8PDwyclJLiER\nUXl5uVKpHBsb0+v1RDQ2NjY+Ps7FPoZhSktLc3JyFh/u9evXX758uX37tqOjIxH5+Pi8efNG\nLpdzR/Py8liWjY2NzczM/Pnzp1arlclk586dS05Olkql3HLgvD7WM8tksqamJu6KOjs7k5KS\nuF2NRhMbG8swzLxCjUZz+PDh6enp9PT08PDwo0ePpqWlcZl1eno6JSVl7969LMuOjo5mZmb2\n9vaePn36d19rAIBlQ7ADgNVQKBRDQ0MXL178z6PPnz/38PDw9PQkoubm5qmpqejo6NTUVC4A\neXh4XLt2LSgo6MiRI4WFhZ8+fVpyuM7OzrCwMDu7f29Z1seEv3z5UqfT6XS6xsZGOzu7yspK\nR0dHg8GQl5enUqm8vb3VavXCPtbziMXi7u7u4eHhtrY2oVDIsiy3gKfRaGQy2a8K6+vr+/r6\n9Hp9WVmZWq1mGIaI2tvbCwoKWltbt2zZUl9ff+XKFT8/v+UsRgIArBUEOwBYMZPJ5OvrS0SP\nHj2yNprNZqFQKBQKvby8lEplXV0dl8PKysrS0tIYhklPT6+urrZYLER06dKlvr6+s2fP9vb2\nikSip0+fLj6ivb391NQUt11SUrJ9+3YfH5/4+Hgi0ul0XV1dLMuyLNvW1vb582ciCggI4JLf\n1q1bTSbTf/bhbNiwISoqSqfTNTU1xcTEiMXiV69eEZFWq42Njf1VYUdHR1RUFMMw7u7u0dHR\nXGNwcDCXZX19fc1m81q80gAAK+Ow3hMAgL+Pu7v71atXz5w5ExERERgYGBERQURubm4Gg4GI\nqqqqiouLBQIBEY2MjNTW1vJ4vNraWiKampqqra1NTk4eGRlxc3OLj4+Pj48/ePDgrVu3UlJS\nuJMrlcqioiIiUqvVu3bt4hqDg4Nv3LgxOTnp6OiYlZWVlZX19u3by5cvE5GTk1NCQkJ+fr51\nei0tLfb29tbd2dnZhX3miouL02q179+/Lykp2bRpk7e3d0NDg5eXl6en58JCbj1vZmaGW6Uj\nIutYDg7/3lHxUB8AWBdYsQOAVdq2bVtlZWViYuLXr1/ntp84ccLFxeX+/ftEpFKppFJpV1eX\nwWAwGAwPHjwoLy9Xq9USicRkMnH9DQaDn5+ftfzChQtGo9FoNFpTHRGJRKKQkJCcnJzx8XEi\nslgsNTU1rq6uRCQWixsaGiYmJogoPz9/4e9Yl+zDBbuBgQE+n09E0dHReXl5cXFxixQKBILW\n1lYiMpvN1h/VzmNnZzf3y3wAAP8DBDsAWD2WZc+fP3/s2LF5Caa4uPjmzZtGo7GsrCw7O9va\nnpCQ0NnZGRgYmJGRIZVKhUJhSEhIV1fX3bt3lxyrurra1dWVz+fz+fzQ0NCZmZmKigoi2rdv\nX2pqqkQikUgk7e3tQUFBC2sX77Nz506LxbJnzx5uNyYmprm5WSaTLVIol8udnJwiIyMVCkVU\nVNTctTorHo/HMMz+/fuXvDQAgLXC4PMCAICV+vbtW2NjY3Jy8uzsrFgsLiwstOZCAIB1hBU7\nAIAV27x5c0NDQ1hYWGRkJMuySHUA8IfAih0AAACAjcCKHQAAAICNQLADAAAAsBEIdgAAAAA2\nAsEOAAAAwEYg2AEAAADYCAQ7AAAAABuBYAcAAABgIxDsAAAAAGwEgh0AAACAjUCwAwAAALAR\nCHYAAAAANuIfc4eQPbcTsCoAAAAASUVORK5CYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "p"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Save Figure\n",
    "fig_file <- file.path('figures', 'classifier_coefficients_scatter.pdf')\n",
    "ggplot2::ggsave(fig_file, plot = p, dpi = 600, width = 4, height = 3)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R [conda env:multiple-myeloma-classifier]",
   "language": "R",
   "name": "conda-env-multiple-myeloma-classifier-r"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.4.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}