openbus_18_haifa_tlv_train_hishmul_boosts.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Imports and config"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T19:43:11.675924Z",
     "start_time": "2018-08-15T19:43:11.269704Z"
    }
   },
   "outputs": [],
   "source": [
    "# Put these at the top of every notebook, to get automatic reloading and inline plotting\n",
    "%reload_ext autoreload\n",
    "%autoreload 2\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T20:33:16.173192Z",
     "start_time": "2018-08-15T20:33:16.001391Z"
    }
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import partridge as ptg\n",
    "import datetime\n",
    "import os\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "#import gtfs_utils as gu\n",
    "#import gtfstk\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T19:43:13.019294Z",
     "start_time": "2018-08-15T19:43:12.503790Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "RendererRegistry.enable('notebook')"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import altair as alt\n",
    "alt.renderers.enable('notebook')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T19:43:16.612411Z",
     "start_time": "2018-08-15T19:43:13.331720Z"
    }
   },
   "outputs": [],
   "source": [
    "import glob\n",
    "rs = (pd.concat([pd.read_pickle(f, compression='gzip') \n",
    "                 for f in glob.glob('data\\\\gtfs_stats_MOD_ffill\\\\2018-07*_route_stats.pkl.gz')] + \n",
    "                [pd.read_pickle(f, compression='gzip') \n",
    "                 for f in glob.glob('data\\\\gtfs_stats_MOD_ffill\\\\2018-08*_route_stats.pkl.gz')])\n",
    "      .assign(date = lambda x: pd.to_datetime(x.date))\n",
    "      .assign(dayofweek = lambda x: x['date'].dt.strftime('%A')))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T19:43:18.158705Z",
     "start_time": "2018-08-15T19:43:17.986872Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(229087, 37)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rs.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T19:43:18.924154Z",
     "start_time": "2018-08-15T19:43:18.689863Z"
    }
   },
   "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>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>route_id</th>\n",
       "      <td>1</td>\n",
       "      <td>10</td>\n",
       "      <td>10000</td>\n",
       "      <td>10001</td>\n",
       "      <td>10002</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>route_short_name</th>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>91</td>\n",
       "      <td>91</td>\n",
       "      <td>92</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>agency_id</th>\n",
       "      <td>25</td>\n",
       "      <td>25</td>\n",
       "      <td>14</td>\n",
       "      <td>14</td>\n",
       "      <td>14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>agency_name</th>\n",
       "      <td>אפיקים</td>\n",
       "      <td>אפיקים</td>\n",
       "      <td>נתיב אקספרס</td>\n",
       "      <td>נתיב אקספרס</td>\n",
       "      <td>נתיב אקספרס</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>route_long_name</th>\n",
       "      <td>ת. רכבת יבנה מערב-יבנה&lt;-&gt;ת. רכבת יבנה מזרח-יבנ...</td>\n",
       "      <td>ת. רכבת יבנה מזרח-יבנה&lt;-&gt;ת. רכבת יבנה מערב-יבנ...</td>\n",
       "      <td>קופת חולים/הנשיא-צפת&lt;-&gt;דרך חטיבת יפתח-צפת-11</td>\n",
       "      <td>האר''י/ישראל ב''ק-צפת&lt;-&gt;שפרינצק/קרן היסוד-צפת-2#</td>\n",
       "      <td>געתון/מירון-כפר ורדים&lt;-&gt;בית ספר כפר ורדים-כפר ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>route_type</th>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_trips</th>\n",
       "      <td>26</td>\n",
       "      <td>35</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_trip_starts</th>\n",
       "      <td>26</td>\n",
       "      <td>35</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_trip_ends</th>\n",
       "      <td>26</td>\n",
       "      <td>35</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>is_loop</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>is_bidirectional</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>start_time</th>\n",
       "      <td>05:10:00</td>\n",
       "      <td>06:05:00</td>\n",
       "      <td>08:00:00</td>\n",
       "      <td>12:00:00</td>\n",
       "      <td>07:30:00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>end_time</th>\n",
       "      <td>21:54:44</td>\n",
       "      <td>22:48:36</td>\n",
       "      <td>08:17:31</td>\n",
       "      <td>12:36:42</td>\n",
       "      <td>07:53:03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max_headway</th>\n",
       "      <td>65</td>\n",
       "      <td>50</td>\n",
       "      <td>NaN</td>\n",
       "      <td>15</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min_headway</th>\n",
       "      <td>20</td>\n",
       "      <td>10</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean_headway</th>\n",
       "      <td>37.8125</td>\n",
       "      <td>28</td>\n",
       "      <td>NaN</td>\n",
       "      <td>5</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>peak_num_trips</th>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>peak_start_time</th>\n",
       "      <td>05:10:00</td>\n",
       "      <td>15:55:00</td>\n",
       "      <td>08:00:00</td>\n",
       "      <td>12:15:00</td>\n",
       "      <td>07:30:00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>peak_end_time</th>\n",
       "      <td>05:29:44</td>\n",
       "      <td>16:03:36</td>\n",
       "      <td>08:17:31</td>\n",
       "      <td>12:21:42</td>\n",
       "      <td>07:53:03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>service_distance</th>\n",
       "      <td>182416</td>\n",
       "      <td>261100</td>\n",
       "      <td>6326</td>\n",
       "      <td>33504</td>\n",
       "      <td>7175</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>service_duration</th>\n",
       "      <td>8.55111</td>\n",
       "      <td>10.85</td>\n",
       "      <td>0.291944</td>\n",
       "      <td>1.44667</td>\n",
       "      <td>0.384167</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>service_speed</th>\n",
       "      <td>21.3324</td>\n",
       "      <td>24.0645</td>\n",
       "      <td>21.6685</td>\n",
       "      <td>23.1594</td>\n",
       "      <td>18.6768</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean_trip_distance</th>\n",
       "      <td>7016</td>\n",
       "      <td>7460</td>\n",
       "      <td>6326</td>\n",
       "      <td>8376</td>\n",
       "      <td>7175</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean_trip_duration</th>\n",
       "      <td>0.328889</td>\n",
       "      <td>0.31</td>\n",
       "      <td>0.291944</td>\n",
       "      <td>0.361667</td>\n",
       "      <td>0.384167</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>start_stop_id</th>\n",
       "      <td>38725</td>\n",
       "      <td>15745</td>\n",
       "      <td>18423</td>\n",
       "      <td>37179</td>\n",
       "      <td>7270</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>end_stop_id</th>\n",
       "      <td>15745</td>\n",
       "      <td>38725</td>\n",
       "      <td>38800</td>\n",
       "      <td>18420</td>\n",
       "      <td>17976</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>start_stop_lat</th>\n",
       "      <td>31.8907</td>\n",
       "      <td>31.8621</td>\n",
       "      <td>32.9605</td>\n",
       "      <td>32.9719</td>\n",
       "      <td>32.9928</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>start_stop_lon</th>\n",
       "      <td>34.7312</td>\n",
       "      <td>34.7442</td>\n",
       "      <td>35.4955</td>\n",
       "      <td>35.4933</td>\n",
       "      <td>35.258</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>end_stop_lat</th>\n",
       "      <td>31.8621</td>\n",
       "      <td>31.8907</td>\n",
       "      <td>32.9684</td>\n",
       "      <td>32.9585</td>\n",
       "      <td>32.9961</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>end_stop_lon</th>\n",
       "      <td>34.7442</td>\n",
       "      <td>34.7312</td>\n",
       "      <td>35.4944</td>\n",
       "      <td>35.4942</td>\n",
       "      <td>35.2771</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_stops</th>\n",
       "      <td>19</td>\n",
       "      <td>24</td>\n",
       "      <td>23</td>\n",
       "      <td>26</td>\n",
       "      <td>16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>start_zone</th>\n",
       "      <td>רחובות</td>\n",
       "      <td>רחובות</td>\n",
       "      <td>גליל עליון ורמת הגולן</td>\n",
       "      <td>גליל עליון ורמת הגולן</td>\n",
       "      <td>נהריה</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>end_zone</th>\n",
       "      <td>רחובות</td>\n",
       "      <td>רחובות</td>\n",
       "      <td>גליל עליון ורמת הגולן</td>\n",
       "      <td>גליל עליון ורמת הגולן</td>\n",
       "      <td>נהריה</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_zones</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>num_zones_missing</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>date</th>\n",
       "      <td>2018-07-01 00:00:00</td>\n",
       "      <td>2018-07-01 00:00:00</td>\n",
       "      <td>2018-07-01 00:00:00</td>\n",
       "      <td>2018-07-01 00:00:00</td>\n",
       "      <td>2018-07-01 00:00:00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>dayofweek</th>\n",
       "      <td>Sunday</td>\n",
       "      <td>Sunday</td>\n",
       "      <td>Sunday</td>\n",
       "      <td>Sunday</td>\n",
       "      <td>Sunday</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                                    0  \\\n",
       "route_id                                                            1   \n",
       "route_short_name                                                    1   \n",
       "agency_id                                                          25   \n",
       "agency_name                                                    אפיקים   \n",
       "route_long_name     ת. רכבת יבנה מערב-יבנה<->ת. רכבת יבנה מזרח-יבנ...   \n",
       "route_type                                                          3   \n",
       "num_trips                                                          26   \n",
       "num_trip_starts                                                    26   \n",
       "num_trip_ends                                                      26   \n",
       "is_loop                                                             0   \n",
       "is_bidirectional                                                    0   \n",
       "start_time                                                   05:10:00   \n",
       "end_time                                                     21:54:44   \n",
       "max_headway                                                        65   \n",
       "min_headway                                                        20   \n",
       "mean_headway                                                  37.8125   \n",
       "peak_num_trips                                                      1   \n",
       "peak_start_time                                              05:10:00   \n",
       "peak_end_time                                                05:29:44   \n",
       "service_distance                                               182416   \n",
       "service_duration                                              8.55111   \n",
       "service_speed                                                 21.3324   \n",
       "mean_trip_distance                                               7016   \n",
       "mean_trip_duration                                           0.328889   \n",
       "start_stop_id                                                   38725   \n",
       "end_stop_id                                                     15745   \n",
       "start_stop_lat                                                31.8907   \n",
       "start_stop_lon                                                34.7312   \n",
       "end_stop_lat                                                  31.8621   \n",
       "end_stop_lon                                                  34.7442   \n",
       "num_stops                                                          19   \n",
       "start_zone                                                     רחובות   \n",
       "end_zone                                                       רחובות   \n",
       "num_zones                                                           1   \n",
       "num_zones_missing                                                   0   \n",
       "date                                              2018-07-01 00:00:00   \n",
       "dayofweek                                                      Sunday   \n",
       "\n",
       "                                                                    1  \\\n",
       "route_id                                                           10   \n",
       "route_short_name                                                    4   \n",
       "agency_id                                                          25   \n",
       "agency_name                                                    אפיקים   \n",
       "route_long_name     ת. רכבת יבנה מזרח-יבנה<->ת. רכבת יבנה מערב-יבנ...   \n",
       "route_type                                                          3   \n",
       "num_trips                                                          35   \n",
       "num_trip_starts                                                    35   \n",
       "num_trip_ends                                                      35   \n",
       "is_loop                                                             0   \n",
       "is_bidirectional                                                    0   \n",
       "start_time                                                   06:05:00   \n",
       "end_time                                                     22:48:36   \n",
       "max_headway                                                        50   \n",
       "min_headway                                                        10   \n",
       "mean_headway                                                       28   \n",
       "peak_num_trips                                                      2   \n",
       "peak_start_time                                              15:55:00   \n",
       "peak_end_time                                                16:03:36   \n",
       "service_distance                                               261100   \n",
       "service_duration                                                10.85   \n",
       "service_speed                                                 24.0645   \n",
       "mean_trip_distance                                               7460   \n",
       "mean_trip_duration                                               0.31   \n",
       "start_stop_id                                                   15745   \n",
       "end_stop_id                                                     38725   \n",
       "start_stop_lat                                                31.8621   \n",
       "start_stop_lon                                                34.7442   \n",
       "end_stop_lat                                                  31.8907   \n",
       "end_stop_lon                                                  34.7312   \n",
       "num_stops                                                          24   \n",
       "start_zone                                                     רחובות   \n",
       "end_zone                                                       רחובות   \n",
       "num_zones                                                           1   \n",
       "num_zones_missing                                                   0   \n",
       "date                                              2018-07-01 00:00:00   \n",
       "dayofweek                                                      Sunday   \n",
       "\n",
       "                                                               2  \\\n",
       "route_id                                                   10000   \n",
       "route_short_name                                              91   \n",
       "agency_id                                                     14   \n",
       "agency_name                                          נתיב אקספרס   \n",
       "route_long_name     קופת חולים/הנשיא-צפת<->דרך חטיבת יפתח-צפת-11   \n",
       "route_type                                                     3   \n",
       "num_trips                                                      1   \n",
       "num_trip_starts                                                1   \n",
       "num_trip_ends                                                  1   \n",
       "is_loop                                                        0   \n",
       "is_bidirectional                                               0   \n",
       "start_time                                              08:00:00   \n",
       "end_time                                                08:17:31   \n",
       "max_headway                                                  NaN   \n",
       "min_headway                                                  NaN   \n",
       "mean_headway                                                 NaN   \n",
       "peak_num_trips                                                 1   \n",
       "peak_start_time                                         08:00:00   \n",
       "peak_end_time                                           08:17:31   \n",
       "service_distance                                            6326   \n",
       "service_duration                                        0.291944   \n",
       "service_speed                                            21.6685   \n",
       "mean_trip_distance                                          6326   \n",
       "mean_trip_duration                                      0.291944   \n",
       "start_stop_id                                              18423   \n",
       "end_stop_id                                                38800   \n",
       "start_stop_lat                                           32.9605   \n",
       "start_stop_lon                                           35.4955   \n",
       "end_stop_lat                                             32.9684   \n",
       "end_stop_lon                                             35.4944   \n",
       "num_stops                                                     23   \n",
       "start_zone                                 גליל עליון ורמת הגולן   \n",
       "end_zone                                   גליל עליון ורמת הגולן   \n",
       "num_zones                                                      1   \n",
       "num_zones_missing                                              0   \n",
       "date                                         2018-07-01 00:00:00   \n",
       "dayofweek                                                 Sunday   \n",
       "\n",
       "                                                                   3  \\\n",
       "route_id                                                       10001   \n",
       "route_short_name                                                  91   \n",
       "agency_id                                                         14   \n",
       "agency_name                                              נתיב אקספרס   \n",
       "route_long_name     האר''י/ישראל ב''ק-צפת<->שפרינצק/קרן היסוד-צפת-2#   \n",
       "route_type                                                         3   \n",
       "num_trips                                                          4   \n",
       "num_trip_starts                                                    4   \n",
       "num_trip_ends                                                      4   \n",
       "is_loop                                                            0   \n",
       "is_bidirectional                                                   0   \n",
       "start_time                                                  12:00:00   \n",
       "end_time                                                    12:36:42   \n",
       "max_headway                                                       15   \n",
       "min_headway                                                        0   \n",
       "mean_headway                                                       5   \n",
       "peak_num_trips                                                     4   \n",
       "peak_start_time                                             12:15:00   \n",
       "peak_end_time                                               12:21:42   \n",
       "service_distance                                               33504   \n",
       "service_duration                                             1.44667   \n",
       "service_speed                                                23.1594   \n",
       "mean_trip_distance                                              8376   \n",
       "mean_trip_duration                                          0.361667   \n",
       "start_stop_id                                                  37179   \n",
       "end_stop_id                                                    18420   \n",
       "start_stop_lat                                               32.9719   \n",
       "start_stop_lon                                               35.4933   \n",
       "end_stop_lat                                                 32.9585   \n",
       "end_stop_lon                                                 35.4942   \n",
       "num_stops                                                         26   \n",
       "start_zone                                     גליל עליון ורמת הגולן   \n",
       "end_zone                                       גליל עליון ורמת הגולן   \n",
       "num_zones                                                          1   \n",
       "num_zones_missing                                                  0   \n",
       "date                                             2018-07-01 00:00:00   \n",
       "dayofweek                                                     Sunday   \n",
       "\n",
       "                                                                    4  \n",
       "route_id                                                        10002  \n",
       "route_short_name                                                   92  \n",
       "agency_id                                                          14  \n",
       "agency_name                                               נתיב אקספרס  \n",
       "route_long_name     געתון/מירון-כפר ורדים<->בית ספר כפר ורדים-כפר ...  \n",
       "route_type                                                          3  \n",
       "num_trips                                                           1  \n",
       "num_trip_starts                                                     1  \n",
       "num_trip_ends                                                       1  \n",
       "is_loop                                                             0  \n",
       "is_bidirectional                                                    0  \n",
       "start_time                                                   07:30:00  \n",
       "end_time                                                     07:53:03  \n",
       "max_headway                                                       NaN  \n",
       "min_headway                                                       NaN  \n",
       "mean_headway                                                      NaN  \n",
       "peak_num_trips                                                      1  \n",
       "peak_start_time                                              07:30:00  \n",
       "peak_end_time                                                07:53:03  \n",
       "service_distance                                                 7175  \n",
       "service_duration                                             0.384167  \n",
       "service_speed                                                 18.6768  \n",
       "mean_trip_distance                                               7175  \n",
       "mean_trip_duration                                           0.384167  \n",
       "start_stop_id                                                    7270  \n",
       "end_stop_id                                                     17976  \n",
       "start_stop_lat                                                32.9928  \n",
       "start_stop_lon                                                 35.258  \n",
       "end_stop_lat                                                  32.9961  \n",
       "end_stop_lon                                                  35.2771  \n",
       "num_stops                                                          16  \n",
       "start_zone                                                      נהריה  \n",
       "end_zone                                                        נהריה  \n",
       "num_zones                                                           1  \n",
       "num_zones_missing                                                   0  \n",
       "date                                              2018-07-01 00:00:00  \n",
       "dayofweek                                                      Sunday  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rs.head().T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T20:59:53.370510Z",
     "start_time": "2018-08-15T20:59:53.120567Z"
    }
   },
   "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>route_id</th>\n",
       "      <th>route_short_name</th>\n",
       "      <th>agency_id</th>\n",
       "      <th>agency_name</th>\n",
       "      <th>route_long_name</th>\n",
       "      <th>route_type</th>\n",
       "      <th>num_trips</th>\n",
       "      <th>num_trip_starts</th>\n",
       "      <th>num_trip_ends</th>\n",
       "      <th>is_loop</th>\n",
       "      <th>...</th>\n",
       "      <th>start_stop_lon</th>\n",
       "      <th>end_stop_lat</th>\n",
       "      <th>end_stop_lon</th>\n",
       "      <th>num_stops</th>\n",
       "      <th>start_zone</th>\n",
       "      <th>end_zone</th>\n",
       "      <th>num_zones</th>\n",
       "      <th>num_zones_missing</th>\n",
       "      <th>date</th>\n",
       "      <th>dayofweek</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>4983</th>\n",
       "      <td>7189</td>\n",
       "      <td>910</td>\n",
       "      <td>3</td>\n",
       "      <td>אגד</td>\n",
       "      <td>ת.מרכזית ת''א ק.7/רציפים-תל אביב יפו&lt;-&gt;ת. מרכז...</td>\n",
       "      <td>3</td>\n",
       "      <td>14</td>\n",
       "      <td>14</td>\n",
       "      <td>14</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>34.780031</td>\n",
       "      <td>32.793368</td>\n",
       "      <td>34.959587</td>\n",
       "      <td>22</td>\n",
       "      <td>גוש דן</td>\n",
       "      <td>סובב חיפה</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>2018-07-01</td>\n",
       "      <td>Sunday</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4985</th>\n",
       "      <td>7190</td>\n",
       "      <td>910</td>\n",
       "      <td>3</td>\n",
       "      <td>אגד</td>\n",
       "      <td>ת.מרכזית ת''א ק.7/רציפים-תל אביב יפו&lt;-&gt;ת. רכבת...</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>34.780031</td>\n",
       "      <td>32.793065</td>\n",
       "      <td>34.957647</td>\n",
       "      <td>22</td>\n",
       "      <td>גוש דן</td>\n",
       "      <td>סובב חיפה</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>2018-07-01</td>\n",
       "      <td>Sunday</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4986</th>\n",
       "      <td>7192</td>\n",
       "      <td>910</td>\n",
       "      <td>3</td>\n",
       "      <td>אגד</td>\n",
       "      <td>ת. מרכזית חוף הכרמל/רציפים בינעירוני-חיפה&lt;-&gt;ת....</td>\n",
       "      <td>3</td>\n",
       "      <td>19</td>\n",
       "      <td>19</td>\n",
       "      <td>18</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>34.958985</td>\n",
       "      <td>32.055512</td>\n",
       "      <td>34.780335</td>\n",
       "      <td>23</td>\n",
       "      <td>סובב חיפה</td>\n",
       "      <td>גוש דן</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>2018-07-01</td>\n",
       "      <td>Sunday</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3 rows × 37 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     route_id route_short_name agency_id agency_name  \\\n",
       "4983     7189              910         3         אגד   \n",
       "4985     7190              910         3         אגד   \n",
       "4986     7192              910         3         אגד   \n",
       "\n",
       "                                        route_long_name  route_type  \\\n",
       "4983  ת.מרכזית ת''א ק.7/רציפים-תל אביב יפו<->ת. מרכז...           3   \n",
       "4985  ת.מרכזית ת''א ק.7/רציפים-תל אביב יפו<->ת. רכבת...           3   \n",
       "4986  ת. מרכזית חוף הכרמל/רציפים בינעירוני-חיפה<->ת....           3   \n",
       "\n",
       "      num_trips  num_trip_starts  num_trip_ends  is_loop    ...      \\\n",
       "4983         14               14             14        0    ...       \n",
       "4985          1                1              0        0    ...       \n",
       "4986         19               19             18        0    ...       \n",
       "\n",
       "      start_stop_lon end_stop_lat end_stop_lon  num_stops  start_zone  \\\n",
       "4983       34.780031    32.793368    34.959587         22      גוש דן   \n",
       "4985       34.780031    32.793065    34.957647         22      גוש דן   \n",
       "4986       34.958985    32.055512    34.780335         23   סובב חיפה   \n",
       "\n",
       "       end_zone  num_zones num_zones_missing       date  dayofweek  \n",
       "4983  סובב חיפה          6                 0 2018-07-01     Sunday  \n",
       "4985  סובב חיפה          6                 0 2018-07-01     Sunday  \n",
       "4986     גוש דן          6                 0 2018-07-01     Sunday  \n",
       "\n",
       "[3 rows x 37 columns]"
      ]
     },
     "execution_count": 71,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "route_sn = ['412', '980', '111', '460', '461', '701', '605']\n",
    "route_sn = ['910']\n",
    "rs[rs.route_short_name.isin(route_sn)].drop_duplicates(subset=['route_short_name', 'route_long_name', 'agency_id', 'start_zone', 'end_zone'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T21:00:18.880124Z",
     "start_time": "2018-08-15T21:00:18.614565Z"
    }
   },
   "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>route_short_name</th>\n",
       "      <th>111</th>\n",
       "      <th>412</th>\n",
       "      <th>460</th>\n",
       "      <th>461</th>\n",
       "      <th>605</th>\n",
       "      <th>701</th>\n",
       "      <th>900</th>\n",
       "      <th>910</th>\n",
       "      <th>980</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>date</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2018-07-01</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>139.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>34.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-02</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>139.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-03</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>139.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-04</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>139.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-05</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>139.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-08</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-09</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-10</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-11</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-12</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-15</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-16</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-17</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-18</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-19</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-22</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-23</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-24</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-25</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-26</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-29</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-30</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-07-31</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-01</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-02</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-05</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-07</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-08</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-09</th>\n",
       "      <td>90.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-12</th>\n",
       "      <td>90.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>42.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-13</th>\n",
       "      <td>90.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>42.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-14</th>\n",
       "      <td>90.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>42.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2018-08-15</th>\n",
       "      <td>90.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>159.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>42.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "route_short_name   111   412   460   461    605   701  900   910   980\n",
       "date                                                                  \n",
       "2018-07-01        90.0  40.0  32.0  60.0  139.0   0.0  0.0  34.0   0.0\n",
       "2018-07-02        90.0  40.0  32.0  60.0  139.0   0.0  0.0  31.0   0.0\n",
       "2018-07-03        90.0  40.0  32.0  60.0  139.0   0.0  0.0  31.0   0.0\n",
       "2018-07-04        90.0  40.0  32.0  60.0  139.0   0.0  0.0  31.0   0.0\n",
       "2018-07-05        90.0  40.0  32.0  60.0  139.0   0.0  0.0  31.0   0.0\n",
       "2018-07-08        90.0  40.0  32.0  60.0  160.0   0.0  0.0  36.0   0.0\n",
       "2018-07-09        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-10        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-11        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-12        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-15        90.0  40.0  32.0  60.0  160.0   0.0  0.0  36.0   0.0\n",
       "2018-07-16        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-17        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-18        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-19        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-22        90.0  40.0  32.0  60.0  160.0   0.0  0.0  36.0   0.0\n",
       "2018-07-23        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-24        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-25        90.0  40.0  32.0  60.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-26        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-29        90.0  40.0  32.0  59.0  160.0   0.0  0.0  36.0   0.0\n",
       "2018-07-30        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-07-31        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-01        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-02        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-05        90.0  40.0  32.0  59.0  160.0   0.0  0.0  36.0   0.0\n",
       "2018-08-07        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-08        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-09        90.0  40.0  32.0  59.0  159.0   0.0  0.0  31.0   0.0\n",
       "2018-08-12        90.0  81.0  32.0  60.0  160.0  42.0  0.0  36.0  42.0\n",
       "2018-08-13        90.0  81.0  32.0  60.0  159.0  42.0  0.0  31.0  42.0\n",
       "2018-08-14        90.0  81.0  32.0  60.0  159.0  42.0  0.0  31.0  42.0\n",
       "2018-08-15        90.0  81.0  32.0  60.0  159.0  42.0  0.0  31.0  42.0"
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "route_ids = ['10953', '10954', '11178', '11179', '11941', '11942',\n",
    "             '15544', '15545', '21331', '21332', '23451', '23452',\n",
    "             '711', '712', '713', '714', '715', '21994',\n",
    "             '15383', '19733', '7189', '7190', '7192']\n",
    "\n",
    "test = rs[rs.route_id.isin(route_ids)]\n",
    "nums = test.groupby(['route_short_name', 'date']).num_trips.sum().to_frame().unstack(level=0).reset_index()\n",
    "nums = nums.assign(dayofweek = lambda x: x['date'].dt.strftime('%A'))\n",
    "nums =nums[~nums.dayofweek.isin(['Friday', 'Saturday'])].set_index('date').drop(columns='dayofweek', axis=0)\n",
    "nums.columns = nums.columns.droplevel(0)\n",
    "\n",
    "nums.fillna(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-08-15T21:00:35.782396Z",
     "start_time": "2018-08-15T21:00:35.282510Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x1ed2e15ae48>"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6kAAAI2CAYAAACymM6RAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzs3Xt8VfWd7//3d++dO5uEkBsEsoMkAUIwiKm2OCpIvXTGYmm99NRHLXVGPUfrcToznk5PR52ix5mfZ9pT608tlnFE20EQO6KgOGixHlqLDSpJCHfMBUggEEJC7jv7e/5IAgFyITfWSvJ6Ph48dvZ3r72+n+1y7b3e6/vdexlrrQAAAAAAcAOP0wUAAAAAANCJkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA2f0wVIUkJCgk1PT3e6DAAAAABwxLZt245ZaxOdrsMNXBFS09PTlZ+f73QZAAAAAOAIY0yp0zW4BdN9AQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAa/QZUo0xLxpjjhpjis5pf9AYs9sYs8MY81SX9h8aY/Z1PHbjcBQNAAAAABidfBewzEuS/n9JL3c2GGMWSrpF0qXW2mZjTFJHe7akb0qaLWmypPeMMVnW2rahLhwAAAAAMPr0OZJqrf1QUvU5zf9N0j9ba5s7ljna0X6LpFettc3W2s8l7ZN0xRDWCwAAAAAYxS5kJLU7WZKuNsb8L0lNkv7OWvsnSamS/thluYMdbYAr1B5v1Du/KFRkTJiSAuOVlO5XUmC8xk2IkDHG6fLQwVqr1uY2NTcE1dzQqqb69tvm+qCaGlrb2+vbb5s6bpsbWhUe5WvfroH27RqfGiOvl6/eu0lba+jMNjy9Hbts44bzt3GozWrilHGnt2timl+RMWFOvxR0YUNWzY1ntmHXffX0PtplX22qD6qlMSj/xEglpXfus36NT4jivdhFrLUKtoT6fh/ucr+pISgjKSYuQjGx4YqOi1BMbITGxUUoJi5c0bERiomLUFi41+mXN6a1tYXU3PV9t8v+2d1+21TfqrZgSPGTz7wXJwX8ivKHO/1SMEoNNKT6JE2Q9EVJX5C0xhhziaTuPllsdyswxtwr6V5JSktLG2AZwIULhazef2mnTh5tlJKkzzaVKRRq/98zyt8eWhO7vPHGxEY4XPHIdzqQdD2o6Xq/pw/I+uDpbdMdj8coIsaniOgwRUT7FOUPV1xytBrrWrT/k6Mq3nJYkuT1eZQwdZyS0vxK7NiuEybFyOPhIHgwQiGrlsaew0dzQ/uBanM3Bz3BllCv6w6P8iki2qfImPZtGxMXKcmqqvyUDnxadXq52MQoJQXat2tyul8JU/0KjxzoRxqk808OdRs0ezix0NwY7OHTvp0vzKOIaJ8iOrarf2KkwiO9qjnaqILN5QoF258cEe07vV07D4Q5iTh4bcFQN++z5+6nZ78vdz4Waut5wxqPad+uHfts5LgwxSZFy1qrhpMtOlJap4btxxRsPX+/D4/yKSY2vCPMtgfXmLjw039Hx7b/7fVxorEnZ04OdXMC4ZwTC+e+L7c29/5NvLBIryKjw05/1sZPipHxGB0/dEolhcdO7+/++MiOfdbffsIpza+IaE4iYvAG+ol+UNJvrLVW0sfGmJCkhI72qV2WmyLpcHcrsNa+IOkFScrLy+vlow0YGp9tKtPhvTW67q5ZmjV/koKtbTp+sF5HS2s7/tWpbMdx2Y7/G2PiIk6f3e8MsFHjxt4Zw1DIqqXLwU37B2D3H3rnHsh2d2BympEiOgJJZ9j0x0eedSDb9QMysksoDYvw9njQaq1V7bFGHS2pO71dd/2xUoW/OyRJ8oV7lJjmV1LamZH02MQomTEWXDsDyfnh4+wD2bOCSJcRsN74wj1nbbPxCVFnb9NonyJiOv8OO/13eLSv1xMITfWtqiqt09Gy9u1aceCk9uZ3fNvESBOSo8+aITFxyrgxOVrT9eTQWaPT3YyAnXsg29vJIeMxZ+2HnSeHOk8odN52buuubb6wnrdDWzCk6sP1p/fXo6W1+uw/uz+JmNxxOxZPInaeHOorfJy1T3c8FuwjkISffi9u32YxcTGKiAlTZJf357Pehzv22bDInt+LO1nbXnd9TYvqTza3/6tpPnO/plmH9p5Qw8mWbgNx5LiwLkH2TIg9HXDjIhTlDx+xJx/PPjl0zuh01+3bcWLhrG3cx8khb5infRvGnPmMTZwyrss+ev42jYjxKSLKJ08vs5BaGoOqKqtr31/LanW0pFb7uzmJ2DlLgpOIGAhjbd/50BiTLmm9tTan4/5/lTTZWvuoMSZL0vuS0iRlS/p3tX8PdXJHe2ZfP5yUl5dn8/PzB/EygN5Vlddp7T/nK/3SBN10b06PH6qtzW2qKq9rPxDuOGCqOdJw+nH/xPYzhglT/PKFj46zu6E2e/50yy4jZC1NvR/c+CK8Zx3IRHY5QD19kNPNB2J4VO+BZCjZkNWJIw2qOn0QXKdj5XWnQ3R4pPf0KHr0+NFxIsJadRzQnj+N64JGq72m24OY09u560HNOScWvGEXb99oqG3R0dLaMwdMJbVqqG2R1B6q4ifFKCndrwnJMTKjY5dVsCXU67T3vkarT++f3YWPbsJm5/3eTg4N+WtsbdOxg6fOei8+UVF/3knExDS/wiJGx4mIUFvHqFg32/SCR6vPPWHQ9X34rH32zHLhUd5eA8nFYkNWTfWtHcG1pT3Idgbak2fuN9a26NxDV2Ok6PHhHSOwZ0KsW/7faGkMngma3ey3FzJaHdnlc7T9vfb8fTTinPdl30U8SXfuScSjJbU6daK540VIE1JilBTwK35SjDxe508o+MK9yrnGfd9INMZss9bmOV2HG/QZUo0xqyQtUPtI6RFJj0l6RdKLkuZKalH7d1J/27H8jyTdLSko6a+tte/0VQQhFcMp2Nqm1/4pX031rfrmI1f0ezS0+fQZw9rTB0y1x5qGqVpneLymhwOZruGj8++zD2hH6lSsUFtI1RUNZ23XY4dOnZ52OCp0Ga0++yDm3JHqs0etL3YgGUrWWtXXdA2utTpaUqem+lanSxtSvnDPeQetXffhbk8sxFzck0NDraUpeF5w7XoScTQ496sM550w6Dq75Jz364t5cshJobaQGmpbTwfYhpNnh9jOgOu2fT48ynfWZ+d5n7U9nCAaqe/F0jknEUtqdaS0To0dJxGdFjU+XHc/9WdOl3EeQuoZFzSSOtwIqRhOW9bs1fbfluurD+YqbfbEIVlna3Nbr6NQI4nHY+QL94zYD8Gh1NYW6nMUaqQwUvvBzQgNJEPJWtvnjICRxOfzjJlA0pdgS5vaehmFGkmM0YgOJG7T1hpSMOiO9/OwCO+IPTk0lKy1am1q621CwEVj1H7iwG0IqWe4b+sAQ6i8uFrbf1uuOQumDFlAleSaKUQYWl6vR94oDv5HG2OMIlx4MILB84V7OZBBt7xhnMxxG2OMK4Mh3Im9F6NWU32r3l9ZrAkp0frS16c7XQ4AAACAC0BIxahkrdXv/n23Gutadf3ds8fkL3wCAAAAIxEhFaPSno+PaN+2o7pi8TQlpvmdLgcAAADABSKkYtSpPd6oD1ft1qSMWF12Q8DpcgAAAAD0AyEVo0ooZPX+SztlJX15aTa/pgcAAACMMIRUjCqfbSrT4b01uvr2LI1PiHK6HAAAAAD9REjFqFFVXqetbx7Q9MsSNfNLKU6XAwAAAGAACKkYFYItbXrv34oVOS5MC+6cycXQAQAAgBGKkIpR4aM39qv6cL0W3TVLkePCnC4HAAAAwAARUjHilRdXq+C3BzVnwRSlzZ7odDkAAAAABoGQihGtqb5V768s1oSUaH3p69OdLgcAAADAIBFSMWJZa/XBr3ersa5V1989W2HhXqdLAgAAADBIhFSMWHu2Vmr/J0d1xeJpSkzzO10OAAAAgCFASMWIVHu8UR++ukeTMmJ12Q0Bp8sBAAAAMEQIqRhxQiGr91/aKSvpy0uz5fFwuRkAAABgtCCkYsT5bFOZDu+t0TV3ZGl8QpTT5QAAAAAYQoRUjChV5XXa+uYBTb8sUTO+mOJ0OQAAAACGGCEVI0awpU2bXixW5LgwLbhzpoxhmi8AAAAw2hBSMWJ89MZ+naio16K7ZilyXJjT5QAAAAAYBoRUjAjlxdUq+O1BzVk4RWmzJzpdDgAAAIBhQkiF6zXVt+r9lcWakBKt+UumO10OAAAAgGFESIWrWWv1wa93q7GuVdffPVu+cK/TJQEAAAAYRoRUuNqerZXa/8lRXbF4mhLT/E6XAwAAAGCYEVLhWrXHG/Xhq3s0KSNWl90QcLocAAAAABcBIRWuFApZvf/STllJX16aLY+Hy80AAAAAYwEhFa702aYyHd5bo2vuyNL4hCinywEAAABwkRBS4TpVZXXa+uYBTb8sUTO+mOJ0OQAAAAAuIkIqXCXY0qZN/1asyHFhWnDnTBnDNF8AAABgLCGkwlU+emO/TlTUa9F3ZilyXJjT5QAAAAC4yAipcI3y4moV/Pag5iycorTsiU6XAwAAAMABhFS4QlN9q95fWawJKdGav2S60+UAAAAAcAghFY6z1uqDX+9WY12rrr97tnzhXqdLAgAAAOAQQioct2drpfZ/clRXLJ6mxDS/0+UAAAAAcBAhFY6qPdaoD1/do0kZsbrshoDT5QAAAABwGCEVjgmFrN57qVhW0peXZsvj4XIzAAAAwFhHSIVjPttUpop9J3XNHVkanxDldDkAAAAAXICQCkdUldVp65sHNH1eomZ8McXpcgAAAAC4BCEVF12wpU2b/q1YkePCtOBbM2UM03wBAAAAtCOk4qL76I39OlFRr0XfmaXIcWFOlwMAAADARQipuKjKi6tV8NuDmrNwitKyJzpdDgAAAACXIaTiommqb9X7K4s1ISVa85dMd7ocAAAAAC5ESMVFYa3VB7/ercZTrbr+7tnyhXudLgkAAACACxFScVHs2Vqp/Z8c1RVfnabENL/T5QAAAABwKUIqhl3tsUb97tU9mpQRq8tuCDhdDgAAAAAXI6RiWIVCVu+9VCxJ+vLSbHk8XG4GAAAAQM8IqRhWn20qU8W+k7rmm1kanxDldDkAAAAAXI6QimFTVVanrW8e0PR5iZpxZYrT5QAAAAAYAQipGBbBljZtenGHIseFacG3ZsoYpvkCAAAA6BshFcPio//YrxOVDVr0nVmKHBfmdDkAAAAARghCKoZceXG1CjYf1KULpygte6LT5QAAAAAYQQipGFJN9a16f2WxJkyK0ZeWTHe6HAAAAAAjDCEVQ8Zaqw9+vUuNp1p1/Xez5Qv3Ol0SAAAAgBGGkIohs3trpfZ/UqUrvjpNiWl+p8sBAAAAMAIRUjEkao816sNX92hSRqwuuyHgdDkAAAAARihCKgYtFLJ676ViSdKXl2bL4+FyMwAAAAAGhpCKQftsU5kq9p3UNd/M0viEKKfLAQAAADCCEVIxKFVlddr65gFNn5eoGVemOF0OAAAAgBGOkIoBC7a0adOLOxQ1LkwLvjVTxjDNFwAAAMDgEFIxYB/9x36dqGzQou9kK3JcmNPlAAAAABgF+gypxpgXjTFHjTFF3Tz2d8YYa4xJ6LhvjDE/N8bsM8YUGGPmDUfRcF5Z8XEVbD6oSxdO0dTseKfLAQAAADBKXMhI6kuSbjq30RgzVdL1ksq6NH9FUmbHv3slPT/4EuE2Tada9f7KnZowKUZfWjLd6XIAAAAAjCJ9hlRr7YeSqrt56P9I+h+SbJe2WyS9bNv9UVKcMWbSkFQKV7DW6oN/36WmU626/rvZ8oV7nS4JAAAAwCgyoO+kGmMWSzpkrd1+zkOpksq73D/Y0YZRYvfWSu3/pEpXLr5EiWl+p8sBAAAAMMr4+vsEY0y0pB9JuqG7h7tps920yRhzr9qnBCstLa2/ZcABtcca9eGrezQpI1Zzr2ebAQAAABh6AxlJnS5pmqTtxpgSSVMkfWKMSVH7yOnULstOkXS4u5VYa1+w1uZZa/MSExMHUAYuplDI6r2XiiVJX16aLY+Hy80AAAAAGHr9Hkm11hZKSuq83xFU86y1x4wxb0r6njHmVUlXSjppra0YqmIvpsoDJ3Xk81qny3CN44dOqWLfSS1aOkvjE6KcLgcAAADAKNVnSDXGrJK0QFKCMeagpMestf/aw+JvS/pzSfskNUj67hDVedGV7jiu/A0lTpfhKjO/mKIZV6Y4XQYAAACAUcxY2+1XRi+qvLw8m5+f73QZZwm2tKktGHK6DFeJiA5zugQAAABgVDLGbLPW5jldhxv0e7rvWOEL93J5FQAAAAC4yAZ0CRoAAAAAAIYDIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAa/QZUo0xLxpjjhpjirq0/W9jzC5jTIEx5j+MMXFdHvuhMWafMWa3MebG4SocAAAAADD6XMhI6kuSbjqnbZOkHGvtpZL2SPqhJBljsiV9U9Lsjuc8Z4zxDlm1AAAAAIBRrc+Qaq39UFL1OW3/aa0Ndtz9o6QpHX/fIulVa22ztfZzSfskXTGE9QIAAAAARrGh+E7q3ZLe6fg7VVJ5l8cOdrQBAAAAANCnQYVUY8yPJAUl/bqzqZvFbA/PvdcYk2+Mya+qqhpMGQAAAACAUWLAIdUY8x1JN0u601rbGUQPSpraZbEpkg5393xr7QvW2jxrbV5iYuJAywAAAAAAjCIDCqnGmJsk/UDSYmttQ5eH3pT0TWNMhDFmmqRMSR8PvkwAAAAAwFjg62sBY8wqSQskJRhjDkp6TO2/5hshaZMxRpL+aK39r9baHcaYNZKK1T4N+AFrbdtwFQ8AAAAAGF3MmZm6zsnLy7P5+flOlwEAAAAAjjDGbLPW5jldhxsMxa/7AgAAAAAwJAipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAADX8DldAAAAAADgfNu2bUvy+XwrJOVo9AwwhiQVBYPBv7r88suPdrcAIRUAAAAAXMjn861ISUmZlZiYeMLj8Vin6xkKoVDIVFVVZVdWVq6QtLi7ZUZLGgcAAACA0SYnMTGxdrQEVEnyeDw2MTHxpNpHh7tf5iLWAwAAAAC4cJ7RFFA7dbymHrMoIRUAAAAA0K3bbrstPT4+PjczM3N2Z9uLL744ISMjY7bH47n8ww8/jO5sr6ys9F555ZVZ0dHRl911111pA+2TkAoAAAAA6Nbdd9997M0339zbtW3u3LmNr7/++r68vLxTXdujo6PtsmXLDv/jP/7jwcH0SUgFAAAAAHTrK1/5yqnExMRg17Z58+Y15ebmNp+77Pjx40M33njjqcjIyNBg+iSkAgAAAABcg0vQAAAAAIDLPbx2+9Q9lXXRfS954bJS/A3/+9bc8qFc51BgJBUAAAAA4BqMpAIAAACAy7lxxHO4EFIBAAAAAN366le/Ou2Pf/yj/8SJE77k5ORL//7v//7wxIkTgw8//HDaiRMnfEuWLMmcNWtWw5YtW/ZKUmpq6pxTp055W1tbzbvvvhv39ttv77n88sub+tMnIRUAAAAA0K233nrr8+7a77rrrpru2g8dOlQ42D75TioAAAAAwDUIqQAAAAAA1yCkAgAAAABcg5AKAAAAAHANQioAAAAAwDUIqQAAAAAA1yCkAgAAAAB6FQwGNWvWrOyFCxdmSNKTTz6ZmJaWlmOMubyiouL0pU2ff/75+KysrOysrKzsyy67bOZHH30U1d++CKkAAAAAgF498cQTyRkZGY2d96+99tpTmzZt2jN58uSWrstlZGQ0//73v9+9Z8+e4h/+8IeH77vvvkB/+yKkAgAAAAB6tH///rB333039p577jnW2XbVVVc1zpgxo+XcZa+//vr6xMTENklauHBhfWVlZXh/+yOkAgAAAAB69MADD0x96qmnDno8/YuPzzzzTMLChQtP9rc/X9+LAAAAAAAc9cYDU3W0OHpI15mU3aCvPVve2yKrVq2KTUhICF599dUN69ev91/oqt966y3/r371q4Q//OEPu/pbFiEVAAAAANCtLVu2jNu0aVNcampqbHNzs6e+vt5zyy23TFu3bt3nPT1n69atUffff39gw4YNe1NSUtr62ychFQAAAADcro8Rz+Hy7LPPHnr22WcPSdL69ev9P/nJT5J7C6h79+4Nv+2226a/+OKLn1966aXNA+mT76QCAAAAAPrliSeeSEpOTr70yJEj4bm5udl33HFHQJL+4R/+YVJNTY3vwQcfDMycOTM7JydnVn/Xbay1Q19xP+Xl5dn8/HynywAAAAAARxhjtllr87q2bd++vSQ3N/dYT88ZybZv356Qm5ub3t1jjKQCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAACgV8FgULNmzcpeuHBhhiSFQiE9+OCDqenp6TmXXHLJ7CeeeCKps33p0qVT09LScrKysrK3bNkS3d++fENdPAAAAABgdHniiSeSMzIyGk+dOuWVpGeeeWbiwYMHw/bv31/k9Xp16NAhnyS99tprsQcOHIgsKSkp2rx5c8z999+fVlBQsKs/fTGSCgAAAADo0f79+8Pefffd2HvuuedYZ9uKFSuSHn/88Qqv1ytJSk1NDUrSunXr4u68887jHo9HixYtqq+trfWVlpaG9ac/QioAAAAAoEcPPPDA1Keeeuqgx3MmPpaXl0e88sorE3JycmZdc801mYWFhRGSVFFREZaent7SudykSZNa+htSme4LAAAAAC73yO8fmbrvxL5+f7+zNxkTMhoev+rx8t6WWbVqVWxCQkLw6quvbli/fr2/s72lpcVERkbaoqKinStXroxbunRp+rZt23Zba89bhzGmX3URUgEAAAAA3dqyZcu4TZs2xaWmpsY2Nzd76uvrPbfccsu05OTklm9961snJOnb3/52zfe+9710SZo8eXJrSUlJeOfzKyoqwtPS0lr70ychFQAAAABcrq8Rz+Hy7LPPHnr22WcPSdL69ev9P/nJT5LXrVv3+f3335/6zjvv+GfMmHH87bff9gcCgWZJWrx4cc1zzz2XdM8991Rv3rw5xu/3twUCAUIqAAAAAGD4LFu2rPLWW2+d9txzzyVHR0eHfvnLX5ZI0u23335yw4YNsYFAICcqKiq0YsWKkv6um5AKAAAAAOjTzTffXHfzzTfXSVJCQkLbBx98sO/cZTwej1555ZWywfTDr/sCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAACgV8FgULNmzcpeuHBhhiSFQiE9+OCDqenp6TmXXHLJ7CeeeCJJkj799NPIuXPnzgwPD5/36KOPJg+krz6vk2qMeVHSzZKOWmtzOtriJa2WlC6pRNLt1toTxhgj6WlJfy6pQdJSa+0nAykMAAAAAOAOTzzxRHJGRkbjqVOnvJL0zDPPTDx48GDY/v37i7xerw4dOuSTpKSkpODTTz9dtnbt2gkD7etCRlJfknTTOW1/L+l9a22mpPc77kvSVyRldvy7V9LzAy0MAAAAAOC8/fv3h7377rux99xzz7HOthUrViQ9/vjjFV6vV5KUmpoa7Ly99tprG8LCwuxA++szpFprP5RUfU7zLZJWdvy9UtLXurS/bNv9UVKcMWbSQIsDAAAAADjrgQcemPrUU08d9HjOxMfy8vKIV155ZUJOTs6sa665JrOwsDBiqPrrc7pvD5KttRWSZK2tMMYkdbSnSirvstzBjraKgZcIAAAAAGPb4f/5o6nNe/dGD+U6IzIzGyY/+b/Ke1tm1apVsQkJCcGrr766Yf369f7O9paWFhMZGWmLiop2rly5Mm7p0qXp27Zt2z0UdQ00pPbEdNPW7TCvMeZetU8JVlpa2hCXAQAAAAAYrC1btozbtGlTXGpqamxzc7Onvr7ec8stt0xLTk5u+da3vnVCkr797W/XfO9730sfqj4HGlKPGGMmdYyiTpJ0tKP9oKSpXZabIulwdyuw1r4g6QVJysvLG/B8ZQAAAAAY7foa8Rwuzz777KFnn332kCStX7/e/5Of/CR53bp1n99///2p77zzjn/GjBnH3377bX8gEGgeqj4HGlLflPQdSf/ccbuuS/v3jDGvSrpS0snOacEAAAAAgNFh2bJllbfeeuu05557Ljk6Ojr0y1/+skSSysrKfF/4whey6+vrvcYYu3z58uSdO3cWxcfHhy503Rdv5xNKAAAgAElEQVRyCZpVkhZISjDGHJT0mNrD6RpjzF9KKpN0W8fib6v98jP71H4Jmu9e+MsEAAAAALjVzTffXHfzzTfXSVJCQkLbBx98sO/cZdLS0oJHjhwpGEw/fYZUa+1/6eGhRd0sayU9MJiCAAAAAABj14VcJxUAAAAAgIuCkAoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAA6NGxY8e8N9100yXTpk2bfckll8x+7733Yo4cOeKdP39+ZiAQyJk/f35mVVWVV5LWr1/v9/v9c2fOnJk9c+bM7L/7u7+b1N/+CKkAAAAAgB7de++9U2+44Ybazz//fEdxcXHx3Llzmx577LFJCxYsqCstLS1asGBB3aOPPprSuXxeXt6pXbt2Fe/atav4X/7lXyr621+f10kFAAAAAIxN1dXVnq1bt/rXrl1bIkmRkZE2MjKybePGjXG/+93vdkvSfffdd/zaa6+dIenQUPTJSCoAAAAAoFu7du2KiI+PD952223ps2bNyr7jjjsCtbW1nuPHj/sCgUCrJAUCgdbq6urTA6CffvrpuBkzZmRfc801mfn5+ZH97ZORVAAAAABwufdf3jm1+tCp6KFcZ3zquIZFd80q722ZYDBodu7cGf3000+XXXfddfXf/e53pz7yyCMpPS0/f/78+tLS0oLY2NjQ6tWrY7/xjW9klJaWFvWnLkZSAQAAAADdSk9Pb0lOTm657rrr6iXpjjvuOLF9+/boiRMnBktLS8MkqbS0NCw+Pj4oSfHx8aHY2NhQx7Ing8Ggqaio6NfgKCOpAAAAAOByfY14Dpe0tLRgSkpKy/bt2yNyc3Ob//M//3P8jBkzmmbMmNG0fPnyiU8++WTl8uXLJ9500001klRWVuabMmVK0OPxaPPmzdGhUEjJycnB/vRJSAUAAAAA9OiZZ54pu/POOy9paWkxaWlpzatWrSppa2vTkiVLpgcCgYTJkye3vPHGG/sl6Ve/+tWEF198Mcnr9drIyMjQyy+/fMDj6d8EXmOtHZYX0h95eXk2Pz/f6TIAAAAAwBHGmG3W2ryubdu3by/Jzc095lRNw2n79u0Jubm56d09xndSAQAAAACuQUgFAAAAALgGIRUAAAAA4BqEVAAAAACAaxBSAQAAAACuQUgFAAAAALgGIRUAAAAA0K3t27dHzJw5M7vz37hx4y5btmxZ0pEjR7zz58/PDAQCOfPnz8+sqqryStKnn34aOXfu3Jnh4eHzHn300eSB9ElIBQAAAAB0Kzc3t3nXrl3Fu3btKi4qKiqOjIwMffOb36x57LHHJi1YsKCutLS0aMGCBXWPPvpoiiQlJSUFn3766bL77rvvyED7JKQCAAAAAPr05ptvjk9LS2vOyspq2bhxY9x99913XJLuu+++4++8884ESUpNTQ1ee+21DWFhYXag/RBSAQAAAAB9WrVqVfytt956XJKOHz/uCwQCrZIUCARaq6urfUPVz5CtCAAAAAAwPN59/mdTj5WXRg/lOhOmBhpu/G9/XX4hyzY1NZn33nsv9qc//enBoayhO4ykAgAAAAB6tXbt2tjs7OyGqVOnBiVp4sSJwdLS0jBJKi0tDYuPjw8OVV+MpAIAAACAy13oiOdwefXVV+Nvv/326tP13HhjzfLlyyc++eSTlcuXL59400031QxVX4RUAAAAAECP6urqPFu2bBm/cuXK0s62H//4xxVLliyZHggEEiZPntzyxhtv7JeksrIy3xe+8IXs+vp6rzHGLl++PHnnzp1F8fHxoQvtj5Dagx+/tUPFh2udLgMAAADAEMqePF6PfXW202WMKH6/P1RTU/NZ17aUlJS2jz76aM+5y6alpQWPHDlSMJj++E4qAAAAAMA1GEntAWdXAAAAAODiYyQVAAAAAOAahFQAAAAAgGsQUgEAAAAArkFIBQAAAAC4BiEVAAAAANCjxx9/PCkzM3N2RkbG7GXLliVJ0pEjR7zz58/PDAQCOfPnz8+sqqrySlIoFNLSpUunpqWl5WRlZWVv2bIlur/9EVIBAAAAAN3605/+FPnyyy8nfvLJJzt37ty5Y+PGjXGFhYURjz322KQFCxbUlZaWFi1YsKDu0UcfTZGk1157LfbAgQORJSUlRc8//3zp/fffn9bfPgmpAAAAAIBuFRYWRs2bN++U3+8PhYWF6aqrrqpbvXp13MaNG+Puu+++45J03333HX/nnXcmSNK6devi7rzzzuMej0eLFi2qr62t9ZWWlob1p09CKgAAAACgW3Pnzm3cunWrv7Ky0ltXV+fZtGlTbHl5efjx48d9gUCgVZICgUBrdXW1T5IqKirC0tPTWzqfP2nSpJb+hlTf0L4EAAAAAMBQq167Z2prZX2/v9/Zm7CUmIb4W7PKe1tm3rx5TQ899FDlddddlxUdHR3Kzs5u8Pl6jpHW2vPajDH9qouRVAAAAABAj77//e8fKy4u3pmfn787Pj6+LTMzs2nixInBzhHS0tLSsPj4+KAkTZ48ubWkpCS887kVFRXhaWlprf3pj5FUAAAAAHC5vkY8h9OhQ4d8qampwb1794Zv2LAh7uOPP971+eefRyxfvnzik08+Wbl8+fKJN910U40kLV68uOa5555Luueee6o3b94c4/f72zqnBV8oQioAAAAAoEeLFy+eXlNT4/P5fPZnP/tZWWJiYtuPf/zjiiVLlkwPBAIJkydPbnnjjTf2S9Ltt99+csOGDbGBQCAnKioqtGLFipL+9kdIBQAAAAD0aNu2bbvPbUtJSWn76KOP9pzb7vF49Morr5QNpj++kwoAAAAAcA1CKgAAAADANQipAAAAAADXIKQCAAAAAFyDkAoAAAAAcA1CKgAAAADANQipAAAAAIAePf7440mZmZmzMzIyZi9btixJkl588cUJGRkZsz0ez+UffvhhdNflf/jDH6akpaXlpKen57z++uvj+9sfIRUAAAAA0K0//elPkS+//HLiJ598snPnzp07Nm7cGFdYWBgxd+7cxtdff31fXl7eqa7Lb9u2LfI3v/lN/O7du3ds3Lhxz1//9V+nBYPBfvVJSAUAAAAAdKuwsDBq3rx5p/x+fygsLExXXXVV3erVq+PmzZvXlJub23zu8mvXro37+te/Xh0VFWVnzpzZEggEmj/44IOY/vRJSAUAAAAAdGvu3LmNW7du9VdWVnrr6uo8mzZtii0vLw/vaflDhw6FT506taXz/uTJk1t6W747vsEUDAAAAAAYfm+88cbUo0ePRve95IVLSkpq+NrXvlbe2zLz5s1reuihhyqvu+66rOjo6FB2dnaDz9dzjLTWntdmjDm/sReMpAIAAAAAevT973//WHFx8c78/Pzd8fHxbZmZmU09LTtlypSzRk4PHz4cPmXKlNb+9MdIKgAAAAC4XF8jnsPp0KFDvtTU1ODevXvDN2zYEPfxxx/v6mnZb3zjGzV33nnnJY8++uiR0tLSsJKSksgFCxbU96c/QioAAAAAoEeLFy+eXlNT4/P5fPZnP/tZWWJiYtvLL78c9/DDD6edOHHCt2TJksxZs2Y1bNmyZW9eXl7T1772teqsrKzZXq9XP/3pT0t7mx7cHdPdnOGLLS8vz+bn5ztdBgAAAAA4whizzVqb17Vt+/btJbm5ucecqmk4bd++PSE3Nze9u8cYSQUAAADQs7agtP+3Utt5VxsZmbwRUtYNTleBXgwqpBpjvi/pryRZSYWSvitpkqRXJcVL+kTSt621LT2uBAAAAIB7ffZr6a3/7nQVQycmSXp4r9NVoBcDDqnGmFRJ/11StrW20RizRtI3Jf25pP9jrX3VGPMLSX8p6fkhqRYAAADAxVWwWpqYId32ktOVDA0Pk0ndbrBbyCcpyhjTKilaUoWk6yR9q+PxlZL+UYRUAAAAYOSpKZNKfy8t/AcpZY7T1WCMGPB1Uq21hyT9i6QytYfTk5K2Saqx1gY7FjsoKXWwRQIAAABwQOHa9ts5tzpbB8aUAYdUY8wESbdImiZpsqQYSV/pZtFufz7YGHOvMSbfGJNfVVU10DIAAAAADAdr26f6Tv2iFD/N6Wowhgw4pEr6sqTPrbVV1tpWSb+RNF9SnDGmcxrxFEmHu3uytfYFa22etTYvMTFxEGUAAAAAGHKVhVLVLunS252uBA57/PHHkzIzM2dnZGTMXrZsWZIk/eEPf4jKzc2dOXPmzOycnJxZmzdvjpakUCikpUuXTk1LS8vJysrK3rJlS3R/+xtMSC2T9EVjTLQxxkhaJKlY0mZJnfMBviNp3SD6AAAAAOCEgtXtPzI0e4nTlcBBf/rTnyJffvnlxE8++WTnzp07d2zcuDGusLAw4uGHH57yox/96PCuXbuKH3nkkcM/+MEPpkrSa6+9FnvgwIHIkpKSoueff770/vvvT+tvn4P5TupWSWvVfpmZwo51vSDpB5L+xhizT9JESf860D4AAAAAOCDU1v591MwbpOh4p6uBgwoLC6PmzZt3yu/3h8LCwnTVVVfVrV69Os4Yo5MnT3olqaamxpucnNwiSevWrYu78847j3s8Hi1atKi+trbWV1paGtafPgf1677W2sckPXZO8wFJVwxmvQAAAAAcVPJ/pVOVTPWF5s6d27hs2bLUyspKb0xMjN20aVNsbm5u/c9//vPyv/iLv8h85JFHpoZCIW3ZsmWXJFVUVISlp6e3dD5/0qRJLaWlpWGBQKD1QvvkIkEAAAAAzlawRooYL2Xd5HQl6FC88wdT60/t6ff3O3sTMy6rIXvW/1fe2zLz5s1reuihhyqvu+66rOjo6FB2dnaDz+fTz3/+88R/+qd/Kl+6dGnNihUrJixdujT9D3/4wx5rz//d3PZvh164wXwnFQAAAMBo09IgFb8pZS+WwqKcrgYu8P3vf/9YcXHxzvz8/N3x8fFtmZmZTa+//vrEu+66q0aS7r777hMFBQUxkjR58uTWkpKS8M7nVlRUhKelpV3wKKrESCoAAACArva8I7XUSXOY6usmfY14DqdDhw75UlNTg3v37g3fsGFD3Mcff7zrhRdeSHr77bf9N998c91bb73lDwQCTZK0ePHimueeey7pnnvuqd68eXOM3+9v689UX4mQCgAAAKCrgjWSf7KU/mdOVwKXWLx48fSamhqfz+ezP/vZz8oSExPbnn/++dK/+Zu/mfq3f/u3JiIiIvSLX/yiVJJuv/32kxs2bIgNBAI5UVFRoRUrVpT0tz9CKgAAAIB29celfe9JX7xf8nidrgYusW3btt3ntt14442nduzYsfPcdo/Ho1deeaVsMP3xnVQAAAAA7Xb8RgoFpUvvcLoSjGGEVAAAAADtCtZISbOllBynK8EYRkgFAAAAIFUfkA5+LF16m9OVYIwjpAIAAACQCl5rv51DSHWRUCgU6t9FRkeAjtcU6ulxQioAAAAw1lkrFa6R0q+WYqc4XQ3OKKqqqoodTUE1FAqZqqqqWElFPS3Dr/sCAAAAY93hT6Tj+6SrHnK6EnQRDAb/qrKyckVlZWWORs8AY0hSUTAY/KueFiCkAgAAAGNdwRrJGyHNWux0Jeji8ssvPyppzG2U0ZLGAQAAAAxEW6tUuFaacZMUFed0NQAhFQAAABjTDnwgNRyT5tzudCWAJEIqAAAAMLYVrJEi46TM652uBJBESAUAAADGruZT0q710uwlki/C6WoASYRUAAAAYOzatUFqbZAuvcPpSoDTCKkAAADAWFWwWopLk6Ze6XQlwGmEVAAAAGAsqjsiHdgszblN8hAL4B783wgAAACMRTt+I9kQv+oL1yGkAgAAAGNRwWppUq6UNNPpSoCzEFIBAACAsaZqj3T4U34wCa5ESAUAAADGmsI1kvFIOd9wuhLgPIRUAAAAYCyxtn2q77RrJX+K09UA5yGkAgAAAGNJ+VappoypvnAtQioAAAAwlhSskXxR0qybna4E6BYhFQAAABgrgi3tl56Z+RdShN/paoBuEVIBAACAsWLfe1LjCab6wtUIqQAAAMBYUbBaip4oTV/odCVAjwipAAAAwFjQdFLa/U77ZWe8YU5XA/SIkAoAAACMBTvfktqameoL1yOkAgAAAGNBwWop/hIp9XKnKwF6RUgFAAAARruTh6TP/2/7KKoxTlcD9IqQCgAAAIx2RWslWWnObU5XAvSJkAoAAACMdgVrpNQ8aeJ0pysB+kRIBQAAAEazIzukI0X8YBJGDEIqAAAAMJoVrJGMV8r5utOVABeEkAoAAACMVqGQVPialPFlKSbB6WqAC0JIBQAAAEar0t9LtYekS293uhLgghFSAQAAgNGqYLUUPk6a8edOVwJcMEIqAAAAMBq1NknFb0qzviqFRztdDXDBCKkAAADAaLT3Xan5JFN9MeIQUgEAAIDRqGCNNC5Zmnat05UA/UJIBQAAAEabhmppz7tSzq2Sx+t0NUC/EFIBAACA0ab4DSnUylRfjEiEVAAAAGC0KXhNSpghTcp1uhKg3wipAAAAwGhyolQq+0P7KKoxTlcD9BshFQAAABhNCl9rv51zm7N1AANESAUAAABGC2ulgtVS2pekCQGnqwEGhJAKAAAAjBYV26Vje/jBJIxohFQAAABgtCh8TfKESdlfc7oSYMAIqQAAAMBoEGprD6lZN0rR8U5XAwwYIRUAAAAYDT7/nXTqCFN9MeIRUgEAAIDRoGCNFBErZd7odCXAoBBSAQAAgJGupUHa+ZaUvVgKi3S6GmBQCKkAAADASLf7banllHTpHU5XAgwaIRUAAAAY6QrWSONTpcBVTlcCDBohFQAAABjJ6o9J+96T5twmeTi8x8jH/8UAAADASFb0G8m28au+GDUIqQAAAMBIVrBaSs6Rkmc7XQkwJAYVUo0xccaYtcaYXcaYncaYLxlj4o0xm4wxeztuJwxVsQAAAAC6OL5fOpTPKCpGlcGOpD4taaO1dqakXEk7Jf29pPettZmS3u+4DwAAAGCoFb4myUg5tzpdCTBkBhxSjTHjJV0j6V8lyVrbYq2tkXSLpJUdi62U9LXBFgkAAADgHNa2T/WddrUUm+p0NcCQ8Q3iuZdIqpL0b8aYXEnbJD0kKdlaWyFJ1toKY0zS4Mu8+KpXrlT1r37tdBkAgFHEGxenyJzZisrJUWROjiKmT5fxDeajGMCYdmibVH1A+rO/6dfTrLUK1derreakQrUn1XbypNpO1rbf1p5UqLZWbTUn1Vbbpe1UfXsoHgW8EyZo2prVTpeBXgzmk9EnaZ6kB621W40xT6sfU3uNMfdKuleS0tLSBlHG8PClTFLUZXOdLgMAMIoEj1ap9q31qln1qiTJREYqctYsRebkKCpntiLnzFF4eroMl5AAcAFC+f+utuZItYXPVig///ywefrv9tvQyY5AWlcntbX1uF4TFiZPbKy8sbHyjh+vsKRkeS4ZJ3nMRXx1w8c7zu90CeiDsQM8I2KMSZH0R2ttesf9q9UeUjMkLegYRZ0k6QNr7Yze1pWXl2fz8/MHVAcAACOJDYXUUlqqpqIdaioqVGPRDjUVF8s2NkqSPDExiszObg+uc9pHXMOmTpUxo+PgEMDZbDB4dojsHL08Wau2kzXto5pnhc2a0+HTtrT0vGKPR16//6yw6Y2NlSd2fMf9WHk7/vaMHy9vbNzp+yYykvccBxhjtllr85yuww0GPJJqra00xpQbY2ZYa3dLWiSpuOPfdyT9c8ftuiGpFACAUcB4PIqYNk0R06Yp9qs3S5JsW5taDhxQY2GRmoqK1LijSCd+/WtVdxyAemJjFTV7tiJzcuRLSebgERghrLUKnarvNWyG6ut7XYcnJqY9WI5vD5sR0y6RNy5WnmC1vPvfkPequ+XNmt8RNmPljWsPpJ5x45iVgRFrwCOpkmSMmStphaRwSQckfVftP8a0RlKapDJJt1lrq3tbDyOpAACczba2qnnvXjUWFampaIcaiwrVvGevFAw6XRqAfjLh4WdGMcd3GdmMO2cUs3O0szNs+v0yYWHdr3TtX0r735f+do/kC7+4LwjDgpHUMwb1aw3W2s8kdfcfctFg1gsAwFhnwsLap/1mZ0sdlz8MNTcrVFfnbGEA+sUzbpw8kZFDu9LmOmnXBmnufyGgYlTiJwUBABghPBER8kREOF0GAKftXC8FG6VL73C6EmBYMFEdAAAAGEkK10hxadLUK52uBBgWhFQAAABgpKirlA580D6Kyo+oYZQipAIAAAAjRdHrkg1Jc253uhJg2BBSAQAAgJGiYLU0aa6UmOV0JcCwIaQCAAAAI0HVbqliOz+YhFGPkAoAAACMBAVrJOORcr7hdCXAsCKkAgAAAG4XCrX/qu8lCyV/stPVAMOKkAoAAAC4XflWqaaMqb4YEwipAAAAgNsVrJbCoqWZf+F0JcCwI6QCAAAAbhZskXb8R3tAjRjndDXAsCOkAgAAAG62b5PUVMNUX4wZhFQAAADAzQpWS9EJ7T+aBIwBhFQAAADArRprpN0bpTm3Sl6f09UAFwUhFQAAAHCrnW9Kbc3Spbc7XQlw0RBSAQAAALcqWCPFT5cmz3O6EuCiIaQCAAAAbnTyoFSypf0Hk4xxuhrgoiGkAgAAAG5UuFaSlS69zelKgIuKkAoAAAC4UcEaacoVUvwlTlcCXFSEVAAAAMBtKoukozv4wSSMSYRUAAAAwG0KVksenzT7605XAlx0hFQAAADATUJt7d9HzfiyFDPR6WqAi46QCgAAALhJ6e+lusNM9cWYRUgFAAAA3KRgtRTul7K+4nQlgCMIqQAAAIBbtDZKxW9K2Yul8GinqwEcQUgFAAAA3GLPRqm5VprDtVExdhFSAQAAALcoWCONS5GmXeN0JYBjCKkAAACAGzRUS3s3/b/27jy6qvO89/j3EZIYxDwIMKOZkQ42HmI7cYyB2I4TEsezHWOD3Nzk9qZpUrddTTPcNkm7elfuzWradLXpzVAGTzEekjjE47IN8RTbwLWxxGyMMZhJaEDzcM5z/9j7CBkzSUjaW5vfZy0WOudsSc/hYe+zn3e/77Nhzs2Q0yfqaEQioyJVRERERCQOyn4NmRY477aoIxGJlIpUEREREZE42LgKRs2GMXOijkQkUipSRURERESiVvEuvP9HOO8WMIs6GpFIqUgVEREREYna248Ef6urr4iKVBERERGRSLnD26tg0uUwdGLU0YhETkWqiIiIiEiU9r0J5dvgvFujjkQkFlSkioiIiIhEaeMq6JMPRV+IOhKRWFCRKiIiIiISlXRrsB51+jXQf1jU0YjEgopUEREREZGovLsG6g7q3qgi7ahIFRERERGJysaHod+Q4EqqiAAqUkVEREREotFcB5t/B0XXQ16/qKMRiQ0VqSIiIiIiUdjyBLTUaaqvyDFUpIqIiIiIRGHjQzB4PEz8eNSRiMSKilQRERERkZ5WewjeeR7OuwVydEou0p72CBERERGRnlb2GHhaU31FjkNFqoiIiIhIT9v4EIyZA4Wzo45EJHZUpIqIiIiI9KTyHbB3va6iipyAilQRERERkZ709irAIHVT1JGIxJKKVBERERGRnuIeTPU9dx4MPifqaERiSUWqiIiIiEhP2bMOKndpqq/ISahIFRERERHpKRsfgtx+MPvzUUciElsqUkVEREREekK6BUofhZmfhX6Do45GJLZUpIqIiIiI9IQdz0FDBZx3a9SRiMSailQRERERkZ6w8SHoPxymfirqSERiTUWqiIiIiEh3azwCW5+A1I2Qmx91NCKxpiJVRERERKS7bVkNrY3q6ityGlSkioiIiIh0t40PwbDJMP5jUUciEnsqUkVEREREutORfbBzLcy5FcyijkYk9nKjDkBEREREJNFKHwG813b1fe43t9KS2R11GF0mL2cin7p+VdRhyEnoSqqIiIiISHfauArOuRBGTo86EpFeQVdSRURERES6y8HNsH8jXPvDqCPpNF11lJ6mK6kiIiIiIt1l4yqwPsGtZ0TktJxxkWpmfczs/5nZ6vDxuWb2mpltN7OHzEw3ghIRERGRs08mA28/DFMXwMDCqKMR6TW64krqN4DN7R7/EPixu08HKoEvdcHvEBERERHpXXa/CtXv696oIh10RkWqmY0HFgG/CB8bsBB4JNxkBXD9mfwOEREREZFe6e1VkFcAsxZFHYlIr3KmjZP+BfgbYFD4eARQ5e6t4eM9wLgz/B3R+ONP4fWfRR2FiIiIiPRW1Xug+AbIL4g6EpFepdNFqpl9Djjo7uvNbH726eNs6if4/q8AXwGYOHFiZ8PoPoPGwriLoo5CRERERHqr8ZfAJ++JOgqRXudMrqReDlxnZp8F+gGDCa6sDjWz3PBq6njgg+N9s7v/DPgZwMUXX3zcQjZSxdcHf0RERERERKTHdHpNqrt/y93Hu/tk4HbgeXdfDLwA3BxuthT47RlHKSIiIiIiImeF7rhP6jeBvzSzHQRrVH/ZDb9DREREREREEuhMGycB4O5rgDXh1zuBS7ri54qIiIiIiMjZpTuupIqIiIiIiIh0iopUERERERERiQ0VqSIiIiIiIhIbKlJFREREREQkNlSkioiIiIiISGyoSBUREREREZHYUJEqIiIiIiIisaEiVc466XQ66hCkG2QyGTKZTNRhSDfQPptMymsyuTvuHnUY0g20z0pPyo06AJGeUF9fz5YtWygtLWXXrl0sXLiQT37yk1GHJWcok8mwa9cuSktL2bx5MyNHjmTx4sX069cv6tDkDFVUVFBaWkpZWRmVlZUsXryYSZMmRR2WnKGmpia2bt1KaWkp77zzDpdccgnXXHMNZhZ1aHIGMpkMe/bsobS0lE2bNmFmXHHFFVx44YXk5upUszerrq6mrKyMsrIyDhw4wM0338ysWbOiDkvOAhaH0a6LL77Y161bF3UYkjCNjY0fOhnKZDIMGzaMQYMGsXv3bq655ho+8YlPRB2mdFAmk2H37t2UlZWxadMm6urqyBzujC4AACAASURBVMvLY8qUKWzfvp1x48Zx55130rdv36hDlQ6qqqqirKyM0tJS9u3bB8D48eOpr6+ntraWO++8k4kTJ0YcpXRUc3Mz27dvp7S0lO3bt9Pa2sqgQYMYNWoUO3fu5PLLL+eqq65SodrLuDt79+5tK2COHDlCnz59mD59OvX19ezevZvBgwczb948LrjgAvr06RN1yHKaampq2LRpE6Wlpbz//vsAjBkzBnfn0KFD3H777cyYMSPiKJPJzNa7+8VRxxEHKlIlUZqbm9m6dStlZWVs376ddDrNkCFDKC4upri4mHPOOYdMJsOjjz7Kpk2buPbaa7nsssuiDltOwd3Zs2dP28lQTU0Nubm5zJgxg+LiYqZPn05+fj5lZWU88sgjTJw4kcWLF5Ofnx916HIKR44cacvrnj17ADjnnHPa9tmhQ4dSU1PDsmXLqK2tZcmSJYwfPz7iqOVUWlpa2LFjB2VlZWzdupWWlhYKCgooKioilUoxYcIEzIwnnniCN954gyuuuIKFCxeqUI05d2f//v1tsxyqqqrIyclh2rRpFBcXM3PmTPr164e788477/DCCy+wd+9ehg4dypVXXsl5552nYjWm6urq2LRpE2VlZezatQuAwsJCiouLSaVSjBgxgoaGBlauXMnBgwf54he/yLRp06INOoFUpB6lIlV6vZaWlrZR+m3bttHa2srAgQPbDqzjxo0jJ+fDy6/T6TQPP/wwW7ZsYdGiRXzsYx+LKHo5EXdn3759bSdD1dXV9OnTh2nTppFKpZgxY8Zxr5a+/fbbPPbYY0yePJk77riDvLy8CKKXk6mtrW0bpd+9ezcAo0ePJpVKUVxczPDhwz/yPdXV1Sxfvpz6+nqWLl3KOeec09Nhyym0trayc+dOSktL2bp1K01NTfTv35+ioiKKi4uZPHnyR47FmUyG1atXs2HDBubPn8/8+fOjCV5O6uDBg23H4sOHD2NmTJkyhVQqxaxZs+jfv/9xv8/d2b59Oy+88AL79u1j+PDhzJ8/n1Qq9ZH/C9Lz2i+Fevfdd3F3RowY0XYsLiwsPO73rFy5kvLycu644w6mTJkSQeTJpSL1KBWp0iu1trZ+aJS+ubmZAQMGtI3ST5w48ZQfgK2traxatYpt27bx+c9/nosuuqiHopcTcXcOHDjQNuWzsrKSnJycD50Mnc5607feeotf//rXTJ06ldtvv12FagzU19ezefPmtnXh7s7IkSNJpVKkUilGjhx5yp9RVVXF8uXLaWxsZOnSpYwdO7YHIpeTSafTvPvuu5SVlbF582YaGxvp168fs2bNIpVKce65557yylkmk+Hxxx/nzTffZOHChcybN6+HopeTKS8vbzsWHzp0CDNj8uTJFBcXM3v2bAoKCk77Z7k7W7ZsYc2aNRw4cICRI0cyf/58ioqKVKz2sBMthcoO7I8ePfqUMxrq6upYsWIFFRUVLF68mHPPPbeHok8+FalHqUiVXiOdTrNz5862k6HsKP3s2bPbRuk7Oo2otbWVX/3qV+zYsYMvfOELXHDBBd0UvZzMoUOH2kbpy8vLMTPOPffctpOhAQMGdPhnbtiwgccff5zp06dz2223qXlHBBoaGtiyZQtlZWXs3LmTTCbD8OHD206GCgsLOzy9s7KykmXLltHS0kJJSQmjR4/upujlRDKZDO+9915bw7L6+nry8/OZNWsWxcXFTJ06tcP7WyaT4Te/+Q0bN27k6quv5vLLL++m6OVkKisr247F+/fvB2DixImkUilmz57NoEGDzujnZzIZNm/ezAsvvEB5eTmFhYUsWLCAWbNmaap3NzqdpVAd/fevra1lxYoVVFVVceedd6qxXRdRkXqUitQTaGhooLa2llGjRkUdylktnU5/6GSooaGBvn37to3ST5ky5YzXt7S0tPDggw+yc+dObrjhBs4///wuil5O5vDhw22j9AcPHgRg0qRJbSdDAwcOPOPfsW7dOlavXs3MmTO55ZZbVKj2gGO7t6bTaYYOHdp2MjR27NgzPhmtqKhg2bJlpNNpSkpKjjslTbrWsd1ba2trycvLY8aMGaRSKaZNm3bGMxYymQyPPfYYpaWlfPrTn+bjH/94F0UvJ5Pt3lpaWsoHH3wAwLhx40ilUhQVFTFkyJAu/52ZTIbS0lLWrFlDRUUFY8eOZcGCBUyfPl3FahfpzFKojqqpqWH58uXU1NRw1113MWHChC6K/uylIvUoFakn8PTTT/Paa69x+eWXM2/ePE0X7EHH696an5/PzJkzKS4uZtq0aV1ebDQ3N/PAAw/w3nvvceONNzJnzpwu/fkSOF731gkTJlBcXExRURGDBw/u8t/5+uuv88QTTzB79mxuvvlmNe3oBs3NzWzbtq1tlD7bvbX9yVBXn3iWl5ezfPly3J277777tKYLS8ecqHtrtmHZjBkzurw5WTqdbmts95nPfIZLL720S3++BGpqatrymu3eOnbs2LbBpGHDhvVIHOl0mo0bN7J27VqqqqoYN24cCxYsYOrUqSpWO6ErlkJ11JEjR1i+fDl1dXXcddddamx3hlSkHqUi9QRqa2t55pln2LhxI0OHDuWzn/2s2m13o2z31uwoffvuralUiunTp3f7QEFzczP3338/u3fv5pZbbqGoqKhbf9/Z4kTdW7Oj9EOHDu32GP74xz/y1FNPUVxczI033qhCtQtku7dmR+mz3VuzJ7kTJkzo9rVmhw4dYvny5eTk5FBSUsKIESO69fedDU7WvTWVSjFz5sxuv71T+8Z2n/vc57j4Yp2vdYUTdW/NNsmJcv9Jp9O8+eabrF27liNHjjBx4kQWLFigtY6nIbsUqrS0lC1btnTJUqiOqq6uZtmyZTQ2NrJkyRI1tjsDKlKPUpF6Crt27WL16tWUl5cza9YsPvOZz3TL1JezUWe7t3anpqYm7rvvPvbu3cutt96qG1Z30vG6t44ZM6atgDle99bu9vLLL/Pss88yZ84cbrjhBjXr6IT23Vu3bNnSNko/e/ZsUqkUkyZN6vF/14MHD7J8+XJyc3O5++67e+wKUNJku7eWlpZSUVFx2t1bu0v7xnbXXXcdF154YY/+/qQ4XvfWkSNHts1yiNuSptbWVjZs2MCLL75ITU0NkydPZsGCBVrveIyeWArVUVVVVSxbtozm5maWLl3KmDFjevT3J4WK1KNUpJ6G1tZWXn31VdauXYuZMX/+fC677DJdjemEE3VvnTp1KsXFxafdvbU7NTY2cu+997Jv3z7dsLoD6urq2Lx5c9sovbszatSotlH6OEzHfPHFF3nuueeYO3cu1113nQrV03Ci7q3ZUfrT6d7a3fbv38+KFSvIz8/n7rvv7pGr80nQld1bu0P7xnbXX389c+fOjTSe3qKxsbGtYVn77q3ZY/HpdG+NWktLC+vXr+fFF1+krq6OqVOnsmDBgrN6KmkUS6E6qqKiguXLl9Pa2srSpUvV2K4TVKQepSK1AyorK3nyySfZtm0bhYWFLFq0SKN7p+lE3Vuzo/Sd6d7anRoaGrj33ns5cOCAblh9EtnuraWlpezcuRN3Z/jw4R86GYqbNWvWsGbNGi688EI+97nPqVA9jpN1b82O0kd9MnSsDz74gJUrV9K/f39KSko04+UEKioq2qbfd0f31q7WvrHdjTfeyHnnnRd1SLHU1NTEtm3bKC0tZceOHR/q3ppKpbqkYVkUmpubeeONN3jppZdoaGhg+vTpLFiw4KyZThqHpVAddfjwYZYtW4a7U1JSErur9XGnIvUoFamdsGXLFp588kmqq6uZO3cuV199deQjznF0vO6t2VH6oqKi2P+b6YbVx3ey7q2pVIoxY8bE/mTo+eef5w9/+AMXX3wxixYtin28PaEnurd2t71797Jy5UoKCgooKSnplkZcvdHxureOHz++bfp93P+d2je2u+mmm0ilUlGHFAvH696abVhWXFzM+PHjE3Nsa2pq4rXXXuOVV16hsbGRWbNmMX/+/EROKXV3Pvjgg7bBpOxSqOnTp7c1LOvppVAdVV5ezrJlyzAzSkpKYjGTqrdQkXqUitQT+OHrP2RLxZYTvm5pY8jeIQzeN5hMnwxVE6qoLayFZHwedFqfxj4UVBQw4PAA+tYFB9HGQY3Uj6infng96fx0xBF2TE5LDqM3jya3MZeDMw/SNKQp6pAiYWmjf1V/BhweQP/K/uR4Dq35rdSPqKduRB3NBc296/++w9D3hzLkgyEcGXOEykmVvSv+ruKQX5dPweFgn81tziVjGRqGNVA/op6GoQ14n+g/Izoivyaf0ZtHk85Ps79oP5n8TNQhRaJPcx8GHB7AgMMD6FcbLKFoKmhq22fTfXvXsdjSRuGWQvrW9OXQjEM0DG+IOqRoZKB/VX8KDhcEx+JMDum8NPXDg7w2DWpK9LHMWo3B+wczeN9gctI51A2vo3p8NS0DWqIO7cw45NXntR2L85rycHMahgTH4vph9Xhu7zoW59XnMXrTaDzHOVB0gNZ+rVGH9CGzhs/im5d8M+owPkJF6lHxmq/Vi3gfp2piFXUj6xj+7nBGvDuCgkMFVJxbQUtBLz9YdlD2ZKjgcAF9a4PCtKmgiYqJFdSPqO91J0PtZfIyHJh9gNGbRlO4tZCDsw7SNPgsKVSPczLUmtdK7eha6kfU0zSwF58MGVRNqMLcGLxvMBhUTjxLCtUTnQwNbaBqYlVwMtTLCtP2mgc1c3DWQQq3FDJ682gOFB0gk3d2FKo5LTkMqBhAQXkBfWv6YhjNA5qpnFBJ/Yj62J0kdoT38SCvmwsZtX0Uh6afRYVqBvpV9wv22coB5KRzSOemqRtZR/2IehoHN54dxy7Ac53q8dXUjKlh8L7BDNo3iAEVA6gbWUf1uGpa+/eu/+O5DbkUlBdQcLiAvMY8HKdxSCNHxh2hfng9mdzee+xqGdDCgaLg/Gn0ptGxLFQl3nQltQu4O2+99RbPPPMMDQ0NXHrppSxYsCD20zHORBy7t3ans+WG1Sfq3lpUVERxcXEk3Vu7k7vz1FNPtd0T+aqrrkrM9LhjtW9YVlFRQU5ODlOmTGlrWNbT3Vu7265du7jvvvsYPnw4S5cujf3ygs46UffW7LrwpK0HO1sa26XTaXbt2tXWsCzbvTXbsCyK7q1xVFdXxyuvvMJrr71GOp3m/PPPZ968ebE+B+ntS6E6at++faxYsYJ+/fpRUlKixnanoCupR6lI7UINDQ0899xzrFu3jkGDBnHttddSVFSUmJPe3tC9tTsl9YbVJ+vemkqleuQea1Fyd37/+9+zbt065s2bx8KFC6MOqcuUl5e3NSxr37012yQnbg3LutrOnTt54IEHGDlyJEuWLEnM+01C99Yz0dDQwMqVKzl48GCiGttlu7dmG5a1796aSqWYOnVq7BqWxUVtbS0vvfQSb7zxBu7O3LlzmTdvXmwKosrKyrY1pvv27QNgwoQJbfcLj1vDsq72wQcfsGLFCgYMGKDGdqegIvUoFandYM+ePaxevZr9+/czYsQI8vPzow7pjGUyGQ4ePIi7M2LEiLYmOYWFhVGH1qOqq6vbCtUob3zelaqrq3tF99bulMlkWL16NRs2bKCwsDARRXlzczOHDx8GYNKkSW2j9AMHDow4sp61Y8cOHnzwQQYMGJCI9+7uHDp0KDHdWzurvr6eFStWUF5enpjPoZqaGmpra8nNzW27rUgcu7fG2ZEjR3jxxRfZsGED7k5hYWHk+0VLSwvl5eUAjBs3rm3G2dlWqO3Zs4eVK1eSm5sbi/fev39/lixZEnUYH6Ei9SgVqd0knU6zfv16duzYEXUoXaawsJDi4uJe0b21O1VVVfHcc8/R1JSMtan9+/dn1qxZvaJ7a3fKZDKsXbu2bZS7t2t/z8u4d2/tbu+88w6vv/46cfi86wrDhw9PXPfWzqirq+PZZ5+lvr4+6lC6RPaqaW/o3hp3VVVVvPrqq1RWVkYdCmbGhAkTKC4uZtiwYVGHE6k9e/bw8ssvk05H36ukX79+3HjjjVGH8REqUo9SkSoiIiIiIhIxFalHJacDioiIiIiIiPR6KlJFREREREQkNlSkioiIiIiISGyoSBUREREREZHYUJEqIiIiIiIisaEiVURERERERGJDRaqIiIiIiIjEhopUERERERERiQ0VqSIiIiIiIhIbKlJFREREREQkNlSkioiIiIiISGyoSBUREREREZHYUJEqIiIiIiIisaEiVURERERERGJDRaqIiIiIiIjEhopUERERERERiQ0VqSIiIiIiIhIbKlJFREREREQkNlSkioiIiIiISGyoSBUREREREZHYMHePOgbM7BDwXtRxJNxIoDzqIKTLKa/Jpdwmk/KaTMprcim3yRTXvE5y91FRBxEHsShSpfuZ2Tp3vzjqOKRrKa/Jpdwmk/KaTMprcim3yaS8xp+m+4qIiIiIiEhsqEgVERERERGR2FCRevb4WdQBSLdQXpNLuU0m5TWZlNfkUm6TSXmNOa1JFRERERERkdjQlVQRERERERGJDRWpIiIiIiIiEhsqUhPEzIZFHYOInB4zG2dm+eHXFnU80nXMbKGZFUQdh4icHjMb0O5rHY8TwswuNLO8qOOQzlGRmgBmNsTM1gPLo45Fuo6ZfdnM/sPMpkYdi3QdM7vNzEqBHwP3AriaAySCmS0Oj8ULgJao45GuYWZfMbNvhF+rgEkQM7vTzN4A/tnMvg46HieBmd1hZm8BnwYyUccjnZMbdQDSJRxoBOaY2RXu/mLUAUnnmVkOcAvwN8A+4FIz2+vujdFGJmfKzD4GfAP4iru/YmabzexCd98QdWzSOWHRkkuQ1+8An3H3P0YblXQFM+sH/BXwVWCAmf3W3XdFG5V0FTO7AfgK8NfAEGCxma1y9/3RRiadFe6z3wNuB+5w91favWYagOhddCW1F8pOEQy/zgH6AA8C/wr8MKq4pGu4ewZ4E/gY8FNgHjA70qCkq0wBXgwL1NFAKVAVcUzSSWaW54EWYBtwP/CemeWb2U1mdk7EIUonmFkfgHBgcJ27jwN+DvxjpIHJGcvmNnQlcL+7rwVqgCMqUHu3cJ89CKwAXjOz/mZ2jZkNUoHa+6hI7WXM7O+BB8xsqZkVhgXNMGCRu/8r0GJmf2Jml0cbqXSEmX3bzC4NvzZgp7tXAY8ABlyhNce9T/u8hnYDE83sYeANgtz+wsx+GG6vqYS9hJl9C/ilmd1tZoOAFwjy+ySwAbgBWGFm3wm31+dtL2Bm3wN+YmY3hU89G/79fYJZLQvC7ZTPXqZdbm8On3oR+Asz+y+CpReTzOwBM/vTcHvluBcws6+Z2Zx2Tz0ODACeAl4nuFq+3My+Em6vvPYSSlQvYmb3AJ8kuLp2FfCdcKS+HMhOaXiFYMT3m9oR48/MxprZowRTe++DYD2Mu7eEU1NagEeBi4ALj/leFTQxdby8Arj7q0AJsBn4rrvfDHwJWGJm4zTSG39mNsvMXgGKgYeBm4Cl7l5DcNL7JHCtu98J3AP8tZmNCAcUJcbCIuYSgsL0z8PP3KEA7l4H/AT4fnhsVj57kePk9i+B3wELCc6FP+vu1wA/A/7MzIYqx/FmZpPMbC3wXeCfs8+7+w5gHVAGfCr8nP0p8FUzG6K89h4qYnqJcIrKBcD33f054B+ABuC/E+Txy2b2PMH0lT8AG7Uj9grVwMPuPhSoCj84MbPcbMHi7s8AuwjWHC8ysz8Ln1dBE18nzCvBGvJBBB+guPu7BINLMyKKVTqmBljl7ne6+++Ax4BPhK+tB77n7nsA3L2UYDR/ZCSRymkLO4B+Evgrd/8N8PfAOcAXs9u4+78RLK+5wcwmmtmiSIKVDjlObv+OILd/AhwCJhN8xgJsJJgJMbjHA5WOqiBYYjEdyJhZSbvXfg38rbsfDB9vIsht/x6NUM6IitReIBy1TQMHCK66AOwguMJ2IcEOej/whLt/ArgRuN7MdGIUc+5eD/w+fHgPwdXxfHdvNbOcdlfDnwK+TXCVPP84P0pi5GR5dfdmgkL12+FamR8B4wjWp0rMuftegv0w6zVgkJn1c/dmd2+C4MTYzP6N4GT3vQhClRM4dhaKmeWEs1Y2c7QofYVgSv4cM2s/gPQjgmUYfyCYUigx0oHcvk4wQ2lk+NovzGwgwbTukQTrGiUmjpNXC2ev3Bv+/Z/A1+zo7WZass0mw4s83yFocHeoB8OWM6QiNYbM7Dprd9uRdlfM/i8w3swuCq+SvkdwgnSeu3/L3X8Ubl8JXODu5T0du5zYsXnNcvea8ID7ErCW4GCLu2fcPWNmo4D/QzA1aZq7/7hHA5eT6mheQ98iOMn9H+HjT7m7Pjxj5iS5rWv3cCHwfvvu22b2BYIT4TRwizpzx07b1ZRjpu6uJlgzPissbN4GjhBcdcPMLgL+J8FncZG7P9yzYctp6Ehuq4BJwN8SDBxmBxa1z8bPsXnNzjRrCJ/+LUHzuu+Hz3u47RKCqb8twJfCCz7SS5hmDMaHmV1FMI13GnBDeHKbHQnMWNDV9x7gQne/LXztJ8Bb7v7LcCphWtNA4+UkeTUIDqbh9N5WCzq+biaY+lkIZNx9i5mN1KBDvJxBXscAze6+zcz6t/uQlZjoYG7/Bfiju/8qLGJ2Elw9NdftSmIlnJ57D8HUzhfc/f7w+T7unjazCcCXgTx3/1b42mpgmbs/Gr6Ou78fyRuQE+pkbn8PrHT3h8ysL9DX3Y9E8w7keE6S1xzCFh7ttr2IYE3xfILP2f3ARKApXKcqvYyupEbMAgPN7HcEi7+/C/yRYHQvuzYxOxI4hKAD3Qgz+244wj+T8Kbx4VRCFagxcJp59fBkdxThPYvd/QDBGreDwHKC9U+oQI2HLsrrfwF54fMqUGOiE7ntG35rATDKzJYRFLZD3P09FajxYmbXENw/8V8JpnouNLNzwkHg7NWVGuBpoNjMvm5mIwj24ToIilMVqPFzBrntQ9A/AHdvUoEaL6fIayY8Fg/NTvF19/XAWwRrVVcAw929TAVq76UiNWLhOU8twb265nvQFOlp4Avh661mlhuubfoFkAH+gmDqw0PAy+6+MqLw5QQ6kNefEKxvm2LBGtS7CDo3f9PdL3H3ssjehHyE8ppcncjtOeEV8jsJGti95e6fVXEaW1cCT3vQ7GodwRW1D7KDwGb2U4IOoVsJBhsuANYQXCV/KpqQ5TQpt8l0qrz+O0EBOy58fA9wNfBtd/+Eu6sXQC+XG3UAZysz+zrBOpf17v6wu/8qfD4HqATeN7O+HjThSAEDgRIP1pvuJ2jE8oPwdYmJTuR1EHB3mFfMbDMw14N7pEpMKK/J1QW5/Tawwt0ronkHcjzt8rrB3VcRrDd8Jlw2swTYamY/J7hCs45g4PeeMK/lZvbfCE6KtTYxZpTbZOpEXguAv8geiwm6956vY3Fy6EpqDwunlN0D3Eawk/3AzErC6WOEI0TvAouyBai7v+nud7t7pQVdygifV4EaE12Q1+y00HUqZOJDeU2urjoWu/uPdVIUH8fJ6/fN7EsEJ7ZzgVHAze5+JfAbgm74Ve5eckxe0ypi4kW5TaYuyGv2c/Y5HYuTRUVqD3N3BxYA33X3RwgWhJ8PXNtum1eAPRZ0iGxr1nHM+gqJkS7Ia2vPRy2norwml47FyXSCvM4BbvNgbdpkYF+4+dsEt3bLTh9UXmNMuU2mLsirPmcTSkVqN8qe0LR7nP33XgdcARCuh9hGsJh/VrjdYGAL0Bxuk221nUEip7wmk/KaXMptMnUgr1uBuWY2DXgO+N/hdiUE69kqw22V15hQbpNJeZWOUJHavfq3f9BuZ9pBcPP3OeHjtQSdeweG2x0BxgOjeyhO6RjlNZmU1+RSbpOpI3ktAIYC/wHkmtkaoBi4y9XVNY6U22RSXuW0qUjtBmZ2mZk9Cvy7mV2TXQeRnTdPMM8+DVxtwW0NNhGMDF3c7sfc7u7LezJuOTnlNZmU1+RSbpOpk3mdAFwSrln7InCru9/m7vujeA9yfMptMimv0hkqUruYmc0nGPV5jGC6wp3AsPbz5sM59m8Q3Cj+b8NvbSK4WTHhNlrUHyPKazIpr8ml3CbTGeS1EdgZvl7v7gd7OHQ5BeU2mZRX6SwVqV3vPOANd78fuA/IA2r96H2d/tHMfgmsB34CXGJm6wluPvxMRDHLqSmvyaS8Jpdym0zKa3Ipt8mkvEqnWNgHQjrJzC4DKtx9W/h4LvA8wQ2GvwpsJhgJehp4H/ga8HfhqBFmNhDIdd2eIlaU12RSXpNLuU0m5TW5lNtkUl6lq+hKaieZ2VAz+z3wLHBruFPh7m8S3MJgEvBVd58PvAxcBdS5+x3uviPb0czda7UjxofymkzKa3Ipt8mkvCaXcptMyqt0NRWpnVdAMAr05+HXV2RfcPfXCW4+/F741PMEHcoqoe2+TmqbHU/KazIpr8ml3CaT8ppcym0yKa/SpVSkdoCZLTGzK81ssLvvBX4GrCJY3H2pmZ0TbtcXeIVgWgPAp4Dh4Xa6r1PMKK/JpLwml3KbTMprcim3yaS8SnfSmtRTMDMDxgAPABngHYIRom+4e3m4zeXArcA6d783fK4Y+Pvwe1uAr7n75p5/B3I8ymsyKa/Jpdwmk/KaXMptMimv0lNyT73J2cvM+rh72swGAXvd/U4L7un0Y4LRohsB3P1lM7sEmGlmQ4Emdy8zs6XAWHffGdmbkI9QXpNJeU0u5TaZlNfkUm6TSXmVnqTpvsdhZrlm9k/AP5nZlcBMgpsM48E9nb4OfDx8LevnwECCBeO7zGycuzdoR4wP5TWZlNfkUm6TSXlNLuU2mZRXiYKK1GOEO9h6YBiwA/gHgmkJC8JRITyYI/0D4HvtvnURwVz7N4E54dx8iQnlNZmU1+RSbpNJeU0u5TaZlFeJiqb7flQG+FG7OfQXAOcCfwf8FLjIgjbZvybYQSe7+y6Cxd9XufsfoglbTkF5vAXAiAAAAyhJREFUTSblNbmU22RSXpNLuU0m5VUioSupH7UeWGVmfcLHLwMT3X050MfM/jzsQjYeSIc7Iu7+W+2Isaa8JpPymlzKbTIpr8ml3CaT8iqRUJF6DHevd/cmd0+HT10NHAq/vhuYbWargQeBDdDW6UxiTHlNJuU1uZTbZFJek0u5TSblVaKi6b4nEI4YOTAaeDx8ugb4NpAC3s3Orw/n4ksvoLwmk/KaXMptMimvyaXcJpPyKj1NV1JPLAPkAeXAeeEo0f8EMu7+khaA91rKazIpr8ml3CaT8ppcym0yKa/So0yDHSdmZpcBr4R/lrn7LyMOSbqA8ppMymtyKbfJpLwml3KbTMqr9CQVqSdhZuOBu4B/dvemqOORrqG8JpPymlzKbTIpr8ml3CaT8io9SUWqiIiIiIiIxIbWpIqIiIiIiEhsqEgVERERERGR2FCRKiIiIiIiIrGhIlVERERERERiQ0WqiIiIiIiIxIaKVBEROWuY2ffM7K9P8vr1ZlbUkzGJiIjIh6lIFREROep6QEWqiIhIhHSfVBERSTQz+w6wBHgfOASsB6qBrwD5wA6CG9TPBVaHr1UDN4U/4t+BUUA98GV339KT8YuIiJxtVKSKiEhimdlFwHLgUiAX2AD8J7DM3Q+H2/wjcMDd/83MlgOr3f2R8LXngD919+1mdinwv9x9Yc+/ExERkbNHbtQBiIiIdKMrgF+7ez2AmT0ePp8Ki9OhwEDg6WO/0cwGAp8AHjaz7NN9uz1iERGRs5yKVBERSbrjTRlaDlzv7m+ZWQkw/zjb5ABV7j63+0ITERGRY6lxkoiIJNkfgBvMrL+ZDQI+Hz4/CNhnZnnA4nbb14Sv4e5HgHfN7BYAC5zfc6GLiIicnbQmVUREEq1d46T3gD3AJqAO+JvwubeBQe5eYmaXAz8HmoCbgQzwU2AskAf8yt1/0ONvQkRE5CyiIlVERERERERiQ9N9RUREREREJDZUpIqIiIiIiEhsqEgVERERERGR2FCRKiIiIiIiIrGhIlVERERERERiQ0WqiIiIiIiIxIaKVBEREREREYkNFakiIiIiIiISG/8fy1JnZ+uODaMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1ed2e18b828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "nums.plot(figsize=(15,10))\n",
    "plt.legend(loc='center left', bbox_to_anchor=(1.0, 0.5))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "לפי הלו\"ז הרשמי של משרד התחבורה, מתוך 7 הקווים שהצהירו שיתוגברו בעקבות השיבושים ברכבות - רק 2 תוגברו באמת (980 ו-412).\n",
    "\n",
    "מה עם כל השאר? לא צריך לעדכן ת'לוז? מאלתרים? מה התדירות החדשה? מתי הנסיעות?\n",
    "\n",
    "למי שמתעניין, מזהי הקווים: 10953, 10954, 11178, 11179, 11941, 11942, 15544, 15545, 21331, 21332, 23451, 23452, 711, 712, 713, 714, 715, "
   ]
  },
  {
   "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.3"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": true
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}