data-cleaning.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import re\n",
    "import numpy as np\n",
    "import networkx as nx\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from Levenshtein import distance as edit_distance\n",
    "from itertools import combinations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "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>country</th>\n",
       "      <th>school</th>\n",
       "      <th>postal_code</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Malaysia</td>\n",
       "      <td>POWIIS</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Philippines</td>\n",
       "      <td>Bangladesh International School, English section</td>\n",
       "      <td>21452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Saudi Arabia</td>\n",
       "      <td>Global International School</td>\n",
       "      <td>23213</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Malaysia</td>\n",
       "      <td>Kolej Tuanku Ja'afar</td>\n",
       "      <td>71700</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>China</td>\n",
       "      <td>the middle school attach to Guangxi university</td>\n",
       "      <td>530000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        country                                            school postal_code\n",
       "0      Malaysia                                            POWIIS           .\n",
       "1   Philippines  Bangladesh International School, English section       21452\n",
       "2  Saudi Arabia                       Global International School       23213\n",
       "3      Malaysia                              Kolej Tuanku Ja'afar       71700\n",
       "4         China    the middle school attach to Guangxi university      530000"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = pd.read_csv('registration-trimmed.csv', encoding = \"ISO-8859-1\")    \n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(175, 3)"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Handle Missing Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "country        0.005714\n",
       "school         0.011429\n",
       "postal_code    0.011429\n",
       "dtype: float64"
      ]
     },
     "execution_count": 43,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Missing Data\n",
    "df.isna().mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "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>country</th>\n",
       "      <th>school</th>\n",
       "      <th>postal_code</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>169</th>\n",
       "      <td>Malaysia</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      country school postal_code\n",
       "56        NaN    NaN         NaN\n",
       "169  Malaysia    NaN         NaN"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df[df['school'].isna()]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "country        0.0\n",
       "school         0.0\n",
       "postal_code    0.0\n",
       "dtype: float64"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = df[~df['school'].isna()]\n",
    "df.isna().mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Convert to string values, standardize format"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "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>country</th>\n",
       "      <th>school</th>\n",
       "      <th>postal_code</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>powiis</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>philippines</td>\n",
       "      <td>bangladesh international school, english section</td>\n",
       "      <td>21452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>saudi arabia</td>\n",
       "      <td>global international school</td>\n",
       "      <td>23213</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>kolej tuanku ja'afar</td>\n",
       "      <td>71700</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>china</td>\n",
       "      <td>the middle school attach to guangxi university</td>\n",
       "      <td>530000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        country                                            school postal_code\n",
       "0      malaysia                                            powiis           .\n",
       "1   philippines  bangladesh international school, english section       21452\n",
       "2  saudi arabia                       global international school       23213\n",
       "3      malaysia                              kolej tuanku ja'afar       71700\n",
       "4         china    the middle school attach to guangxi university      530000"
      ]
     },
     "execution_count": 46,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "for col in df.columns:\n",
    "    df[col] = df[col].astype(str).str.strip()\n",
    "    df[col] = df[col].str.lower()\n",
    "    df[col] = df[col].apply(lambda text: re.sub('\\s+', ' ', text))\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(173, 3)"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Inspect Country"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "brunei darussalam         2\n",
       "china                     2\n",
       "egypt                     1\n",
       "hong kong (s.a.r.)        1\n",
       "indonesia                 1\n",
       "iraq                      1\n",
       "japan                     2\n",
       "kuwait                    1\n",
       "malaysia                148\n",
       "maldives                  2\n",
       "oman                      1\n",
       "philippines               1\n",
       "qatar                     1\n",
       "saudi arabia              1\n",
       "singapore                 6\n",
       "sudan                     1\n",
       "united arab emirates      1\n",
       "Name: country, dtype: int64"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['country'].value_counts().sort_index()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Inspect Postal Code"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1      12\n",
       "3       1\n",
       "4       2\n",
       "5     140\n",
       "6       9\n",
       "8       2\n",
       "9       1\n",
       "10      1\n",
       "11      1\n",
       "12      2\n",
       "17      1\n",
       "45      1\n",
       "Name: postal_code, dtype: int64"
      ]
     },
     "execution_count": 49,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['postal_code'].str.len().value_counts().sort_index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "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>country</th>\n",
       "      <th>school</th>\n",
       "      <th>postal_code</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>powiis</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>china</td>\n",
       "      <td>the middle school attach to guangxi university</td>\n",
       "      <td>530000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>qatar</td>\n",
       "      <td>the international school of choueifat</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>oman</td>\n",
       "      <td>the sultan's school</td>\n",
       "      <td>don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>singapore</td>\n",
       "      <td>global indian international school</td>\n",
       "      <td>459663</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>china</td>\n",
       "      <td>ningbo binhai internationalcooperation school</td>\n",
       "      <td>315800</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's international school kuala lumpur</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>singapore</td>\n",
       "      <td>nan chiau high school</td>\n",
       "      <td>545079</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>singapore</td>\n",
       "      <td>nan chiau high school</td>\n",
       "      <td>545079</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>singapore</td>\n",
       "      <td>nan chiau high school</td>\n",
       "      <td>545079</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>iraq</td>\n",
       "      <td>almtamezen-harthea</td>\n",
       "      <td>non</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's international school kuala lumpur</td>\n",
       "      <td>i don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's college</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>brunei darussalam</td>\n",
       "      <td>chung hua middle school</td>\n",
       "      <td>ka1131</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>kuen cheng high school</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>seri omega</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>64</th>\n",
       "      <td>japan</td>\n",
       "      <td>international budo university</td>\n",
       "      <td>299-5295</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sri kuala lumpur</td>\n",
       "      <td>i dont know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>japan</td>\n",
       "      <td>???????????????????????????????????????</td>\n",
       "      <td>169-0075</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>hong kong (s.a.r.)</td>\n",
       "      <td>homeschooled</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>singapore</td>\n",
       "      <td>evergreen secondary school</td>\n",
       "      <td>738489</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>singapore</td>\n",
       "      <td>evergreen secondary school</td>\n",
       "      <td>738489</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk tengku mahmud 2</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>brunei darussalam</td>\n",
       "      <td>chung hua middle school, kuala belait</td>\n",
       "      <td>[they changed it so i can't remember anymore]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>inti</td>\n",
       "      <td>dont know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>st.johns international</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smjk sin min</td>\n",
       "      <td>8000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sekolah menengah kebangsaan ibrahim</td>\n",
       "      <td>8000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>142</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk tengku ampuan rahimah</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>143</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smjk sam tet</td>\n",
       "      <td>i don??t know it</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>145</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk bandar utama damansara 3</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>146</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sri emas international school</td>\n",
       "      <td>i don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>148</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sekolah menengah stella maris</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                country                                          school  \\\n",
       "0              malaysia                                          powiis   \n",
       "4                 china  the middle school attach to guangxi university   \n",
       "6                 qatar           the international school of choueifat   \n",
       "9                  oman                             the sultan's school   \n",
       "10            singapore              global indian international school   \n",
       "12                china   ningbo binhai internationalcooperation school   \n",
       "20             malaysia      taylor's international school kuala lumpur   \n",
       "21            singapore                           nan chiau high school   \n",
       "22            singapore                           nan chiau high school   \n",
       "23            singapore                           nan chiau high school   \n",
       "30                 iraq                              almtamezen-harthea   \n",
       "33             malaysia      taylor's international school kuala lumpur   \n",
       "45             malaysia                                taylor's college   \n",
       "51    brunei darussalam                         chung hua middle school   \n",
       "58             malaysia                          kuen cheng high school   \n",
       "61             malaysia                                      seri omega   \n",
       "64                japan                   international budo university   \n",
       "65             malaysia                                sri kuala lumpur   \n",
       "69                japan         ???????????????????????????????????????   \n",
       "71   hong kong (s.a.r.)                                    homeschooled   \n",
       "75            singapore                      evergreen secondary school   \n",
       "76            singapore                      evergreen secondary school   \n",
       "77             malaysia                             smk tengku mahmud 2   \n",
       "78    brunei darussalam           chung hua middle school, kuala belait   \n",
       "82             malaysia                                            inti   \n",
       "83             malaysia                          st.johns international   \n",
       "104            malaysia                                    smjk sin min   \n",
       "125            malaysia             sekolah menengah kebangsaan ibrahim   \n",
       "142            malaysia                       smk tengku ampuan rahimah   \n",
       "143            malaysia                                    smjk sam tet   \n",
       "145            malaysia                    smk bandar utama damansara 3   \n",
       "146            malaysia                   sri emas international school   \n",
       "148            malaysia                   sekolah menengah stella maris   \n",
       "\n",
       "                                       postal_code  \n",
       "0                                                .  \n",
       "4                                           530000  \n",
       "6                                                0  \n",
       "9                                       don't know  \n",
       "10                                          459663  \n",
       "12                                          315800  \n",
       "20                                               -  \n",
       "21                                          545079  \n",
       "22                                          545079  \n",
       "23                                          545079  \n",
       "30                                             non  \n",
       "33                                    i don't know  \n",
       "45                                               .  \n",
       "51                                          ka1131  \n",
       "58                                               -  \n",
       "61                                               -  \n",
       "64                                        299-5295  \n",
       "65                                     i dont know  \n",
       "69                                        169-0075  \n",
       "71                                               0  \n",
       "75                                          738489  \n",
       "76                                          738489  \n",
       "77                                               -  \n",
       "78   [they changed it so i can't remember anymore]  \n",
       "82                                       dont know  \n",
       "83                                               -  \n",
       "104                                           8000  \n",
       "125                                           8000  \n",
       "142                                              -  \n",
       "143                              i don??t know it  \n",
       "145                                              -  \n",
       "146                                   i don't know  \n",
       "148                                              -  "
      ]
     },
     "execution_count": 50,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df[df['postal_code'].str.len() != 5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "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>country</th>\n",
       "      <th>school</th>\n",
       "      <th>postal_code</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>powiis</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>qatar</td>\n",
       "      <td>the international school of choueifat</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>oman</td>\n",
       "      <td>the sultan's school</td>\n",
       "      <td>don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's international school kuala lumpur</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>iraq</td>\n",
       "      <td>almtamezen-harthea</td>\n",
       "      <td>non</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's international school kuala lumpur</td>\n",
       "      <td>i don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>taylor's college</td>\n",
       "      <td>.</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>brunei darussalam</td>\n",
       "      <td>chung hua middle school</td>\n",
       "      <td>ka1131</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>kuen cheng high school</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>seri omega</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sri kuala lumpur</td>\n",
       "      <td>i dont know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>hong kong (s.a.r.)</td>\n",
       "      <td>homeschooled</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk tengku mahmud 2</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>brunei darussalam</td>\n",
       "      <td>chung hua middle school, kuala belait</td>\n",
       "      <td>[they changed it so i can't remember anymore]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>inti</td>\n",
       "      <td>dont know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>st.johns international</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smjk sin min</td>\n",
       "      <td>8000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sekolah menengah kebangsaan ibrahim</td>\n",
       "      <td>8000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>142</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk tengku ampuan rahimah</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>143</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smjk sam tet</td>\n",
       "      <td>i don??t know it</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>145</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>smk bandar utama damansara 3</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>146</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sri emas international school</td>\n",
       "      <td>i don't know</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>148</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>sekolah menengah stella maris</td>\n",
       "      <td>-</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                country                                      school  \\\n",
       "0              malaysia                                      powiis   \n",
       "6                 qatar       the international school of choueifat   \n",
       "9                  oman                         the sultan's school   \n",
       "20             malaysia  taylor's international school kuala lumpur   \n",
       "30                 iraq                          almtamezen-harthea   \n",
       "33             malaysia  taylor's international school kuala lumpur   \n",
       "45             malaysia                            taylor's college   \n",
       "51    brunei darussalam                     chung hua middle school   \n",
       "58             malaysia                      kuen cheng high school   \n",
       "61             malaysia                                  seri omega   \n",
       "65             malaysia                            sri kuala lumpur   \n",
       "71   hong kong (s.a.r.)                                homeschooled   \n",
       "77             malaysia                         smk tengku mahmud 2   \n",
       "78    brunei darussalam       chung hua middle school, kuala belait   \n",
       "82             malaysia                                        inti   \n",
       "83             malaysia                      st.johns international   \n",
       "104            malaysia                                smjk sin min   \n",
       "125            malaysia         sekolah menengah kebangsaan ibrahim   \n",
       "142            malaysia                   smk tengku ampuan rahimah   \n",
       "143            malaysia                                smjk sam tet   \n",
       "145            malaysia                smk bandar utama damansara 3   \n",
       "146            malaysia               sri emas international school   \n",
       "148            malaysia               sekolah menengah stella maris   \n",
       "\n",
       "                                       postal_code  \n",
       "0                                                .  \n",
       "6                                                0  \n",
       "9                                       don't know  \n",
       "20                                               -  \n",
       "30                                             non  \n",
       "33                                    i don't know  \n",
       "45                                               .  \n",
       "51                                          ka1131  \n",
       "58                                               -  \n",
       "61                                               -  \n",
       "65                                     i dont know  \n",
       "71                                               0  \n",
       "77                                               -  \n",
       "78   [they changed it so i can't remember anymore]  \n",
       "82                                       dont know  \n",
       "83                                               -  \n",
       "104                                           8000  \n",
       "125                                           8000  \n",
       "142                                              -  \n",
       "143                              i don??t know it  \n",
       "145                                              -  \n",
       "146                                   i don't know  \n",
       "148                                              -  "
      ]
     },
     "execution_count": 51,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#The only acceptable format are \n",
    "# empty\n",
    "# 5 digits \n",
    "# 6 digits\n",
    "# 3 digits follows by hypen and 4 digits\n",
    "def is_valid_postal_code(t):\n",
    "    return re.match('^\\d{5}$', t) is not None or \\\n",
    "           re.match('^\\d{6}$', t) is not None or \\\n",
    "           re.match('^\\d{3}-\\d{4}$', t) is not None or \\\n",
    "           t == ''\n",
    "\n",
    "df[~df['postal_code'].apply(is_valid_postal_code)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.0"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Fix invalid postal code, convert to missing values\n",
    "df['postal_code'].loc[~df['postal_code'].apply(is_valid_postal_code)] = ''\n",
    "df['postal_code'].apply(is_valid_postal_code).mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0     23\n",
       "5    140\n",
       "6      8\n",
       "8      2\n",
       "Name: postal_code, dtype: int64"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Check again\n",
    "df['postal_code'].str.len().value_counts().sort_index()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Identify school groups"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "smjk yu hua                              7\n",
       "smk gajah berang                         4\n",
       "catholic high school                     4\n",
       "nan chiau high school                    3\n",
       "smjk jit sin                             3\n",
       "                                        ..\n",
       "smk taman setiawangsa                    1\n",
       "smk tun habab                            1\n",
       "seri omega                               1\n",
       "chong hwa independent high school        1\n",
       "center for higher secondary education    1\n",
       "Name: school, Length: 129, dtype: int64"
      ]
     },
     "execution_count": 54,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['school'].value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "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>country</th>\n",
       "      <th>postal_code</th>\n",
       "      <th>schools</th>\n",
       "      <th>n_schools</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>brunei darussalam</td>\n",
       "      <td></td>\n",
       "      <td>{chung hua middle school, chung hua middle sch...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>malaysia</td>\n",
       "      <td></td>\n",
       "      <td>{smk tengku mahmud 2, sri emas international s...</td>\n",
       "      <td>16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>11600</td>\n",
       "      <td>{penang free school, han chiang high school}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>14000</td>\n",
       "      <td>{jit sin, smjk jit sin, high school bukit mert...</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>26600</td>\n",
       "      <td>{sekolah sains sultan haji ahmad shah, pekan, ...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>41200</td>\n",
       "      <td>{smjk chung hwa, smk bukit tinggi, smk tengku ...</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>43000</td>\n",
       "      <td>{smjk yu hua kajang, smk bandar baru sungai lo...</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>46000</td>\n",
       "      <td>{catholic high school, smjk katholik pj}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>47100</td>\n",
       "      <td>{help international school, smk bandar puchong...</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>47400</td>\n",
       "      <td>{smk damansara utama, smk damansara jaya}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>47500</td>\n",
       "      <td>{sri kuala lumpur secondary school, sri kuala ...</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>47800</td>\n",
       "      <td>{nobel international school, bandar utama dama...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>52000</td>\n",
       "      <td>{smk jinjang, chong hwa independent high school}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>44</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>52100</td>\n",
       "      <td>{smk kepong baru, sekolah sri bestari, smk kep...</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>56000</td>\n",
       "      <td>{smk seri mutiara, smk taman connaught, oxburg...</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>56100</td>\n",
       "      <td>{smk seri bintang seletan, sri sempurna intern...</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>75200</td>\n",
       "      <td>{smk gajah berang, catholic high school}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>75250</td>\n",
       "      <td>{smk tun haji abdul malek, oakrich internation...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>81400</td>\n",
       "      <td>{seri omega private and international school, ...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>malaysia</td>\n",
       "      <td>98000</td>\n",
       "      <td>{smk chung hua (cf) miri, smk lutong}</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>maldives</td>\n",
       "      <td>20159</td>\n",
       "      <td>{centre for higher secondary education (chse),...</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              country postal_code  \\\n",
       "0   brunei darussalam               \n",
       "10           malaysia               \n",
       "12           malaysia       11600   \n",
       "16           malaysia       14000   \n",
       "18           malaysia       26600   \n",
       "25           malaysia       41200   \n",
       "26           malaysia       43000   \n",
       "29           malaysia       46000   \n",
       "31           malaysia       47100   \n",
       "35           malaysia       47400   \n",
       "36           malaysia       47500   \n",
       "38           malaysia       47800   \n",
       "43           malaysia       52000   \n",
       "44           malaysia       52100   \n",
       "50           malaysia       56000   \n",
       "51           malaysia       56100   \n",
       "60           malaysia       75200   \n",
       "61           malaysia       75250   \n",
       "65           malaysia       81400   \n",
       "78           malaysia       98000   \n",
       "79           maldives       20159   \n",
       "\n",
       "                                              schools  n_schools  \n",
       "0   {chung hua middle school, chung hua middle sch...          2  \n",
       "10  {smk tengku mahmud 2, sri emas international s...         16  \n",
       "12       {penang free school, han chiang high school}          2  \n",
       "16  {jit sin, smjk jit sin, high school bukit mert...          3  \n",
       "18  {sekolah sains sultan haji ahmad shah, pekan, ...          2  \n",
       "25  {smjk chung hwa, smk bukit tinggi, smk tengku ...          4  \n",
       "26  {smjk yu hua kajang, smk bandar baru sungai lo...          4  \n",
       "29           {catholic high school, smjk katholik pj}          2  \n",
       "31  {help international school, smk bandar puchong...          5  \n",
       "35          {smk damansara utama, smk damansara jaya}          2  \n",
       "36  {sri kuala lumpur secondary school, sri kuala ...          4  \n",
       "38  {nobel international school, bandar utama dama...          2  \n",
       "43   {smk jinjang, chong hwa independent high school}          2  \n",
       "44  {smk kepong baru, sekolah sri bestari, smk kep...          3  \n",
       "50  {smk seri mutiara, smk taman connaught, oxburg...          3  \n",
       "51  {smk seri bintang seletan, sri sempurna intern...          5  \n",
       "60           {smk gajah berang, catholic high school}          2  \n",
       "61  {smk tun haji abdul malek, oakrich internation...          2  \n",
       "65  {seri omega private and international school, ...          2  \n",
       "78              {smk chung hua (cf) miri, smk lutong}          2  \n",
       "79  {centre for higher secondary education (chse),...          2  "
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "school_grps = df.groupby(['country', 'postal_code']) \\\n",
    "                  .apply(lambda x: set(x['school'])).reset_index() \\\n",
    "                  .rename(columns={0:'schools'})\n",
    "school_grps['n_schools'] = school_grps['schools'].apply(len)\n",
    "school_grps = school_grps[school_grps['n_schools'] > 1]\n",
    "school_grps"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create School Mapper "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [],
   "source": [
    "#################\n",
    "# Configuration\n",
    "################\n",
    "EDIT_DISTANCE_FACTOR = 1/3\n",
    "MAX_WORD_DIFFERENCES = 3\n",
    "\n",
    "def get_school_mapper(df, country, postal_code):\n",
    "    schools_same_area = set(df.query(\"country==@country and postal_code==@postal_code\")['school'])\n",
    "\n",
    "    G = nx.Graph()\n",
    "\n",
    "    for s1, s2 in combinations(schools_same_area, 2):\n",
    "        ######################\n",
    "        # partial matching\n",
    "        ######################\n",
    "        if s1 in s2 or s2 in s1:\n",
    "            G.add_edge(s1, s2)\n",
    "\n",
    "        ##############################\n",
    "        # Check character differences\n",
    "        ##############################\n",
    "        s1_num_words = len(s1.split(' '))\n",
    "        s2_num_words = len(s2.split(' '))\n",
    "        diff_in_words = abs(s1_num_words - s2_num_words)\n",
    "        if edit_distance(s1, s2) <= (min(len(s1), len(s2)) * EDIT_DISTANCE_FACTOR) and \\\n",
    "           diff_in_words <= MAX_WORD_DIFFERENCES:\n",
    "            G.add_edge(s1, s2)\n",
    "\n",
    "    ###########################\n",
    "    # Create homogenous groups\n",
    "    ###########################\n",
    "    similar_schools = list(nx.connected_components(G))\n",
    "\n",
    "    ################\n",
    "    # Create lookup\n",
    "    ################\n",
    "    school_mapper = {}\n",
    "    for schools in similar_schools:\n",
    "        ######################################\n",
    "        # Set main school as the shortest one\n",
    "        ######################################\n",
    "        main_school = min(schools, key=lambda x: len(x))\n",
    "        for school in schools:\n",
    "            school_mapper[school] = main_school\n",
    "\n",
    "    return school_mapper"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'sri kuala lumpur': 'sri kuala lumpur',\n",
       " 'sri kuala lumpur international school': 'sri kuala lumpur',\n",
       " 'sri kuala lumpur secondary school': 'sri kuala lumpur'}"
      ]
     },
     "execution_count": 57,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "get_school_mapper(df, 'malaysia', '47500')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'smk seri bintang selatan': 'smk seri bintang utara',\n",
       " 'smk seri bintang seletan': 'smk seri bintang utara',\n",
       " 'smk seri bintang utara': 'smk seri bintang utara',\n",
       " 'sri sempurna': 'sri sempurna',\n",
       " 'sri sempurna international school': 'sri sempurna'}"
      ]
     },
     "execution_count": 58,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "get_school_mapper(df, 'malaysia', '56100')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Inspect all mappers (Modify manually if necessary)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{('brunei darussalam',\n",
       "  ''): {'chung hua middle school': 'chung hua middle school', 'chung hua middle school, kuala belait': 'chung hua middle school'},\n",
       " ('malaysia', '14000'): {'jit sin': 'jit sin', 'smjk jit sin': 'jit sin'},\n",
       " ('malaysia',\n",
       "  '26600'): {'sekolah sains sultan haji ahmad shah': 'sekolah sains sultan haji ahmad shah', 'sekolah sains sultan haji ahmad shah, pekan': 'sekolah sains sultan haji ahmad shah'},\n",
       " ('malaysia', '43000'): {'smjk yu hua': 'smjk yu hua',\n",
       "  'smjk yu hua kajang': 'smjk yu hua'},\n",
       " ('malaysia',\n",
       "  '47100'): {'axcel international school': 'help international school', 'help international school': 'help international school'},\n",
       " ('malaysia', '47400'): {'smk damansara jaya': 'smk damansara jaya',\n",
       "  'smk damansara utama': 'smk damansara jaya'},\n",
       " ('malaysia', '47500'): {'sri kuala lumpur': 'sri kuala lumpur',\n",
       "  'sri kuala lumpur international school': 'sri kuala lumpur',\n",
       "  'sri kuala lumpur secondary school': 'sri kuala lumpur'},\n",
       " ('malaysia', '52100'): {'smk kepong': 'smk kepong',\n",
       "  'smk kepong baru': 'smk kepong'},\n",
       " ('malaysia', '56100'): {'smk seri bintang selatan': 'smk seri bintang utara',\n",
       "  'smk seri bintang seletan': 'smk seri bintang utara',\n",
       "  'smk seri bintang utara': 'smk seri bintang utara',\n",
       "  'sri sempurna': 'sri sempurna',\n",
       "  'sri sempurna international school': 'sri sempurna'},\n",
       " ('maldives',\n",
       "  '20159'): {'center for higher secondary education': 'center for higher secondary education', 'centre for higher secondary education (chse)': 'center for higher secondary education'}}"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "school_grps = df.groupby(['country', 'postal_code']) \\\n",
    "                  .apply(lambda x: set(x['school'])).reset_index() \\\n",
    "                  .rename(columns={0:'schools'})\n",
    "school_grps['n_schools'] = school_grps['schools'].apply(len)\n",
    "school_grps = school_grps[school_grps['n_schools'] > 1]\n",
    "\n",
    "school_mappers = {}\n",
    "for row in school_grps[['country', 'postal_code']].itertuples(index=False):\n",
    "    country, postal_code = row[0], row[1]\n",
    "    mapper = get_school_mapper(df, country, postal_code)\n",
    "    if len(mapper) > 0:\n",
    "        school_mappers[(country, postal_code)] = mapper\n",
    "    \n",
    "school_mappers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "# modify if necessary:\n",
    "school_mappers = {('brunei darussalam',\n",
    "  ''): {'chung hua middle school, kuala belait': 'chung hua middle school'},\n",
    " ('malaysia', '14000'): {'smjk jit sin': 'jit sin'},\n",
    " ('malaysia',\n",
    "  '26600'): {'sekolah sains sultan haji ahmad shah, pekan': 'sekolah sains sultan haji ahmad shah'},\n",
    " ('malaysia', '43000'): {'smjk yu hua kajang': 'smjk yu hua'},\n",
    " ('malaysia',\n",
    "  '47100'): {'axcel international school': 'help international school'},\n",
    " ('malaysia', '47400'): {'smk damansara utama': 'smk damansara jaya'},\n",
    " ('malaysia',\n",
    "  '47500'): {'sri kuala lumpur international school': 'sri kuala lumpur',\n",
    "  'sri kuala lumpur secondary school': 'sri kuala lumpur'},\n",
    " ('malaysia', '52100'): {'smk kepong baru': 'smk kepong'},\n",
    " ('malaysia', '56100'): {'smk seri bintang selatan': 'smk seri bintang utara',\n",
    "  'smk seri bintang seletan': 'smk seri bintang utara',\n",
    "  'sri sempurna international school': 'sri sempurna'},\n",
    " ('maldives',\n",
    "  '20159'): {'centre for higher secondary education (chse)': 'center for higher secondary education'}}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Fix school names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Changed row 78-th school: \"chung hua middle school, kuala belait\" -> \"chung hua middle school\"\n",
      "Changed row 13-th school: \"smjk jit sin\" -> \"jit sin\"\n",
      "Changed row 141-th school: \"smjk jit sin\" -> \"jit sin\"\n",
      "Changed row 171-th school: \"sekolah sains sultan haji ahmad shah, pekan\" -> \"sekolah sains sultan haji ahmad shah\"\n",
      "Changed row 113-th school: \"smjk yu hua kajang\" -> \"smjk yu hua\"\n",
      "Changed row 41-th school: \"axcel international school\" -> \"help international school\"\n",
      "Changed row 120-th school: \"smk damansara utama\" -> \"smk damansara jaya\"\n",
      "Changed row 42-th school: \"sri kuala lumpur international school\" -> \"sri kuala lumpur\"\n",
      "Changed row 67-th school: \"sri kuala lumpur secondary school\" -> \"sri kuala lumpur\"\n",
      "Changed row 129-th school: \"smk kepong baru\" -> \"smk kepong\"\n",
      "Changed row 57-th school: \"smk seri bintang selatan\" -> \"smk seri bintang utara\"\n",
      "Changed row 62-th school: \"sri sempurna international school\" -> \"sri sempurna\"\n",
      "Changed row 84-th school: \"sri sempurna international school\" -> \"sri sempurna\"\n",
      "Changed row 132-th school: \"smk seri bintang seletan\" -> \"smk seri bintang utara\"\n",
      "Changed row 133-th school: \"smk seri bintang selatan\" -> \"smk seri bintang utara\"\n",
      "Changed row 16-th school: \"centre for higher secondary education (chse)\" -> \"center for higher secondary education\"\n"
     ]
    }
   ],
   "source": [
    "df_clean = df.copy()\n",
    "\n",
    "def change_school(row, mapper):\n",
    "    school = row['school']\n",
    "    index = row.name\n",
    "    if school in mapper and school != mapper[school]:\n",
    "        print(f'Changed row {index}-th school: \"{school}\" -> \"{mapper[school]}\"')\n",
    "        return mapper[school]\n",
    "    else:\n",
    "        return school\n",
    "\n",
    "for (country, postal_code), school_mapper in school_mappers.items():\n",
    "    df_qry = df_clean.query(\"country==@country and postal_code==@postal_code\")\n",
    "    new_schools = df_qry.apply(lambda row: change_school(row, school_mapper), axis=1)\n",
    "    df_clean['school'].loc[df_qry.index] = new_schools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "smjk yu hua                          8\n",
       "smk seri bintang utara               4\n",
       "sri kuala lumpur                     4\n",
       "catholic high school                 4\n",
       "smk gajah berang                     4\n",
       "                                    ..\n",
       "smk tun habab                        1\n",
       "smk bukit mewah                      1\n",
       "chong hwa independent high school    1\n",
       "smjk jit sin                         1\n",
       "wesley methodist ipoh                1\n",
       "Name: school, Length: 117, dtype: int64"
      ]
     },
     "execution_count": 62,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_clean['school'].value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyEAAAJOCAYAAACgOwjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3XmUZWV57/Hvj24BAQUVAk2itAoOzENDZFBRiddIO6OdBKWJUYK5aib0EgdEVLQDN3YQJ/QaVHBi0IuQCAZBRcYGegIVI9AxoihXISCD0P3cP85beihPTT3squr+ftaqdc55h+d99q6qXuepd+/TqSokSZIkqSsbTXYCkiRJkjYsFiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSxi3JF5K8Yx3E/UCST67tuJKmJosQSZLWQJJ7+r5WJbmv7/Xha3mtU5L8MMndSW5M8qfD+vdNsjjJvUmuTrLrKLFmJ/lKkjuS3JVkaZI/W5v5StJILEIkSVoDVbXF0Bfwn8CL+trOXMvL/Tfwx8CWwFHAx5LsA5DkkcD/BU4DHgOcBXw5ycwRYn0e+D7weGBr4M+BO9ZyvpI0kEWIJEnrUJJHJvlwkp8k+a8kJyV5ROt7QZL/SPLuJL9IcnOSV44Uq6reUVU3VdWqqroMuAp4Ruv+I+D+qvpIVT0A/G/gUcBBA3IKMAf4l6q6r6oerKprq+qivjEHJ7my7ZL857Bdkq2TXNh2ZL6TZIe+ec9Ocl2bd2WSffv6npDkX9ux3pRk/gjnbPN22dcvktyZ5Kokjxn7bEuaLixCJElat94N7A7sBuwDHAy8ta9/NrAxsB293Y1PJ3niWEGTbAHsDdzQmnYBlgz1V9UqYHlrf5iqKnoFzMeTvCrJHwyLvSNwPnAS8LiW9w19Q/4M+AfgscBP2jGS5PeArwIfaPM+Bvxrki3bvLPo7b7MajE+mOTAAYf3OmAm8Pv0dmneCPx6rHMiafqwCJEkad06HHhXVd1RVbcD7wVe09f/EPDuqvp1Vf078O/AYaMFbDsZnwQuq6pLW/MWwF3Dht5FbzdkkJcCV9MrIP4zyaIke7W+1wBfrapzquqhqvp5VS3pm/ulqrquqh4EPgfs2dpfAiyuqi+1eacD/wX8cZKdgD2At1XVA1W1CPj0sHMx5EFgG+DJLc41VfWr0c6JpOnFIkSSpHWkFQvbASv6mlfQ+wv/kJ9X1f3D+rcfI/QpwA7Aq/va7gEePWzco4G7BwVoRdFbqurpLcebgHNb9+OBH46y/k/7nt9LrwCi5b1i2Nih492e3rHeN6BvuP8DfBM4u13CdmKSGaPkI2masQiRJGkdaZc9/ZRewTDkCcCP+15vnWTTYf23jRQzyQfo3efxx1V1T1/XDfR2GobGbQTsysMvoxopz58B/wTMTrI58CPgyWPNG+A2Hn6s8NvjvQ3Ypt1AP7xveD4PVNVxVfU04FnAK4E/WY18JE1RFiGSJK1bnwfeleRx7Z6JtwNn9PU/Anhnko2TPJfeDebnDAqU5N30Lnl6flXdOaz768AjkxydZBPgb4FfAZeNEOvkJDsnmdHu2TgaWN4ue/oMMDfJy1r/Nkl2H8exngfsleSwJDOTHEGv0Pga8B/AUuC9STZJsjcwH/idTxBLckjLbSN6nwj2ELByHOtLmiYsQiRJWreOA26ktyOxGPgO8I99/bfSe5P9U+BTwJ9X1c3Dg7TC4jjgScAtff8Xyd8BtMucXkKvmLiT3s7BS6vqoRHyejS9ouEuegXCNsDLW6wftlhvA34JLGLADe7DtXteXkyv0Pp/9G4on1tVd7ZdoVcBO7dj/SLwlqr69oBQv0/v44bvpndz/b8CXxprfUnTR3r/JkiSpK4leQFwalXtONm5SFKX3AmRJEmS1CmLEEmSJEmd8nIsSZIkSZ1yJ0SSJElSp2ZOdgKSRpfE7UpJkjRd3FFV24w1yJ0QSZIkSWvLivEMsgiRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEEmSJEmdsgiRJEmS1Ck/olea4jbebkdmzV842WlIkqRpasWCuZOdwu9wJ0SSJElSpyxCJEmSJHXKIkSSJElSpyxCtF5Jsn2Ss9vzI5Ocug7Xmp1k+bqKL0mStL7yxnStV6rqNuCwyc5DkiRJI3MnRFNOks2TXJBkSZLlSea19luTnJjkiiSLkuyd5MIkP0xydBszcHciyaFt3tZ9bRsl+UGSbfpe/0eSrZOcnuSwvrH3jJDujCSfSHJDkouSPLKNvzTJnPZ86yS39uX37STXta8D1s5ZkyRJmj4sQjQVvQC4rar2qKpdga/19f2oqvYHvg2cTm/X4xnACSMFS/Iy4FjghVV1x1B7Va0CzgAOb02HAEv6x4zDTsCHq2oX4E7gFWOM/xnwR1W1NzAPOGWEnI9qhdailffeNYF0JEmSpj6LEE1Fy4BDkixI8syq6n8Xfl7fmKuq6u6q+jlwf5KtBsR6DvC/gEOr6pcD+j8FHNGevxb4lwnmektVLW7PrwVmjzH+EcAnkiwDzgJ2HjSoqk6rqjlVNWfGZltOMCVJkqSpzSJEU05V3QTsQ6/QeH+S4/q6H2iPq/qeD70edI/TzcCjgKeMsNaPgNuTPBf4Q+DfWtdDtN+PJAE2HiHd/hxW9uXwm/nApn1j/ha4HdgDmDNKXEmSpPWWRYimnCTbA/dW1RnAycDeaxBuBfBy4DNJdhlhzCfpXZb1papa2dpupVcIAbyE3g7GRPTP779RfkvgJ+1SsNcAMyYYV5IkadqzCNFUtBtwdZLFwNuB905wfj3sRdX36d33cVaSJw8Yfx6wBQ+/FOsTwLOTXE1vh+RXE8zhZOANSS4Htu5r/wgwP8mV9HZnJhpXkiRp2ktVjT1KmiaS7AP8U1U9ewJz5gAfrKpnrrvMVt8ms3aqWfMXTnYakiRpmlqxYG6Xy11bVXPGGuROiNYbrZj4PPDPE5hzLHAO8A/rKi9JkiQ9nDsh0hSXxF9SSZI0XbgTIkmSJGnqsQiRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEEmSJEmdsgiRJEmS1KmZk52ApNFtvN2OzJq/cLLTkDRFrVgwd7JTkKQJcydEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQjZgSQ5Ocv46XuOTSXYe0H5rkq0HtB+d5IgxYu6Z5IVrM8+JSLJVkr+arPUlSZKmO4sQrTNJZlTV66rqxvHOqaqPVdVnxhi2JzBpRQiwFTChIiQ9/r5JkiRhETLtJdk8yQVJliRZnmRea781yYlJrkiyKMneSS5M8sMkRw+Is2+S65M8aVj7LkmuTrI4ydIkO7X2V/e1fzzJjNZ+T5ITklwF7J/k0iRzRkj/LS3G1Ul2bPOPT3JMe35pkgWt/6Ykz0yyMXACMK+tPS/Jfkkub/lfnuSpbf6RSc5N8rUkP0jyj33H9Rct5qVJPpHk1AHn5De5tNfLk8wGPgA8ua1/UpItklyc5Loky5K8pI2fneS7ST4CXAc8PslH2/fjhiTvHs/3WJIkaX3jR/ROfy8AbquqQwGSbNnX96Oq2j/JB4HTgQOBTYEbgI8NDUpyAPAh4CVV9Z/D4h8N/HNVndkKgBlJng7MAw6sqgfbm+zDgc8AmwPLq+q4Fnu03P+7qvZrl18tBAZ9zuTMNuaFwLuq6pAkxwFzquqNbY1HA8+qqoeSHAKcCLyizd8T2At4APh+kg8BK4F3AnsDdwPfAJaMlugwxwK7VtWebf2ZwMuq6r/bJWZXJjmvjX0q8OdV9Vdt7Nur6hetaLs4ye5VtXT4AkmOAo4CmPHobSaQmiRJ0tRnETL9LQNOTrIAOL+qvt3Xd17fmC2q6m7g7iT3J9mq9T0dOA14flXdNiD+FcDbk/wBcG5V/SDJ84B9gGtakfFI4Gdt/ErgnHHm/vm+xw+OMObc9ngtMHuEMVsCn267NAU8oq/v4qq6CyDJjcAOwNbAN6vqF639LOAp48x5kAAnJnkWsAr4fWDb1reiqq7sG/uqVmDMBGYBOwO/U4RU1Wn0vi9sMmunWoPcJEmSphwvx5rmquomegXBMuD9bZdgyAPtcVXf86HXQwXoT4D76e0WDIr/OeDFwH3AhUmeS+9N96eras/29dSqOr5Nub+qVo43/RGe9xvKeyUjF83vAS6pql2BF9Hb7Rk+vz/GqNszfR7i4b8jm44w7nBgG2Cftjtye9/YXw0NSvJE4BjgeVW1O3DBKDElSZLWWxYh01yS7YF7q+oM4GR6lxhNxJ3AofT+kn/wgPhPAm6uqlPo7azsDlwMHJbk99qYxybZYTXSn9f3eMUE5t0NPKrv9ZbAj9vzI8cx/2rg2Uke0y6lesUI426lnc8kewNPHGX9n7VL055Db7dlkEfTK0ruSrIt8MfjyFWSJGm94+VY099uwElJVgEPAm+YaICquj3Ji4B/S/Laqrqqr3se8OokDwI/BU5o9zS8A7govU98ehD4n8CKCS69SbuBfSPgTycw7xLg2CSLgfcD/0jvcqy/o3d/x6iq6sdJTgSuAm4DbgTuGjD0HOCIts41wE1t/v9L8p0ky4F/AxYAX02yCFgMfG+EdZckuZ7ePTk3A9+ZwDFLkiStN1Ll5eba8CTZoqruaTshXwY+VVVfnuy8Btlk1k41a/7CyU5D0hS1YsGgz/SQpElzbVWN9Mmov+HlWNpQHd92OJYDtwBfmeR8JEmSNhjuhEhTXBJ/SSVJ0nThTogkSZKkqcciRJIkSVKnLEIkSZIkdcoiRJIkSVKnLEIkSZIkdcoiRJIkSVKnLEIkSZIkdcoiRJIkSVKnLEIkSZIkdcoiRJIkSVKnLEIkSZIkdcoiRJIkSVKnZk52ApJGt/F2OzJr/sLJTkPSFLViwdzJTkGSJsydEEmSJEmdsgiRJEmS1CmLEEmSJEmdsghRJ5KckOSQAe2XJpmzFuIfmeTUCYw/OMn5a7quJEmSJs4b07XOJZlRVcdNdh5TVTs/Kyc7D0mSpK64E6I1kmTzJBckWZJkeZJ5rf3WJMcluQx4ZZLTkxw2SpyNknw6yXvb63v6+g5Lcnp7/qIkVyW5Psm/J9l2QKwxxwwbf3ySY/peL08yu319L8knW9uZSQ5J8p0kP0iyX9/8zyb5Rmt/fWt/2G5LklOTHDno/Ix9piVJktYf7oRoTb0AuK2qDgVIsmVf3/1VdVBrf8EoMWYCZwLLq+p9Y6x3GfCMqqokrwPeCvz9aowZrx3pFQlHAdcAfwYcBLwYeBvw0jZud+AZwObA9UkuGEfs35yf4ZIc1dZkxqO3Wc3UJUmSpiaLEK2pZcDJSRYA51fVt/v6vjjOGB8HvjSOAgTgD4AvJpkFbAzcsppjxuuWqloGkOQG4OJW3CwDZveN+79VdR9wX5JLgP2AO8eIPeL5qarTgNMANpm1U61B/pIkSVOOl2NpjVTVTcA+9IqR9yfpv/fjV+MMcznwnCSb9ofue97f/iHg1KraDfjLYX0TGdPvIR7+u9A//oG+56v6Xq/i4UX88EKhxogL4z8/kiRJ6xWLEK2RJNsD91bVGcDJwN6rEeb/AP8KnJVk6I397UmenmQj4GV9Y7cEftyezx8h3njG9LuVlneSvYEnjj/133hJkk2TPA44mN6lWyuAnZNs0i5Te95qxJUkSVrveDmW1tRuwElJVgEPAm9YnSBV9U/tjfpnkxwOHAucD/wIWA5s0YYeT69Y+TFwJYMLhvGM6XcOcESSxfSKh5tW4xCuBi4AngC8p6puA0jyJWAp8APg+tWIK0mStN5JlZebS2siyfHAPVV18rqIv8msnWrW/IXrIrSk9cCKBXMnOwVJ6ndtVY35f8B5OZYkSZKkTrkTIk1xSfwllSRJ04U7IZIkSZKmHosQSZIkSZ2yCJEkSZLUKYsQSZIkSZ2yCJEkSZLUKYsQSZIkSZ2yCJEkSZLUKYsQSZIkSZ2yCJEkSZLUKYsQSZIkSZ2yCJEkSZLUKYsQSZIkSZ2yCJEkSZLUqZmTnYCk0W283Y7Mmr9wstOQNEWtWDB3slOQpAlzJ0SSJElSpyxCJEmSJHXKIkTjkuTgJAf0vT49yWETmD87yfL2fE6SUyYw99YkWw9of3GSY8eR9/njXWuc+dyzluIcn+SYtRFLkiRpOvGeEI3XwcA9wOVrGqiqFgGL1kKc84Dz1jSOJEmSuuVOyAYsyRFJliZZkuSzre1FSa5Kcn2Sf0+ybZLZwNHA3yZZnOSZLcSzklye5OahXZH0nJRkeZJlSeYNWPc3uxNJtkjyL23s0iSvGCHdNyW5ro17Wpt7ZJJT2/MnJ7kyyTVJThi2W7FFkrOTfC/JmUkyIKc3J7mx5fCFsXJL8r523q5Msm1r2yHJxW3sxUmeMFq7JEnShsoiZAOVZBfg7cBzq2oP4K9b12XAM6pqL+ALwFur6lbgY8AHq2rPqvp2GzsLOAiYC3ygtb0c2BPYAzgEOCnJrFFSeSdwV1XtVlW7A98YYdwdVbU38FFg0CVM/wz8c1XtC9w2rG8v4G+AnYEnAQcOmH8ssFfL4egxctscuLKdt28Br2/tpwKfaWPPBE4Zo31ESY5KsijJopX33jXWcEmSpGnFImTD9Vzg7Kq6A6CqftHa/wC4MMky4C3ALqPE+EpVraqqG4FtW9tBwOeramVV3Q58E9h3lBiHAB8eelFVvxxh3Lnt8Vpg9oD+/YGz2vPPDeu7uqr+q6pWAYtHmL8UODPJq4GHxsjt18DQfSb9+ezft/Zn6Z2L0dpHVFWnVdWcqpozY7MtxxouSZI0rViEbLgC1ID2DwGnVtVuwF8Cm44S44Fh8fof1zSPkdZaycTvZerPc6T5h9IrOPYBrk0yc5TcHqyqofbR8hnpuMZzvJIkSesti5AN18XAq5I8DiDJY1v7lsCP2/P5fePvBh41jrjfAuYlmZFkG+BZwNWjjL8IeOPQiySPGV/6v+NKYOiejT+ZyMQkGwGPr6pLgLcCWwFbrEZul/etfTi9S9tGa5ckSdogWYRsoKrqBuB9wDeTLAH+qXUdD5yV5NvAHX1Tvgq8bNiN6YN8md6lTUvo3UPx1qr66Sjj3ws8pt3IvgR4zmodUO+ej79LcjW9e1UmciPFDOCMdgna9fTufblzNXJ7M/DnSZYCr+G399mM1C5JkrRBym+vKpGmrySbAfdVVSX5E+BPq+olk53X2rDJrJ1q1vyFk52GpClqxYK5k52CJPW7tqrmjDXI/ydE64t9gFPbx+/eCbx2kvORJEnSCNwJkaa4JP6SSpKk6WJcOyHeEyJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUzMnOwFJo9t4ux2ZNX/hZKchaYpasWDuZKcgSRPmTogkSZKkTlmESJIkSeqURYgkSZKkTlmEaJ1IcnCS89dB3E8m2XmMMbcm2XqMMUcmOXXtZidJkqTx8MZ0TStV9brJzgEgyYyqWjnZeUiSJE1H7oRoREk2T3JBkiVJlieZ19pvTXJikiuSLEqyd5ILk/wwydED4uyb5PokTxrWvlmSLyVZmuSLSa5KMqf1fbTFviHJu/vmXDrWmOZNSa5LsizJ00Y4xMcn+VqS7yd5V98ar05ydZLFST6eZEZrvyfJCUmuAvZPclySa9q5OS1J+nJc0GLclOSZYx2vJEnShsQiRKN5AXBbVe1RVbsCX+vr+1FV7Q98GzgdOAx4BnBCf4AkBwAfA15SVTcPi/9XwC+ranfgPcA+fX1vr6o5wO7As5PsPiC/0cbcUVV7Ax8Fjhnh+PYDDgf2BF6ZZE6SpwPzgAOrak9gZRsDsDmwvKr+sKouA06tqn3buXkk0P85mTOraj/gb4ChAme0432YJEe1AmvRynvvGmmYJEnStGQRotEsAw5pf9V/ZlX1vxs+r2/MVVV1d1X9HLg/yVat7+nAacCLquo/B8Q/CPgCQFUtB5b29b0qyXXA9cAuwKD7QEYbc257vBaYPcLxfb2q/l9V3dfGHwQ8j15xcE2Sxe310A7OSuCcvvnPabsZy4DnthxGW3+0432YqjqtquZU1ZwZm2050jBJkqRpyXtCNKKquinJPsALgfcnuaiqhnY6HmiPq/qeD70e+rn6CbApsBdw24AlMmjdJE+kt3uxb1X9MsnpLc5ExgzltJKRf85rwOsAn66qfxgw/v6h+0CSbAp8BJhTVT9Kcvw41h94vJIkSRsad0I0oiTbA/dW1RnAycDeEwxxJ3AocGKSgwf0Xwa8qq21M7Bba3808CvgriTbAn88YO54xozlj5I8NskjgZcC3wEuBg5L8nstr8cm2WHA3KGC444kW9C7HG0sIx2vJEnSBsWdEI1mN+CkJKuAB4E3TDRAVd2e5EXAvyV5bVVd1df9EeDTSZbSu6RqKXBXVf0gyfXADcDN9IqDYWFryRhjxuMy4LPAjsDnqmoRQJJ3ABcl2Yjecf9PYMWwBO5M8gl6l6PdClwzjvUGHu9q5C1JkjStpWr4FSlSN9qnTj2iqu5P8mR6uxBPqapfjzJnGfDiqrqlqzzXltU5XoBNZu1Us+Yv7CRHSdPPigVzxx4kSd25tn1w0KjcCdFk2gy4JMkj6N0v8YYxCpCvA8umYwHSTOh4h/z6p//hmwxJkrRecSdEmuKS+EsqSZKmi3HthHhjuiRJkqROWYRIkiRJ6pRFiCRJkqROWYRIkiRJ6pRFiCRJkqROWYRIkiRJ6pRFiCRJkqROWYRIkiRJ6pRFiCRJkqROWYRIkiRJ6pRFiCRJkqROWYRIkiRJ6tTMyU5A0ug23m5HZs1fONlpSFPOigVzJzsFSdJqcidEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQqRRJLk8yfZJzm6v90zywonGWDfZSZIkTU8WIdIoquqAqrqtqg5rTXsCEypCquqAtZ+ZJEnS9GURIo0iyT1JZidZnmRj4ARgXpLFSeYNG7tLkqtb39IkOw3FaI8HJ7k0ydlJvpfkzCTp/qgkSZImlx/RK41TVf06yXHAnKp644AhRwP/XFVntoJlxoAxewG7ALcB3wEOBC4bPijJUcBRADMevc1aOgJJkqSpwZ0Qae25Anhbkv8F7FBV9w0Yc3VV/VdVrQIWA7MHBaqq06pqTlXNmbHZlusuY0mSpElgESKtJVX1OeDFwH3AhUmeO2DYA33PV+JupCRJ2gBZhEgTczfwqEEdSZ4E3FxVpwDnAbt3mZgkSdJ0YREija6Gvb4E2HnQjenAPGB5ksXA04DPdJGgJEnSdOOlINIIkjwO+EVV3QrsClBVvwD2HTS+qt4PvH9A+xbt8VLg0r72QTe3S5IkrfdSNfwPvZKSbE+vYPhQVX1oknPxl1SSJE0X11bVnLEGuRMiDVBVtwFPmew8JEmS1kfeEyJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjplESJJkiSpUxYhkiRJkjo1c7ITkDS6jbfbkVnzF052GtKUs2LB3MlOQZK0mtwJkSRJktQpixBJkiRJnbIIkSRJktQpi5BpJMnsJMunQB63Jtm6w/WOT3LMCH2Xj2P+Ws03yelJDlsLcabE91OSJKlrFiGa1qrqgMnOQZIkSRNjETJNJXlSkuuT7JvkyCSn9vWdn+Tg9vz5Sa5Icl2Ss5Js0dpvTfLu1r4sydMGrDEjycmtf2mSN/V1v2n43CSbJ/lUkmtabi9p7UcmOTfJ15L8IMk/9q3xF0luSnJpkk/0H8cwO7cxNyd5c9/8e9rjRkk+kuSGdvz/Omy34nfyHXasuyS5Osnidqw7tfYj2uslST7bN+VZSS5v+RzWxibJSUmWt3XmjdYuSZK0ofIjeqehJE8FvgD8eVUtTrLLCOO2Bt4BHFJVv0ryv4C/A05oQ+6oqr2T/BVwDPC6YSGOAp4I7FVVDyV5bF/foLlvB75RVa9NshVwdZJ/b+P3BPYCHgC+n+RDwErgncDewN3AN4AlIxz204DnAI9q8z9aVQ/29b8cmA3sBvwe8F3gU2Pk2+9o4J+r6swkGwMz2nl9O3BgVd0x7PhnAQe1vM4Dzm457AnsAWwNXJPkW8ABI7SPKMlR9M4/Mx69zWhDJUmSph13QqafbYD/C7y6qhaPMfYZwM7Ad5IsBuYDO/T1n9ser6X3Bn64Q4CPVdVDAFX1izHmPh84tq11KbAp8ITWd3FV3VVV9wM3tjz2A75ZVb9oBcVZoxzLBVX1QFXdAfwM2HZY/0HAWVW1qqp+ClwyrH+sY70CeFsr1HaoqvuA5wJntzWHH/9X2lo39uVyEPD5qlpZVbcD3wT2HaV9RFV1WlXNqao5MzbbcrShkiRJ0447IdPPXcCPgAOBG1rbQzy8oNy0PQb4elX96QixHmiPKxn8sxCgJjA3wCuq6vsPC5L8Yd/4/jkZIfZo642U71ixRj3WqvpckquAQ4ELk7yO8R1//9oj5TCR45QkSVrvuRMy/fwaeClwRJI/a223Anu2+yIeT2+HAeBK4MAkOwIk2SzJUyaw1kXA0UlmtvmPHWP8hfTuvUgbv9cY468Gnp3kMW2NV0wgt+EuA17RzsG2wMETmZzkScDNVXUKvcurdgcuBl6V5HFtzFjH/y1gXruXZhvgWfSOcaR2SZKkDZI7IdNQu79jLvD1JL+i96b5FmAZsBy4ro37eZIjgc8n2aRNfwdw0ziX+iTwFGBpkgeBTwAj3TgO8B5gYRsfesXR3FGO48dJTgSuAm6jd5nWXePMbbhzgOfRO/6bWsyJxJoHvLod50+BE6rqF0neB3wzyUrgeuDIUWJ8Gdif3n0tBby1qn6aZKT22RPIT5Ikab2RqpGuNpHWvSRbVNU9bSfky8CnqurLaxjrcfR2Gg5s94dMa5vM2qlmzV842WlIU86KBSP+jUOSNHmurao5Yw1yJ0ST7fgkh9ApXH2JAAAgAElEQVS7j+Ui4CtrEOv89qlcGwPvWR8KEEmSpPWROyHSFJfEX1JJkjRdjGsnxBvTJUmSJHXKIkSSJElSpyxCJEmSJHXKIkSSJElSpyxCJEmSJHXKIkSSJElSpyxCJEmSJHXKIkSSJElSpyxCJEmSJHXKIkSSJElSpyxCJEmSJHXKIkSSJElSp2ZOdgKSRrfxdjsya/7CyU5DmnJWLJg72SlIklaTOyGSJEmSOmURIkmSJKlTFiGSJEmSOmURomkhyQlJDpnsPCRJkrTmvDFdU16SGVV13GTnMR5JAqSqVk12LpIkSVOVOyGaVEk2T3JBkiVJlieZ19pvTXJcksuAVyY5PclhA+a/OcmNSZYm+UJfzE8luSbJ9Ule0tqPTPKVJF9NckuSNyb5uzbmyiSPbeMuTbIwyeUtp/1a+/FJjulbe3mS2e3ru0k+AlwHPD7JPUne147ryiTbtjkvSnJVW/Pfh9olSZI2JBYhmmwvAG6rqj2qalfga31991fVQVX1hVHmHwvsVVW7A0e3trcD36iqfYHnACcl2bz17Qr8GbAf8D7g3qraC7gCOKIv7uZVdQDwV8CnxnEcTwU+U1V7VdUKYHPgyqraA/gW8Po27jLgGW3NLwBvHRQsyVFJFiVZtPLeu8axvCRJ0vRhEaLJtgw4JMmCJM+sqv533F8cx/ylwJlJXg081NqeDxybZDFwKbAp8ITWd0lV3V1VPwfuAr7al8fsvrifB6iqbwGPTrLVGHmsqKor+17/Gji/Pb+2L/YfABcmWQa8BdhlULCqOq2q5lTVnBmbbTnG0pIkSdOLRYgmVVXdBOxDrwh4f5L+ez9+NY4QhwIfbjGuTTITCPCKqtqzfT2hqr7bxj/QN3dV3+tVPPweqRqeKr0ip/93ZtNRcn2wqoZirOyL/SHg1KraDfjLYTEkSZI2CBYhmlRJtqd3SdQZwMnA3hOYuxHw+Kq6hN5lTVsBWwAXAm9qN4mTZK/VSG3o3pSDgLvaDs2tQ/kl2Rt44mrE3RL4cXs+fzXmS5IkTXt+OpYm22707tlYBTwIvGECc2cAZyTZkt7uxwer6s4k7wEWAktbIXIrMHeCef0yyeXAo4HXtrZzgCPaZV7XADdNMCbA8cBZSX4MXMnqFTKSJEnTWn57xYgk6H06FnBMVS2a7FwANpm1U82av3Cy05CmnBULJvq3BUlSB66tqjljDfJyLEmSJEmdcidEmuKS+EsqSZKmC3dCJEmSJE09FiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSJKlTFiGSJEmSOmURIkmSJKlTFiGSJEmSOjVzshOQNLqNt9uRWfMXTnYa0pSzYsHcyU5BkrSa3AmRJEmS1CmLEEmSJEmdsgiRJEmS1CmLEK2RJKcnOWxA+/ZJzl5La1yaZM6A9jlJThlj7uwky9dGHn0xb02y9VqIc2SSU9dGTpIkSdOJN6Zrnaiq24DfKU7W8hqLgEXrcg1JkiStfe6ErGfaX/6/m+QTSW5IclGSR7a+1ye5JsmSJOck2ay1n57klCSXJ7l50M5GG3dEkqVt/mf7up41fG7/DkR7/u0k17WvA1r7wUnO74t/apIjRzi0Vya5OslNSZ45fH6SbZJ8vcX/eJIVfbsVMwadj2HH9soky9uxfau1zUhycpJl7bjf1DflTW2tZUme1sY/NslX2tgrk+w+WrskSdKGyiJk/bQT8OGq2gW4E3hFaz+3qvatqj2A7wJ/0TdnFnAQMBf4wPCASXYB3g48t83/6/HOBX4G/FFV7Q3MA0a9hGoEM6tqP+BvgHcN6H8X8I22xpeBJ/T1jXQ++h0H/I92bC9ubUcBTwT2qqrdgTP7xt/R1voocExrezdwfRv7NuAzY7SPKMlRSRYlWbTy3rvGGi5JkjStWISsn26pqsXt+bXA7PZ817YjsQw4HNilb85XqmpVVd0IbDsg5nOBs6vqDoCq+sUE5j4C+ERb9yxg59U4pnMHHE+/g4AvtNy+Bvyyr2+k89HvO8DpSV4PzGhthwAfq6qHWtz+Yx6Uz0HAZ9vYbwCPS7LlKO0jqqrTqmpOVc2ZsdmoQyVJkqYd7wlZPz3Q93wlMHT50enAS6tqSbvs6eAR5mRAzAA1jvUGzf1b4HZgD3qF7/2t/SEeXghvOkL8/jVWMvjndtC6g/LrPx+/UVVHJ/lD4FBgcZI9Gd8x9+czKIcapV2SJGmD5E7IhuVRwE+SPILeTshEXAy8KsnjoHefwwTmbgn8pKpWAa/htzsNK4Cdk2zSdgaeN8Gc+l0GvKrl9nzgMROZnOTJVXVVVR0H3AE8HrgIODrJzDZmrGP+Fu28JjmY3iVb/z1KuyRJ0gbJnZANyzuBq+i9+V9GrygZl6q6Icn7gG8mWQlcDxw5zukfAc5J8krgEuBXLeaPknwJWAr8oMVcXe8GPp9kHvBN4CfA3cAW45x/UpKd6O1aXAwsAZYDTwGWJnkQ+AQw2kfqHg/8S5KlwL3A/DHaJUmSNkip8qoQTX9JNgFWVtVDSfYHPlpVe052XmvDJrN2qlnzF052GtKUs2LB3MlOQZL0u66tqt/5/92GcydE64snAF9KshHwa+D1k5yPJEmSRuBOiDTFJfGXVJIkTRfj2gnxxnRJkiRJnbIIkSRJktQpixBJkiRJnbIIkSRJktQpixBJkiRJnbIIkSRJktQpixBJkiRJnbIIkSRJktQpixBJkiRJnbIIkSRJktQpixBJkiRJnbIIkSRJktQpixBJkiRJnZo52QlIGt3G2+3IrPkLJzsNacpZsWDuZKcgSVpN7oRIkiRJ6pRFiCRJkqROrdMiJMnRSY5Yl2tMhiSfTLLzas5929rOZxxrXppkzoD2MY8jyUtX91gnavi5SXL5Oljj9CSHrWGMe9ZSLscnOWZtxJIkSZpO1lkRkmRmVX2sqj6zrtaYDElmVNXrqurG1QzReREyknEex0uBCRUhSVb3XqOHnZuqOmA140iSJGkKG7UISbJ5kguSLEmyPMm81r5Pkm8muTbJhUlmtfZLk5yY5JvAX4/0l94kOyS5OMnS9viE1n56ko8muSTJzUmeneRTSb6b5PS++c9PckWS65KclWSL1v7CJN9LclmSU5Kc39r3S3J5kuvb41MH5HRwkm8l+XKSG5N8LMlGre+eJCckuQrYf2hnIckbkvxjX4wjk3yoPf9KOz83JDmqtX0AeGSSxUnObG2vTnJ1a/t4khkDcjsuyTXte3BakvSd7wVt/k1JntnaH5nkC+38fhF45Ajf39/skLRjfF/7Xl+ZZNskBwAvBk5q+T25fX2tHdu3kzyt73v3T0kuARa07/2n2ho3J3lz37rjPTf3tMckOakd/7K+n8ODW/yz2/f9zL5zM/CcjSTJm9v3fWmSL7S2LZL8S1tzaZJX9I1/2LlqbSP9XA9slyRJ2lCNtRPyAuC2qtqjqnYFvpbkEcCHgMOqah/gU8D7+uZsVVXPrqr/PUrcU4HPVNXuwJnAKX19jwGeC/wt8FXgg8AuwG5J9kyyNfAO4JCq2htYBPxdkk2BjwN/XFUHAdv0xfwe8Kyq2gs4DjhxhLz2A/4e2A14MvDy1r45sLyq/rCqLusbf3bfGIB5wBfb89e28zMHeHOSx1XVscB9VbVnVR2e5OltzoFVtSewEjh80Pmqqn3b9+CRQP9Hwsysqv2AvwHe1dreANzbzu/7gH1GON5+mwNXVtUewLeA11fV5cB5wFtazj8ETgPe1I7tGOAjfTGeQu/78vft9dOA/0HvvL6r/eyM69wMy+3lwJ7AHsAh9IqiWa1vr3bsOwNPAg4cxzkb5Fhgr3bOjm5t7wTuqqrdWvs3RjpXQ2sy+Od6tJ/3gZIclWRRkkUr771rrOGSJEnTylhFyDLgkPbX9mdW1V3AU4Fdga8nWUyvIPiDvjlfHBBnuP2Bz7XnnwUO6uv7alVVW/v2qlpWVauAG4DZwDPoveH8Tlt/PrADvTe8N1fVLS3O5/tibgmclWQ5vy1qBrm6qm6uqpVt/lBeK4Fzhg+uqp8DNyd5RpLH0Ts332ndb06yBLgSeDyw04D1nkevQLimHcvz6L2RHu45Sa5Ksoxegdaf/7nt8Vp65wfgWcAZLcelwNIRjrffr4HzB8T6jfR2nA6gdy4X0yv6ZvUNOauduyEXVNUDVXUH8DNg29Y+nnPT7yDg81W1sqpuB74J7Nv6rq6q/2o/I4v78h7tnA2yFDgzyauBh1rbIcCHhwZU1S/b05HO1Ug/16P9vA9UVadV1ZyqmjNjsy3HGi5JkjStjHrtflXdlGQf4IXA+5NcBHwZuKGq9h9h2q9WI4/qe/5Ae1zV93zo9Ux6BcHXq+pP+wMk2WuU+O8BLqmqlyWZDVw6jjz6X98/7M11vy8Cr6K32/LlqqokB9N7A7t/Vd2b5FJg0wFzA3y6qv5hpMTbDs9HgDlV9aMkxw+LNXSOVvLw7+fwYxnLg634GxRryEbAnW3XZpDh3/v+799KYOYEzk2/0S6lGrTGWOdskEPpFW8vBt6ZZJe27qDzOJ5zxQhzR2uXJEnaIIx1T8j29C7rOQM4Gdgb+D6wTZL925hHtDdsE3E58Cft+eHAZaOMHe5K4MAkO7b1N0vyFHpFwJNakQG9y5yGbAn8uD0/cpTY+yV5Ynr3gswbZ17n0rt5+0/57S7QlsAv25vsp9HbvRnyYN9lSRcDhyX5vXYsj02yw7D4Q2+e72g7EeP5ZKdv0S7rSrIrsPs45ozkbuBRAFX138AtSV7ZYifJHhOMN95z0+9bwLwkM5JsQ69YuHqUNSZ0ztr3+/FVdQnwVmArYAvgIuCNfeMeM/qhjfhzvSY/75IkSeudsS7H2g24ul1683bgvVX1a3pv6ha0S2oW07tEZyLeDPx5kqXAa4C/Hu/EdgnUkcDn2/wrgadV1X3AX9G7b+Uy4HZg6GL6f6S3k/Md4Hdu/O5zBfABYDlwC71dn7Hy+SVwI7BDVQ29Mf4avb/IL6W3C3Nl35TTgKVJzmyfTPUO4KI29us8/PImqupO4BP0Lk/7CnDNWDkBHwW2aDHfyuhv2MfyBeAt6d3U/2R6b6L/on3vbwBeMsF44zo3w+Z8md7lUkvo3Zfx1qr66UgLrMY5mwGc0S7duh74YIvxXuAx7eb2JcBzxogz0s/1av+8S5IkrY/y26tKpr8kW1TVPe2TkD4M/KCqPjjOuQcDx1TVWDcwS53aZNZONWv+wslOQ5pyVizwn2tJmoKurarf+f/phlvf/sf017ddmxvoXfbz8UnOR5IkSdIw69VOiLQ+SuIvqSRJmi42yJ0QSZIkSVOcRYgkSZKkTlmESJIkSeqURYgkSZKkTlmESJIkSeqURYgkSZKkTlmESJIkSeqURYgkSZKkTlmESJIkSeqURYgkSZKkTlmESJIkSeqURYgkSZKkTlmESJIkSerUzMlOQNLoNt5uR2bNXzjZaUhTzooFcyc7BUnSanInRJIkSVKnLEIkSZIkdcoiRJIkSVKnLEI0LSXZPsnZ7fmRSU4dY/zxSY5ZB3n8Zu0kRyc5Ym2vIUmStL7xxnRNS1V1G3DYZOfRr6o+Ntk5SJIkTQfuhGjSJNk8yQVJliRZnmRea781yYlJrkiyKMneSS5M8sMkR7cxs5MsHxDz0DZv6wFL7pzk0iQ3J3nzoDhJjklyfHv++iTXtPzOSbLZGMfzm92WkeYmOT3JKUkub3lMqUJKkiSpCxYhmkwvAG6rqj2qalfga319P6qq/YFvA6fT2/V4BnDCSMGSvAw4FnhhVd0xYMjTgP8B7Ae8K8kjxsjv3Krat6r2AL4L/MX4DmvMubOAg4C5wAdGOJajWgG2aOW9d01gWUmSpKnPy7E0mZYBJydZAJxfVd/u6zuvb8wWVXU3cHeS+5NsNSDWc4A5wPOr6r9HWO+CqnoAeCDJz4Btx8hv1yTvBbYCtgAuHN9hjTn3K1W1CrgxycAcquo04DSATWbtVBNYV5IkacpzJ0STpqpuAvahV2i8P8lxfd0PtMdVfc+HXg8qnm8GHgU8ZZQl++OsbHEe4uG/B5v2PT8deGNV7Qa8e1jfWEab259HJhBTkiRpvWARokmTZHvg3qo6AzgZ2HsNwq0AXg58JskuE5h3O/B7SR6XZBN6l0gNeRTwk3bZ1uETzGdN5kqSJK3XvBxLk2k34KQkq4AHgTdMcP7DLlOqqu8nORw4K8mLquqHYwaoejDJCcBVwC3A9/q639naV9DbrXnUBHJbk7mSJEnrtVR5ubmmnyT7AP9UVc+e7FzWtU1m7VSz5i+c7DSkKWfFgrljD5Ikde3aqpoz1iCLEE07SeYAnwOOrapzJzufdS2Jv6SSJGm6GFcR4uVYmnaqahGj34AuSZKkKcwb0yVJkiR1yiJEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQiRJkiR1yiJEkiRJUqcsQiRJkiR1yiJEkiRJUqdmTnYCkka38XY7Mmv+wslOQ5pyViyYO9kpSJJWkzshkiRJkjplESJJkiSpUxYhkiRJkjplEaINQpKDk5y/lmOekOSQAe2vTHJDklVJ5vS175dkcftakuRlazMfSZKk6cIb06XVVFXHjdC1HHg58PEB7XOq6qEks4AlSb5aVQ+tyzwlSZKmGndCNG0k2TzJBW0XYXmSea391iQnJrkiyaIkeye5MMkPkxw9IM6+Sa5P8qRh7UcmObXv9fltB2VGktPbmsuS/G3rPz3JYcPjV9V3q+r7A9rv7Ss4NgVqzc6IJEnS9OROiKaTFwC3VdWhAEm27Ov7UVXtn+SDwOnAgfTe6N8AfGxoUJIDgA8BL6mq/xznunsCv19Vu7YYW63uAST5Q+BTwA7Aa0baBUlyFHAUwIxHb7O6y0mSJE1J7oRoOlkGHJJkQZJnVtVdfX3n9Y25qv5/e3ceZllV3n3/+6ObeVThhcKBRm0HFGigQVHBRtGotCJRQ4xR0LwQYpSgF/ExMVFMnDryKApGBURUEKKAiKCAQUaRoSe6AUEj0IogSARkHu/nj7NKDkUNXXRzThX1/VxXXbVr7TXce50q2PdZa5+uuqOqfg/c25U0vBA4HHjjOBIQgGuAZyc5NMnrgD8+3guoqour6kXA9sA/JVljhHqHV9Xsqpo9ba31h6siSZI0aZmEaNKoql8A29FJND6dpPuZjPva94e7jgd/HlzxuxG4F9hmhCEe5NF/E2u0cW8FtgbOAf4eOPJxX0RTVT8H7gJevKJ9SZIkTTZux9KkkWRT4A9VdUySO4G9x9nFbcDfAGcmuauqzhly/jrgvUlWAZ4O7NDG3RC4v6pOTPIrOtu9Hk/8m9PZNvZgks2A57cxJUmSphSTEE0mWwKfTfIw8ADwd+PtoKpuSvJG4EdJ3lNVF3ed/ilwLZ2VlsuBha386cDXW3IC8E/dXQ4do3307qHARsBpSRZX1Z8BrwA+nOQBOis0762qW8Z7DZIkSZNdqvyAHunxSPID4HNVdfYTOc7qAzNrYK9DnsghpElp2by5/Q5BkvRYC6pq9liVfCZEehySHAWsBVzQ71gkSZImG1dCpAkuiX+kkiRpsnAlRJIkSdLEYxIiSZIkqadMQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPWUSIkmSJKmnTEIkSZIk9ZRJiCRJkqSemt7vACSNbrVNnsvAXof0OwxNEMvmze13CJIkrTBXQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhE1SSGUkuH2ebg5IcOI76s5N8cYw6GyR573jieLySzEnysq6f90vyrpU8xrjndZg+xjXPY/R158roR5IkaTLxwfQprKrmA/PHqLYB8F7gP5e33yQBUlUPjzOkOcCdwIUtvq+Ms70kSZImAVdCJrZpSY5IckWSM5OsCZDkOUlOT7IgyflJXjC0YZJzkhyS5MIklyfZYZg6c5Kc2o4PSnJUa3dNkv1btc8Az0myOMlnW91/THJpkiVJPt7KZiT5eZL/BBYCz0xyZ5JPJrksyUVJNm5135jk4iSLkvx3ko2TzAD2Az7Qxtqpe8UhyazWx5Ik30vylK7rnJfkkiS/SLJTVzznJ1nYvv60wjKcJANJzmtjX97Vz+ta+8uSnNXVZIth5ookH2ztL09ywFjlkiRJU5FJyMQ2E/hSVb0IuA14Sys/HHh/VW0HHMjIqxRrV9XL6KxkHLUc470A+DNgB+BjSVYFPgz8qqpmVdU/Jnlti2sHYBawXZKdW/vnA9+sqm2qahmwNnBRVW0NnAfs0+pdALy0qrYBjgc+VFXXAV8BPt/GOn9IbN8E/k9VbQUsBT7WdW56Ve0AHNBVfjPwmqraFtgTGHXbGfBXwBlVNQvYGlicZCPgCOAt7RreNtpcJdkOeDfwEuClwD5JthmpfLRgkuybZH6S+Q/dffsYoUuSJE0ubsea2K6tqsXteAEwI8k6wMuA73Z2PQGw+gjtjwOoqvOSrJdkg6q6bZTxTquq+4D7ktwMbDxMnde2r0Xt53XoJCW/BpZV1UVdde8HTu2K/zXt+BnAfyUZAFYDrh0lJpKsD2xQVee2om8A3+2qclLXGDPa8arAYUlmAQ8BzxttDOBS4KiWeJ1cVYuTzAHOq6prAarqD131h5urVwDfq6q7WtwnATsBGaF8ESOoqsPpJJusPjCzxohdkiRpUjEJmdju6zp+CFiTzurVbe0d+7EMvXkd62Z26HjD/X4E+HRVffVRhZ3tVHcNqftAVQ2O2d3focDnquqUdqN/0BhxjWUw7u4xPgDcRGdVYxXg3tE6aInazsBuwLfa1rPbGHnOhpurjFB3pHJJkqQpye1Yk0xV/RG4NsnboPMQeJKtR6i+Z6vzCuD2qno8+3ruANbt+vkM4D1tRYYkT0/y/42zz/WB37bjvUYZC4AW962Dz2kA7wTOHVpvmDFubA/HvxOYNlrlJJsBN1fVEcDXgG2BnwGvTLJ5q/PUMcY8D3hzkrWSrA3sAZw/SrkkSdKU5ErI5PQO4MtJ/oXOtqPjgcuGqXdrkguB9YD3PJ6Bqup/k/w0nY+1/VF7LuSFwM/adrA7gb+msxqwvA6is53st8BFwOat/AfACUl2B94/pM1ewFeSrAVcQ+cZi9H8J3BiS9bO5rGrNEPNAf4xyQN0ruldVfX7JPsCJyVZhfacyUgdVNXCJEcDl7SiI6tqEcBI5ZIkSVNRHtktoyeTJOcAB7aP4dUktvrAzBrY65B+h6EJYtm8uf0OQZKk0SyoqtljVXI7liRJkqSeciVEmuCS+EcqSZImC1dCJEmSJE08JiGSJEmSesokRJIkSVJPmYRIkiRJ6imTEEmSJEk9ZRIiSZIkqadMQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeqp6f0OQNLoVtvkuQzsdUi/w9AEsWze3H6HIEnSCnMlRJIkSVJPmYRIkiRJ6imTEEmSJEk9ZRKiUSWZk+TUcbY5J8nsJyqmfkvywyQb9DsOSZKkycoH0yUgyfSqenB56lbVG57oeCRJkp7MXAmZApKsneS0JJcluTzJnq38uiSfSvKzJPOTbJvkjCS/SrLfMP1sn2RRkmcPKV8zyfFJliT5L2DNrnNfbn1fkeTjXeVjjp1knSRnJVmYZGmS3Vv5jCQ/T3JE6/fMJGu2c/snubLFcnwr2yHJhS32C5M8v5XvneS7SX4AnDnSeMPMw3VJNmzHJydZ0OLYt5X9TZLPd9XfJ8nnRqovSZI01bgSMjW8DrihqnYDSLJ+17nfVNWO7ab5aODlwBrAFcBXBisleRlwKLB7Vf16SP9/B9xdVVsl2QpY2HXuI1X1hyTTgLOSbFVVS5Zz7HuBParqj+2m/6Ikp7S2M4G3V9U+Sb4DvAU4BvgwsHlV3de1ZeoqYOeqejDJrsCnWn2AHYGtWozThxuvqmqUuX1Pa7smcGmSE4HjgSVJPlRVDwDvBv52pPpV9b9DO20Jyr4A09bbaJThJUmSJh+TkKlhKXBwknnAqVV1fte5U7rqrFNVdwB3JLm36yb+hcDhwGur6oZh+t8Z+CJAVS1JsqTr3F+0G+rpwACwBTB4fqyx7wI+lWRn4GHg6cDGrc21VbW4HS8AZrTjJcCxSU4GTm5l6wPfSDITKGDVrvh+XFV/aMcZYbzfDXPNg/ZPskc7fiYws6ouSvITYG6SnwOrVtXSkeoDj0lCqupwOnPO6gMzR0uCJEmSJh23Y00BVfULYDs6N/ufTvLRrtP3te8Pdx0P/jyYpN5IZ1Vim9GGGVqQZHPgQODVVbUVcBqdlY7lHfsdwEbAdlU1C7ipq313/Ye6Yt0N+BKd613QVjf+HTi7ql4MvHFIDHd1HY823mMkmQPsCuxYVVsDi7rqHwnsTWcV5OvLUV+SJGnKMAmZApJsSme71DHAwcC24+ziNjo3959qN9JDnUfnBp4kLwa2auXr0bnJvz3JxsDrxznu+sDNVfVAkl2AzUarnGQV4JlVdTbwIWADYJ3Wz29btb1X1nit/q1VdXeSFwAvHTxRVRfTWen4K+C4sepLkiRNJW7Hmhq2BD6b5GHgATrPcIxLVd2U5I3Aj5K8p91kD/oy8PW2DWsxcElrc1mSRXSe8bgG+Ok4hz0W+EGS+a3fq8aoPw04pj3zEuDzVXVbkv+gsx3rg8BPVtJ4BZwO7Neu+2rgoiF1vgPMqqpb289j1ZckSZoSMvozt5K6tQfsbwY2aQ+dj1b3VDqJ0FkrMubqAzNrYK9DVqQLPYksmze33yFIkjSaBVU15r8X53YsaXyuAI4cLQFJskGSXwD3rGgCIkmS9GTkSog0wSXxj1SSJE0WroRIkiRJmnhMQkT7UuYAACAASURBVCRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPWUSIkmSJKmnTEIkSZIk9ZRJiCRJkqSeMgmRJEmS1FMmIZIkSZJ6anq/A5A0utU2eS4Dex3S7zA0QSybN7ffIUiStMJcCZEkSZLUUyYhkiRJknrKJGSCSnJAkrWWo95+Sd7Vg3hmJLl8hHNHJ5nzBIx54eNsNyPJX63seLr63zvJYSupr+uSbLgy+pIkSZosTEImrgOAMZOQqvpKVX2zB/H0TJJpAFX1ssfZxQzgCUtCJEmStGJMQvosydpJTktyWZLLk+yZZH9gU+DsJGeP0f6gJAe241lJLkqyJMn3kjylle+f5MpWfnzXuEcluTTJoiS7t/K9k3w/yelJrk7ysa7hpiU5IskVSc5MsmYrvx24v7X/TNdYB48Q77eS/CTJL5Ps08rnJDk7ybeBpa3szvb9v5K8oauPo5O8pa14nJ9kYfsaTFo+A+yUZHGSDySZluSz7VqXJPnb5XkdWvn2SS5s5ZckWbc12bTN0S+T/EdXP29PsrT1MW+sckmSpKnIT8fqv9cBN1TVbgBJ1q+q25N8ENilqm4ZR1/fBN5fVecm+TfgY3RWVD4MbF5V9yXZoNX9CPCTqnpPK7skyX+3czsALwbuBi5NchpwCzATeHtV7ZPkO8BbgGOq6h9a7E8F9gBeUFXVNdZQWwEvBdYGFrX+/zRuVV07pP7xwJ7AD5OsBrwa+DsgwGuq6t4kM4HjgNnteg+sqrktrn2B26tq+ySrAz9NcuaQcR7zOrSx/gvYs6ouTbIecE+rPwvYBrgPuDrJocBDwDxgO+BW4MwkbwYuGa68qk4eYX4GY94XYNp6G41UTZIkaVJyJaT/lgK7JpmXZKequv3xdJJkfWCDqjq3FX0D2LkdLwGOTfLXwIOt7LXAh5MsBs4B1gCe1c79uKr+t6ruAU4CXtHKr62qxe14AZ1tT93+CNwLHJnkz+kkMcP5flXd0xKss+kkHwCXDJOAAPwIeFVLIF4PnNdiWxU4IslS4LvAFiOM91rgXe1aLwaeRieh6jbc6/B84MaquhSgqv5YVYPzd1ZV3V5V9wJXApsB2wPnVNXvW71j6bwGI5WPqKoOr6rZVTV72lrrj1ZVkiRp0jEJ6bOq+gWdd8iXAp9O8tEnYJjdgC+1cRYkmU5nFeEtVTWrfT2rqn4+GNbQMNv3+7rKHmLISlq7wd4BOBF4M3D6CPGM1P9dw1bu3OifA/wZnRWR49upDwA3AVvTWQFZbYTxQmeFaPBaN6+qM4eMMdzrkGFiHTTcXGSU8SVJktSYhPRZkk2Bu6vqGOBgYNt26g5g3REbDtHeub81yU6t6J3AuUlWAZ5ZVWcDHwI2ANYBzgDenyQtjm26untNkqe2Zz7eDPx0Oa9lHWD9qvohnW1gs0aounuSNZI8DZgDXLoc3R8PvBvYqcUOsD6dlYqH2/VOa+VD5+4M4O+SrNrifF6StYfEPtzrcBWdZz+2b3XWbQncSC4GXplkw3Qern87cO4o5ZIkSVOSz4T035bAZ5M8DDxA51kHgMOBHyW5sap2SXIk8JWqmj9MH4Pv1u8FfCWdj/a9hs5N+zTgmLZdK8Dnq+q2JP8OHAIsaYnIdcDgP8V8AfAt4LnAt6tqfpIZy3Et6wLfT7JGG+sDI9S7BDiNzvavf6+qG5I8b4y+z6TzzMspVXV/K/tP4MQkb6OzrWtwJWUJ8GCSy4CjgS/Q2Tq2sF3r7+kkV90e8zpU1f3tAfVDW0J2D7DrSAFW1Y1J/qnFEuCHVfV9gJHKJUmSpqJUjbTbRJNBeyB6YVV9fSX1tzcwu6retzL6G6b/g4A7q+oxn5yl4a0+MLMG9jqk32Foglg2b+7YlSRJ6p8FVTV7rEpux5rE2mrGS4BT+h2LJEmStLxcCZEmuCT+kUqSpMnClRBJkiRJE49JiCRJkqSeMgmRJEmS1FMmIZIkSZJ6yiREkiRJUk+ZhEiSJEnqKZMQSZIkST1lEiJJkiSpp0xCJEmSJPWUSYgkSZKknjIJkSRJktRTJiGSJEmSesokRJIkSVJPTe93AJJGt9omz2Vgr0P6HYYmiGXz5vY7BEmSVpgrIZIkSZJ6yiREkiRJUk+ZhExxSTZI8t4ejHNh+75pkhPG2fa6JBsOU/6mJB8eo+2cJKeOL9ox47lzJfVzUJIDV0ZfkiRJk4lJiDYAnvAkpKpe1r7fUFVvXUl9nlJVn1kZfUmSJKl3TEL0GeA5SRYn+WySgSTntZ8vT7ITPPrd/yRvTXJ0Oz46yReTXJjkmiTDJhiD7ZPMSHJ5O947yUlJTk/yyyT/MUqc70+yMMnSJC/oan9YO35OkouSXJrk34asVqyT5IQkVyU5NkmGiW//JFcmWZLk+Fa2TpKvtzGXJHlLV/1PJrmsjblxK9ssyVmt7llJnjVauSRJ0lRlEqIPA7+qqllV9Y/AXwFnVNUsYGtg8XL0MQC8AphLJ6kZj1nAnsCWwJ5JnjlCvVuqalvgy8BwW5i+AHyhqrYHbhhybhvgAGAL4NnAy4dp/2Fgm6raCtivlf0rcHtVbdnKf9LK1wYuqqqtgfOAfVr5YcA3W91jgS+OUT6iJPsmmZ9k/kN33z5WdUmSpEnFJERDXQq8O8lBwJZVdcdytDm5qh6uqiuBjcc53llVdXtV3QtcCWw2Qr2T2vcFwIxhzu8IfLcdf3vIuUuq6vqqephOUjVc+yXAsUn+Gniwle0KfGmwQlXd2g7vBwafM+mOZ8eusb9FJzEbrXxEVXV4Vc2uqtnT1lp/rOqSJEmTikmIHqWqzgN2Bn4LfCvJuwZPdVVbY0iz+7qOH7PVaQzdbR9i5H+75r7lqLMiY+xGJ+HYDliQZDqda6lh6j5QVYPlo8UzXNvRyiVJkqYEkxDdAaw7+EOSzYCbq+oI4GvAtu3UTUlemGQVYI/ehzmmi4DBZzb+cjwN2zU9s6rOBj5E52H9dYAzgfd11XvKGF1d2DX2O4ALxiiXJEmakkxCpriq+l/gp+0h9M8Cc4DFSRbRuan/Qqv6YTpbkH4C3Ph4hloJ4Y7mAOCDSS6h84zKeB6kmAYck2QpsAj4fFXdBnwCeEqbm8uAXcboZ386W9mWAO8E/mGMckmSpCkpj+wqkZ4YSZ4GLKyqkZ73WBljrAXcU1WV5C+Bt1fV7k/UeL20+sDMGtjrkH6HoQli2by5/Q5BkqTRLKiq2WNVGu/eemlckmwKnAMc/AQPtR1wWPv43duA9zzB40mSJOlxciVEmuCS+EcqSZImi+VaCfGZEEmSJEk9ZRIiSZIkqadMQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPWUSIkmSJKmnTEIkSZIk9ZRJiCRJkqSeMgmRJEmS1FPT+x2ApNGttslzGdjrkH6HoQli2by5/Q5BkqQV5kqIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJERPKknmJDl1jDp7JzmsVzFJkiTp0UxCJEmSJPWUSYgmnCRrJzktyWVJLk+yZyu/Lsmnkvwsyfwk2yY5I8mvkuw3TD/bJ1mU5NmjjLVb62/DJBslOTHJpe3r5a3OQUm+leQnSX6ZZJ9WniSfbTEu7YpzTpJzkpyQ5KokxyZJO/eGVnZBki+OtWojSZL0ZORH9Goieh1wQ1XtBpBk/a5zv6mqHZN8HjgaeDmwBnAF8JXBSkleBhwK7F5Vvx5ukCR7AB8E3lBVtyb5NvD5qrogybOAM4AXtupbAS8F1gYWJTkN2BGYBWwNbAhcmuS8Vn8b4EXADcBPgZcnmQ98Fdi5qq5NctxIE5BkX2BfgGnrbTTWfEmSJE0qJiGaiJYCByeZB5xaVed3nTulq846VXUHcEeSe5Ns0M69EDgceG1V3TDCGLsAs1udP7ayXYEt2qIFwHpJ1m3H36+qe4B7kpwN7AC8Ajiuqh4CbkpyLrA98Efgkqq6HiDJYmAGcCdwTVVd2/o8jpZoDFVVh7drYPWBmTXSREmSJE1GbsfShFNVvwC2o5NofDrJR7tO39e+P9x1PPjzYFJ9I3AvndWIkVwDrAs8r6tsFWDHqprVvp7ekhyAoYlAAWFk3bE91GIbrb4kSdKUYRKiCSfJpsDdVXUMcDCw7Ti7uA3YDfhUkjkj1FkG/DnwzSQvamVnAu/rimNWV/3dk6yR5GnAHOBS4DxgzyTTkmwE7AxcMkpcVwHPTjKj/bznOK5JkiTpScMkRBPRlsAlbRvTR4BPjLeDqroJeCPwpSQvGaHO1cA7gO8meQ6wPzA7yZIkVwLdD7tfApwGXAT8e9vm9T1gCXAZ8BPgQ1X1u1Fiugd4L3B6kguAm4Dbx3ttkiRJk12q3G4ujSbJQcCdVXXwSuhrnaq6s31a1peAX1bV50drs/rAzBrY65AVHVpPEsvmze13CJIkjWZBVc0eq5IPpku9tU+SvYDVgEV0Pi1rVPf/7n+88ZQkSU8qroRIE1wS/0glSdJksVwrIT4TIkmSJKmnTEIkSZIk9ZRJiCRJkqSeMgmRJEmS1FMmIZIkSZJ6yiREkiRJUk+ZhEiSJEnqKZMQSZIkST1lEiJJkiSpp0xCJEmSJPWUSYgkSZKknjIJkSRJktRT0/sdgKTRrbbJcxnY65B+h6EJYtm8uf0OQZKkFeZKiCRJkqSeMgmRJEmS1FMmIZIkSZJ66kmThCQ5J8nscdS/LsmGY9S5c8Uje0yfF67sPpdjzD/NTZJ/Xon9HpTkwHG2+WGSDdrxuObiiXg9hhljTpJTl7PujCSXj7P/o5O89fFFJ0mS9OTwpElCJouqelmfQ1hpScjjUVVvqKrb2nG/50KSJEl9MGGTkCRrJzktyWVJLk+yZyvfLsm5SRYkOSPJwJB2qyT5RpJPtJ/fnmRp62PeCGOd3Pq7Ism+Q859ssVwUZKNh2n7yiSL29eiJOsmWSfJWUkWtrF376p/Z/s+p61QnJDkqiTHJkk795kkVyZZkuTg5RzzUe/gJzksyd5D2n0GWLO1O3a0a09y51jX3mzRruOaJPuPNafdK1AjrWyM9/VoqwtfTnJ2i+OVSY5K8vMkR3e1/XKS+a3fj3eVv669BhcAfz5CTC9KckmbuyVJZrZT05Ic0fo8M8marf4+SS5tsZ6YZK2u7nZOcmGL1VURSZI05UzYJAR4HXBDVW1dVS8GTk+yKnAo8Naq2g44CvhkV5vpwLHAL6rqX5JsCswDXgXMArZP8uZhxnpP6282sH+Sp7XytYGLqmpr4Dxgn2HaHgj8fVXNAnYC7gHuBfaoqm2BXYD/O5hgDLENcACwBfBs4OVJngrsAbyoqrYCPrGcY46pqj4M3FNVs6rqHSvh2gFeAPwZsAPwsfYajdbv8ng8MT2Fzuv8AeAHwOeBFwFbJpnV6nykqmYDWwGvTLJVkjWAI4A30pnLTUaIaT/gC23OZwPXt/KZwJeq6kXAbcBbWvlJVbV9i/XnwN909TUAvAKYC3xmuMGS7NsSpvkP3X37iBMlSZI0GU3kJGQpsGuSeUl2qqrbgecDLwZ+nGQx8C/AM7rafBW4vKoGE5PtgXOq6vdV9SCdBGXnYcbaP8llwEXAM+ncWALcDwyuLiwAZgzT9qfA59oqwAZtnACfSrIE+G/g6cBwKwmXVNX1VfUwsLj1/0c6ScyRSf4cuHs5x3y8VuTaAU6rqvuq6hbgZh65zpH6faJi+kFVFZ3fm5uqammb1yu66v1FkoXAIjoJyhZ0kqhrq+qXrf0xI8T0M+Cfk/wfYLOqGkz8rq2qxcPE9OIk5ydZCryjjTfo5Kp6uKquZPjfC6rq8KqaXVWzp621/gghSZIkTU4TNgmpql8A29G5qfx0ko/Subm/or2TP6uqtqyq13Y1uxDYpb27Tas/qiRzgF2BHdu71ouAwfYPtBtTgIcY5h93rKrPAP8/sCZwUZIX0Lnp3AjYrr1zflNXn93u6zp+CJjeEoodgBOBNwOnL+eYD/Lo13O48VbqtY90DWP0+0TFNBjHw0NierjFtDmdFaRXtxWm07r6LcZQVd8G3kRn1emMJK8aMu7QmI4G3ldVWwIf59HX391mzN9RSZKkJ5sJm4S0rVR3V9UxwMHAtsDVwEZJdmx1Vk3S/Q7z14AfAt9NMh24mM62mw2TTAPeDpw7ZKj1gVur6u52M//Sccb5nPau+zxgPp131tcHbq6qB5LsAmw2jv7WAdavqh/S2ao1a5g6w425jM7zGasnWR949QhDPNC1ZWqFrn0UK9LvExXTesBdwO3tWZLXt/KrgM2TPKf9/PbhGid5NnBNVX0ROIXOlq7RrAvc2Ob6HWPUlSRJmlJGend7ItgS+GySh4EHgL+rqvvbg7xfbDfa04FD6Gy5AaCqPtfOfYvOzd8/AWfTecf5h1X1/SHjnA7s17ZOXU1nC9B4HNASjYeAK4Ef0bkB/UGS+XS2WV01jv7WBb7fVnNC5xmHMcesqvuSfAdYAvySzgrCcA4HlrRtSe9hxa59JGPN6WgrDyv6egyrqi5LsojO78o1dLa0UVX3toffT0tyC3ABnS1/Q+0J/HWSB4DfAf9GJ7EZyb/SSYKX0VnNW3dlXIckSdKTQR7Z3SI9sdpq1M3AJlX1QL/jmSxWH5hZA3sd0u8wNEEsmze33yFIkjSaBe2DgEY1Ybdj6UnpCuBIExBJkqSpzZUQaYJL4h+pJEmaLFwJkSRJkjTxmIRIkiRJ6imTEEmSJEk9ZRIiSZIkqadMQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPWUSIkmSJKmnpvc7AEmjW22T5zKw1yH9DkMTxLJ5c/sdgiRJK8yVEEmSJEk9ZRIiSZIkqadMQiRJkiT1lEnIJJHk6CRv7fGYd/ZyvCFjz07yxRHOXZdkw2HKD0pyYDt+wuZrpPEfRz97JzlsZcQkSZI0mfhguiakqpoPzO93HJIkSVr5XAmZgJK8K8mSJJcl+VbXqZ2TXJjkmsF3+ZPMSXJqV9vDkuzdjq9L8vEkC5MsTfKCVr5Rkh+38q8mWTbSO/tJPtniuCjJxq3sUasMgysmSdZJclbXeLuP0OedSeYlWZDkv5PskOScdl1vGnpdSZ6W5Mwki5J8FUhXXx9JcnWS/waeP8J42yU5t413RpKBYeq8Lcnl7VrPa2XTkhzcrmVJkvd3NXn/MPP61CQnt7oXJdlqtHJJkqSpyiRkgknyIuAjwKuqamvgH7pODwCvAOYCn1nOLm+pqm2BLwMHtrKPAT9p5d8DnjVC27WBi1oc5wH7jDHWvcAerd9dgP+bJMPUWxs4p6q2A+4APgG8BtgD+Ldh6n8MuKCqtgFOGYw3yXbAXwLbAH8ObD+0YZJVgUOBt7bxjgI+OcwYHwX+rF3rm1rZvsDmwDZVtRVwbFf94eb148CiVvefgW+OUT6iJPsmmZ9k/kN33z5WdUmSpEnF7VgTz6uAE6rqFoCq+kPXuZOr6mHgysFVieVwUvu+gM6NOnQSmT1a/6cnuXWEtvcDg6ssC+gkCqMJ8KkkOwMPA08HNgZ+N0y/p7fjpcB9VfVAkqXAjGH63Xkw9qo6rSvenYDvVdXdAElOGabt84EXAz9u+dA04MZh6v0UODrJd3hkznYFvlJVD7axu1+Lkeb1La3uT9oKzvqjlI+oqg4HDgdYfWBmjVZXkiRpsjEJmXgCjHTTed+QegAP8ugVrTVGaPMQj7zew61ODOeBqhqMpbv9n8ZsKx2rtfJ3ABsB27Wk4rph4hna78ODMVbVw0lG+p0caU7GukEPcEVV7ThaparaL8lLgN2AxUlmsXyvxVjzWqOUS5IkTUlux5p4zgL+IsnToPM8wRj1lwFbJFm9vbv+6uUY4wLgL1r/rwWeMs4YrwO2a8e7A6u24/WBm1sCsguw2Tj7Hcl5dBIckryeR+I9D9gjyZpJ1gXeOEzbq4GNkuzY2q/atrw9SpLnVNXFVfVR4BbgmcCZwH6DidFyvBbdcc6hs2Xrj6OUS5IkTUmuhEwwVXVFkk8C5yZ5CFgE7D1K/d+0LURLgF+2+mP5OHBckj2Bc+lsT7pjHGEeAXw/ySV0kqa7WvmxwA+SzAcWA1eNo8/liXdhi/fXAFW1MMl/tbGWAecPbVhV97eH6L/YkrTpwCHAFUOqfjbJTDqrFmcBlwGXA88DliR5oF33aB+pexDw9SRLgLuBvcYolyRJmpLyyK4YTRVJVgceqqoH2wrBl6tqVr/j0vBWH5hZA3sd0u8wNEEsmze33yFIkjSaBVU1e6xKroRMTc8CvpNkFToPiY/1qVeSJEnSSuNKiDTBJfGPVJIkTRbLtRLig+mSJEmSesokRJIkSVJPmYRIkiRJ6imTEEmSJEk9ZRIiSZIkqadMQiRJkiT1lEmIJEmSpJ4yCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPTW93wFIGt1qmzyXgb0O6XcYmiCWzZvb7xAkSVphroRIkiRJ6imTEEmSJEk9ZRIiSZIkqadMQiapJDOSXN7D8Y5O8tZhyjdNcsITPPbeSQ4b4dwPk2wwRvtzksxeifEclOTAldTXnSujH0mSpMnEB9O1QqrqBuAxyUkPx39Dv8aWJEnS4+NKyOQ2LckRSa5IcmaSNQGS7JPk0iSXJTkxyVqt/OgkX0xyYZJrhlvZaPXelWRJa/+trlM7D23bvSLTjs9PsrB9vayVz2mrESckuSrJsUnSzr2hlV3QYjt1hGvdNMnpSX6Z5D+6Yr0uyYbt+F9bXz9OctyQ1Yq3JbkkyS+S7DTMNQ8kOS/J4iSXD9ZJ8rp2LZclOauryRbtmq5Jsn9XPx9s7S9PcsBY5ZIkSVORKyGT20zg7VW1T5LvAG8BjgFOqqojAJJ8Avgb4NDWZgB4BfAC4BTgUVupkrwI+Ajw8qq6JclTu06P2ha4GXhNVd2bZCZwHDC4DWob4EXADcBPgZcnmQ98Fdi5qq5Nctwo1zqr9XEfcHWSQ6vqN11xz27Xvw2d3+uFwIKu9tOraockbwA+Buw6pP+/As6oqk8mmQaslWQj4Iiu+Lrn4gXALsC6LZ4vA1sB7wZeAgS4OMm5dJL9x5RX1aKRLjbJvsC+ANPW22iUaZEkSZp8TEImt2uranE7XgDMaMcvbsnHBsA6wBldbU6uqoeBK5NsPEyfrwJOqKpbAKrqD+NouypwWJJZwEPA87rOXVJV1wMkWdxivRO4pqqubXWOo914D+Osqrq9tb8S2Az4Tdf5VwDfr6p7Wp0fDGl/UvvePU/dLgWOSrJqu87FSeYA5w3GN2QuTquq+4D7ktwMbNxi+F5V3dViOAnYiU7iMVz5iElIVR0OHA6w+sDMGqmeJEnSZOR2rMntvq7jh3gkqTwaeF9VbQl8HFhjhDYZps8AI930jtX2A8BNwNZ0VkBWGyPW4foYyUjXOlo8w7Ufri1VdR6wM/Bb4FtJ3sXyz8VY1zOe65QkSXrSMwl5cloXuLG9q/+OcbY9C/iLJE8DGLIFaSzrAze21ZJ3AtPGqH8V8OwkM9rPe44v1Ee5AHhjkjWSrAPsNp7GSTYDbm7b2L4GbAv8DHhlks1bnbHm4jzgzUnWSrI2sAdw/ijlkiRJU5LbsZ6c/hW4GFgGLKWTlCyXqroiySeBc5M8RGfL0N7L2fw/gROTvA04G7hrjLHuSfJe4PQktwCXLG+cw/R1aZJTgMvoXPd84PZxdDEH+MckD9DZJvauqvp9ezbjpCSr0J55GSWGhUmO5pHrOHLwuY+RyiVJkqaiVLndXP2TZJ2qurN9WtaXgF9W1edXsK+16Kw+7FtVC1dmvP2w+sDMGtjrkH6HoQli2by5/Q5BkqTRLKiqMf99Nrdjqd/2aQ+qX0FnO9dXV6Cvw1tfC4ETnwwJiCRJ0pORKyHSBJfEP1JJkjRZuBIiSZIkaeIxCZEkSZLUUyYhkiRJknrKJESSJElST5mESJIkSeopkxBJkiRJPeW/mC5NfHcCV/c7iAlgQ+CWfgcxATgPj3AuOpyHRzgXHc7DI5yLjl7Ow2bLU8kkRJr4rl6ez9t+sksy33lwHro5Fx3OwyOciw7n4RHORcdEnAe3Y0mSJEnqKZMQSZIkST1lEiJNfIf3O4AJwnnocB4e4Vx0OA+PcC46nIdHOBcdE24eUlX9jkGSJEnSFOJKiCRJkqSeMgmRJEmS1FMmIdIEleR1Sa5O8j9JPtzvePolyVFJbk5yeb9j6ackz0xydpKfJ7kiyT/0O6Z+SLJGkkuSXNbm4eP9jqnfkkxLsijJqf2OpV+SXJdkaZLFSeb3O55+SrJBkhOSXNX+e7Fjv2PqtSTPb78Lg19/THJAv+PqhyQfaP+tvDzJcUnW6HdMg3wmRJqAkkwDfgG8BrgeuBR4e1Vd2dfA+iDJznT+wcZvVtWL+x1PvyQZAAaqamGSdYEFwJun2u9EkgBrV9WdSVYFLgD+oaou6nNofZPkg8BsYL2qmtvvePohyXXA7Kqa8v8oXZJvAOdX1ZFJVgPWqqrb+h1Xv7T/n/4WeElVLet3PL2U5Ol0/hu5RVXdk+Q7wA+r6uj+RtbhSog0Me0A/E9VXVNV9wPHA7v3Oaa+qKrzgD/0O45+q6obq2phO74D+Dnw9P5G1XvVcWf7cdX2NWXfTUvyDGA34Mh+x6L+S7IesDPwNYCqun8qJyDNq4FfTbUEpMt0YM0k04G1gBv6HM+fmIRIE9PTgd90/Xw9U/CGU8NLMgPYBri4v5H0R9t+tBi4GfhxVU3JeWgOAT4EPNzvQPqsgDOTLEiyb7+D6aNnA78Hvt626B2ZZO1+B9Vnfwkc1+8g+qGqfgscDPwauBG4varO7G9UjzAJkSamDFM2Zd/t1SOSrAOcCBxQVX/sdzz9bP3qZQAAAgZJREFUUFUPVdUs4BnADkmm5Da9JHOBm6tqQb9jmQBeXlXbAq8H/r5t45yKpgPbAl+uqm2Au4Cp/EzhasCbgO/2O5Z+SPIUOrsoNgc2BdZO8tf9jeoRJiHSxHQ98Myun5/BBFpCVX+0ZyBOBI6tqpP6HU+/tW0m5wCv63Mo/fJy4E3teYjjgVclOaa/IfVHVd3Qvt8MfI/Oltap6Hrg+q7VwRPoJCVT1euBhVV1U78D6ZNdgWur6vdV9QBwEvCyPsf0JyYh0sR0KTAzyebtnZy/BE7pc0zqo/ZA9teAn1fV5/odT78k2SjJBu14TTr/k72qv1H1R1X9U1U9o6pm0PlvxE+qasK8y9krSdZuH9ZA23r0WmBKfppeVf0O+E2S57eiVwNT6sMrhng7U3QrVvNr4KVJ1mr/D3k1necJJ4Tp/Q5A0mNV1YNJ3gecAUwDjqqqK/ocVl8kOQ6YA2yY5HrgY1X1tf5G1RcvB94JLG3PQwD8c1X9sI8x9cMA8I32iTerAN+pqin70bQCYGPge517LKYD366q0/sbUl+9Hzi2vYF1DfDuPsfTF0nWovMJk3/b71j6paouTnICsBB4EFgEHN7fqB7hR/RKkiRJ6im3Y0mSJEnqKZMQSZIkST1lEiJJkiSpp0xCJEmSJPWUSYgkSZKknjIJkSRJktRTJiGSJEmSeur/AU2Xzwecl+wLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x18c06c23f60>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "top = 20\n",
    "fig, ax = plt.subplots(figsize=(10, 10))\n",
    "top_schools = df_clean['school'].value_counts()[:top].plot.barh(ax=ax)\n",
    "ax.invert_yaxis()\n",
    "ax.set_facecolor('black')\n",
    "ax.set_title(f'Top {top} Schools')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "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.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}