climate_analysis.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib import style\n",
    "style.use('fivethirtyeight')\n",
    "import matplotlib.pyplot as plt\n",
    "# https://matplotlib.org/3.1.1/gallery/style_sheets/fivethirtyeight.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import datetime as dt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Reflect Tables into SQLAlchemy ORM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Python SQL toolkit and Object Relational Mapper\n",
    "import sqlalchemy\n",
    "from sqlalchemy.ext.automap import automap_base\n",
    "from sqlalchemy.orm import Session\n",
    "from sqlalchemy import create_engine, func"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "engine = create_engine(\"sqlite:///hawaii.sqlite\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# reflect an existing database into a new model\n",
    "Base =  automap_base()\n",
    "# reflect the tables\n",
    "Base.prepare(engine,reflect=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['measurement', 'station']"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# We can view all of the classes that automap found\n",
    "Base.classes.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Save references to each table\n",
    "Measurement = Base.classes.measurement\n",
    "Station =  Base.classes.station"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create our session (link) from Python to the DB\n",
    "session = Session(engine)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Exploratory Climate Analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='date'>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAEjCAYAAADpH9ynAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABEPklEQVR4nO3dd3gU1foH8G86IYUNSUhCGoKLVCliIKB4FWmCiFwF1J+oVAEbCpemIlgCIoo0L4goIHopgtIMKDWBhKaABAgLgVBSN8lmk03d3fn9EbJsrzOzO7vv53l4HrI7c+bM2Zl555w5c46XTCZjQAghhPDA29kZIIQQ4jko6BBCCOENBR1CCCG8oaBDCCGENxR0CCGE8IaCDiGEEN5Q0CGEEMIbVoNO586dIRKJDP6NHDmSzc0QQggRKF82Ezt06BBUKpXm74KCAvzrX//C8OHD2dwMIYQQgWI16EREROj8vXHjRoSEhFDQIYQQAoDDZzoMw2Djxo0YNWoUmjZtytVmCCGECAhnQefQoUPIzc3Fyy+/zNUmCCGECAxnQWf9+vXo3r07HnzwQa42AQCQSCScpu9uqLxsR2VmOyoz23hSeXESdIqLi7F371688sorXCRPCCFEoDgJOps2bUJAQABGjBjBRfKEEEIEivWgwzAMNmzYgBEjRiAkJITt5AkhhAgY60EnLS0NOTk51LRGCCHEAKvv6QBA3759IZPJ2E6WCIxSqYRCoXB2NljRpEkTlJeXOzsbgmJNmQUFBcHXl/VLEHFx9IvruVGhxDW5Er1a+CPIj4ams4dSqURFRQVEIhG8vLycnR2HBQQEoEmTJs7OhqBYKjOGYSCTyRASEkKBx8PQr60lo7AWI/aVoFrFQNzMF2nDWqCJr/AvmnxTKBRuE3AIN7y8vCASiSCXy9GsWTNnZ4fwiG7ltbyZLkO1igEASMqV+OlqlZNzJFwUcIgldIx4Jgo6Wq7KlTp/H8qrcVJOCCHEPVHQIYQQwhsKOoQ4mUgkwm+//Wb18rm5uRCJRPj77785yU9KSgqSk5M5SZsQCjqEOFl2djYGDRpk9fJxcXHIzs5G586dATS8GycSiVBSUmLTdk0FrzfffBN79uyxKS1CrEW91wixU11dHfz9/R1OJyoqyqblfXx8bF7HFsHBwZylTQjVdAi5a8iQIZg2bRpmzpyJxMREJCYm4oMPPoBarQbQMB17SkoKpk6dioSEBEyYMAEAcOLECTz11FOIiYlB+/bt8e6770Iul2vSZRgGy5cvR/fu3dGiRQt06NAB8+fP13yv3bzWWPvYunUrBg0ahKioKDz88MM4ePCgZnntGkpubi6efvppAECbNm0gEokwefJkAMCff/6JwYMHIzExEa1atcKIESOQnZ2tSadLly4AgMcffxwikQhDhgwBYNi8plar8fnnn6Njx45o0aIFevfurVMTaszPb7/9huHDhyMmJgaPPvooDh06xMKvQtwN1XTMYBhn58C9iL6/w+v2ZK/F2rzO1q1b8cILL+CPP/5AVlYW3n77bYSHh+Odd94BAKxatQrTp0/H4cOHwTAMsrKyMGLECMyaNQvLly9HWVkZZs+ejTfeeAMbNmwAACxYsADfffcdPv30U/Tp0wdSqRTnz583m4958+bh008/RceOHfHtt9/ixRdfxF9//YWWLVvqLBcXF4cNGzZgzJgxyMzMRFhYmOalTIVCgddffx2dOnVCdXU1vvjiC4wePRonTpyAv78/Dh48iCeeeAK//PILOnXqZLLW9s0332D58uX48ssv0a1bN2zevBkvv/wyDh8+rDN1ySeffIIFCxZgyZIlWLhwIcaOHYt//vmHak5EBwUdQrRERUXh888/h5eXF9q2bYurV69i9erVmqDTu3dvvP3225rlJ02ahGeffRZvvvmm5rMlS5agb9++KC4uRmBgIFatWoWUlBTNhIatW7dGUlKS2XyMHTsWzz77LABg0aJFOHjwINatW4f3339fZzkfHx+EhYUBACIjIxEeHq757plnntFZduXKlYiPj8eZM2eQnJysWbZ58+Zmm+tWrFiBN954A88//zwAYO7cuTh+/DhWrFiBNWvWaJabMmUKBg8eDACYM2cOtm7din/++Yc6JRAdFHTMoHfXPE+PHj10XlpMSkrCp59+qmku69atm87y586dQ05ODnbs2KH5jLlbRb5+/Tp8fHxQW1uLxx57zKZ8PPzww5r/e3t746GHHsLly5dtSuP69ev49NNPcfr0aZSUlECtVkOtVuP27dtWpyGXy5Gfn49evXrpfJ6cnIz9+/frfNaxY0fN/6OjowE0zK1FiDYKOmZQ8xrRFxQUpPO3Wq3GmDFjMGXKFINlY2JikJWVxVfWDIwePRoxMTFYunQpYmJi4Ovri549e6Kuro6V9PVHFPDz8zP4jqGTiOihoEN4Y88zFr6dOXMGDMNoLpqnTp1CdHQ0QkNDjS7fpUsXXLp0Ca1btzb6/QMPPICAgAAcOXIEbdq0sTofp0+f1tSOGIbBX3/9ZdBc1qjxWYxKpdJ8VlpaiuzsbCxevBh9+/YFAJw9exZKpdLsevpCQ0MRExODzMxMndpaRkYGHnjgAav3h5BGFHQI0VJQUIBZs2Zh/PjxuHjxIpYtW6Z5nmPM22+/jf79+2PatGl49dVXERISgitXriA1NRVLly5FSEgIXn/9dcyfPx/+/v7o06cPSktLcfbsWYwbN85kuuvWrcP999+PDh06YO3atbh16xbGjh1rdNn4+Hh4eXlh3759GDx4MJo0aQKRSITw8HBs2LABcXFxyMvLw4cffqgzonNkZCQCAwNx4MABJCQkICAgwOjgm2+++SZSUlLQpk0bdO3aFZs3b0ZGRgYOHz5sdbkawzAMVGqHkiACRF2mzaBnOp7n+eefh1qtRr9+/fDWW2/h5ZdfxqRJk0wu36lTJ+zduxc3b97E0KFD8cgjj2DBggWIjIzULDNv3jy88847WLx4MZKSkjBmzBjk5eWZzce8efOwcuVKPPLIIzhw4AB+/PFHxMYarym2bNkSs2fPxieffAKxWIwZM2bA29sb69atQ1ZWFpKTkzFjxgzMnTsXAQEBmvV8fX2xaNEibNy4Ee3atcOLL75oNP3XX38db775JubNm4fk5GTs2bMHGzZs0Om5ZqsqpRpZZUpsu16FD07RXEWexEsmkwm60VUikUAsFrOSln6X3qEJTfBjv3ATSwsTm+VlSnl5uSCHqx8yZAg6dOiAxYsX63xeU1PD23w6ubm56NKlCw4dOmTQaUFILJWZpLwelfUMTt0swX/OqXF6RAvc38zP5PLujo/z0lVQTYcQwrvKet173d25NKK7p6CgQwghhDfUkYCQu1xhkMvExETIZDJnZ4MQzlBNhxBCCG9YDzoFBQV4/fXX0aZNG0RFRaFnz55IT09nezOEEDdA7456Hlab12QyGQYOHIhevXphy5YtCA8PR25urk73UeIZtF+wJMQohkF5Hb2o42lYDTrLli1DdHQ0Vq9erfmsVatWbG6CCEBQUBBkMhlEIhEFHmIcwyCnsAxrrjEA6BjxJKwGnT179qBfv3547bXXkJaWhujoaIwZMwYTJkygi48H8fX1RUhIiM6cMkIml8tNDoNDjLNUZn9crcSaawyk9XRd8DSsvhzaODz6lClTMHz4cPzzzz+YOXMm5s2bh4kTJxpdRyKRsLV5hz2c3lTn73+FK7G4PTuDIxJC7tE/195oVYdX4pQmliZCYuklV1ZrOmq1Gt26dcO8efMANAyGmJOTg7Vr15oMOo6+hcvqm7zpuiMSBAcFQyymEQk8HZWZ7SyWmd65FhEeAbE4hONcuS5POsZY7b0WFRVlMPJs27ZtbZq/gxBCiPtiNej06tULV69e1fns6tWriI+PZ3MzhBBCBIrVoDNlyhScOnUKX3zxBXJycvDrr79izZo1GD9+PJubIYQQIlCsBp3u3btj06ZN2LFjB5KTk/Hxxx9jzpw5FHQIIYQA4GDstYEDB2LgwIFsJ0sIIcQN0NhrhBBCeENBxwwaFooQQthFQYcQQghvKOiYQQN0EMIPalXwHBR0CCGE8IaCjhl090UIP6hVwXNQ0CGEEMIbCjpm0N0XIYSwi4KOGdS8Rgg/6FzzHBR0CCGE8IaCDiHE6agp23NQ0DGDTgRCCGEXBR0zqJ2ZEELYRUGHEOJ0dIPnOSjoEEII4Q0FHUKI09HzU89BQYcQQghvKOgQQgjhDQUdM6jKTwg/qCOB56CgYwadCIQQwi4KOoQQp6NWBc/BatBJSUmBSCTS+de2bVs2N0GIy0svqMWIfVJMTS9DaY3K2dkhxKX4sp2gWCzG7t27NX/7+PiwvQlCXFaNksHoP0pQqWxonA3w9sKXvUXOzZQAUFO252A96Pj6+iIqKortZAkRhF251ZqAAwDrshUUdAjRwnrQuXHjBtq3bw8/Pz/06NEDH374IVq1amVyeYlE4vA22UijQVOdvxSVlZBISllK23WwV16ew9oyyy3wARBg17ruxvx+655rUqkUEkkBtxlyce5ynIjFYrPfsxp0evTogVWrVkEsFkMqlWLx4sUYMGAAMjMz0bx5c7syaIlEInE4DY30Ozp/BgUHQywOZydtF8FqeXkIW8qsBaMArsp0PvPE8rZYZnrnWmREBMTiEI5z5bo86bxkNej0799f5+8ePXqga9eu+Omnn/DGG2+wuSlCXBL1wiLEPE67TAcHB6Ndu3bIycnhcjOEEIGjjgSeg9OgU1NTA4lEQh0LCCGEAGC5ee3999/HoEGDEBcXp3mmU1VVhRdeeIHNzRDisryofY0Qs1gNOnl5eRg/fjxKSkoQERGBHj164I8//kBCQgKbmyGEuBmK1Z6D1aCzbt06NpMjhBDiZmjsNTPo4SYh/KBzzXNQ0CGEEMIbCjpmUDszIYSwi4IOISyiGxX7ULl5Dgo6ZlA7MyGEsIuCDiHE6egGz3NQ0DGDqvyEEMIuCjpm0N0XsRXdqBBiHgUdQojTUbD2HBR0CCGE8IaCjhl090UIP6gp23NQ0DGDTgRiKy8aZpoQsyjoEEII4Q0FHUKI01H90HNQ0CGERXTxtA81ZXsOCjqEEEJ4Q0GHEEIIbyjoEEII4Q0FHUJYRD2m7UPF5jko6BBCnI46EngOToPOkiVLIBKJMGPGDC43QwghRCA4CzqnTp3C+vXr0bFjR642QQghRGA4CTrl5eWYMGECli9fDpFIxMUmCCHE4xRXqzDmYAke+a0Im69VOTs7duEk6Lzzzjt45pln8Nhjj3GRPCHEzVBHAut8eb4CO3NrcKG0HlPSyiCtUTk7SzbzZTvB9evXIycnB6tXr7ZqeYlE4vA22UijQVOdvxSVlZBISllK23WwV16ew9oyKyzyARBg17ruxvx+655rxVIpJJICbjPk4qw5Tr65eK/cVAzwxbGbGJeg5DJbNhOLxWa/ZzXoSCQSLFiwAL///jv8/f2tWsdSBq3ZpqNpaKTf0fkzKDgYYnE4O2m7CFbLy0PYUmZ/e1cBV8p0PvPE8rZYZnrnWkREBMTiEI5z5bqsPsb0yi20eTjE4lCOcsUNVoPOyZMnUVJSguTkZM1nKpUKx48fx7p165CXl4eAgAAzKRBCCHFnrAadIUOGoFu3bjqfTZ06FW3atMG7775rde2HEEKIe2I16IhEIoPeak2bNkVYWBg6dOjA5qZ4wdAba8RG9EDcPlRunoNGJCCEOB3d33kO1nuv6duzZw/Xm+AMjaNFCCHsopqOGdS8RmxFNyqEmEdBhxBCCG8o6BBCCOENBR0zqKmEEELYRUHHDHqmQ2xF9ymEmEdBhxBCCG8o6BBCnI5qiJ6Dgg4hxOmoJdtzUNDh0I7rVRixT4p5p8pRq6LTyhPQHTvhkxCfO3M+IoGnulGhxGuHG4a4P5hXi6imPpjSMdjJuSKEEOeimg5HFp2t0Pl7zslyJ+WEEOKuhPhaBwUdMxyZClZWq2YxJ4S4NwFeO12CEJvXKOiYcaq4HpLyemdngxC3J8BrJ7ETBR0L5p2WOzsLhBBiFDWvuaG9N2ucnQUiIEK8CBDCJwo6hBBCeENBhxDidFRB9BwUdAghTkcdCTwHBR1CWORF9+yEmEVBhxBCCG8o6BBCCOENq0Hn22+/Re/evREfH4/4+Hj0798f+/btY3MThLg06jJNiHmsBp2WLVti/vz5OHLkCA4dOoS+ffvipZdewoULF9jcDCGEEIFiNegMGTIE/fv3R+vWrXH//ffjgw8+QHBwME6dOsXmZggRlKxSGkqJkEacPdNRqVT45ZdfoFAokJSUxNVmCHF5/XYXocSBwWMJMUWIA36yPp9OVlYWBgwYgJqaGgQFBeHHH39Ex44dTS4vkUgc3iYbaTRoylr6CoU/9IuXvXw6xlXyISTWllmB1AdAgM5nNSpg/tFbePM+z6rxmC8z3XNNKpVCIingNkMuzrpjTLfcSktLIZEUcpMhO4nFYrPfsx50xGIx0tLSUF5ejp07d2Ly5MnYvXs3OnToYFcGLZFIJA6noZF+x+jH9qQflFsClOqO28ZaPh3Aanl5CFvK7KJfNXC51ODzw+VNsEzciuWcuS6LZaZ3rkVEREAsDuE4V67L6mNMr9zCw5tDLA7lKFfcYD3o+Pv7o3Xr1gCAbt264a+//sKqVauwYsUKtjdFiGDcrKTmNcI+al4zQq1Wo66ujuvNEBYp6tXIr7p3kWzq642WQT5OzBFxVWqGwf+uVqG0Vo0xbYMQ6k+v/hHzWA06H330EQYMGIDY2FhUVlZi27ZtSE9Px5YtW9jcDOHYsYI6jPyzRPN3/9gAbB0Q4cQcCYenvaaz4IwcS/+pBABsv16Ng0+3cHKOPIsQ3wtjNegUFhZi4sSJKCoqQmhoKDp27Iht27ahX79+bG6GEOIiGgMOAPwlrcfFsnp0CPNzYo48i8c3r33zzTdsJkecRIDHMXER0hq1s7NAXBw1wBKLhFiFJ8QTCPHcpKBDDDBU1yGEcISCDrFIgDdThBAXRUGHGBDiw0lXIcTmDr7RkECejfP3dIgb8KArab9dRSi++zDcC8AfQyNtWt9zSso+DMNg6O9SZ2eDOBEFHUK03FKoUFR9rweWmmp9rDojrcclmdLZ2SBORM1rxABdZ++hmgu7iqupac3TUdAhFnnShZeeZxHCLQo6xIAnX3j1d92DHmcRwgu3faZzq1KJyvp7b5zEBfnQYIR28uTrrrP2nWEYZJcrIfL3RnRTGmzVk9UoGeRUKBHnJoPuum3QmXmiHHtv3pvP5scnmmNoYqATcyQcHlzRcbiWx1aQmpouw09Xq9DEB1j3r+Z4KoGOXU8kr1Nj8N5iZJW5T9DxmFt/T76QOsqTmphcoXntank9frpaBaBh1tHJaWX8Z4K4hPXZCmSVNfT2u60w7IQhxKZwt63peNB10mMUVauQUdgwN5MXgMhAbyRHBZhfyUHOOI5OFetOa11eJ7wrC8M0zLUD3CtDL0+6e2HJ99kKZ2eBdW4bdPQJ8Y7AWVy1qC6U1uOVQ/emgn68ZQB2DGQ36LjCuHNKN3g5KEPmjaQf8jR/P9EyANsH0pxMbBNiHHfb5jV7fosu4TQPiDGuclzzcSnWvzmxdd/ZuAioBBJzjL1z01h+BuXoKgeRmxHizbRgazo1SgYLz8pRVe6LuJoKhPp749UHgkwuL8Dfxmlc9UB2NCDYwxlNQipX/QH0HMqrtXpZijmkkWCDTkW9+u6shf7ATTnCA3SDDt1ZscdVi5KL39gVLvdKmgeNWEmI1znBNq9V67VByOsdP1MFcoPJOVctBl6a1/T+dsY5rTSyo3eM9FxyNlt+D77KsbBKhRH7pOiwOR9fna/gaavOI8RrlnCDjt6ZWa8GsmX3ev3oH+RC/HFchaveTblittjIk8pIR4JVWZUspMwfZ51u31ysxMG8WuRVqTH/jBw5chpc1NW4TdABgA9PyzX/d9ULpRDtyq2xvBAPeLlxcIGbE2M1HaH1aHPW+04NTe7af7t3bUeI1zlWg86XX36Jxx9/HPHx8WjTpg1GjRqFixcvsrkJDWOn4L5bNfj0Lzme3SfFbzdq9JYX1knragqqnN+8w0fTlyu8HGqspiN0Arw2CoIQW3BYDTrp6ekYN24c9u3bh507d8LX1xfDhw9HWRn7b1SbKuzF5yps6lWjk6YD+XF3qbdco7ajjY+A4CoXS6Hc0daqGKy5WInt+a7x+oEQL8rujtXea9u3b9f5e/Xq1UhISEBmZiYGDx7M5qZsDhB08FnPWFmV1Di/SxUftVVXPU5cMeYYK6vXDpdCoWQA6I0TJpSoKTBCLFZOn+lUVlZCrVZDJBJxuRnCg3oXbPJx1+Y1Yw7esa/2zjeFsQdShGjh9D2dWbNmoXPnzkhKSjK5jEQisSvtmxXeAJpYvXx+QQEkKvPPJWprm8BYHLYnjwqFP/SL1959ZZulfOQV+wDQHV6mtLQEEkkhh7my7E6Jbr4UiirWy1StDoR2OMu5dg1Nfaz/7fJKDMuukbVplJT4AvDX+Sy7XOkyx0+jwiLT+6qvSlEJiaTUZPlIpVJIJAUs5Kqpzl9yuRwSiZSFdLln7Petrzd+TWpUWlLq9PNSn1gsNvs9Z0Fnzpw5yMzMRGpqKnx8TA/JbSmDppQX1wHniq1ePio6GuLWTc0u43+xCFDUG3xuTx6DckuAUt3nIPbuK5skEonFfFzwrQKydZ/DNW8eDrE4lMusWSQJqAYu3Rt7LSg4CGJxAqvb8MrMA7RqdW3atEH+jWtW/3ZX9fKozdo0whVy4KZhrytXOH60RXlVAVese14bFBQMsTjcZPlERERALA5xPFPpd3T+DA0NhVgc5ni6HDN1XvqdK2gYatyE5uHNnX5e2oqToDN79mxs374du3btQqtWrbjYhMu2vbsDoZQtJ81rDg61I5Ci452LtFISF8D6M52ZM2di27Zt2LlzJ9q2bct28hq2PlQWyoWUmOaM39AZz3RMjfd2o8K1XnS0aUQCJ0Uddz/thXhdYzXoTJ8+HT/99BPWrl0LkUiEwsJCFBYWorJSWG9TE0Ou8kBdGzcdCew7i3PkStaCAmPiSvLZ33KjnwuB0F5uJdxhtXlt7dq1AIBnnnlG5/OZM2di9uzZbG7K5ghvzeKmTnZPY6wUXKFouM5Cz+2FBs3nXlaEtv9kyrDmEveTbclqnd9tXZstQX//bef0vnPBeyVWueLNoCWsBh2ZTMZmcrxzhQsrMY3r7szXjdRULG1CqWZYDzhCmWFTCKeLEPLoCCFeswQ59trtSiUG7rWtG6RVNR37suN2jJWDK14H2c6SPftYxeN7KS74ExBiM0EGHa5Oc2p2Ns0V7qg8pSOBJ3KBw8uoyno1Zp2Q4YU/S3CswPVe0BXi8SnISdzsKWdrLliu1WLuPK4QYKzBek3Hjm0IpKiInT4/W4H/XmxoPj2cVwvJC9EI9nOde3WhnKvaXKf0bGBPm7c1v41aiL8gTywV+ZniOiw+K8eJQu7uBrn+dazpNKCPDhl2GCt5hmHw6/VqLD4rR66Tuosvu3Cv5221isH/rlY5JR/uRJA1Ha5Q85p9skrrMXBPMZQMkHK2AgeGRqJbhL/lFR3EdtOCsfRcqvnCpTLDfW/PH7KrMC1DBgD470UFLoyMRqCvc8ugst61LhIudkhYRZg1HTvWseYEoaDTwNYu07NPlmsmHlMzwPS7Fwqu8XG+UfMae/berLZp+Wlax1FJrRqpt2xbH3D/30eINW1hBh2Orjb0TMc+mXpNamek9ZhzUsb6drg+wbi6mSENXk8zPU6bNaWYV0VnqDvwmOY1614OdWwb50vqkFWmRL9Y60bedVVsdJmuZXGi0RolAzUYg9ECWG9es/IzbRRyrCevc6y0guxoWnOFbvVccrX8WEOQQYercnakI8GRvFr8e78USgZoEegNcTNBFi1r2PqNdt6oxuS0MtSpGfSO4j+YWzqpqUmWHVyd02z/PK5WsXW1/FhDkFdGu5pBrFjGkQvIlLQyzXONomo1iqrr7E/Mydg4kNm6iEw4WqqpNR3J123Gs6e3mTn23DXyec7T+GXEHdAzHS2OnNJ3qlhsT3IHLP1G5prpXOE9HT7jgCtMGe7pXK05y9XyYw1BBh17WPVyqAveSKoZBmsvVeI/mTJcLDOcYM5VCfBcMMrSO2F8HjOudniymR9r0nKXY4pNQmxeE2TQcaR5TalmMO14GRI25eHf+6Uo0xq51xWDzn8vKjA9sxxrLinQf3cxKuq5v9tloxh46c7MR1XHAnqhmD9CvKsnhoQZdBw4+NILavF9dhXkdQwO3KnF+ux7IwSXuNjQ8QAw52S55v8KJYPvLxuOaFyjZKCoV6NKqUa1koHKBaInHxcIV7gG8XnEuML+cuVInuuNa0a4Icig44gPTulOhPXRmYa/s2XCaLq6ozB8yDH+SClif8xHy435iNmYh903axzahi3vnijVDOqMXHnd+QKp7Uwx+x1GTJWdO9/pH8mvRXE1+89F2a6IutpPIMRjQpBBh4ty3nhFWGMqbbiiwIDdxZieITMYXp/P43CyiRf+eKnpuMAJt+ICf7PiuntL3nIey9JervYTCPGY8Jwu0xZ+HG8XuIBZS1Jej7eOyQAAJ43caTt6Mbb2OK5XM9iaY3xoEra7M7sqH1eIfE7C9vWugIMeoB7887gsYQYdB44kU2suE8BdVqOFf1eY/Z6v80xp5oGGvXmYf7oc1SoGXrB8wXCF64mQblZcXS0HzyKFWBPQZukcWHyuAlM7BkMUIJxGK0EGHXuYO/biNubxlg82VFroweYK10F77wu+u6yA3MqRfF1hP/kMOgK/flpUI4BX3VzhmNO3/ooCb3cOcXY2rCac8KiF7R++kscph9lgKbd8Na+ZzYMTt80nXoOO0ArHRrUqbnbwVFEdxh0uxSdn5Jxtw5nmnZZbXsiFeE5Nx/2ONZO4uA7amqbdQceW38kFbjv5fKbj7odwtdbN3w0jk7bZU9IKpRrDUqWovhtsvLyAud1D7c0iYQHrNZ1jx45h9OjRaN++PUQiETZt2sT2JuziLiesVW9uO1rTMbIRo3PsmMmNvXkQ2u/kQ81rrNHuhTmfpbv3327UaAIO0PAMxBEucJ8jeKwHHYVCgQ4dOmDhwoUIDAxkO3kAzumRwjAMtlyrwudn5ci7+65MSY0KS85VYMMVwxc2nYmLnmMpf1fg7WNlOs+TzNVK+Ghec4ULQKsQ/hoL3H3uHu2azo4btk/Yxge+fwF3/MlZP2MGDBiAAQMGAACmTJnCdvIA7LvYvJshw9H8Wpwvte8l0FUXFZh7d3SAdZcV+GdkNAbvleJKOf9zt1vaf0cvxqaO8/VXqtBO5IfJHYPNLgc4UNOx4SRzhaCT1MIf67L5uelwxwuQtmoLz1bVDDDvVDl25lYjOSoAi3s1Q5CfIB9LezSP+sV+deDuaa7WcDQF1WrMOVHulIBjDS5rgrO1ysFs0LEzfXNNdq6Ii9ya+v1cbZAmtoOgykKCxwpr8fWFSlyvUOGnq1XYZuIdMeLanN6RQCKR2LyOQgkATVnPiykNedTd3rdGxkCzLh3rSOuAEacDoX/5LpfJoKjxBuBjct28vDuQVJu+RFnKR2GhDwDTE6Y1rl9p5ncoKyuDRFJkdjvGqNWG+2xKRUUFJJISm7dhctsqw2037qupMiuwUFbG0rJEKvUF4G/weW1tnV3ni76TMm/sLPRF2yA1XopV2v1cqsiGfW+Ul5dvcp06pUpr/wyPqy3XdIPM28dl6O2t/7qD5euCbWWom55UKoVEUmDD+tYzlq/6+iawpm7AxnHBFrFYbPZ7pwcdSxk0prJeDWTmc5Ab48RiMZB+h510rBRapUL1ScODu5lIBFmFEigzPUBiXGwsxLFNjH4nkUgs5iOKUQBXZSa/b1y/vM7079C8eXOIxXb0Ejp+x+rqw54iX3zTPx6h/uxU2L1P5gF6XWrFYrHZMmsBBSCRWZW+tb9/eFUFkGv4IN3X3x9icbxVaZiSX6XCm8cKoGaAfcVAYnQkxrcPtiutFhaOE2NatowBLpUa/c7bx+deGVl5vhmUqRXr2XTN0UsvMiICYjH778SYOsb8zhVY9QKTPddRZ3F60LGHK7Tlc82RJjJbVy2tUaFG1TBMPwPLbeuNbG1eqVUxWHdZAS+vhjw28fHCKw8E6SxjbPBQcz75S47Pe4lsW8kEeZ3t7UV8Nq+x0Zz15bkKnSk8pmeW2x102OYCg6MTHggy6HgCRwKrreuOPVKGw1pDyw9v5XivQ2N5qFYyOs+EmvnrBp0iO0YZXnNJwUrQmZ4hs2s9S0MSsYmN5135LjzDLZ9zEynVDK7JlYgK9BHUEDLugPWgU1lZiZycHACAWq3G7du3cf78eYSFhSE+3rGmgUZCHcSvqFqFFoGmn8VoM/umu4UCsLV89Be39uS3tfeafrr6iyxx8B0Ke10rV2KtHc/oAOC2kakmTJHVqrH6UiUCfbwwsX0wmvja9kNpF9/tSiW+z1YgLsgXrzzQFN5OOCnYDhF8hZx6NYNhqVJkFNYhsok3dgyMQKfmfprva1UM5p4sx6E8x6YIIcaxHuL//vtv9O3bF3379kV1dTVSUlLQt29ffPbZZ6xtQ6gjGC/7x/pBRfms6egHOGubOcy9N2IsD/pL618nM4vYn5vGGkfz+ZlA7IUDJUj5uwIfnpbjrePGp4Qwp7HlUalm0H9PMZacr8S0DJnDLzy6Cr4qOjtvVCOjsOFYK65R4/1T5Qbfr72swDW582uFQr3BNof1ms6jjz4KmUzGdrJuYUVWJT5JambVso6MpG1r2DGo6VixTmZhLbJlpruMWxN0vPWWcvf3UBovdEBDT6w1fW1bv7F8fr9Vg/yqe79Syt8VmNlV+EO7tA7lp7VffzqOw3qzlk44avsNAbGeIBsz3TD4GzC1j5uvVWlGRDC5roPNa9Zc/Jecq8Dbx2U25UG/BqW/jLNiTuptYTSjNIYZLuadcQV8BR1zzkrN17b/ktbj56tVkLng1PZCIcyg4wFRx9Q+ltcx+MfCqAoON69ZsY6lZYzWdPSiiv52Lb0cyJV9t4QRdGKsfB5ojrESHr5Pih+yFbwMs2NuC42bd2Znhzkny81+v+NGNSanleGxnUWo42DE6qJqFSYeLcXz+6WcTIXuCpx/a2EHD4g5jj3TsXlIaP1mLvMnU45cCYmF0RiMZUE/UNlTw/IEpsqhX5xtL2Ja63BeLQ7n1aJLuB+6RRi+lGrMJolCM3utLcw9L2z8ypk1uetGRrc2JrdShf9dq8KYtkGWF7bB7BPl+OV6Q/Pf+dISBNrY2UQIBFnT8QR8vqdj68W/+y+FuFlpqYnPMBeWajr0nkYDU70H2bj8mJtPxtJdfqMbFUpMTZfZtX3zQefu9AN2pew4Wa1a51mZJX+xXBO5Vq7UBBwAKKxW40aF+zWlCjLoeELzmiM/DB/Na/bkwVKXaUdaK4qrVRj5hxSdthTgy/PC7s2lsPBy7oxM64KDvop6Nf68Y7qnXkmNdb/8UgfK15rmNfsHi3Wsq//qS7ZNWc/2U506D7nrEmbQcXYGeOBQTcfh93Ts37a5PBj0XtNbyJGXA1dfVGD/7VrcVqiw4IwcV2T2jSaur07F4PtbvngzvcziszS2LDXRtf6ymd6C1viBpdGwLQVFc8z9xo3HnS3vHJVrDWHhaPmk2PiiL9sxwpFZaOW2DuXhRIIMOp7AkXeR9NfclVuNxE15aLUpD0dKDB9Gc9HM5YWGtvm3jpVhcloZciuUFtN15LT5Qu/um63azpLzFViV64+Nkio8tbfY6iGCuLD9ejW25VTZvf7PEvPrWrtnjtz0WfNMx5aLUq7WM5h6nmsKjTXzOhWDWhUDlZpBtZIxWuNqXMYcR8rVkVoSwzCov5v30hoVqpTcBjDqSOCi2HpLh2EYvJlehvK744qlXPXH+F7mt8XGcCReACanleHQ3XcgLpXV44fHm+ssY9h7zeHNarDxgisALDp7L3hV1DPYmsP+w2NbjD9i/zskrtB4Y03QsaWmrp2e1VOHs9RjRc0w2JNbjUlHy1CpdTPyXOtArOkbpqmxrc9WYHqmDIE+Xvhv3zA8lWB8mClHajqO1B4O3KnFc3/ojtZ+6/9iEMLRXEWCrOl4xDMdlprX6tSATGsgy5J6L4MLrX55snGfI6tTawIOAJwtqddpCgGMBDs7N/yfTJl9K8L2Wp21zz2s4eozgarUDDZJFFhzsZK1u19zqWie6diQ3gsHSjAsVYrTxXW8d7lnGGByum7AAYBtOdU4cvfYr1czmJ4pQ70akNczeNNMBwxHhjJybF3DzzZftb9GbQnVdKzw6V/szNdurc/PylHjwG2/9gFo7ERUM9CZQ4WLkQGMPZdQ6V1xDIOdfRtec8nwWYW1Kdl6KWXz2NP/HVzNnJPlWH23bHffrMHOQRENXziQZ3NBvvG3sOWGK79KjfyqWrx6qBTDWhmfzoMrDEyPTP7fSwo8HtsExdVqaM3wjhIzL5WydaNpDTXD4MvzlUj5W44+0YZd8T88Leds9HFBBh2+8T221WcOjlzcePwVVqlwtsTw4beK0Z0CTv+APV7IzUtp+u3O+tVsZ3TesXWbbNay1TA3FR/7bC3e1VrB/Gh+LfIUKrQMcizHZptu735nz2FwW6HCqix+pg1vpD+cjjHGWqjqVAz8jdxtOHJo/SSpwupLlUgM9sWKR0SIDzZ/aV90tkLTdGxs7MEqJYMbFUq0CmE/RAi0ec2Fbw8d9NnfjteqvABcKK1Hzx2FGPWn4cyaSgtdl7nyt1Q3AOr/jGw+07E2KVtrdfplxdVDdUcwDIOMwlpkFNay2oSn3zxqD2ue6bhFz+G75W7smK4y0RnFkcva7JPluFGhwpH8WqtGa9d+VmkKVy/pCjLouKNLZfWYeKQUn1txMFji5QXMPVmu8yxHm/6J4Ei13hZX5bpdWg27TPOTD51t2nhfzXbz2hVZPbbnVKHEitkhrfXRaTkG75Vi8F4pPjp97ybG0byzca9nzXs6XB8HfB5mSr2dadnU2+T8PddZGtX6hyuGz2M2SRRI+DEPHTcXILPQulHVre6YYSNqXnOS7TlVuFmpQoi/FwJ8vDAjoxzVLN3qewE4Yma4foNnK6xs1TJj7xLI69T47UY1EoJ9WZ3Ey9qkbL7AsVhYJ4pqMerPEtSqgOhAb2Q+G4Vm/o5tgGEYfH3h3vO0ry9U4qMeofDy8mLtYutIDvWPPW0MGprf2AzAzmZ4g2e89BgG+Pd+KWvbfXJ3EVKfioSvtxdqlAzeOd7YmUGlM5GiOVw9b6Sg4wQMw2CsA11fHWVwcecp6ui3gUvKlUjYlM/Pxk2w+ZmO/t9e9ne8ePe4DLV3r68F1WqsuVSJ9x4MsS+xu4y13DR2WLD2nuastA7HjDzXY+MwMdfDrLRWjef2l+BgHrfzG/HZyU2/zH1NtC1dVngZ/e3sdbq4Hh+fkWN8+yCDzgz6zdymcNUCQkHHCU5yPFmZpVqx/sHtjm2sVj/TMfOdl5Hv2XyemKM3rtaBO7V4p7ODQcdITaKx44jKigh7vqQOj+8q5qwJSl5vOuW/rLwYOorroKZNP8jK6xgk7yiEj7cX/LyB8ABvbBsQgZ/v+JlIwX5fX6jE1xcq0THMvsu8L0dRh4KOE/zvGnd94AHLd6Tad18Mw7B6hyU0ZjtTGfmMzY4ExtJ29HmGfieRxs/8YflOWlKuRN+dxRa3ccKBm6YFZ/h9/cAYDmYkMOm3G7rTZigZBpe0huuJCmy45fu9mLtLcVaZfcMDcdW85o43uR7P0t14hdazlZJaNbZZ0fVTaKx/pmPbFYjr65W97yo1MvbMpLH242hAa1zfHUc+Zp2XF1JvVeMT/Xf89JvbXLgnLgUdN8J1m7KlY2Vq+r3nSWy+YS9Etv4UNVrVBYZhWL1r9vIC6uy8njd2jTY2/ljjNNmOdtRQMa4/ioLLYBi8YWT0Af3mRVPPeFwBV81rLrzLxF6Wev+cKr7Xds7GuxdCZmvQiAu+93JkKctTFmcU1qHVT/Z1rBh3t2OKsSa01RcrIatVI8+GuWKMUTGM2V6R5J6TxXWQWnFDl1upwtDfLTdpOgNXHQko6PCsuFrFeRPNM/sMXwg1pbDa84KOtEaF1Rcr8eftGuzJtW2qau3z0NQQKM6w/Xo16tUNb5HrO5hXa3cw01alZHCswD2nUGZbWa31x0a6i5apoLpMr127FsuWLUNhYSHatWuHlJQU9O7dm4tNCY74fwXOzgIA4GJZPerVDF4+WOrsrHBix41qfG/k83o1g3/tLMZthX3tWH9L6zCyTVMApt8sd5ZhqVJNUxoXBu+V4hmexzcjziOY3mvbt2/HrFmzsGTJEvTq1Qtr167F888/j8zMTMTHx7O9OWKn3r8WOTsLnBN9fwcHhkaiRaA39t6swfnSeijVjN0BBwC+uajAozENAySeZnm6YkdxGXAa6ffGchcrs2ybNdQT+HPUDuYlk8lYvV3r168fOnbsiGXLlmk+6969O5555hnMmzePte2Ivr/DWlrEfXVq7ofoQG+z0zQTQgwVv9ISfhzUdliNZXV1dTh79iyeeOIJnc+feOIJnDhxgs1NYWgCVfOJZV4AhrUKxL/va5g4K9TBYWaI++JrDEKh4CLgACw3r5WUlEClUiEyMlLn88jISBQVGW/OkUgkdm1rRiwgqwhAehm7g8P7gEF0EwZ3aqiPBQDEN1FjzYM1yKvxxrjzwgv0YajG6RtV+OWOH7qFqvBstBIfXjGcP8ReXUNVOCvnb4ICXy8GSkZ4V0cfMEgKUyOD5fOVTW4xujVLHglT2X1tFovFZr/npCOB/suJDMOYfGHRUgbN2d0OOPzPVQw/bXz6V20tAr0xvFUg3ukcgpZBPjhdXIeiahV25dbgZFEt/tWyCZRqBi/c3xTtRH54Pa0Mqbe4a7/uEOaLiybeFJ7RJQSbrzUMCAoATX29jD60/uihUPSO9sf8M3K7ehU92NwP50uNDz0y/cEQJEf7o3uEP8LujorbJrEO/9p1r3tn13A/DIxvgjc6BWPLtSq8l9EwkGBSpD/6tgzAxPZBiGzijX23azD/tBw3K1V4uW1TtArxRVp+LfbcbCjfkW0Cca1ciTNaw6D83K85dlyvxhatF1dHtg7E6r5h8PLywudn5RbnHZrQLggtAr3xcAt/hIXV4pZChUcfCEa3Uplm/Kn3u4di87UqSMobfosX7m+KwfFNMOaQbgcLf2/giVjDoDs0sQnePiYz6Hr9YHM/BPp6mX17//OezXC+tB4/ShpGqOge4YeEYF8UVqs0z2e2PBmOsABvrMtWQKVm8JK4KW4pVPjgVLlOD6nmAd7oHxeAzdd0X/TtGu6nmVOpS7gfzhmZXynEzwsVZoan0d6nf7cOxLzTpkcV8PbSvXh/8nAo3ugUAoZhcE2uxDcXFfjusgJtm/niyt0y/+6xME2Xb2tEB3qjwEyvy8RgHzzeMsDoaMuNVvcNw4HbNdiSUw0fL2BShyCcKKzDP6X1uP5iDAJ8vLDkfAVS7h5jvVr44y9pHfTfMAj197LYi7FLuB++7RuG7dercb1CibggH8jrGXxrZPJBWzzeMkBndl42tAn1waD4QEzvEqI579nG6jOduro6xMTE4LvvvsPw4cM1n0+fPh0XL17E3r172dqUhkQicShweRoqL9tRmdmOysw2nlRerIYyf39/dO3aFYcOHdL5/NChQ+jZsyebmyKEECJArDevTZ06FZMmTcJDDz2Enj17Yt26dSgoKMBrr73G9qYIIYQIDOtBZ8SIESgtLcXixYtRWFiI9u3bY8uWLUhISGB7U4QQQgSGk44E48ePx/jx47lImhBCiIBRv2BCCCG8oaBDCCGEN6wPg0MIIYSYQjUdQgghvKGgQwghhDcUdAghhPCGgg4hhBDeUNAhhBDCG4tB58svv8Tjjz+O+Ph4tGnTBqNGjcLFixd1lmEYBikpKWjXrh2io6MxZMgQXLp0SWeZH374AUOHDkVCQgJEIhFyc3ONbu/AgQPo378/YmJikJCQgGHDhlnciaysLDz11FOIjo5G+/btsWjRIjCMbqe8rVu34pFHHkFMTAzatm2LiRMnorCw0OF9/+STT/Dwww+jZcuWSExMxLBhwzBt2jSXLrOamhpMnjwZvXv3RkREBIYMGWKwTEFBAcaPH4+HH34YzZs3x+TJk82mqW3t2rV48MEHERUVhcceewzHjx/X+b6oqAiTJ09Gu3btEBMTg4ceegi9e/fmvLzS0tIgEomM/vv111/N7pOlY8ze8rLmGKusrMSMGTPQoUMHREdHo0ePHhg5ciQvxxiXZaYtIyMD4eHhSE5Otlhmx44dw+jRo9G+fXuIRCJs2rTJYBl3LbPJkycbTbdly5YOl5k1+84Gi0EnPT0d48aNw759+7Bz5074+vpi+PDhKCu7NxT5119/jZUrV2LRokU4ePAgIiMj8eyzz6Ki4t7Q81VVVXjiiScwa9Ysk9vavXs3xo4di1GjRuHo0aP4448/8H//939m8yeXy/Hss8+iRYsWOHjwIBYuXIjly5djxYoVmmUyMzMxadIkvPDCC8jIyMCmTZtw+fJlTJgwweF9F4vF+OKLL3D8+HGkpqYiMTERGzZswMiRI122zFQqFZo0aYKJEydiwIABRpepra1F8+bN8c4776BHjx5m09PWOF35e++9h6NHjyIpKQnPP/88bt26BaDhwH7ppZeQk5ODTZs24ejRo6iurkZ+fj5+/fVXTsurZ8+eyM7O1vn37rvvIjg4GE8++aTJfbLmGLO3vKw5xubOnYv9+/fjv//9L06cOIH33nsPf/zxBzp37sz5McZlmTWSyWR4/fXX8dhjj1lVZgqFAh06dMDChQsRGGh8WhN3LbOFCxcapN2qVSudUf3tLTNr9p0VMpmMseXf7du3GW9vb+bnn39mZDIZU1ZWxkRFRTHvv/++Zpn8/HwmODiY+eqrrwzWP3ToEAOAOXfunM7nJSUlTFxcHPP111/blJ8lS5YwISEhTH5+vuazuXPnMjExMUxZWRkjk8mYjz/+mImLi9NZb8WKFUxQUJBD+27s382bNxkAzC+//OKyZab9b8KECUyfPn3MLjNw4EDmhRdesCq9hx56iBkzZozOZ61bt2amTZvGyGQy5vTp0wwAJi0tTfN9aWkpExERwSxbtozT8jL27/7772deeeUVh48xe8vLmmOsffv2zH/+8x+d5Xr37s1MmDCB82OMjzIbOnQoM2vWLGbmzJlM+/btbSqvoKAgZuXKlQafu3uZNf5LTU1lADD79u1zqMxs3XdH/tn8TKeyshJqtRoikQgAkJubi8LCQp0pqgMDA9G7d2+bpqg+e/Ysbt++DX9/f/Tt2xdt27bFs88+i3Pnzpld7+TJk0hOTtaJ3v369UN+fr6m2tuzZ08UFhbi999/B8MwKCkpwfbt29G/f38b9txw3/XV1dVh/fr1CA0NRefOnU2u5+wy44o105XX1jZMOtWkyb0J0by9vREQEICMjAwA3JWXvrS0NFy9ehWvvvqq2eWsOcbYYuwY69WrF1JTU3H79m0AwIkTJ3DhwgX069fP5HpCKbO1a9eiqKgIM2bMsDtPxrhzmWlbv3492rdv7/DUMVztuzE2B51Zs2ahc+fOSEpKAgDNcxFbpqg25saNGwCATz/9FO+99x62bNmCli1bYujQocjPzze5XlFRkdFtN34HAElJSVi7di0mTpyIyMhItGnTBgzD4JtvvrE6f4DhvjdKTU1FbGwsoqKisGrVKuzYsQMtWrQwuZ6zy4wr1kxX3rZtW8THx2PBggUoKytDXV0dli5dijt37mjKhavy0rd+/Xp06tQJ3bp1M7ucNccYW4wdY4sWLULnzp3RqVMnzTO4jz76CIMGDTK5nhDKLCsrC4sWLcKaNWvg48PuNNbuWmbaysvL8dtvv2HMmDF2560RV/tujE1BZ86cOcjMzMTGjRsNDhJbpqg2Rq1umAd2+vTpeOaZZ9C1a1d8/fXXaNasGTZv3gyg4e4lNjYWsbGxeO6558xuW/vzy5cvY9asWZgxYwYOHz6MX375BYWFhXjnnXcAAMePH9ekGxsbiy1btti0748++ijS0tKwf/9+9OvXD6+++ioKCgpcuswcZa7MzO2Xn58fNm7ciOvXr+O+++5DTEwM0tLS0L9/f/j4+HBaXtpKS0uxa9cug7tPe48xSxw5xlavXo0TJ07g559/xuHDh/HZZ5/hgw8+wJ9//ml2PVP5doUyq62txbhx4/Dxxx+jVatWRrdnTZmZ4o5lpm/Lli1QqVQYPXq05jNHyszU9u3dd1Osntpg9uzZ2L59O3bt2qVzkERFRQFoiMRxcXGaz6VSqUHUNKcxnQceeOBe5nx90bp1a00VecuWLVAqG+ZVb2yeadGihUEklkqlAO5F7S+//BLdu3fHW2+9BQDo1KkTmjZtisGDB+ODDz5At27dkJaWpllfP9+m9r1RUFAQWrdujdatW+Phhx9G9+7dsWHDBpSVlblkmbHBWJkFBATAx8fH6O+hvV9du3ZFeno6ysvLUV9fj4iICPTr1w+1tbU4f/48Z+Wl7eeff4a3tzeef/55nc/tPcYssfcYq66uxoIFC/DDDz9g8ODBABqO33/++QfLly/HgQMHOD3GtLFZZgUFBbh8+TKmTp2KqVOnAmi4iWIYBuHh4di6dSuSk5PNlpkp7lpm+tavX49hw4YhLCxM85ml48wULvbdFKtqOjNnzsS2bduwc+dOtG3bVue7xMREREVF6UxRXVNTg4yMDJvaGbt27YqAgABIJBLNZ2q1GtevX0d8fDwAICEhQXNxb+wimJSUhIyMDNTU1GjWO3ToEGJiYpCYmAig4SDUv5tp/JthGAQGBmrSbd26NUJCQqzad1PUajVSU1NdtszYYKzMbJ2uvFmzZoiIiMC1a9dw5swZ3Lp1i9Py0rZx40YMHz4czZo10/nc3mPMEnuPsfr6etTX1xs9fiUSCefHmDY2y6xly5Y4fvw40tLSNP/Gjh2L1q1bIy0tDUlJSWbLzBx3LTNtp0+fxoULFwya1uwtMy723RSLNZ3p06dj8+bN+PHHHyESiTRtf0FBQQgODoaXlxcmT56MJUuWQCwW4/7778cXX3yBoKAgnWpjYWEhCgsLcfXqVQBAdnY2ysvLER8fj7CwMISGhuK1117DwoULERsbi4SEBKxZswbl5eUYOXKkyfw999xzWLRoEaZMmYLp06fj6tWrWLp0Kf7zn/9oqoWDBg3C22+/je+++w79+vVDQUEBZs+ejS5dumguzvbsu1wux7JlyzBo0CBERUWhpKQE3377LW7evImioiL8/PPPLllmQEOTY11dHUpKSqBQKHD+/HkAwIMPPqhZpvEzuVwOLy8vnD9/Hv7+/mjXrp3JdK2ZrvzXX39F8+bNkZCQgKysLEyaNAm+vr7YsGEDp+XVKCMjA5cvX8bSpUvNllEja44xe8vL0jEWGhqKPn36YP78+QgKCkJ8fDyOHTuGjRs3wtfXF1u2bBFkmfn5+aFDhw4660RERCAgIMDgc32VlZXIyckB0HCTdfv2bZw/fx5hYWGIj4932zLTtn79erRp0waPPPKIVWlbKjNr950Vlrq3ATD6b+bMmTrd7WbOnMlERUUxAQEBTO/evZnjx4/rpDNz5kyj6Wh33SsuLmbeeustpkWLFkxISAjTp08f5vDhwxa74B07doxJTk5mAgICmKioKGbWrFkGXQwXLVrEtGvXjgkMDGSioqKY5557jsnKynJo3/Py8pghQ4Yw0dHRjL+/PxMdHc0MHjxYEGUWHx9vNG1L+x8fH28x7S+++IKJj49n/P39mS5dujB79uzR+X7hwoVMbGws4+fnx8TFxfFaXjKZjBk9ejTzwAMP2NTN05pjzJ7ysmbfs7OzmRdffJGJiYlhmjRpwojFYrcpM/28WNNleteuXUbzrN1N3Z3L7NatW0xQUBAzf/58q9O1psys2Xc2/tF8OoQQQnhDY68RQgjhDQUdQgghvKGgQwghhDcUdAghhPCGgg4hhBDeUNAhhBDCGwo6hNhp06ZNZifXI4QYoqBDCM9SU1ORkpLi7GwQ4hQUdAjh2b59+7Bo0SJnZ4MQp6CgQwghhDcUdAixwqlTpzBgwABERUWhU6dO+OqrrzRznTTau3cvRo0ahfbt26NFixbo1KkT5s2bp5ktFQAmT56M77//HgAgEok0/7SfC/3yyy/o168fYmJikJCQgFGjRuHy5cv87CghHLN6Ph1CPNXly5cxfPhwhISEYPr06fD398cPP/yAoKAgneV+/PFH+Pj4YOLEiRCJRDhx4gSWL1+OO3fuYO3atQCA1157DXfu3MHRo0exevVqzboREREAgKVLl+Kjjz7C008/jdGjR0OhUGDt2rUYOHAgjhw5YnLCM0KEggb8JMSCl19+GampqTh58iTuu+8+AA2TW3Xv3h1yuRznzp1DYmIiqqqq0LRpU511Fy9ejM8++wwXLlxAbGwsAGDatGn4/vvvIZPJdJa9desWunXrhvfeew+zZ8/WfF5QUICkpCQMGzYMK1as4HZnCeEYNa8RYoZKpcKBAwcwaNAgTcABGmom+nMWNQYctVqN8vJylJSUoHfv3mAYBufOnbO4rV27dkGpVOLf//43SkpKNP/8/PzQo0cPHD16lN2dI8QJqHmNEDOkUimqqqogFosNvrv//vt1/r506RI+/PBDpKeno7q6Wue78vJyi9u6du0agIYZJI3Rr0URIkQUdAgxo7GzgP7MjdrfAQ1B5emnn0ZgYCA++OAD3HfffQgMDEReXh6mTJkCtVptcVuNy2zbtg2+voanprc3NUwQ4aOgQ4gZkZGRaNq0Ka5cuWLwXWPNBADS0tIglUqxe/dunSmEteecb2QsgAHQNN/FxcWZneKaECGjWydCzPDx8cETTzyB1NRUXL9+XfO5VCrF1q1bdZYDdGs/arUaK1euNEizsZlMvyPBsGHD4Ovri5SUFKM1I6lU6tC+EOIKqKZDiAVz5szBwYMHMXjwYIwfPx5+fn744YcfEB8fr3lW06tXLzRv3hyTJ0/GpEmT4Ovri507d6KystIgvW7dugEAZsyYgSeffBK+vr4YNGgQWrVqhfnz52Pu3Ll48skn8fTTTyMsLAy3bt3C/v370aNHD3z11Ve87jshbKMu04RY4cSJE3j//fdx/vx5REZGYty4cYiMjMQbb7yh6TJ9+vRpzTJBQUEYNmwYxo4diz59+mDlypV46aWXADT0iJszZw527NiB4uJiTe+2xMREAA3D5Cxfvhznzp2DUqlETEwMevXqhfHjx6N79+7OLAZCHEZBhxBCCG/omQ4hhBDeUNAhhBDCGwo6hBBCeENBhxBCCG8o6BBCCOENBR1CCCG8oaBDCCGENxR0CCGE8IaCDiGEEN5Q0CGEEMKb/wc0kvr9kW1bXgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Design a query to retrieve the last 12 months of precipitation data and plot the results. \n",
    "# Starting from the last data point in the database. \n",
    "# prev_year = dt.date(2017, 8, 23)\n",
    "\n",
    "# Calculate the date one year from the last date in data set.\n",
    "prev_year = dt.date(2017, 8, 23) - dt.timedelta(days=365)\n",
    "\n",
    "# Perform a query to retrieve the data and precipitation scores\n",
    "results = []\n",
    "# Save the query results as a Pandas DataFrame and set the index to the date column\n",
    "# results = session.query(Measurement.date, Measurement.prcp)\n",
    "# print(results.all())\n",
    "\n",
    "# Sort the dataframe by date\n",
    "#results = session.query(Measurement.date, Measurement.prcp).filter(Measurement.date >= prev_year)\n",
    "results = session.query(Measurement.date, Measurement.prcp).filter(Measurement.date >= prev_year).all()\n",
    "#print(results)\n",
    "\n",
    "#This is to add it to a dataFrame\n",
    "df = pd.DataFrame(results, columns=['date','precipitation'])\n",
    "df.set_index(df['date'], inplace=True)\n",
    "\n",
    "# This will sort the index\n",
    "df = df.sort_index()\n",
    "\n",
    "# we'll convert the DataFrame to strings, and then we'll set our index to \"False.\"\n",
    "# This will allow us to print the DataFrame without the index.\n",
    "#print(df.to_string(index=False))\n",
    "\n",
    "# Use Pandas Plotting with Matplotlib to plot the data\n",
    "df.plot()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>precipitation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>2021.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>0.177279</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>0.461190</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>0.020000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>0.130000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>6.700000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       precipitation\n",
       "count    2021.000000\n",
       "mean        0.177279\n",
       "std         0.461190\n",
       "min         0.000000\n",
       "25%         0.000000\n",
       "50%         0.020000\n",
       "75%         0.130000\n",
       "max         6.700000"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Use Pandas to calcualte the summary statistics for the precipitation data\n",
    "df.describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[(9,)]"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# How many stations are available in this dataset?\n",
    "# session.query()\n",
    "# session.query(func.count(Station.station))\n",
    "session.query(func.count(Station.station)).all()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('USC00519281', 2772),\n",
       " ('USC00519397', 2724),\n",
       " ('USC00513117', 2709),\n",
       " ('USC00519523', 2669),\n",
       " ('USC00516128', 2612),\n",
       " ('USC00514830', 2202),\n",
       " ('USC00511918', 1979),\n",
       " ('USC00517948', 1372),\n",
       " ('USC00518838', 511)]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# What are the most active stations?\n",
    "# List the stations and the counts in descending order.\n",
    "#session.query(Measurement.station, func.count(Measurement.station))\n",
    "\n",
    "#Group by\n",
    "#session.query(Measurement.station, func.count(Measurement.station)).\\\n",
    "    #group_by(Measurement.station)\n",
    "    \n",
    "# ORDER BY\n",
    "# session.query(Measurement.station, func.count(Measurement.station)).\\\n",
    "# group_by(Measurement.station).order_by(func.count(Measurement.station).desc())\n",
    "\n",
    "# ALL RESULTS FOR OUR QUERY\n",
    "session.query(Measurement.station, func.count(Measurement.station)).\\\n",
    "group_by(Measurement.station).order_by(func.count(Measurement.station).desc()).all()\n",
    "\n",
    "\n",
    "# In the left column is the station ID, and on the right are the counts for each station.\n",
    "# The counts indicate which stations are most active.\n",
    "# We can also see which stations are the least active."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[(54.0, 85.0, 71.66378066378067)]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Using the station id from the previous query, calculate the lowest temperature recorded, \n",
    "# highest temperature recorded, and average temperature most active station?\n",
    "\n",
    "# calculate the minimum, maximum, and average temperatures\n",
    "#session.query(func.min(Measurement.tobs), func.max(Measurement.tobs), func.avg(Measurement.tobs))\n",
    "\n",
    "# With the minimum, maximum, and average in our query, we now need to add one filter\n",
    "# If you look at the outcome of the previous query, you can see that the most active station is USC00519281\n",
    "# session.query(func.min(Measurement.tobs), func.max(Measurement.tobs), func.avg(Measurement.tobs)).\\\n",
    "# filter(Measurement.station == 'USC00519281')\n",
    "\n",
    "#add the .all() function to return our results as a list.\n",
    "session.query(func.min(Measurement.tobs), func.max(Measurement.tobs), func.avg(Measurement.tobs)).\\\n",
    "filter(Measurement.station == 'USC00519281').all()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " tobs\n",
      " 77.0\n",
      " 77.0\n",
      " 80.0\n",
      " 80.0\n",
      " 75.0\n",
      " 73.0\n",
      " 78.0\n",
      " 77.0\n",
      " 78.0\n",
      " 80.0\n",
      " 80.0\n",
      " 78.0\n",
      " 78.0\n",
      " 78.0\n",
      " 73.0\n",
      " 74.0\n",
      " 80.0\n",
      " 79.0\n",
      " 77.0\n",
      " 80.0\n",
      " 76.0\n",
      " 79.0\n",
      " 75.0\n",
      " 79.0\n",
      " 78.0\n",
      " 79.0\n",
      " 78.0\n",
      " 78.0\n",
      " 76.0\n",
      " 74.0\n",
      " 77.0\n",
      " 78.0\n",
      " 79.0\n",
      " 79.0\n",
      " 77.0\n",
      " 80.0\n",
      " 78.0\n",
      " 78.0\n",
      " 78.0\n",
      " 77.0\n",
      " 79.0\n",
      " 79.0\n",
      " 79.0\n",
      " 79.0\n",
      " 75.0\n",
      " 76.0\n",
      " 73.0\n",
      " 72.0\n",
      " 71.0\n",
      " 77.0\n",
      " 79.0\n",
      " 78.0\n",
      " 79.0\n",
      " 77.0\n",
      " 79.0\n",
      " 77.0\n",
      " 78.0\n",
      " 78.0\n",
      " 78.0\n",
      " 78.0\n",
      " 77.0\n",
      " 74.0\n",
      " 75.0\n",
      " 76.0\n",
      " 73.0\n",
      " 76.0\n",
      " 74.0\n",
      " 77.0\n",
      " 76.0\n",
      " 76.0\n",
      " 74.0\n",
      " 75.0\n",
      " 75.0\n",
      " 75.0\n",
      " 75.0\n",
      " 71.0\n",
      " 63.0\n",
      " 70.0\n",
      " 68.0\n",
      " 67.0\n",
      " 77.0\n",
      " 74.0\n",
      " 77.0\n",
      " 76.0\n",
      " 76.0\n",
      " 75.0\n",
      " 76.0\n",
      " 75.0\n",
      " 73.0\n",
      " 75.0\n",
      " 73.0\n",
      " 75.0\n",
      " 74.0\n",
      " 75.0\n",
      " 74.0\n",
      " 75.0\n",
      " 73.0\n",
      " 75.0\n",
      " 73.0\n",
      " 73.0\n",
      " 74.0\n",
      " 70.0\n",
      " 72.0\n",
      " 70.0\n",
      " 67.0\n",
      " 67.0\n",
      " 69.0\n",
      " 70.0\n",
      " 68.0\n",
      " 69.0\n",
      " 69.0\n",
      " 66.0\n",
      " 65.0\n",
      " 68.0\n",
      " 62.0\n",
      " 75.0\n",
      " 70.0\n",
      " 69.0\n",
      " 76.0\n",
      " 76.0\n",
      " 74.0\n",
      " 73.0\n",
      " 71.0\n",
      " 74.0\n",
      " 74.0\n",
      " 72.0\n",
      " 71.0\n",
      " 72.0\n",
      " 74.0\n",
      " 69.0\n",
      " 67.0\n",
      " 72.0\n",
      " 70.0\n",
      " 64.0\n",
      " 63.0\n",
      " 63.0\n",
      " 62.0\n",
      " 70.0\n",
      " 70.0\n",
      " 62.0\n",
      " 62.0\n",
      " 63.0\n",
      " 65.0\n",
      " 69.0\n",
      " 77.0\n",
      " 70.0\n",
      " 74.0\n",
      " 69.0\n",
      " 72.0\n",
      " 71.0\n",
      " 69.0\n",
      " 71.0\n",
      " 71.0\n",
      " 72.0\n",
      " 72.0\n",
      " 69.0\n",
      " 70.0\n",
      " 66.0\n",
      " 65.0\n",
      " 69.0\n",
      " 68.0\n",
      " 68.0\n",
      " 68.0\n",
      " 59.0\n",
      " 60.0\n",
      " 70.0\n",
      " 73.0\n",
      " 75.0\n",
      " 64.0\n",
      " 59.0\n",
      " 59.0\n",
      " 62.0\n",
      " 68.0\n",
      " 70.0\n",
      " 73.0\n",
      " 79.0\n",
      " 75.0\n",
      " 65.0\n",
      " 70.0\n",
      " 74.0\n",
      " 70.0\n",
      " 70.0\n",
      " 71.0\n",
      " 71.0\n",
      " 71.0\n",
      " 69.0\n",
      " 61.0\n",
      " 67.0\n",
      " 65.0\n",
      " 72.0\n",
      " 71.0\n",
      " 73.0\n",
      " 72.0\n",
      " 77.0\n",
      " 73.0\n",
      " 67.0\n",
      " 62.0\n",
      " 64.0\n",
      " 67.0\n",
      " 66.0\n",
      " 81.0\n",
      " 69.0\n",
      " 66.0\n",
      " 67.0\n",
      " 69.0\n",
      " 66.0\n",
      " 68.0\n",
      " 65.0\n",
      " 74.0\n",
      " 69.0\n",
      " 72.0\n",
      " 73.0\n",
      " 72.0\n",
      " 71.0\n",
      " 76.0\n",
      " 77.0\n",
      " 76.0\n",
      " 74.0\n",
      " 68.0\n",
      " 73.0\n",
      " 71.0\n",
      " 74.0\n",
      " 75.0\n",
      " 70.0\n",
      " 67.0\n",
      " 71.0\n",
      " 67.0\n",
      " 74.0\n",
      " 77.0\n",
      " 78.0\n",
      " 67.0\n",
      " 70.0\n",
      " 69.0\n",
      " 69.0\n",
      " 74.0\n",
      " 78.0\n",
      " 71.0\n",
      " 67.0\n",
      " 68.0\n",
      " 67.0\n",
      " 76.0\n",
      " 69.0\n",
      " 72.0\n",
      " 76.0\n",
      " 68.0\n",
      " 72.0\n",
      " 74.0\n",
      " 70.0\n",
      " 67.0\n",
      " 72.0\n",
      " 60.0\n",
      " 65.0\n",
      " 75.0\n",
      " 70.0\n",
      " 75.0\n",
      " 70.0\n",
      " 79.0\n",
      " 75.0\n",
      " 70.0\n",
      " 67.0\n",
      " 74.0\n",
      " 70.0\n",
      " 75.0\n",
      " 76.0\n",
      " 77.0\n",
      " 74.0\n",
      " 74.0\n",
      " 74.0\n",
      " 69.0\n",
      " 68.0\n",
      " 76.0\n",
      " 74.0\n",
      " 71.0\n",
      " 71.0\n",
      " 74.0\n",
      " 74.0\n",
      " 74.0\n",
      " 74.0\n",
      " 80.0\n",
      " 74.0\n",
      " 72.0\n",
      " 75.0\n",
      " 80.0\n",
      " 76.0\n",
      " 76.0\n",
      " 77.0\n",
      " 75.0\n",
      " 75.0\n",
      " 75.0\n",
      " 75.0\n",
      " 72.0\n",
      " 74.0\n",
      " 74.0\n",
      " 74.0\n",
      " 76.0\n",
      " 74.0\n",
      " 75.0\n",
      " 73.0\n",
      " 79.0\n",
      " 75.0\n",
      " 72.0\n",
      " 72.0\n",
      " 74.0\n",
      " 72.0\n",
      " 72.0\n",
      " 77.0\n",
      " 71.0\n",
      " 73.0\n",
      " 76.0\n",
      " 77.0\n",
      " 76.0\n",
      " 76.0\n",
      " 79.0\n",
      " 81.0\n",
      " 76.0\n",
      " 78.0\n",
      " 77.0\n",
      " 74.0\n",
      " 75.0\n",
      " 78.0\n",
      " 78.0\n",
      " 69.0\n",
      " 72.0\n",
      " 74.0\n",
      " 74.0\n",
      " 76.0\n",
      " 80.0\n",
      " 80.0\n",
      " 76.0\n",
      " 76.0\n",
      " 76.0\n",
      " 77.0\n",
      " 77.0\n",
      " 77.0\n",
      " 82.0\n",
      " 75.0\n",
      " 77.0\n",
      " 75.0\n",
      " 76.0\n",
      " 81.0\n",
      " 82.0\n",
      " 81.0\n",
      " 76.0\n",
      " 77.0\n",
      " 82.0\n",
      " 83.0\n",
      " 77.0\n",
      " 77.0\n",
      " 77.0\n",
      " 76.0\n",
      " 76.0\n",
      " 79.0\n"
     ]
    }
   ],
   "source": [
    "# Choose the station with the highest number of temperature observations.\n",
    "# Query the last 12 months of temperature observation data for this station and plot the results as a histogram\n",
    "#session.query(Measurement.tobs)\n",
    "\n",
    "# filter out all the stations except the most active station\n",
    "# results = session.query(Measurement.tobs).\\\n",
    "# filter(Measurement.station == 'USC00519281')\n",
    "\n",
    "# apply another filter to consider only the most recent year.\n",
    "results = session.query(Measurement.tobs).\\\n",
    "filter(Measurement.station == 'USC00519281').\\\n",
    "filter(Measurement.date >= prev_year).all()\n",
    "\n",
    "#Print results\n",
    "# print(results)\n",
    "\n",
    "# CONVERT TO DF\n",
    "df = pd.DataFrame(results, columns=['tobs'])\n",
    "print(df.to_string(index=False))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEQCAYAAAD2/KAsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAjGUlEQVR4nO3df1QU5f4H8PeKoYDoenMdVARLVhGTBBQM85iYIFIiSqmZdckfhVqSioqVptlBEbpiF/eqaJraLUUsb7+8FXQNQ+TquWlltoSiFiyCLoLBMWC+f3TcrysiLOzyLLvv1zmc484zO/PZp4F3z8wzswq9Xi+DiIiojXUQXQAREdknBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMIBaSavVii7BqrA/jLE/jLE/jNl7fzCAiIhICAYQEREJwQAiIiIhGEBERCRER9EFEBGJVltbi+vXr7f5fjt37oyKioo23685ubi4oGPHlkUJA4iI7FptbS0qKyuhVCqhUCjadN+dOnVC586d23Sf5iTLMvR6PVxdXVsUQjwFR0R27fr160LCxxYoFAoolcoWjx4ZQERk9xg+LdeavuMpOCJqM8NznIGcX9tkX/qYPm2yH2o5joCIiEgIBhARkZ2IiIhAfHy86DIMeAqOiOgOlO+0zanCm0w9ZRgREQEfHx9s2LDBQhVZHkdAREQkBAOIiKidiY2NxdGjR7Ft2zYolUoolUoUFRXh6NGjGDt2LCRJglqtRkJCAm7cuGH03traWixbtgyenp7w9PTEa6+9hvr6ekP7oUOHEBwcDDc3N/Tr1w8TJkxAaWmpRT4HA4iIqJ1Zt24dAgMDMWPGDJw9exZnz57FPffcgyeeeAK+vr44cuQI3n77bRw4cACrV682eu/+/ftRX1+PL774Ahs3bsSuXbuwefNmAIBOp8OsWbMwffp05OXl4dNPP8W0adMs9jl4DYiIqJ3p1q0b7rnnHjg7O0OSJADAG2+8AUmSkJKSgg4dOmDgwIFYtWoVXn75ZbzyyitwdnYGAEiShKSkJCgUCgwYMAAFBQXYvHkzFixYgOLiYvzxxx+IjIyEh4cHAMDHx8din4MjICIiG3D27FkMHz4cHTr8/5/1hx56CDdu3EBhYaFh2bBhw4xuHg0MDMRvv/2Ga9euYciQIXjkkUcQHByMmTNnYvv27SgrK7NYzQwgIiIbIMtyo08laO7TChwcHHDw4EFkZmZi8ODB2L17N/z9/XH69GlzlmrAACIiaoccHR1RV1dneO3t7Y38/HyjCQW5ublwdHTEfffdZ1h24sQJyLJseJ2fn49evXqha9euAP4Mq8DAQCxfvhzZ2dno1asXDh48aJHPwAAiImqHPDw8cOLECRQVFaG8vByzZs1CSUkJFi9ejLNnz+Lw4cNYvXo15syZY7j+AwAlJSVYvnw5tFotPvroI2zatAnz5s0D8GcYbdiwASdPnsTFixfx6aef4tdff8XAgQMt8hk4CYGIqB168cUXERsbixEjRqC6uhrfffcd9u/fj5UrV2LUqFHo1q0boqOjsXLlSqP3PfHEE6ivr8fYsWOhUCgwc+ZMQwB17doVeXl52Lp1KyoqKtCnTx/Ex8dj6tSpFvkMCr1eLze9GjVGq9VCrVaLLsNqsD+MsT+MteXTBZr7ZIGKigp069bNwtXcWU1NTbv+PqCbWtqHPAVHRERCMICIiEgIBhAREQkhNIBKSkrwwgsvoH///pAkCUFBQcjJyTG0y7KMxMREeHt7w83NDREREThz5ozAiomIyFyEBZBer0dYWBhkWca+ffuQl5eHpKQkqFQqwzqpqalIS0vD+vXrkZWVBZVKhaioKFRWVooqm4iIzETYNOxNmzbBzc0NW7ZsMSzr16+f4d+yLEOj0SAuLg6RkZEAAI1GA7VajYyMDMTExLR1yURko+72FAG6u1tvajWVsBHQJ598goCAAMTExMDLywsPP/wwtm7davgwRUVF0Ol0CAkJMbzHyckJwcHByMvLE1U2EdkYFxcX6PX6Vv0htVeyLEOv18PFxaVF7xc2Ajp//jy2b9+OefPmIS4uDqdPn8ayZcsAAHPnzoVOpwMAo1NyN18XFxc3ul2tVmu5oq1on9aM/WGM/XEr56ZXMRNT+/3atWtGD/KkptXX1+PGjRuNPrC0qXvghAVQfX09/Pz8sGrVKgDAgw8+iMLCQqSnp2Pu3LmG9W4fFjc1VG7rm/54o6Ex9ocx9sdtctruRtT20O/2fnwIi3tJkho8X2jAgAG4dOmSoR1Ag2/iKysrazAqIiKi9kdYAI0YMQIFBQVGywoKCtC3b18AgKenJyRJQnZ2tqG9pqYGubm5CAoKatNaiYjI/IQF0Lx585Cfn4/k5GQUFhbiww8/xNatWzF79mwAf556i42NxcaNG3Ho0CH8+OOPmDdvHlxcXBAdHS2qbCIiMhNh14D8/f2xd+9erFmzBhs2bIC7uztWrFhhCCAAWLhwIaqrqxEfHw+9Xo+AgABkZmbC1dVVVNlERGQmQr+OISwsDGFhYY22KxQKJCQkICEhoQ2rIiKitsA5h0REJAQDiIiIhGAAERGREAwgIiISQugkBCISry2/JpvoVhwBERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMICIiEgIBhAREQnBp2ETkU1qy6d862P6tNm+bAlHQEREJAQDiIiIhGAAERGREAwgIiISggFERERCCAugxMREKJVKo58BAwYY2mVZRmJiIry9veHm5oaIiAicOXNGVLlERGRmQkdAarUaZ8+eNfx8++23hrbU1FSkpaVh/fr1yMrKgkqlQlRUFCorKwVWTERE5iI0gDp27AhJkgw/PXr0APDn6Eej0SAuLg6RkZHw8fGBRqNBVVUVMjIyRJZMRERmIjSAzp8/j0GDBsHX1xfPPfcczp8/DwAoKiqCTqdDSEiIYV0nJycEBwcjLy9PULVERGROwp6EMGzYMGzevBlqtRplZWXYsGEDQkNDcezYMeh0OgCASqUyeo9KpUJxcfFdt6vVai1WszXt05qxP4xZf384iy6g3WvNf2PrPz5aTq1W37VdWACNGzfO6PWwYcMwdOhQvPfeexg+fDgAQKFQGK0jy3KDZbdr6gObm1arbfN9WjP2h7F20R85bffIGlvV0v/G7eL4sCCrmYbdpUsXeHt7o7CwEJIkAQBKS0uN1ikrK2swKiIiovbJagKopqYGWq0WkiTB09MTkiQhOzvbqD03NxdBQUECqyQiInMRdgru1Vdfxfjx4+Hu7m64BvT7779j+vTpUCgUiI2NRUpKCtRqNby8vJCcnAwXFxdER0eLKpmIiMxIWAD99ttvmD17NsrLy9GjRw8MGzYMX3zxBTw8PAAACxcuRHV1NeLj46HX6xEQEIDMzEy4urqKKpmIiMxIWADt2LHjru0KhQIJCQlISEhoo4qIiKgtWc01ICIisi8MICIiEoIBREREQjCAiIhICAYQEREJwQAiIiIhTA4gvV5vgTKIiMjemBxAAwcOxDPPPINPPvkEf/zxhyVqIiIiO2ByAM2dOxcnTpzA008/jYEDB2LJkiXIz8+3RG1ERGTDTA6gN954A99//z0+/PBDhIWFYd++fQgLC4O/vz+SkpIMXypHRER0Ny2ahKBQKDB69GhoNBr8/PPP2LZtG9RqNTZs2AB/f3+Eh4dj586dvF5ERESNavUsuM6dO2PKlCl4+eWXER4eDlmWcezYMbz88ssYNGgQ4uPjce3aNXPUSkRENqRVDyP95Zdf8MEHH2D//v0oKipCz549sWDBAkyfPh2Ojo7YuXMn0tPT8euvv+K9994zV81ERGQDTA6g8vJyHDhwAPv27cPJkyfh6OiICRMmICkpCWPHjkWHDv8/qFq7di0kSUJiYqJZiyYiovbP5ADy9vZGbW0tAgMD8dZbbyEqKgrdunVrdH21Wo0ePXq0qkgiIrI9JgfQwoULMX36dPTv379Z648fPx7jx483uTAiIrJtJgfQq6++aok6iIjIzpg8C2737t2YOXNmo+3PPPMMJxwQEVGTTA6g7du3Q5KkRtvd3NyQnp7eqqKIiMj2mRxAv/zyCwYPHtxo+6BBg1BQUNCqooiIyPaZHEAKhQLl5eWNtl+5cgX19fWtKoqIiGyfyQH04IMPYv/+/aipqWnQVl1djf3798PX19csxRERke0yOYAWLVoErVaLsLAwfPTRR9BqtSgoKMBHH32E8PBwaLVaLFq0yBK1EhGRDTE5gMaMGYPNmzfj3LlziImJQVBQEAIDA/HXv/4VhYWFePvtt/Hoo4+aXEhKSgqUSiXi4+MNy2RZRmJiIry9veHm5oaIiAicOXPG5G0TEZH1adGz4KZNm4aIiAhkZWXh/PnzkGUZ9913H0JCQuDq6mry9vLz87Fr164GkxtSU1ORlpaGtLQ0qNVqJCUlISoqCvn5+S3aD1FrKN/5tQXvcgZyTH+fPqZPC/ZF1L60+GGkrq6uiIyMbHUBFRUVmDNnDt5++20kJSUZlsuyDI1Gg7i4OMN+NBoN1Go1MjIyEBMT0+p9ExGROC0OoMrKSly6dAlXr16FLMsN2keOHNms7dwMmNGjRxsFUFFREXQ6HUJCQgzLnJycEBwcjLy8PAYQEVE7Z3IA6fV6LF26FAcPHkRdXR2AP0crCoXC6N9Xrlxpclu7du1CYWEhtmzZ0qBNp9MBAFQqldFylUqF4uJiU8smIiIrY3IAxcXF4eOPP8acOXMwcuRIKJXKFu1Yq9VizZo1+Oyzz+Do6NjoejeD7aZbw66x7bY1Efu0ZrbbH85ttqeWXW8iUVpzzNvu78uf34ZwNyYH0Jdffonnn38eb775ZouLAoDjx4+jvLwcDz30kGFZXV0dvv32W+zYsQPHjh0DAJSWlsLd3d2wTllZWYNR0a2a+sDmptVq23yf1sym+6MFkwnIPrT0mLfp35dmMDmAHB0dm/1VDHcTEREBPz8/o2Xz589H//79sWjRInh5eUGSJGRnZ8Pf3x8AUFNTg9zcXKxZs6bV+yciIrFMDqDIyEh88cUXeO6551q1Y6VS2eD0nbOzM7p37w4fHx8AQGxsLFJSUqBWq+Hl5YXk5GS4uLggOjq6VfsmIiLxTA6gF198EbNmzcILL7yAWbNmoW/fvnBwcGiw3t1OkzXXwoULUV1djfj4eOj1egQEBCAzM5P3ABER2QCFXq9vOIf6Lrp37w6FQtHkZIDmzIKzBfZ+Dvd2ttwfnBhAjWnpjcO2/PvSHCaPgJYuXXrX4CEiImoOkwMoISHBEnUQEZGdMflhpLeqq6vDlStXUFtba656iIjITrQogE6ePIlJkyahd+/e8PLywtGjRwEA5eXlePLJJ/Gf//zHrEUSEZHtMTmAjh8/jgkTJuDcuXOYNm2a0XPg7r33XlRVVWH37t1mLZKIiGyPyQH0xhtvoH///sjLy8PKlSsbtI8aNQr//e9/zVIcERHZLpMD6OTJk3j66afRuXPnO86G69Onj+FBokRERI0xOYA6dOiADh0af5tOp4OTk1OriiIiIttncgANHToUn3/++R3bbty4gf379yMwMLDVhRERkW0zOYAWLVqEI0eOYMGCBTh9+jQAoKSkBF9++SUmTpyIc+fOYfHixWYvlIiIbIvJN6KOGTMGW7ZsQXx8PN577z0Afz40VJZldOvWDenp6Rg+fLjZCyUiItvSoq/kjo6OxoQJE5CVlYXCwkLU19fjvvvuw9ixY9GlSxdz10hERDaoRQEE/PnVCY899pg5ayEiIjticgBdvHixWev17dvX5GKIiMh+mBxAvr6+zXoatr18HQMREbWMyQH097//vUEA1dXVoaioCO+//z569uyJ2bNnm61AIiKyTSYH0IwZMxpti4uLQ0hICKqqqlpVFBER2b5WfR3D7bp06YIZM2Zg8+bN5twsERHZILMGEADcc889KC4uNvdmiYjIxpg1gE6fPo1//OMfGDhwoDk3S0RENshss+AqKipw7do1dOnSBWlpaWYpjoiIbJfJATRy5MgGAaRQKKBUKnH//fdjypQpUCqV5qqPiIhslMkBpNFoLFEHERHZGbNPQiAiImoOk0dA69evN3knCoUCS5cuNVq2bds2vPPOO4ZH+3h7e2PJkiUICwsDAMiyjHXr1mHXrl3Q6/UICAhAcnIyBg0aZPL+iYjI+pgcQOvWrTNcA5Jl2ajtbstvD6DevXtj9erV6N+/P+rr6/HPf/4TM2bMwNdff40HHngAqampSEtLQ1paGtRqNZKSkhAVFYX8/Hy4urqaWjYREVkZk0/B/fTTT3jggQcQHR2NrKwsXLhwARcuXMBXX32FKVOm4IEHHsBPP/2Eq1evGn7u9Fy4iIgIjBs3Dvfffz+8vLzw2muvoUuXLsjPz4csy9BoNIiLi0NkZCR8fHyg0WhQVVWFjIwMs3xwIiISy+QAWrp0KTw9PbF161b4+fnB1dUVrq6u8Pf3x7Zt2+Dh4dFgtNOUuro6HDhwANevX0dgYCCKioqg0+kQEhJiWMfJyQnBwcHIy8sztWQiIrJCJp+Cy87OxqpVqxptHzNmDFavXt2sbf3www8IDQ1FTU0NXFxcsGfPHgwePNgQMiqVymh9lUrV5FMWtFpts/ZtTiL2ac1stz+cRRdAVqo1x7zt/r4AarX6ru0mB5CDgwO+//77RttPnTqFDh2aN7BSq9X45ptvUFFRgUOHDiE2NhYff/yxof32+41kWW7yqyCa+sDmptVq23yf1sym+yPnV9EVkJVq6TFv078vzWDyKbiJEyfi3XffRUpKCiorKw3LKysrkZycjD179iAyMrJZ23J0dMT9998PPz8/rFq1CkOGDMHmzZshSRIAoLS01Gj9srKyBqMiIiJqn0weAa1duxbnzp3D2rVrkZiYiJ49e0KhUECn06Gurg4jR47E2rVrW1RMfX09bty4AU9PT0iShOzsbPj7+wMAampqkJubizVr1rRo20REZF1MDiBXV1ccOnQIn332GQ4fPoxLly5BlmWEhoYiNDQU4eHhzdrO66+/jtDQUPTp08cwuy0nJwf79u2DQqFAbGwsUlJSoFar4eXlheTkZLi4uCA6OtrkD0lEZEnKd1p6eta5Rad29TF9Wrg/62JyAN0UHh7e7LC5E51Oh7lz56K0tBRdu3bF4MGDkZGRgbFjxwIAFi5ciOrqasTHxxtuRM3MzOQ9QERENkKh1+vlpldr6OLFizh69CguX76MqKgouLu7o7a2FlevXkX37t3RsWOLs61dsfeLiLez5f5o+f/lEpmXXY+AVqxYga1bt6Kurg4KhQK+vr5wd3fH77//Dn9/fyxfvhzz5883d61ERGRDTJ4Ft2nTJmg0GsyfPx8ffvih0WN3unbtioiICKOp1ERERHdicgDt2rULTz75JFavXo0hQ4Y0aB88eDB++eUXsxRHRES2y+QAunTpEoKDgxttd3V1RUVFRauKIiIi22fyNaC//OUvKCkpabT9hx9+QK9evVpVFFFzcFIAUftm8ggoNDQUu3btQnl5eYO27777Dnv27EFERIRZiiMiIttlcgCtWLECHTp0QHBwMF5//XUoFArs3bsXzz33HMaNG4fevXsjPj7eErUSEZENMTmAJEnC119/jfHjx+Nf//oXZFnG/v378eWXX2Lq1Kn497//DaVSaYFSiYjIlph0DeiPP/7A8ePH4ebmhtTUVKSmpqKsrAz19fXo0aNHs5+CTUREZFJiODg4YNKkScjKyjIs69GjB3r27MnwISIik5iUGh06dICHhweqqqosVQ8REdkJk4ct8+bNw86dO3H58mVL1ENERHbC5PuAqqqq4OLiAn9/f0RERKBfv35wcnIyWkehUOCll14yW5FERGR7TH4advfu3ZveqEKBK1eutLio9sSWn/7cEm3ZH7wRleyVXT0Ne/HixXj66afh5+eH7777DrIs4+rVq+jatavdfO0CERGZV7PSY8eOHRgxYgT8/Pzg4eGBK1euwM/PDwcPHsTo0aMtXSMREdmgFs+dvvVrGIiIiEzFm3eIiEgIBhAREQnR7BkE58+fx4kTJwAA165dA/DnjKcuXbrccf2AgAAzlEdERLaqWdOwu3fvDoVCYbRMluUGy25dzmnY9onTsIksz66mYaelpVm6DiIisjPNCqCnnnrK0nUQEZGd4SQEIiISQlgAvfXWWxgzZgz69u2L/v37Y+rUqfjxxx+N1pFlGYmJifD29oabmxsiIiJw5swZQRUTEZE5CQugnJwczJo1C4cPH8ahQ4fQsWNHTJo0CVevXjWsk5qairS0NKxfvx5ZWVlQqVSIiopCZWWlqLKJiMhMhD3ILTMz0+j1li1b4OHhgWPHjiE8PByyLEOj0SAuLg6RkZEAAI1GA7VajYyMDMTExIgom4iIzMRqrgFVVVWhvr4eSqUSAFBUVASdToeQkBDDOk5OTggODkZeXp6gKomIyFys5lHWy5cvx5AhQxAYGAgA0Ol0AACVSmW0nkqlQnFxcaPb0Wq1livSivZpzdquP5zbaD9E1qW9/M1p6p5AqwigFStW4NixY/j888/h4OBg1NbcG2BvauubQnkjqrE27Y8c3ohK9slW/uYIPwWXkJCAAwcO4NChQ+jXr59huSRJAIDS0lKj9cvKyhqMioiIqP0RGkDLli1DRkYGDh06hAEDBhi1eXp6QpIkZGdnG5bV1NQgNzcXQUFBbV0qERGZmbBTcEuWLMEHH3yAPXv2QKlUGq75uLi4oEuXLlAoFIiNjUVKSgrUajW8vLyQnJwMFxcXREdHiyqbiIjMRFgApaenA4BhivVNy5YtQ0JCAgBg4cKFqK6uRnx8PPR6PQICApCZmQlXV9c2r5eIiMyrWU/DpsZxEoIxPg2byPJs5WnYwichEBGRfWIAERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMICIiEiIjqILINsyPMcZyPlVdBlE1A5wBEREREIIDaCjR49i2rRpGDRoEJRKJfbu3WvULssyEhMT4e3tDTc3N0RERODMmTOCqiUiInMSGkDXr1+Hj48P1q1bBycnpwbtqampSEtLw/r165GVlQWVSoWoqChUVlYKqJaIiMxJaACFhoZi5cqViIyMRIcOxqXIsgyNRoO4uDhERkbCx8cHGo0GVVVVyMjIEFQxERGZi9VeAyoqKoJOp0NISIhhmZOTE4KDg5GXlyewMiIiMgerDSCdTgcAUKlURstVKhVKS0tFlERERGZk9dOwFQqF0WtZlhssu5VWq7V0SVaxT+vlLLoAIpvXXv7mqNXqu7ZbbQBJkgQAKC0thbu7u2F5WVlZg1HRrZr6wOam1WrbfJ9WjfcAEVmcrfzNsdpTcJ6enpAkCdnZ2YZlNTU1yM3NRVBQkMDKiIjIHISOgKqqqlBYWAgAqK+vx6VLl3Dq1Cl0794dffv2RWxsLFJSUqBWq+Hl5YXk5GS4uLggOjpaZNlERGQGCr1eL4va+TfffIPHH3+8wfLp06dDo9FAlmWsW7cOO3fuhF6vR0BAAJKTk+Hj4yOg2jvjKThjynd4Co7I0vQxfUSXYBZCA8gWMICMMYCILM9WAshqrwEREZFts9pZcEREdGdteabBkqMtjoCIiEgIBhAREQnBU3CC8GI9Edk7joCIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREJyGfYuWTY125nfgEBG1AEdAREQkBAOIiIiEYAAREZEQDCAiIhKCAUREREIwgIiISAgGEBERCcEAIiIiIRhAREQkBAOIiIiEYAAREZEQDCAiIhKiXQRQeno6fH19IUkSRo8ejW+//VZ0SURE1EpWH0CZmZlYvnw5Fi9ejCNHjiAwMBBPPPEELl68KLo0IiJqBasPoLS0NDz11FN49tlnMXDgQGzYsAGSJGHHjh2iSyMiolaw6u8DunHjBv73v//hxRdfNFoeEhKCvLw8s+9PH9PH7NskIqI7s+oRUHl5Oerq6qBSqYyWq1QqlJaWCqqKiIjMwaoD6CaFQmH0WpblBsuIiKh9seoAuvfee+Hg4NBgtFNWVtZgVERERO2LVQeQo6Mjhg4diuzsbKPl2dnZCAoKElQVERGZg1VPQgCA+fPn4/nnn0dAQACCgoKwY8cOlJSUICYmRnRpRETUClY9AgKAyZMnIzExERs2bMCoUaNw7Ngx7Nu3Dx4eHm1WQ0lJCV544QX0798fkiQhKCgIOTk5hnZZlpGYmAhvb2+4ubkhIiICZ86cabP62lpT/REbGwulUmn08+ijjwqs2HKGDBnS4LMqlUo8+eSTAOzv2GiqP+zp2ACAuro6rF271nAjva+vL9auXYva2lrDOvZ2jNzK6kdAADB79mzMnj1byL71ej3CwsIwYsQI7Nu3D/feey+KioqMrkGlpqYiLS0NaWlpUKvVSEpKQlRUFPLz8+Hq6iqkbktpTn8AwCOPPIItW7YYXjs6OrZ1qW0iOzsbdXV1htclJSV45JFHMGnSJAD2dWwATfcHYD/HBgBs3LgR6enp0Gg08PHxwQ8//IDY2Fg4Ojpi6dKlAOzvGLlVuwggkTZt2gQ3NzejX5h+/foZ/i3LMjQaDeLi4hAZGQkA0Gg0UKvVyMjIsLlThU31x02dOnWCJEltWJkYPXr0MHq9e/duuLq6YtKkSXZ3bAB374+b7OXYAIDjx49j/PjxCA8PBwB4enoiPDwcJ06cAGB/fz9uZ/Wn4ET75JNPEBAQgJiYGHh5eeHhhx/G1q1bIcsyAKCoqAg6nQ4hISGG9zg5OSE4ONgiN8uK1lR/3JSbmwsvLy8EBATgpZdewuXLlwVV3HZkWcbu3bsxdepUODs7292xcbvb++Mmezo2RowYgZycHPz8888AgJ9++gnffPMNxo0bB8D+/n7cjiOgJpw/fx7bt2/HvHnzEBcXh9OnT2PZsmUAgLlz50Kn0wHAHW+WLS4ubvN6La2p/gCARx99FI8//jg8PT1x4cIFrF27FhMnTsTXX3+NTp06iSzforKzs1FUVISZM2cCgN0dG7e7vT8A+zs24uLiUFVVhaCgIDg4OKC2thZLliwxXFKw92OEAdSE+vp6+Pn5YdWqVQCABx98EIWFhUhPTzf8wQXs52bZ5vTHlClTDOsPHjwYQ4cOxZAhQ3D48GFMnDhRSN1tYdeuXfD394evr6/Rcns5Nm53p/6wt2MjMzMT77//PtLT0+Ht7Y3Tp09j+fLl8PDwwDPPPGNYz16PEZ6Ca4IkSRg4cKDRsgEDBuDSpUuGdgB2c7NsU/1xJ7169ULv3r1RWFho6fKEuXz5Mj799FM8++yzhmX2dmzc6k79cSe2fmysXLkSCxYswJQpUzB48GBMmzYN8+fPx9/+9jcA9n2MAAygJo0YMQIFBQVGywoKCtC3b18Af15UlCTJ6GbZmpoa5Obm2uTNsk31x52Ul5ejuLjYpi887927F506dcLkyZMNy+zt2LjVnfrjTmz92Pj999/h4OBgtMzBwQH19fUA7PsYAXgKrknz5s1DaGgokpOTMXnyZJw6dQpbt27Fa6+9BuDPoXNsbCxSUlKgVqvh5eWF5ORkuLi4IDo6WnD15tdUf1RVVWHdunWYOHEiJEnChQsXsGbNGqhUKjz22GOCq7cMWZbx7rvvYvLkyUbTZu3t2Lipsf6wx2Nj/Pjx2LhxIzw9PeHt7Y1Tp04hLS0N06ZNA2C/x8hNCr1eLze9mn07fPgw1qxZg4KCAri7u2POnDl4/vnnDedoZVnGunXrsHPnTuj1egQEBCA5ORk+Pj6CK7eMu/VHdXU1ZsyYgVOnTqGiogKSJGHUqFF45ZVX4O7uLrp0izhy5AgmTpyIr776CgEBAUZt9nZsAI33hz0eG5WVlXjzzTfx8ccfo6ysDJIkYcqUKVi6dCk6d+4MwD6PkZsYQEREJASvARERkRAMICIiEoIBREREQjCAiIhICAYQEREJwQAiIiIhGEBERCQEA4iIiIRgABERkRD/B3hy5f7MpxBrAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# create an histogram with 12 bins\n",
    "# Bin\" refers to each rectangular column in the plot\n",
    "df.plot.hist(bins=12)\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEQCAYAAAD2/KAsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAkmElEQVR4nO3deVhU9f4H8Pc4yurokOIIIpgyofBACKZcvBQq4YKC4oaZmUsWakYpqGmWREGKmXURMTSNNJXESvNmdYVMw+W6pGnRKIpK7DYsKikwvz+8zs9RWQZmOLO8X88zz+N8v2fOfM6XL7w9c5YRKZVKFYiIiFpZG6ELICIi88QAIiIiQTCAiIhIEAwgIiISBAOIiIgEwQAiIiJBMICIiEgQDCAiIhIEA0jPFAqF0CUYDI7FHRyHOzgOd5jzODCAiIhIEAwgIiISBAOIiIgEwQAiIiJBtBW6AENSU1OD69ev63SdVlZWKC8v1+k6DZGtrS3atuV0IqKm41+M/6mpqUFlZSWkUilEIpHO1mtpaQkrKyudrc8QqVQqKJVKSCQShhARNRk/gvuf69ev6zx8zIVIJIJUKtX53iMRmTYG0D0YPs3HsSMibfHzEiJqlPSTfD2s1QY4qI/16pZyWjehSzBZ3AMiIiJBMIDMVEhICKKjo4Uug4jMGD+Ca4R+PnqoX3N290NCQuDu7o6VK1fqoSIiIv3gHhAREQmCAWTkIiMjcejQIXz88ceQSqWQSqXIy8vDoUOHMGTIEMhkMsjlcixevBi3bt3SeG1NTQ0WLlwIFxcXuLi44I033kBdXZ26/+uvv4a/vz+6du2KHj16YMSIESguLm7tTSQiE8UAMnIJCQno378/Jk+ejJycHOTk5KBdu3YYP348vLy8cODAAXz00UfYuXMnli9frvHa9PR01NXV4fvvv8cHH3yAzZs3Y+3atQCAoqIizJgxA5MmTcKRI0ewd+9eRERECLGJRGSieAzIyHXs2BHt2rWDjY0NZDIZAODtt9+GTCbDqlWr0KZNG7i5ueHNN9/Eq6++iiVLlsDGxgYAIJPJsGLFCohEIjz22GM4f/481q5di7lz56KgoAC3b99GWFgYnJ2dAQDu7u6CbScRmR5B94AOHTqEiIgI9OnTB1KpFFu2bNHoj4yMVH+sdPcRFBSksczff/+N6Oho9OzZE46OjoiIiEB+vuFfW6BPOTk5eOKJJ9Cmzf//eP/xj3/g1q1byM3NVbf169dP4wLS/v37488//0RFRQU8PT0RGBgIf39/TJkyBRs2bEBpaWmrbgcRmTZBA+j69etwd3dHQkICrK2tH7pMYGCg+qOlnJwcpKena/QvXrwYu3fvxoYNG7B3715UVlZi4sSJqK2tbY1NMEgqlareOxM09Y4FYrEYu3btQkZGBjw8PJCWlgYfHx+cOXNGl6USkRkTNICCg4OxbNkyhIWFafxv/V6WlpaQyWTqh52dnbqvvLwcaWlpiI2NxaBBg+Dt7Y2UlBScPXsWWVlZrbQVwrOwsNAI3N69e+PYsWMaJxRkZ2fDwsICjz76qLrt+PHjUKlU6ufHjh2Dg4MDOnToAOBOWPXv3x+LFi1CZmYmHBwcsGvXrlbYIiIyBwZ/EkJ2djZcXV3h6+uLefPmoaSkRN136tQp3L59G4MHD1a3OTk5wc3NDUeOHBGiXEE4Ozvj+PHjyMvLQ1lZGWbMmIHCwkLMnz8fOTk52LdvH5YvX44XXnhBffwHAAoLC7Fo0SIoFAp89dVX+PDDDzF79mwAd8Jo5cqVOHHiBK5cuYK9e/ciPz8fbm5uQm0mEZkYgz4JISgoCKNGjYKLiwsuX76MuLg4hIaGIisrC5aWliguLoZYLEanTp00Xmdvb9/g6cIKheKBNisrK1haWup8G7RVXV2t9WtmzZqFefPmwc/PDzdv3sTRo0exZcsWxMbGIiAgAB06dEB4eDhiYmLU66+rq0N4eDhu3bqFIUOGQCQSYdKkSZg+fTqqq6thaWmJn3/+GSkpKaioqICjoyOioqIQFhZWb40VFRWNnqb9sLE3R8Y3DjaNL2KiWuNnZXzzoWnkcnmD/SKlUqlqcIlW0q1bN6xYsQKTJ0+ud5mCggJ4enpi48aNCA0NRXp6Ol566SWUlpZqHNsYOXIk5HI5Vq9e3eT3Ly8vR8eOHVu0DQ9TXV1t8t8HdFdjY6hQKBqdkObAGMehte8IYkj0fTNSY5wPumLwH8Hdy8HBAY6Ojuozubp06YLa2lqUlZVpLFdaWgp7e3shSiQioiYyqgAqKytDQUGB+noXb29vtGvXDpmZmepl8vPzkZOTgwEDBghVJhERNYGgx4CqqqrUezN1dXW4evUqTp8+DTs7O9jZ2SEhIQGhoaGQyWS4fPkyYmNjYW9vj5EjRwK4cxHmlClTsGzZMtjb28POzg5LliyBh4cHAgMDBdwyIiJqjKABdPLkSYwaNUr9PD4+HvHx8Zg0aRLef/99nDt3Dtu2bUN5eTlkMhkCAgLwySefQCKRqF/z7rvvQiwWY9q0aaiursaTTz6JdevWQSwWC7FJRETURIIGUEBAAJRKZb39GRkZja7DysoKK1eu5FcREBEZGaM6BqRv916USdrh2BGRthhA/2NrawulUsk/pM2gUqmgVCpha2srdClEZEQM+kLU1tS2bVtIJBJUVFTodL0VFRXqW9uYMolEgrZtOZ2IqOn4F+Mebdu21fnFqMXFxejevbtO10lEZAr4ERwREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREgmAAERGRIBhAREQkCAYQEREJggFERESCYAAREZEgGEBERCQIBhAREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREgmAAERGRIBhAREQkCAYQEREJQtAAOnToECIiItCnTx9IpVJs2bJF3Xf79m28+eab8Pf3h6OjI9zc3DBz5kxcuXJFYx0hISGQSqUaj+nTp7f2phARkZYEDaDr16/D3d0dCQkJsLa21ui7ceMGfvnlFyxYsAA//vgjtm7divz8fIwbNw41NTUay06ePBk5OTnqx+rVq1tzM4iIqBnaCvnmwcHBCA4OBgDMnj1bo69jx4748ssvNdpWr14NPz8/5OTkwMPDQ91uY2MDmUym93qJiEh3jOoYUGVlJQBAKpVqtO/cuRM9e/aEn58fli5dql6OiIgMl6B7QNq4desWli5dimHDhqFbt27q9vHjx6N79+7o2rUrfv/9dyxfvhy//vrrA3tP91IoFK1QsXDvZ8g4FncY3zjYCF2AYFrjZ2V886Fp5HJ5g/1GEUA1NTWYNWsWysvL8fnnn2v0Pf/88+p/e3h4oEePHhgyZAhOnToFb2/vh66vsUHRJYVC0arvZ8g4FncY5TgczBe6AsHo+2dllPNBRwz+I7iamhrMmDEDZ8+exVdffYVHHnmkweX79u0LsViM3NzcVqqQiIiaw6D3gG7fvo3p06fjt99+w549e5p0osHZs2dRW1vLkxKIiAycoAFUVVWl3lOpq6vD1atXcfr0adjZ2cHBwQFTp07FyZMn8fnnn0MkEqGoqAgA0KFDB1hbW+PixYvYsWMHgoOD8cgjjyAnJwdLly6Fl5cX/Pz8hNw0IiJqhKABdPLkSYwaNUr9PD4+HvHx8Zg0aRIWLVqEvXv3AgACAwM1XpeUlITJkyejXbt2+PHHH7Fu3Tpcv34d3bp1Q3BwMBYtWgSxWNyam0JERFoSNIACAgKgVCrr7W+oDwCcnJzUIUVERMbF4E9CICIi08QAIiIiQTCAiIhIEAwgIiISBAOIiIgEwQAiIiJBMICIiEgQDCAiIhIEA4iIiATBACIiIkEwgIiISBAMICIiEgQDiIiIBMEAIiIiQWgdQI19RQIREVFTaB1Abm5ueO655/DNN9/g9u3b+qiJiIjMgNYBNGvWLBw/fhzPPvss3NzcsGDBAhw7dkwftRERkQnTOoDefvtt/Prrr/jyyy8xdOhQ7NixA0OHDoWPjw9WrFiBS5cu6aFMIiIyNc06CUEkEuGpp55CcnIy/vjjD3z88ceQy+VYuXIlfHx8MHz4cGzatInHi4iIqF4tPgvOysoKY8eOxauvvorhw4dDpVLh8OHDePXVV9GnTx9ER0ejoqJCF7USEZEJaduSF1+4cAHbt29Heno68vLy0KVLF8ydOxeTJk2ChYUFNm3ahNTUVOTn52Pr1q26qpmIiEyA1gFUVlaGnTt3YseOHThx4gQsLCwwYsQIrFixAkOGDEGbNv+/UxUXFweZTIb4+HidFk1ERMZP6wDq3bs3ampq0L9/f7z//vsYM2YMOnbsWO/ycrkcnTt3blGRRERkerQOoFdeeQWTJk1Cr169mrT8sGHDMGzYMK0LIyIi06Z1AC1dulQfdRARkZnR+iy4tLQ0TJkypd7+5557rsknHBw6dAgRERHo06cPpFIptmzZotGvUqkQHx+P3r17o2vXrggJCcFvv/2msczff/+N6Oho9OzZE46OjoiIiEB+fr62m0VERK1M6wDasGEDZDJZvf1du3ZFampqk9Z1/fp1uLu7IyEhAdbW1g/0r1mzBklJSXjvvfewf/9+2NvbY8yYMaisrFQvs3jxYuzevRsbNmzA3r17UVlZiYkTJ6K2tlbbTSMiolakdQBduHABHh4e9fb36dMH58+fb9K6goODsWzZMoSFhWmcPQfc2ftJTk5GVFQUwsLC4O7ujuTkZFRVVeGLL74AAJSXlyMtLQ2xsbEYNGgQvL29kZKSgrNnzyIrK0vbTSMiolakdQCJRCKUlZXV23/t2jXU1dW1qCgAyMvLQ1FREQYPHqxus7a2hr+/P44cOQIAOHXqFG7fvq2xjJOTE9zc3NTLEBGRYdL6JITHH38c6enpmDt3LqysrDT6bt68ifT0dHh5ebW4sKKiIgCAvb29Rru9vT0KCgoAAMXFxRCLxejUqdMDyxQXF9e7boVC0eL6tNHa72fIOBZ3GN842AhdgGBa42dlfPOhaeRyeYP9WgfQa6+9hrFjx2Lo0KF47bXX4O7uDpFIhLNnz2L16tVQKBTYvn17swu+n0gk0niuUqkeaLtfY8s0Nii6pFAoWvX9DBnH4g6jHIeD5ntij75/VkY5H3RE6wAaNGgQ1q5di5iYGEybNk3drlKpIJFI8NFHHyEoKKjFhd090aG4uBhOTk7q9tLSUvVeUZcuXVBbW4uysjKNi11LS0vh7+/f4hqI7pJ+oss/wDZm/Qed6K5m3QsuIiICISEh2L9/Py5dugSVSoVHH30UgwcPhkQi0UlhLi4ukMlkyMzMhI+PDwCguroa2dnZiI2NBQB4e3ujXbt2yMzMxPjx4wEA+fn5yMnJwYABA3RSBxER6Uezb0YqkUgQFhbWojevqqpCbm4uAKCurg5Xr17F6dOnYWdnh+7duyMyMhKrVq2CXC6Hq6srEhMTYWtri3HjxgEAOnbsiClTpmDZsmWwt7eHnZ0dlixZAg8PDwQGBraoNiIi0q9mB1BlZSWuXr2Kv/76CyqV6oH+gQMHNrqOkydPYtSoUern8fHxiI+Px6RJk5CcnIxXXnkFN2/eRHR0NJRKJXx9fZGRkaGxl/Xuu+9CLBZj2rRpqK6uxpNPPol169ZBLBY3d9OIiKgViJRK5YPp0QClUomYmBjs2rVLfbHnvQf97/772rVruq/WCJnzAcb7GfNY6PYYEBkT5bRuel2/Mf9etJTWe0BRUVHYs2cPXnjhBQwcOBBSqVQPZRERkanTOoB++OEHvPjii3jnnXf0UQ8REZkJre+EYGFh0eSvYiAiIqqP1gEUFhaG77//Xh+1EBGRGdE6gF5++WUUFhbipZdewrFjx1BYWIiSkpIHHkRERA3R+hiQr68vRCIRTp06hR07dtS7HM+CIyKihmgdQDExMY3ei42IiKgxWgfQ4sWL9VEHERGZGa2PAd2rtrYW165dQ01Nja7qISIiM9GsADpx4gRGjx4NR0dHuLq64tChQwCAsrIyTJgwAT/++KNOiyQiItOjdQAdPXoUI0aMwMWLFxEREaFxH7hOnTqhqqoKaWlpOi2SiIhMj9YB9Pbbb6NXr144cuQIli1b9kB/QEAA/vvf/+qkOCIiMl1aB9CJEyfw7LPPwsrK6qFnw3Xr1k39ddpERET10TqA2rRpgzZt6n9ZUVERrK2tW1QUERGZPq0DyNvbG99+++1D+27duoX09HT079+/xYUREZFp0zqAXnvtNRw4cABz587FmTNnAACFhYX44YcfEBoaiosXL2L+/Pk6L5SIiEyL1heiDho0CCkpKYiOjsbWrVsBAJGRkVCpVOjYsSNSU1PxxBNP6LxQIiIyLc36Su5x48ZhxIgR2L9/P3Jzc1FXV4dHH30UQ4YMQfv27XVdIxERmaBmBRAA2NjYYOTIkbqshYiIzIjWAXTlypUmLde9e3etiyEiIvOhdQB5eXk16W7Y/DoGIiJqiNYB9K9//euBAKqtrUVeXh62bduGLl26YObMmTorkIiITJPWATR58uR6+6KiojB48GBUVVW1qCgiIjJ9Lfo6hvu1b98ekydPxtq1a3W5WiIiMkE6DSAAaNeuHQoKCnS9WiIiMjE6DaAzZ85g3bp1cHNz08n6PD09IZVKH3hMmDABwJ0LYO/vCwoK0sl7ExGRfunsLLjy8nJUVFSgffv2SEpK0klxmZmZqK2tVT8vLCxEYGAgRo8erW4LDAxESkqK+rmFhYVO3puIiPRL6wAaOHDgAwEkEokglUrRs2dPjB07FlKpVCfFde7cWeN5WloaJBKJRgBZWlpCJpPp5P2IiKj1aB1AycnJ+qijUSqVCmlpaZg4cSJsbGzU7dnZ2XB1dUXHjh0xcOBAvPHGG7C3txekRiIiarpm34qntWVmZiIvLw9TpkxRtwUFBWHUqFFwcXHB5cuXERcXh9DQUGRlZcHS0lLAaomIqDEipVKp0uYF7733nvZvIhIhJiZG69fda+rUqbhy5Qr2799f7zIFBQXw9PTExo0bERoaWu9yCoWiRbWQ+XnioE3jC5FJOvbPG0KXYLTkcnmD/VrvASUkJKiPAalUmtnVUHtLAqikpAR79+5FYmJig8s5ODjA0dERubm5DS7X2KDokkKhaNX3M2RGPRYH84WugASi7zlr1L8XLaR1AP3+++8YP348evfujcjISLi6ugK4M4jJycnIyclBenq6Tk8M2LJlCywtLREeHt7gcmVlZSgoKOBJCURERkDr64BiYmLg4uKC9evXo2/fvpBIJJBIJPDx8cHHH38MZ2fnFn/cdi+VSoVPP/0U4eHhkEgk6vaqqiosXboUR48eRV5eHn766SdERETA3t6eXxNBRGQEtA6gzMxMBAYG1ts/aNAgZGZmtqQmDT/99BNyc3MxdepUjXaxWIxz587hmWeeQb9+/dR7Y999951GUBERkWHS+iM4sViMX3/9td7+06dPo00b3d1g4cknn4RSqXyg3draGhkZGTp7HyIial1aJ0VoaCg+/fRTrFq1CpWVler2yspKJCYm4rPPPkNYWJhOiyQiItOj9R5QXFwcLl68iLi4OMTHx6NLly4QiUQoKipCbW0tBg4ciLi4OH3USkREJkTrAJJIJPj666/x73//G/v27cPVq1ehUqkQHByM4OBgDB8+XB91EhGRiWn2nRCGDx/OsCEiomZrdgBduXIFhw4dQklJCcaMGQMnJyfU1NTgr7/+gp2dHdq2NZq7/BARkQCalRKvv/461q9fj9raWohEInh5ecHJyQk3btyAj48PFi1ahDlz5ui6ViIiMiFanwX34YcfIjk5GXPmzMGXX36pcdudDh06ICQkBHv27NFpkUREZHq0DqDNmzdjwoQJWL58OTw9PR/o9/DwwIULF3RSHBERmS6tA+jq1avw9/evt18ikaC8vLxFRRERkenTOoAeeeQRFBYW1tt/9uxZODg4tKgoIiIyfVoHUHBwMDZv3oyysrIH+n755Rd89tlnCAkJ0UlxRERkurQOoNdffx1t2rSBv78/3nrrLYhEImzZsgXTp0/H008/DUdHR0RHR+ujViIiMiFaB5BMJkNWVhaGDRuG3bt3Q6VSIT09HT/88AMmTpyI7777DlKpVA+lEhGRKdHqOqDbt2/j6NGj6Nq1K9asWYM1a9agtLQUdXV16Ny5s07vgk1ERKZNq8QQi8UYPXo09u/fr27r3LkzunTpwvAhIiKtaJUabdq0gbOzM6qqqvRVDxERmQmtd1tmz56NTZs2oaSkRB/1EBGRmdD6XnBVVVWwtbWFj48PQkJC0KNHD1hbW2ssIxKJMG/ePJ0VSUREpkfrAHrrrbfU/96+fftDl2EAERFRY5oUQPPnz8ezzz6Lvn374pdffoFKpcJff/2FDh068GsXiIioWZqUHhs3boSfnx/69u0LZ2dnXLt2DX379sWuXbvw1FNP6btGIiIyQc0+d/rer2EgIiLSFi/eISIiQTCAiIhIEE0+g+DSpUs4fvw4AKCiogIAoFAo0L59+4cu7+vrq4PyiIjIVImUSmWjB3Ps7OwgEok02lQq1QNt97Zfu3ZNd1UaMYVCAblcLnQZBsGYx0L6Sb7QJZBAlNO66XX9xvx70VJN2gNKSkrSdx0PFR8fj/fee0+jrUuXLvjjjz8A3Am7hIQEbN68GUqlEr6+vkhMTESfPn2EKJeIiLTQpAB65pln9F1HveRyOfbs2aN+LhaL1f9es2YNkpKSkJSUBLlcjhUrVmDMmDE4duwYJBKJEOUSEVETGfxJCG3btoVMJlM/OnfuDODO3k9ycjKioqIQFhYGd3d3JCcno6qqCl988YXAVRMRUWMMPoAuXbqEPn36wMvLC9OnT8elS5cAAHl5eSgqKsLgwYPVy1pbW8Pf3x9HjhwRqFoiImoqg76PTr9+/bB27VrI5XKUlpZi5cqVCA4OxuHDh1FUVAQAsLe313iNvb09CgoKGlyvQqHQW82G8H6GzHjHwkboAkggrTFnjff3omGNnVxh0AH09NNPazzv168fvL29sXXrVjzxxBMA0OSz8+7VmmecmPMZLvcz6rE4yLPgzJW+56xR/160kMF/BHev9u3bo3fv3sjNzYVMJgMAFBcXayxTWlr6wF4REREZHqMKoOrqaigUCshkMri4uEAmkyEzM1OjPzs7GwMGDBCwSiIiagqD/ghu6dKlGDZsGJycnNTHgG7cuIFJkyZBJBIhMjISq1atglwuh6urKxITE2Fra4tx48YJXToRETXCoAPozz//xMyZM1FWVobOnTujX79++P777+Hs7AwAeOWVV3Dz5k1ER0erL0TNyMjgNUBEREagSbfioeYz5wOM9zPmseCteMwXb8WjP0Z1DIiIiEwHA4iIiATBACIiIkEwgIiISBAGfRYcGZaWH4i34R0FiEiNe0BERCQIBhAREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREguB1QEREDdD/jWgN+/o4fd6MlXtAREQkCAYQEREJggFERESCYAAREZEgGEBERCQIBhAREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREgmAAERGRIAw6gN5//30MGjQI3bt3R69evTBx4kScO3dOY5nIyEhIpVKNR1BQkEAVExFRUxn03bAPHjyIGTNmwMfHByqVCu+++y5Gjx6NI0eOwM7OTr1cYGAgUlJS1M8tLCyEKJeIiLRg0AGUkZGh8TwlJQXOzs44fPgwhg8frm63tLSETCZr7fKIiKgFDPojuPtVVVWhrq4OUqlUoz07Oxuurq7w9fXFvHnzUFJSIkyBRETUZCKlUqkSuoimev7553HhwgVkZWVBLBYDAHbu3Alra2u4uLjg8uXLiIuLQ11dHbKysmBpafnQ9SgUitYs22Q8cdBG6BKIqJUd++eNZr9WLpc32G80AfT6668jIyMD3377LXr06FHvcgUFBfD09MTGjRsRGhraegXWQ6FQNPpDMBb6/2ZIIjI0+vxGVIM+BnTX4sWLkZGRgd27dzcYPgDg4OAAR0dH5Obmtk5xRETULAYfQAsXLkRGRgb27NmDxx57rNHly8rKUFBQwJMSiIgMnEEH0IIFC7B9+3Z89tlnkEqlKCoqAgDY2tqiffv2qKqqQkJCAkJDQyGTyXD58mXExsbC3t4eI0eOFLh6IiJqiEEHUGpqKgAgLCxMo33hwoVYvHgxxGIxzp07h23btqG8vBwymQwBAQH45JNPIJFIhCiZiIiayKADSKlUNthvbW39wLVCRERkHIzqOiAiIjIdDCAiIhIEA4iIiATBACIiIkEwgIiISBAMICIiEoRBn4ZtaJp3LzQb4CDvoUZEdD/uARERkSAYQEREJAgGEBERCYIBREREgmAAERGRIBhAREQkCAYQEREJggFERESCYAAREZEgGEBERCQIBhAREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREgmAAERGRIBhAREQkCJMJoNTUVHh5eUEmk+Gpp57Czz//LHRJRETUAJMIoIyMDCxatAjz58/HgQMH0L9/f4wfPx5XrlwRujQiIqqHSQRQUlISnnnmGUydOhVubm5YuXIlZDIZNm7cKHRpRERUj7ZCF9BSt27dwqlTp/Dyyy9rtA8ePBhHjhzR6Xspp3XT6fqIiMyZ0e8BlZWVoba2Fvb29hrt9vb2KC4uFqgqIiJqjNEH0F0ikUjjuUqleqCNiIgMh9EHUKdOnSAWix/Y2yktLX1gr4iIiAyH0QeQhYUFvL29kZmZqdGemZmJAQMGCFQVERE1xuhPQgCAOXPm4MUXX4Svry8GDBiAjRs3orCwENOmTRO6NCIiqofR7wEBQHh4OOLj47Fy5UoEBATg8OHD2LFjB5ydnVuthsLCQrz00kvo1asXZDIZBgwYgIMHD6r7VSoV4uPj0bt3b3Tt2hUhISH47bffWq2+1tLYOERGRkIqlWo8goKCBKxY9zw9PR/YRqlUigkTJgAwn7nQ2DiYw1wAgNraWsTFxakvlPfy8kJcXBxqamrUy5jLnLifSewBAcDMmTMxc+ZMQd5bqVRi6NCh8PPzw44dO9CpUyfk5eVpHINas2YNkpKSkJSUBLlcjhUrVmDMmDE4duwYJBKJIHXrWlPGAQACAwORkpKifm5hYdHapepVZmYmamtr1c8LCwsRGBiI0aNHAzCPuQA0Pg6A6c8FAPjggw+QmpqK5ORkuLu74+zZs4iMjISFhQViYmIAmM+cuJ/JBJCQPvzwQ3Tt2lXjF6lHjx7qf6tUKiQnJyMqKgphYWEAgOTkZMjlcnzxxRcm81FhY+Nwl6WlJWQyWStW1ro6d+6s8TwtLQ0SiQSjR482m7kANDwOd5n6XACAo0ePYtiwYRg+fDgAwMXFBcOHD8fx48cBmM/fh4cxiY/ghPbNN9/A19cX06ZNg6urK/75z39i/fr1UKlUAIC8vDwUFRVh8ODB6tdYW1vD399f5xfLCqmxcbgrOzsbrq6u8PX1xbx581BSUiJQxfqnUqmQlpaGiRMnwsbGxmzmwv3uH4e7zGEu+Pn54eDBg/jjjz8AAL///jt++uknPP300wDM5+/Dw3APSAcuXbqEDRs2YPbs2YiKisKZM2ewcOFCAMCsWbNQVFQEAA+9WLagoKDV69WXxsYBAIKCgjBq1Ci4uLjg8uXLiIuLQ2hoKLKysmBpaSlk+XqRmZmJvLw8TJkyBQDMZi7c7/5xAMxnLkRFRaGqqgoDBgyAWCxGTU0NFixYoD5kYK5zAmAA6URdXR369u2LN998EwDw+OOPIzc3F6mpqeo/vIDpXyzblHEYO3asenkPDw94e3vD09MT+/btQ2hoqCB169PmzZvh4+MDLy8vjXZTnwv3e9g4mMtcyMjIwLZt25CamorevXvjzJkzWLRoEZydnfHcc8+plzO3OQHwIzidkMlkcHNz02h77LHHcPXqVXU/AJO/WLaxcXgYBwcHODo6Ijc3V9/ltbqSkhLs3bsXU6dOVbeZy1y418PG4WFMdS4sW7YMc+fOxdixY+Hh4YGIiAjMmTMHq1evBmCec+IuBpAO+Pn54fz58xpt58+fR/fu3QHcOegok8k0Lpatrq5Gdna2SV0s29g4PExZWRkKCgpM8kD0li1bYGlpifDwcHWbucyFez1sHB7GVOfCjRs3IBaLNdrEYjHq6uoAmOecuIsfwenA7NmzERwcjMTERISHh+P06dNYv3493njjDQB3dq0jIyOxatUqyOVyuLq6IjExEba2thg3bpzA1etOY+NQVVWFhIQEhIaGQiaT4fLly4iNjYW9vT1GjhwpcPW6pVKp8OmnnyI8PFzjNFpzmQt31TcO5jQXhg0bhg8++AAuLi7o3bs3Tp8+jaSkJERERAAwvzlxL5FSqVQ1vhg1Zt++fYiNjcX58+fh5OSEF154AS+++KL6M1yVSoWEhARs2rQJSqUSvr6+SExMhLu7u8CV61ZD43Dz5k1MnjwZp0+fRnl5OWQyGQICArBkyRI4OTkJXbpOHThwAKGhofjPf/4DX19fjT5zmQtA/eNgTnOhsrIS77zzDvbs2YPS0lLIZDKMHTsWMTExsLKyAmBec+JeDCAiIhIEjwEREZEgGEBERCQIBhAREQmCAURERIJgABERkSAYQEREJAgGEBERCYIBREREgmAAERGRIP4Ppw8keQyuN9EAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df.plot.hist(bins=5)\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEQCAYAAAD2/KAsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAhhklEQVR4nO3df1RUZf4H8PcAoYCD4+o4WCAWjCCsJJBgmEfF/EEUiGJqah1LM8SSVFDKLM0O/oA2dXFWRVfS2hLCMvvh2krHMCRWT2mt0RhBajAIOsioHAXu94+O83VEhIGBZ5x5v87hHOd57tz7mYfLvH3m/hiZXq+XQERE1MUcRBdARET2iQFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAygDtJqtaJLsCocD1McD1McD1P2Ph4MICIiEoIBREREQjCAiIhICAYQEREJ4SS6gK7S0NCAy5cvW3y93bt3R21trcXXa03c3Nzg5GQ3uwoRdRG7eFdpaGhAXV0dFAoFZDKZRdfdrVs3dO/e3aLrtCaSJEGv10MulzOEiMii7OIjuMuXL3dK+NgDmUwGhULRKbNHIrJvdhFAABg+HcCxI6LOwM9UiMgmKf55rkPP18++z0KVUEvsZgZERETWhQFkg6Kjo5GcnCy6DCKiO7Lrj+A6OkU3V3um9NHR0QgICMD69es7oSIiInE4AyIiIiEYQFYsISEBR44cwbZt26BQKKBQKFBeXo4jR45gzJgxUKlUUKvVSE1NxbVr10ye29DQgKVLl8Lb2xve3t547bXX0NTUZOzft28fIiIi4OHhgQEDBuCxxx5DVVVVV79EIrJjDCArtmbNGoSFhWHGjBkoKSlBSUkJ7rnnHkyZMgVBQUE4fPgwNm3ahI8++ggrV640eW5OTg6amppw8OBBvPPOO8jOzsbmzZsBADqdDs899xymT5+OoqIifP7555g2bZqIl0hEdsyujwFZu549e+Kee+6Bq6srVCoVAODNN9+ESqVCRkYGHBwc4Ofnh9dffx0vv/wyXn31Vbi6ugIAVCoV1q1bB5lMhoEDB+L06dPYvHkzFixYgIqKCly/fh2xsbHo378/ACAgIEDY6yQi+8QZ0F2mpKQEQ4cOhYPD///qHn74YVy7dg2lpaXGtoceesjkAtKwsDD88ccfuHTpEgYPHoxRo0YhIiICs2bNwvbt21FdXd2lr4OISFgApaWlGY9r3PgZOHCgsV+SJKSlpcHf3x8eHh6Ijo7GqVOnRJVrNSRJavHOBG29Y4GjoyP27t2LvLw8BAYGYteuXQgJCcHJkyctWSoR0R0JnQGp1WrjsY2SkhJ8++23xr4NGzYgMzMTa9euxaFDh6BUKhEXF4e6ujqBFXc9Z2dnNDY2Gh/7+/ujuLjY5ISCwsJCODs74/777ze2HTt2DJIkGR8XFxejX79+cHd3B/BnWIWFhWHZsmXIz89Hv379sHfv3i54RUREfxIaQE5OTlCpVMafPn36APjzf/kajQZJSUmIjY1FQEAANBoNDAYDcnNzRZbc5fr3749jx46hvLwcNTU1eO6551BZWYnFixejpKQEBw4cwMqVKzF37lzj8R8AqKysxLJly6DVavHJJ59g48aNmD9/PoA/w2j9+vU4fvw4zpw5g88//xznzp2Dn5+fqJdJRHZIaACVlZVh0KBBCAoKwrPPPouysjIAQHl5OXQ6HSIjI43Luri4ICIiAkVFRYKqFePFF1+Es7Mzhg0bBh8fH1y/fh05OTk4ceIERowYgQULFmDy5MlYsWKFyfOmTJmCpqYmjBkzBi+99BJmzZplDCB3d3cUFRVh6tSpCA0NxfLly5GcnIypU6eKeIlEZKdker1ean0xyzt48CAMBgPUajWqq6uxfv16aLVaHD16FFqtFuPHj8fJkyfh5eVlfE5iYiIqKiqQl5fX4nq1Wm2ztu7du0OpVHbK67AX58+fR319vegyiNpsaIFr6wvdQfEjVyxUif1Sq9V37Bd2GvbYsWNNHj/00EMYMmQI3n//fQwdOhRA84PqdzoAf8PtXnBtbW2nfWlcfX29TX8h3Q3u7u4m/xloiVarbXWnsyccD1NdOh4FHbvVVlfUae/7h9Wcht2jRw/4+/ujtLTUeM3LrVfmV1dXcyZDRGQjrCaA6uvrodVqoVKp4O3tDZVKhfz8fJP+wsJChIeHC6ySiIgsRdhHcMuXL8eECRPg6elpPAZ05coVTJ8+HTKZDAkJCcjIyIBarYavry/S09Ph5uaG+Ph4USUTEZEFCQugP/74A3PmzEFNTQ369OmDhx56CAcPHjTeGmbhwoW4evUqkpOTodfrERoairy8PMjlclElExGRBQkLoB07dtyxXyaTITU1FampqRbZXltOYKDbu/mCViIiS7GaY0Cdyc3NDXq9nm+k7SBJEvR6Pdzc3ESXQkQ2xi7uhu3k5AS5XI5Lly5ZfN2XLl0y3t7GVsnlcjg52cWuQkRdyG7eVZycnNCzZ0+Lr7eqqqpN18cQEZEpu/gIjoiIrA8DiIiIhGAAERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCMICIiEgIBhAREQnBACIiIiEYQEREJAQDiIiIhGAAERGREAwgIiISggFERERCWE0AZWRkQKFQIDk52dgmSRLS0tLg7+8PDw8PREdH49SpUwKrJCIiS7GKACouLkZ2djYCAwNN2jds2IDMzEysXbsWhw4dglKpRFxcHOrq6gRVSkREliI8gGprazF37lxs2rQJCoXC2C5JEjQaDZKSkhAbG4uAgABoNBoYDAbk5uaKK5iIiCxCeADdCJiRI0eatJeXl0On0yEyMtLY5uLigoiICBQVFXV1mUREZGFOIjeenZ2N0tJSbNmypVmfTqcDACiVSpN2pVKJioqKLqmPiIg6j7AA0mq1WLVqFb744gs4Ozu3uJxMJjN5LElSs7Zb19vVRGzTmnE8TLVnPIYWuHZom8WPXOnQ8ztT1+0fHRvDrqrTlv9e1Gr1HfuFBdB3332HmpoaPPzww8a2xsZGfPvtt9ixYweOHj0KAKiqqoKnp6dxmerq6mazopu19oItTavVdvk2rRnHw1S7x6PgXIe2a62/gy7dP+6CMbT3vxdhARQdHY3g4GCTtsTERPj4+GDRokXw9fWFSqVCfn4+QkJCAAD19fUoLCzEqlWrRJRMREQWJCyAFAqFyVlvAODq6opevXohICAAAJCQkICMjAyo1Wr4+voiPT0dbm5uiI+PF1AxERFZktCTEFqzcOFCXL16FcnJydDr9QgNDUVeXh7kcrno0ohsmuKfHfv4CgD0s++zQCVky6wqgD777DOTxzKZDKmpqUhNTRVUERERdRbh1wEREZF9YgAREZEQDCAiIhKCAUREREIwgIiISAgGEBERCWF2AOn1+k4og4iI7I3ZAeTn54enn34an332Ga5fv94ZNRERkR0wO4Cef/55HDt2DDNnzoSfnx+WLFmC4uLizqiNiIhsmNkB9Oabb+LHH3/Exx9/jPHjx2PPnj0YP348QkJCsG7dOpSVlXVCmUREZGvadRKCTCbDyJEjodFo8Msvv2Dbtm1Qq9VYv349QkJCEBUVhZ07d/J4ERERtajDZ8F1794dkydPxssvv4yoqChIkoSjR4/i5ZdfxqBBg5CcnIxLly5ZolYiIrIhHboZ6a+//ooPP/wQOTk5KC8vR9++fbFgwQJMnz4dzs7O2LlzJ7KysnDu3Dm8//77lqqZiIhsgNkBVFNTg48++gh79uzB8ePH4ezsjMceewzr1q3DmDFj4ODw/5Oq1atXQ6VSIS0tzaJFExHR3c/sAPL390dDQwPCwsLw9ttvIy4uDj179mxxebVajT59+nSoSCIisj1mB9DChQsxffp0+Pj4tGn5CRMmYMKECWYXRkREts3sAFq+fHln1EFERHbG7LPgdu3ahVmzZrXY//TTT/OEAyIiapXZAbR9+3aoVKoW+z08PJCVldWhooiIyPaZHUC//vorAgMDW+wfNGgQTp8+3aGiiIjI9pkdQDKZDDU1NS32X7hwAU1NTR0qioiIbJ/ZAfTggw8iJycH9fX1zfquXr2KnJwcBAUFWaQ4IiKyXWYH0KJFi6DVajF+/Hh88skn0Gq1OH36ND755BNERUVBq9Vi0aJFnVErERHZELNPwx49ejQ2b96MlJQUzJ4929guSRLkcjk2bdqERx991KJFEhGR7WnXveCmTZuG6OhoHDp0CGVlZZAkCffffz8iIyMhl8stXSMR3YUU/zx3m1ZXoOB27WSP2n0zUrlcjtjYWEvWQkREdqTdAVRXV4ezZ8/i4sWLkCSpWf/w4cM7VBgREdk2swNIr9cjJSUFe/fuRWNjI4A/j//IZDKTf1+4cMGylRIRkU0xO4CSkpKwf/9+zJ07F8OHD4dCoeiEsoiIyNaZHUBfffUV5s2bh7feeqsz6iEiIjth9nVAzs7Obf4qBiIiopaYHUCxsbE4ePBgZ9RCRER2xOwAevHFF1FZWYkXXngBxcXFqKysxPnz55v9tGbbtm2IiIiAl5cXvLy8MHbsWBw4cMDYL0kS0tLS4O/vDw8PD0RHR+PUqVPmlktERFbK7GNAoaGhkMlk+P7777Fnz54Wl2vtLLh7770XK1euhI+PD5qamvCvf/0LM2bMwNdff42//vWv2LBhAzIzM5GZmQm1Wo1169YhLi4OxcXFvNiViMgGmB1AKSkpxlOuOyI6Otrk8WuvvYbt27ejuLgYgYGB0Gg0SEpKMl7sqtFooFarkZuba3ILICIiujuZHUCpqakWL6KxsREff/wxLl++jLCwMJSXl0On0yEyMtK4jIuLCyIiIlBUVMQAIiKyAe2+EwLwZ3DU1tbC3d0dTk7mr+qnn37CuHHjUF9fDzc3N+zevRuBgYEoKioCACiVSpPllUolKioqOlIyERFZiXYF0PHjx7Fq1SoUFhbi+vXr2Lt3L0aOHImamhokJCQgMTERI0eObHU9arUa33zzDWpra7Fv3z4kJCRg//79xv5bP+q7+Y4LLdFqte15SR0iYpvWjONhqn3j4Spgm5bbvi3oqv3Ylv9e1Gr1HfvNDqDvvvsOMTExUKlUmDZtGt59911jX+/evWEwGLBr1642BZCzszMeeOABAEBwcDCOHz+OzZs3Y8mSJQCAqqoqeHp6Gpevrq5uNiu6VWsv2NK0Wm2Xb9OacTxMtXs8OnjH6A7/DnjH6i7Zj+3978Xs07DffPNN+Pj4oKioCCtWrGjWP2LECPz3v/9tVzFNTU24du0avL29oVKpkJ+fb+yrr69HYWEhwsPD27VuIiKyLmbPgI4fP47ly5eje/fuuHLlSrP+++67DzqdrtX1vPHGGxg3bhzuu+8+GAwG5ObmoqCgAHv27IFMJkNCQgIyMjKgVqvh6+uL9PR0uLm5IT4+3tySiYjICpkdQA4ODnBwaHnipNPp4OLi0up6dDodnn/+eVRVVcHd3R2BgYHIzc3FmDFjAAALFy7E1atXkZycDL1ej9DQUOTl5fEaICIiG2F2AA0ZMgRffvkl5s2b16zv2rVryMnJQVhYWKvr0Wg0d+yXyWRITU3tlNO+iYhIPLOPAS1atAiHDx/GggULcPLkSQBAZWUlvvrqK8TExOC3337D4sWLLV4oERHZFrNnQKNHj8aWLVuQnJyM999/HwCQkJAASZLQs2dPZGVlYejQoRYvlIiIbEu7rgOKj4/HY489hkOHDqG0tBRNTU24//77MWbMGPTo0cPSNRIRkQ1q950QXF1d8fjjj1uyFiIisiNmB9CZM2fatJyXl5fZxRARkf0wO4CCgoLadDfs1r6OgYiI7JvZAfT3v/+9WQA1NjaivLwcH3zwAfr27Ys5c+ZYrEAiIrJNZgfQjBkzWuxLSkpCZGQkDAZDh4oiIiLbZ/Z1QHfSo0cPzJgxA5s3b7bkaomIyAZZNIAA4J577uF39hARUassGkAnT57EP/7xD/j5+VlytUREZIMsdhZcbW0tLl26hB49eiAzM9MixRERke0yO4CGDx/eLIBkMhkUCgUeeOABTJ48GQqFwlL1ERGRjTI7gFq7izUREVFbWPwkBCIiorYwewa0du1aszcik8mQkpJi9vOIiMh2mR1Aa9asMR4DkiTJpO9O7QwgIiK6mdkB9PPPP2PKlCnw9/dHQkICfH19AQBarRYajQYlJSXIycmBSqWyeLFERGQ7zD4GlJKSAm9vb2zduhXBwcGQy+WQy+UICQnBtm3b0L9/f852iIioVWYHUH5+PkaNGtVi/+jRo5Gfn9+RmoiIyA6YHUCOjo748ccfW+w/ceIEHBx4ch0REd2Z2UkRExODd999FxkZGairqzO219XVIT09Hbt370ZsbKxFiyQiIttj9kkIq1evxm+//YbVq1cjLS0Nffv2hUwmg06nQ2NjI4YPH47Vq1d3Rq1ERGRDzA4guVyOffv24YsvvsCBAwdw9uxZSJKEcePGYdy4cYiKiuqMOomIyMaYHUA3REVFMWyIiKjd2h1AZ86cwZEjR3D+/HnExcXB09MTDQ0NuHjxInr16gUnp3avmoiI7EC7UuKVV17B1q1b0djYCJlMhqCgIHh6euLKlSsICQnBsmXLkJiYaOlaiYjIhph9FtzGjRuh0WiQmJiIjz/+2OS2O+7u7oiOjsb+/fstWiQREdkes2dA2dnZePLJJ7Fy5UpcuHChWX9gYCAOHTpkkeKI7kaKf5676ZErUHCuxWWJ7JnZM6CzZ88iIiKixX65XI7a2toOFUVERLbP7AD6y1/+gsrKyhb7f/rpJ/Tr169DRRERke0zO4DGjRuH7Oxs1NTUNOv74YcfsHv3bkRHR1ukOCIisl1mB9Arr7wCBwcHRERE4I033oBMJsN7772HZ599FmPHjsW9996L5OTkzqiViIhsiNkBpFKp8PXXX2PChAn49NNPIUkScnJy8NVXX2Hq1Kn497//DYVC0ep63n77bYwePRpeXl7w8fHB1KlT8b///c9kGUmSkJaWBn9/f3h4eCA6OhqnTp0yt2QiIrJCZgXQ9evXceTIEdTW1mLDhg347bffoNVqUVJSgrKyMmzatAm9e/du07oKCgrw3HPP4cCBA9i3bx+cnJwwceJEXLx40bjMhg0bkJmZibVr1+LQoUNQKpWIi4szuQkqERHdncwKIEdHR0ycONHkNOs+ffqgb9++Zn8FQ15eHmbOnImAgAAEBgZiy5YtqK6uxtGjRwH8OfvRaDRISkpCbGwsAgICoNFoYDAYkJuba9a2iIjI+piVGg4ODujfvz8MBoPFCzEYDGhqajJ+fFdeXg6dTofIyEjjMi4uLoiIiEBRUZHFt09ERF3L7AtR58+fj40bN2LmzJlQKpUWK2TZsmUYPHgwwsLCAAA6nQ4Amm1DqVSioqKixfVotVqL1dRWIrZpzTgerqILsMDvQPxrEK2r9mNb/ntRq9V37Dc7gAwGA9zc3BASEoLo6GgMGDAALi4uJsvIZDK89NJLbV7nK6+8gqNHj+LLL7+Eo6Njs3XdTJKkZm03a+0FW5pWq+3ybVozjges4s4HHf4dWMFrEK0r9mN7/3sxO4DeeOMN478//PDD2y5jTgClpqYiLy8Pn376KQYMGGBsV6lUAICqqip4enoa26urqy068yIiIjHaFECLFy/GzJkzERwcjB9++AGSJOHixYtwd3fv0NcuLF26FHl5edi/fz8GDhxo0uft7Q2VSoX8/HyEhIQAAOrr61FYWIhVq1a1e5tERGQd2pQeO3bswLBhwxAcHIz+/fvjwoULCA4Oxt69ezFy5Mh2bXjJkiX48MMPsXv3bigUCuMxHzc3N/To0QMymQwJCQnIyMiAWq2Gr68v0tPT4ebmhvj4+HZtk6gtTG8mSkSdpd3Tl5u/hqE9srKyAACxsbEm7UuXLkVqaioAYOHChbh69SqSk5Oh1+sRGhqKvLw8yOXyDm2biIjEE/a1pXq9vtVlZDIZUlNTjYFERES2w+xb8RAREVlCm2dAZWVlOHbsGADg0qVLAP48hbBHjx63XT40NNQC5RERka2S6fX6Vg/m9OrVq83X49xov923pdoiez+P/1a2MB48CYEsRT/7vjv228LfS0e0aQaUmZnZ2XUQEZGdaVMAPfXUU51dBxER2RmehEBEREIwgIiISAgGEBERCcEAIiIiIRhAREQkBAOIiIiEYAAREZEQDCAiIhKCAUREREIwgIiISAhh3wdE1JKO3gy0tRtAEpF14AyIiIiEYAAREZEQDCAiIhKCAUREREIwgIiISAgGEBERCcEAIiIiIRhAREQkBAOIiIiEYAAREZEQDCAiIhKC94Ijixpa4AoUdOxebh3V0XvJEVHX4AyIiIiEYAAREZEQDCAiIhKCAUREREIIPQnhyJEj2LRpE3744QdUVFQgMzMTM2bMMPZLkoQ1a9YgOzsber0eoaGhSE9Px6BBgwRWTUTUNq2fEHPnk3Zs/csVhc6ALl++jICAAKxZswYuLi7N+jds2IDMzEysXbsWhw4dglKpRFxcHOrq6gRUS0REliQ0gMaNG4cVK1YgNjYWDg6mpUiSBI1Gg6SkJMTGxiIgIAAajQYGgwG5ubmCKiYiIkux2mNA5eXl0Ol0iIyMNLa5uLggIiICRUVFAisjIiJLsNoLUXU6HQBAqVSatCuVSlRUVLT4PK1W26l1Wcs2rZer6AKIbMbd/t6iVqvv2G+1AXSDTCYzeSxJUrO2m7X2gi1Nq9V2+TatmuC7IBDZElt/b7Haj+BUKhUAoKqqyqS9urq62ayIiIjuPlYbQN7e3lCpVMjPzze21dfXo7CwEOHh4QIrIyIiSxD6EZzBYEBpaSkAoKmpCWfPnsWJEyfQq1cveHl5ISEhARkZGVCr1fD19UV6ejrc3NwQHx8vsmwiIrIAmV6vl0Rt/JtvvsETTzzRrH369OnQaDTGC1F37txpciFqQECAgGpvj8eATPFO1ESWY+sXogoNIFvAADLFACKyHFsPIKs9BkRERLaNAUREREIwgIiISAgGEBERCcEAIiIiIRhAREQkBAOIiIiEYAAREZEQDCAiIhKCAUREREJY/fcB2RNL3MbG1m/dQUS2gzMgIiISggFERERCMICIiEgIBhAREQnBkxBu0r6TAFyBgj+fxxMAiIjajjMgIiISggFERERCMICIiEgIHgOyMR29mJXHsYishyUuTu+oznxP4AyIiIiEYAAREZEQDCAiIhKCx4DIhDV85kxE9oEzICIiEoIBREREQjCAiIhICAYQEREJwQAiIiIhGEBERCQEA4iIiIRgABERkRB3xYWoWVlZ2LhxI3Q6Hfz9/ZGWloaIiAjRZTXDiziJiNrO6mdAeXl5WLZsGRYvXozDhw8jLCwMU6ZMwZkzZ0SXRkREHWD1AZSZmYmnnnoKzzzzDPz8/LB+/XqoVCrs2LFDdGlERNQBVv0R3LVr1/D999/jxRdfNGmPjIxEUVGRxbfH78IhIuo6Vj0DqqmpQWNjI5RKpUm7UqlEVVWVoKqIiMgSrDqAbpDJZCaPJUlq1kZERHcXqw6g3r17w9HRsdlsp7q6utmsiIiI7i5WHUDOzs4YMmQI8vPzTdrz8/MRHh4uqCoiIrIEqz4JAQASExMxb948hIaGIjw8HDt27EBlZSVmz54tujQiIuoAq54BAcCkSZOQlpaG9evXY8SIETh69Cj27NmD/v37d1kNlZWVeOGFF+Dj4wOVSoXw8HAUFBQY+yVJQlpaGvz9/eHh4YHo6GicOnWqy+rraq2NR0JCAhQKhcnPo48+KrDizjN48OBmr1WhUODJJ58EYH/7RmvjYU/7BgA0NjZi9erVCAoKgkqlQlBQEFavXo2GhgbjMva2j9zM6mdAADBnzhzMmTNHyLb1ej3Gjx+PYcOGYc+ePejduzfKy8tNjkFt2LABmZmZyMzMhFqtxrp16xAXF4fi4mLI5XIhdXeWtowHAIwaNQpbtmwxPnZ2du7qUrtEfn4+GhsbjY8rKysxatQoTJw4EYB97RtA6+MB2M++AQDvvPMOsrKyoNFoEBAQgJ9++gkJCQlwdnZGSkoKAPvbR252VwSQSBs3boSHh4fJH8yAAQOM/5YkCRqNBklJSYiNjQUAaDQaqNVq5Obm2txHha2Nxw3dunWDSqXqwsrE6NOnj8njXbt2QS6XY+LEiXa3bwB3Ho8b7GXfAIDvvvsOEyZMQFRUFADA29sbUVFROHbsGAD7e/+4ldV/BCfaZ599htDQUMyePRu+vr545JFHsHXrVkiSBAAoLy+HTqdDZGSk8TkuLi6IiIjolItlRWttPG4oLCyEr68vQkND8dJLL+H8+fOCKu46kiRh165dmDp1KlxdXe1u37jVreNxgz3tG8OGDUNBQQF++eUXAMDPP/+Mb775BmPHjgVgf+8ft+IMqBVlZWXYvn075s+fj6SkJJw8eRJLly4FADz//PPQ6XQAcNuLZSsqKrq83s7W2ngAwKOPPoonnngC3t7e+P3337F69WrExMTg66+/Rrdu3USW36ny8/NRXl6OWbNmAYDd7Ru3unU8APvbN5KSkmAwGBAeHg5HR0c0NDRgyZIlxkMK9r6PMIBa0dTUhODgYLz++usAgAcffBClpaXIysoyvuEC9nOxbFvGY/LkycblAwMDMWTIEAwePBgHDhxATEyMkLq7QnZ2NkJCQhAUFGTSbi/7xq1uNx72tm/k5eXhgw8+QFZWFvz9/XHy5EksW7YM/fv3x9NPP21czl73EX4E1wqVSgU/Pz+TtoEDB+Ls2bPGfgB2c7Fsa+NxO/369cO9996L0tLSzi5PmPPnz+Pzzz/HM888Y2yzt33jZrcbj9ux9X1jxYoVWLBgASZPnozAwEBMmzYNiYmJ+Nvf/gbAvvcRgAHUqmHDhuH06dMmbadPn4aXlxeAPw8qqlQqk4tl6+vrUVhYaJMXy7Y2HrdTU1ODiooKmz7w/N5776Fbt26YNGmSsc3e9o2b3W48bsfW940rV67A0dHRpM3R0RFNTU0A7HsfAfgRXKvmz5+PcePGIT09HZMmTcKJEyewdetWvPbaawD+nDonJCQgIyMDarUavr6+SE9Ph5ubG+Lj4wVXb3mtjYfBYMCaNWsQExMDlUqF33//HatWrYJSqcTjjz8uuPrOIUkS3n33XUyaNMnktFl72zduaGk87HHfmDBhAt555x14e3vD398fJ06cQGZmJqZNmwbAfveRG2R6vV5qfTH7duDAAaxatQqnT5+Gp6cn5s6di3nz5hk/o5UkCWvWrMHOnTuh1+sRGhqK9PR0BAQECK68c9xpPK5evYoZM2bgxIkTqK2thUqlwogRI/Dqq6/C09NTdOmd4vDhw4iJicF//vMfhIaGmvTZ274BtDwe9rhv1NXV4a233sL+/ftRXV0NlUqFyZMnIyUlBd27dwdgn/vIDQwgIiISgseAiIhICAYQEREJwQAiIiIhGEBERCQEA4iIiIRgABERkRAMICIiEoIBREREQjCAiIhIiP8DqnUvgTxsTuQAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df.plot.hist(bins=20)\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[71.0, 77.21989528795811, 83.0]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Write a function called `calc_temps` that will accept start date and end date in the format '%Y-%m-%d' \n",
    "# and return the minimum, average, and maximum temperatures for that range of dates\n",
    "def calc_temps(start=None,end=None):\n",
    "    sel = [func.min(Measurement.tobs), func.avg(Measurement.tobs), func.max(Measurement.tobs)]\n",
    "\n",
    "    if not end:\n",
    "        results = session.query(*sel\n",
    "        ).filter(Measurement.date >= start).all()\n",
    "        temps = list(np.ravel(results))\n",
    "        return temps\n",
    "\n",
    "    results = session.query(*sel).filter(Measurement.date >= start\n",
    "    ).filter(Measurement.date <= end).all()\n",
    "    temps = list(np.ravel(results))\n",
    "    return temps\n",
    "\n",
    "calc_temps(\"2017-06-01\", \"2017-06-31\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sqlalchemy import extract\n",
    "\n",
    "# sel = [func.min(Measurement.tobs), func.avg(Measurement.tobs), func.max(Measurement.tobs)]\n",
    "start= \"2011-06-01\"\n",
    "end = \"2011-06-31\"\n",
    "juneeee = session.query(func.min(Measurement.tobs), func.avg(Measurement.tobs), func.max(Measurement.tobs)).filter(Measurement.date >= start).filter(Measurement.date <= end).all()\n",
    "\n",
    "\n",
    "print(juneeee)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Challenge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}