genome-wide.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import researchpy as rp\n",
    "from scipy import stats\n",
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [],
   "source": [
    "df = pd.read_csv('/home/masha/FreePSI/FreePSI/example/psi_freePSI.summary', sep = '\\t')\n",
    "df1 = pd.read_csv('/home/masha/FreePSI/FreePSI/example/TruePSI', sep = '\\t')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Chromosome</th>\n",
       "      <th>Strand</th>\n",
       "      <th>Exon starting site</th>\n",
       "      <th>Exon ending site</th>\n",
       "      <th>PSI</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5022492</td>\n",
       "      <td>5022693</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5025008</td>\n",
       "      <td>5025049</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5026279</td>\n",
       "      <td>5026431</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5026431</td>\n",
       "      <td>5026630</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5027934</td>\n",
       "      <td>5028225</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5032052</td>\n",
       "      <td>5032217</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5033407</td>\n",
       "      <td>5033443</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5034581</td>\n",
       "      <td>5036782</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>5040584</td>\n",
       "      <td>5040668</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6059028</td>\n",
       "      <td>6059063</td>\n",
       "      <td>0.400740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6060339</td>\n",
       "      <td>6060413</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6060413</td>\n",
       "      <td>6060499</td>\n",
       "      <td>0.400740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6062523</td>\n",
       "      <td>6062591</td>\n",
       "      <td>0.400740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6066866</td>\n",
       "      <td>6067049</td>\n",
       "      <td>0.400740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6076003</td>\n",
       "      <td>6076171</td>\n",
       "      <td>0.400740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6076171</td>\n",
       "      <td>6076305</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6111072</td>\n",
       "      <td>6111129</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6111129</td>\n",
       "      <td>6111252</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6112069</td>\n",
       "      <td>6112239</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6112738</td>\n",
       "      <td>6112855</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6116223</td>\n",
       "      <td>6116287</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6116461</td>\n",
       "      <td>6116623</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6116893</td>\n",
       "      <td>6117018</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6117127</td>\n",
       "      <td>6117251</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6117803</td>\n",
       "      <td>6118027</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6118255</td>\n",
       "      <td>6118402</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6118782</td>\n",
       "      <td>6118908</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6119023</td>\n",
       "      <td>6119240</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6119719</td>\n",
       "      <td>6120001</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>chr21</td>\n",
       "      <td>+</td>\n",
       "      <td>6120486</td>\n",
       "      <td>6120718</td>\n",
       "      <td>0.245082</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2676</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46265924</td>\n",
       "      <td>46266166</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2677</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46266981</td>\n",
       "      <td>46267142</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2678</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46270400</td>\n",
       "      <td>46270563</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2679</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46272560</td>\n",
       "      <td>46272829</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2680</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46273387</td>\n",
       "      <td>46273585</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2681</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46275185</td>\n",
       "      <td>46275325</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2682</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46277526</td>\n",
       "      <td>46277717</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2683</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46279992</td>\n",
       "      <td>46280137</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2684</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46280496</td>\n",
       "      <td>46280575</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2685</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46283614</td>\n",
       "      <td>46283838</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2686</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46284067</td>\n",
       "      <td>46285394</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2687</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46300429</td>\n",
       "      <td>46300794</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2688</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46301129</td>\n",
       "      <td>46301713</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2689</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46301713</td>\n",
       "      <td>46301996</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2690</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46301996</td>\n",
       "      <td>46302154</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2691</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46302484</td>\n",
       "      <td>46302576</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2692</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46311455</td>\n",
       "      <td>46311567</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2693</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46314715</td>\n",
       "      <td>46314880</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2694</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46314880</td>\n",
       "      <td>46315049</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2695</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46315473</td>\n",
       "      <td>46315547</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2696</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46317207</td>\n",
       "      <td>46317268</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2697</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46317268</td>\n",
       "      <td>46317348</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2698</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46318011</td>\n",
       "      <td>46318220</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2699</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46320636</td>\n",
       "      <td>46320825</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2700</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46322638</td>\n",
       "      <td>46323136</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2701</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46323136</td>\n",
       "      <td>46323857</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2702</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46323857</td>\n",
       "      <td>46323899</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2703</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46598617</td>\n",
       "      <td>46599503</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2704</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46602277</td>\n",
       "      <td>46602416</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2705</th>\n",
       "      <td>chr21</td>\n",
       "      <td>-</td>\n",
       "      <td>46605012</td>\n",
       "      <td>46605122</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2706 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     Chromosome Strand  Exon starting site  Exon ending site       PSI\n",
       "0         chr21      +             5022492           5022693  1.000000\n",
       "1         chr21      +             5025008           5025049  1.000000\n",
       "2         chr21      +             5026279           5026431  1.000000\n",
       "3         chr21      +             5026431           5026630  1.000000\n",
       "4         chr21      +             5027934           5028225  1.000000\n",
       "5         chr21      +             5032052           5032217  1.000000\n",
       "6         chr21      +             5033407           5033443  0.000000\n",
       "7         chr21      +             5034581           5036782  1.000000\n",
       "8         chr21      +             5040584           5040668  0.000000\n",
       "9         chr21      +             6059028           6059063  0.400740\n",
       "10        chr21      +             6060339           6060413  0.000000\n",
       "11        chr21      +             6060413           6060499  0.400740\n",
       "12        chr21      +             6062523           6062591  0.400740\n",
       "13        chr21      +             6066866           6067049  0.400740\n",
       "14        chr21      +             6076003           6076171  0.400740\n",
       "15        chr21      +             6076171           6076305  1.000000\n",
       "16        chr21      +             6111072           6111129  1.000000\n",
       "17        chr21      +             6111129           6111252  0.245082\n",
       "18        chr21      +             6112069           6112239  0.245082\n",
       "19        chr21      +             6112738           6112855  0.245082\n",
       "20        chr21      +             6116223           6116287  0.245082\n",
       "21        chr21      +             6116461           6116623  0.245082\n",
       "22        chr21      +             6116893           6117018  0.245082\n",
       "23        chr21      +             6117127           6117251  0.245082\n",
       "24        chr21      +             6117803           6118027  0.245082\n",
       "25        chr21      +             6118255           6118402  0.245082\n",
       "26        chr21      +             6118782           6118908  0.245082\n",
       "27        chr21      +             6119023           6119240  0.245082\n",
       "28        chr21      +             6119719           6120001  0.245082\n",
       "29        chr21      +             6120486           6120718  0.245082\n",
       "...         ...    ...                 ...               ...       ...\n",
       "2676      chr21      -            46265924          46266166  0.000000\n",
       "2677      chr21      -            46266981          46267142  0.000000\n",
       "2678      chr21      -            46270400          46270563  0.000000\n",
       "2679      chr21      -            46272560          46272829  0.000000\n",
       "2680      chr21      -            46273387          46273585  0.000000\n",
       "2681      chr21      -            46275185          46275325  0.000000\n",
       "2682      chr21      -            46277526          46277717  0.000000\n",
       "2683      chr21      -            46279992          46280137  0.000000\n",
       "2684      chr21      -            46280496          46280575  0.000000\n",
       "2685      chr21      -            46283614          46283838  0.000000\n",
       "2686      chr21      -            46284067          46285394  0.000000\n",
       "2687      chr21      -            46300429          46300794  1.000000\n",
       "2688      chr21      -            46301129          46301713  1.000000\n",
       "2689      chr21      -            46301713          46301996  1.000000\n",
       "2690      chr21      -            46301996          46302154  1.000000\n",
       "2691      chr21      -            46302484          46302576  1.000000\n",
       "2692      chr21      -            46311455          46311567  1.000000\n",
       "2693      chr21      -            46314715          46314880  1.000000\n",
       "2694      chr21      -            46314880          46315049  1.000000\n",
       "2695      chr21      -            46315473          46315547  1.000000\n",
       "2696      chr21      -            46317207          46317268  1.000000\n",
       "2697      chr21      -            46317268          46317348  1.000000\n",
       "2698      chr21      -            46318011          46318220  1.000000\n",
       "2699      chr21      -            46320636          46320825  1.000000\n",
       "2700      chr21      -            46322638          46323136  1.000000\n",
       "2701      chr21      -            46323136          46323857  1.000000\n",
       "2702      chr21      -            46323857          46323899  1.000000\n",
       "2703      chr21      -            46598617          46599503  0.000000\n",
       "2704      chr21      -            46602277          46602416  0.000000\n",
       "2705      chr21      -            46605012          46605122  0.000000\n",
       "\n",
       "[2706 rows x 5 columns]"
      ]
     },
     "execution_count": 66,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [],
   "source": [
    "x, y = pd.Series(df['PSI'], name=\"Estimated PSI by freePSI\"), pd.Series(df1['PSI'], name=\"True PSI in the simulated date\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f2548309908>"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xl8VPW5+PHPM0kmK9nYIYEkLFKxKooLioiKLdYWaxerdtO6oF2wentv9Xdbfr3c9qftbWuxWoVS22pbbdXW4nVHBETRgriBsoYtbGEJkH0yM8/vjzMTJpPJZJLMTCbheb9eeWXOmTPnfE+Weeac7/f7PKKqGGOMMbFw9XYDjDHG9B0WNIwxxsTMgoYxxpiYWdAwxhgTMwsaxhhjYmZBwxhjTMwsaBhjjImZBQ1jjDExs6BhjDEmZum93YB4GzRokJaVlfV2M4wxpk95++23D6rq4M6263dBo6ysjDVr1vR2M4wxpk8RkR2xbGe3p4wxxsTMgoYxxpiYWdAwxhgTMwsaxhhjYmZBwxhjTMx6dfSUiDwMfBqoVtVTIjwvwHzgU0ADcJ2qrk1EW65Z8AarttW0Lk8pL+Kx2ecl4lBJtWxDNQtWVLKq8lDE5++YMY45M8YnuVXxN/XuJVQdbW5dLinIZOVdMxJ2vPuWbGLRym3Ue3zkutO4cWo5c2aMb/1576ppoLQoh9nTKpg+YUib15bd+Wy7/W2/5/Jut2XmvcvYsL++dXnC0FxeuH16p22FxPzdRzueSYx4/01FI71ZuU9EpgF1wCMdBI1PAd/BCRrnAPNV9Zxo+5w8ebJ2dcht+D9OUF8PHMs2VDN38Xp2Hm6Iul1fDxzhASMoUYHjviWbmL90Cy4Bl4Bfna9Zpw7j7Z1HyUgTsjPSaGzx0eJT5s2a2Bo4Iv1zB3Xnnzw8YAQFA0dHbb3t4rGs2now7n/30Y7Xl//GUlm8/qZE5G1VndzZdr16e0pVVwCHo2xyBU5AUVV9EygUkeHxbkekf5xo6/uKBSsqyUiTTrdbtHJbElqTOJECRrT1PbVo5TZcAukuFy5xBb7D4vf3kZEm5LjTEXG+Z6QJC1ZUJqQdQMSAEbq+o7YuWrktIX/30Y5n+odU79MYCewKWa4KrGtDRG4WkTUisubAgQNJa1yq21XTQHZGWqfb1Xt8SWhN/1Hv8eEKi8UuAZ9f2/28szPSqKqJfqWXSB21NVG/82QfzyRfqgeNSB+T291PU9WFqjpZVScPHtzpLPgTRmlRDo0tnf+z5ro7DyzmuFx3Gv6wv0K/QppL2v28G1t8lBTlJLF1bXXU1kT9zpN9PJN8qR40qoDSkOUSYE+8DzKlvKhL6/uK2dMqaPF13md149TyJLQmcUoKMru0vqdunFqOX8Hr9+NXf+C706fR4lMaPF5Une8tPmX2tIqEtAOcvoto6ztq641TyxPydx/teKZ/SPWgsRj4mjjOBY6q6t54H+Sx2ee1+0fp653gANMnDGHerIlMqRjY4TZ9vRMcYOVdM9oFiESOnpozYzy3XTyW7Iw0vH7nFtRtF4/l3qvPYN6siQwZkMXRxhaGDMhq0wkOHXdMdnekywu3T28XOEJHT3XU1jkzxifk7z7a8UxixPtvqjO9PXrqMWA6MAjYD/xfIANAVR8KDLm9H5iJM+T2elWNOjSqO6OnjDHmRBfr6Klenaehqtd08rwC30pSc4wx5oR0tLEl5m37XWp0Y4wxsWn0+DhY10yLzx/zayxoGGPMCcbj9XO43kODx0vlgbouzSWyoGGMMScIv1+pafBwrMnLobpmHn59Gy+s29dumHQ0FjSMMeYEcKyphZp6D/XNXp54u4rH/rWTphbntlRpUTYxle3DgoYxxvRrjR4fh+qbaWrxseSjaha9VsnBOg8ABdkZfH3KaD596nBOujO2/VnQMMaYfqjF5/Rb1Dd7eXfXER5ctpXN1XUAZKQJn5s0ki+fM5q8rHQyY0g3FGRBwxhj+pHQfoudh+pZuKKS17ceL41w0UmDufGCcoYXZJOR5qIo101eZuyhwIKGMcb0E8F+i8N1Hh55cweL39uDL9DLffLwfL45fQwnj8gnzSUU5rjJz3IyMneFBQ1jjOnjgv0WdU1e/vHObv701g7qm53kmcMLsrjpggouHD+INJeLguwMCrIzcIWnI46RBQ1jjOmjgv0WdU0tLN90kN++Vsneo00A5Gam8ZVzRnPlpJFkZqSRl5lOca6btG4GiyALGsYY08f4/cqRxhaONrawfvdRHly+lfV7jgFOiv5Zp43ga+eOpiAng9zMdIpy3LjT45Of1oKGMcb0IcF+i6qaBha9to1XNx4vPHf+mIHcPK2C0uIcsjLSKM51k9WFkVGxsKBhjDF9QFOLkyfqcJ2HP7+1g7+/s7u1Xs7YIXl8c/oYTi8tJCPNRXGum9wujIjqCgsaxhiTwoL9FkcbPDzz/l7++MZ2jjV5ARiU5+bGqeXMOHkoGWmubo+I6goLGsYYk4KC/RZHGjy8sfUgC5ZXsqumEYCsDBfXnD2KL55ZQrY7nYLsDAp7MCKqKyxoGGNMijnW1MKR+hY+2nuUB5dX8u6uIwC4BC47ZTjXn19Gca6bvKx0inPcpKclrwirBQ1jjEkRwX6L3TWNPPz6Nl5av59gAtqzyoqYPa2CisF55LjTKcrNIDM9vp3csbCgYYwxvSzYb3Gwtpm/rt7F39bsosnrZKAtG5jDLReO4ezyYtzpLgbmZpLtTn6wCLKgYYwxvSTYb3G43sML6/by+9e3c6jeyUBblJPB9eeXcdkpw8lMT6MoN4MBWRm93GILGsYY0yuC8y3eqjzEQysqqTxQD4A73cUXzyzhmrNLycvMoDDHSfuRyBFRXWFBwxhjkiiYJ2rT/loWLK/krW2HW5+79OSh3HB+GUMLshmQ5czk7mnaj3izoGGMMR1YtqGaBSsq2VXTQGlRDrOnVTB9wpBu7SvYb1FV08Af3tjOs+/vbS2zelpJAbdOH8P4oQPIy0ynKNdNRhJHRHWFBQ1jjIlg2YZq5i5eT0aaUJidQXVtE3MXr2cedClwBOtbHKht5sm3q/jLv3bS4HEy0JYUZXPzBRWcP3Yg2e70bqf9iGdw64wFDWOMiWDBikoy0oQct/M2meNOp8HjZcGKypjfkI81tXCorpmXP9zPote2UV3bDEB+VjpfmzKaz5w2gpxAsOhu2o94BbdY9WrQEJGZwHwgDVikqveEPT8K+CNQGNjmTlV9LukNNcaccHbVNFCY3Xa0UnZGGlU1DZ2+NthvsWb7YR5cXsnGfbUApLuEKyeN5CvnjqIwxx2XtB/xCG5d0WtBQ0TSgAeAS4EqYLWILFbVD0M2+wHwN1V9UEROBp4DypLeWGNMyrtvySYWrdxGvcdHrjuNG6eWM2fG+G7vr7Qoh+raptY3Y4DGFh8lRTkdvsbjdfotNu+vZeFrlby2+WDrc9PGDeKmaRWUFOXENe1HT4Jbd/TmlcbZwBZVrQQQkceBK4DQoKFAfuBxAbAnqS00xvQJ9y3ZxPylW3AJpLucN/f5S7e0Pt+dYDJ7WgVzF6+nweMlOyONxhYfLT5l9rSKdtsG+y2qahp5ZNV2/vnuHryBXu4JwwZw64Vj+HhJQULSfnQnuPVEbwaNkcCukOUq4JywbX4EvCQi3wFygRnJaZoxpjPx/mTfE4tWbgsEDOfN2CXg9ft5YNlWvH6NGEw6a+v0CUOYh3P7p6qmgZIOOpiPNbVQfayJp9bu5k9v7qA2kIF2yIBMbrqggosmDCY3UDUvEWk/uhLc4qE3g0ak6zINW74G+IOq/kJEpgCPisgpqupvsyORm4GbAUaNGpWQxhpjjov2yb43Ake9x0ewMJ3Pr3j9/sBw1kDASHPerIPBZNHKbTG1c/qEIR32CzR6fBysa+KVj6pZ+Fole44Eyqy607j2nFF8/owS58oi193mKiDeYg1u8dKbQaMKKA1ZLqH97acbgJkAqrpKRLKAQUB16EaquhBYCDB58uTwwGOMibOOPtnH+mbcVRV3Pos/bF2aS1qvcHLdzidsVaXF13ZLv4LX52+9JeQSJ8h0d5hqcL7F2zsO8+CyrXyw+1jrfj9z6gi+ft5oBuVlUZibQX6S0n5EC27x1ptBYzUwTkTKgd3A1cC1YdvsBC4B/iAiHwOygAMYY3pV6Cf7oOCbcbxFChjgXFEEr3DOHl3Iv3YcweePtCX4VFvf7PwK7jTp8jDVYL/Fpv11LHqtklc2HP/sem5FMbOnVVA+KC/l0n7EW68FDVX1isi3gRdxhtM+rKrrRWQesEZVFwP/BvxWRG7HuXV1naralUSKu/3xtSx+fx8+vzpF7k8dxr1Xn9HbzTJxFPxkHzr4x6/O+niLHAYc6S4XXr+f9Xtrue3isdz7ymZUnQA2MNeNS4T9tc34FfzqD3yHQXmZXRqmeqypharDDTz65g6efLuqtczqmMG53HrhGM4sK07ZtB/x1qvzNAJzLp4LWzc35PGHwPnJbpfpvtsfX8s/3t3buuzza2B5rQWOfuTGqeXMX7oFr9+PS2h9M75xannS2+ISONbkZdHKbW0CxrCCbMC5nXS0yYvXT+vtrL+9XUV22MzrSMNUGz0+9tc28vQ7e/jD69s50tgCOPv/xtRyPnHyUGf4bI4bd/ilVz/VYdAQkaj/4aq6Nv7NMX3d4vf3ARB6Za7qrL/36l5qlIm7YL9Fb4+eamw5fjvsWJMXwQleB+o8NLd4GZyfTV5WBr+86vQ2VxCrKg9HHabq8fo5VNfMqxudfo8dh5xgkpXu4ktnlXLVWaUU5bi7nfajL4t2pfGLKM8pcHGc22L6AZ8/8t3DjtabvmvOjPFJCRIuot+iChX6V3as2c8Il/DDy09ud8upo2GqN00t52BdM+/uOsKDr27h7Z1OmVUBZp4yjOvPL2N4QTYD8xI7IiqVdXjWqnpRMhti+oc0l0QMEP39Pq9JnMp7Lu+wMzwSEUCdAFKUmxmxjyJ8mOrIwmy+eu5oCnMzmPv0Op5ft681AJ0xqpBbLhzDhGH5SR0Rlaqi3Z46C9ilqvsCy18DPg/sAH6kqoc7eq05cc06dRj/eHcv4cMVZp06rHcaZLokmdlSu6LynstbH8caQASoqmno8JyCXw0eL1WHG/nLv3bw+JO7aGpx9j66OIfZF1YwpWIghTluCuKU9qOvi3Z9tYDADGwRmQbcA3wHOB1nTsQXEt460+c4nd02eqqvmXr3EqqONrcuu4DMdFdCs6V2V6cBI3CVUZyTQV5meodDa88bO4iDdc08/c5ufvf6Ng7VOWVWC7IzuO680Xz61BEU5bpPiBFRXREtaKSFXE18CVioqk8BT4nIu4lvmumr7r36DOv07kPCAwY4b8zbDtRTPjg3YdlSE6k4J4O8rAxUtd3Q2vrmFu5/dQuHGzz85tWtbDlQB0BGmvD5M0q49pxRDM3PougEGhHVFVGDhoikq6oXZ4LdzTG+zhjTh4QHjCA/kYeh9vYtrECXRURTKga2SaXxg3+ua80Aq6r4FXx++GD3Ue7423utr7t4whBunFpO2aDcPjkiKplzo6K9+T8GLBeRg0Aj8BqAiIwFjiakNcaYlBKeLTXeBX/Cr3JKCjJZeVfHeUmXbagmM11o8rYPG1eePrzdG2XpCicDbFZ6Gk1eHwfrPBwNzLUAmDgin29OH8OpJYU9KoTUm5I9Nyra6KmfiMgrwHDgpZCZ2C6cvg1jTIqK1yfP8GypPS34E3qVsu9II+Hv/VVHmym781mmlBfx2Ozz2r1+wYpKRhblcOBYE8eaj8/RGJiTHvH8bji/jB8uXsfBumZqm7ytNbmLc9185+KxXHzSEIry3H16RFRowAhfn4jbxNFGT2UB5wJjgSEi8jtV9arqpvg3wxgTL1395FlSkBnxFlW6wLxZE9sEg54U/Am/Sqmqaexw21XbarhmwRvtAkfw+KMH5bWuU9U2Vw/g5Ik6XN/MjpoGmlr8HG100pWniTBz4lBumzGeIfmZ/TpHVKJEuxb7I9CCc1vqMuBk4LZkNMoY031dnZW/8q4ZMd8m6knBn/CrlM6s2lZD2Z3Pti6nC5xVPrDT4x9rauG1TQd44NUtfLjXKbOa5hKuOH0EX5tSRmlRNoU2Iqrbov32TlbVjwOIyO+AfyWnScaYnujOrPxo/QihelLwJ9JVSld4Fd6qPERJcU7E4zd6fLy/+wi/eXUryzcdT4Y9dewgbp5WzoRh+RTlusmIY9W8E1G0oNF6vRfISJuE5hhjeiqRs/J7UvAn/CqlK+lBgvw4t8xCj3/D+WWMLMpm7j/X8fS7u1sz0Ka7hLKBuXxu0kjOLh/Y50ZEpapoQeM0ETkWeCxAdmBZAFXV/I5faozpLYmeld/dgj/hVynlg3PZeqC+28f3+5Xq2ib+/NZO7njiPY4Fyqy6BAqzMxiY68arfn71ymaKc919bq5Jqoo2esrCsjmh9Jc6IKk6Kz/SVUowmeA1C95g1baamPd1pMHDs+/v5cHlW1s71LMz0ijMziAtDfKzjvdZdGV0l+lcTD1SIjIVGKeqvxeRQcAAVd2W2KYZkzz9rQ5Iqs7K7+gqJXyU1Ni7nm03HBecam1LPtzPr5du5r0qZ7qYS+BTHx/ON84v59t/WUtBVtscUbGO7jKx6TRoiMj/BSYDJwG/B9zAn7DiSKYfsTogidXVq7gtd1/eLnC4gEtOHspNj6xpnRF+dnkxt1w4htNKCijMcTN6YG6g3+R4Z3eso7tMbGK50rgSmASsBVDVPSIyIKGtMibJrA5I4nT3Km7L3Zfj8yu7axr47Wvb+NuaXbz04X4AKgblMvvCCi46aUibEVE9Gd3VV3U0z6akIDMhx4slaHhUVUVEAUQkNyEtMaYXWR2QxOnOVZyqUtPQwl/e2sHvVm6jpsEZzFmc6+b688r47KSRDB6Q2W5EVE9Gd/VVXZlnEw+xBI2/icgCoFBEbgK+Afw2Ia0xppdYHZDE6epVXH2zl+fX7eWBV7ey7aAzuioz3cVVk0v4yrmjKSnKiZojqruju/qyRAWISDoNGqr6cxG5FDiG068xV1VfTnjLjEmiVB1xlKpiyXR735JNLFp5fLyMaturjfCruGavj9XbDzN/yWZWbz8+kuoTJw/lpmnljB+aT35WuqX9iCCZmYdFwz9aRdpIZDTO6KklIpKDU2ujNiEt6qHJkyfrmjVrersZxvRboTmkQvsNQvNU3bdkE/OXbsEl4PNpxEl8way0Pr+yaX8t9y/dwvPr9rYmFTy9tIBvTh/LWWXFVjUvilh+H7EQkbdVdXJn28UyeuomnFoaxcAYYCTwEE6NDWPMCSaWTLeLVm7DJZDucpHuAo/Xjy/wATV4FffLL01i79FGFr22jb+8tZPGFidrbUlRNrdcOIZPThxKcW6m9St1oqeZh7sqlj6NbwFnA28BqOpmETmxbhgaY1rtqmmgtrGFo03HZ3MXZKXj8R6/nqj3+AgteudOd+FXP14/bP1/n+JYYwuLVm5j4fJKDtQ5Hbj5Wel8bUoZ15xdyuABWVY1L0Y9yTzcHbEEjWZV9QTvI4pIOh0XzuoSEZkJzMeZs7NIVe+JsM1VwI8Cx3xPVa+Nx7GNOZEE+xfqPT5y3WncOLWcOTPGt9kmlvviyzZUs/9YU2t+p6CjTV4yQt7kc93ObZLQiwS/Qk6Gi2ff38P8Vzazaf/xMqtXThrJN6aWUzYw13JEdVFpUQ7bDtZR2+TF4/PjTnMxICud8pD08fEUSyhfLiL/Byf31KXAE8AzPT2wiKQBD3A87fo1InJy2DbjgLuA81V1IvDdnh7XmBNNsH+hscX59N/Y4mP+0i3ct+R4aZzgffHq2qY2FfmWbahut014wAg6WOdp3f7GqeX4Fbx+P3710+Lz4fUpRTluvvWXd1oDxvTxg/nTDecw9zMTmTAs3wJGN0ypKOZAnQePz49LwOPzc6DOw5SK4oQcL5agcSdwAPgAmA08B/wgDsc+G9iiqpWq6gEeB64I2+Ym4AFVrQFQ1WqMMV0S2r/gElfgO21GNoXeFxdxvmekCQtWVLbbJppgoJkzYzy3XTyWrHQXHq8CAgI7A3miTh4+gAeuncT8qydxdnkxeX2wzGqqWFV5mIKsdHx+pdmr+PxKQVY6qyoPJ+R4UX9TgauBP6rqV4j/3IyRwK6Q5SrgnLBtxgfa8TrOLawfqeoLcW6HMf1aeP8COPma6j3Hy6XGcl88lnoYh+ubWztgrz13NE1eP4+s2kFds5OBdlh+FjdPq2DWacMpsk7uuNhcXUttk5eMNBciztDm2iYvm6sTM8A1atBQVZ+IDBYRd+BqIJ4i/bWEX/emA+OA6UAJ8JqInKKqR9rsSORmnBFejBo1Ks7NNKZv66h/Idd9/FZQLBX5gtsUZmdwJKy8alBds48dh+p4dNV2frNsK3uPNrW24cvnjuar545iWEG2FUKKI4/Xj18Vn09b58JIYH0ixPKb2w68LiI/FJE7gl9xOHYVUBqyXALsibDNP1W1JZBVdyNOEGlDVReq6mRVnTx48OA4NM2Y/iO8f8H57qwPmj2tghaf0uDxoup8D8/ZFNxmYJ476vH213r44T/Xs/doEy6Bz54+gidvncIdl46ntDjXAkac+dWPT50PAorz3aegmpigEcuNxD2BLxcQz0SFq4FxIlIO7AauBsJHRj0NXAP8IZCSfTxQiTEmZsFRUtFGT8WSsyl0m2jFk4LpQaZUDOTbF41l0ujCmOuCm65ziQuXOAEidNa9SGKCc4e/SRF5VFW/ChxR1fnxPnCghOy3gRdx+iseVtX1IjIPWKOqiwPPfUJEPgR8wL+r6qF4t8WY/m7OjPHththCR0Nxz+1wP8G8TuGZa8Pd+6XTuORjQ8nP6n5NcBMbd7qLNI/gcklrn4bfrwmb59JhGpHAG/VlwGKcPoU2fRCqmpiu+R6yNCLGxCY01YdLwOPV1k7FznJvHWv0cOp/RU5BJzgT+CztR3Jcs/BNth+q41jj8Xka+dnplA3M47GbO/4AEC4eaUQeAl4AKoC3aRs0NLDeGNNHhQ7F9Xj9bUahdFTzornFx+L39vDrpVs63O9nTx9uASOJgjVEhhWkJ6WGSIfXL6p6n6p+DOe2UYWqlod8WcAwpo+r9xwfUeULu+MQvC8erIXh8ysrNx/kmt++yb8/+T47DzfgEhhRkNX6aTLNJa1JCE3yTJ8whHmzJjJkQBZHG1sYMiCry8kKuyKW1Oi3JuTIxpheFWkobjgnA+0x7n15My+s29d6NTJ5dBFzLhnHWWXFZLttFndvS2YNERvSYMwJ6sap5cxfugWvP/rQzFm/fp2mwJj/0QNz+PZFY7ns48NtFvcJyn7rxpygQofien3eNn0aoXermrx+inIy+Mb55Vx7ziiKc91WCOkE1mkRpsCw2D8H8z+lOhs9ZUz33PbY2yx+b1+b4OFOd/HFM0u4eVoFJUU5lvajH4tbESZgGLBaRNYCDwMvaizl/owxfYKq8s7OGvbXetoEjBkfG8Kci8fxsRH5NovbtIqlI/wHIvJD4BPA9cD9IvI34HequjXRDTTGdE8sNTR2HKrn3pc3sfi9Pa1lVj8+Mp/bZoxn6thBlqrctBNTn4aqqojsA/YBXqAIeFJEXlbV/0hkA40xXRc6cS+0hgY4fRlHGzw8uHwrj6zaQUMg2+2Iwiy+OX0snz19BHk2k9t0IJYa4XOArwMHgUU4qTxaxElsshmwoGFMigmduAfOjG+v38/CFVspzHVz/9ItVNc6ZVYHZKXz9Smjue78cgZaJ7fpRCxXGoOAz6nqjtCVquoXkU8nplnGmJ4Ir6GhqqgqdR5l7j/XA5DuEj47aSTfumgso4tzbBa3iUksfRpzReQMEbkCJ33I66q6NvDcR4luoDGm60In7vn8frx+be2zAJg2bhC3zRjHaSWFpFsnd58XS333eInl9tQPgauAvwdW/V5EnlDVHyekRcaYHrtxajn3LtmM1+drMyJqUK6bH195ChdNGEJmunVy9wfLNlRz3R9Wty5X1TSyqvIQf7jurIQEjljmaXwETFLVpsByNrA2kJcq5dg8DZPKZt67jA37j9eimDA0lxdunx6Xfd/++FoWv78Pn18RnPxRwasLAS4cP4j7rz3DOrn7mbI7n+3wue33XB7zfuI5T2M7kAU0BZYzARtqa0wXhQcMgA3765l577IeB47w+haKM6s7zSXcMLWcGy8oZ8iArB4dwxiIXoTp1zh/e83AehF5ObB8KbAyOc0zpv8IDxidrY+VqvJ0BwWR/H7lzpkTrJPbxE20K43gPZ63gX+ErF+WsNYYY7pkS3UtP3thIx3dZFawgGHiqsOgoap/TGZDjDGxO1jXzH2vbObx1bvweNtnqQ2W/bRcUf1fnttFnaf930CeOzGj4mysnTFJMmFobpfWR9Lk8bJg+VYu/eVyHlm1A4/Xz8BcNxOH57XZLji+Zdapw7rdXtM3rJt3WbsAked2sW7eZQk5Xqejp/oaGz1lUll3R0/5/cpzH+zl5y9tZPuhBgCy0l1cffYobrmwgmEF2W1GT3VW49uYcLGOnrKgYUwvO2Xu821uL4R/Snx7x2HueX4Dq7c71QkEmHnKML57yTjGDxvQYdqPZE74Mn1fj4fcisgz0GH/Gqo6q5ttM8YEhAcMgDqPn1PmPs/z372Qn72wgWc/2Ns63+KMUYXccelJTBkzMGp/xbIN1cxdvJ6MNKEwO4Pq2ibmLl7PPLDAYXok2uipnyetFcb0cx2lKY/UgQlO4Lj03uU0tTjPjyrOYc7F4/jMacPJjCFd+YIVlWSkCTlu5188x51Og8fLghWVFjRMj0QbPbU8dFlEMoBTgN2qWp3ohhnTX0RLUx5NU4ufguwMbphaznXnlZGfHftM7l01DRSGbZ+dkUZVTUOX229MqA5HT4nIQyIyMfC4AHgPeAR4R0SuicfBRWSmiGwUkS0icmeU7b4gIioind5vMybVhKYpd4kr8N1ZH801Z4/i+dumMueScV0KGAClRTk0tvjarGts8VFSlNPl9hsTKtrtqQtU9ZbA4+uBTar6WREZBjwPPNaTA4tIGvAAzgzzKpySsotV9cOw7QYAc4Ca9yI/AAAgAElEQVS3enI8c2KJlI8nPys9YvW6RAtPUw5OfYt6j4/sdGj0tn9Ndjr8vytP6XZti9nTKpi7eD0NHi/ZGU7G2xafMntaRbf2Z1JbMgc9RJun4Ql5fCnwNICq7ovTsc8Gtqhqpap6gMeBKyJs99/Azzie+8qYqDpK4Hasycv8pVu4b8mmpLRj2YZqyu58Fp9fafYqjS2+1k//Xp/iAjz+9kEhN0P46MeX96gY0vQJQ5g3ayJDBmRxtLGFIQOymDdrovVn9EPBQQ/VtU1tBj0s25CYXoRoVxpHAkWW9gDnAzcAiEg6kB2HY48EdoUsVwHnhG4gIpOAUlX9XxH5XhyOaU5wwdtCib7aCE9XHSoYOHyB4e7DC7L49kVj+cKZJTF1csdq+oQhFiROAMke9BAtaMwG7gOGAd8NucK4BOg4F2/sIn2Mah3iGygney9wXac7ErkZuBlg1KhRcWia6a+Ct4USbcGKyk63yc1M4/rzyrnpggoKcixduemeZA96iDZ6ahMwU0QGqerBkPUvAi/G4dhVQGnIcgnOVU3QAJzRWssCl+nDgMUiMktV28zeU9WFwEJwJvfFoW2mn2r2KvlZsVQE6JldnfzDfuHMEm67ZCylxbGnEDEmktKiHKprm1qvNCCxgx6ijZ76jIgcAN4XkSoROS/Ox14NjBORchFxA1cDi4NPqupRVR2kqmWqWga8CbQLGMZ01Y1TyxN+jNJO/mH/5wunWsAwcTF7WgUtPqXB40XV+Z7IQQ/ROsJ/gjOCagTweeDueB5YVb3At3GuWj4C/qaq60VknojYbHPTbZ1VK0t0f8bB2mayMjr+18pMo0ed3MaESvagh2jX6V5V3QCgqm8Fhr7Glao+BzwXtm5uB9tOj/fxjYmnphYfv1u5jQUrtnIs0jhanICx8Sexl+A0JhbJHPQQLWgMEZE7OlpW1V8mrlnG9Ewyawz4/X6eeW8vP395I7sONwJOR+RXp4zmlmljKM5zx/2YxvSWaEHjtzid0R0tG5Oy1s27rNPssfHwr22HuPv5Dbyz8wjgjM769KkjuH3GOMoH53XyamP6nmijp/4rmQ0xJt4SVYQGYPvBen76wgZeWLevdZz4WWVFfH/mBM4cXWR9FqbfSvzYQ2P6kZp6D79eupk/v7WT5kCZ1fJBufzbpeO57OPDrbyq6RWRsiB0NiCkuyxoGBMDj9fPH9/YzoPLt3K43smwU5STwewLx/D1KWVku+M3k9uYrugobU7Znc8mJHBY0DAmClXl+XX7+PmLG6k86JRpdae7uOasUr5z8TgGDcjs5RZG1lH9DmN6Klrlvjs6eg5s9JTp/97ZWcM9z2/grW2HW9d9cuJQvveJkxg3NHXHhESr32GBw/RUtCuN4H/FScBZHJ+t/RlgRSIbZUxv2nW4gZ+/tJFn3tvTWmb1tJICvj9zAlPGDEz5Tu7Q+h3gjOjy+v1JSdRo+r9OR0+JyEvAGapaG1j+EfBEUlpnTBIda2zhwWVb+eOq7TQEkhqWFGUzcfgAXt9ykK8+/K8+casnWv0OY3oqlj6NUbStreEByhLSGmN6gdfr57HVO/n10i1U1zYDgYJNF1TgafHx4IrKPnWrJ9ftFF0KHcjlV2e96X+233N5yo2eehT4l4j8Ayd1+ZU4ZV+N6dNUlWUbD3DPCxvYuK8WgHSX8IUzS/jujPEMK8ji1B+92Odu9dw4tZz5S7fg9ftxiRMw/JqcRI2mdyQqQETSadBQ1Z+IyPPABYFV16vqO4ltljGJ9eHeo9zz3AZWbG7N+s/08YP5j5kncfKIgtZ1ffFWTzCY2egpkwixDrnNAY6p6u9FZLCIlKvqtkQ2zJhE2H+siV+8tJGn1u7GF+jlPnl4Pv/+yZOYftLgdp3cffVWz5wZ4y1ImIToNGiIyP8FJuOMovo9kAH8CacErOkjkll4vidm3ruMDfvrW5cnDM3lhdun93i/DR4vC1dUsui1bdQ1Oxloh+ZnMueScVx1ZgkZ6ZGDgN3qMaatWK40rgQmAWsBVHVPItKkm8QJFp7PSJM2hefnQUoFjvCAAbBhfz0z713W7cDh9/v5yXMbeHTVDjw+J+1HZrqLGy8oZ/a0CvKzo2egtVs9xrQVS9DwqKqKiAKIiJUb62OSXXi+u8IDRmfrO7Nq60Hu+vsHbD90vPRqXmYaA7IyOGt0cacBI8hu9RhzXCxB428isgAoFJGbgG8AixLbLBMv4enBAT4+siChhed725bqOn76wgZe/nB/67pcdxrDCrLIznD6KEID5tS7l1B1tLl125KCTFbeNSPp7Tamu5I55LbTijSq+nPgSeApnH6Nuap6X0JaY+IqUsAA+GD30YQWnu8th+uamfvPdVw2f0VrwEh3CaVFWZQPyiXHnY6ItAmY4QEDoOpoM1PvXpL09hvTHdESFiZCLB3hP1XV7wMvR1hnUlikgBGUyMLz3TVhaG7EW1EThka/I9rc4uPh17dz/9LNrUNhMwLzLbYeqONQvafNqKjQgBkeMII6Wm/MiS6W2peXRliXuOo2JikSWXi+u164fXq7ABFt9JTf72fxu7v5xL0r+OkLG6j3+BBgYK6bEUVZvL71EFPHDqLFpzR4vKg631MxYBrTV0TLcnsr8E2gQkTeD3lqAPB6ohtmEivVAkZQrKOk1mw/zN3PbeDtnTWt63Iz0xhRkE1WhjN8tsHjZVXlYebNmsiCFZVU1TRQksLDjY3pC6LdnvoL8DxwN3BnyPpaVT0c+SUmleS5XRFvUeW5Y7nATE07Dtbzsxc38twHe1vLrJ45uoidh+oZMiATl+v4uQX7LqZPGNJhkCgpyIx4K6qkIDXrZBjT2zp891DVo6q6XVWvUdUdQCNO7qk8ERmVtBaabls377J2ASLP7Upo7exEqWlo5sf/+yGf/NUKng0EjNHFOfzq6tP5683nMnbIAJq8bQNkLJ39K++a0S5A2Ogp05d0NEoqUaOnRFWjbyDyGeCXwAigGhgNfKSqExPSoh6aPHmyrlmzprebYeKkucXHn97ayW+WbeFQnZNsuTA7g5unVXDdeWXkZDoXy6ETGIPDalt8mpJ9N8akIhF5W1Und7ZdLPM0fgycCyxR1UkichFwTU8bCCAiM4H5QBqwSFXvCXv+DuBGwAscAL4RuOox/Zzfr7z84X7+56WNbKmuAyAjTfjSWaXMuXgcQ/Kz2mw/fcIQ5oH1XRiTYLEEjRZVPSQiLhFxqeqrIvLTnh5YRNKAB3BGZ1UBq0Vksap+GLLZO8BkVW0IdMz/DPhST49tUtt7u47w0xc28MbWQ63rZnxsKN/75HgmDMvv8HXR+i6MMfERS9A4IiJ5OCVe/ywi1Tif/HvqbGCLqlYCiMjjwBVAa9BQ1VdDtn8T+EocjmtSVFVNA794aRP/fHd3a5nVU0bm8x+fnMDUsYNwuVK7zKoxJ4JYgsYVQBNwO/BloACYF4djjwR2hSxXAedE2f4GnNFcpp852uDhoRWVPPLG9tbJeSMKs5hz8Tg+f8bIDjPQGmOSL5YiTPUAIpIPPBPHY0f62BixV15EvoKTnv3CDp6/GbgZYNQoG9jVV3i8Pv66ehf3v7qF/cecYa8DstK57rwybrqggvzsjF5uoTEmXCxpRGbjXFk0An6cN3sFejqltgooDVkuAfZEOP4M4D+BC1U1Ym4HVV0ILARn9FQP22USzO9XXt1Yzf+8uJENIWVWr5w0ku/OGMfIfpYTy5j+JJbbU98DJqrqwU637JrVwDgRKQd2A1cD14ZuICKTgAXATFWtjvPxTZKpKh/uPcbPXtjI8k0HWtdfMG4Q//7Jk/j4yIJ2lfOMMakllqCxFYh7Dm1V9YrIt4EXcYbcPqyq60VkHrBGVRcD/wPkAU8E3kx2quqseLfFJN6eI43MX7KZp9ZW4Q30cp80bAD/dul4LvnYUNKsk9uYPiGWoHEX8IaIvAW03h5S1Tk9PbiqPgc8F7Zubshjm5bbx9U2tfC7ldt4eOU2jjU5g+6GDMjkm9PH8KWzSsl2x1qm3hiTCmL5j10ALAU+wOnTMKZTHq+Pe57fwJ/f2kmz93iZ1a+fV8YtF1ZQnGu5nYzpi2IJGl5VvSPhLTH9gt+vvFF5iB8+vY5tB4/XxsjLTGNAZgbnVQy0gGFMHxZL0Hg1MKT1GdrenrJMt6aVqrJpfx2/eGkjL4WVWR2an0mOO71dmVVjTN8TS9AIjmi6K2RdPIbcmn7iQG0Tv166hcdX78ITuBWV7hKG5WdSmONuHRHVn+uSG3OiiGVyX3kyGmL6nvrmFh5ZtYOFKyqpaWgBoDjXzc3TKnjlw/0cbui4zKoxpm+KVrnvYlVdKiKfi/S8qv49cc0yqczr8/PsB3v51ZLNrf0Wmekurj6rlG9dNJYh+VlMGDqAuYvX0+DxtklVbmVWjenbol1pXIgzauozEZ5TwILGCcbvV97eUcPPXtzA6u1OmVUBPjFxKHdcOp7xQwe0XllYqnJj+qdYijCVq+q2ztalCivClBiVB+r45cubeO6Dva0ZaE8vLeDfLj2J88YOssl5xvRx8SzC9BRwRti6J4Ezu9Mw07ccqmvmwWVb+fNbO2lscTLQlhRlM+eSscw6bQRZGenct2QTi1Zuo97jI9edxo1Ty5kzY3wvt9wYkwjR+jQmABOBgrB+jXwgK/KrTH/R0OzlsdU7eWhZJQfqnJHW+VnpfGNqOdefX0ZBthuA+5ZsYv7SLbgE0l1OZ/f8pVsALHAY0w9Fu9I4Cfg0UEjbfo1a4KZENsr0Hq/Pz5KP9vOLlzaxOaTM6ufPKOFbF4+lpDCbO/76Dovf34fPf/zWpjvDqXnhEvD6/Sxauc2ChjH9UIdBQ1X/CfxTRKao6qoktsn0Ar9fea/qCD9/aSOvbzleZvWikwZz+4zxTBxZQJpLuP3xtfzj3b3tXt/Y4iM7JHAEiykZY/qXWPo0rhSR9Tj1NF4ATgO+q6p/SmjLTFKoKjsPNzD/lc388909rVcPJw8fwB2XnsS08YNxp7tat1/8/j4AgtMvIo2j8KszEzzI+jyM6T9iCRqfUNX/EJErcQonfRF4FbCg0cfV1Hv47WuV/HHVduqbnSuD4QVZ3Dp9DF84o4SczPZ/HqG3pMIFO8rTXMKNU505odbnYUz/EkvQCNbc/BTwmKoetkI5fVujx8tTa3fzwKtb2Hu0CYDczDS+dm4ZN1xQzqC8jhMKprmkNXB0NFo7O11aA8KildsCAcO5WrE+D2P6tliCxjMisgHn9tQ3RWQw0JTYZplEaPH5Wb7pAPe+vIn1e44BThCYddoIvn3xWMoH5uKKMN9i6t1LqDrattJutOk9dZ7jGfTrPT5C7m4B1udhTF8WS+6pO0Xkp8AxVfWJSANwReKbZuLF51fW7znKr5ZsZumG41Vzzx8zkO/OGM+kUYWkp7nave6k/3yW5h6+t+e6nRQiobEovM/DGNN3tH+nCBCR/whZnKGqPgBVrQd6XLXPJJ6qsutwPT/4xwd8/sE32gQMgDe2HuKs8uI2AeOaBW9QduezlN3Z84ABcOPUcvzq3JLyqz/wndY+D2NM39Jh0ACuDnl8V9hzMxPQFhNHh+ubue+VLXz616/z2OpdtPja309SoOLOZ1uXK+58llXbanp87JKC430ic2aM57aLx5KdkYbX76RHv+3isdafYUwfFe32lHTwONKySRENzV6eeX8PD7y6lZ2HndoVWRkumloiV+r1Q4dzL7qjpCCTlXe1Le0+Z8b4mIPE7Y+vbZ04mOYSZp06jHuvDs9iY4zpLdGChnbwONKy6WUer583Kw9x75JNvLPzCOB0OM88ZRjfvngsn5q/ssPX9jRg5LldrJt3WY/2Ae2Dl8+vgeW1FjiMSRHRgsZpInIM56oiO/CYwLLlnkoRXp+fjftquW/pZl5av781mk8eXcRtl4zjrPJisjIS1+kcr4AB7ScOgjNKa/H7+7j36g5eZIxJqmhpRGx4Swrz+5W9RxtZuKKSv67eRVOgzOrogTl8+6KxfOrjw8kNmZznwrkV1ROZabDxJ5f3cC8d62jiYLQJhcaY5IplnoZJMTUNHh57aye/W7mNQ/UeAIpyMrj+/HK+fM4oinPdhE/ArLzncirufLbbgSNdEhswoO3EwfD1xpjUYEGjD2nweHlh3T7uf3ULlQecMqvudBdfPLOEm6ZVUFqU0+4NNpj36ViTt9vHTRfYcndiAwbArFOH8Y9397abODjr1GEJP7YxJja9GjREZCYwH0gDFqnqPWHPZwKP4BR8OgR8SVW3J7udva2pxcfanTXMX7KZt7Ydbl0/42ND+NZFY5k4oqBNUsGg+5Zs4pdLNvf4+F6FMf/nuYSPaHL22b3RU5FGXQE8/e7e1n6eQXlufv6F0064krPLNlSzYEUlu2oaKI1D2d2e7C/ebTGOspCh80Hb70nMB71Oy70mioikAZuAS3ESIa4GrlHVD0O2+SZwqqreIiJXA1eq6pei7bc/lXtt8fnZWl3HA8u28Oz7x8usnlpSwHcuGsv54waR4+447kf6Q4qXK08fnjIjmroyZDjHncZvrj3jhHmjWrahmrmL15ORJmRnOLPzW3zKvFkTu/Uz6Mn+4t0W44j2f96VwBFruddok/sS7Wxgi6pWqqoHeJz26UmuAP4YePwkcImcANkSfX5ld00DP31hA5978A2eec8JGCMLs/nvKybyyPVnM+PkoVEDRiKIHB/ZFBzplApCR12FtjGSBo+PBSsqk9Sy3rdgRSUZaUKOOx0R53tGmnT7Z9CT/cW7LaZ39ObtqZHArpDlKuCcjrZRVa+IHAUGAgdDNxKRm4GbAUaNGpWo9iacqlLT4OGpt3fz29cqqa49Xmb1q1NG85VzRzNkQFZKdAyn0oimrralqqYhQS1JPbtqGijMzmizLjsjrds/g57sL95tMb2jN4NGpHe+8P/+WLZBVRcCC8G5PdXzpiVfbVMLyzYe4P5Xt7BxXy0A6S7hykkjufGCcsoH5UXst+gtqRC4gjoaddWRkqKcBLYmtZQW5VBd29TmqrSxxdftn0FP9hfvtpje0ZvvQlVAachyCbCno21EJB0oAA7TjzR6fLxZeZBv/nkt33nsndaAMW38IP50wzn81xUTOWlYfkoEDNXjKdFTaURTsC3B9kXrpstxpzF7WkWSWtb7Zk+roMWnNHi8qDrfW3za7Z9BT/YX77aY3tGbHeHpOB3hlwC7cTrCr1XV9SHbfAv4eEhH+OdU9apo++0rHeHNXh/bD9bz29cqefqdPXgDn5QnDBvAty4ay7Rxg8nPTm8336Kr4tUZHvw0n6r5oGz0VMeCI5aqahooiePoqe7sL95tMY54jJ6KtSO814IGgIh8CvgVzpDbh1X1JyIyD1ijqotFJAt4FJiEc4VxtapG7TVL9aDh9fmprm3mT2/u4NE3d1AbmD8xND+Tmy6oYNZpIxiYl5lSt3+MMf1frEGjV+dpqOpzwHNh6+aGPG7CqUne5/n9Tif3sx/sZcHySnYfaQScYkTXnjOKL587iuEF2WSmW/YWY0zqshnhCaaq1DZ7WbX1IPcv3coHu48CTgbaz5w6ghsuKGfM4Lw2eaKMMSZV2TtVAtU3e/lo7zEeWr6VJR8dr5p3bkUxt144htNLi+LSb2GMMcliQSMBmlp87DrcwO9f38aTa3fjCWSgHTM4l1unj2X6SYMpynFbv4Uxps+xoBFHHq+fA3VNPLGmij++sZ2ahhYABua6+cb5ZXx2UgmDB2SmxPBZY4zpDgsaceDzK4frm1nyUTUPLd/KjkOBMqvpLq46q5SvnDuK0qJcst3WyW2M6dssaPSA368cbWxh7c4afrNsK2/vqAGcaeyfnDiMGy8oZ9zQARSEpU4wxpi+yoJGNwRHRG3ZX8ei1yp5ft2+1glkZ4wq5NbpYzhjVBFFOW5c1m9hjOlHLGh0UX2zl901jfzlXzt4fPUumlqcTu5Rxc7s1ks+NoTiXOu3MMb0TxY0YtTU4uNgXTOL393D717fxqE6p8xqQXYG1503misnjWRofrb1Wxhj+jULGp3weP3UNHh4bfMBHlpeyZbqOgAy0oTPn1HCV6eMprQ4h/ws67cwxvR/FjQ64Auk/Vi/+ygPLt/Km5XHk+tePGEIN11QzknD8inMzrB+C2PMCcOCRhhVZ0TU9oP1/P6N7Tzz3p7WMqunjMjn1uljOKusmKJcNxlp1m9hjDmxWNAIUdvUwr6jTTyxZhd/fmsn9R4fACMKs7j5ggpmnDyUQXmZZGVYv4Ux5sRkQQOnENLBuiZeWr+fha9Vsv+YU2Y1LzOdr547is+fUcrQgkwGWL+FMeYEd0IHjWavj8P1HtZsP8yDy7by4V6nal6aS7ji9BF8fUoZo4pzKMzJsKSCxhjDCRo0vD4/hxs8bN5fy29XbGPZpgOtz00dO4ibp5UzYXg+xTlu0q3fwhhjWp1QQcPnV440eNhzpJFH39zBP97ZTYvP6eUePzSPW6eP4ezygQzMdVu/hTHGRHBCBA1V5Vijl4N1TfzjnT08smo7xwJlVocMyOSGqeV8cuIwBua5rd/CGGOi6PdBo7aphcN1HlZsPhAoaO+UWc3OSOPac0r54pmlDM3Psn4LY4yJQb8NGo0eH4fqm1m3+ygPLd/Ku7uOl1m9/OPD+fp5ZYwamGP9FsYY0wX9Lmiowr6jTew4VM/vVm7j5Q/3t2agPbusiNkXjuFjw/Mptn4LY4zpsn4XNFp8fu5/dTNPrKmiOVBmtXxQLrdcWMGUikEU57nJy+x3p22MMUnR7949tx6o409v7gSgKCeD688v5/KPD6c41239FsYY00P9Lmh4/UpmuosvTi7h6rNKGZKfZf0WxhgTJ70SNESkGPgrUAZsB65S1ZqwbU4HHgTyAR/wE1X9a2f7LsjO4JFvnE1pcQ4D89xkplu/hTHGxEtvffy+E3hFVccBrwSWwzUAX1PVicBM4FciUtjZjkcWZnPKyAJGFGZbwDDGmDjrraBxBfDHwOM/Ap8N30BVN6nq5sDjPUA1MLizHae5hFzr6DbGmIToraAxVFX3AgS+D4m2sYicDbiBrR08f7OIrBGRNQcOHIi0iTHGmDhI2EdyEVkCDIvw1H92cT/DgUeBr6uqP9I2qroQWAgwefJkjbSNMcaYnktY0FDVGR09JyL7RWS4qu4NBIXqDrbLB54FfqCqbyaoqcYYY2LUW7enFgNfDzz+OvDP8A1ExA38A3hEVZ9IYtuMMcZ0oLeCxj3ApSKyGbg0sIyITBaRRYFtrgKmAdeJyLuBr9N7p7nGGGMARLV/dQFMnjxZ16xZ09vNMMaYPkVE3lbVyZ1tZ9OkjTHGxMyChjHGmJhZ0DDGGBOzftenISIHgB1x3OUg4GAc99db7DxSi51HarHzgNGq2mnWjX4XNOJNRNbE0jmU6uw8UoudR2qx84id3Z4yxhgTMwsaxhhjYmZBo3MLe7sBcWLnkVrsPFKLnUeMrE/DGGNMzOxKwxhjTMwsaIQRkWIReVlENge+F0XY5nQRWSUi60XkfRH5Um+0NRIRmSkiG0Vki4i0q4goIpki8tfA82+JSFnyW9m5GM7jDhH5MPDzf0VERvdGOzvT2XmEbPcFEVERSckRPLGch4hcFfidrBeRvyS7jbGI4e9qlIi8KiLvBP62PtUb7YxGRB4WkWoRWdfB8yIi9wXO8X0ROSOuDVBV+wr5An4G3Bl4fCfw0wjbjAfGBR6PAPYChSnQ9jScQlUVOEWr3gNODtvmm8BDgcdXA3/t7XZ38zwuAnICj2/tq+cR2G4AsAJ4E5jc2+3u5u9jHPAOUBRYHtLb7e7meSwEbg08PhnY3tvtjnAe04AzgHUdPP8p4HlAgHOBt+J5fLvSaC9hpWiT4Gxgi6pWqqoHeBznfEKFnt+TwCUiIklsYyw6PQ9VfVVVGwKLbwIlSW5jLGL5fQD8N86HlaZkNq4LYjmPm4AHVLUGQFUj1sjpZbGchwL5gccFwJ4kti8mqroCOBxlkytwSkqoOnWICgN1i+LCgkZ7cS1Fm2QjgV0hy1WBdRG3UVUvcBQYmJTWxS6W8wh1A84nq1TT6XmIyCSgVFX/N5kN66JYfh/jgfEi8rqIvCkiM5PWutjFch4/Ar4iIlXAc8B3ktO0uOrq/0+XJKxyXypLZinaJIt0xRA+PC6WbXpbzG0Uka8Ak4ELE9qi7ol6HiLiAu4FrktWg7oplt9HOs4tquk4V32vicgpqnokwW3riljO4xrgD6r6CxGZAjwaOI9U+P+OVUL/x0/IoKH9txRtFVAaslxC+8vr4DZVIpKOcwke7VK3N8RyHojIDJxAf6GqNiepbV3R2XkMAE4BlgXuEA4DFovILFVNpaIwsf5dvamqLcA2EdmIE0RWJ6eJMYnlPG4AZgKo6ioRycLJ55SKt9s6EtP/T3fZ7an2+nIp2tXAOBEpD7TxapzzCRV6fl8Almqg9yyFdHoegds6C4BZKXr/HDo5D1U9qqqDVLVMVctw+mZSLWBAbH9XT+MMTkBEBuHcrqpMais7F8t57AQuARCRjwFZwIGktrLnFgNfC4yiOhc4GrzlHhe9PRIg1b5w7u+/AmwOfC8OrJ8MLAo8/grQArwb8nV6b7ddj4+c2ITTx/KfgXXzcN6MwPkneALYAvwLqOjtNnfzPJYA+0N+/ot7u83dOY+wbZeRgqOnYvx9CPBL4EPgA+Dq3m5zN8/jZOB1nJFV7wKf6O02RziHx3BGbLbgXFXcANwC3BLyu3ggcI4fxPtvymaEG2OMiZndnjLGGBMzCxrGGGNiZkHDGGNMzCxoGGOMiZkFDWOMMTGzoGF6REQGisi7ga99IrI7ZNkdx+PMEJGjgf1+JCL/GVifJyKPi8gHIrJORF4TkRwRSReRdrORRSRNRF7r4rG/ISLDQparRKSw52fV5hgvisiAOOwn4nmHbVMsIrd0Y98/FgjxFrgAAASaSURBVJHv9tbxTWo4IWeEm/hR1UPA6QAi8iOgTlV/HrpNICGiaM9TMbyqqp8VkTzgfRH5X2AWsFNVrw4cawLO+PWO2usDLujicb8BrAX2da/ZnVPVTyZq3xEU44zrfyiJx0yl45sesCsNkxAiMjbwyf8hnDfc0tBPoCJytYgsCjweKiJ/F5E1IvKvwCzWDqlqXWCfY4DhwO6Q5zaok8qio3a1fhIOXL28Ejj2RhF5JML2X8IJin8Nu3r6rhyvuTA+sG2eiPwhcA7viMhnIuxvpIisDOxrnYicF1hfJSKFIT+3h8WpS/GIiHxSRN4QkU0SqLcR/qlfRDaISEnYsfJFZKmIrA2089OBp+4BTgq04Z7AtncG2v2+iMwN2cfcwM/mZZy0IJF+pmPEqc2yGifhX9yOb1JQb89utK/+84XzhvG9wOOxgB84K7CcDhwJ2fZqjs+w/ytwbuBxGRHqBAAzgKcDjwfjpHs4CTgTJ83DGzgpxsdGOl7IflrXB/ZZgxN40nDSTJwb4TUrCZnxjzMLN1hzYQ7H65P8jMBMaKAIZ+ZxVti+vg98P/A4DcgL2Wdh4OfWgjMz2YUzK/m3gW0+DzwZePxj4Lsh+92Ak2Mo9PwygAGBx0OAzSG/m3dDXvsp4Dc4M4ldwAvAeTipxN8DsnFylG0LPWbI658Drg08vi1ex+/tv2f7ivxlt6dMIm1V1VgS1s3A+eQZXC4SkWxVbQzb7iIReQcnGP23qm4EEJEK4BOB/awRJ119rHmP3tRAXh4ReRcnaMWSgPLvge9v47zpEWjDZXK8IlwWMAoneAStBhaIkwjvaVV9L8K+t6jqh4E2fYiTMgWclBB3xXJSAQL8VESm4vzMSsXJCxXuE8BlOEWUAPJwckcNAp4K/B4aReSZDo4zBQheVT0K/Fecjv9GTGdpksqChkmk+pDHftqmbM4KeSzA2eoUxonmVVWNVBSrFngKeCrQf3IZTu6dWIRmx/UR+/9E8HWhrxHgs6raYW0VVV0qItOBy4E/i8jdqvrnKG3yhyz7Q47lpe3t5dCfZ9DXcK4QzlBVrzg1IiJtJ8CPVfV3bVaKfI/YUmprB9v16PgmNVmfhkkKdTrBa0RknDh1JK4MeXoJ8K3ggoicHut+RWRqcCSTiGQCHwN2xKfVrWpx0ph35kWc21XBtk0K30CcWub7VHUh8Aeg3TYx2o5zay5YCKw0wjYFQHXgDftSjhfiCT+fF4EbRCQ3sL+SwBXBCuBzIpIlTimATxPZm8BVgcdfjuPxTQqyKw2TTN/HuV+9EycbamZg/beAB0Xkepy/yVcJCSKdGBd4LTgfgp7BSWefFr9m83tgkYg04tzn78h/Ab8SkQ8CbdlC+3KilwB3iEgLUIeTMbk7nsCpMPcOTrbiSLfjHgWeEZE1OAMHgiWK9wcGHXwAPKuqdwZGnb0Z+DnW4vRR/EtE/oHTr7EdJ4hEMgfnqukOnJIBcTk+cLBbPxmTUJbl1hhjTMzs9pQxxpiYWdAwxhgTMwsaxhhjYmZBwxhjTMwsaBhjjImZBQ1jjDExs6BhjDEmZhY0jDHG/P8RDQAD4/2o8AiK1AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.regplot(y, x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "SpearmanrResult(correlation=0.8380276496703185, pvalue=0.0)"
      ]
     },
     "execution_count": 68,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stats.spearmanr(df['PSI'], df1['PSI'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.8442792843244691, 0.0)"
      ]
     },
     "execution_count": 69,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stats.pearsonr(df['PSI'], df1['PSI'])"
   ]
  }
 ],
 "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.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}