47_KNNplots.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'\\nparse and plot outputs from 46a/b_KNN\\n\\nShaina Lu\\nZador & Gillis Labs\\nMay 2021\\n'"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\"\"\"\n",
    "parse and plot outputs from 46a/b_KNN\n",
    "\n",
    "Shaina Lu\n",
    "Zador & Gillis Labs\n",
    "May 2021\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from __future__ import division\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import scipy as sp\n",
    "from scipy import stats\n",
    "from sklearn import metrics\n",
    "from sklearn.linear_model import Lasso\n",
    "from sklearn.linear_model import LinearRegression\n",
    "from sklearn.preprocessing import StandardScaler\n",
    "from sklearn.model_selection import train_test_split #stratify train/test split\n",
    "import random"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def set_style():\n",
    "    #function for matplotlib formatting\n",
    "    plt.style.use(['seaborn-white','seaborn-notebook'])\n",
    "    plt.rcParams['figure.figsize'] = [6,4]\n",
    "    plt.rcParams['axes.spines.top'] = False       #remove top line\n",
    "    plt.rcParams['axes.spines.right'] = False     #remove right line\n",
    "    plt.rcParams['axes.linewidth'] = 2.0          #set weight of axes\n",
    "    plt.rcParams['axes.titlesize'] = 20           #set font size of title\n",
    "    plt.rcParams['axes.labelsize'] = 18           #set font size of x,y labels\n",
    "    plt.rcParams['axes.labelpad'] = 14            #space between labels and axes\n",
    "    plt.rcParams['xtick.labelsize'] = 14          #set x label size\n",
    "    plt.rcParams['ytick.labelsize'] = 14          #set y label size\n",
    "    plt.rcParams['legend.fontsize'] = 16          #set legend font size\n",
    "    \n",
    "set_style()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ST KNN all areas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "STTRAINFILE = \"/home/slu/spatial/data/multiclass_knn/STknnAUROC_train_051121.npy\"\n",
    "STTESTFILE = \"/home/slu/spatial/data/multiclass_knn/STknnAUROC_test_051121.npy\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "STtrain = np.load(STTRAINFILE)\n",
    "STtest = np.load(STTESTFILE)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAE7CAYAAADtredSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3dd1RU1/428GdQkICggtj1IuhAjFS7ggWxxBLAXrHXaDSxERMhURRN9GpEI4oitqtGEKyxNxQLKHbsEoGoFKUISj3vH77Mz3FABpw5tOezlms5p373kPC4z9lnH4kgCAKIiIhEolHSBRARUcXC4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVJVLugAiotLu9OnT8PPzw927d5Gbm4smTZpg2LBhcHFxUfoYMTEx6Nq1a4Hre/XqhZUrVyosj46Oxrp163DhwgUkJiaiWrVqaNOmDb799luYmprKbZuVlYUzZ87g9OnTuHnzJv7991/k5OSgUaNGcHR0xLhx41C1alXlG64mEj5ASkRUsO3bt2PRokWoXr06evXqBU1NTRw9ehQvXrzA2LFjMW/ePKWOkxc85ubmcHR0VFjftGlT9OzZU27ZnTt34Orqijdv3qBt27Zo1qwZXrx4gePHj0NTUxObN2+GtbW1bPvHjx+jV69e0NHRQZs2bWBqaor09HScP38ez549g7GxMXbu3AkDA4PP+1I+l0BERPmKjo4WmjdvLrRu3VqIjo6WLU9KShIcHR0FqVQqXLt2TeljSaVSYd68eUqf38nJSZBKpcLmzZvlll+7dk1o1qyZ0L17dyEzM1O2/MWLF8L27duFtLQ0ue0zMjKEiRMnClKpVFi4cKHS51cX3uMhKuNiYmJgZmYGNzc3PHv2DN999x3atGkDGxsbjB07Fg8ePAAAvHr1CgsWLICdnR0sLCzQv39/XLp0Kd9jZmdnY8eOHRg0aBBsbW1hZWUFZ2dnbN++Hbm5uQrb7927F9OnT0fXrl1haWkJW1tbDBkyBPv27cv3+CNHjoSZmRmys7Ph4+OD7t27o3nz5ujUqRN+//13ZGZmqu4L+gyBgYHIzMzE8OHD0aBBA9nyatWqYdKkSQCAXbt2qeXc0dHRiIyMhKGhIVxdXeXW2djYoGvXroiKikJISIhsee3atTF8+HDo6OjIba+lpSWr98qVK2qptyh4j4eonIiNjcXAgQNhamoKFxcXxMbG4vjx4xg5ciR2796N8ePHo2rVqvj666+RnJyMw4cPY8KECTh69Cjq1asnO05WVhYmT56M8+fPo3HjxujTpw+qVKmCy5cvY9GiRbhx4wZ+//13uXP/8ssvaNKkCVq1agUjIyMkJSXh7NmzmDt3Lp4+fYqZM2fmW/OsWbNw9epV2Nvbo1OnTjh37hw2btyIV69ewcvLS63flzLygtne3l5hXceOHeW2UVZcXBx27dqFpKQkVK9eHdbW1jA3N1fYLj4+HgBQv359aGgo9hHygvDixYtwcHAo9LyVK7//dV+pUqUi1asODB6icuLKlSuYOXMmpkyZIlu2du1arF69GoMGDULPnj3xyy+/yH6JtW/fHvPmzYO/vz/mz58v28fHxwfnz5/HiBEjMH/+fNkvqpycHCxYsACBgYHo0aOH3H2KgwcPolGjRnL1ZGZmYsKECfD19cXQoUNRu3ZthZqjo6Nx8OBBVK9eHQDw/fffw8nJCcHBwfjhhx9gZGRUaLsjIyNx4sSJInxTwKhRo6Cvr1/odk+fPgUAGBsbK6yrVasWdHR08OLFC7x9+xZffPGFUue+cOECLly4ILesdevWWLZsmdw/AGrUqAEA+PfffyEIAiQSidw+MTExAIAnT54odd7AwEAA+Yeo6Er6Wh8RfZ68ewddunQRsrOz5dbFxsYKUqlUsLKyElJTU+XWZWdnC82aNRNGjBghW5aTkyO0bt1a6NChg5CVlaVwruTkZMHMzEz47rvvlKrt6NGjglQqFYKCguSWjxgxQpBKpcKFCxcU9vnjjz8EqVQqnDp1SqlzBAYGClKptEh/Prxf8ylfffWVIJVK8/0uBEEQ7OzsBKlUKrx8+bLQYyUkJAirVq0Sbt++LSQnJwvJycnClStXhJEjRwpSqVTo1q2bwr2Z7t27C1KpVNiyZYvc8uvXrwvNmjUTpFKpMGDAgELPfeLECcHMzEzo2LGjkJSUVOj26sYeD1E58eWXXypcRqlVqxaA9/9i/3gYbaVKlWBoaIiXL1/Klj19+hRJSUkwNjbGunXr8j2Ptra2wr+y//33X/j6+uLixYt4/vw53r17J7f+w3N8qHnz5grL6tatCwBITk7Od5+P9evXD/369VNqW1UT/v+g4I97I/kxNDTEjBkz5Ja1atUKfn5+GDZsGG7cuIE9e/Zg1KhRsvULFy7E+PHjsXjxYpw+fRrm5uZ4+fIljh07BlNTU9y/f7/QS2fXrl3D7Nmz8cUXX8Db2xvVqlUrRktVi8FTCDMzMwDA/fv3S7gSok/T09NTWJZ3XT+/dXnrs7OzZZ+TkpIAAFFRUVizZk2B50pLS5P9PTo6GgMGDEBKSgpatmwJOzs7VK1aFZUqVUJsbCyCgoIKHCyQ3+WuvF+k+Q1iEFvVqlXx+vVrpKamyi59fSjve/icZ2MqV66MgQMH4saNGwgPD5cLnjZt2mDPnj1Yt24dwsLCEBYWhjp16mDKlCkwNzfH1KlTPzk0OiIiAhMmTIBEIoGvry8sLS2LXacqMXiISCYvoLp16/bJ4PnQ5s2bkZSUBC8vL4Wex8GDBxEUFKTyOj+kzns8jRs3xuvXrxEVFaUQPHFxcUhPT0edOnWUvr9TkLxjp6enK6wzNzfHH3/8obB89erVAAALC4t8jxkeHo6JEydCQ0MDmzZtknvep6SVueDx8fHBypUrMXz4cLi7uwN4391ds2YNdu/ejZSUFFhZWcHd3R1NmzaV7ZecnAxPT0+cOnUKAODg4IAFCxYo9R8fUUVhYmICfX19XL9+HVlZWdDU1Cx0n3/++QcA0L17d4V1YgzdjYyMVDok87i4uCj1/37btm1x7do1hISEwMbGRm7duXPnZNt8rhs3bgAAGjZsqNT2mZmZCA4OhoaGBnr37q2w/uLFi5g6dSo0NTWxcePGUtPTyVOmnuO5fv06/vrrL9nlrzy+vr7w8/PDggULEBAQAAMDA4wZMwZv3ryRbTNr1izcvXsXvr6+2LhxI+7evYu5c+eK3QSiUq1y5coYMWIE4uPj4enpqXCvBnj/L/1Hjx7JPtevXx+AYsiEhIQgICBAvQXj/T2e+/fvF+nPh8/kFHZsLS0t7NixQzaKDHj/D9n169cDAIYMGSK3T2pqKh4/foy4uDi55Tdu3Mj3kuPFixfh7+8PAPjmm2/k1qWnpyMnJ0duWVZWFjw8PBAbG4uhQ4cqjCY8f/48Jk+eDC0tLfj7+5e60AHKUI8nNTUVs2fPxuLFi/Hnn3/KlguCgK1bt2LixIno0aMHAGDZsmVo164dDh48iCFDhuDx48cICQnB//73P9ja2gIAfv31VwwfPhxPnjyBiYlJibSJqDSaOnUq7t27h127duH06dNo27YtateujcTERPzzzz+4du0avv/+ezRp0gQAMGzYMOzduxczZsxA9+7dUbt2bTx8+BAhISH4+uuvcfjw4RJuUfE1bNgQc+fOhaenJ/r375/vlDkf94SOHz+OH3/8ES4uLli6dKls+fLly/Hw4UO0bt0aderUAfD+3nHec0AzZsyQ/X7Kc/nyZfz8889o164d6tatizdv3uDs2bOIjY1F586dFabrefLkCaZOnYqMjAx06tQJJ0+exMmTJxXaNX36dJV8P8VVZoJnwYIF6NGjB9q1aycXPDExMYiPj0eHDh1ky7S1tdGqVStERERgyJAhiIiIgI6OjtwPtUWLFtDR0UFERASDh+gDmpqa+PPPP7Fv3z4EBQXhzJkzSE9PR40aNdCgQQPMmDEDffv2lW1vbm6OrVu3YtWqVTh37hyys7Nhbm6ONWvWQE9Pr0wHD/B+loX69evDz88PwcHBEAQBpqammDlzZpEmCf3mm29w4sQJ3L59GyEhIcjKykLNmjXx9ddfY8SIEWjZsqXCPsbGxrC1tUVYWBgSExOhra0Nc3NzTJs2Dc7OzgoPlsbHxyMjIwMAcPToURw9ejTfWko6eMrEJKF//fUXdu3ahV27dkFLSwsjR45E06ZN4e7ujmvXrmHo0KE4ffq03MNXP/74I+Li4rBp0yb4+Phgz549CsnftWtXDBo0SDaVRH44qo2ISLVKfY/nyZMn+O9//4sdO3ZAS0urwO0KG0ef33ohn6eBiYhIvUp98Fy/fh2vX7+W69rn5OQgLCwMu3btwsGDBwG872LmPXgGAImJiahZsyYAoGbNmkhMTJQLGkEQ8Pr1axgaGorYGiIiKvXB4+joqPB0848//ghjY2NMmjQJjRs3hpGREUJDQ2WjNzIyMhAeHi4btWZjY4P09HRERETI7vNEREQgPT1d4cYgERGpV6kPHn19fYXx9jo6OqhWrRqkUikAwNXVFT4+PjAxMZFN9aGjo4M+ffoAAExNTWFvbw8PDw8sWrQIgiDAw8MDXbp04cACIiKRlfrgUcaECROQkZGBhQsXIjk5GVZWVvDz85ObxmL58uXw9PTE2LFjAbx/gDTvAVQiIhJPmRjVVpI4qo2ISLXK1MwFRERU9jF4iIhIVAweIiISFYOHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETF4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETF4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETF4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRlfrg2bFjB/r27QtbW1vY2tpi8ODBOHPmjGy9IAjw9vaGnZ0dLC0tMXLkSDx8+FDuGMnJyZgzZw5atGiBFi1aYM6cOUhJSRG5JUREBJSB4KlduzZmz56NoKAgBAYGom3btvj2229x7949AICvry/8/PywYMECBAQEwMDAAGPGjMGbN29kx5g1axbu3r0LX19fbNy4EXfv3sXcuXNLqklERBWbUAa1atVK2Llzp5Cbmyt06NBB+PPPP2Xr3r59K1hbWws7d+4UBEEQHj16JEilUiE8PFy2TVhYmCCVSoXHjx8Xei6pVCpIpVLVN4KIqIIq9T2eD+Xk5ODQoUNIT0+HjY0NYmJiEB8fjw4dOsi20dbWRqtWrRAREQEAiIiIgI6ODmxtbWXbtGjRAjo6OrJtiIhIPJVLugBl3L9/H0OGDEFGRgZ0dHSwZs0amJmZ4dq1awCAmjVrym1vaGiIuLg4AEBCQgIMDAwgkUhk6yUSCQwMDJCQkCBeI4iICEAZCZ7GjRsjODgYKSkpOHbsGObNm4dt27bJ1n8YKvnJb70gCIXuR0REqqd08Lx+/RqvXr2CqampbFl0dDT8/f2RlJQEZ2dn2Nvbq6VILS0t/Oc//wEAWFhY4NatW/D398eUKVMAAPHx8ahbt65s+8TERFkvqGbNmkhMTJQLGkEQ8Pr1axgaGqqlXiIiKpjS93gWL16MefPmyT6npaVh+PDh2LFjBw4dOoRJkyYhLCxMLUV+LDc3F5mZmWjQoAGMjIwQGhoqW5eRkYHw8HDY2NgAAGxsbJCeni53PyciIkJ2n4iIiMSldPBcv34dHTt2lH0+fPgw4uLisGHDBoSEhMDU1BQbN25UeYHLly9HeHg4YmJicP/+faxYsQJXrlxB3759IZFI4Orqig0bNuDYsWN48OAB3NzcoKOjgz59+gAATE1NYW9vDw8PD1y/fh0RERHw8PBAly5dYGJiovJ6iYjo05S+1JaYmCh3OSskJATNmzeXhZGLiws2b96s8gITEhIwZ84cxMfHQ09PD2ZmZvD19ZVd1pswYQIyMjKwcOFCJCcnw8rKCn5+fqhatarsGMuXL4enpyfGjh0LAHBwcIC7u7vKayUiosIpHTyVK1dGRkaG7POVK1fg4uIi+6ynp4ekpCTVVgdg6dKln1wvkUgwffp0TJ8+vcBtqlevjuXLl6u6NCIiKgalL7UZGxvj6NGjEAQBJ0+eRHJyMtq1aydb/+LFC1SrVk0tRRIRUfmhdI9n+PDhcHNzQ6tWrfDu3Ts0bNhQLnjCwsJgZmamliKJiKj8UDp4nJ2dIZFIcOLECVStWhWTJ0+GpqYmgPdDrd+8eYOhQ4eqrVAiIiofJIIgCCVdRGmW14u7f/9+CVdCRFQ+lKm52oiIqOwr0pQ52dnZOHHiBG7cuIGUlBTk5ubKrZdIJFiyZIlKCyQiovJF6UttSUlJcHV1xcOHD2XTz+Ttmvd3iUSCyMhItRYsNl5qIyJSLaUvta1atQpPnjyBp6cnjh8/DkEQsGnTJhw+fBi9e/eGhYUFLl++rM5aiYioHFA6eM6ePQtnZ2f0799fNiuAhoYGTExMsHz5cmhra2PFihVqK5SIiMoHpYMnPj4eFhYWAN7PYgAAmZmZsvVdu3bFqVOnVFweERGVN0oHT/Xq1fH27VsAgK6uLipXroznz5/L1mtqaiIlJUX1FRIRUblSpClzHj169H4nDQ00a9YMQUFByMzMxNu3bxEcHIyGDRuqrVAiIioflA6eDh064OjRo7LLa6NHj8aNGzfQunVrtG/fHrdv38aoUaPUVigREZUPSg+nFgQBWVlZ0NLSki07duwY9u/fDw0NDfTs2RO9evVSW6ElhcOpiYhUi1PmFILBQ0SkWsWaMiczMxMvX76UG9VGRESkjCIFz507d+Dq6gpbW1t07twZV69eBfD+7aSjRo1CaGioWookIqLyQ+ngiYyMxPDhwxEdHQ0nJye5dYaGhsjIyEBQUJDKCyQiovJF6eD5448/UKtWLRw8eBCzZs3Cx7eG2rZti5s3b6q8QCIiKl+UDp6rV69i4MCB0NXVhUQiUVhfr149xMXFqbQ4IiIqf5QOnoyMDOjp6RW4/s2bNyopiIiIyjel38fTqFEj3Llzp8D1ly5dQpMmTVRSVHli7HaoRM4btbR3iZyXiKgwSvd4+vTpg3379smNXMu75Obn54eQkBCFQQdEREQfU7rHM3bsWFy4cAHjxo2DiYkJJBIJvLy88OrVKyQkJKB9+/YYNmyYOmslIqJyQOkej5aWFjZv3ox58+ahSpUqqFKlCqKiolCjRg3MmTMH69evh4ZGsZ5HJSKiCkTpHg/w/j08o0ePxujRo9VUDhERlXdKdVHS0tLg6OgIf39/NZdDRETlnVLBo6uri6SkJOjq6qq7HiIiKueUviljZWWFW7duqbMWIiKqAJQOntmzZ+PIkSMIDAxUmC6HiIhIWUoPLvDy8oK+vj5+/vln/P7772jUqBG0tbXltpFIJNiyZYvKiyQiovJD6eCJiYkBANStWxcAkJCQoJ6KiIioXFM6eE6dOqXOOoiIqILgE59ERCQqBg8REYmqSDMXPHv2DP7+/rhx4wZSUlKQm5srt14ikeDEiRMqLZCIiMoXpXs89+/fh4uLC/bs2YOsrCxER0dDR0cHGRkZiI2NRaVKlWQDD4iIiAqidPCsXr0ampqa2Ldvn2zqnPnz5+P8+fNYuHAhUlJS4OHhoa46iYionCjSq68HDx4seyXChwYNGoSOHTti+fLlKi+QiIjKF6WDJy0tDQ0bNgQAaGpqAgDS09Nl621tbXHt2jUVl0dEROWN0sFTs2ZN2UOjVatWxRdffIGoqCjZ+pSUFOTk5Ki8QCIiKl+UHtVmbm4uN0lo69atsXXrVlhaWiI3Nxfbt2+Hubm5WookIqLyQ+keT9++fZGUlIR3794BAGbMmIHU1FS4urpi9OjRSE1Nxffff6/yAtevX4/+/fvD1tYWbdu2xeTJk/HgwQO5bQRBgLe3N+zs7GBpaYmRI0fi4cOHctskJydjzpw5aNGiBVq0aIE5c+YgJSVF5fUSEdGnSYTPmGr6+fPnOH78OCpVqoSOHTvK7gGp0rhx49CrVy9YWFhAEASsXr0a169fx6FDh1C9enUAwIYNG7Bu3TosXboUjRs3xtq1a3H16lUcOXIEVatWBQCMHz8ez58/x6JFiyCRSPDzzz+jYcOG8PHx+eT5zczMALwfTl4cxm6HirXf54pa2rtEzktEVJjPCp6SkJaWhpYtW2Lt2rVwcHCAIAiwt7fH8OHDMWXKFADAu3fv0K5dO8ybNw9DhgzB48eP0atXL/zvf/9DixYtAADh4eEYPnw4/v77b5iYmBR4PgYPEZFqlbkpc9LS0pCbmwt9fX0A72fNjo+PR4cOHWTbaGtro1WrVoiIiAAAREREQEdHB7a2trJtWrRoAR0dHdk2REQkjjIXPIsXL8aXX34JGxsbAEB8fDyA96PuPmRoaCgbhZeQkAADAwO5548kEgkMDAz4egciIpEVaa62kubl5YWrV69i586dqFSpkty6jx9q/Vh+6wVBKHQ/IiJSrTLT41myZAkOHTqELVu2yA1iMDIyAvB/PZ88iYmJsl5QzZo1kZiYKPfKbkEQ8Pr1axgaGopQPRER5SkTwePp6YmDBw9iy5YtMDU1lVvXoEEDGBkZITQ0VLYsIyMD4eHhsstxNjY2SE9Pl7ufExERgfT0dNk2REQkjlJ/qe3XX3/Fvn37sHbtWujr68t6Njo6OtDV1YVEIoGrqyt8fHxgYmICY2NjrFu3Djo6OujTpw8AwNTUFPb29vDw8MCiRYsgCAI8PDzQpUuXT45oIyIi1SvWcOq3b98iKSkJ+e1ar149lRSWJ28488emTZuG6dOnA3h/2WzNmjXYvXs3kpOTYWVlBXd3d0ilUtn2SUlJ8PT0lL3C28HBAe7u7rLRcYWdn8OpiYhUQ+ngyc3NxcaNG7Ft27ZPjgSLjIxUWXGlAYOHiEi1lL7Utnz5cvj5+aFp06bo0aOHbNYAIiKiolA6ePbv3w97e3v4+vqqsx4iIirnlB7VlpKSgq5du6qzFiIiqgCUDh6pVKrwrAwREVFRKR0806ZNw65du/D8+XN11kNEROWc0vd4bt++jXr16qFXr17o1q0bGjRoAA0N+dySSCT49ttvVV4kERGVH0oHz5o1a2R/379/f77bMHiIiKgwSgfPyZMn1VkHERFVEEoHT/369dVZBxERVRBlYpJQIiIqPwrs8axZswYSiQRTpkyBhoaG3D2egvAeDxERFabAudrMzc0hkUhw48YNaGlpwdzcvPCDSSScq+0jnKuNiEhegT2evMEEWlpacp+JiIg+R4HB8/FgAg4uICIiVeDgAiIiElWR30B669Yt3Lx5E8nJycjNzZVbx8EFRERUGKWD5927d5g2bRouXLgAQRAgkUhkbyDN+zuDh4iICqP0pba1a9fiwoULmDx5MrZu3QpBELB06VL4+vqiZcuWsLCwwKFDJTOCi4iIyg6lg+fo0aPo2bMnZsyYgaZNmwIAateuDXt7e2zevBlZWVkICgpSW6FERFQ+KB08z58/R6tWrQAAlSpVAgBkZWUBACpXrozevXuzx0NERIVSOnh0dXWRk5Mj+7uGhgbi4uJk6/X09JCQkKD6ComIqFxROngaNWqEqKgoAO97PE2aNMHRo0cBAIIg4Pjx46hbt65aiiQiovJD6eBp164djh49Kuv1DB48GCEhIXB0dET37t0RGhqK/v37q61QIiIqH5QeTj1x4kQ4OTnJhlAPHz4cmZmZ2L9/PzQ0NPD9999jwoQJaiuUiIjKhwInCaX3OEkoEZFqKXWpLS0tDY6OjvD391dzOUREVN4pFTy6urpISkqCrq6uuushIqJyTunBBVZWVrh165Y6ayEiogpA6eCZPXs2jhw5gsDAQPC2EBERFZfSo9q8vLygr6+Pn3/+Gb///jsaNWoEbW1tuW0kEgm2bNmi8iKJiKj8UDp4YmJiAED2kChnKSAiouJQOnhOnTqlzjqIiKiC4BtIiYhIVEV+A2lmZiYuX76M6OhoAEDDhg3RunVrVKlSReXFERFR+VOk4AkODoaXlxdSUlLk3j6qr6+PefPmoV+/fmopkoiIyg+lg+fw4cNwc3NDvXr1MG7cOJiamkIQBDx+/Bi7du3CTz/9BG1tbfTq1Uud9RIRURmn9Fxt33zzDbKzs/HXX3+hatWqcutSU1MxcOBAaGpq4sCBA2optKRwrjYiItVSenDB06dP0a9fP4XQAd6/BK5fv374559/VFocERGVP0oHj5GR0SdnLNDQ0EDNmjVVUhQREZVfSgePi4sLgoKCkJaWprDuzZs3CAwM5OACIiIqVIGDC8LCwuQ+t2zZEqdPn0bfvn0xbNgwmJiYQCKR4NGjR9i5cydq1KiBFi1aqL1gIiIq2wocXGBubg6JRCK37MNN89Z9vCwyMlIddZYYDi4gIlKtAns8Xl5eYtbxSWFhYdi0aRPu3LmDuLg4eHl5yV3WEwQBa9aswe7du5GSkgIrKyu4u7ujadOmsm2Sk5Ph6ekpm/rHwcEBCxYsgL6+vujtISKqyAoMHhcXFzHr+KT09HRIpVI4Oztj3rx5Cut9fX3h5+eHpUuXonHjxli7di3GjBmDI0eOyEbhzZo1C8+fP4evry8kEgl+/vlnzJ07Fz4+PmI3h4ioQisTc7V16tQJP/zwA3r27AkNDfmSBUHA1q1bMXHiRPTo0QNSqRTLli1DWloaDh48CAB4/PgxQkJCsHDhQtja2sLGxga//vorTp8+jSdPnpREk4iIKqwyETyfEhMTg/j4eHTo0EG2TFtbG61atUJERAQAICIiAjo6OrC1tZVt06JFC+jo6Mi2ISIicZT54ImPjwcAhWeIDA0NZe8MSkhIgIGBgdxgCYlEAgMDA75XiIhIZGU+ePJ8PAJPmfWCIBS6HxERqVaZDx4jIyMA/9fzyZOYmCjrBdWsWROJiYlyQ78FQcDr169haGgoXrFERFT2g6dBgwYwMjJCaGiobFlGRgbCw/qJAiUAABgBSURBVMNhY2MDALCxsUF6errc/ZyIiAikp6fLtiEiInEU+UVwJSEtLQ3Pnj0DAOTm5uLff/9FZGQkqlWrhnr16sHV1RU+Pj4wMTGBsbEx1q1bBx0dHfTp0wcAYGpqCnt7e3h4eGDRokUQBAEeHh7o0qULTExMSrJpREQVjtKvRShJly9fhqurq8JyFxcXLF26VO4B0uTkZNkDpFKpVLZtUlKSwgOk7u7uhT5AypkLiIhUq0wET0li8BARqVaZv8dDRERlC4OHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETF4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRlYk3kFLRldR7gAC+C4iIPo09HiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETFKXNI5Upquh5O1UNUNjB4qNzg/HREZQMvtRERkagYPEREJCoGDxERiYrBQ0REomLwEBGRqBg8REQkKgYPERGJis/xEJVhfFiXyiIGD5EKlOTDq0RlDS+1ERGRqNjjIaIi4/RE9DnY4yEiIlGxx0NEZQoHVJR9DB4iolKuvIVthQueHTt2YNOmTYiPj0fTpk0xf/58tGzZsqTLIqJSjiMXVadC3eM5fPgwlixZgsmTJyM4OBg2NjaYMGEC/v3335IujYiowqhQwbN582a4uLhg0KBBMDU1xYIFC2BkZISdO3eWdGlERBVGhQmezMxM3LlzBx06dJBb3qFDB0RERJRQVUREFU+Fucfz+vVr5OTkoGbNmnLLDQ0NERoaWuj+ZmZmxTpvlWLtRURU8syCfvis/e/fv5/v8grT48kjkUiUWkZEROpRYXo8NWrUQKVKlRAfHy+3PDExUaEX9KGCEpuIiIqnwvR4tLS08NVXXylcVgsNDYWNjU0JVUVEVPFUmB4PAIwZMwZz586FpaUlbG1tsXPnTsTFxWHIkCElXRoRUYVRoYKnV69eeP36NdatW4e4uDhIpVJs2LAB9evXL+nSiIgqDIkgCEJJF0FERBVHhbnHQ0REpQODh4iIRMXg+Uw7duyAg4MDLCws0K9fP4SHhxe47eXLl2FmZqbw5/HjxyJW/PmK0mbg/awRf/zxBxwcHNC8eXN07twZW7duFala1ShKm93c3PL9OVtbW4tY8ecp6s/4wIEDcHJygpWVFTp06IDZs2crPLpQ2hW1zTt27MDXX38NS0tL9OjRA8HBwSJVqhphYWGYPHky7O3tYWZmhr179xa6z/379zFixAhYWlrC3t4ea9asQbHu1ghUbIcOHRKaNWsm7N69W3j06JGwcOFCwdraWoiNjc13+0uXLglSqVR4+PChEBcXJ/uTnZ0tcuXFV9Q2C4IgTJs2Tejfv79w/vx5ITo6Wrh+/bpw6dIlEav+PEVtc0pKitzPNy4uTujatavg5uYmcuXFU9T2hoeHC+bm5sLmzZuFZ8+eCREREYKzs7Pg6uoqcuXFV9Q279ixQ7C2thYOHDggPHv2TDh48KBgbW0tnDx5UuTKi+/MmTPCihUrhL///luwtLQUAgMDP7l9amqq0L59e+G7774T7t+/Lxw5ckSwtrYWNm3aVORzM3g+w4ABA4SffvpJblm3bt2E5cuX57t9XvAkJiaKUZ5aFLXNISEhgq2tbYVq88fCw8MFqVQqXL16VR3lqVxR27tx40ahc+fOcssCAgIEa2trtdWoakVt8+DBg4XFixfLLfPy8hKGDBmithrVydrautDg2bFjh2BjYyO8fftWtmzt2rWCnZ2dkJubW6Tz8VJbMX3OpKMDBgyAnZ0dRo0ahUuXLqmzTJUqTptPnDgBCwsL+Pv7o2PHjujevTs8PT2RlpYmRsmfTRWTy+7ZswdNmzaFra2tOkpUqeK019bWFvHx8Th16hQEQcCrV69w+PBhdOzYUYySP1tx2pyZmYkqVeRnYqxSpQpu3bqFrKwstdVakq5fv46WLVtCW1tbtszOzg5xcXGIiYkp0rEYPMX0qUlHC7q2bWRkhF9++QWrV6+Gt7c3GjdujNGjRyMsLEyMkj9bcdocHR2Nq1ev4t69e/D29saCBQsQEhKCH3/8UYySP1tx2vyh1NRUHDlyBAMHDlRXiSpVnPba2NhgxYoVmD17Npo3b4527dpBEAQsW7ZMjJI/W3HabGdnh8DAQNy8eROCIODWrVsICAhAVlYWXr9+LUbZoktISIChoaHcsrzvLCEhoUjHqlAPkKpDUSYdNTExgYmJieyzjY0NYmNjsWnTJrRq1UptNapaUdosCAIkEglWrFgBPT09AMCCBQswbtw4JCQkfHKevNKkuJPL7t+/Hzk5OXByclJHWWpTlPY+evQInp6emDp1Kuzs7BAfH4/ffvsN7u7u+O2339RdqsoUpc1Tp05FfHw8hg4dCkEQYGhoCGdnZ2zcuBGVKlVSd6kl5uPvQ/j/AwuKOtEyezzFVNxJRz9mZWWFf/75R9XlqUVx2mxkZITatWvLQgcATE1NAaBMvPn1c3/Of/31F7p3747q1aurq0SVKk57169fD0tLS4wfPx7m5uawt7eHh4cH9u3bh+fPn4tR9mcpTpu1tbXh5eWF69ev49SpUzhz5gzq168PXV1d1KhRQ4yyRVezZk2Fnk1iYiIAKPSECsPgKSZVTToaGRkJIyMjVZenFsVps62tLeLi4uTu6URFRQFAmZiq6HN+zjdv3sS9e/cwaNAgdZaoUsVp77t37xT+lZ/3WSgDE6N8zs9YU1MTderUQaVKlXD48GF06dIFGhrl89eqtbU1wsPDkZGRIVsWGhqKWrVqoUGDBkU6VqVffvnlFxXXV2FUrVoV3t7eMDIygra2Nv7880+Eh4djyZIl0NfXx9y5c3H8+HF069YNAODv74+EhARoaGggPj4e/v7+2LNnD9zc3GS9gNKuqG1u3LgxAgMDcfv2bTRp0gRRUVFYuHAhWrVqhf79+5dwa5RT1DbnWb16NdLT0zFv3rwy9c6norb33bt32LRpE2rUqIHq1avj0aNHWLx4MYyMjDB+/PgSbo1yitrmp0+f4uzZs9DW1sazZ8+wdOlS3LhxAytXroS+vn4Jt0Y5aWlpePz4MRISErBnzx5IpVLo6ekhKysLenp6WLFiBdavXw8XFxcAgLGxMXbv3o3IyEiYmJjg6tWrWLZsGSZNmlTkgTO8x/MZCpt09OPLDFlZWVi2bBlevnwJbW1tNGnSBBs2bECnTp1KovxiKWqbdXV1sXnzZnh6emLAgAHQ19eHo6MjZs2aVRLlF0tR2wwAb968weHDhzF16tQyFTpA0dvbr18/pKWlYceOHVi2bBn09PTQpk0bzJkzpyTKL5aitjk3Nxf+/v54+vQpKleujDZt2mDnzp1F/pd/Sbp9+zZcXV1ln729veHt7Q0XFxcsXboU8fHxiI6Olq3X09ODn58fFi5ciP79+6NatWoYO3YsxowZU+Rzc5JQIiISVfm8GElERKUWg4eIiETF4CEiIlExeIiISFQMHiIiEhWDh4iIRMXgISIiUTF4iIhIVAweIiISFYOHiIhExeAhIiJRMXiIiEhUDB4iIhIVg4eIiETF4CEiIlExeIiISFQMHiIAbm5uMDMzK+kyiCoEBg+VapGRkfD29kZMTExJl1JqeHt748SJE8XaNycnB3Z2djAzM8Off/5Z4HYODg5wcHAocH1eUH/4c/H29oaZmZnsj7m5OVq3bo1Ro0bh5MmTBR4rNzcXQUFBcHV1RZs2bdC8eXN07twZc+bMQWRk5Cfbc/HiRcyYMQOdOnVC8+bNYWNjg379+mHlypV48eLFJ/elksPgoVItMjISa9asQWxsrFrPs2jRIty8eVOt51CVNWvWFDt4zp07h/j4eDRq1Ah79+6FOt58/9133+G3337DkiVLMGzYMDx48ABTp07F/v37FbZNT0/H+PHj4ebmhnfv3mHChAnw8PBA7969cf78efTv3x//+9//FPbLzc3Fzz//jNGjR+PWrVvo06cPPDw8MHv2bJibm2PHjh1wcXFRedtINSqXdAFEqpSTk4PMzEx88cUXRdpPU1NTTRWVLgEBAWjYsCHc3NwwdepUXL58GW3btlXpOTp27AgLCwvZ5549e8LJyQm+vr745ptv5Lb18PDAhQsXMHnyZHz//fdy68aNG4fRo0dj4cKFMDY2Rvv27WXrvL29sWfPHvTp0wdeXl7Q0tKS29fNzQ1r1qxRabtIddjjoVLL29sbP/74IwDA1dVVdgnHzc0NALB3716YmZkhNDQUa9euhaOjIywtLfH3338DAM6fP4+ZM2eia9eusLS0RMuWLTF27FhcuXJF4Vz53ePJW5aamgoPDw+0a9cOFhYWGDJkCG7cuKFUG5KSkrBkyRI4OjrCwsICbdq0Qb9+/bBx40aFbQ8fPoyhQ4fCxsYGVlZWGDhwII4cOSJbHxMTI6sxKChI7rKWMhITE3H27Fk4OzujU6dOMDQ0REBAgFL7fg5zc3PUqFEDUVFRcsvv3buH/fv3w8rKCjNnzlTYz8DAACtWrAAALF++XLY8MTERmzZtQv369bFkyRKF0AEAfX19zJ8/X7UNIZVhj4dKrW7duiE+Ph67d+/G5MmTYWJiAgBo1KiR3HbLli1DdnY2Bg0aBF1dXTRu3BjA+1/OycnJcHZ2Rp06dfDy5Uvs2bMHo0ePxtatW9GyZUul6hg3bhwMDAzw7bffIikpCZs3b8bEiRNx8uRJVK1a9ZP7zpgxA+Hh4Rg8eDDMzc3x9u1bPHnyBFeuXMH48eNl261cuRI+Pj6wt7fHjBkzoKGhgePHj2PGjBlwd3fH8OHDYWBggN9++w1z585Fy5YtMWjQoKJ8nQgKCkJ2djacnJxQuXJl9O3bF7t27UJqair09PSKdKyiSE5ORkpKCgwMDOSWHzt2DAAwcOBASCSSfPdt2rQprK2tERERgdjYWNSvXx9nzpxBRkYGnJycUKVKFbXVTerD4KFSy9zcHNbW1ti9ezfat2+PNm3a5Lvdu3fvEBwcrHB5bdGiRdDR0ZFbNmTIEPTu3Rvr169XOniaNWuGX375RfbZ1NQUM2fOxMGDBzFkyJAC90tNTcWlS5cwdOhQuLu7F7jdnTt34OPjg0mTJuGHH36QLXd1dcXUqVOxYsUKODk5oWrVqnBycsLcuXPRsGFDODk5KVV/nr1796JVq1Zo2LAhAMDFxQX+/v44cOAAhg0bVqRjfcqbN2/w6tUr5OTk4NmzZ1i1ahVycnIULrM9fPgQwPvv91O++uorRERE4MGDB6hfv75svy+//FJlNZO4GDxU5g0dOjTfezofhk5aWhoyMzOhoaEBKysrpS+VAcDo0aPlPufdE/nnn38+uV+VKlWgpaWFmzdvIiYmBg0aNMh3uwMHDkAikcDZ2RmvXr2SW+fg4ICTJ0/i+vXrsLOzU7rmj127dg2PHz/GuHHjZMvMzc3x5ZdfIjAwUKXB8/H3paWlhQkTJihcTnvz5g0AFNrbyutVpqamyu1XWG+TSi8GD5V5eZfWPvbs2TOsXLkS58+fR0pKity6gi7t5Cevh5CnRo0aAN7fv/kULS0tzJ8/H4sXL0bXrl3RpEkTtG3bFo6OjmjXrp1su8ePH0MQBHz99dcFHishIUHpevMTEBAATU1NNGvWTC4w7ezs4Ovri3v37sHc3LzIx83ve3R3d0fjxo3x9u1bXL58Gdu2bUNKSgoqV5b/dfNxoBQkLS0NwP8FVN5+ecup7GHwUJmnra2tsCwtLQ3Dhw/H27dvMWrUKEilUujq6kJDQwPr16/HpUuXlD5+pUqV8l2uzFDkoUOHomvXrjh79iyuXLmCo0ePYvv27ejVqxdWrlwpO45EIoGvr2+B52rSpInS9X4sLS0Nf//9N7KysuDs7JzvNoGBgfjpp59kn6tUqfLJYE1PTweQ/3dvaWkpG9XWtWtX1KxZEytWrMCXX36JoUOHyrZr2rQpjh07hrt37+Krr74q8Fx37twBAEilUtl+AHD37l1069atwP2o9GLwUKlWlJ7Jhy5evIi4uDgsWbIE/fv3l1u3atUqVZSmtFq1amHgwIEYOHAgcnJyMHfuXBw8eBBjxoyBpaUljI2NERISgnr16sHU1FTl5//777+Rnp6OH374Af/5z38U1m/btg379+/HnDlzZCPEGjRogCdPnuDVq1cKgwIA4MmTJ9DV1ZX1/j5lzJgxCAgIwKpVq9C3b19Zj6V79+5Yu3YtAgICMGDAgHx/1o8ePUJERAS++uor1K9fHwDQuXNnVKlSBfv27cOUKVPyHdVGpRuHU1OplnefJjk5uUj75fUcPu6VnD9/vkj3dz7H27dv8fbtW4W68oY/57Up76b7f//7X+Tk5CgcJzExUe6zjo5OoZf5PhQQEIDq1atj3Lhx6Nmzp8KfAQMGICkpSW52AUdHRwCAv7+/wvFCQkLw8OFDdOnSBRoahf8K0dTUxKRJk5CUlIStW7fKlpubm6NPnz64fv06vL29FfZLSkrCnDlzAACzZs2SLTc0NMS4ceMQGxuLn376CZmZmQr7vnnzBkuWLCm0NioZ7PFQqWZhYQENDQ34+PggOTkZOjo6aNCgAaysrD65X4sWLWBkZIRly5YhNjYWderUQWRkJPbt2wepVIoHDx6ovfaoqCiMGDEC3bp1Q9OmTaGvr48nT55g586daNCggWxUnaWlJaZPnw5vb284OzujR48eqF27NuLi4nDnzh2cO3cOt2/flh3X2toaFy9exIYNG1CvXj1IJBL07t073xoeP36MiIgI9OvXT+EeSx4HBwdoamoiICBAdp+pX79+OHDgANavX4+7d++iXbt20NbWxt27dxEcHAwjIyO5MCiMk5MT1q5dC39/f4wcOVJ2v+bXX39FQkIC1q5di9DQUHTr1g3VqlVDVFQU9u7di9evX8Pd3R0dOnSQO9706dMRHx+PPXv24OrVq+jduzcaNWqErKws3Lt3D0eOHIGmpiaf5SmlGDxUqtWrVw9LliyBr68vfv31V2RlZcHFxaXQ4NHX18fGjRvx+++/Y/v27cjOzkbz5s3h6+uLgIAAUYKnTp066N+/Py5fvowTJ04gMzMTtWvXxsCBAzFhwgS5kXjTpk1D8+bNsW3bNmzduhXp6ekwNDRE06ZNFX55enh4YOHChfDx8ZHdYC8oePIeEP3UvZBq1aqhTZs2CA0NxfPnz1G3bl1oamrCz88P/v7+OHToELy9vZGdnY06depg8ODBmDx5MmrVqqX0d1G5cmVMnDgRHh4e2LJlC6ZNmwbg/UABPz8/BAcHIzg4GOvXr5e1vUOHDhg7dmy+w6Y1NDTg6emJXr16YdeuXdi3bx9evXoFTU1NNG7cGMOGDVPpSD1SLYmgjsmaiIiICsB7PEREJCoGDxERiYrBQ0REomLwEBGRqBg8REQkKgYPERGJisFDRESi4gOkRKVUdHQ0NmzYgLCwMDx//hxaWlowMjKChYUFXFxcEBwcjKCgIKWONW3aNEyfPl3NFRMph8FDVArdunULI0eOROXKleHs7IwmTZrg3bt3iIqKwunTp6Grq4vBgwfLvV4BAObOnQsTExNMnjxZbrmyr8cmEgNnLiAqhSZPnozTp08jODhYYcqY3NxcxMfHo3bt2gr7mZmZoXXr1ti2bZtYpRIVGe/xEJVCUVFRqF69eoHzlOUXOkRlBYOHqBRq1KgRkpKScOzYsZIuhUjleI+HqBSaMmUKQkNDMX36dBgbG8PW1hYWFhZo06aNWl4WRyQm3uMhKqXu37+PzZs349y5c3Ivg2vRogWWLVuGhg0bKuzDezxUFjB4iMqA2NhYhIWFYc+ePQgPD4dUKkVgYKDCa58ZPFQW8B4PURlQv359ODs7Y/v27bC1tcWDBw9w8+bNki6LqFgYPERliEQikb19NS4uroSrISoeBg9RKXThwgVkZ2crLH/37h0uXLgAABxkQGUWR7URlUJeXl5ISkqCg4MDpFIptLW18eLFCxw4cABRUVFwdnbmbARUZjF4iEohNzc3nDx5ElevXsXRo0eRmpoKPT09SKVSTJgwAf369SvpEomKjaPaiIhIVLzHQ0REomLwEBGRqBg8REQkKgYPERGJisFDRESiYvAQEZGoGDxERCQqBg8REYmKwUNERKJi8BARkaj+H+08he1V8w/fAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(STtrain)\n",
    "plt.xlabel(\"train set AUROC\\nST\")\n",
    "plt.ylabel(\"brain areas\")\n",
    "plt.ylim([-2,410])\n",
    "plt.text(0.75, 400, \"mean = %.3f\" %np.mean(STtrain),size=20)\n",
    "plt.savefig(\"051221_1.pdf\", bbox_inches='tight')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAE2CAYAAABRM/SMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deVxO6f8/8NddSu6UpbIvCZUZlYosFWMZ2clOI9tYZoYx9mUkQzNhzGdmxIhIloylZM0Y60wYS2SNkInC6C7tadP5/dG383OruEv3afF6Ph4eH53tfl81n16uc13nOjJBEAQQERFJRKOsCyAiog8Lg4eIiCTF4CEiIkkxeIiISFIMHiIikhSDh4iIJMXgISIiSVUp6wKIiMq706dPw9fXF+Hh4cjNzUWLFi0wevRoODs7q3yNmJgYdO/evcj9ffr0wc8//1zovqCgIPj7+yMyMhIaGhr46KOPMGHCBHTt2rXQ46Ojo+Hj44NLly7h2bNn0NXVRePGjTF48GA4OztDW1tb5brVQcYHSImIirZjxw4sX74cNWvWRJ8+faClpYVjx47hv//+w4QJEzB//nyVrpMfPObm5ujRo0eB/S1btkSvXr0KbF+5ciV8fX1Rr149ODk5ITs7G8HBwUhMTISbmxs+++wzpeNv3LiBsWPHIiMjA46OjmjZsiVSU1Nx+vRpPH/+HA4ODti0aRNkMlnJviGlQSAiokJFR0cLrVu3Fuzs7ITo6Ghxe2JiotCjRw/B1NRUuHr1qsrXMjU1FebPn6/y51+5ckUwNTUVevToISQmJipdy87OTmjdurVSXYIgCJMmTRJMTU2Fffv2KW1PS0sT+vTpI5iamgqXLl1SuQZ14BjPO5iZmcHMzKysyyAqUkxMDMzMzLBgwQI8fvwYX3/9Ndq3bw9ra2tMmDAB9+7dAwC8ePECbm5ucHBwgIWFBYYMGYILFy4Ues2cnBz4+/tj+PDhsLGxgZWVFQYNGoQdO3YgNze3wPH79u3D9OnT0b17d1haWsLGxgYjR47EgQMHCr3+mDFjYGZmhpycHHh7e6Nnz55o3bo1unTpgh9//BFZWVml9w16D4GBgcjKyoKLiwsaNWokbq9RowamTJkCANi1a5faPj//2lOnTkWNGjXE7Y0aNcLo0aORlZWFffv2KZ0THR0NAOjWrZvSdrlcjo4dOwLI+2+hLDF4iCqJJ0+eYNiwYYiLi4OzszMcHBxw/vx5jBkzBlFRURg+fDhu3ryJ3r17o3fv3oiIiMCkSZPw9OlTpetkZ2djypQpWLZsGZKTk9GvXz8MHz4cubm5WL58eaG3lpYuXYonT56gXbt2GDt2LPr27YunT59i3rx5+OWXX4qsefbs2dixYwdsbW0xatQo6OjoYNOmTXB3dy/1709J5Aezo6NjgX2dO3dWOkZVsbGx2LVrF7y9vbFr1y7cvXu3VD+/ZcuWAIAzZ84obX/58iUuXLgAuVyONm3aFKvm0sbJBUSVxKVLl/DNN9/giy++ELetW7cOa9aswfDhw9GrVy8sXboUGhp5/97s1KkT5s+fDz8/PyxatEg8x9vbG2fPnsVnn32GRYsWQVNTEwDw6tUruLm5ITAwEE5OTkrjFIcPH0aTJk2U6snKysKkSZPg4+ODUaNGoW7dugVqjo6OxuHDh1GzZk0AwMyZMzFw4EDs378fs2bNgpGR0TvbfefOHZw4caIY3ylg7Nix0NfXf+dx//77LwDA2Ni4wL46depALpfjv//+w8uXL1GtWjWVPvvcuXM4d+6c0jY7OzusXLkSDRo0ELelp6fj+fPnkMvlqFOnToHrNG3aFAAQFRWltH3GjBm4evUqFi5ciKNHj6JFixZITU3FmTNnkJOTg19++aXQn4WUGDxElUTDhg0xefJkpW3Ozs5Ys2YNsrKyMG/ePDF0AKB///749ttvcefOHXFbbm4uduzYASMjIyxcuFAMHQDQ1NTEggULsG/fPhw6dEgpeN4MHQDQ1taGi4sLLly4gH/++QeDBg0qcMycOXPE0AHybgf1798f69atw61bt4qctfW6O3fuYO3ate887nXOzs4qBU9qaioAQE9Pr9D91atXR3p6OlJSUt4ZPNWqVcOXX36JHj16oHHjxgCAiIgIeHl54eLFixg3bhz2798PuVwOAEhJSXnrZ+dvT05OVtrevHlzBAYGYtasWTh9+jROnz4NANDS0oKrq2uZ93YABg9RpdGqVSuloAAg/kvZ2NgY1atXV9qnqakJAwMDPH/+XNz277//IjExEcbGxli/fn2hn6Ojo4OHDx8qbXv69Cl8fHzwzz//4NmzZ8jIyFDa//pnvK5169YFttWvXx8AkJSUVOg5bxo8eDAGDx6s0rGlTfi/ScGqzBAzMDDAjBkzlLa1a9cOvr6+GD16NK5fv469e/di7Nixxarhzc8ODw/HV199hdq1a8Pf3x+tWrVCSkoKDh48iF9++QUnT55EQEBAkYEmBQYPUSVR2C+SKlWqFLkvf39OTo74dWJiIoC82zdv60WkpaWJf4+OjsbQoUORnJyMtm3bwsHBAdWrV4empiaePHmCoKCgIicLFNbryA/PwiYxSK169epISEhASkoKatWqVWB//vfhzVAvjipVqmDYsGG4fv06QkNDxeDJ/5nl93zeVFiPKCcnB9988w1evHiBPXv2iLcqdXV1MXnyZMTFxWHr1q3w8/PD9OnTS1zz+2LwEJEo/5fYp59+qvLtqy1btiAxMRGenp4Feh6HDx9GUFBQqdf5OnWO8TRr1gwJCQmIiooqEDyxsbFIT09HvXr1VB7fKUr+tdPT08VtcrkcdevWxfPnzxEbG1tgnOfRo0cAlMefHj58iEePHuHjjz8udHysffv22Lp1K27fvv1e9b6vChc83t7e+Pnnn+Hi4oIlS5YAyOvurl27Frt370ZycjKsrKywZMkScXYHkNdt9/DwwKlTpwDkTTV0c3NT6T++kjBecEQt141a0Vct1yUCABMTE+jr6+PatWvIzs6GlpbWO8/J/wXYs2fPAvsuXbpU6jW+SZ1jPB06dMDVq1cREhICa2trpX1///23eMz7un79OgCIYz+vf/6BAwcQEhKCIUOGvPPz83uWCQkJhX5O/jRqVX6u6lShplNfu3YNe/bsKfBcjY+PD3x9feHm5oaAgADUrl0b48ePFwcGgbxpm+Hh4fDx8cGmTZsQHh6OefPmSd0EonKtSpUq+Oyzz6BQKODh4VFgrAbI+5f+gwcPxK8bNmwIoGDIhISEICAgQL0FI2+MJyIiolh/Xn8m513X1tbWhr+/P2JiYsTtSUlJ2LBhAwBg5MiRSuekpKQgMjISsbGxStuvX79e6C3Hf/75B35+fgCAAQMGKO3Lv7a3t7fSmFdMTAx27twJbW1tpV6mqakp9PX18fTpU+zdu1fpWsnJyfD19QVQOmH5PipMjyclJQVz5szB999/j99++03cLggCtm3bhsmTJ8PJyQlA3hITHTt2xOHDhzFy5EhERkYiJCQEO3fuhI2NDQDgu+++g4uLCx4+fAgTE5MyaRNRefTll1/i7t272LVrF06fPo0OHTqgbt26iI+Px6NHj3D16lXMnDkTLVq0AACMHj0a+/btw4wZM9CzZ0/UrVsX9+/fR0hICHr37o3g4OAyblHJNW7cGPPmzYOHhweGDBlS6JI5b/aEjh8/joULF8LZ2RkrVqwQt69evRr379+HnZ0d6tWrByBvVlv+czgzZswQfz/ls7Gxwfjx47FlyxYMGDCg0CVzXg9RbW1tLFq0CAsXLsTixYtx5MgRtGrVCsnJyTh16hRevHiBNm3aYNiwYer6lqmkwgSPm5sbnJyc0LFjR6XgiYmJgUKhgL29vbhNR0cH7dq1Q1hYGEaOHImwsDDI5XKlH6qtrS3kcjnCwsIYPESv0dLSwm+//YYDBw4gKCgIZ86cQXp6OmrVqoVGjRphxowZ6N+/v3i8ubk5tm3bhl9++QV///03cnJyYG5ujrVr10JPT69CBw+Qt8pCw4YN4evri/3790MQBDRv3hzffPNNsRYJHTBgAE6cOIFbt24hJCQE2dnZMDQ0RO/evfHZZ5+hbdu2hZ63YMECmJmZYceOHdizZw9kMhk+/vhjTJw4sdDp5s7OzmjUqBG2bt2Ka9eu4fLly9DW1kazZs0wfvx4jB07louEqmLPnj3YtWsXdu3aBW1tbYwZMwYtW7bEkiVLcPXqVYwaNQqnT59Wevhq4cKFiI2NxebNm+Ht7Y29e/fi5MmTStft3r07hg8fLi59UZj823oRERHFqpljPEREhSv3PZ6HDx/if//7H/z9/d+a0u+aR1/YfkEQynaFViKiD1C5D55r164hISFBqWv/6tUrXL58Gbt27cLhw4cBAAqFQnzwDADi4+NhaGgIADA0NER8fLxS0AiCgISEBBgYGEjYGiIiKvfB06NHjwJPNy9cuBDGxsaYMmUKmjVrBiMjI5w/fx6WlpYAgMzMTISGhoqz1qytrZGeno6wsDBxnCcsLAzp6ekFBgaJiEi9yn3w6OvrF5hvL5fLUaNGDZiamgIAXF1d4e3tDRMTE3GpD7lcjn79+gHIW7vI0dER7u7uWL58OQRBgLu7O7p27cqJBUREEiv3waOKSZMmITMzE8uWLUNSUhKsrKzg6+urtIzF6tWr4eHhgQkTJgDIe4A0/wFUIiKSToWY1VaWOKuNiKh0VaiVC4iIqOJj8BARkaQYPEREJCkGDxERSYrBQ0REkmLwEBGRpBg8REQkKQYPERFJisFDRESSYvAQEZGkGDxERCQpBg8REUmKwUNERJJi8BARkaQYPEREJCkGDxERSYrBQ0REkmLwEBGRpBg8REQkKQYPERFJisFDRESSYvAQEZGkGDxERCQpBg8REUmKwUNERJJi8BARkaQYPEREJCkGDxERSYrBQ0REkmLwEBGRpBg8REQkKQYPERFJisFDRESSYvAQEZGkGDxERCQpBg8REUmKwUNERJJi8BARkaQYPEREJCkGDxERSarcB4+/vz/69+8PGxsb2NjYYMSIEThz5oy4XxAEeHl5wcHBAZaWlhgzZgzu37+vdI2kpCTMnTsXtra2sLW1xdy5c5GcnCxxS4iICKgAwVO3bl3MmTMHQUFBCAwMRIcOHfDVV1/h7t27AAAfHx/4+vrCzc0NAQEBqF27NsaPH4/U1FTxGrNnz0Z4eDh8fHywadMmhIeHY968eWXVJCKiD5rKwZOQkIDIyEilbdHR0Vi+fDlmz56NkJCQUi8OAHr06IEuXbqgadOmaNasGWbOnAldXV1cu3YNgiBg27ZtmDx5MpycnGBqaoqVK1ciLS0Nhw8fBgBERkYiJCQEy5Ytg42NDaytrfHdd9/h9OnTePjwoVpqJiKioqkcPN9//z3mz58vfp2WlgYXFxf4+/vjyJEjmDJlCi5fvqyWIvO9evUKR44cQXp6OqytrRETEwOFQgF7e3vxGB0dHbRr1w5hYWEAgLCwMMjlctjY2IjH2NraQi6Xi8cQEZF0VA6ea9euoXPnzuLXwcHBiI2NxcaNGxESEoLmzZtj06ZNaikyIiIC1tbWsLCwgLu7O9auXQszMzMoFAoAgKGhodLxBgYGiIuLAwDExcWhdu3akMlk4n6ZTIbatWuLxxARkXSqqHpgfHw86tevL34dEhKC1q1bi2Hk7OyMLVu2lH6FAJo1a4b9+/cjOTkZf/75J+bPn4/t27eL+18PlcIUtl8QhHeeR0REpU/lHk+VKlWQmZkpfn3p0iW0a9dO/FpPTw+JiYmlW93/0dbWRtOmTWFhYYHZs2ejVatW8PPzg5GREQCIPZ988fHxYi/I0NAQ8fHxEARB3C8IAhISEmBgYKCWeomIqGgqB4+xsTGOHTsGQRBw8uRJJCUloWPHjuL+//77DzVq1FBLkW/Kzc1FVlYWGjVqBCMjI5w/f17cl5mZidDQUFhbWwMArK2tkZ6erjSeExYWJo4TERGRtFS+1ebi4oIFCxagXbt2yMjIQOPGjZWC5/LlyzAzMyv1AlevXo1PPvkE9erVE2erXbp0CRs2bIBMJoOrqyu8vb1hYmICY2NjrF+/HnK5HP369QMANG/eHI6OjnB3d8fy5cshCALc3d3RtWtXmJiYlHq9RET0dioHz6BBgyCTyXDixAlUr14dU6dOhZaWFoC8qdapqakYNWpUqRcYFxeHuXPnQqFQQE9PD2ZmZvDx8YGjoyMAYNKkScjMzMSyZcuQlJQEKysr+Pr6onr16uI1Vq9eDQ8PD0yYMAEA0K1bNyxZsqTUayUioneTCa8PflAB+b24iIiIYp1nvOCIOspB1Iq+arkuEZFUyv3KBUREVLmofKsNAHJycnDixAlcv34dycnJyM3NVdovk8nwww8/lGqBRERUuagcPImJiXB1dcX9+/fFZ2Dy79Ll/53BQ0RE76LyrbZffvkFDx8+hIeHB44fPw5BELB582YEBwejb9++sLCwwMWLF9VZKxERVQIqB89ff/2FQYMGYciQIeKMMQ0NDZiYmGD16tXQ0dHBTz/9pLZCiYioclA5eBQKBSwsLADkrWIAAFlZWeL+7t2749SpU6VcHhERVTYqB0/NmjXx8uVLAICuri6qVKmCZ8+eifu1tLT4cjUiInqnYi2Z8+DBg7yTNDTw0UcfISgoCFlZWXj58iX279+Pxo0bq61QIiKqHFQOHnt7exw7dky8vTZu3Dhcv34ddnZ26NSpE27duoWxY8eqrVAiIqocVJ5OPXXqVEycOBHa2toAgD59+qBKlSo4ePAgNDQ00KtXL/Tp00dthRIRUeWgcvDIZDIxdPL17NkTPXv2LPWiiIio8irRkjlZWVl4/vy50qw2IiIiVRQreG7fvg1XV1fY2Njgk08+wZUrVwDkvXht7NixSu/FISIiKozKwXPnzh24uLggOjoaAwcOVNpnYGCAzMxMBAUFlXqBRERUuagcPL/++ivq1KmDw4cPY/bs2XjzbQodOnTAjRs3Sr1AIiKqXFQOnitXrmDYsGHQ1dWFTCYrsL9BgwaIjY0t1eKIiKjyUTl4MjMzoaenV+T+1NTUUimIiIgqN5WDp0mTJrh9+3aR+y9cuIAWLVqUSlFERFR5qRw8/fr1w4EDB5RmruXfcvP19UVISEiBSQdERERvUvkB0gkTJuDcuXOYOHEiTExMIJPJ4OnpiRcvXiAuLg6dOnXC6NGj1VkrERFVAir3eLS1tbFlyxbMnz8fVatWRdWqVREVFYVatWph7ty52LBhAzQ0SvQ8KhERfUBU7vEAee/hGTduHMaNG6emcoiIqLJTqYuSlpaGHj16wM/PT83lEBFRZadS8Ojq6iIxMRG6urrqroeIiCo5lQdlrKyscPPmTXXWQkREHwCVg2fOnDn4448/EBgYWGC5HCIiIlWpPLnA09MT+vr6WLx4MX788Uc0adIEOjo6SsfIZDJs3bq11IskIqLKQ+XgiYmJAQDUr18fABAXF6eeioiIqFJTOXhOnTqlzjqIiOgDwSc+iYhIUgweIiKSVLFWLnj8+DH8/Pxw/fp1JCcnIzc3V2m/TCbDiRMnSrVAIiKqXFTu8URERMDZ2Rl79+5FdnY2oqOjIZfLkZmZiSdPnkBTU1OceEBERFQUlYNnzZo10NLSwoEDB8SlcxYtWoSzZ89i2bJlSE5Ohru7u7rqJCKiSqJYr74eMWKE+EqE1w0fPhydO3fG6tWrS71AIiKqXFQOnrS0NDRu3BgAoKWlBQBIT08X99vY2ODq1aulXB4REVU2KgePoaGh+NBo9erVUa1aNURFRYn7k5OT8erVq1IvkIiIKheVZ7WZm5srLRJqZ2eHbdu2wdLSErm5udixYwfMzc3VUiQREVUeKvd4+vfvj8TERGRkZAAAZsyYgZSUFLi6umLcuHFISUnBzJkz1VYoERFVDjLhPZaafvbsGY4fPw5NTU107txZHAOqTMzMzADkTScvDuMFR9RRDqJW9FXLdYmIpFKsB0jfVL9+fbi6upZWLURE9AEo90vmbNiwAUOGDIGNjQ06dOiAqVOn4t69e0rHCIIALy8vODg4wNLSEmPGjMH9+/eVjklKSsLcuXNha2sLW1tbzJ07F8nJyVI2hYiIUAGC59KlSxg9ejR27dqFrVu3QlNTE+PHj0diYqJ4jI+PD3x9feHm5oaAgADUrl0b48ePR2pqqnjM7NmzER4eDh8fH2zatAnh4eGYN29eWTSJiOiD9l5jPGUhLS0Nbdu2xbp169CtWzcIggBHR0e4uLjgiy++AABkZGSgY8eOmD9/PkaOHInIyEj06dMHO3fuhK2tLQAgNDQULi4uOHr0KExMTIr8PI7xEBGVrnLf43lTWloacnNzoa+vDyDvBXUKhQL29vbiMTo6OmjXrh3CwsIAAGFhYZDL5bCxsRGPsbW1hVwuF48hIiJpVLjg+f7779GqVStYW1sDABQKBYC8B1xfZ2BgID7wGhcXh9q1ayst9SOTyVC7dm2+SZWISGLvNatNap6enrhy5Qp+//13aGpqKu17c/24NxW2XxCEd55HRESlq0TB8/LlSyQmJqKw4aEGDRq8d1GF+eGHHxAcHIytW7cqPS9kZGQEIK/n8/prGeLj48VekKGhIeLj45WCRhAEJCQkwMDAQC31EhFR4VQOntzcXGzatAnbt29/6+2pO3fulEphr/Pw8EBwcDC2b9+O5s2bK+1r1KgRjIyMcP78eVhaWgIAMjMzERoaKs5as7a2Rnp6OsLCwsRxnrCwMKSnp4u37IiISBoqB8/q1avh6+uLli1bwsnJCTVr1lRnXaLvvvsOBw4cwLp166Cvry+O6cjlcujq6kImk8HV1RXe3t4wMTGBsbEx1q9fD7lcjn79+gEAmjdvDkdHR7i7u2P58uUQBAHu7u7o2rXrW2e0ERFR6VN5OrWDgwNatWoFHx8fddekJH8685umTZuG6dOnA8i7bbZ27Vrs3r0bSUlJsLKywpIlS2Bqaioen5iYCA8PD5w6dQoA0K1bNyxZskScHfeuz+d0aiKi0qFyjyc5ORndu3dXZy2FUuUXvkwmw/Tp08UgKkzNmjX5ojoionJA5enUpqam4m0uIiKiklI5eKZNm4Zdu3bh2bNn6qyHiIgqOZVvtd26dQsNGjRAnz598Omnn6JRo0bQ0FDOLZlMhq+++qrUiyQiospD5eBZu3at+PeDBw8WegyDh4iI3kXl4Dl58qQ66yAiog+EysHTsGFDddZBREQfiAq3SCgREVVsRfZ41q5dC5lMhi+++AIaGhpKYzxF4RgPERG9S5ErF5ibm0Mmk+H69evQ1taGubn5uy8mk6llrbayxJULiIhKV5E9nvzJBNra2kpfExERvY8ig+fNyQScXEBERKWBkwuIiEhSxX4R3M2bN3Hjxg0kJSUhNzdXaR8nFxAR0buoHDwZGRmYNm0azp07J77JM39eQv7fGTxERPQuKt9qW7duHc6dO4epU6di27ZtEAQBK1asgI+PD9q2bQsLCwscOaKemVxERFR5qBw8x44dQ69evTBjxgy0bNkSAFC3bl04Ojpiy5YtyM7ORlBQkNoKJSKiykHl4Hn27BnatWsHANDU1AQAZGdnAwCqVKmCvn37ssdDRETvpHLw6Orq4tWrV+LfNTQ0EBsbK+7X09NDXFxc6VdIRESVisrB06RJE0RFRQHI6/G0aNECx44dAwAIgoDjx4+jfv36aimSiIgqD5WDp2PHjjh27JjY6xkxYgRCQkLQo0cP9OzZE+fPn8eQIUPUVigREVUOKk+nnjx5MgYOHChOoXZxcUFWVhYOHjwIDQ0NzJw5E5MmTVJboUREVDkUuUgo5eEioUREpUulW21paWno0aMH/Pz81FwOERFVdioFj66uLhITE6Grq6vueoiIqJJTeXKBlZUVbt68qc5aiIjoA6By8MyZMwd//PEHAgMDwWEhIiIqKZVntXl6ekJfXx+LFy/Gjz/+iCZNmkBHR0fpGJlMhq1bt5Z6kUREVHmoHDwxMTEAID4kylUKiIioJFQOnlOnTqmzDiIi+kDwDaRERCSpYr+BNCsrCxcvXkR0dDQAoHHjxrCzs0PVqlVLvTgiIqp8ihU8+/fvh6enJ5KTk5XePqqvr4/58+dj8ODBaimSiIgqD5WDJzg4GAsWLECDBg0wceJENG/eHIIgIDIyErt27cK3334LHR0d9OnTR531EhFRBafyWm0DBgxATk4O9uzZg+rVqyvtS0lJwbBhw6ClpYVDhw6ppdCywrXaiIhKl8qTC/79918MHjy4QOgAeS+BGzx4MB49elSqxRERUeWjcvAYGRm9dcUCDQ0NGBoalkpRRERUeakcPM7OzggKCkJaWlqBfampqQgMDOTkAiIieqciJxdcvnxZ6eu2bdvi9OnT6N+/P0aPHg0TExPIZDI8ePAAv//+O2rVqgVbW1u1F0xERBVbkZMLzM3NIZPJlLa9fmj+vje33blzRx11lhlOLiAiKl1F9ng8PT2lrIOIiD4QRQaPs7OzlHW81eXLl7F582bcvn0bsbGx8PT0VBpPEgQBa9euxe7du5GcnAwrKyssWbIELVu2FI9JSkqCh4eHuOZct27d4ObmBn19fcnbQ0T0IasQa7Wlp6fD1NRUfEj1TT4+PvD19YWbmxsCAgJQu3ZtjB8/HqmpqeIxs2fPRnh4OHx8fLBp0yaEh4dj3rx5UjaDiIhQQYKnS5cumDVrFnr16gUNDeWSBUHAtm3bMHnyZDg5OcHU1BQrV65EWloaDh8+DACIjIxESEgIli1bBhsbG1hbW+O7777D6dOn8fDhw7JoEhHRB6tCBM/bxMTEQKFQwN7eXtymo6ODdu3aISwsDAAQFhYGuVwOGxsb8RhbW1vI5XLxGCIikkaFDx6FQgEABR5eNTAwEF9WFxcXh9q1ayvN0pPJZKhduzZfaEdEJLEKHzz53pz6rcp+QRDeeR4REZWuCh88RkZGAP5/zydffHy82AsyNDREfHy80jNHgiAgISEBBgYG0hVLREQVP3gaNWoEIyMjnD9/XtyWmZmJ0NBQWFtbAwCsra2Rnp6uNJ4TFhaG9PR08RgiIpJGsd9AWhbS0tLw+PFjAEBubpQLdr4AABhzSURBVC6ePn2KO3fuoEaNGmjQoAFcXV3h7e0NExMTGBsbY/369ZDL5ejXrx8AoHnz5nB0dIS7uzuWL18OQRDg7u6Orl27wsTEpCybRkT0wVH5fTxl6eLFi3B1dS2w3dnZGStWrFB6gDQpKUl8gNTU1FQ8NjExscADpEuWLHnnA6RcMoeIqHRViOApSwweIqLSVeHHeIiIqGJh8BARkaQYPEREJCkGDxERSYrBQ0REkmLwEBGRpBg8REQkKQYPERFJisFDRESSYvAQEZGkGDxERCQpBg8REUmKwUNERJJi8BARkaQYPEREJCkGDxERSYrBQ0REkmLwEBGRpBg8REQkKQYPERFJisFDRESSYvAQEZGkGDxERCQpBg8REUmKwUNERJJi8BARkaQYPEREJCkGDxERSapKWRdAxWO84EipXzNqRd9SvyYRUVHY4yEiIkkxeIiISFIMHiIikhSDh4iIJMXgISIiSTF4iIhIUgweIiKSFIOHiIgkxeAhIiJJMXiIiEhSXDKH1LIMD8CleIiocOzxEBGRpD644PH390e3bt1gYWGBwYMHIzQ0tKxLIiL6oHxQwRMcHIwffvgBU6dOxf79+2FtbY1Jkybh6dOnZV0aEdEHQyYIglDWRUhl2LBhMDMzg4eHh7itZ8+ecHJywuzZsws9x8zMDAAQERFRrM9S17gJceyIqKL7YCYXZGVl4fbt25gwYYLSdnt7e4SFhZVRVVQSnAxR8b4HfI8Uve6DCZ6EhAS8evUKhoaGStsNDAxw/vz5d56f3/NRVdViHU3lgVnQrLIuQWXq+u9LXd8DddRbkX5eH6qi7hR9UGM8ACCTyVTaRkRE6vHB9Hhq1aoFTU1NKBQKpe3x8fEFekGvK+7YDhERvd0H0+PR1tbGxx9/XOC22vnz52FtbV1GVRERfXg+mB4PAIwfPx7z5s2DpaUlbGxs8PvvvyM2NhYjR44s69KIiD4YH1Tw9OnTBwkJCVi/fj1iY2NhamqKjRs3omHDhmVdGhHRB+ODeo6HiIjK3gczxkNEROUDg4eIiCTF4JFYcRYpvXjxIszMzAr8iYyMlLDi4inuIqxZWVn49ddf0a1bN7Ru3RqffPIJtm3bJlG1xVec9i1YsKDQn1+bNm0krLh4ivvzO3ToEAYOHAgrKyvY29tjzpw5BR5ZKC+K2zZ/f3/07t0blpaWcHJywv79+yWqtPguX76MqVOnwtHREWZmZti3b987z4mIiMBnn30GS0tLODo6Yu3atZBs5EUgyRw5ckT46KOPhN27dwsPHjwQli1bJrRp00Z48uRJocdfuHBBMDU1Fe7fvy/ExsaKf3JyciSuXDXFbZ8gCMK0adOEIUOGCGfPnhWio6OFa9euCRcuXJCwatUVt33JyclKP7fY2Fihe/fuwoIFCySuXDXFbV9oaKhgbm4ubNmyRXj8+LEQFhYmDBo0SHB1dZW48ncrbtv8/f2FNm3aCIcOHRIeP34sHD58WGjTpo1w8uRJiStXzZkzZ4SffvpJOHr0qGBpaSkEBga+9fiUlBShU6dOwtdffy1EREQIf/zxh9CmTRth8+bNktTL4JHQ0KFDhW+//VZp26effiqsXr260OPzgyc+Pl6K8t5bcdsXEhIi2NjYVNr2vSk0NFQwNTUVrly5oo7y3ltx27dp0ybhk08+UdoWEBAgtGnTRm01llRx2zZixAjh+++/V9rm6ekpjBw5Um01lpY2bdq8M3j8/f0Fa2tr4eXLl+K2devWCQ4ODkJubq66SxR4q00i+YuU2tvbK21XZZHSoUOHwsHBAWPHjsWFCxfUWWaJlaR9J06cgIWFBfz8/NC5c2f07NkTHh4eSEtLk6LkYnmfn1++vXv3omXLlrCxsVFHie+lJO2zsbGBQqHAqVOnIAgCXrx4geDgYHTu3FmKklVWkrZlZWWhalXlFeaqVq2KmzdvIjs7W221SuXatWto27YtdHR0xG0ODg6IjY1FTEyM2j+fwSORty1SWtQ9cSMjIyxduhRr1qyBl5cXmjVrhnHjxuHy5ctSlFwsJWlfdHQ0rly5grt378LLywtubm4ICQnBwoULpSi5WErSvtelpKTgjz/+wLBhw9RV4nspSfusra3x008/Yc6cOWjdujU6duwIQRCwcuVKKUpWWUna5uDggMDAQNy4cQOCIODmzZsICAhAdnY2EhISpChbreLi4mBgYKC0Lf/7ExcXp/bP/6AeIC0PirNIqYmJCUxMTMSvra2t8eTJE2zevBnt2rVTW43vozjtEwQBMpkMP/30E/T09AAAbm5umDhxIuLi4t66hl5ZKekiswcPHsSrV68wcOBAdZRVaorTvgcPHsDDwwNffvklHBwcoFAosGrVKixZsgSrVq1Sd6nFVpy2ffnll1AoFBg1ahQEQYCBgQEGDRqETZs2QVNTU92lSuLNtgv/N7FAikWT2eORSEkXKX2TlZUVHj16VNrlvbeStM/IyAh169YVQwcAmjdvDgDl7q2w7/vz27NnD3r27ImaNWuqq8T3UpL2bdiwAZaWlvj8889hbm4OR0dHuLu748CBA3j27JkUZaukJG3T0dGBp6cnrl27hlOnTuHMmTNo2LAhdHV1UatWLSnKVitDQ8MCPZv4+HgAKNATUgcGj0RKa5HSO3fuwMjIqLTLe28laZ+NjQ1iY2OVxnSioqIAoNwtY/Q+P78bN27g7t27GD58uDpLfC8laV9GRkaBf/3nfy2UowVR3udnp6WlhXr16kFTUxPBwcHo2rUrNDQq/q/NNm3aIDQ0FJmZmeK28+fPo06dOmjUqJHaP19z6dKlS9X+KQQAqF69Ory8vGBkZAQdHR389ttvCA0NxQ8//AB9fX3MmzcPx48fx6effgoA8PPzQ1xcHDQ0NKBQKODn54e9e/diwYIFYs+gPClu+5o1a4bAwEDcunULLVq0QFRUFJYtW4Z27dphyJAhZdyagorbvnxr1qxBeno65s+fX67f/VTc9mVkZGDz5s2oVasWatasiQcPHuD777+HkZERPv/88zJujbLitu3ff//FX3/9BR0dHTx+/BgrVqzA9evX8fPPP0NfX7+MW1NQWloaIiMjERcXh71798LU1BR6enrIzs6Gnp4efvrpJ2zYsAHOzs4AAGNjY+zevRt37tyBiYkJrly5gpUrV2LKlCmSTH7hGI+E3rVI6Zu3J7Kzs7Fy5Uo8f/4cOjo6aNGiBTZu3IguXbqURfnvVNz26erqYsuWLfDw8MDQoUOhr6+PHj16YPbs2WVR/jsVt30AkJqaiuDgYHz55ZflOnSA4rdv8ODBSEtLg7+/P1auXAk9PT20b98ec+fOLYvy36q4bcvNzYWfnx/+/fdfVKlSBe3bt8fvv/8uSW+gJG7dugVXV1fxay8vL3h5ecHZ2RkrVqyAQqFAdHS0uF9PTw++vr5YtmwZhgwZgho1amDChAkYP368JPVykVAiIpJUxb9ZSUREFQqDh4iIJMXgISIiSTF4iIhIUgweIiKSFIOHiIgkxeAhIiJJMXiIiEhSDB4iIpIUg4eIiCTF4CEiIkkxeIiISFIMHiIikhSDh4iIJMXgISIiSTF4iIhIUgweIiKSFIOHKpQ7d+7Ay8sLMTExknyen58f9u3bJ8lnFeXEiRPw8vIq8fmrVq2CmZkZevbsiaJeOOzl5QUzMzNcvHix0P0XL16EmZmZUh0xMTEwMzNT+mNpaYm+fftizZo1ePnyZZE13b17F/Pnz0fXrl3RunVr2NnZYcyYMQgMDMSrV6+KPO/58+dYtWoV+vfvD2tra7Ru3RrdunXDnDlz8M8//6j4HaGyVqWsCyAqjjt37mDt2rWws7NDo0aN1P5527ZtQ8OGDTF48GC1f1ZRTpw4gaCgIEyfPr3Y5+bk5ODgwYNo0qQJHj16hMuXL8POzq5U67O3t8fAgQMBAAkJCTh27BjWrVuHa9euwdfXt8DxO3fuhIeHB/T19eHs7AwTExMkJyfjzz//xKJFi3DkyBGsXbsWcrlc6bwzZ85g1qxZyMrKQq9evTBixAhUrVoVT548wcmTJzFu3Dhs3LgRXbp0KdX2Uelj8BBVYmfOnIFCoYCfnx9mz56NwMDAUg8eY2NjMXgAYMyYMRg+fDjOnTuHW7duoXXr1uK+c+fOYdmyZTA1NYWfnx9q164t7ps4cSJ+/vlneHt7Y+nSpVi1apW47/79+5gxYwZq1KiBLVu2oHnz5ko1zJgxAwcPHoS2tnapto3Ug7faqMLw8vLCwoULAQCurq7i7Z0FCxaIx2RlZcHb2xt9+/aFhYUF2rZti6lTpyI8PFzpWoIgwM/PT7xlY2NjAycnJyxatAjZ2dkAADMzMzx58gSXLl1Sup30rtt8+/fvx9ChQ9G2bVu0adMG3bt3x+zZs/HixQul46KiojB37lw4ODiIt4xWrlyJ9PR08ZgxY8YgKChIrCf/j6q3/wICAtC4cWN06NAB/fv3x7Fjx5CamqrSuSWlqakphtujR4+U9v3000/i/74eOvm++eYbWFlZ4cCBA4iIiBC3r1mzBhkZGfDw8CgQOgAgk8kwcOBAdOzYsTSbQmrCHg9VGJ9++ikUCgV2796NqVOnwsTEBADQpEkTAEB2djYmTpyIsLAwDBw4EC4uLkhNTcWePXswatQo7NixAxYWFgCA3377DWvWrEHXrl0xcuRIaGpqIiYmBqdOnUJWVha0tLSwatUqeHp6olatWpg6dapYR2G/MPMdOHAA8+fPR9u2bfH1119DR0cHT58+xd9//434+Hjx3Fu3bmHs2LHQ19fHiBEjULduXdy9exfbt29HWFgYtm/fDi0tLUydOhW5ubkIDQ1V6gHY2Ni88/ulUCgQEhKCL774AjKZDM7OzvDz88ORI0cwYsSI4v8AiiE6OhoAUKNGDaVtt2/fho2NDVq2bFnoeTKZDEOHDsX169dx/PhxmJmZITMzE2fOnEH9+vXRuXNntdZN0mDwUIVhbm6ONm3aYPfu3ejUqRPat2+vtN/f3x+XLl3Cpk2b4OjoKG4fPXo0+vXrh1WrVmH79u0A8sZNmjdvDm9vb6VrzJkzR/z7wIED8euvv8LQ0FDpVtLbHD9+HLq6uti6dSuqVPn///f65ptvlI5btGgRjIyMEBAQgOrVq4vbO3bsiGnTpuHQoUMYPHgw7O3tcejQIYSGhqpcQ76goCC8evVKPM/c3BytWrVCQEBAqQZPZmam2JtLSEjA0aNHcfz4cdSrV0/ptt79+/cBAB9//PFbr/fRRx8BAO7duwcgr2eYlZUFc3PzUquZyhZvtVGlcfDgQZiYmODjjz/GixcvxD9ZWVno1KkTrly5goyMDABA9erV8fz5c4SGhpZqDXp6esjIyMCZM2eKnEEWERGBiIgI9OvXD1lZWUq12traQi6X49y5c+9dS2BgINq2bYvGjRuL25ydnXHjxg0xBEpDQEAAOnbsiI4dO6JPnz7w8vJC+/bt4efnpzTmkn+L7/WgLYyenh4AICUlpVjnUcXBHg9VGpGRkcjIyHjrff6EhATUr18fs2bNwldffQUXFxfUqVMHdnZ2+OSTT+Dk5PReA9RTpkzB5cuX8dVXX6FmzZqws7ND586d0bt3b/EXZ2RkJIC8MauipknHxcWVuAYACA0NRVRUFAYMGKA0zmJlZQUNDQ0EBASI42XFIZPJCmzr3r07PvvsM7x69QqPHj3Cpk2b8N9//xX4Pua3/11jTPn78wMo/7y0tLRi10vlE4OHKg1BEGBqavrWX6j5YyzW1tY4fvw4zp49i4sXL+LixYs4fPgw1q9fj507d6JmzZolqsHY2BjBwcH4559/8M8//+DSpUtYvHgx1qxZA39/f3E8CgAmTJigdEvwdfr6+iX6/HwBAQEA8gbl16xZU2D/wYMHMWfOHGhpaQEAqlatCgBij/BN+c/k5B/3unr16qFTp04AAEdHR3Tu3BkDBgzArFmzsGvXLjGs8sd1bt++/dba8yeCmJqaAsj7nmpra+Pu3btvPY8qDgYPVSiF/Ys7X9OmTZGQkIAOHTpAQ+Pdd5F1dXXh5OQEJycnAHljRMuWLUNAQAA+//zzEteora2NLl26iM+T/PXXX5g8eTK2bNkCd3d3NG3aFACgoaEh/sJ+m7e1uTCpqak4duwY7O3tMXz48AL7IyIi8Ntvv+HUqVNi2/OfiYqMjCz0OZj8Xtrrt+2K0qRJE0yYMAHr1q3D4cOH0b9/f/HcVq1aISwsDA8ePECLFi0KnCsIghiaPXr0AJAXdl26dBH/oeDg4KDKt4HKMY7xUIWS/1BhUlJSgX2DBg2CQqHAli1bCj339dtXb05tBv7/oPfr19bV1UViYqLK9RV23fzB8vzrfvTRRzA1NcWuXbvE2V+vy8nJUfrM/DarWkdwcDDS09MxcuRI9OrVq8CfyZMno1q1aggMDBTPsbe3h1wux549ewrcCktNTcXevXshl8tVnq48btw46OnpYe3atUorEcyaNQsAMHfuXCQkJBQ4z8vLC9euXcPAgQOVJhPkzxBcvHgxHj58WOhnHjp0iKsXVBCaS5cuXVrWRRCpSi6XY/v27Xj06BE0NDQQFRWFly9fol69erCwsMCNGzfw+++/48aNG1AoFLh//z6OHz+OX375BSdOnBBXIOjatStCQ0Px+PFjPHr0CH/99Rf+97//ISMjA4sXL4ahoSGAvNs+Z8+eRU5ODhQKBR48eIBGjRqJt6jeNHr0aBw5cgQxMTF4/PgxLly4gBUrViAhIQFz5syBsbExZDIZrK2tceDAAfj7+yMuLg5Pnz5FaGgoDhw4gCVLlqBevXpo1aoVgLzAOn78OJ4+fYqMjAxERkaiWrVqSlOVX/fdd98hOTkZy5YtK7ROLS0thIeH48yZMxg6dCiqV68OHR0d1KxZE/v378ehQ4eQlJSEx48f48SJE3B3d0dMTAzc3Nxga2srXic5ORnbtm2DpaVlgV5S1apVkZaWhhMnTqBJkyZiiDRt2hQ1a9ZEUFAQAgMDkZCQgGfPnuHChQtYuXIlDh48CHt7e6xatUqpdgMDA7Rq1QoHDx6Ev78/oqKiEBsbi3v37uHEiRP44YcfsHPnTvTt21fsUVL5JROKmnpDVE4FBQXBx8cHjx8/RnZ2NpydnbFixQoAeb2FnTt34sCBA+LtoTp16sDCwgLOzs7ibZqNGzfir7/+wsOHD5GSkgIDAwNYWVlhypQpStN94+PjsXTpUly8eBHJyckQBAEnT54scrmePXv24OjRo7h37x6SkpJQs2ZNtGrVChMnTkSHDh2Ujn3y5Ak2bNiAs2fPIjY2Frq6umjYsCHs7e0xevRo1K9fHwCQm5uLH3/8EUeOHIFCoUBubi48PT0LXcbnwYMH6Nu3L3r27PnW9d0OHz6M2bNnY+bMmUrPKJ07dw5+fn64ceMGUlJSoKenB0tLS4wfP77AbcGYmBh0794dLi4uWLJkSYHPSEhIQLdu3WBkZISjR49CU1NT3BceHo4tW7bg0qVLiI+PR7Vq1WBubo5BgwZh0KBBSse+7vnz5/Dz88PZs2cRExOD7Oxs1KlTB7a2thg6dGiBKfZUPjF4iIhIUhzjISIiSTF4iIhIUgweIiKSFIOHiIgkxeAhIiJJMXiIiEhSDB4iIpIU12ojKqeio6OxceNGXL58Gc+ePYO2tjaMjIzEh2H3798vvp30XaZNm4bp06eruWIi1TB4iMqhmzdvYsyYMahSpQoGDRqEFi1aICMjA1FRUTh9+jR0dXUxYsSIAmunzZs3DyYmJkqrEQB5r80mKi+4cgFROTR16lScPn0a+/fvF9dsy5ebmwuFQoG6desWOM/MzAx2dnbim1aJyiOO8RCVQ1FRUeI6b2/S0NAoNHSIKgoGD1E51KRJEyQmJuLPP/8s61KISh3HeIjKoS+++ALnz5/H9OnTYWxsDBsbG1hYWKB9+/Zo3rx5WZdH9F44xkNUTkVERGDLli34+++/ER8fL263tbXFypUrC30bKMd4qCJg8BBVAE+ePMHly5exd+9ehIaGwtTUFIGBgdDW1lY6jsFDFQHHeIgqgIYNG2LQoEHYsWMHbGxscO/ePdy4caOsyyIqEQYPUQUik8lgZWUFAIiNjS3jaohKhsFDVA6dO3cOOTk5BbZnZGTg3LlzAMBJBlRhcVYbUTnk6emJxMREdOvWDaamptDR0cF///2HQ4cOISoqCoMGDeJqBFRhMXiIyqEFCxbg5MmTuHLlCo4dO4aUlBTo6enB1NQUkyZNwuDBg8u6RKIS46w2IiKSFMd4iIhIUgweIiKSFIOHiIgkxeAhIiJJMXiIiEhSDB4iIpIUg4eIiCTF4CEiIkkxeIiISFIMHiIiktT/A5jzTuxJpy6BAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(STtest)\n",
    "plt.xlabel(\"test set AUROC\\nST\")\n",
    "plt.ylabel(\"brain areas\")\n",
    "plt.text(0.75, 400, \"mean = %.3f\" %np.mean(STtest),size=20)\n",
    "plt.xlim([0.45,1.05])\n",
    "plt.ylim([-2,420])\n",
    "plt.savefig(\"051221_2.pdf\", bbox_inches='tight')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ABA KNN all areas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "ABATRAINFILE = \"/home/slu/spatial/data/multiclass_knn/ABAknnAUROC_train_051221.npy\"\n",
    "ABATESTFILE = \"/home/slu/spatial/data/multiclass_knn/ABAknnAUROC_test_051221.npy\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "ABAtrain = np.load(ABATRAINFILE)\n",
    "ABAtest = np.load(ABATESTFILE)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgEAAAE8CAYAAABZ6mKPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de1zO9/8/8MfVSevkkMaIT0LlUCnHlFMoc1glGSKnsWbzYWPWZmGWQz4zJraI5Mx0pIVhs2U5RQ4j5hQV00kHRcf37w+/rq9L0VWuq6t6P+63W7eb6/V+vd/X81Wfz67H9Xq/3u+3RBAEAURERCQ6aqougIiIiFSDIYCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiERKQ9UFEBER1QW///47goODce3aNZSVlaFDhw6YMGEC3NzcqnWcoqIi7N+/HxEREUhOTkZRURFatmwJe3t7TJ06Fa1bt650v4iICOzatQu3b9+GmpoaOnfujGnTpmHQoEGV9i8tLcWOHTsQFhaGe/fuQVtbG9bW1vjoo49ga2srV60S3iyIiIjEbufOnfj222/RpEkTDB8+HJqamjhy5Aj+/fdfTJs2DV988YVcxykpKcGkSZNw4cIFmJqaom/fvtDS0sKVK1dw7tw56OvrY+/evejQoYPMfv7+/ggODkbLli3h7OyM4uJixMTEIDs7G76+vpg4caJMf0EQMGfOHBw5cgTt2rXDoEGDkJOTg0OHDqGwsBDr1q3DkCFDqi5YICIiErHk5GSha9euQq9evYTk5GRpe3Z2tjBkyBDBzMxMuHDhglzHiomJEczMzITJkycLpaWlMtt++OEHwczMTPDx8ZFpP3/+vGBmZiYMGTJEyM7OlqmrV69eQteuXWXqEgRBOHjwoGBmZia8//77wrNnz6Ttly5dErp06SL06dNHyMvLq7JergkgIqoHUlJSYG5uDh8fH9y/fx///e9/0bt3b9jY2GDatGn4559/AABZWVnw9fWFg4MDLC0t4e7ujtOnT1d6zJKSEuzatQtjx46Fra0trK2t4erqip07d6KsrKxC//DwcMyePRuDBw+GlZUVbG1tMW7cOERFRVV6/EmTJsHc3BwlJSUIDAyEk5MTunbtigEDBuB///sfioqKFPcLegNhYWEoKiqCp6cnjI2Npe2NGzfGhx9+CADYu3evXMdKTk4GAAwcOBBqarIfsYMHDwYAPH78WKa9/Nje3t5o3LixtN3Y2BgTJkxAUVERwsPDZfbZs2cPAGDu3Llo1KiRtN3KygrDhw9HVlYWjhw5UmW9DAFERPVIamoqPDw8kJGRATc3Nzg4OCAuLg6TJk1CUlISxo4diytXruDdd9/Fu+++ixs3bmDGjBl48OCBzHGKi4vx4YcfYunSpcjNzcXIkSMxduxYlJWV4dtvv610+nvJkiVITU1Fz549MXnyZIwYMQIPHjzAggULsHbt2lfWPG/ePOzcuRPdu3fH+PHjoa2tjc2bN2Px4sUK//3URHlI6tevX4Vt/fv3l+lTlY4dOwIA/vzzzwpB6vfffwcA2NnZvdH7FxUVISEhAW+99RZ69OjxRjVzYSARUT1y9uxZzJ07Fx999JG0bcOGDVi3bh3Gjh2LYcOGYcmSJdJvoX379sUXX3yBkJAQfPXVV9J9AgMDcfLkSUycOBFfffUV1NXVATxfbObr64uwsDA4OzvLnFeOjo5G27ZtZeopKirCjBkzEBQUhPHjx6NFixYVak5OTkZ0dDSaNGkCAPj000/h4uKCyMhIfPbZZzAyMqpy3ImJiTh27Fg1flPA5MmTYWBgUGW/u3fvAgBMTEwqbHv77beho6ODf//9F0+fPsVbb7312mMNHDgQTk5O+PXXXzFq1CjY2dlBU1MTV69exYULFzBp0iSZ8/sFBQV49OgRdHR08Pbbb1c43n/+8x8AQFJSkrTt3r17KC0tRZs2baChUfFjvLJ9XoUhgIioHmndujVmzpwp0+bm5oZ169ahqKgICxYskJmGHjVqFBYuXIjExERpW1lZGXbu3AkjIyN8+eWX0gAAAOrq6vDx8UF4eDgOHjwoEwJeDgAAoKWlBU9PT5w+fRqnTp2Cq6trhT7z58+XBgAA0NHRwahRo7Bhwwb8/fffr1z9/qLExESsX7++yn4vcnNzkysEPHnyBACgr69f6XY9PT0UFBQgLy+vyhAgkUiwbt06bNiwAT/++CNu3bol3WZnZ4eRI0fK/L7z8vJe+97l7bm5uRX20dPTk3ufV2EIICKqRzp16iTzIQJA+g3SxMSkwgeDuro6DA0N8ejRI2nb3bt3kZ2dDRMTE/z000+Vvo+2tjbu3Lkj0/bgwQMEBQXh1KlTePjwIZ49eyaz/cX3eFHXrl0rtL3zzjsAgJycnEr3edno0aMxevRoufoqmvD/L6KTSCRV9i0sLMSCBQvw559/YtGiRRg8eDDeeustnD9/HsuWLcPEiROxdu1a+Vbuv0Ce965JvQwBRET1SGXfGMunhF/1bVJDQwMlJSXS19nZ2QCeTxe/7tt1fn6+9N/JyckYM2YMcnNz0aNHDzg4OEBPTw/q6upITU1FRETEKxf6VfZtvDzIVLYAsbbp6enh8ePHyMvLQ9OmTStsL/89vOqb94s2bdqEw4cPY+HChRg3bpy0fcCAAWjRogVcXFywfPlyaQgo/5uVf7t/WWUzBeX/Lp/BeFlVMxsvYgiogrm5OQDgxo0bKq6EiEgxyj8chg4dKvcU+9atW5GdnY0VK1ZU+EYeHR2NiIgIhdf5ImWuCWjXrh0eP36MpKSkCiEgLS0NBQUFaNmyZZWnAoD/W/zXu3fvCtssLCzQpEkTpKam4vHjx2jatCl0dHTQokULPHr0CGlpaRXWBdy7dw+A7HqF//znP1BXV0dycjJKSkoqrAuobJ9XqRch4Ny5c9iyZQuuXr2KtLS0Cv8jFAQB69evx759+5Cbmwtra2ssWrRIukoTeD7l5Ofnh99++w0A4OjoCF9fX7n+B0JE1JCYmprCwMAAFy9eRHFxMTQ1Navcp/yDxcnJqcK2s2fPKrzGlylzTUCfPn1w4cIFxMbGwsbGRmbbn3/+Ke0jj+LiYgDPL9V8WVFRkfRb+ou/8z59+iAqKgqxsbFwd3ev8v21tLRgY2OD+Ph4xMfHV6itOjXXi0sECwoKYGZmhoULF0JbW7vC9qCgIAQHB8PX1xehoaFo1qwZpk6dKjNVMm/ePFy7dg1BQUHYvHkzrl27hgULFtTmMIiI6gQNDQ1MnDgR6enp8PPzq3BuH3j+DfjFRW3lt7p9+QM/NjYWoaGhyi0Yz9cE3Lhxo1o/L17zX9WxtbS0sGvXLqSkpEjbc3JysHHjRgCQmdoHnk/T3759G2lpaTLt3bt3BwBs3LixwumRgIAAlJSUwNLSUubUQvmxAwMDZdZIpKSkYPfu3dDS0qow+zJ+/HgAwNq1a1FYWChtv3z5MmJiYtCsWTM4OztXOfZ6MRMwYMAADBgwAADw5ZdfymwTBAHbt2/HzJkzpQP29/eHnZ0doqOjMW7cONy+fRuxsbHYvXu39H7K33zzDTw9PXHnzh2YmprW7oCIiFRs1qxZuH79Ovbu3Yvff/8dffr0QYsWLZCZmYl79+7hwoUL+PTTT6W3t50wYQLCw8MxZ84cODk5oUWLFrh58yZiY2Px7rvvIiYmRsUjqrk2bdpgwYIF8PPzg7u7e6W3DX55huDo0aP48ssv4ebmhpUrV0rbP/roI/z+++84deoUhg0bhn79+kFbWxsXLlzA5cuXoa2tjYULF8ocy9bWFlOnTsXWrVvx3nvvVXrb4JcDzYgRI/Drr7/iyJEjcHV1xaBBg5CdnY1Dhw5J7/UgzxqGehECXiclJQXp6emwt7eXtmlra6Nnz55ISEjAuHHjkJCQAB0dHZkHKnTv3h06OjpISEhgCCAi0dHU1MSPP/6IqKgoRERE4MSJEygoKEDTpk1hbGyMOXPmYNSoUdL+FhYW2L59O9auXYs///wTJSUlsLCwwPr166Gvr1+vQwDw/O6GrVu3RnBwMCIjIyEIAtq3b4+5c+dW6wFCLVq0QEREBIKCgnDixAmEh4dDEAQYGRlh9OjR+OCDD9C+ffsK+/n4+MDc3Bw7d+7Ezz//DIlEgi5dumD69OmVXkIpkUjw/fffY+fOnQgLC8POnTvRqFEj9OjRo2E/QMjGxga+vr7SqZELFy5g/Pjx+P3339GqVStpvy+//BJpaWnYsmULAgMDsX//fhw/flzmWIMHD8bYsWOlt4WsDBcGEhFRQ1Uv1gTIo6rrISvbLghCta69JCIiakjqfQgov91kenq6THtmZiaaN28OAGjevDkyMzPx4qSHIAh4/PgxDA0Na69YIiKiOqTehwBjY2MYGRkhLi5O2lZYWIj4+HjpQg4bGxsUFBQgISFB2ichIQEFBQUVFnsQERGJRb1YGJifn4/79+8DeH53qQcPHiAxMRGNGzdGq1at4OXlhcDAQJiamkpvg6mjo4ORI0cCANq3b49+/fph8eLF+PbbbyEIAhYvXoxBgwZxUSAREYlWvVgYeObMGXh5eVVoL78048WbBeXk5EhvFmRmZibtm52dXeFmQYsWLaryRhJcGEhERA1VvQgBqsQQQFS3mPj8ouoSal3SyhGqLoEaqHq/JoCIiIhqhiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikGAKIiIhEql7cMZCI6hYxXqtP1BBxJoCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQaRAgoLS3F2rVr4ejoCEtLSzg6OmLNmjUoKSmR9hEEAQEBAXBwcICVlRUmTZqEmzdvqrBqIiIi1WoQISAoKAi7d+/G119/jUOHDmHhwoXYvXs3Nm7cKNMnODgYvr6+CA0NRbNmzTB16lQ8efJEhZUTERGpToMIAQkJCRg0aBAcHR1hbGyMwYMHw9HREZcvXwbwfBZg+/btmDlzJpydnWFmZgZ/f3/k5+cjOjpaxdUTERGpRoMIAd27d8eZM2dw+/ZtAMCtW7dw+vRp9O/fHwCQkpKC9PR02NvbS/fR1tZGz549kZCQoJKaiYiIVE1D1QUowowZM5Cfn48RI0ZAXV0dJSUl8Pb2hqenJwAgPT0dANC8eXOZ/QwNDZGWllbr9RIREdUFDSIExMTEIDIyEqtXr0aHDh2QmJiI5cuXw9jYGB4eHtJ+EolEhVUSERHVLQ0iBKxatQrTpk3DiBEjAADm5uZ48OABNm3aBA8PDxgZGQF4PiPwzjvvSPfLzMysMDtAREQkFg1iTcCzZ8+grq4u06auro6ysjIAgLGxMYyMjBAXFyfdXlhYiPj4eNjY2NRqrURERHVFg5gJGDRoEDZt2gRjY2Pp6YCtW7fC1dUVwPPTAF5eXggMDISpqSlMTEzw008/QUdHByNHjlRx9URERKrRIELA119/jR9++AHffPMNMjMzYWRkhLFjx+Ljjz+W9pkxYwYKCwuxdOlS5OTkwNraGsHBwdDT01Nh5URERKojEQRBUHURdZm5uTkA4MaNGyquhKjuMPH5RdUliErSyhGqLoEaqAaxJoCIiIiqjyGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEik5L5PwOPHj5GVlYX27dtL25KTkxESEoLs7Gy4urqiX79+SimSiIiIFE/uELBs2TIkJSUhNDQUAJCfnw9PT0/pU/gOHTqEbdu2oWfPnsqplIiIiBRK7hBw8eJFvPfee9LXMTExSEtLw6ZNm9CpUydMmzYNmzdvZgggqkW8aQ8RvQm51wRkZmbKPIEvNjYWXbt2Rf/+/WFkZAQ3Nzdcu3ZNKUUSERGR4skdAjQ0NFBYWCh9ffbsWZlv/fr6+sjOzlZsdURERKQ0cocAExMTHDlyBIIg4Pjx48jJyYGdnZ10+7///ovGjRsrpUgiIiJSPLnXBHh6esLHxwc9e/bEs2fP0KZNG5kQcO7cOenDdoiIiKjukzsEuLq6QiKR4NixY9DT04O3tzc0NTUBPL988MmTJxg/frzSCiUiIiLF4qOEq8BHCVNdxqsDxIGPEiZl4R0DiYiIREru0wEAUFJSgmPHjuHSpUvIzc1FWVmZzHaJRILly5crtEAiIiJSDrlDQHZ2Nry8vHDz5k0IggCJRILyMwnl/2YIICIiqj/kPh2wdu1a3LlzB35+fjh69CgEQcCWLVsQExODESNGwNLSEmfOnFFmrURERKRAcoeAP/74A66urnB3d4eent7zndXUYGpqiu+++w7a2tpYvXq10golIiIixZI7BKSnp8PS0hLA87sHAkBRUZF0++DBg/Hbb78puDwiIiJSFrlDQJMmTfD06VMAgK6uLjQ0NPDw4UPpdk1NTeTm5iq+QiIiIlKKat02+NatW893UlND586dERERgaKiIjx9+hSRkZFo06aN0golIiIixZI7BNjb2+PIkSPSUwBTpkzBpUuX0KtXL/Tt2xd///03Jk+erLRCiYiISLHkvkTQ29sb06dPh5aWFgBg+PDh0NDQwIEDB6CmpoZhw4Zh+PDhSiuUiIiIFEvuECCRSKQBoJyTkxOcnJwUXhQREREpX41uG1xUVIRHjx7JXB1ARERE9Uu1QsDVq1fh5eUFW1tbDBw4EOfPnwcAZGZmYvLkyYiLi1NKkURERKR4coeAxMREeHp6Ijk5GS4uLjLbDA0NUVhYiIiICIUXSERERMohdwj44Ycf8PbbbyM6Ohrz5s3Dy08g7tOnDy5fvqzwAomIiEg55A4B58+fh4eHB3R1dSGRSCpsb9WqFdLS0hRaHBERESmP3CGgsLAQ+vr6r9z+5MkThRREREREtUPuENC2bVtcvXr1ldtPnz6NDh06KKSomkhLS8MXX3yBPn36wNLSEsOHD8fZs2el2wVBQEBAABwcHGBlZYVJkybh5s2bKquXiIhI1eQOASNHjkRUVJTMFQDlpwWCg4MRGxtbYcFgbcnNzcX48eMhCAI2bdqEmJgY+Pr6wtDQUNonKCgIwcHB8PX1RWhoKJo1a4apU6dyBoOIiERL7psFTZs2DX/99RemT58OU1NTSCQSrFixAllZWcjIyEDfvn0xYcIEZdb6Sps3b4aRkRFWrVolbXvxOQaCIGD79u2YOXMmnJ2dAQD+/v6ws7NDdHQ0xo0bV+s1ExERqZrcMwFaWlrYunUrvvjiCzRq1AiNGjVCUlISmjZtis8//xwbN26EmlqN7j30xo4dOwZra2vMnTsXdnZ2cHFxwc6dO6VXMKSkpCA9PR329vbSfbS1tdGzZ08kJCSopGYiIiJVk3smAAA0NDQwZcoUTJkyRUnl1ExycjJ2796NKVOmYObMmUhMTISfnx8AYOLEiUhPTwcANG/eXGY/Q0NDXtFACmHi84uqSyAiqja5vrrn5+djyJAhCAkJUXI5NSMIArp06YJ58+ahc+fOcHd3x6RJk7Br1y6ZfpVd2khERCRWcoUAXV1dZGdnQ1dXV9n11IiRkRHat28v02ZqaoqHDx9KtwOQzgiUy8zMrDA7QEREJBZyn8S3trbGlStXlFlLjdna2uLu3bsybUlJSWjVqhUAwNjYGEZGRjJXNhQWFiI+Ph42Nja1WisREVFdIXcImD9/Pg4fPoywsLAKtwxWtcmTJ+PSpUv46aefcO/ePRw6dAg7duyAp6cngOenAby8vLBp0yb8+uuv+Oeff+Dj4wMdHR2MHDlSxdUTERGphkSQ8xPdy8sLDx48QGpqKho3boy2bdtCW1tb9mASCbZt26aUQqty4sQJfP/997h79y5atWoFT09PTJo0SboOQBAErF+/Hvv27UNOTg6sra2xaNEimJmZvfa45ubmAIAbN24ofQxUf3FhIClT0soRqi6BGii5Q4Cjo6NcB/ztt9/eqKC6hiGA5MEQQMrEEEDKIvclgg3tw52IiEjsVHN3HyIiIlI5hgAiIiKRqtYdA+/fv4+QkBBcunQJubm5KCsrk9kukUhw7NgxhRZIREREyiH3TMCNGzfg5uaG/fv3o7i4GMnJydDR0UFhYSFSU1Ohrq6Od955R5m1EhERkQLJHQLWrVsHTU1NREVFSW8f/NVXX+HkyZNYunQpcnNzsXjxYmXVSURERAomdwg4f/483n//feljhF80duxY9O/fH999953CCyQiIiLlkDsE5Ofno02bNgAATU1NAEBBQYF0u62tLS5cuKDg8oiIiEhZ5A4BzZs3R0ZGBgBAT08Pb731FpKSkqTbc3NzUVpaqvACiYiISDnkvjrAwsJC5gFCvXr1wvbt22FlZYWysjLs3LkTFhYWSimSiIiIFE/umYBRo0YhOzsbz549AwDMmTMHeXl58PLywpQpU5CXl4dPP/1UaYUSERGRYsn97IDKPHz4EEePHoW6ujr69+8vXTPQkPDZASQPPjuAlInPDiBlqdbNgl72zjvvwMvLS1G1EBERUS3ibYOJiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikanSJ4NOnT5GdnY3KbjHQqlWrNy6KiIiIlE/uEFBWVobNmzdjx44d0mcIVCYxMVEhhREREZFyyR0CvvvuOwQHB6Njx45wdnZGkyZNlFkXERERKZncIeDAgQPo168fgoKClFkPERER1RK5Fwbm5uZi8ODByqyFiIiIapHcIcDMzAzp6enKrIWIiIhqkdwh4JNPPsHevXvx8OFDZdZDREREtUTuNQF///03WrVqheHDh2Po0KEwNjaGmppshpBIJPj4448VXiQREREpntwhYP369dJ/HzhwoNI+DAFERET1h9wh4Pjx48qsg4iIiGqZ3CGgdevWyqyDiIiIahmfHUBERCRSr5wJWL9+PSQSCT766COoqanJrAl4Fa4JICIiqj8kQmVPAQJgYWEBiUSCS5cuQUtLCxYWFlUfTCKpE88OCAwMxJo1a+Dp6YlFixYBAARBwPr167Fv3z7k5ubC2toaixYtQseOHV97LHNzcwDAjRs3lF431V8mPr+ougRqwJJWjlB1CdRAvXImoHwhoJaWlszruu7ixYv4+eefpR/e5YKCghAcHIyVK1eiXbt22LBhA6ZOnYrDhw9DT09PRdUSERGpzitDwMsLAevDwsC8vDzMnz8fy5Ytw48//ihtFwQB27dvx8yZM+Hs7AwA8Pf3h52dHaKjozFu3DhVlUxERKQyDWphoK+vL5ydnWFnZyfTnpKSgvT0dNjb20vbtLW10bNnTyQkJNR2mURERHWC3JcIlrty5QouX76MnJwclJWVyWxT5cLAn3/+Gffv38eqVasqbCt/5kHz5s1l2g0NDZGWllYr9ZHy8bw8EVH1yB0Cnj17hk8++QR//fUXBEGARCJB+ZrC8n+rKgTcuXMH33//PXbt2iVdw1AZiURSi1URERHVbXKHgA0bNuCvv/6Ct7c37Ozs4OXlhZUrV8LQ0BCbNm3Cs2fP4O/vr8xaX+nixYt4/PgxRo0aJW0rLS3FuXPnsHfvXkRHRwN4PiPwzjvvSPtkZmZWmB0gIiISC7nXBBw5cgTDhg3DnDlzpJfVtWjRAv369cPWrVtRXFyMiIgIpRX6OkOGDMHBgwcRGRkp/enatStGjBiByMhItGvXDkZGRoiLi5PuU1hYiPj4eNjY2KikZiIiIlWTeybg4cOHmDJlCgBAXV0dAFBcXPz8IBoaGDFiBPbs2YN58+YpvsoqGBgYwMDAQKZNR0cHjRs3hpmZGQDAy8sLgYGBMDU1hYmJCX766Sfo6Ohg5MiRtV4vERFRXSB3CNDV1UVpaan032pqajKL6vT19ZGRkaH4ChVkxowZKCwsxNKlS5GTkwNra2sEBwfzHgFERCRacoeAtm3bIikpCcDzmYAOHTrgyJEjGDNmDARBwNGjR2XOt6vajh07ZF5LJBLMnj0bs2fPVlFFREREdYvcawLs7Oxw5MgR6WzA+++/j9jYWAwZMgROTk6Ii4uDu7u70golIiIixZJ7JmDmzJlwcXGRXhbo6emJoqIiHDhwAGpqavj0008xY8YMpRVKREREivXKBwjRc3yAUP3BmwVRQ8UHCJGyyHU6ID8/H0OGDEFISIiSyyEiIqLaIlcI0NXVRXZ2NnR1dZVdDxEREdUSuRcGWltb48qVK8qshYiIiGqR3CFg/vz5OHz4MMLCwsBlBERERPWf3FcHrFixAgYGBvj666/xv//9D23btoW2trZMH4lEgm3btim8SCIiIlI8uUNASkoKAEhvCFSX7w5IREREVZM7BPz222/KrIOIiIhqmdxrAoiIiKhhkXsmoFxRURHOnDmD5ORkAECbNm3Qq1cvNGrUSOHFERERkfJUKwRERkZixYoVyM3NlV4hIJFIYGBggC+++AKjR49WSpFERGKmqrth8k6FDZ/cISAmJgY+Pj5o1aoVpk+fjvbt20MQBNy+fRt79+7FwoULoa2tjeHDhyuzXiIiIlIQuZ8d8N5776GkpAQ///wz9PT0ZLbl5eXBw8MDmpqaOHjwoFIKVRU+O6D+4LMDiBSLMwENn9wLA+/evYvRo0dXCAAAoK+vj9GjR+PevXsKLY6IiIiUR+4QYGRk9No7BaqpqaF58+YKKYqIiIiUT+4Q4ObmhoiICOTn51fY9uTJE4SFhXFhIBERUT3yyoWB586dk3ndo0cP/P777xg1ahQmTJgAU1NTSCQS3Lp1C3v27EHTpk3RvXt3pRdc33BVLxER1VWvDAGTJk2CRCKRaSs/HfDdd99Jt5W3PXjwANOmTUNiYqKyaiUiIiIFemUIWLFiRW3WQURERLXslSHAzc2tNusgIiKiWsZnBxAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiVS1HiVMJA8+yIeIqH7gTAAREZFIMQQQERGJFEMAERGRSDWIELBx40a4u7vD1tYWffr0gbe3N/7550IBRvgAABtUSURBVB+ZPoIgICAgAA4ODrCyssKkSZNw8+ZNFVVMRESkeg0iBJw9exYTJkzA3r17sW3bNqirq2Pq1KnIzs6W9gkKCkJwcDB8fX0RGhqKZs2aYerUqXjy5IkKKyciIlKdBhECtmzZAnd3d5iZmcHc3ByrVq1CVlYWLly4AOD5LMD27dsxc+ZMODs7w8zMDP7+/sjPz0d0dLSKqyciIlKNBhECXpafn4+ysjIYGBgAAFJSUpCeng57e3tpH21tbfTs2RMJCQmqKpOIiEilGmQIWLZsGTp16gQbGxsAQHp6OgCgefPmMv0MDQ2RkZFR6/URERHVBQ3uZkErVqzA+fPnsWfPHqirq8tsk0gkKqqKiIio7mlQMwHLly/HL7/8gm3btqFNmzbSdiMjIwD/NyNQLjMzs8LsABERkVg0mBDg5+eH6OhobNu2De3bt5fZZmxsDCMjI8TFxUnbCgsLER8fLz1lQEREJDYN4nTAN998g6ioKGzYsAEGBgbSb/w6OjrQ1dWFRCKBl5cXAgMDYWpqChMTE/z000/Q0dHByJEjVVw9ERGRajSIELB7924AwJQpU2TaP/nkE8yePRsAMGPGDBQWFmLp0qXIycmBtbU1goODoaenV9vlEhER1QkNIgTcuHGjyj4SiQSzZ8+WhgIiIiKxaxAhgCri43yJiKgqDWZhIBEREVUPQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYmU6ELArl274OjoCEtLS4wePRrx8fGqLomIiEglRBUCYmJisHz5cnh7eyMyMhI2NjaYMWMGHjx4oOrSiIiIap2oQsDWrVvh5uaGsWPHon379vD19YWRkRH27Nmj6tKIiIhqnYaqC6gtRUVFuHr1KqZNmybTbm9vj4SEhCr3Nzc3r9H7NqrRXkREqmce8dkb7X/jxg0FVULKIpqZgMePH6O0tBTNmzeXaTc0NER6erqKqiIiIlId0cwElJNIJHK1lWOSJSKihko0MwFNmzaFurp6hW/9mZmZFWYHiIiIxEA0IUBLSwtdunRBXFycTHtcXBxsbGxUVBUREZHqiOp0wNSpU7FgwQJYWVnB1tYWe/bsQVpaGsaNG6fq0oiIiGqdqELA8OHD8fjxY/z0009IS0uDmZkZNm3ahNatW6u6NCIiolonEQRBUHURREREVPtEsyaAiIiIZDEEEBERiRRDwBuqzgOJzpw5A3Nz8wo/t2/frsWK31x1H8JUVFSEH374AY6OjujatSsGDhyI7du311K1b6464/Xx8an0b9ytW7darPjNVfdvfPDgQbi4uMDa2hr29vaYP39+vbsJV3XHvGvXLrz77ruwsrKCs7MzIiMja6lSxTh37hy8vb3Rr18/mJubIzw8vMp9bty4gYkTJ8LKygr9+vXD+vXrwTPK9ZxANfbLL78InTt3Fvbt2yfcunVLWLp0qdCtWzchNTW10v6nT58WzMzMhJs3bwppaWnSn5KSklquvOaqO2ZBEIRPPvlEcHd3F06ePCkkJycLFy9eFE6fPl2LVddcdcebm5sr87dNS0sTBg8eLPj4+NRy5TVX3THHx8cLFhYWwtatW4X79+8LCQkJgqurq+Dl5VXLlddcdce8a9cuoVu3bsLBgweF+/fvC9HR0UK3bt2E48eP13LlNXfixAlh9erVwqFDhwQrKyshLCzstf3z8vKEvn37Cv/973+FGzduCIcPHxa6desmbNmypZYqJmVgCHgDY8aMERYuXCjTNnToUOG7776rtH95CMjMzKyN8pSiumOOjY0VbG1t6+2Yqzvel8XHxwtmZmbC+fPnlVGeUlR3zJs3bxYGDhwo0xYaGip069ZNaTUqWnXH/P777wvLli2TaVuxYoUwbtw4pdWoTN26dasyBOzatUuwsbERnj59Km3bsGGD4ODgIJSVlSm7RFISng6oofIHEtnb28u0y/NAojFjxsDBwQGTJ0/G6dOnlVmmQtVkzMeOHYOlpSVCQkLQv39/ODk5wc/PD/n5+bVR8ht5k79xuf3796Njx46wtbVVRokKV5Mx29raIj09Hb/99hsEQUBWVhZiYmLQv3//2ij5jdVkzEVFRWjUSPbxYI0aNcKVK1dQXFystFpV6eLFi+jRowe0tbWlbQ4ODkhLS0NKSooKK6M3wRBQQzV5IJGRkRGWLFmCdevWISAgAO3atcOUKVNw7ty52ij5jdVkzMnJyTh//jyuX7+OgIAA+Pr6IjY2Fl9++WVtlPxG3vShU3l5eTh8+DA8PDyUVaLC1WTMNjY2WL16NebPn4+uXbvCzs4OgiDA39+/Nkp+YzUZs4ODA8LCwnD58mUIgoArV64gNDQUxcXFePz4cW2UXesyMjJgaGgo01b+O8vIyFBFSaQAorpZkDJU54FEpqamMDU1lb62sbFBamoqtmzZgp49eyqtRkWrzpgFQYBEIsHq1auhr68PAPD19cX06dORkZFRL57bUN2HTpU7cOAASktL4eLiooyylKo6Y7516xb8/Pwwa9YsODg4ID09HatWrcKiRYuwatUqZZeqMNUZ86xZs5Ceno7x48dDEAQYGhrC1dUVmzdvhrq6urJLVZmXfx/C/18UKM//H6hu4kxADSnqgUTW1ta4d++eostTipqM2cjICC1atJAGAABo3749AODBgwfKK1YB3vRv/PPPP8PJyQlNmjRRVokKV5Mxb9y4EVZWVvjggw9gYWGBfv36YfHixYiKisLDhw9ro+w3UpMxa2trY8WKFbh48SJ+++03nDhxAq1bt4auri6aNm1aG2XXuubNm1f4xp+ZmQkAFWYIqP5gCKghRT2QKDExEUZGRoouTylqMmZbW1ukpaXJrAFISkoCgDp/u+Y3+RtfvnwZ169fx9ixY5VZosLVZMzPnj2r8O23/LVQDy4fe5O/s6amJlq2bAl1dXXExMRg0KBBUFNrmP9Z7datG+Lj41FYWChti4uLw9tvvw1jY2MVVkZvQn3JkiVLVF1EfaWnp4eAgAAYGRlBW1sbP/74I+Lj47F8+XIYGBhgwYIFOHr0KIYOHQoACAkJQUZGBtTU1JCeno6QkBDs378fPj4+0m/HdV11x9yuXTuEhYXh77//RocOHZCUlISlS5eiZ8+ecHd3V/Foqlbd8ZZbt24dCgoK8MUXX9S7qdLqjvnZs2fYsmULmjZtiiZNmuDWrVtYtmwZjIyM8MEHH6h4NPKp7pjv3r2LP/74A9ra2rh//z5WrlyJS5cuYc2aNTAwMFDxaOSTn5+P27dvIyMjA/v374eZmRn09fVRXFwMfX19rF69Ghs3boSbmxsAwMTEBPv27UNiYiJMTU1x/vx5+Pv748MPP6w3C1+pIq4JeANVPZDo5anQ4uJi+Pv749GjR9DW1kaHDh2wadMmDBgwQBXl10h1x6yrq4utW7fCz88PY8aMgYGBAYYMGYJ58+apovxqq+54AeDJkyeIiYnBrFmz6l0AAKo/5tGjRyM/Px+7du2Cv78/9PX10bt3b3z++eeqKL9GqjvmsrIyhISE4O7du9DQ0EDv3r2xZ8+eevWN+O+//4aXl5f0dUBAAAICAuDm5oaVK1ciPT0dycnJ0u36+voIDg7G0qVL4e7ujsaNG2PatGmYOnWqKsonBeEDhIiIiESqYZ68IiIioioxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBAREYkUQwAREZFIMQQQERGJFEMAEQAfHx+Ym5urugwiolrFEEB1WmJiIgICApCSkqLqUuqMgIAAHDt2rEb7lpaWwsHBAebm5vjxxx9f2c/R0RGOjo6v3F4eml78uwQEBMDc3Fz6Y2FhgV69emHy5Mk4fvz4K49VVlaGiIgIeHl5oXfv3ujatSsGDhyIzz//HImJia8dz6lTpzBnzhwMGDAAXbt2hY2NDUaPHo01a9bg33//fe2+RMQQQHVcYmIi1q9fj9TUVKW+z7fffovLly8r9T0UZf369TUOAX/++SfS09PRtm1bhIeHQxAEBVcH/Pe//8WqVauwfPlyTJgwAf/88w9mzZqFAwcOVOhbUFCADz74AD4+Pnj27BlmzJiBxYsXY8SIETh58iTc3d2xe/fuCvuVlZXh66+/xpQpU3DlyhWMHDkSixcvxvz582FhYYFdu3bBzc1N4WMjamg0VF0AkSKVlpaiqKgIb731VrX209TUVFJFdUtoaCjatGkDHx8fzJo1C2fOnEGfPn0U+h79+/eHpaWl9PWwYcPg4uKCoKAgvPfeezJ9Fy9ejL/++gve3t749NNPZbZNnz4dU6ZMwdKlS2FiYoK+fftKtwUEBGD//v0YOXIkVqxYAS0tLZl9fXx8sH79eoWOi6gh4kwA1VkBAQH48ssvAQBeXl7SaWYfHx8AQHh4OMzNzREXF4cNGzZgyJAhsLKywqFDhwAAJ0+exNy5czF48GBYWVmhR48emDZtGs6ePVvhvSpbE1DelpeXh8WLF8POzg6WlpYYN24cLl26JNcYsrOzsXz5cgwZMgSWlpbo3bs3Ro8ejc2bN1foGxMTg/Hjx8PGxgbW1tbw8PDA4cOHpdtTUlKkNUZERMhMvcsjMzMTf/zxB1xdXTFgwAAYGhoiNDRUrn3fhIWFBZo2bYqkpCSZ9uvXr+PAgQOwtrbG3LlzK+zXrFkzrF69GgDw3XffSdszMzOxZcsWtG7dGsuXL68QAADAwMAAX331lWIHQtQAcSaA6qyhQ4ciPT0d+/btg7e3N0xNTQEAbdu2lenn7++PkpISjB07Frq6umjXrh2A5x+UOTk5cHV1RcuWLfHo0SPs378fU6ZMwfbt29GjRw+56pg+fTqaNWuGjz/+GNnZ2di6dStmzpyJ48ePQ09P77X7zpkzB/Hx8Xj//fdhYWGBp0+f4s6dOzh79iw++OADab81a9YgMDAQ/fr1w5w5c6CmpoajR49izpw5WLRoETw9PdGsWTOsWrUKCxYsQI8ePTB27Njq/DoRERGBkpISuLi4QENDA6NGjcLevXuRl5cHfX39ah2rOnJycpCbm4tmzZrJtP/6668AAA8PD0gkkkr37dixI7p164aEhASkpqaidevWOHHiBAoLC+Hi4oJGjRoprW4iMWAIoDrLwsIC3bp1w759+9C3b1/07t270n7Pnj1DZGRkhVMA3377LXR0dGTaxo0bhxEjRmDjxo1yh4DOnTtjyZIl0tft27fH3LlzER0djXHjxr1yv7y8PJw+fRrjx4/HokWLXtnv6tWrCAwMxIcffojPPvtM2u7l5YVZs2Zh9erVcHFxgZ6eHlxcXLBgwQK0adMGLi4uctVfLjw8HD179kSbNm0AAG5ubggJCcHBgwcxYcKEah3rdZ48eYKsrCyUlpbi/v37WLt2LUpLSyucCrh58yaA57/f1+nSpQsSEhLwzz//oHXr1tL9OnXqpLCaicSKIYDqvfHjx1e6BuDFAJCfn4+ioiKoqanB2tpa7ul8AJgyZYrM6/Jz6Pfu3Xvtfo0aNYKWlhYuX76MlJQUGBsbV9rv4MGDkEgkcHV1RVZWlsw2R0dHHD9+HBcvXoSDg4PcNb/swoULuH37NqZPny5ts7CwQKdOnRAWFqbQEPDy70tLSwszZsyoMOX/5MkTAKhyFqJ8tiUvL09mv6pmYYioagwBVO+VT/+/7P79+1izZg1OnjyJ3NxcmW2vmn6uTPk353JNmzYF8Px8/+toaWnhq6++wrJlyzB48GB06NABffr0wZAhQ2BnZyftd/v2bQiCgHffffeVx8rIyJC73sqEhoZCU1MTnTt3lgkvDg4OCAoKwvXr12FhYVHt41b2e1y0aBHatWuHp0+f4syZM9ixYwdyc3OhoSH7n5uXP9xfJT8/H8D/hYXy/crbiajmGAKo3tPW1q7Qlp+fD09PTzx9+hSTJ0+GmZkZdHV1oaamho0bN+L06dNyH19dXb3Sdnkurxs/fjwGDx6MP/74A2fPnsWRI0ewc+dODB8+HGvWrJEeRyKRICgo6JXv1aFDB7nrfVl+fj4OHTqE4uJiuLq6VtonLCwMCxculL5u1KjRa0NOQUEBgMp/91ZWVtKrAwYPHozmzZtj9erV6NSpE8aPHy/t17FjR/z666+4du0aunTp8sr3unr1KgDAzMxMuh8AXLt2DUOHDn3lfkRUNYYAqtOq8439RadOnUJaWhqWL18Od3d3mW1r165VRGlye/vtt+Hh4QEPDw+UlpZiwYIFiI6OxtSpU2FlZQUTExPExsaiVatWaN++vcLf/9ChQygoKMBnn32G//znPxW279ixAwcOHMDnn38uXWlvbGyMO3fuICsrq8KCPgC4c+cOdHV1pbMirzN16lSEhoZi7dq1GDVqlPSbvJOTEzZs2IDQ0FCMGTOm0r/1rVu3kJCQgC5duqB169YAgIEDB6JRo0aIiorCRx99VOnVAUQkH14iSHVa+Xn9nJycau1X/o365W/rJ0+erNZ6gDfx9OlTPH36tEJd5Zf0lY+pfMHc999/j9LS0grHyczMlHmto6NT5amIF4WGhqJJkyaYPn06hg0bVuFnzJgxyM7Olrmr35AhQwAAISEhFY4XGxuLmzdvYtCgQVBTq/o/IZqamvjwww+RnZ2N7du3S9stLCwwcuRIXLx4EQEBARX2y87Oxueffw4AmDdvnrTd0NAQ06dPR2pqKhYuXIiioqIK+z558gTLly+vsjYiseNMANVplpaWUFNTQ2BgIHJycqCjowNjY2NYW1u/dr/u3bvDyMgI/v7+SE1NRcuWLZGYmIioqCiYmZnhn3/+UXrtSUlJmDhxIoYOHYqOHTvCwMAAd+7cwZ49e2BsbCy9OsHKygqzZ89GQEAAXF1d4ezsjBYtWiAtLQ1Xr17Fn3/+ib///lt63G7duuHUqVPYtGkTWrVqBYlEghEjRlRaw+3bt5GQkIDRo0dXOCdfztHREZqamggNDZWuSxg9ejQOHjyIjRs34tq1a7Czs4O2tjauXbuGyMhIGBkZyXwwV8XFxQUbNmxASEgIJk2aJD2//8033yAjIwMbNmxAXFwchg4disaNGyMpKQnh4eF4/PgxFi1aBHt7e5njzZ49G+np6di/fz/Onz+PESNGoG3btiguLsb169dx+PBhaGpq8l4BRFVgCKA6rVWrVli+fDmCgoLwzTffoLi4GG5ublWGAAMDA2zevBn/+9//sHPnTpSUlKBr164ICgpCaGhorYSAli1bwt3dHWfOnMGxY8dQVFSEFi1awMPDAzNmzJC5ouGTTz5B165dsWPHDmzfvh0FBQUwNDREx44dK3yQLV68GEuXLkVgYKB0cdyrQkD5zYBed+68cePG6N27N+Li4vDw4UO888470NTURHBwMEJCQvDLL78gICAAJSUlaNmyJd5//314e3vj7bfflvt3oaGhgZkzZ2Lx4sXYtm0bPvnkEwDPF/kFBwcjMjISkZGR2Lhxo3Ts9vb2mDZtWqWXAqqpqcHPzw/Dhw/H3r17ERUVhaysLGhqaqJdu3aYMGGCQq94IGqoJIIybh5OREREdR7XBBAREYkUQwAREZFIMQQQERGJFEMAERGRSDEEEBERiRRDABERkUgxBBDVATk5ObC0tIS5uTmioqIq7ePo6Ahzc3Ppj4WFBezt7TFx4kT88ssvrz3+rl27YG5uju7du1e4iyERiRdvFkRUBxw8eBDFxcUwNjZGaGgoXFxcKu3XsmVLfPbZZwCAsrIyPHr0CBEREfjss8+Qnp5e4TG+5cLCwtC2bVvcv38fhw8fhpubm7KGQkT1CGcCiOqA0NBQ9O7dG5MnT8a5c+dw//79Svvp6+vDxcUFLi4ucHNzg7e3N/bu3QtNTU2Eh4dXus/169dx9epVfPzxx+jcuTPCwsKUORQiqkcYAohU7OrVq0hMTISbmxtGjRoFDQ2Nan1QN27cGFpaWtDU1Kx0e2hoKHR0dODk5AQ3NzecO3cOSUlJCqqeiOozhgAiFXvxQ7pp06YYOHAgIiMjUVZWVqFvaWkpsrKykJWVhYyMDCQmJmLhwoXIz8/HuHHjKvQvKirCwYMHMWzYMOjo6GDkyJGvnTUgInHhmgAiFSosLMQvv/wCZ2dn6WOTXV1dcfToUcTGxmLAgAEy/e/cuQM7OzuZtkaNGmHp0qXw8PCocPxjx44hOzsbrq6uAIBmzZphwIABiIiIwJw5c6SPXCYicWIIIFKhX3/9FTk5OdIPaQAYOHAgDA0NERYWViEEtG7dGn5+fgAAQRCQlpaG8PBwLFmyBBoaGnB3d5fpHxoaitatW6NXr17SNjc3Nxw7dgyxsbEYOHCg8gZHRHUeQwCRCoWGhqJZs2Zo2bIl7t27J23v27cvDh8+jKysLDRr1kzarqOjg759+8ocY9SoUXBzc4Ofnx8cHR3RtGlTAEBqaipOnToFDw8PmYWG7dq1g56eHkJDQxkCiESOIYBIRZKTk3HmzBkIggBnZ+dK+xw4cOCVl/2V09DQQJ8+fbB9+3ZcvnxZOnsQHh6OsrIy7Nu3D/v27auw34kTJ5CZmQlDQ8M3HgsR1U8MAUQqEh4eDkEQ4OfnB319/Qrb165di7CwsCpDAACUlJQAAPLz8wE8P1UQERGBTp06wdvbu0L/jIwMfPvtt4iKisK0adPebCBEVG8xBBCpQFlZGSIiImBmZlbpgj4AuHXrFgICAnD58mVYWVm98liFhYWIjY0FAHTu3BkAEBcXh9TUVEyaNAnDhg2rdL+tW7ciLCyMIYBIxHiJIJEKnDx5Eg8fPnzlaQAAcHJyAvB83UC5vLw8REVFISoqCpGRkdi0aRPGjBmD5ORkjB07FiYmJjL7DB069LXHv3XrFi5evKiAERFRfcSZACIVkOdD2szMDCYmJoiJicFXX30FAPj333+xYMECaZ+33noL7du3x+LFi6X3CcjOzsaxY8fQpUsXGBsbv/L4Tk5OCA4ORlhYGLp166aIYRFRPSMRBEFQdRFERERU+3g6gIiISKQYAoiIiESKIYCIiEikGAKIiIhEiiGAiIhIpBgCiIiIRIohgIiISKQYAoiIiESKIYCIiEikGAKIiIhE6v8BQ3C/CSWX8UAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(ABAtrain)\n",
    "plt.xlabel(\"train set AUROC\\nABA\")\n",
    "plt.ylabel(\"brain areas\")\n",
    "#plt.ylim([-2,175])\n",
    "plt.text(0.95, 100, \"mean = %.3f\" %np.mean(ABAtrain),size=20)\n",
    "plt.savefig(\"051321_1.pdf\", bbox_inches='tight')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaUAAAE3CAYAAAAKfF7mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3dd1QUV/8G8GdpEpoFiQ0NglnQKMWO2DViLAEUe8QWDTHx1dhjDyGxJL5JxBY1gBpsARFE1EAswYZg7CIRFQEbC9KVIszvD37M6wrqUpZd4fmc4wl7Z2bne+GEhztz945EEAQBREREakBD1QUQERGVYCgREZHaYCgREZHaYCgREZHaYCgREZHaYCgREZHaYCgREZHaeCtCKSoqCu7u7ujRowcsLS2xf/9+ue2CIMDLywvdu3eHtbU1xo8fj1u3bsntk5GRgXnz5qFDhw7o0KED5s2bh8zMzOrsBhERvcFbEUpPnz6FVCrF4sWLoaurW2r71q1b4e3tjaVLl8Lf3x8NGjTApEmTkJ2dLe4zZ84c3LhxA1u3bsW2bdtw48YNzJ8/vzq7QUREbyK8ZWxtbYWAgADxdVFRkeDg4CBs3LhRbHv27Jlga2sr7N69WxAEQYiLixOkUqkQHR0t7hMVFSVIpVLh9u3b1Vc8ERG91lsxUnqdpKQkyGQyODg4iG26urro1KkTLl68CAC4ePEi9PT00L59e3GfDh06QE9PT9yHiIhU760PJZlMBgBo2LChXLuxsTFSUlIAACkpKWjQoAEkEom4XSKRoEGDBuI+r2JpaQlLS8sqrpqIiMry1odSiRcDR9HtgiC88TgiIqo+b30omZiYAPjfiKlEamqqOHpq2LAhUlNTIbywILogCEhLS4OxsXH1FUtERK/11oeSqakpTExMcObMGbEtLy8P0dHRsLOzAwDY2dnh6dOncvePLl68iKdPn4r7EBGR6mmpugBF5OTkICEhAQBQVFSEBw8eICYmBnXr1kXTpk3h5uaGzZs3w9zcHGZmZti0aRP09PQwZMgQAICFhQV69OiB5cuX49tvv4UgCFi+fDn69OkDc3NzVXaNiIheIBEE9X/IX2RkJNzc3Eq1u7i4YNWqVRAEAevXr8fevXuRkZEBGxsbLFu2DFKpVNw3PT0dnp6eOHbsGACgb9++WLZsGYyMjF577pJJDrGxsVXYIyIiKstbEUqqxFAiIqo+b/09JSIiqjkYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDbeimWGiIhqo+PHj8Pb2xs3btxAUVERWrVqhbFjx8LFxaXc7yUIAg4cOID9+/cjNjYWubm5aNiwIdq1a4dZs2ahZcuWcvsnJiZi06ZNOH36NFJTU1G3bl106dIFX3zxBSwsLEq9/8KFCxEYGPjK84eGhpZ53MsYSkREauj333/Ht99+i3r16uHjjz+GtrY2jh49ioULF+Lff//FggULFH6vvLw8zJw5E8ePH0fLli0xZMgQ6OvrIzk5GdHR0YiPj5cLpevXr8PNzQ3Z2dno2rUrBg0ahEePHuHPP//E8ePH4ePjA1tb2zLP5ebmVubybfXr11esWFU98vZtIZVKBalUquoyiKgWSUxMFNq2bSt07txZSExMFNvT09OF/v37C1KpVPjnn38Ufr8VK1YIUqlUWLt2rVBYWFhqe35+vtxrJycnQSqVCj4+PnLt//zzj9CmTRthwIABpY5ZsGCBIJVK5eqtCI6UlMxs4SGVnDd+1WCVnJdqj6SkJPTr1w8uLi6YPn06fvzxR0RGRiI/Px92dnZYuHAhpFIpnjx5gp9++gnHjx9HRkYGpFIp5s2bh65du5Z6z+fPn2Pv3r0ICgpCXFwcCgsL0bJlS7i6umLs2LHQ0JC/Db5//34cP34cN27cgEwmg5aWFqRSKcaMGQMnJ6dS7z9+/HicP38e169fx7Zt27B//348ePAAxsbGGDJkCGbOnAkdHR2lfc8UFRAQgPz8fEydOhWmpqZie926dfHZZ59h8eLF2LNnj0KP3klISMCePXvQrl07fPXVV2U+2FRbW1v8OjExETExMTA2Ni61ELadnR369euHo0ePIiIiAn379q1EL8vGUCKiSrl//z5GjBgBCwsLuLi44P79+wgLC8P48eOxd+9efPrppzAwMMBHH32EjIwMhIaGYurUqTh69CiaNm0qvk9BQQHc3d1x6tQp8RJTnTp1EBkZiW+//RaXL1/GDz/8IHfuFStWoFWrVujUqRNMTEyQnp6OkydPYv78+bh79y5mzZpVZs1z5szBhQsX0KNHD/Tq1Qt///03tm3bhidPnmDlypVK/X4p4ty5cwCAHj16lNrWs2dPuX3eJCQkBEVFRXBxcUF2djaOHTuGR48eoV69eujatSvee+89uf1LHpjarFmzUn8EABBD8uzZs2WG0t9//43s7GxoamrivffeQ9euXWFgYKBQrQBDiYgq6fz585g1axY+//xzsW3Dhg1Yt24dRo4ciYEDB2LFihXiL7hu3bphwYIF8PX1xaJFi8RjNm/ejFOnTuGTTz7BokWLoKmpCQAoLCzE0qVLERAQAEdHR/Tv3188JiQkBC1atJCrp2SEsXXrVowZMwaNGjUqVXNiYiJCQkJQr149AMBXX30FJycnHDhwALNnzxafaP06MTExCA8PL8d3CpgwYcIbH5cDAHfv3gUAmJmZldr27rvvQk9PD48ePcKzZ8/wzjvvvPa9rl69CgDIyspC//79kZ6eLm6TSCQYM2YMlixZIn6/S+79PHjwAIIglBpZJSUlAQDu3LlT5vm++eYbudf6+vqYM2cOxo0b99o6SzCUiKhSmjVrhmnTpsm1ubi4YN26dcjPz8f8+fPl/uIeOnQoFi9ejJiYGLGtqKgIv//+O0xMTPD111+LvyABQFNTEwsXLsT+/ftx8OBBuVB6OZAAQEdHB+PGjcO5c+dw9uxZODs7l9pn7ty5YiABgJ6eHoYOHYoNGzbg2rVr6NOnzxv7HRMTg/Xr179xvxe5uLgoFErZ2dkAAENDwzK3GxgY4OnTp8jKynpjKD158gQAsG7dOtjb22PBggVo1qwZrly5guXLl2PXrl1o0KABZsyYAQBo2bIlzMzMEB8fj507d8pdwrt8+TL++usvAEBmZqbceTp27IiePXvC1tYWxsbGSE5ORlhYGNavXw8PDw9oaWlh1KhRb+w7Q4mIKqV169ZyIQIU/zUPFP+l//KlG01NTRgbG+Px48di2927d5Geni4+Obosurq6pf46f/DgAbZu3YqzZ8/i4cOHyM3Nldv+4jle1LZt21JtTZo0AQBkZGSUeczLhg0bhmHDhim0b1UT/v8xeGXdH3pZYWEhAMDExAQbNmyArq4uAMDe3h7r1q2Di4sLfHx88Nlnn4n30zw8PPDpp5/iu+++w/Hjx2FlZYXHjx/jzz//hIWFBWJjY0v9zF1dXeVeN2/eHJMnT0bLli3h7u6On376Ca6urqWOexlDiYgqpay/5rW0tF65rWT78+fPxdcll5Ti4+NfO/rIyckRv05MTISrqysyMzPRsWNHdO/eHQYGBtDU1MT9+/cRGBiI/Pz8Mt+nrNFKyS/LoqKiV56/uhgYGCAtLQ1ZWVllTqUu+T4ocq+mbt26AIrvT5UEUgkrKyuYmpoiISEBd+7cgZWVFQCgS5cu+OOPP7Bp0yZERUUhKioKjRs3xueffw4rKytMnz4dDRo0UKgvffr0QaNGjfD48WPExcWJD059FYYSEalcSXh9+OGHCl8S8/HxQXp6OlauXFlqxBISEvLaD3JWBWXeU2rZsiXS0tIQHx9fKpSSk5Px9OlTNG7c+I2X7kre69SpU6/8A6GknpdHmVZWVvjll19K7b9u3ToAQLt27d547hINGjTA48eP8ezZszfuy1AiIpUzNzeHkZERLl26hIKCArkpyq9y7949AMCAAQNKbTt//nyV1/gyZd5T6tq1K/755x9ERESUmvb9999/i/soomvXrti5cydu3bpValt+fr74fWzWrNkb3ys/Px8HDhyAhoYGBg9W7GMnWVlZuHPnDiQSiULn4Np3RKRyWlpa+OSTTyCTyeDp6Vnqr3ageIQQFxcnvi75BfdyAEVERMDf31+5BaP4nlJsbGy5/r34maM3vbeOjg78/PzE2W5A8f2uX3/9FQAwevRouWOysrJw+/ZtJCcny7X37NkTzZs3x6lTp3D69Gm5bRs2bEBWVhY6d+4sN+Pw6dOn4r2oEgUFBVi+fDnu37+PMWPGyE0ykclkYri9KCcnBwsXLkReXh66deum0KxGjpSISC1Mnz4dN2/exJ49e3D8+HF07doVjRo1QmpqKu7du4d//vkHX331FVq1agUAGDt2LPbv34+ZM2diwIABaNSoEW7duoWIiAh89NFHCA0NVXGPKq558+aYP38+PD09MXz4cAwaNEhcZujRo0eYPHlyqRFUWFgYvv76a7i4uGDVqlViu46ODlatWoUpU6Zg6tSp6N+/P5o1a4arV68iKioKDRo0gIeHh9x7RUZGYsmSJbC3t0eTJk2QnZ2NkydP4v79++jdu3epJY7u3LkDNzc32NnZwdzcXJzIcubMGchkMjRv3hyenp4K9Z2hRERqQVtbGxs3bkRQUBACAwNx4sQJPH36FPXr14epqSlmzpyJoUOHivtbWVlhx44d+Pnnn/H333/j+fPnsLKywvr162FoaPhWhxJQvPpEs2bN4O3tjQMHDkAQBFhYWGDWrFnlXpC1Y8eOCAgIwPr16xEZGYljx47B2NgYo0aNwvTp09G4cWO5/c3MzNC+fXtERUUhNTUVurq6sLKywpdffglnZ+dSH6pt0aIFRo0ahatXr+LYsWPIysqCrq4uWrZsiXHjxmH8+PEKf4BWIpTMLaQylcwUiY2NrdDxXGaIiEhxvKdERERqg6FERERqg6FERERqg6FERERqg6FERERqg6FERERqg6FERERqg6FERERqg6FERERqg6FERERqo0asfVdYWAgvLy8EBwdDJpPBxMQEQ4cOxYwZM8SHjQmCgPXr12Pv3r3IzMyEjY0Nli1bhvfff1/F1VNNwOWkiKpGjRgpbd26Fbt27cKSJUtw+PBhLF68GLt27RKXeC/Zx9vbG0uXLoW/vz8aNGiASZMmITs7W4WVExHRi2pEKF28eBF9+vRB3759YWpqin79+qFv3764cuUKgOJR0o4dOzBt2jQ4OjpCKpVi9erVyMnJQUhIiIqrJyKiEjUilDp06IDIyEjcvn0bABAXF4dz586hZ8+eAICkpCTIZDI4ODiIx+jq6qJTp064ePGiSmomIqLSasQ9palTpyInJweDBw+GpqYmnj9/Dnd3d4wbNw5A8VMRAaBhw4ZyxxkbG5d6SiNVHu+vEFFF1YhQCg0NxYEDB7B27Vq0atUKMTEx+P7772FqaooRI0aI+0kkEhVWSUREb1IjQmnNmjWYPHkyBg8u/kvZ0tISDx48wJYtWzBixAjxufAymQxNmjQRj0tNTS01eiIiItWpEfeUcnNzoampKdemqamJoqIiAICpqSlMTExw5swZcXteXh6io6NLPeeeiIhUp0aMlPr06YMtW7bA1NRUvHzn4+MDZ2dnAMWX7dzc3LB582aYm5vDzMwMmzZtgp6eHoYMGaLi6omIqESNCKUlS5bgl19+wTfffIPU1FSYmJhg5MiR+OKLL8R9pk6diry8PHh4eCAjIwM2Njbw9vaGgYGBCisnIqIX1YhQMjAwwOLFi7F48eJX7iORSDBjxgzMmDGjGisjIqLyqBH3lIiIqGZgKBERkdpgKBERkdpgKBERkdqoERMdiGorLulENQ1HSkREpDYYSkREpDYYSkREpDYYSkREpDYUnuiQlpaGJ0+ewMLCQmxLTEyEr68v0tPT4ezsjB49eiilSCIiqh0UDqXvvvsO8fHx8Pf3BwDk5ORg3Lhx4kPyDh8+jO3bt6NTp07KqZSIiGo8hS/fXbp0SXy8OFD8YL3k5GRs2bIFERERsLCwwLZt25RSJBER1Q4Kh1JqaqrcA/IiIiLQtm1b9OzZEyYmJnBxccGNGzeUUiQREdUOCoeSlpYW8vLyxNfnz5+Xu1RnaGiI9PT0qq2OiIhqFYVDyczMDEePHoUgCPjrr7+QkZEBe3t7cfujR49Qt25dpRRJRES1g8ITHcaNG4eFCxeiU6dOyM3NRfPmzeVCKSoqCpaWlkopkoiIageFQ8nZ2RkSiQTh4eEwMDCAu7s7tLW1ARRPF8/OzsaYMWOUVigREdV85VqQ1cnJCU5OTqXa69evj/3791dZUUREVDtxRQciIlIb5RopPX/+HOHh4bh8+TIyMzNRVFQkt10ikeD777+v0gKJiKj2UDiU0tPT4ebmhlu3bkEQBEgkEgiCAADi1wwlIiKqDIUv3/3888+4c+cOPD09ERYWBkEQ8NtvvyE0NBSDBw9Gu3btEBkZqcxaiYiohlM4lE6ePAlnZ2cMHz4cBgYGxQdraMDc3Bw//vgjdHV1sXbtWqUVSkRENZ/CoSSTydCuXTsAxas7AEB+fr64vV+/fjh27FgVl0dERLWJwqFUr149PHv2DACgr68PLS0tPHz4UNyura2NzMzMqq+QiIhqjXItMxQXF1d8kIYG2rRpg8DAQOTn5+PZs2c4cOAAmjdvrrRCiYio5lM4lBwcHHD06FHxkt3EiRNx+fJldO7cGd26dcO1a9cwYcIEpRVKREQ1n8JTwt3d3TFlyhTo6OgAAAYNGgQtLS0EBwdDQ0MDAwcOxKBBg5RWKBER1XwKh5JEIhEDqcSAAQMwYMCAKi+KiIhqpwotM5Sfn4/Hjx/Lzb4jIiKqrHKF0vXr1+Hm5ob27dujd+/euHDhAoDip9JOmDABZ86cUUqRRERUOygcSjExMRg3bhwSExNLrRRubGyMvLw8BAYGVnmBRERUeygcSr/88gveffddhISEYM6cOeK6dyW6du2KK1euVHmBikpOTsaCBQvQtWtXtGvXDoMGDcL58+fF7YIgwMvLC927d4e1tTXGjx+PW7duqaxeIiIqTeFQunDhAkaMGAF9fX1IJJJS25s2bYrk5OQqLU5RmZmZGDNmDARBwJYtWxAaGoqlS5fC2NhY3Gfr1q3w9vbG0qVL4e/vjwYNGmDSpEnIzs5WSc1ERFSawrPv8vLyYGho+Mrtqvzlvm3bNpiYmGDNmjVi24sf5BUEATt27MC0adPg6OgIAFi9ejXs7e0REhKC0aNHV3vNRERUmsIjpRYtWuD69euv3H7u3Dm0atWqSooqr/DwcNjY2GDWrFmwt7eHk5MTfv/9d/ESY1JSEmQyGRwcHMRjdHV10alTJ1y8eFElNRMRUWkKh9KQIUMQFBQkN8Ou5DKet7c3IiIiynxUenVITEzErl270Lx5c/z2229wc3PD2rVr4efnB6B4MVkAaNiwodxxxsbGSElJqfZ6iYiobApfvps8eTJOnz6NKVOmwNzcHBKJBCtXrsSTJ0+QkpKCbt26YezYscqs9ZUEQUDbtm0xZ84cAECbNm1w7949+Pn54ZNPPhH3K+teGBERqQ+FR0o6Ojrw8fHBggULUKdOHdSpUwfx8fGoX78+5s2bh19//RUaGhX6LG6lmZiYwMLCQq7N3NxcXMXcxMQEwP9GTCVSU1NLjZ6IiEh1FB4pAcXPUZo4cSImTpyopHIqpn379rh7965cW3x8PJo2bQoAMDU1hYmJCc6cOQNra2sAxRM3oqOjMX/+/Gqvl4iIyqbQ0CYnJwf9+/eHr6+vksupmAkTJuDy5cvYtGkT7t27h8OHD2Pnzp0YN24cgOLLdm5ubtiyZQv+/PNP/Pvvv1i4cCH09PQwZMgQFVdPREQlFBop6evrIz09Hfr6+squp0Ksra2xYcMG/Pe//8XGjRvRtGlTzJw5U+4e19SpU5GXlwcPDw9kZGTAxsYG3t7e4qPdiYhI9RS+fGdjY4OrV69ixIgRyqynwnr37o3evXu/crtEIsGMGTMwY8aM6iuKiIjKReGZCXPnzsWRI0cQEBBQaokhIiKiqqDwSGnlypUwMjLCkiVL8MMPP6BFixbQ1dWV20cikWD79u1VXiQREdUOCodSUlISAKBJkyYAwA+dEhFRlVM4lI4dO6bMOoiIiCr25FkiIiJlYCgREZHaKNeKDgkJCfD19cXly5eRmZmJoqIiue0SiQTh4eFVWiAREdUeCodSbGwsxo4di/z8fLRs2RKJiYl4//33kZaWhpSUFLRo0QKNGjVSZq1EpCbMFh5S2bnjVw1W2blJ+RS+fLdu3Tpoa2sjKChIXG5o0aJFOHXqFDw8PJCZmYnly5crq04iIqoFFB4pXbhwAaNGjYK5uTnS0tLkto0cORLR0dH48ccfsXnz5iovkspPlX/JEhFVlMIjpZycHPER49ra2gCAp0+fitvbt2+Pf/75p4rLIyKi2kThkVLDhg3FD8waGBjgnXfeQXx8vLg9MzMThYWFVV4gkaI4OiR6+ykcSlZWVrh69ar4unPnztixYwesra1RVFSE33//HVZWVkopkoiIageFL98NHToU6enpyM3NBQDMnDkTWVlZcHNzw8SJE5GVlYWvvvpKaYUSEVHNJxEqseT3w4cPERYWBk1NTfTs2VO851STWFpaAiieEl8RvKREVLU4JbxmK9eHZ1/WpEkTuLm5VVUtRERUy3GZISIiUhsMJSIiUhsMJSIiUhsMJSIiUhsMJSIiUhsMJSIiUhsVmhL+7NkzpKeno6yPODVt2rTSRRERUe2kcCgVFRVh27Zt2Llzp7gGXlliYmKqpDAiIqp9FA6lH3/8Ed7e3nj//ffh6OiIevXqKbMuIiKqhRQOpeDgYPTo0QNbt25VZj1ERFSLKTzRITMzE/369VNmLUREVMspHEpSqRQymUyZtRARUS2ncCh9+eWX2LNnDx4+fKjMeoiIqBZT+J7StWvX0LRpUwwaNAgffvghTE1NoaEhn2kSiQRffPFFlRdJRES1g8KhtH79evHr4ODgMvdhKBERUWUoHEp//fWXMusgIlKIqh6cyYcLVg+FQ6lZs2bKrIOIiKhmrn23efNmWFpawsPDQ2wTBAFeXl7o3r07rK2tMX78eNy6dUuFVRIR0cteOVJav349JBIJPv/8c2hoaMjdU3oVdbindOnSJezbtw+WlpZy7Vu3boW3tzdWrVqFli1bYsOGDZg0aRKOHDkCAwMDFVVLREQvemMoTZ06FTo6Om9FKGVlZWHu3Ln47rvvsHHjRrFdEATs2LED06ZNg6OjIwBg9erVsLe3R0hICEaPHq2qkomI6AWvDKWSiQ06Ojpyr9XZ0qVL4ejoCHt7e7lQSkpKgkwmg4ODg9imq6uLTp064eLFiwwlIiI18cpQenlig7pPdNi3bx8SEhKwZs2aUttKVqJo2LChXLuxsTGSk5OrpT4iInqzCj1PSd3cuXMH//3vf+Hn5yeO7MoikUiqsSoiIiqvcofS1atXceXKFWRkZKCoqEhum6ruKV26dAlpaWkYOnSo2FZYWIioqCjs2bMHISEhAIpHTE2aNBH3SU1NLTV6IiIi1VE4lHJzc/Hll1/i9OnTEAQBEolEfPJsydeqCqX+/fujbdu2cm1ff/01zMzM8Nlnn6Fly5YwMTHBmTNnYG1tDQDIy8tDdHQ05s+fX+31EhFR2RQOpQ0bNuD06dNwd3eHvb093NzcsGrVKhgbG2PLli3Izc3F6tWrlVnrKxkZGcHIyEiuTU9PD3Xr1oVUKgUAuLm5YfPmzTA3N4eZmRk2bdoEPT09DBkyRBUlExFRGRQOpaNHj2LgwIGYOXMm0tLSAACNGjWCvb097O3t4erqisDAQMyZM0dpxVbG1KlTkZeXBw8PD2RkZMDGxgbe3t78jBIRkRpROJQePnyIiRMnAgA0NTUBAAUFBcVvoqWFwYMHY/fu3WoTSjt37pR7LZFIMGPGDMyYMUNFFRER0ZsovMyQvr4+CgsLxa81NDTkplMbGhoiJSWl6iskIqJaQ+FQatGiBeLj4wEUj5RatWqFo0ePAiheMSEsLExuZhsREVF5KRxK9vb2OHr0qDhaGjVqFCIiItC/f38MGDAAZ86cwfDhw5VWKBER1XwK31OaNm0anJycxGng48aNQ35+PoKDg6GhoYGvvvoKU6dOVVqhRERU80mEkpShMpWsNh4bG1uh41X1QDIiqlp8yF/1UOjyXU5ODvr37w9fX18ll0NERLWZQqGkr6+P9PR06OvrK7seIiKqxRSe6GBjY4OrV68qsxYiIqrlFA6luXPn4siRIwgICABvQxERkTIoPPtu5cqVMDIywpIlS/DDDz+gRYsW0NXVldtHIpFg+/btVV4kERHVDgqHUlJSEgCIH5Dl6g1ERFTVFA6lY8eOKbMOIiIixe8pERERKVu5nzybn5+PyMhIJCYmAgCaN2+Ozp07o06dOlVeHBER1S7lCqUDBw5g5cqVyMzMlHvqrJGRERYsWIBhw4YppUgiIqodFA6l0NBQLFy4EE2bNsWUKVNgYWEBQRBw+/Zt7NmzB4sXL4auri4GDRqkzHqJiKgGU3jtu48//hjPnz/Hvn37Sj2tNSsrCyNGjIC2tjYOHjyolEJVhWvfERHAte+qi8ITHe7evYthw4aV+fhwQ0NDDBs2DPfu3avS4oiIqHZROJRMTExeu5KDhoYGGjZsWCVFERFR7aRwKLm4uCAwMBA5OTmltmVnZyMgIIATHYiIqFJeOdEhKipK7nXHjh1x/PhxDB06FGPHjoW5uTkkEgni4uKwe/du1K9fHx06dFB6wUREVHO9cqKDlZUVJBKJXNuLu5Zse7ktJiZGGXWqDCc6EBHAiQ7V5ZUjpZUrV1ZnHURERK8OJRcXl+qsg4iIiGvfERGR+mAoERGR2mAoERGR2mAoERGR2mAoERGR2mAoERGR2mAoERGR2mAoERGR2ij349DV0a+//oo///wTd+/ehY6ODmxtbTF79mxIpVJxH0EQsH79euzduxeZmZmwsbHBsmXL8P7776uwciJ6W6hyybDatMRRjRgpnT9/HmPHjsWePXuwfft2aGpqYtKkSUhPTxf32bp1K7y9vXaZ/G4AABsxSURBVLF06VL4+/ujQYMGmDRpErKzs1VYORERvahGhNJvv/2G4cOHQyqVwtLSEmvWrMGTJ0/wzz//ACgeJe3YsQPTpk2Do6MjpFIpVq9ejZycHISEhKi4eiIiKlEjQullOTk5KCoqgpGREQAgKSkJMpkMDg4O4j66urro1KkTLl68qKoyiYjoJTUylL777ju0bt0adnZ2AACZTAYApZ6Ma2xsjJSUlGqvj4iIylYjJjq8aOXKlbhw4QJ2794NTU1NuW0vPx+KiIjUS40aKX3//fc4dOgQtm/fjubNm4vtJiYmAP43YiqRmppaavRERESqU2NCydPTEyEhIdi+fTssLCzktpmamsLExARnzpwR2/Ly8hAdHS1e4iMiItWrEZfvvvnmGwQFBWHDhg0wMjISR0R6enrQ19eHRCKBm5sbNm/eDHNzc5iZmWHTpk3Q09PDkCFDVFw9ERGVqBGhtGvXLgDAxIkT5dq//PJLzJgxAwAwdepU5OXlwcPDAxkZGbCxsYG3tzcMDAyqu1wiInoFiSAIgqqLUGeWlpYAgNjY2Aodr8pPgRNRzcAVHYiIiFSAoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGqDoURERGpDS9UFEBHR65ktPKSS88avGlzt5+RIiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1AZDiYiI1EatCyU/Pz/07dsX7dq1w7BhwxAdHa3qkoiI6P/VqlAKDQ3F999/D3d3dxw4cAB2dnaYOnUqHjx4oOrSiIgItSyUfHx84OLigpEjR8LCwgJLly6FiYkJdu/ererSiIgItWiV8Pz8fFy/fh2TJ0+Wa3dwcMDFixffeLylpWWFzlunQkcREameZeDsSh0fGxtb7mNqzUgpLS0NhYWFaNiwoVy7sbExZDKZiqoiIqIX1ZqRUgmJRKJQW4mKJD0REVVMrRkp1a9fH5qamqVGRampqaVGT0REpBq1JpR0dHTwwQcf4MyZM3LtZ86cgZ2dnYqqIiKiF9Wqy3eTJk3C/PnzYW1tjfbt22P37t1ITk7G6NGjVV0aERGhloXSoEGDkJaWhk2bNiE5ORlSqRRbtmxBs2bNVF0aEREBkAiCIKi6CCIiIqAW3VMiIiL1x1CqpPKspRcZGQlLS8tS/27fvl2NFVdeedcPzM/Pxy+//IK+ffuibdu26N27N3bs2FFN1VZeefq7cOHCMn/Gtra21Vhx5ZX3Z3zw4EE4OTnBxsYGDg4OmDt37lv3+b/y9tnPzw8fffQRrK2t4ejoiAMHDlRTpVUjKioK7u7u6NGjBywtLbF///43HhMbG4tPPvkE1tbW6NGjB9avX48qv9gmUIUdOnRIaNOmjbB3714hLi5O8PDwEGxtbYX79++Xuf+5c+cEqVQq3Lp1S0hOThb/PX/+vJorr7jy9lkQBOHLL78Uhg8fLpw6dUpITEwULl26JJw7d64aq6648vY3MzNT7mebnJws9OvXT1i4cGE1V15x5e1zdHS0YGVlJfj4+AgJCQnCxYsXBWdnZ8HNza2aK6+48vbZz89PsLW1FQ4ePCgkJCQIISEhgq2trfDXX39Vc+UVd+LECWHt2rXC4cOHBWtrayEgIOC1+2dlZQndunUT/vOf/wixsbHCkSNHBFtbW+G3336r0roYSpXg6uoqLF68WK7tww8/FH788ccy9y8JpdTU1OooTynK2+eIiAihffv2b22fy9vfl0VHRwtSqVS4cOGCMspTivL2edu2bULv3r3l2vz9/QVbW1ul1VjVytvnUaNGCd99951c28qVK4XRo0crrUZlsrW1fWMo+fn5CXZ2dsKzZ8/Etg0bNgjdu3cXioqKqqwWXr6roJK19BwcHOTaFVlLz9XVFd27d8eECRNw7tw5ZZZZpSrS5/DwcLRr1w6+vr7o2bMnBgwYAE9PT+Tk5FRHyZVSmZ9xiT/++APvv/8+2rdvr4wSq1xF+ty+fXvIZDIcO3YMgiDgyZMnCA0NRc+ePauj5EqrSJ/z8/NRp478ypZ16tTB1atXUVBQoLRaVenSpUvo2LEjdHV1xbbu3bsjOTkZSUlJVXYehlIFVWQtPRMTE6xYsQLr1q2Dl5cXWrZsiYkTJyIqKqo6Sq60ivQ5MTERFy5cwM2bN+Hl5YWlS5ciIiICX3/9dXWUXCmVXS8xKysLR44cwYgRI5RVYpWrSJ/t7Oywdu1azJ07F23btoW9vT0EQcDq1auro+RKq0ifu3fvjoCAAFy5cgWCIODq1avw9/dHQUEB0tLSqqPsapeSkgJjY2O5tpLvWUpKSpWdp1Z9TkkZyrOWnrm5OczNzcXXdnZ2uH//Pn777Td06tRJaTVWtfL0WRAESCQSrF27FoaGhgCApUuXYsqUKUhJSXkrlngq73qJJYKDg1FYWAgnJydllKVU5elzXFwcPD09MX36dHTv3h0ymQxr1qzBsmXLsGbNGmWXWmXK0+fp06dDJpNhzJgxEAQBxsbGcHZ2xrZt26CpqansUlXm5e+H8P+THBT5/0FRHClVUFWtpWdjY4N79+5VdXlKUZE+m5iYoFGjRmIgAYCFhQUAqP3DFSv7M963bx8GDBiAevXqKavEKleRPv/666+wtrbGp59+CisrK/To0QPLly9HUFAQHj58WB1lV0pF+qyrq4uVK1fi0qVLOHbsGE6cOIFmzZpBX18f9evXr46yq13Dhg1LjYhSU1MBoNQIqjIYShVUVWvpxcTEwMTEpKrLU4qK9Ll9+/ZITk6Wu4cUHx8PAGq/kkZlfsZXrlzBzZs3MXLkSGWWWOUq0ufc3NxSo4OS18Jb8Nn8yvyctbW10bhxY2hqaiI0NBR9+vSBhkbN/LVqa2uL6Oho5OXliW1nzpzBu+++C1NT0yo7j+aKFStWVNm71TIGBgbw8vKCiYkJdHV1sXHjRkRHR+P777+HkZER5s+fj7CwMHz44YcAAF9fX6SkpEBDQwMymQy+vr74448/sHDhQnH0oO7K2+eWLVsiICAA165dQ6tWrRAfHw8PDw906tQJw4cPV3Fv3qy8/S2xbt06PH36FAsWLKjSSxvVobx9zs3NxW+//Yb69eujXr16iIuLw3fffQcTExN8+umnKu6NYsrb57t37+LkyZPQ1dVFQkICVq1ahcuXL+Onn36CkZGRinujmJycHNy+fRspKSn4448/IJVKYWhoiIKCAhgaGmLt2rX49ddf4eLiAgAwMzPD3r17ERMTA3Nzc1y4cAGrV6/GZ599VqUTeXhPqRLetJbey5cuCgoKsHr1ajx+/Bi6urpo1aoVtmzZgl69eqmi/Aopb5/19fXh4+MDT09PuLq6wsjICP3798ecOXNUUX65lbe/AJCdnY3Q0FBMnz79rQskoPx9HjZsGHJycuDn54fVq1fD0NAQXbp0wbx581RRfoWUt89FRUXw9fXF3bt3oaWlhS5dumD37t1VOmJQtmvXrsHNzU187eXlBS8vL7i4uGDVqlWQyWRITEwUtxsaGsLb2xseHh4YPnw46tati8mTJ2PSpElVWhfXviMiIrVRMy9+EhHRW4mhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhREREaoOhRG+VmJgYeHl5ISkpqVrO5+vri/3791fLuV4lPDwcXl5eFT5+zZo1sLS0xIABA/CqB017eXnB0tISkZGRZW6PjIyEpaWlXB1JSUmwtLSU+2dtbY3Bgwdj3bp1ePbs2StrunnzJhYsWIA+ffqgbdu26Ny5M8aPH4+AgAAUFha+8rjHjx9jzZo1GDp0KOzs7NC2bVv07dsXc+fOxdmzZxX8jpA601J1AUTlERMTg/Xr16Nz584wNTVV+vl27NiBZs2aYdiwYUo/16uEh4cjMDAQM2bMKPexz58/R3BwMFq0aIF79+4hKioKnTt3rtL6HBwc4OTkBABIS0vD0aNHsWHDBly6dAne3t6l9t+1axc8PT1hZGQEFxcXmJubIzMzE3/++ScWLVqEQ4cOYf369dDT05M77sSJE5g9ezby8/MxcOBAjBo1CnXq1MH9+/fx119/YeLEidiyZQt69epVpf2j6sVQIqrBTpw4AZlMBl9fX8yZMwcBAQFVHkpmZmZiKAHA+PHjMXLkSJw+fRrXrl1D27ZtxW2nT5+Gh4cHpFIpfH190aBBA3HblClT8NNPP2Hz5s1YsWIF1qxZI267desWZs6cibp168LHxwcWFhZyNcycORPBwcHQ0dGp0r5R9ePlO3preHl54euvvwYAuLm5iZeMFi5cKO6Tn5+PzZs3Y/DgwWjXrh06duwId3d33LhxQ+69BEGAr6+veBmoffv2cHR0xKJFi1BQUAAAsLS0xP3793H+/Hm5S1RvunR44MABuLq6omPHjrC1tUW/fv0wZ84cPHnyRG6/+Ph4zJs3D927dxcvQ61evRpPnz4V9xk/fjwCAwPFekr+KXpJ0d/fH82bN0fXrl0xdOhQHD16FNnZ2QodW1Gamppi8N27d09u29q1a8X/vhhIJWbNmgUbGxsEBQUhNjZWbF+3bh1yc3Ph6elZKpAAQCKRwMnJCfb29lXZFVIBjpTorfHhhx9CJpNh7969cHd3h7m5OQCgRYsWAICCggJMmTIFFy9ehJOTE8aNG4fs7Gzs27cPY8aMwe+//4527doBADZu3Ih169ahT58+GD16NDQ1NZGUlIRjx44hPz8f2traWLNmDVauXIn69evD3d1drKOsX6YlgoKCsGDBAnTs2BH/+c9/oKuriwcPHuDvv/9GamqqeOy1a9cwYcIEGBkZYdSoUWjUqBFu3ryJnTt34uLFi9i5cye0tbXh7u6OoqIiREdHy40c2rdv/8bvl0wmQ0REBD7//HNIJBK4uLjA19cXhw4dwqhRo8r/AyiHxMREAEDdunXl2q5fv4727dvj/fffL/M4iUQCV1dXXL58GWFhYbC0tEReXh5OnDiBJk2aoGfPnkqtm1SPoURvDSsrK9ja2mLv3r3o1q0bunTpIrfdz88P58+fx7Zt29CjRw+xfezYsRgyZAjWrFmDnTt3Aii+T2NhYYHNmzfLvcfcuXPFr52cnPDLL7+gYcOGcpenXicsLAz6+vrYvn07tLT+97/XrFmz5PZbtGgRTExM4O/vDwMDA7Hd3t4eX375JQ4ePIhhw4bBwcEBBw8eRHR0tMI1lAgMDERhYaF4nJWVFVq3bg1/f/8qDaW8vDxxFJiWlobDhw8jLCwMjRs3lrtUeOvWLQDABx988Nr3a9OmDQDg33//BVA8oszPz4eVlVWV1Uzqi5fvqMYIDg6Gubk5PvjgAzx58kT8l5+fj27duuHChQvIzc0FABgYGODx48eIjo6u0hoMDQ2Rm5uLEydOvHKmW2xsLGJjYzFkyBDk5+fL1dqhQwfo6enh9OnTla4lICAAHTt2RPPmzcU2FxcXXLlyRQyIquDv7w97e3vY29tj0KBB8PLyQpcuXeDr6yt3j6fksuGLIVwWQ0NDAEBWVla5jqOagSMlqjFu376N3Nzc195XSEtLQ5MmTTB79mx88cUXGDduHN5991107twZvXv3hqOjY6Vuln/22WeIiorCF198gXr16qFz587o2bMnPvroI/GX6u3btwEU3yN71VTvlJSUCtcAANHR0YiPj8fHH38sd1/HxsYGGhoa8Pf3F+/PlYdEIinV1q9fP3zyyScoLCzEvXv3sG3bNjx69KjU97Gk/2+6p1WyvSScSo7Lyckpd7309mEoUY0hCAKkUulrf9mW3NOxs7NDWFgYTp06hcjISERGRiIkJASbNm3Crl27UK9evQrVYGZmhtDQUJw9exZnz57F+fPnsWTJEqxbtw5+fn7i/S8AmDx5stxlxhcZGRlV6Pwl/P39ARRPEFi3bl2p7cHBwZg7dy60tbUBAHXq1AEAcST5spLPHJXs96LGjRujW7duAIAePXqgZ8+e+PjjjzF79mzs2bNHDLKS+0jXr19/be0lk1KkUimA4u+pjo4Obt68+drjqGZgKNFbpay/1Eu89957SEtLQ9euXaGh8eYr0/r6+nB0dISjoyOA4ntSHh4e8Pf3x6efflrhGnV0dNCrVy/x8zInT57EtGnT4OPjg+XLl+O9994DAGhoaIi/zF/ndX0uS3Z2No4ePQoHBweMHDmy1PbY2Fhs3LgRx44dE/te8pmv27dvl/k5n5LR3YuXAl+lRYsWmDx5MjZs2ICQkBAMHTpUPLZ169a4ePEi4uLi0KpVq1LHCoIgBmr//v0BFAdhr169xD8iunfvrsi3gd5SvKdEb5WSD1RmZGSU2ubs7AyZTAYfH58yj33xktjL07OB/92Af/G99fX1kZ6ernB9Zb1vyY37kvdt06YNpFIp9uzZI85Se9Hz58/lzlnSZ0XrCA0NxdOnTzF69GgMHDiw1L9p06bhnXfeQUBAgHiMg4MD9PT0sG/fvlKX17Kzs/HHH39AT09P4SnXEydOhKGhIdavXy+3QsPs2bMBAPPmzUNaWlqp47y8vHDp0iU4OTnJTWwomcm4ZMkS3Llzp8xzHjx4kKs61ACaK1asWKHqIogUpaenh507d+LevXvQ0NBAfHw8nj17hsaNG6Ndu3a4cuUKdu/ejStXrkAmk+HWrVsICwvDzz//jPDwcHFlhj59+iA6OhoJCQm4d+8eTp48if/+97/Izc3FkiVL0LBhQwDFl5JOnTqF58+fQyaTIS4uDqampuJlr5eNHTsWhw4dQlJSEhISEnDu3DmsWrUKaWlpmDt3LszMzCCRSGBnZ4egoCD4+fkhJSUFDx48QHR0NIKCgrBs2TI0btwYrVu3BlAcZmFhYXjw4AFyc3Nx+/ZtvPPOO3LTrV/0zTffIDMzEx4eHmXWqa2tjRs3buDEiRNwdXWFgYEBdHV1Ua9ePRw4cAAHDx5ERkYGEhISEB4ejuXLlyMpKQlLly5Fhw4dxPfJzMzEjh07YG1tXWp0VadOHeTk5CA8PBwtWrQQA+a9995DvXr1EBgYiICAAKSlpeHhw4c4d+4cVq9ejeDgYDg4OGDNmjVytRsbG6N169YIDg6Gn58f4uPjkZycjH///Rfh4eH4/vvvsWvXLgwePFgcidLbSSK8aooQkZoKDAzE1q1bkZCQgIKCAri4uGDVqlUAikcZu3btQlBQkHjJ6d1330W7du3g4uIiXvrZsmULTp48iTt37iArKwvGxsawsbHBZ599JjdlOTU1FStWrEBkZCQyMzMhCAL++uuvVy5xtG/fPhw+fBj//vsvMjIyUK9ePbRu3RpTpkxB165d5fa9f/8+fv31V5w6dQrJycnQ19dHs2bN4ODggLFjx6JJkyYAgKKiIvzwww84dOgQZDIZioqKsHLlyjKXPoqLi8PgwYMxYMCA166XFxISgjlz5uCrr76S+wzW6dOn4evriytXriArKwuGhoawtrbGpEmTSl1qTEpKQr9+/TBu3DgsW7as1DnS0tLQt29fmJiY4PDhw9DU1BS33bhxAz4+Pjh//jxSU1PxzjvvwMrKCs7OznB2dpbb90WPHz+Gr68vTp06haSkJBQUFODdd99Fhw4d4OrqWupjAvT2YSgREZHa4D0lIiJSGwwlIiJSGwwlIiJSGwwlIiJSGwwlIiJSGwwlIiJSGwwlIiJSGwwlIjWQkZGBdu3awdLSEkFBQWXu07dvX7mnz1pZWcHBwQGffPIJDh069Nr39/Pzg6WlJTp06CAurkqkjrggK5EaOHjwIAoKCmBqagp/f/9XPtCvcePG4vpxRUVFePz4MQIDAzF79mzIZDJMnDixzOMCAgLQokULJCQk4MiRI3BxcVFWV4gqhSMlIjXg7++PLl26YMKECYiKikJCQkKZ+xkaGsLJyQlOTk5wcXGBu7s79uzZA21tbezfv7/MY27evInr16/jiy++QJs2beQWYiVSNwwlIhW7fv06YmJi4OLigqFDh0JLS6tcwVG3bl3o6Oi8cpFYf39/6OnpYcCAAXBxcUFUVBTi4+OrqHqiqsVQIlKxF0Ojfv366N27Nw4cOICioqJS+xYWFoqPTk9JSUFMTAwWL16MnJwcjB49utT++fn5OHjwIAYOHAg9PT0MGTLktaMqIlXjPSUiFcrLy8OhQ4fg6OgoPjfJ2dkZYWFhiIiIKPVIiDt37pR6plGdOnXg4eGBESNGlHr/8PBwpKenw9nZGUDxk3d79eqFwMBAzJw585WrcROpCkOJSIX+/PNPZGRkiKEBAL1794axsTECAgJKhVKzZs3g6ekJoPgprcnJydi/fz9WrFgBLS0tDB8+XG5/f39/NGvWDJ07dxbbXFxcEB4ejoiICPTu3Vt5nSOqAIYSkQr5+/ujQYMGaNy4Me7duye2d+vWDUeOHMGTJ0/QoEEDsV1PT6/Uc42GDh0KFxcXeHp6om/fvqhfvz6A4uc1nT17FiNGjJCbONGyZUsYGBjA39+foURqh6FEpCKJiYmIjIyEIAhwdHQsc5/g4OBXTvMuoaWlha5du2LHjh24cuWKOLrav38/ioqKsHfvXuzdu7fUcSdOnEBqaiqMjY0r3ReiqsJQIlKR/fv3QxAEeHp6wtDQsNT2n3/+GQEBAW8MJaD4ibsAkJOTA6D40l5gYCBat24t92TZEikpKfj2228RFBSEyZMnV64jRFWIoUSkAkVFRQgMDIRUKi1zggJQ/GhzLy8vXLlyBdbW1q98r7y8PERERAAA2rRpAwA4c+YM7t+/j/Hjx2PgwIFlHufj44OAgACGEqkVTgknUoFTp07h4cOHr7xsBwADBgwAUHzfqURWVhaCgoIQFBSEAwcOYMuWLXB1dUViYiJGjhwJMzMzuWM+/PDD175/XFwcLl26VAU9IqoaHCkRqYAioSGVSmFmZobQ0FAsWrQIAPDo0SPMnz9f3Oedd96BhYUFli9fLn5OKT09HeHh4fjggw9gamr6yvcfMGAAvL29ERAQAFtb26roFlGlSQRBEFRdBBEREcDLd0REpEYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDYYSkREpDb+D1F4fJnNtgv/AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(ABAtest)\n",
    "plt.xlabel(\"test set AUROC\\nABA\")\n",
    "plt.ylabel(\"brain areas\")\n",
    "plt.text(0.8, 87, \"mean = %.3f\" %np.mean(ABAtest),size=20)\n",
    "#plt.xlim([0.45,1.05])\n",
    "plt.ylim([0,100])\n",
    "plt.savefig(\"051321_2.pdf\", bbox_inches='tight')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}