{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# IS608 Project -- Sreejaya Vasudevannair\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Project Title: US Motor Vehicle Fatality Analysis\n", "The goal of this project is to create choropleth map of US Motor Vehicle Fatality on US State level. \n", "The National Highway Traffic Safety Administration(NHTSA) has a dataset called **[Fatality Reporting System(FARS)](ftp://ftp.nhtsa.dot.gov/fars) ** to record every fatal traffic crash in the US. FARS is a nationwide census providing NHTSA, Congress and the American public yearly data regarding fatal injuries suffered in motor vehicle traffic crashes. \n", "Collected data include \n", "1.\tCrash characteristics\n", "2.\tVehicle characteristics\n", "3.\tPerson characteristics\n", "4.\tWeather and time of the day.\n", "\n", "### Data Source #### \n", "**[Download Raw Data from FTP Site](ftp://ftp.nhtsa.dot.gov/fars) ** \n", "and http://www-fars.nhtsa.dot.gov/Trends/TrendsGeneral.aspx.\n", " \n", "The files are saved as .dbf. I converted them to comma delimited vectors(CSV). \n", "\n", "The project provide a general idea on the Fatal Traffic crashes in US. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "%matplotlib inline\n", "import numpy as np\n", "from numpy import *\n", "import random\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "\n", "import plotly\n", "from plotly.graph_objs import Scatter, Layout\n", "import plotly.plotly as py\n", "import plotly.graph_objs as go" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "** Fatality trend from 1994 to 2014 **" ] }, { "cell_type": "code", "execution_count": 55, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "#df=pd.read_csv('C:/DataAnalytics/MSDASemester3/608/Project/Data_proj/CSV/USAFatalCrashes1994_2014.csv')\n", "df=pd.read_csv('https://raw.githubusercontent.com/srajeev1/IS608Project/master/data/USAFatalCrashes1994_2014.csv')\n", "plotly.offline.init_notebook_mode() # run at the start of every notebook\n", "plotly.offline.iplot({\n", "\"data\": [{\n", " \"x\": df.Year,\n", " \"y\": df.Total\n", "}],\n", " \n", "\"layout\": go.Layout(\n", " title='Total Crashes By Year',\n", " xaxis=dict(\n", " title='Year',\n", " titlefont=dict(\n", " family='Courier New, monospace',\n", " size=18,\n", " color='#7f7f7f'\n", " )\n", " ),\n", " yaxis=dict(\n", " title='Count',\n", " titlefont=dict(\n", " family='Courier New, monospace',\n", " size=18,\n", " color='#7f7f7f'\n", " )\n", " )\n", ")\n", "})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The Fatality Rate has declining from 2005." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Top 5 US States based of Traffic fatality Rates.###" ] }, { "cell_type": "code", "execution_count": 57, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Anaconda2\\lib\\site-packages\\ipykernel\\__main__.py:8: FutureWarning:\n", "\n", "sort(columns=....) is deprecated, use sort_values(by=.....)\n", "\n" ] }, { "data": { "text/html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
CaliforniaTexasFloridaGeorgiaNorth Carolina
199442323187268714251431
199541923183280514881448
199639893742275315731494
199736883513278515771483
199834943586282515681596
\n", "
" ], "text/plain": [ " California Texas Florida Georgia North Carolina\n", "1994 4232 3187 2687 1425 1431\n", "1995 4192 3183 2805 1488 1448\n", "1996 3989 3742 2753 1573 1494\n", "1997 3688 3513 2785 1577 1483\n", "1998 3494 3586 2825 1568 1596" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" } ], "source": [ "FatalitybySTATE_Data=pd.read_csv('https://raw.githubusercontent.com/srajeev1/IS608Project/master/data/USAFatalities_FatalityRatesbySTATE1994_2014.csv')\n", "FatalitybySTATE_Data =FatalitybySTATE_Data[pd.notnull(FatalitybySTATE_Data['State'])]\n", "\n", "#find the total fatality for 1994-2014\n", "FatalitybySTATE_df =FatalitybySTATE_Data[FatalitybySTATE_Data.columns[0:23]]\n", "FatalitybySTATE_df =FatalitybySTATE_df.drop(FatalitybySTATE_df.columns[1],axis=1)\n", "FatalitybySTATE_df['Total']=FatalitybySTATE_df.sum(axis=1)\n", "FatalitybySTATE_df=FatalitybySTATE_df.sort(['Total'],ascending=[0])\n", "\n", "#top 5 States\n", "FatalitybySTATE_df[0:1].drop(FatalitybySTATE_df.columns[0],axis=1)\n", "FatalitybySTATE_dt=FatalitybySTATE_df[1:6].transpose()\n", "FatalitybySTATE_dt=FatalitybySTATE_dt[1:22]\n", "FatalitybySTATE_dt.columns=['California', 'Texas', 'Florida', 'Georgia', 'North Carolina']\n", "FatalitybySTATE_dt.head(5)" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmQAAAEACAYAAAAdnGptAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdcleX/x/HXhRNxgwPcAzN3uUe5ylGucuvP1LZ9czbc\nSZqW5izThpaaqaRljhJLTa00yYHm3htRAZG9zvX74wKExMmB+wCf5+NxHpxzn3WBt/A+1/hcSmuN\nEEIIIYSwjpPVDRBCCCGEyO4kkAkhhBBCWEwCmRBCCCGExSSQCSGEEEJYTAKZEEIIIYTFJJAJIYQQ\nQljsvgKZUiqHUmqfUmpdwm0vpdTFhGP7lFLtkz12tFLqhFLqqFKqTbLjdZVS/ybcN8f+34oQQggh\nROZ0vz1kQ4HDQGLRMg3M1Fo/lnDZAKCUqgb0BKoB7YB5SimV8Jz5wEtaa0/AUynVzl7fhBBCCCFE\nZnbPQKaUKg08AywAEsOVSnY9uc7Acq11rNb6LHASaKiUcgcKaK19Ex63BOiSxrYLIYQQQmQJ99ND\nNgt4B7AlO6aBwUqp/UqphUqpwgnHPYCLyR53ESiVyvFLCceFEEIIIbK9uwYypVQH4KrWeh8pe8Tm\nAxWAOoA/MCPdWiiEEEIIkcXlvMf9TYBOSqlngLxAQaXUEq31C4kPUEotANYl3LwElEn2/NKYnrFL\nCdeTH7+U2hsqpWRzTSGEEEJkGlrr1KZxPRB1v5uLK6WaA29rrTsqpdy11v4Jx4cD9bXWfRIm9S8D\nGmCGJDcBlbXWWim1CxgC+AI/A59orX1SeR8tG56L++Hl5YWXl5fVzRCZhJwv4n7JuSIehFLKLoHs\nXj1kKd6TW6sspymlaifcPgO8BqC1PqyU+h6zIjMOeCNZunoDWAQ4A7+kFsaEEEIIIbKj+w5kWuut\nwNaE6/3u8rgpwJRUju8Baj5wC4UQQgghsjip1C8yrRYtWljdBJGJyPki7pecK8IK9z2HLKPIHDIh\nhBBCZBb2mkMmPWRCCCGEEBaTQCaEEEIIYTEJZEIIIYQQFpNAJoQQdhQfD/v2QViY1S0RQmQmD1KH\nTAghxB3ExsLSpfDhh2Czgb8/VKoEDRtCgwbma7VqkFN+6wohUiGrLIUQIg2io2HRIvjoIxPAxo+H\n5s0hJgYOHABfX9i1y3y9eBEeeyxlSCtTBlSa12cJIaxir1WWEsiEEOIhRETAggXw8cdQqxaMGweN\nG9/9OTduwO7dt0Larl3meGI4a9AA6teHwoXTv/1CCPuQQCaEEBYIC4P582HmTGjUyASxunUf7rW0\nhgsXUvai7d0LpUqlDGm1a0Pu3Pb9PoQQ9iGBTAghMlBICHz6KXzyCbRqBWPHQs102AwuLg4OH04Z\n0k6eNO+VPKRVrixDnUI4AglkQgiRAQIDYfZs0yv27LMwZgw88kjGtiE8HPbsSRnSQkNNMGvQAFq3\nhieflIAmhBWydCDr109TsiSUKGEuiddLlgRXV3CSYh1CZFpRUSZQVKpkhuYcVUAAzJgBCxdC164w\nahRUrGh1q265csX8HH19YeVKKFvWtLdWLatbJkT2kqUD2TffaAICzC+cgABSXA8JATe3lCHtv8Et\n8XrRohLehHAEwcHw88+wZg389psZbjt7FgoWhCeeuHWpUsX6Xp5Ll8xE/SVLoG9feOcdE3YcWWws\nfPklTJwIHTvCpEng7m51q4TIHrJ0ILtbm2Jj4erVlCEtteB25YqZfFus2J0DW8WKprvf6j8AQmRF\n58+bAPbTT/DPP9CyJXTpAh06mP+XNhscPQp//HHrEhkJzZrdCmh16mRc3a6zZ2HqVPD2hoED4e23\nM1+ouXEDpkwxvXrDhsFbb0G+fFa3SoisLdsGsgcRE2PC250C2/795pfVu+9Ct25SsFGItNDa1N36\n6ScTxM6fN701nTvD00+Di8u9X+P8efjzz1sB7fx5s5LxiSdMUGvY0P4B4+RJE2LWrIHXXoPhw01g\nzMzOnIHRo+Gvv2DyZPi//5PRAiHSiwQyO7DZ4JdfzKfiy5fNp8mBA8HZOUPeXohMLy7OBKjEEKaU\n6QXr0gWaNEn7h5zAQBMq/vjDvM+BA2aOVGIPWtOmZmrCwzh82ASxjRvhzTdh8OCHfy1HtXMnjBhh\nPpzOmAEtWljdIiGyHglkdrZjB0ybBn//bX45v/FG1vvlLIQ9hIebELNmjZkXVr686QXr0gVq1Ejf\nKQAREWaVYWIP2q5dZn5X8nloZcrc/TX8/Eyv0fbtpjfsjTfMXLasSmv4/nuzKKFWLfN7LqNXiQqR\nlUkgSydHjpgJvT/9BP37m0+X9/oFL0RWd/UqrFtnQtjWrWbosHNn6NTJ2gnvcXEmYCWfh+biciuc\nNWsGjz5qQqKvL3zwgamU//bbZnjyfoZRs4qoKFNHbepU6NMHJkwwq9aFEGkjgSydXbwIc+bA11+b\nScjvvGM+/QuRXZw4cWtS/sGD0KaN6QV75hnH3dpHazh27FY4+/NPuHnT9OJdvQojR8JLL0HevFa3\n1DrXrsH775vFC6NGmRGBPHmsbpUQmZcEsgxy44YpCPnJJ2Z7lJEjzaduWZkpshqbzfQeJYawoCDT\nA9ali6lMn1n/aF+8aHq+mzeX7YeSO3LELGg6fNj0mnXtKr/XhHgYEsgyWFSUqUv08cdmBdbIkWYF\nmaxcEplZVBRs2WKGI9euNXOpunQxw5ENGsj5nR1s3mwWNOXPbyb+N2xodYuEyFwkkFkkPh5Wrzaf\nKMPCzFBm376Zt/dAZD8BAWYy/rp1JozVqnWrPIVM9s6e4uPNB85x48wWTB9+aIZ5hRD3Zq9Adl+f\nf5VSOZRS+5RS6xJuF1VK/aaUOq6U+lUpVTjZY0crpU4opY4qpdokO15XKfVvwn1z0tpwq+TIYWqW\n+frCZ5+Z1UuVKsH06WauihCORmv4919T4qFxYxO6fHzg+efh1Ckz1+rddyWMZWc5cpiSP8eOmfOg\nbl0zvywkxOqWCZF93O+AxFDgMJDYdTUK+E1rXQXYnHAbpVQ1oCdQDWgHzFMqaVbCfOAlrbUn4KmU\namefb8EaSpl5NT4+sH497N1rKv+PHm2KzgphpZgYs0XRkCHmvOzY0ZyXEyeaye3ffw/9+pltyIRI\nlD8/eHmZem9Xr5pwNn++Wc0qhEhf9wxkSqnSwDPAAiAxXHUCFidcXwx0SbjeGViutY7VWp8FTgIN\nlVLuQAGttW/C45Yke06mV6cOLFtmtocJDYVq1cyS+hMnrG6ZyE4CA+Hbb6F7dyheHN57z2wRtnat\nqdz+ySemYr5MbBf3UqqUWWHu4wOrVplh7V9+Mb2tQoj0cT89ZLOAdwBbsmMltNYBCdcDgBIJ1z2A\ni8kedxEolcrxSwnHs5QKFWDuXNPtX7KkqVTerZsJakLYm9ZmL8hp00zNrYoV4ccfoX17cw7u3Alj\nxkDNmrJ6TjycOnVg0yYzZ3bECFP6ZP9+q1slRNZ0141NlFIdgKta631KqRapPUZrrZVSdv3c5OXl\nlXS9RYsWtMhk+30UK2bq/Lzzjtnkt1s3M89s7Fho3drq1onMLDbW1NZat85cIiPNcOTo0WYIPTvX\n1xLpQylzjrVrB19+aULZwIHw0UdWt0wIa2zdupWtW7fa/XXvuspSKTUF6AfEAXmBgsCPQH2ghdb6\nSsJw5O9a66pKqVEAWuuPEp7vA0wAziU85tGE472B5lrr11N5T4deZfkwYmNNEcbx4+Gpp2DmTChQ\nwOpWicwiONgMHa1bZ74mzgnr1Mn0YEjvl8hIwcFmD9ERI+Dll61ujRDWy/CyF0qp5sDbWuuOSqlp\nQKDWempCCCustR6VMKl/GdAAMyS5Caic0Iu2CxgC+AI/A59orX1SeZ8sF8gS3bxpfolt2QKLFpnl\n5UKk5tw5M/y4bp0p1vrkkyaEdehg5vcIYaWjR80wuY+PWZEpRHZmVSB7S2vdSSlVFPgeKAucBXpo\nrW8kPG4M8CKmV22o1npjwvG6wCLAGfhFaz3kDu+TZQNZonXrzKT/vn1h0iQZZhK37Ntn5oT99psp\n0Nqpkxnmzk57LorMYdUqMy1j927ZE1Nkb1IYNpO7dg1ef91Mvv72W3jsMatbJKyitek1nTrVbGMz\nbBi8+qqpmi+EI3vnHVPj7uefTS0zIbIjCWRZgNbw3XdmGHPYMFOcM+ddl1mIrCQ+3gxLTp0K4eHm\n379PH9n1QWQecXGmlMoTT5gad0JkRxLIspALF8yqpfBws32Jp6fVLRLpKTISFi82uzsULy77oorM\nLSAA6tUzBWQ7dLC6NUJkvAzdOkmkrzJl4NdfTe9I48Ywb54UYMyKgoPN9kUVK5ohnkWLYMcOs4ek\nhDGRWZUoYVaRv/ii2YpLCPFwpIfMwRw9Ci+8AEWKmErZsqIu87t4EWbNgm++MZP0334batSwulVC\n2NfcubBggfmQkS+f1a0RIuNID1kWVbWq+YXWtKmZ6L98ufSWZVaHD5uh6Fq1zL/h/v2mV0zCmMiK\n/vc/qF4dBg2S31lCPAzpIXNgu3eb3rKaNc0wpiwtzxz++stM1N+1CwYPhjfegKJFrW6VEOkvPBwa\nNTLh7PXbyn4LkTXJpP5sIjLSbLnk7Q1ffQXPPGN1i0RqbDYzL2zqVPD3N8OSAwaAs7PVLbuzy6GX\n2XN5D3v89+B3xY+I2AiUUigUTsop6Xryr07K6b6O3ek1Eh+bK0cuGpRqQNtKbXEv4G71j0LY0YkT\npod/3Tpo2NDq1giR/iSQZTNbt5o/8G3bwowZkD+/1S0SADExsGwZfPyxKVcxciR07epY5Uu01lwK\nvZQUvvb472Gv/15i42Op61GXuu51eazkYxTMUxCNRmuNRmPTtqTryb/atO2hjyW+bmRsJH9e+JNN\npzdRtlBZ2lZqS7vK7Whapil5ckrdj8xuzRrTO7xnj9nbV4isTAJZNnTzpqlXtm2bKY/RtKnVLcq+\nQkNNj+WsWWbe38iRpqK+1ftKaq25ePOiCV7JAphN26jrbsJXYggrW6gsyuIGx9ni8L3ky8aTG/E5\n5cORa0d4styTtKvcjraV2lK5aGXL2ygezpgx4OsLGzdK0ViRtUkgy8bWrDHzM/r3h/ffl0Ki6c3v\nih9v/foWZ4LP4JbHnZv+7pz71x1Pd3eef9qdxjU8cC/gjnt+d1zzueKkMmatjNaa8yHnk3q8EkOY\nUuq28FW6YOlMEWwCIwLZdHoTPqd82HhyI865nGlbqS1tK7WlVYVWFMhTwOomivsUHw/t2kH9+qbc\nixBZlQSybO7qVbMf5qlTZuul2rWtblHWExgRyPjfx/PDkR94s/pEjvzcmnVb/Wn8tD8Nn/InJvdl\n/MP8zSXUfA2NDqVE/hK453fHo4AH7vndk8Ja4lePAh4UdylODqf77zbQWnMu5Nxtw445VI6k0JUY\nwEoVKJUpwte9aK359+q/bDy5kY2nNrLr0i7qutdNGt6sXbJ2hoVf8XCuXTNFY+fMMXuzCpEVSSDL\nZALCAliyfwkvPf4SRZ3ts+ROazN0+fbb8NZbZl85GRpIu3hbPF/u+ZIJWyfQrEgPojZM5J/tRXnt\nNRgyxFTXv5OouCiuhF1JCmjJv14Ou5x0OygyCLd8bikDW2KIK+BOCZcSKeZ97fXfS+4cuVOEr8fd\nH8ejgEeWCF/3IzwmnK1nt7Lx1EZ8TvoQEh1Cm0ptaFepHU9XepriLnf5hxGW8fU1Ffz/+kt2IRFZ\nkwSyTCQ2Ppanvn2K2PhYTgSdYFjDYQxtNJT8ue0zM//cOVPvKjraBLRKlezystnS9nPbGfzLEGJu\nFgafT7BdrsXw4ab8iD2LXcbGx3I1/Cr+Yf5cDr18W4C7EnaFkvlLphh2lNWIKZ0JPpMUzn4/+zue\nRT3N8GbltjQu3ZhcOXJZ3USR4IsvTOHYv/8GFxerWyOEfUkgy0Te/vVtDl07xPre6zkdfBqvbV5s\nPr2Z0c1G81q918ibM2+a38Nmg08/hQ8+MJdXX7V+gnlmcvHmRYb9/C6bjv+J0+aPqZOrB2+NULRv\nL9saZQYx8THsvLCTjafM8OapoFO0rNCSdpXa0bZyW8oXLm91E7M1rc2HxthYWLpUfjeJrEUCWSbx\n/aHvGbVpFLtf3Z1iqPJAwAHGbRmH3xU/JjSfQP86/cnplPZaCUeOQL9+Zqn5woXg4ZHml8zSouKi\nGLNuJvP8ZqD+eYPnS4zi3WEuMicvkwsIC+C307+ZgHZyI675XOlTow8v1H6BcoXLWd28bCkyEpo0\nMXteDh5sdWuEsB8JZJnA4WuHab6oOb/+36885v5Yqo/ZeWEnY7aM4XLoZSa2mEj36t3TPFE5NhYm\nTzYTaXv2NKUyqlZN00tmOTab5sMf1zNlz3DiLtfglbIzGftGRdxlVDDLsWkb/1z6hyX7l+B9yJta\nJWrRv3Z/ulbrardpA+L+nD4NjRvD6tUmnAmRFUggc3A3o2/S4KsGjGo2igF1Btz1sVprNp/ZzJjN\nY4iJj2Fyq8k84/lMmidrBwTA/PnmUq8eDB/uGLWyrBQTA7O/Pc4Hu4cRne80r5efw4cvtZXNkLOJ\n6Lho1h1fx+L9i/nj3B90qdqF/rX707x8c1mxmUF+/tmsEN+zB0qUsLo1QqSdBDIHprWm28puFMtX\njM87fP5Az1tzbA3jtoyjUN5CTGk1heblm6e5PVFR8N13poipk5MJZn36ZK/6ZYGBMGd+KDP3fEBM\n9YX0rziaT/sNJm+u3FY3TVgkICyA7/79jsX7FxMSFcILtV+gf+3+VCoqq2LS24QJpsD1pk2OtauF\nEA9DApkDm/bXNH448gPbB2x/qG1g4m3xLD+4nAlbJ1C5aGUmt5pMPY96aW6X1vDbbyaY+fnBoEHm\nkpW3Njl6FGbNtrF0/3fw1CiertiGz7t/SMn8Ja1umnAgflf8WOS3iOUHl1PFtQr9a/enR/UeFMxT\n0OqmZUnx8aYURo0aZtsxITIzCWQOasuZLfT9sS++L/tSplCZNL1WTHwMX+/7mknbJ9GodCMmtZxE\ntWLV7NLOQ4dg9mxYtQq6dTO9ZtXs89KW0xq2bDHBc8fZPbh0G4xr8Vg+7/QpjUo3srp5woHFxsey\n4eQGFvktYsuZLTxb5VkG1B5AqwqtHqiQr7i3wEAzleLjj83vICEyKwlkDuhCyAUaLGjAd89/R6sK\nrez2upGxkXz2z2d8vONj2lVuh1dzLyoUqWCX1756FT7/3Mwzq10bRoyAp5/OnPPMoqNh+XITxKJy\nXKNE77EcU2uZ0moyAx8bKHOExAO5HnGd5f8uZ/H+xQSEB/B/Nf+P/nX6U9VNVsjYy5490L49bN8u\nC49E5iWBzMFEx0Xz5KIn6fpoV95t+m66vMfN6JvM2jmLT30/pWf1nox7cpzdioUmhpmZM01Ns+HD\noW9fyHuXEmlRcVEcvX6Uw9cOc+jqIcJjw/Es6omnqydVXKtQpmCZDOlVuHbNhMp586Bm7Tgq957P\nyqsT6VuzL14tvCict3C6t0FkbQevHmSx32KW/ruUcoXKMaDOAHpW70kR5yJWNy3TW7gQZsyAXbug\ngGxVKjIhCWQOZtD6QQSEB/BDjx/SfSub6xHXmfrnVL72+5qXH3uZd5u+i2s+V7u8ttawebPpZdq9\n28wxG/hKFDdyHOPQtUMcunqIQ9cOcfjaYS7cvEDFIhWpXqw61YtVxyW3CycCT3Ai6ATHA48TGBlI\nxSIV8SxqAlrysOae3z3NP6fDh82w68qVZsijSZ/fmXl0CMVdivNJu0+oXry6XX4mQiSKs8Xx66lf\nWbx/MRtPbqRNpTYMqDOANpXa2KWOYHb1yisQEgLe3pmzd15kbxkSyJRSeYFtQB4gJ7BKa+2llPIC\nXgauJTx0jNZ6Q8JzRgMvAvHAEK31rwnH6wKLgLzAL1rroXd4z0wXyBb5LeKjPz/C9xXfDJ0EfOnm\nJT7Y/gErD69kaMOhDGs0jAJ5Hv4jZlRcFMeu3wpeu84cZve5Q4RwgcK2ijSoUJ0mlatTrVg1qhev\njmdRz7tuTxMeE86p4FMcDzzOicATHA9K+Bp4nIjYCDxdPVOEtSquVfB09cTV2TUprGlt5pocP57y\ncuyY6RkbNAg6/t95Ptr7Nr6XfJnZdibPVX0u2+zvKKwTHBmM9yFvFu9fzNkbZ5OGNGsUr2F10zKd\nqCho1sz0yg8fbnVrhHgwGdZDppTKp7WOUErlBP4EhgLtgFCt9cz/PLYasAyoD5QCNgGeWmutlPIF\n3tRa+yqlfgE+0Vr7pPJ+mSqQ7fXfS9ulbdk2YJvdJtw/qFNBp/Da5sVvp35jZNORDKo/6K7bMSUG\nr8PXDpvwlRDA/tvjlRi8itg8+XpBLj77zKyKGj4c2rZN25ZCN6JuJPWmnQg8weGrxznkf4IzN4+j\nbU7ki/KEwCqEn/fEKbgKFQt5Ut3dkxqeBalSBapUgfKVI/l0z3Tm7JrD4AaDeafpO+TLJQXFRMY7\ndv0Yi/cv5tsD31LCpQTdqnWjYJ6C5FA5yOGUI8VXJ+V027E7fb3fx+bLlS/FB5nM6Nw5aNgQvv8e\nnnzS6tYIcXc2bePP83/yy4lfmPr01IwdslRK5QP+AAYB7YEwrfWM/zxmNGDTWk9NuO0DeAHngC1a\n60cTjvcCWmitX0/lfTJNIAuKDKLel/WY+tRUulfvbnVzOHj1ION/H8/uy7t578n36F2zN6eDTycN\nMyYONZ4POZ8UvKoVq2YCWPHqVC5amdw57lyXKzraDCnMnGkKrA4bZrZpcna+v/bFxsKZM7f3dh0/\nDkFBULkyeFbRlKlyHZeyx8HtBOF5jnMp0vSunQw6SYHcBZJ61rae3crj7o8zo80M2Q5HOIR4Wzxb\nTv7Ggd+WcqZMfmKVJl7Hm4vt3l9t2nbfj03+NSwmDK01VVyr3HbxLOqZpp7zjLRxo9laafduZNcM\n4XC01uy6tAvvg96sPLySos5FzXzu5uMyrIfMCdgLVALmaq1HK6UmAAOBEGA38JbW+oZS6lPgb631\ndwnPXQBsAM4CH2mtn044/gTwrta6Yyrvl2ogy8yf/IT9XAy5aIZAg05QxbUKLcq3sLpJQtzi5wev\nvgqXLpnx9r594YUXoGbNdH/rwIhAjgcev3UJOp40XaBw3sK3BbVHXB+hQpEKd/0QZoUPPgAfH/j9\nd8h15xkRQmQIrTV7/ffifcib7w99j3MuZ3pW70nP6j15tNijgAWT+pVShYDVwGDM3LHE+WOTAHet\n9Uv2CmQTJkxIut2iRQtatGiR+A0/3HcpsgQ5B4TDCg8HLy9YsgQ+/BAGDjRVib/91lzc3KB/f+jd\nO8P3C7JpGxdvXkwZ1hIuF29epEyhMiakFU0Z2EoVLGVJqRibDTp3hkqVzKIdITKa1pqDVw/ifcgb\n70Pe2LSNXtV70bNGT2oWr8m2bdvYunVr0uPff//9jF9lqZQaD0QkH6pUSpUH1mmtayqlRiV8Mx8l\n3OcDTMAMWf6ebMiyN9D8QYYs5Y+xkHNAOKRffoH//Q+eeAKmT4fixVPeHx8PW7easLZmjZm9/sIL\n0KnT3evKZICY+BhOB59ONayFRIdQuWjlVMOaaz5XtDbDsTHxMcTGxxITH5PqJdZ2l/vu8Lyb4TEs\nXBRL4ydiqFAphvoe9Rn42EBLf1Yi6zt6/SjeB00IC48Np0e1HvSs0ZO67nXvOkqXUass3YC4hOFI\nZ2Aj8BGwV2t9JeExw4H6Wus+ySb1N+DWpP7KCZP6dwFDAF/gZx5wUr/8MRZyDgiHcvmymUi5b5+p\nrPzUU/d+TlgYrF4Nixeb53XrZsJZkyYOV+8hNDo0qYRN8suxwGNExkYSZ4vDSTmRO0ducuXIRe4c\nuVO95HJK/b6k5zilft/1q7lZ8Hluhg3OxepL82hZviWz282W8iLCrk4Hn04KYVfDr9K9Wnd61ehF\nw9IN77uHOKMCWU1gMZADcAK8tdYfKKWWAHUADZwBXtNaByQ8Zwym7EUcMFRrvTHheGLZC2dM2Ysh\nd3hPCWQiVXIOCIcQHw9ffGF2yH7tNRg79v5XtiR34QJ8950JZ7GxZoVMv35QsaL922xHWmui46PJ\n5ZQr3Qs/L1kCkyfDpj9DeMmnOzmccuDdzVv2GBVpcj7kPN8f+h7vQ96cDzlP10e70rN6T5qVbfZQ\n53S2KwybVf8YOzk5cfLkSSpWrMigQYMoVaoU48aNA2D+/Pl4eXkRGRnJuXPnKFLEflXBv/vuO5Ys\nWcLGjRvt9prpLaueAyITOXDATNrPmdOEsup2KD6stdlDaPFiWLECHn3U9Jp17w6FCqX99TO5N94w\nRaCXrYhl4u7B/HXhL9b3Xi8rq8UD8Q/1Z+XhlXgf8ubo9aM8V/U5etXoRYvyLdLc6yqBzMEsW7aM\nmTNncuzYMQoUKECdOnUYO3YsTZs2vevzkgey5GJjYylUqBC+vr7UqCGFJsHxzwGRhUVEwMSJ8PXX\nMGWKqc2QlkJ8dxITAxs2mK6hTZvMRo/9+5sNZnNmz6G6+Hh47z1YuhRWrtT8FT+bj3d8zOqeq2lY\nuqHVzRP3cOz6Ma5HXMcltwsuuVxSfE3v4edr4df44cgPrDi4gv0B++n0SCd6Vu/JUxWfsuvqYglk\nDmTmzJlMnTqVL774grZt25I7d258fHzYvn07U6dOvetz7xTILly4QLly5YiNjSVHjgfrQk38OWW1\nUiGOfA6ILMzHx3TTNG5sivBl1CrJwEBT+G/JElM1tU8f03NWu3bGvL+D+ekn0zk5eTKUfHIdL659\nkXnPzHOIGpDidlpr5vrOZdL2SVQuWpnw2HDCY8JTfM3plDMpnOXPnf+2wOaSK+X1/Lnzpxrskr8G\nwIYTG/A+5I3vJV/ae7anZ/WetKvc7q4F09NCApmDCAkJoXTp0ixatIiuXbvedr+vry9Dhw7l6NGj\nODs707VrV2bOnEmuhAI7yQPZgAEDKFOmDC+88AJ16tQhMjISFxcXGjZsyKZNm9ixYwdDhw7lxIkT\nVKlShTlz5tC4cWPAlAdp1qwZv//+O35+fhw4cABPT0/mz5/PjBkzuHbtGn379mXu3LkALFq0iIUL\nF/LHH39wP9dTAAAgAElEQVQAMHToUFavXk1ISAienp7Mnj2bZs2aZdBP8f446jkgsqgrV8yk/X/+\nMZP227Sxri3Hjt0qoVGkiAlmffpAyZLWtckCx47B88+bNRAvj/Oj+4+dGFRvEKOajcpyH0Azs+i4\naAb9PIg9/nv4qedPVChS4bbHaK2JiotKNagl/xoWE5b6fcmuh8WEJV2PiY+hZfmW9Kzek2erPJsh\nu7fYK5ChtXaoi2nS7e503GobNmzQOXPm1PHx8anev2fPHr1r1y4dHx+vz549qx999FE9e/bspPuV\nUvrUqVNaa60HDBigx48fr7XW+uzZs1oplfS6gYGBunDhwnrp0qU6Pj5eL1++XBcpUkQHBQVprbVu\n3ry5LleunD58+LCOj4/XMTExWimlO3bsqENCQvT58+d1sWLFtI+Pj9Za62+++UY3a9YsqR1Lly7V\nQUFBOj4+Xs+YMUOXLFlSR0dH2/8HlgaOeg6ILCY+XuvPP9e6WDGtR4/WOjzc6hbdEh+v9ZYtWg8Y\noHXhwlq3b6/1ihVaR0RY3bIMExqqdffuWterp7XvkUv68S8e1/1X99fRcY71+yq7unTzkm74VUPd\n7ftuOjQ61OrmZIiEv01pzj8ZX/UvnShln8uDCgwMxM3NDac7zCd5/PHHadCgAU5OTpQrV45XX32V\nbdu23fH1dEIPUOLXRD///DOPPPIIffv2xcnJiV69elG1alXWrl2b8P0rBgwYwKOPPoqTk1NSD9yo\nUaMoWLAgZcqUoWXLlvj5+aX6vn379qVIkSI4OTkxYsQIoqOjOXbs2AP/PITI1A4eNPXEliyBLVvM\nfLF8DrQ/qpMTtGwJ33wDFy+anQAWLgQPD3N99WqIjLS6lekqf34zkturF3Rq6cH75bcTEh3C098+\nTWBEoNXNy9b+vvg3Db5qQMcqHfm+2/dJQ4ji/mSZQKa1fS4PytXVlevXr2Oz2VK9//jx43To0AF3\nd3cKFSrE2LFjCQx88F8aly9fpmzZsimOlStXjsuXLyfdLlOmzG3PK5lsSCNfvnyEh4en+vrTp0+n\nWrVqFC5cmCJFihASEsL169cfuJ1CZEqRkTBmDLRqZYYD//gDHH0xjYuLCWG//gpHjpiCs3Pnmk0g\ne/WCVavMDgJZkFLw1luwbBm80t+Fhmd/oIFHQxotbMTxwONWNy9b+nrf13Ra3onPO3zO2CfHyhDy\nQ8gygcwqjRs3Jk+ePKxevTrV+wcNGkS1atU4efIkISEhTJ48+Y7h7W5KlSrFuXPnUhw7d+4cpUqV\nSrr9sP8B/vjjDz7++GNWrlzJjRs3CA4OplChQjJfS2QPv/5qwtfp06asxWuvpc8KyvRUsiQMGgSb\nN8OJE9C6NXz5pek569bNdCmFhVndSrtr2RJ8fWH1j06c+XIaQx5/lye+eYKtZ7da3bRsIzY+lsG/\nDGbqX1PZPnA7Hap0sLpJmVYm+63jeAoVKsTEiRP53//+x5o1a4iIiCA2NpYNGzYwcuRIwsLCKFCg\nAPny5ePo0aPMnz//jq91twD0zDPPcPz4cZYvX05cXBze3t4cPXqUDh063NfzE+9P7TGhoaHkzJkT\nNzc3YmJimDhxIjdv3ryP716ITCwgwPQwvfaa6VlasSJrTJIvVgxeecUEzdOn4ZlnYNEiE86ee84U\no81C/7/LlIHt28HVFea9/Aof1V9Gj5U9+GbfN1Y3Lcu7HnGdNkvbcPrGaXa9vIuqblWtblKmJoHM\nDkaMGMHMmTP54IMPKF68OGXLlmXevHk899xzTJ8+nWXLllGwYEFeffVVevXqlaIn67/X73Rf0aJF\nWb9+PTNmzMDNzY3p06ezfv16ihYtmurj73Q78Vjy6+3ataNdu3ZUqVKF8uXL4+zsfNvwqBBZhs0G\nCxZAzZpQujQcOmTqfWVFrq6mZtqGDaZ0RpcuJniWLm320lyyBG7csLqVaZYnj6nT+/bbMLJ7a8aV\n2sYHf3zAmM1jsOkHH5EQ9+Z3xY/6X9WncenGrO21lsJ5C1vdpExPyl6ITEPOAZFmhw+bHrHYWDOk\nV6uW1S2yRkgIrFtn5plt2WIWMnTrBp07Q7IPeZnR7t3mW+nU6xp7PJ/Do6A7i7sszpDyB9mF90Fv\n3tzwJp898xk9qvewujmWkzpkItuRc0A8tKgoU1H088/h/fdNKHvAgstZVmgo/PwzrFxpdgdo1Mhs\n29SlC7i5Wd26h3L9ulnXEK+icB34EufDT7C291pK5s8CQ9IWirfFM27LOFYcWsHqnqupU7KO1U1y\nCBLIRLYj54C4LzYbnD9vhiIPHjSX7duhYUOYPdvMpRKpCwszw5urVpkdCurXN91Nzz2XcTsU2Elc\nHIwbB8tXaNp9OImNAV+zrvc6apaoaXXTMqUbUTfo80MfIuMi+b7b9xRzKWZ1kxyGBDKR7cg5IFLQ\n2kzMTwxdiZdDh8ym3DVqmM2/a9SAxx6DOvJp/oFERJhQtmoV/PKL+Rl262ZK5bu7W926+/bjj6ZD\n9Ln3lvNT1FAWd1lMe88sOmcwnRy5doQu3l1oW6ktM9rMIFeOXFY3yaFIIBPZjpwD2VhQUMoer8Tr\nYAJX8kv16mZ7IWE/kZFm1eaqVbB+PTRoAJMmma+ZwJEjJkdWab0D3/JdGfvkWN5s8KbVzcoU1h1b\nx0trX2LqU1MZ+NhAq5vjkLJ0IDt67ShOyokcTjlwUk44KSfKFS4nf4yzOQlk2UBYmJl4nzx8HTxo\n5jkl9nYlvxQv/nBbbIiHFx1tVme+/74ZBp48Gao6frmD0FAYMABOBZ4h/LlnaV/lKWa2nUlOp5xW\nN80haa2Z8scU5u+ez6oeq2hUupHVTXJYWTqQVfm0CvG2eGzahk3biNfxXBxxUf4YZ3MSyLKQ6Giz\nU/R/hxuvXDF/3P/b41W2rAQvRxMZCZ9+Ch9/bFZnTphgioI5MK1h2jSYOe8GZUb0oESxXKzouoIC\neQpY3TSHEhYTxsA1A7kQcoEfe/6IRwGZd3k3WTqQyZClSI2cAw4qMhICA83StsDA26+ndl9kJHh6\nppznVaMGVKokqx8zm+BgE8q++AIGDoTRo039Mwe2eTP06RdLxf8NJrzoDtb3WU/ZQlJ7EeBM8Bk6\nr+hMPY96zHt2Hnlz5rW6SQ5PApnIduQcSGdam3Gdu4Wr1IKWzWb+ALu5ma/3c71gwcy3PZG4u8uX\nzbyylSth+HAYNszst+mgzp+H57tq4urN4mqlGazp/RP1S9W3ulmW2nx6M31/7MvYJ8wcO9mP8v5I\nIBPZjpwDdnbzpikKumGDqT914YIpef6g4SpfPhlOFLecOAHjx5tSI+PGwcsvQ+7cVrcqVVFR8Oab\nsPHsWsKfeomvOn1O12pdrW5WhtNaM2fXHD768yOWd11OywotrW5SpiKBzEHkz58/6VNEeHg4efPm\nJUfCkMuXX35J7969rWxeluKo50CmoTX8+68JYD4+pqR5o0Zm26C2baFyZRPIhLCHvXthzBgT0CZN\nMpVaHbRX9Kuv4N1Ze8nxf515+8n/MbLpyGzTOxQVF8Xr619n35V9rOm1hvKFy1vdpDSJjDT5PyNn\nPkggc0AVKlRg4cKFtGrVyuqmZEmZ4RxwODdumN6vxBCWN68JYO3aQcuWDj2kJLKI33+HUaPMQo4P\nPzTnngOGHV9f6PLCJWw9O1LS3UZd97rUKVmH2iVrU7tEbQrlLWR1E+3u0s1LPOf9HOULl+ebzt/g\nkjtz/j64fBl++gl++AF27DA7oxUs+OCd/Q/bkSuBzAElD2Q2m41p06axYMECbty4QevWrfn8888p\nUqQIgwYN4tq1a6xatQqAkSNHsmfPHjZt2kRwcDD9+vXD19eXuLg4mjZtyueff06pUqUAWLRoEZMm\nTeLatWu4ubnxwQcf0KdPHyu/7QyTGc4By9ls4Od3K4D5+UGzZiaEtW9vesEc8I+hyOK0Nn8xx4wx\npUo++ggaN7a6Vbe5ehV69I4mquheegzx43jIfvYH7OffgH8p5lLMBLQSJqDVKVmH8oXLZ9qetB0X\ndtB9ZXferP8mo5qNynTfx5kzpujvjz+aSjnPPgtdu5rO/jx5zFqTB5kOGxhoPq8+SIBLnLHh5CSB\nzOEkD2Rz5szB29ubVatWUaxYMQYPHszNmzdZtmwZkZGR1KlThzFjxlCxYkWef/559u/fj4eHB0FB\nQWzbto327dsTFxfHiy++SGxsLKtXryY8PBwPDw92796Np6cnAQEBBAYGUq1aNau/9QyRGc4BSwQG\nmqKdPj6wcaOpUp/YC9a8OTg7W91CIYy4OFPDzMsLHn/c1DCrXt3qVqUQGwtDhsAff5gauOXLmz0c\nTwWfwu+KH/uvmJDmd8WP0JjQpIBWu6QJadWLVcc5l2P/n1uwdwFjNo/hm87f8GyVZ61uzn07csQE\nsB9+gIsXTbWV55+H1q3TPk1RazOt9n4WiSe/rhRERmZAIFNK5QW2AXmAnMAqrbWXUqoo4A2UA84C\nPbTWNxKeMxp4EYgHhmitf004XhdYBOQFftFaD73Dez5UIFPv2yfd6wkP/wc/eSCrVq0ac+fOTRq+\n9Pf3p1y5ckRFReHk5ISvry/t2rWjYMGCTJ06lZ49e6b6mn5+frRq1YqgoCDCw8MpXbo0CxcupH37\n9jhnsz+0EsgSxMfDnj2mF2zDBvPxsEULE8DatYOKFa1uoRB3FxUFn30GU6earg0vLyhXzupWJdEa\nPvnEdOT9+OOdO/OuR1xPCmiJIe144HEqFqmY1IuWGNYyYmPzqLgogiKDki6BEYEpbgdFBnE25Czn\nbpxjTa81POL2SLq3KS20Np38P/xg/h1u3jTbqnbtajr+czpATd+ICHBxyaAeMqVUPq11hFIqJ/An\nMBToClzXWk9TSo0EimitRymlqgHLgPpAKWAT4Km11kopX+BNrbWvUuoX4BOttU8q75clesjy5ctH\nrly5cEo2iTU6OppTp07hnrAPXL169bh+/TqnT59OelxERATDhw9n48aNBAcHAxAWFkZcXBxKKX79\n9VemT5/OP//8Q9OmTZkxYwaPPOLY/6nsJTOcA+nm6lXT++XjY3rDihe/1Qv2xBMyGV9kTiEhMH06\nzJsHL7xghjSLOc6m1T//bKr7f/IJ3O/6rOi4aI5cP5KiJ21/wH5yOeUyvWgl6iT1plVxrZLqTgGR\nsZEERt4eplIErKjbj8fZ4nDN50pR56JJF1fnlLeLOhelXeV2FMxT0L4/LDux2eDvv28NRyplAtjz\nz5uduhxxXUiGzyFTSuUD/gAGAUuA5lrrAKVUSWCr1rpqQu+YTWs9NeE5PoAXcA7YorV+NOF4L6CF\n1vr1VN4nSwSyqlWr8s0339D4Dh+tPvvsM7744gvy589Pp06dGDVqFACTJk1iy5YteHt7U7x4cfz8\n/Hj88ceJi4u7LdyNHTsWX19ftm/fniHfn9UywzlgN3FxZpZx4lywEyegVatbvWBlpYilyEKuXIEP\nPoDly2HoUFPHrIBjVM8/cAA6dTI1b9977+GmYGqtuRR66bYhz0uhl3jU7VFccrukCF5a6xTBKrVQ\nldp9+XLly3RzwcD8utu+3fSErV4NRYveCmG1ajn+tFd7BbJ7dvgppZyAvUAlYG5CD1cJrXVAwkMC\ngBIJ1z2Av5M9/SKmpyw24XqiSwnHs6zXX3+dMWPGsHjxYsqWLcu1a9fYuXMnnTp14vjx44wfP55t\n27bh7OxMgwYNaN++PbVr1yYsLAxnZ2cKFSpEUFAQ77//ftJrXr16lZ07d/LUU0/h7OyMi4tLUokN\nkY60NhMWrl41Qy1RUWbFWPKv9zr2oM+JiTG/idq3N1XQmzRx2FpOQqRZyZIwd64JYu+9Z3ZxGDMG\nXnvN8t7fWrVMj03nzma3r6+/NpO/H4RSitIFS1O6YGk6VOmQdDwsJox/A/4lOj46RbBy9Dlo9hAd\nbRaA//gjrF1rRqy7djWLcrPJoM9t7hnItNY2oI5SqhCwWilV4z/3a6WUXbstvLy8kq63aNGCFi1a\n2PPlM8TQoUPRWtOmTRsuX75M8eLF6dWrF88++yz9+vVj1KhR1KxZE4ApU6bQr18/9uzZw7Bhw+jT\npw9ubm6UKlWKESNGsHbtWgBsNhuzZs2if//+KKV47LHHmD9/vpXfZtYTHm72VDxwIOUlTx7w8DC/\nifPmNbeTf/3vsQIFzNDLne6/17HcuR2zb16I9FSpEnz3nZk4NHYszJoFEydCnz6WbqlVsiRs3WqG\nL1u2NAtGS5S417PuLX/u/DQu43irTdNLeLjp8P/hB9P5X6OGCWHvvedQUwjvaevWrWzdutXur/tA\nqyyVUuOBCOAVzJDjFaWUO/B7wpDlKACt9UcJj/cBJmCGLH9PNmTZGzPkmaWGLEX6sus5YLPBuXOw\nf3/K4HXxotnculYtc6ldG2rWNHO2hBAZa/t2U8MsTx7TjWLxMKbNZtYffPstrFtnAoW4u5AQ87P6\n8Uezh2jDhmYosksXE3SzggyZQ6aUcgPitNY3lFLOwEbgI6AFEKi1npoQwgr/Z1J/A25N6q+c0Iu2\nCxgC+AI/kwUn9Yv09dDnwM2bt3q9EgPYv/+a8hCJwSsxfHl6Qq5c9m+8EOLh2GwwaJD5f7thAxQu\nbHWL+O47M7q6eLGZVSBSCgszIWzFCtOz2Ly56Qnr2NHMD8tqMiqQ1QQWAzkAJ8Bba/1BQtmL74Gy\n3F72Ygym7EUcMFRrvTHheGLZC2dM2Yshd3hPCWQiVfc8B2w2OHXqVm9XYvgKCDC1jpKHr1q1suZv\nBiGyIq1hxAjYts2sMnZzs7pF/PUXdOtmproNHmx1a6wXGWny8ooVZkF406Zmt6zOnc1n36xMKvWL\nbCfFOXDzpplnkny48eBBM2/rv8GrcmVL558IIexAa7NZ+U8/mdngCeWDrHT6tOn1adkSZs92jLpY\nGSkmBn77zYSw9etNrd9evcyQpKur1a3LOBLIRLajlEK/+CLs2gVnz5q5XbVr3wpeNWtm/Y9iQmR3\nkyebscLNm6FMGatbQ0gI9OhhSjN4e2f9X0FxcaajcsUKU6KialXo2dP0FjpARrZEhpW9EMKh1K0L\n//ufCV8y10uI7GfsWHBxgSefND1llSpZ2pxChUwB2aFDTXWa9euhQgVLm2R3NpvZtHvFCli50uTg\nXr1g714piWhP0kMmMg05B4QQSb74whST/e03001jMa1NKbUpU0xZhyZNrG5R2mgNu3ebEPb992Yt\nRa9epjescmWrW+dYpIdMCCFE9vXaa+DsbHaw8PEx0xYspJSZ3F+pkinpMHu2KZ+WmWhtFqB7e5sg\n5uRkQlhizTCRvqSHTGQacg4IIW6zcqVJQuvWQf36VrcGMKGmY0fo39/ULXP0rX+OHbsVwsLDTS9Y\nr17w2GOO33ZHYK8eMikFng7Onj2Lk5MTNpvtoZ7/4Ycf8sorr9zx/vLly7N58+aHbZ4QQmQd3bvD\nV1/Bs8/Cn39a3RrATHHdtcuUf+jd25SEcDRnz8LUqWZlZIsWEBgICxbAmTMwbZo5LmEsY0kPWRqV\nL1+eq1evJu0pqZRi48aNNG3a9LYNwe0l+Sbm2YmjngNCCAfw22/Qt6/ZoLx1a6tbA5gg9uKLJuSs\nWWOf7ZYeVkQE7NljJuf/9BOcOGGKtfbqZdZHSGWghyc9ZA5CKcX69esJDQ0lNDSUmzdv4p6Gtb/x\n8fF2bJ0QQmQTTz8Nq1aZLqmff7a6NYCZ4rZsGbRrZ7YM+vffjHlfreHkSVi6FN58E+rVMyUa33rL\n7A43YQL4+5t1ES1bShhzFBLIMsDly5fp1KkTrq6ueHp6smDBgqT7vLy86NatG/369aNQoUIsWrQI\nLy8v+vXrl/SYb7/9lnLlyuHm5saUKVNSvLavry+NGzemSJEieHh4MHjwYGJjYzPsexNCCIfx5JOm\n7sSLL5qljg5AKTOPbMoUs/7gl1/s/x43b5qybJMnm7lrxYuboLVmjSnB8ckncP06+PrCp5+agChV\ngxyPBDI7uNcwWq9evShbtiz+/v6sWrWKMWPG8Pvvvyfdv3btWrp3705ISAh9+/ZFJRu4P3z4MG+8\n8Qbfffcdly9fJjAwkIsXLybdnzNnTubMmUNgYCA7d+5k8+bNzJs3z/7fpBBCZAYNGpjJW2++aTad\ndBB9+piA9NJLJiA97OwLmw0OH4avv4ZXXjHz1dzdTa9XcDAMGAD79sGFC2a9w1tvmRIczs52/XZE\nOsg6ZS/sNfvwAf+XaK3p0qULORP2zGjZsiWzZs1Kuv/ChQvs2LGDDRs2kDt3bmrXrs3LL7/MkiVL\naNmyJQBNmjShU6dOAOTNmzdFwFu1ahUdO3akWbNmAEyaNIm5c+cm3f/4448nXS9Xrhyvvvoq27Zt\nY+jQoQ/4jQshRBZRp47pMmrTxkyeussiqYzUpImZw9Whg1nZOGfOvbdbCgoyCwR27oS//za9XG5u\n0KiRubz2mqn4kTt3xnwPIv1knUBm0WRvpRRr1qxJMcH+7NmzSdcvX75M0aJFcXFxSTpWtmxZdu/e\nnXS7dOnSd3z9y5cvp7g/X758uCbbJOz48eOMGDGCPXv2EBERQVxcHPXq1UvrtyWEEJlbtWqwdSs8\n9ZSZXT9kiNUtAswQ4o4dZjL9s8+achOFC5v74uLMPLO//7518fc31TwaNTLVPRo2NEOSIuvJHIHs\nIctHOAIPDw+CgoIICwsjf/78AJw/fz5FyFJ36d3z8PDgyJEjSbcjIiIIDAxMuj1o0CDq1q2Lt7c3\nLi4uzJ49mx8cZO6EEEJYqnJls/Fi69amp2zUKKtbBJjtltatg+HDTa9Zp04mfO3ZY7YlatQImjaF\nt982uVIm3WcPjhvILlwwy5h/+810PWdSZcqUoUmTJowePZrp06dz7Ngxvv76a5YtW3Zfz+/atSuN\nGjXir7/+on79+rz33nsp6puFhYVRoEAB8uXLx9GjR5k/fz7F5eOTEEIY5crB9u2mpywiAt5/3yEK\nbOXMaSbYL15saoKNGWOmvyX2lonsxzEn9VetakoEb9xoPtn884/VLXpgyXu9li9fztmzZ/Hw8OD5\n559n4sSJSUOcSqnbesiSH6tevTqfffYZffr0wcPDg6JFi1KmTJmkx06fPp1ly5ZRsGBBXn31VXr1\n6nXXHjchhMh2PDxMT9natabbyYHqGfbvbybkt2kjYSy7c8zCsLt3m0CWrKiqFAUVcg4IIdIkONjU\nfKhb1+wEng6Fu0X2Y6/CsI4ZyDJRpX6RceQcEEKk2c2bZpljpUpmryCZoCXSSCr1CyGEEA+qYEHY\nsAEuXTJbLUkhbeEgJJAJIYTIXlxczHyy8HDo1g2io61ukRASyIQQQmRDefOa7ZXy5DF1JyIirG6R\nyOYkkAkhhMiecuc2u3+XKAHt20NoqNUtEtmYBDIhhBDZV86csGgRPPooPP003LhhdYtENiWBTAgh\nRPbm5ATz50PjxtCypZnwL0QGu2cgU0qVUUr9rpQ6pJQ6qJQaknDcSyl1USm1L+HSPtlzRiulTiil\njiql2iQ7Xlcp9W/CfXPS51sSQgghHpBSMHMmPPec2a/oscfgnXdMgXKZXyYywD3rkCmlSgIltdZ+\nSqn8wB6gC9ADCNVaz/zP46sBy4D6QClgE+CptdZKKV/gTa21r1LqF+ATrbXPf54vdcjSYNCgQZQq\nVYpx48ZZ3RS7k3NACJEhYmPNDjGbNpnL3r1mh+/Wrc0WTPXqmaFOIbCwMKxS6idgLtAUCNNaz/jP\n/aMBm9Z6asJtH8ALOAds0Vo/mnC8F9BCa/36f56fKQPZihUrmDVrFocOHcLFxYUKFSrQv39/Bg0a\nZHXTsgxHPweEEFlUWJjZDzMxoJ0/Dy1amHD21FPwyCMOsT+msIYlhWGVUuWBx4C/Ew4NVkrtV0ot\nVEol7sLlAVxM9rSLmJ6y/x6/lHA805sxYwbDhg1j5MiRBAQEEBAQwOeff85ff/1FTEyM3d4nLi7O\nbq8lhBDiPuXPD888Y4Y0DxyAY8egZ0/Yt89sQlmmjNmU8ttv4fJlq1srMlJ8vN1e6r4DWcJw5Spg\nqNY6DJgPVADqAP7AjLs8/YF4eXklXbZu3Wqvl00XISEhTJgwgfnz5/P888/j4uICQJ06dVi6dCm5\nc+cmOjqat99+m3LlylGyZEkGDRpEVFRU0mt89dVXeHp64urqSufOnfH390+6z8nJiXnz5uHp6ckj\njzwCwLRp0/Dw8KB06dIsWLAAJycnTp8+DcCAAQMYP348AMHBwXTo0IHixYtTtGhROnbsyCWZrCqE\nEGlTogT07g0LF8K5c/D779CoEfz0E9SoAdWrw9ChsG6d2apJZA02Gxw9ytYxY/Bq2BCv0qXxSvib\nbxda63tegFzARmDYHe4vD/ybcH0UMCrZfT5AQ6AkcCTZ8d7A56m8lk7NnY5bbcOGDTpnzpw6Pj7+\njo8ZNmyY7ty5sw4ODtahoaG6Y8eOevTo0VprrTdv3qzd3Nz0vn37dHR0tB48eLB+8sknk56rlNJt\n2rTRwcHBOioqSm/YsEGXLFlSHz58WEdEROi+fftqpZQ+deqU1lrrAQMG6PHjx2uttQ4MDNQ//vij\njoyM1KGhobp79+66S5cu6fjTSF+Oeg4IIUSSuDitfX21njJF69attXZx0bpJE63Hj9d62zato6Ot\nbqG4Hzab1ufOab1qldYjR2rdqpXWBQtqXaGC1j16aP3xx1pv3ar1zZuJf5vuK0/d7XI/k/oVsBgI\n1FoPT3bcXWvtn3B9OFBfa90n2aT+Btya1F9Za62VUruAIYAv8DNZYFL/0qVLeeedd1L0ajVp0oQj\nR44QHR2Nj48PzzzzDAcOHKBixYoA7Ny5k759+3L69GleeuklihUrxkcffQRAeHg4RYoU4eTJk5Qt\nWxYnJye2bNlCixYtAHjxxRdxd3dn8uTJAJw6dQpPT09OnjxJxYoVGThwIKVLl2bSpEm3tdXPz49W\nrUbWtB0AACAASURBVFoRFBSUzj+V9OGo54AQQtxRZCT89det+WfHj0OzZrfmn9WoYcpuCGtdv24W\nciS/aA0NGpgFHfXrm8UcxYrd9lR7zSG7n2UiTYH/Aw4opfYlHBsD9FZK1QE0cAZ4DUBrfVgp9T1w\nGIgD3kiWsN4AFgHOwC//DWNpoew0tKkTgs/9cnV15fr169hsNpwS/lPt2LEDgDJlyhAQEEBERAR1\n69a99R5aY7PZAPD396devXpJ97m4uODq6sqlS5coW7Zs0usk8vf3p0GDBkm3S5cufce2RUREMHz4\ncDZu3EhwcDAAYWFhJonLBFQhhEh/zs63whdAYKAZ4ty0ydQ+CwmBWrXAw+PWxd095fW8ea39HrKa\nsDCzctbX91b4CgyEunVNABswAD77zMwNzMC/lfcMZFrrP0l9rtmGuzxnCjAlleN7gJoP0sD79aBB\nyl4aN25Mnjx5+Omnn3j++edvu9/NzQ1nZ2cOHz6Mu7v7bfd7eHhw9uzZpNvh4eEEBgZSqtSt9Q7J\nw5O7uzsXLlxIup38+n8fP2PGDI4fP46vry/FixfHz8+Pxx9/XAKZEEJYxdXVbGjerZu5fe6cWSRw\n+TL4+8OJE7Bt263b/v5mM/Tkge2/oc3DA0qWNPtypqfoaAgONuElKOj2S2rHg4NNoCxc2FyKFEn5\n9W7HChWCXLnS1uaYGLMQ459/bgWwM2egZk0Tvjp0gPffhypVLO+plEIqaVS4cGEmTJjAG2+8gdaa\nNm3a4OLiwoEDBwgPD8fJyYlXXnmFYcOGMXfuXIoVK8alS5c4dOgQbdq0oXfv3vTu3Zs+ffpQtWpV\nxowZQ6NGjZJ6x/6rR48evPjii/Tr14+yZcveNjSpb83FIywsDGdnZwoVKkRQUBDvv/9+uv88hBBC\nPIBy5czlTmw2E2wSA9rly+Zy9KjpaUu8feUKFCx499Dm7m6CW+Jr3i1IpXY8JgaKFk394uoKpUvf\nfrxIEYiKMltS3bhhAlryr4GBcPJkymOJ10NCTJi7W4j7b6DLlw8OH74VwA4dgkqVzJBjo0YweLAZ\nJs6dO+P+je+TBDI7eOeddyhVqhTTpk3jhRdewMXFhYoVKzJt2jSaNGlCo0aNmDhxIo0aNeL69euU\nKlWKN954gzZt2tC6dWsmTZpE165dCQ4OpmnTpqxYsSLptf/bk9WuXTuGDBlCy5YtyZEjB+PGjePb\nb78lT8InI6VU0nOGDRtGnz59cHNzo1SpUowYMYK1a9dm3A9GCCFE2jg5gZubudSqdefH2Wwm3CQG\ntMQAd+gQ/PbbrTB35Yp5zcQQlVq4Kl8+9eP58z/8EF7Jkg/+HK3Nhu93CnM3bsDZs+Dnd+tYWJip\nC1e/PvTqZXZcsOdKyHT0wIVh01tmm9RvtSNHjlCzZk1iYmKS5rBlVXIOCCFEGtlsJlTJtBW7saQw\nrHAMq1evJjo6muDgYEaOHEmnTp2yfBgTQghhB05OEsYclPwVz4S+/PJLSpQoQeXKlcmVKxfz58+3\nuklCCCGESAMZshSZhpwDQgghHI0MWQohhBBCZBESyIQQQgghLCaBTAghhBDCYhLIhBBCCCEsJoFM\nCCGEEMJiEsgyEScnJ06fPm11M1I1YMAAxo8fD8Aff/xB1apVLW6REEIIkXlIILOD8uXLU6JECSIi\nIpKOLViwgJYtWz70a7Zo0YKFCxemqV3Hjx+ne/fuFCtWjMKFC1O7dm1mzZqFzWZL0+umJvmWTU88\n8QRHjx61+3sIIYQQWZUEMjux2WzMmTMnza+jtcZms922h+WDOnXqFA0bNqRcuXIcPHiQGzdusHLl\nSvbs2UNoaOgDv979hDipESaEEEI8HAlkdqCU4u2332b69OmEhISk+pgdO3ZQv359ChcuTIMGDdi5\nc2fSfS1atGDcuHE0a9aM/2/v3uOsKO88j3+eOvdb9+nTDU3TNKACihAvCIhRY2tMxFFnTEajWScJ\nG2dncnFRk11jmE3UzLxmTBayWd1M9rW5jEnMsPGSnVGMEAU7MYmCaAAVkYtAQ3Nr6D597reqZ/84\n1YfTzeHSFzh09+/9etWrnnrqcuo0RfW3n6p6KhAI8NnPfpZXX32Ve+65h1AoxOLFi0vLvvTSS8yY\nMYO6ujruueee4+7TQw89xFVXXcXSpUtpbGwEYMaMGTz55JPU1tYCcPvtt9PU1EQ4HOaaa65h8+bN\npfUXLVrEF7/4Rf7sz/6MYDBIW1sb7733Hq2trdTV1TF79myef/75ip/d1tZGS0tLaXrq1KksW7aM\niy++mHA4zJ133kk2mwUgGo1y8803M378eCKRCLfccgsdHR0n+5ELIYQQo4oEsmEyd+5cWltbWbp0\n6THzurq6uOmmm7jvvvvo6uriK1/5CjfddBPd3d2lZZ588kl++MMfkkgkeOKJJ7j66qv5/ve/Tzwe\n57HHHist98ILL7B+/Xo2bdrEU089xapVqyruz+rVq7nttttOuM833XQT27dvp7Ozkzlz5nDXXXf1\nmb98+XK+8Y1vkEgkmDdvHrfccgsLFy6ks7OTxx9/nLvuuoutW7ee9GejlOLpp59m1apV7Ny5k02b\nNvHEE08AxZa3u+++m/b2dtrb2/H5fCcMmkIIIcRoJIFsmCil+Na3vsXjjz/O4cOH+8x74YUXOP/8\n87nrrrswDIM777yTCy64gOeee6607qJFi5g5cyaGYeB0OoHKlwAffPBBampqaGlp4dprr2XDhg0V\n9+fIkSM0NTWdcJ8XLVpEIBDA5XLx0EMPsXHjxj6XM2+99VauuOIKADZs2EAymeTBBx/E6XRy7bXX\ncvPNN7N8+fJT+vksXryYCRMmUFdXxy233FLa70gkwic+8Qm8Xi/BYJAlS5bw29/+9pS2KYQQQowW\nzmrvwHBpU23Dsp1W3TrodWfNmsXNN9/Mo48+ysyZM0v1+/btY/LkyX2WnTJlCvv27StNl1/i61Xp\nPrIJEyaUyn6/n0QiUXFf6uvr+2y/P8uyWLJkCc888wydnZ0YRjGbHz58mFAohFKK5ubmPt+h/z72\n/w4nUr7fPp+vtF4qleL+++9n1apVpRbDRCKB1nrI99EJIYQQI8WoCWRDCVLD6ZFHHmHOnDl89atf\nLdU1Nzfzq1/9qs9yu3fv5sYbbyxN9w8fQw0j119/Pc8++yyLFi2qOP8Xv/gFzz33HKtXr2bKlClE\no1EikUifVrnyfZg4cSJ79uzpE5R2797dp3uLwezzsmXL2Lp1K+vWrWP8+PFs2LCBOXPmSCATQggx\npsgly2F23nnncccdd/R54vLGG29k69atLF++nEKhwC9/+Uu2bNnCzTffXFqm/+XJxsZGduzYccLP\nOtFTjY888gh//OMfeeCBBzh48CAA27dv5zOf+Qw9PT0kEgk8Hg+RSIRkMsmSJUtOuO0FCxbg9/v5\nzne+Qz6fp62tjRUrVnDnnXeWlh/MU5aJRAKfz0dtbS1dXV088sgjA96GEEIIMdJJIDsNvvnNb5JK\npUotPPX19axYsYJly5bR0NDA0qVLWbFiBZFIpLRO/9age++9l2eeeYZIJMJ9991X8XPK+/7q79xz\nz+W1115j165dzJo1i3A4zG233ca8efMIhUJ89rOfZcqUKTQ3NzN79myuuOKKPtvqv22Xy8Xzzz/P\niy++yLhx47jnnnv4+c9/zowZMyouf6LWrfJl77vvPtLpNA0NDXz4wx/mxhtvlJYxIYQQY4462/qO\nUkrpSvuklJJ+rsY4OQaEEEKcbezfTUNuSZAWMiGEEEKIKjtpIFNKtSilXlFKvauUekcptdiujyil\nXlJKbVVK/UYpFS5b5+tKqW1KqS1KqY+X1V+mlHrbnjf0bu2FEEIIIUaBU2khywP3a61nAQuALyul\nZgIPAi9prWcAq+1plFIXAncAFwILgX9WR28K+gFwt9Z6OjBdKbVwWL+NEEIIIcQIdNJAprU+oLXe\nYJcTwHtAM/DnwE/txX4K3GqX/wJYrrXOa613AduBy5VSTUBIa73OXu5nZesIIYQQQoxZA7qHTCk1\nFbgUWAs0aq0P2rMOAo12eSKwt2y1vRQDXP/6DrteCCGEEGJMO+WOYZVSQeBZ4F6tdby8awKttVZK\nDdvjbw8//HCp3NraSmtr63BtWgghhBBi0Nra2mhraxv27Z5StxdKKRewAnhRa/09u24L0Kq1PmBf\njnxFa32BUupBAK31o/ZyK4GHgN32MjPt+k8D12itv9Dvs6TbC1GRHANCCCHONsPV7cVJW8jsG/J/\nDGzuDWO254DPAd+2x/9WVv+vSqnvUrwkOR1YZ7eixZRSlwPrgM8Ajw1kZ6XDUCGEEEKMRidtIVNK\nXQX8DtgE9C78dYqh6ilgMrAL+JTWOmqvswT4PFCgeIlzlV1/GfAE4AN+rbVeXOHzKraQCSGEEEKc\nbYarhWzE9NQvhBBCCHG2kZ76hRBCCCFGCQlkQgghhBBVJoFMCCGEEKLKJJAJIYQQQlSZBDIhhBBC\niCqTQCaEEEIIUWUSyIQQQgghqkwCmRBCCCFElUkgE0IIIYSoMglkQgghhBBVJoFMCCGEEKLKnNXe\nASGEGE0ypskHmQy1TifjXC7chvzdK4Q4OQlkQggxSFprdqTTrI3HeT0WY20sxjvJJJM8HuKmyeF8\nHr9hMM7lYpzbzTiXi/Fl5dJQNu1zOKr9tYQQVaC01tXehz6UUvps2ychhACI5vO8URa+Xo/F8Dkc\nXB4KsaCmhstrargsFMJvhyqtNdFCgc58ns58nkO5XKncWV4um3YpddzAVinMBRwOlFJV/skIMXYp\npdBaD/k/oQQyIYSooGBZvJtK8bodvNbGYrRnMswpC1+Xh0JM8nqH7TO11sRNs2JgO16Y00Cjy8Vk\nr5fJHg9TvF4me71M8XhKdUGnXAwR4nSRQCaEEMNoXzZbavVaG4vxZiJBs9tdCl8LamqYHQjgOsvu\nCUuaJgdyOfZkMuzOZmnPZNidydCezZbGfsPoE9Km2EGtt268241xkla23KEc3S9107WyCzNpMvVb\nUwnODp6hbynE2UsCmRBCDFLaNHkrkehz6TFhmn3C1/xQiDqXq9q7OmRaaw7n88eEtN2ZDO12OVYo\n0NIvpE12eJjytkntb1PoNXEy29LUXVdH3Q11WBmL9n9sp+GTDZzzrXNwj3dX+2sKUTUSyIQYRSyt\n2ZxM8rueHn4XjfKHWIwmt5uFkQgLIxHmh0I4z7KWmZHA0ppYocD+XI718XgpfG1OpbjQ7y+FrwU1\nNUzz+cbsvVgp02RPNsvurTGiq7pQa+LUvJbhyCSDt+bBK5eZHLjYSXPQxxSvl1l+P59w1OH/bicH\nnzzI5K9NZtLiSRgeOUbF2COBTIgRrGBZbEgkSgHs1Z4e6l0uPlJby0fCYT5cU8PebJaVXV2s7Oqi\nPZvl+ro6FkYi3BCJ0OzxVPsrnFGm1vQUCnTl8xyxx139xkcq1EULBQIOB+NcLi617/1aUFPDnGBQ\nnmYEzJRJtC1K16ouulZ2UeguELkhQt0NdUQ+FsHdWGz5MrXmQC5XalF7Ixbj6c5OfIbBokQdVy9L\n4ng/y3nfOY+GTzaM2WArxiYJZEKMIFnL4o1YrBTAXovFmOL1lgLY1bW1NJ0gZO3LZvmNHc5e6u6m\n2eMptZ5dWVuLZwS1nllasyuT4VAuR1ehUDFIlertcqxQIOR0EnE6qXe5iDidRPqN612uY+rCTudZ\nd89XNWmtSb6bpGtlF92ruom9HiM4J0hkYYTIDRGClwRRxqn9XtFasy4e5+lDh3i6s5NL34K/ftwi\n3ODhou/NoGZuzWn+NkKcHSSQibNCPpqn57c9dK/pptBVIDQ3RGheiOClQRy+sdsCkTRNXuvpKQWw\n9fE4MwOBUgC7qraW+kHen2Rqzfp4nBePHGFlVxfvpVJcEw6XAtq5Pt8wf5vBi+bzbEom2ZRIlMbv\nJJNEXC6a3O5TDldhpxOHtLoMSr4rT/fL3cVWsFVdGC6DyMJiK1jddXU4a4b+BGZvOHtq/0E6/+Ug\nf/mjAqnWABf+03lcNK1OWszEqCaBbIQ5kM2yLZ2moDWm1lgUWwpMe2xR/EXbp9y7zCCXB3ApVRwM\nA/dxyi6litN22aUU7rJy+bSRssj+MUGyLUrslR5SW1LUfLiGuuvqcDW4iK+PE3sjRmpzCt8MHzXz\nawjNK4a0wOwAhnN0tlZ05/P8oSyAvZNMcmko1OcSZOg0dT1wJJ/n5e7u0uXNkMNRCmet4XCpT6xy\n2tTkDubIdmTJdmTJ7c+hnApHyIGzxokj5OhTdtY4T3h/UMGy2J5Os7Ff+OoqFJgdCHBRIMDFwSAX\nBQJ8KBik9jR2w2DlLcykiZkwsZJWcZyxcAQdOGodOGudOGucKMfoDAna1MTeiJVawZLvJqm9urbU\nCuabfnrvldNa8/rebt75h500/d84a+50UXPvBG6b3MisQEDCmRh1JJCdpbTWtGezvBWP81YiwZ/s\ncdayuMDvx20YGIBDKQzAUKpv2R6XlikvD3B5DeQti7zWpSFXNl1ezlsWuQplnTGZ/I7F+etNLnzT\nYuo22DEd3pwDb14K22eC8hTDWtjp5Fyfj2k+H9OUhxk7YMK7Jv4NGTJvJsnsyRC8OEhoXoiaeTWE\n5ofwTTv7b6TOR/Nk27NoS+MMO3GGnRzxWvw+djSA7chkWFBTUwpg80OhqtyjpLVm44EYv32vk43b\nujncnuLShJdZMRctRwzcB02yHVnyh/I465x4mj14mj24m9xoS2PGTMy4SSFWwIyXlWMmKHDUOFAh\nB7mAIuWHbq/FYY/FAU8BI+Sgts5NfZ2XpnovkxsCNNV7cdY4j4a8GgfOUDHcaVOXglNp3L9cFqqO\nN7//urqgcQQcxQBmD4bXwEyYFHoKmD3F7+QIFMNZKaTVFv9tT7XOEXSc8uW90/nvrQvFcN39m2Ir\nWPfL3XiaPURuiBBZGKH2qtqq3Wyf/CDFm1/dSnJtnCf+RvH2Qhe3N47n9nHjJJyJUeOMBTKl1E+A\nm4BDWusP2XUPA38NdNqLLdFav2jP+zrwecAEFmutf2PXXwY8AXiBX2ut7z3O542YQGZpzbZ0uhS6\n3orH+VMigdcwmBMKcWkwyJxgkDmhEC0ez4g4+VgFi8RbCbrXdBNdEyX2Wgz/TD/h68LUXVdH7ZW1\nOALFoGFpTaEs2HUVCuxIp48OmUyp7Hc4uLDgYd4HTmZugaZ3TfwbMxgJi5q5oVJLWs28GjzNZ+6G\ndW1qsvuzZHdnybRnyLZnyezOFMt2nWVqjBY3aUOTixYwYiaeFBT8CmodeOtchCJuXOGyX+DlQ4U6\nR61jwK2Fx7RqdRwtl0/rgi6GrGY3RpObAw2ad2ty/NGfomuc4qJpYT5yfgMfbaw/YUtV3rJ4P5Uq\ntnbF47zTnWT7oQQ6bnIRfi40PUw3PZyTdzEx58KRtPqGuLiJGTMpxAt9g17MLH4fS+Pwl4WmgNEn\nRPUPVX2m+weusnUNj3HS/2va0qWAVugpUIjaQa1suhTeyqbL68ykWWw9rBDcHEEHWKAL+tghX6Gu\nwmDlrZMugwnKqXCGnYQ/Gi6GsI9Hzuj/oVMR/X2U7ffvIKlNXv1agP8zKUbQ4eD2ceMknIkR70wG\nsquBBPCzskD2EBDXWn+337IXAv8KzAOagZeB6VprrZRaB9yjtV6nlPo18JjWemWFzzsrA1nesngv\nlTra8pVIsDGRoMHlKoWuS4NBLg0GmTCCnoDTlib5TrIUwKK/i+Jt8RL+qB3APlKLKzy0vpi0/YRW\n/5C2I53m8L40Le+azN9xNKgZboVjToC6+TVMWhAmPL8WV2SQ91slTTJ7MmR2ZYjvThPdmSa5O022\nPYu1J4dxsEC+1iA50SA6weBQI+wfr2kfr/mg3mTvOAtXnZN6t5uZfj/XhMN8JBxmltePTphHf5FH\nC32HnsKJ63oKGF7j+KGt1kkhVugTuvKH8jjrj7ZqeZo9uCe6j5abi2Vn2Fnxl5vWmq3pdOnS5h96\nergkGOSGSITr6+qIFQp97vd6P5VissfDRfalxouCQS4OBpk8xD8utC6GEuVSI/qXsDY1hVhZSCsP\nbQkT5VAoZ7/BVaFuIPOdCsNllMoYjJifobY0B39xkJ1LdlJzZQ2xvxvP074enunsLIWzT40fz6xA\noNq7KsSAnNFLlkqpqcDz/QJZQmu9rN9yXwcsrfW37emVwMPAbmCN1nqmXX8n0Kq1/kKFz6p6IMuY\nJm8nk6VWr7cSCTYnk0zxeoutXqEQc4JBLgkGR1zHkVpr0tvTRNdEiyHslSjOWmepBSzcGi496n6m\n9BQKfFAW0vZvT1B4M4l/Y5ZJ75rM2AapeoP4xR6MOX7Cc0NMvbwOX9DFkf1puj9IkbBDlrknh7E3\nj3tfgeA+E1dKc7gRDoyHQ42QbHKQaXZgNbtQLR48LR4igeLN5fVOJw0uV7HsctHgchE6Te8J1Lqs\nhaZfcOv95e4IOkohq/eyouEavktPadPkdz09rOzqYk13NxGXqxS8LgoEmBUIVLz/TIihMJMme5bu\nYe9je5n4hYlM+loLb+oUT3d28nRnJyEJZ2KEORsC2X8EeoD1wFe11lGl1OPA61rrX9jL/Qh4EdgF\nPKq1/phdfzXwgNb6lgqfdUYDWbxQYIPd4tUbvran05zv9zPHbvGaEwpxUSAw6PfBWTmL3P4c+SN5\nlEthuA2Upzg2PAbKbde5T0+LQWZP5mgAWxNFa13scfujdYSvDeOdPHzv4htuWctiZzLNzk1RDq/t\nIf9mEv/GDJHtJmjI+xTxJoNMs4NCswvV4sbd4sE/1UftFC/1TT4aPMXAJeFCiLNHZm+GnUt20v1y\nN+f8/TlMWDQBbcBau4+z3nD2qXHjuKuxkel+f7V3WYiKhiuQDfZRpx8A37LLfw8sA+4e6s70evjh\nh0vl1tZWWltbh7Q9S2v2ZLNsS6XYlk6zLZ1mezrNllSKjmyWDwUCzAmFuKq2lv88aRKzA4FT6tdJ\na02hq0B234nv6SlEC7gnuHHVu4r3hmQtrJyFzmqsnIWVtdC5o5dxSiGtLKwdU+cxSiHumDo77OWP\n5ImuiVKIFghfGyZ8XZgpfzfltD9lNZw8hsEFoQAXXBmAK5tL9VbWwspbOIPy0mQhRiLvJC8zfzaT\n2Bsxtt+/nY7/1cF53z2PK66t44raWpaedx5rYzGe6uzkyj/9iatra/na5MnMr5H+zUR1tbW10dbW\nNuzbHVQL2fHmKaUeBNBaP2rPWwk8RPGS5Stllyw/DVwznJcsLa3psLuW2JZOl8LX9nSaDzIZ6p1O\npvv9TLefApzu83G+388Mn6/iK2msrEV2X5bcvsohK9tRnKc8qu89PWWXmHqn3ePdp/Q0lraKoaw8\npJWPK4W4inX22BFyEG4NE5gdqPrTYEIIcTxaazqf6eSDBz4geEmQc//7ufinHW0RS5omP96/n2V7\n9nCuz8cDLS0sjERGzB+WYnSr9iXLJq31frt8PzBPa/0fym7qn8/Rm/qn2Tf1rwUWA+uAFxjETf1a\na/blcn1aunpD1450mrDT2SdwTff7mebxcg4evFmwUhZmyiyNzbhZDF0VWrYKPcVWreOFLM/EYrn3\niUMhhBBDY2ZM9n5vL3uW7mHC5yYw5RtT+jxUlLcsfnnoEN/ZswcFPDB5MneMGyfveRVVdSafslwO\nXAM0AAcptni1ApcAGtgJ/K3W+qC9/BKK3V4UgHu11qvs+t5uL3wUu71YfJzP09t+fZCOaJr9PRk6\nezIcieWIxnMkEgVqsopxpoOGvJNw3qAmpwjkFJ4MqLTuE7isVLHVyPAaGH4Dh9/Rdxx04GkaWquW\nEEKI4ZU7mGPnN3Zy+N8PM/WbU2n626Y+3cRorVnZ1cW329vZlcnw1ZYWPt/UREDuExVVMKo7hn3s\nsjacfgfegAN/0EUo5CIcchOpceMPuiqHq0pjn4HhNSRYCSHECJTYlGD7V7aTfDtJ3UeLDyLVXV+H\nd8rRB5HWxmJ8u72d3/f08KWJE7mnuZkG95l9UlyMbaM6kJ1t+ySEEKJ60rvSRFdH6V7dTffqbhwh\nRymcha8N425w834qxdI9e3i2s5O/amzkK5MmMfUseq+rGL0kkAkhhBhztLY7s17dTffL3fS82oPv\nXF+xM+vr60jP9/F4935+tH8/CyMRHpg8mYuDwWrvthjFJJAJIYQY86y8RfyNeCmgxd+ME5oTwndd\nLWs+VOCfxnUyOxziay0tXBMOy5OZYthJIBNCCCH6MZMmPb/vofvl4uXN9I40sbleXvhQjr0LXHzu\n2qncOn4cDglmYphIIBNCCCFOInc4R/SVKF0vd9Px0hHS0Tyb5xpM+XgDN9w6mfA0eT2TGBoJZEII\nIcQApXelWbtiH5tfPMjE13N4A06aP15P48fqCV8Xxj3u9DyhqS2NNoudf+tC38FV78LwSF9qI5UE\nMiGEEGIINsbj/PiVXURXd3PTO26a1ufxnePFN80HJsVX3eWtYwJUxaFC0CofsEA5FcqliuPewaEo\n9BQIzAoQmhsiNK84+Gf6+/S9Js5eEsiEEEKIYbA7k+G7e/bwi44D3L2/luszQSJeF/U+NyGP89gQ\n1W8wXMYJ5yunAoPjPlBgpkwSf0oQXx8n9kaM+Btxsh1ZgpcEqZlXUwpqvmm+s7pfTa012b1ZUltS\npN5LFcdbU0ff0+w2ij9Ll/0zc5eVXar4TmbX0WUGu0552O0zLivj4NhlBvmzlUAmhBBCDKMj+Tzf\n7+jg1Z4e9mWzdGSzZLVmotvNRI+HZns80e2m2eMplSd6PMP+loBCT4H4m3Hib8RLQa0QLRC6rBjO\neoOaZ7LnjD85auUs0tvSxcC1JUXyvSSpLSnS76dxBB34L/Djn+nHf4Ef3wwfhtcotiDa72ruLfef\ntvJ2OVdWPpV1cv3WL9iXhwu61NJZXtd/3LsMihMHtkqBzqmYv3G+BDIhhBDidEqaJvuzWTpyEUI6\nFQAAB7ZJREFUOfZls+zL5ejIZvuWczk8ShVDW1lI6x/gJrjduIbw3s3codzRkPZGMaRhUbzMOfdo\nUHM3Ds99cPlovm9rl13OtGfwTvH2CV69Q/m7R0cabZ04sB0v0NVcWiOBTAghhKg2rTXRQqEUzvYd\nJ8Adyuepczr7hjaPh4jTSdDhIORwECwbQnZ90OHAbxjHtIRprcl2ZEutaL1jR8BRuhctNLc4uOoq\nByWtNdk92T6Bq7dsJsyjYassePmm+TDccn9bL7lkKYQQQowgptYcyuVKoa03rHUXCsQLBRKmWRri\n/cpZyyJQHtb6h7fesmFQ36FpeKdA7aYsvk1ZnJsyqEYX7jkB/JcF8WbB2pYphq/3UzhrnH1buuyy\np/nMXw4diSSQCSGEEGOEqTXJfkEtYZp9gtwx8+xxMlfAuyNPwzt5Jmw26fFYfNCiiZ7rIHuem2Cd\nmwaX65ihvt90yOGQgFaBBDIhhBBCDErOsjiSz3M4ny+NTzQcKRTIWtYpBbfyodKl1tFGApkQQggh\nzpiMaXKkUDh+aOs33ZnPo7XG73DgMQw8ShXHhoHXHpfXDde0z77nzmcY+BwOfIaBcRpDoQQyIYQQ\nQpzVUqZJyjTJak3Wso4Op2k6Yw8p0yRtWaTtabdS+O1w5jOMiuWTzi8LeL1lv2FwYTA4LIHMORw/\ncCGEEEKI/vwOB/5h7qNtoLTWZOxw1hvUUpZF+gTllL18LJcrlcvry8vDRVrIhBBCCCEGabguWUpH\nIkIIIYQQVSaBTAghhBCiyiSQCSGEEEJUmQQyIYQQQogqO2kgU0r9RCl1UCn1dlldRCn1klJqq1Lq\nN0qpcNm8ryultimltiilPl5Wf5lS6m173v8c/q8ihBBCCDEynUoL2b8AC/vVPQi8pLWeAay2p1FK\nXQjcAVxor/PP6mgXvT8A7tZaTwemK6X6b1OIAWlra6v2LogRRI4XcarkWBHVcNJAprV+FejuV/3n\nwE/t8k+BW+3yXwDLtdZ5rfUuYDtwuVKqCQhprdfZy/2sbB0hBkVOmmIg5HgRp0qOFVENg72HrFFr\nfdAuHwQa7fJEYG/ZcnuB5gr1HXa9EEIIIcSYN+Sb+u1eXKUnVyGEEEKIQTqlnvqVUlOB57XWH7Kn\ntwCtWusD9uXIV7TWFyilHgTQWj9qL7cSeAjYbS8z067/NHCN1voLFT5Lwp0QQgghRoxqvsvyOeBz\nwLft8b+V1f+rUuq7FC9JTgfWaa21UiqmlLocWAd8Bnis0oaH40sJIYQQQowkJw1kSqnlwDVAg1Jq\nD/BN4FHgKaXU3cAu4FMAWuvNSqmngM1AAfhS2YspvwQ8AfiAX2utVw7vVxFCCCGEGJnOupeLCyGE\nEEKMNae9p/7jdCx7sVLqNaXUJqXUc0qpkF3vVkr9i12/QSl1TYXtPVe+LTF6DNexopS6Qym1USn1\njlLq0Wp8F3F6KaValFKvKKXetf+dF9v1A+60umy+nFtGqeE8XuT8MroN9Fix619RSsWVUo8fZ5un\ndG45E69OqtSx7I+AB7TWFwH/D/ivdv1/Aiy7/mPAsrKOZVFKfRKII091jlZDOlYAlFL1wHeA67TW\ns4EJSqnrzsTOizMqD9yvtZ4FLAC+rJSaycA6rS6d/+TcMuoNx/Gi5PwyJgzoWAEywH8D/kuljQ3k\n3HLaA9lxOpadbtcDvAz8pV2eCbxir9cJRIG5AEqpIHA/8A+A3Pg/Cg31WFFKzQPOBbZprY/Yy60u\nW0eMElrrA1rrDXY5AbxH8UGigXRaPR/k3DIWDNPxcjlyfhn1BnqsaK1TWus/ANn+2xrouaVaLxd/\nVyn1F3b5dqDFLm8E/lwp5VBKnQNcBkyy5/09sBRIndE9FdU20GNlG3C+UmqKUspJ8T9NS/+NitHD\n7pbnUmAtA+u0eqJdlnPLGDLE40XOL2PIKR4rvSq1gA3o3FKtQPZ54EtKqfVAEMjZ9T+heOCvB/4H\n8EfAVEpdApyrtf535C/YsWZAx4rWOgp8Efgl8DtgJ2Ce6Z0WZ4b9F+izwL1a63j5vFPotFrJuWVs\nGeLxgpxfxo6hHiuDObcMth+yIdFavw/cAKCUmgHcZNebwFd6l1NK/QHYCrQCc5VSOynu83il1Bqt\ntVy7H+UGcaygtV4BrLDr/4ZiFyxilFFKuSieMH+ute7tC/GgUmpCWafVh+z6Dvq2ZEyiGOgXIOeW\nMWEYjpcOkPPLWDDAY+V4BnxuqUoLmVJqnD02KN4M9wN72qeUCtjljwF5rfUWrfX/1lo3a63PAa4C\ntsoJc2wY6LFiT4+3x3UU/5r9URV2XZxG9sM+PwY2a62/Vzart9NqOLbT6jtV8encczjaabWcW8aA\n4Tpe7G3J+WUUG8SxUlq1fGIw55bT3kKmju1Y9iEgqJT6sr3Is1rrJ+xyI7BSKWVR/Ov1M5U2iTwJ\nNSoN47HyPaXUxXb5Ea319tO/9+IMuxL4K2CTUupPdt3XGVyn1b3k3DJ6DefxIueX0W1AxwqAUmoX\nEALc9j3PH+9tIOhdhFM4t0jHsEIIIYQQVVatm/qFEEIIIYRNApkQQgghRJVJIBNCCCGEqDIJZEII\nIYQQVSaBTAghhBCiyiSQCSGEEEJUmQQyIYQQQogqk0AmhBBCCFFl/x8C39XCMsXrKgAAAABJRU5E\nrkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "py.sign_in('sreejaya', '18o7jxvx8l')\n", "\n", "plots = FatalitybySTATE_dt.plot(subplots=False, figsize=(10, 4),)\n", "plt.legend(loc='center left')\n", "\n", "plt .show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "When we look at the top5 states we could not see any significant reduction in the traffic accidents." ] }, { "cell_type": "markdown", "metadata": {}, "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create choropleth map of the rate of accidents per State " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 58, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Anaconda2\\lib\\site-packages\\ipykernel\\__main__.py:15: SettingWithCopyWarning:\n", "\n", "\n", "A value is trying to be set on a copy of a slice from a DataFrame.\n", "Try using .loc[row_indexer,col_indexer] = value instead\n", "\n", "See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy\n", "\n" ] } ], "source": [ "#Collected State-County Data from Bokeh sample Data. Joined this fime with accidents file to get the state names.\n", "import plotly.plotly as py\n", "USAccidents_Data=pd.read_csv('https://raw.githubusercontent.com/srajeev1/IS608Project/master/data/USAAccidents2014.csv')\n", "\n", "#USAccidents_Data.head(10)\n", "\n", "Counties=pd.read_csv('https://raw.githubusercontent.com/srajeev1/IS608Project/master/data/US_Counties.csv')\n", "#Counties.head(10)\n", "\n", "Accidents_county = pd.merge(left=USAccidents_Data,right=Counties, left_on=['STATE','COUNTY'], right_on=['STATE num','COUNTY num'])\n", "\n", "Accidents_countydf =Accidents_county[['State Abbr.','STATE','COUNTY','VE_TOTAL','LATITUDE','LONGITUD','State-County','Geographic Name']]\n", "#Accidents_countydf.head(10)\n", "Accidents_countydf1 =Accidents_county[['State Abbr.','VE_TOTAL']]\n", "Accidents_countydf1['TotalAccidents'] =Accidents_countydf1.groupby(['State Abbr.'])['VE_TOTAL'].transform(np.sum)\n", "\n", "Accidents_countydf1=Accidents_countydf1.drop_duplicates(['State Abbr.','TotalAccidents'])\n", "Accidents_countydf1.head(50)\n", "\n", "for col in Accidents_countydf1.columns:\n", " Accidents_countydf1[col] = Accidents_countydf1[col].astype(str)\n", "\n", "scl = [[0.0, 'rgb(242,240,247)'],[0.2, 'rgb(218,218,235)'],[0.4, 'rgb(188,189,220)'],\\\n", " [0.6, 'rgb(158,154,200)'],[0.8, 'rgb(117,107,177)'],[1.0, 'rgb(84,39,143)']]\n", "\n", "\n", "\n", "data = [ dict(\n", " type='choropleth',\n", " colorscale = scl,\n", " autocolorscale = False,\n", " locations = Accidents_countydf1['State Abbr.'],\n", " z = Accidents_countydf1['TotalAccidents'].astype(float),\n", " locationmode = 'USA-states',\n", " text = Accidents_countydf1['State Abbr.'],\n", " marker = dict(\n", " line = dict (\n", " color = 'rgb(255,255,255)',\n", " width = 2\n", " )\n", " ),\n", " colorbar = dict(\n", " title = \"US Traffic crashes\"\n", " )\n", " ) ]\n", "\n", "layout = dict(\n", " title = '2014 US Traffic Crashes By States)',\n", " geo = dict(\n", " scope='usa',\n", " projection=dict( type='albers usa' ),\n", " showlakes = True,\n", " lakecolor = 'rgb(255, 255, 255)',\n", " ),\n", " )\n", " \n", "fig = dict( data=data, layout=layout )\n", "py.iplot(fig, validate=False, filename='d3-2014 US Traffic Crashes By States')\n", "\n", "url = py.plot( fig, filename='d3-2014 US Traffic Crashes By States' )\n" ] }, { "cell_type": "code", "execution_count": 59, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 59, "metadata": {}, "output_type": "execute_result" } ], "source": [ "\n", "from IPython.display import HTML\n", "HTML('')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### What causes most of the accidents?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1. Alcohol\n" ] }, { "cell_type": "code", "execution_count": 26, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from bokeh.plotting import *\n", "dfAlcohol=pd.read_csv('C:/DataAnalytics/MSDASemester3/608/Project/Data_proj/CSV/USAAlcohol1994-2014.csv')\n", "dfAlcohol.head(5)\n", "percents =dfAlcohol['Year']" ] }, { "cell_type": "code", "execution_count": 60, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import plotly.plotly as py\n", "from plotly.graph_objs import *\n", "\n", "fig = {\n", " 'data': [\n", " {\n", " 'labels': ['without Alcohol-21005','with Alcohol-11906'],\n", " 'values': [21005,11906],\n", " 'type': 'pie',\n", " 'name': '2010-Total(32999)',\n", " 'marker': {'colors': ['rgb(56, 75, 126)']},\n", " 'domain': {'x': [0, .48],\n", " 'y': [0, .49]},\n", " 'hoverinfo':'label+percent+name',\n", " 'textinfo':'none'\n", " },\n", " {\n", " 'labels': ['without Alcohol-20848','with Alcohol-11527'],\n", " 'values': [20848,11527],\n", " 'type': 'pie',\n", " 'name': '2011-Total(32479)',\n", " 'marker': {'colors': ['rgb(56, 75, 126)']},\n", " 'domain': {'x': [.52, 1],\n", " 'y': [0, .49]},\n", " 'hoverinfo':'label+percent+name',\n", " 'textinfo':'none'\n", " },\n", " {\n", " 'labels': ['without Alcohol-21563','with Alcohol-12118'],\n", " 'values': [21563,12118],\n", " 'type': 'pie',\n", " 'name': '2012-Total(33782)',\n", " 'marker': {'colors': ['rgb(56, 75, 126)']},\n", " 'domain': {'x': [0, .48],\n", " 'y': [.51, 1]},\n", " 'hoverinfo':'label+percent+name',\n", " 'textinfo':'none'\n", " },\n", " {\n", " 'labels': ['without Alcohol-20843','with Alcohol-11941'],\n", " 'values': [20843,11941],\n", " 'type': 'pie',\n", " 'name': '2013-Total(32894)',\n", " 'marker': {'colors': ['rgb(56, 75, 126)']},\n", " 'domain': {'x': [.52, 1],\n", " 'y': [.51, 1]},\n", " 'hoverinfo':'label+percent+name',\n", " 'textinfo':'none'\n", " },\n", " {\n", " 'labels': ['without Alcohol-20856','with Alcohol-11731'],\n", " 'values': [20856,11731],\n", " 'type': 'pie',\n", " 'name': '2014-Total(32675)',\n", " 'marker': {'colors': ['rgb(56, 75, 126)']},\n", " 'domain': {'x': [2, .48],\n", " 'y': [2, .49]},\n", " 'hoverinfo':'label+percent+name',\n", " 'textinfo':'none'\n", " }\n", " ],\n", " 'layout': {'title': '2010-2014 Alcohol Relatd Crashes ',\n", " 'showlegend': True}\n", "}\n", "\n", "#fig = dict(data=data, layout=layout)\n", "\n", "#py.iplot(fig, validate=False, filename='d3-Alcohol Related Accidents from 2010-2014')\n", "url = py.plot(fig, filename='d3-Alcohol Related Accidents from 2010-2014')" ] }, { "cell_type": "code", "execution_count": 61, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 61, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import HTML\n", "HTML('')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "** More than 35% of accidents are Alcohol related. The rate has no change in years. **" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 2. Weather Related ###" ] }, { "cell_type": "code", "execution_count": 62, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 62, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAFTCAYAAADRKgwqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcZFV9/vHPA6OCASEEREBZVPaAIAouiK0Ygj8VUREw\nimhQI6i4RYNbGIkSl7jFhaggi1ECBhfcEEQGUIIji4JsIooCsigKiAuCPr8/zqmZmqZnenr63Orp\nvs/79erXVJ2qut9bPV3fuvfcc75HtomIiH5ZZaZ3ICIiRi/JPyKih5L8IyJ6KMk/IqKHkvwjInoo\nyT8iooeWmfwlPUTSWZIuk/RDSYfW9vmSrpd0cf156tBr3iTpaklXStpjqH0nSZfWxz401H4/SSfV\n9vMlbdLFG42IiMUmO/K/G3it7W2BxwCvkLQ1YOD9tnesP18HkLQNsB+wDbAn8DFJqts6CjjI9ubA\n5pL2rO0HAbfW9g8A7274/iIiYgLLTP62b7L9/Xr7TuAKYKP6sCZ4yTOBE23fbfta4MfALpI2ANa0\nvbA+7wRg73p7L+D4evsUYPcVfC8REbGclrvPX9KmwI7A+bXpVZJ+IOkYSWvXtg2B64dedj3ly2J8\n+w0s/hLZCLgOwPY9wO2S1pna24iIiKmYtzxPkrQG8L/Aq23fKeko4Ij68L8B76N033RGUupQRERM\nke2JemkmP/KXdB9Kd8x/2/5i3dgtroCjgZ3r028AHjL08gdTjvhvqLfHtw9es3GNNQ9Yy/avl/Im\nlvvn8MMPn9LzV/QncVbOGImz8sZInNHFWJbJRvsIOAa43PYHh9o3GHras4BL6+1Tgf0l3VfSZsDm\nwELbNwF3SNqlbvMA4EtDrzmw3t4HOHOZexwREdM2WbfP44EXAJdIuri2vRl4nqQdKKN+fgr8E4Dt\nyyWdDFwO3AMc4sVfP4cAxwGrA1+zfVptPwb4tKSrgVuB/Vu8sYiIWLplJn/b32bis4OvL+M1RwJH\nTtB+IbDdBO13AftOuqdTNDY21nqTiTOLYiTOyhsjcVaOGJqsX2hlIcmzZV8jIlYGkvCKXvCNiIi5\nJ8k/IqKHkvwjInooyT8iooeS/CMieijJPyKih5L8IyJ6KMk/IqKHkvwjInooyT8iooeS/CMieijJ\nPyKih5L8IyJ6KMk/IqKHkvwjInooyT8iooeS/CMieijJPyKih5L8IyJ6KMk/IqKHkvwjInooyT8i\nooeS/CMieijJPyKih5L8IyJ6aN5M78B0SVqh19luvCcREbPHrE/+xVQT+Yp9YUREzBXp9omI6KEk\n/4iIHkryj4jooST/iIgeSvKPiOihJP+IiB5K8o+I6KEk/4iIHlpm8pf0EElnSbpM0g8lHVrb15F0\nhqQfSTpd0tpDr3mTpKslXSlpj6H2nSRdWh/70FD7/SSdVNvPl7RJF280IiIWm+zI/27gtba3BR4D\nvELS1sBhwBm2twDOrPeRtA2wH7ANsCfwMS2uv3AUcJDtzYHNJe1Z2w8Cbq3tHwDe3ezdRUTEhJaZ\n/G3fZPv79fadwBXARsBewPH1accDe9fbzwROtH237WuBHwO7SNoAWNP2wvq8E4ZeM7ytU4Ddp/um\nIiJi2Za7z1/SpsCOwHeB9W3fXB+6GVi/3t4QuH7oZddTvizGt99Q26n/Xgdg+x7gdknrTOVNRETE\n1CxXYTdJa1COyl9t+7fDlTRtW9JISmTOnz9/0e2xsTHGxsZGETYiYlZYsGABCxYsWK7narLSxpLu\nA3wF+LrtD9a2K4Ex2zfVLp2zbG8l6TAA2++qzzsNOBz4WX3O1rX9ecButg+uz5lv+3xJ84Abba83\nwX54on0tX0RTr+qZks4RMddJwvaEZYwnG+0j4Bjg8kHir04FDqy3DwS+ONS+v6T7StoM2BxYaPsm\n4A5Ju9RtHgB8aYJt7UO5gBwRER1a5pG/pF2Bc4BLWHx4/SZgIXAysDFwLbCv7dvqa94M/CNwD6Wb\n6Bu1fSfgOGB14Gu2B8NG7wd8mnI94VZg/3qxePy+5Mg/ImIKlnXkP2m3z8oiyT8iYmpWuNsnIiLm\npiT/iIgeSvKPiOihJP+IiB5K8o+I6KEk/4iIHkryj4jooST/iIgeSvKPiOihJP+IiB5K8o+I6KEk\n/4iIHkryj4jooST/iIgeSvKPiOihJP+IiB5K8o+I6KEk/4iIHkryj4jooST/iIgeSvKPiOihJP+I\niB5K8o+I6KEk/4iIHkryj4jooST/iIgeSvKPiOihJP+IiB5K8o+I6KEk/4iIHkryj4jooST/iIge\nSvKPiOihJP+IiB5K8o+I6KFJk7+kT0m6WdKlQ23zJV0v6eL689Shx94k6WpJV0raY6h9J0mX1sc+\nNNR+P0kn1fbzJW3S8g1GRMS9Lc+R/7HAnuPaDLzf9o715+sAkrYB9gO2qa/5mCTV1xwFHGR7c2Bz\nSYNtHgTcWts/ALx7Wu8oIiImNWnyt30u8JsJHtIEbc8ETrR9t+1rgR8Du0jaAFjT9sL6vBOAvevt\nvYDj6+1TgN2Xf/cjImJFTKfP/1WSfiDpGElr17YNgeuHnnM9sNEE7TfUduq/1wHYvge4XdI609iv\niIiYxIom/6OAzYAdgBuB9zXbo4iI6Ny8FXmR7VsGtyUdDXy53r0BeMjQUx9MOeK/od4e3z54zcbA\nLyTNA9ay/euJ4s6fP3/R7bGxMcbGxlZk9yMi5qQFCxawYMGC5XqubE/+JGlT4Mu2t6v3N7B9Y739\nWuDRtv+hXvD9LLAzpTvnm8DDbVvSd4FDgYXAV4H/tH2apEOA7WwfLGl/YG/b+0+wD55oX8v15Mnf\nw7hXsTzvOyJiNpOE7Ymuz05+5C/pROCJwLqSrgMOB8Yk7UDJuj8F/gnA9uWSTgYuB+4BDhnK2IcA\nxwGrA1+zfVptPwb4tKSrgVuBeyX+iIhoa7mO/FcGOfKPiJiaZR35Z4ZvREQPJflHRPRQkn9ERA8l\n+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflH\nRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0\nUJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCS\nf0REDyX5R0T00KTJX9KnJN0s6dKhtnUknSHpR5JOl7T20GNvknS1pCsl7THUvpOkS+tjHxpqv5+k\nk2r7+ZI2afkGIyLi3pbnyP9YYM9xbYcBZ9jeAjiz3kfSNsB+wDb1NR+TpPqao4CDbG8ObC5psM2D\ngFtr+weAd0/j/URExHKYNPnbPhf4zbjmvYDj6+3jgb3r7WcCJ9q+2/a1wI+BXSRtAKxpe2F93glD\nrxne1inA7ivwPiIiYgpWtM9/fds319s3A+vX2xsC1w8973pgownab6jt1H+vA7B9D3C7pHVWcL8i\nImI5zJvuBmxbklvszGTmz5+/6PbY2BhjY2OjCBsRMSssWLCABQsWLNdzZU+etyVtCnzZ9nb1/pXA\nmO2bapfOWba3knQYgO131eedBhwO/Kw+Z+va/jxgN9sH1+fMt32+pHnAjbbXm2AfPNG+lksKU/3u\nEcvzviMiZjNJ2NZEj61ot8+pwIH19oHAF4fa95d0X0mbAZsDC23fBNwhaZd6AfgA4EsTbGsfygXk\niIjo0KRH/pJOBJ4IrEvp3/9XSuI+GdgYuBbY1/Zt9flvBv4RuAd4te1v1PadgOOA1YGv2T60tt8P\n+DSwI3ArsH+9WDx+P3LkHxExBcs68l+ubp+VQZJ/RMTUdNHtExERs1iSf0REDyX5R0T0UJJ/REQP\nJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5\nR0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE\n9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQkn9ERA8l+UdE9FCSf0REDyX5R0T0UJJ/REQPJflHRPRQ\nkn9ERA9NK/lLulbSJZIulrSwtq0j6QxJP5J0uqS1h57/JklXS7pS0h5D7TtJurQ+9qHp7FNERExu\nukf+BsZs72h759p2GHCG7S2AM+t9JG0D7AdsA+wJfEyS6muOAg6yvTmwuaQ9p7lfERGxDC26fTTu\n/l7A8fX28cDe9fYzgRNt3237WuDHwC6SNgDWtL2wPu+EoddEREQHWhz5f1PSBZJeWtvWt31zvX0z\nsH69vSFw/dBrrwc2mqD9htq+UpG0Qj8RESujedN8/eNt3yhpPeAMSVcOP2jbkjzNGIvMnz9/0e2x\nsTHGxsZabXo5TfWtJPlHxOgsWLCABQsWLNdzZbfJzZIOB+4EXkq5DnBT7dI5y/ZWkg4DsP2u+vzT\ngMOBn9XnbF3bnwc80fbLx23fE+1rObqeelKe6vseVZyIiFYkYXvCo9AV7vaRdH9Ja9bbfwXsAVwK\nnAocWJ92IPDFevtUYH9J95W0GbA5sND2TcAdknapF4APGHpNRER0YDrdPusDX6j92vOAz9g+XdIF\nwMmSDgKuBfYFsH25pJOBy4F7gEOGDuUPAY4DVge+Zvu0aexXRERMolm3T9fS7RMRMTWddPtERMTs\nleQfEdFDSf4RET2U5B8R0UNJ/hERPZTkHxHRQ0n+ERE9lOQfEdFDSf4RET2U5B8R0UNJ/hERPZTk\nHxHRQ0n+ERE9lOQfEdFDSf4RET2U5B8R0UNJ/hERPZTkHxHRQ0n+ERE9lOQfEdFDSf4RET00b6Z3\nIJYkacqvsd3BnkTEXJbkv1KaSjKf+pdFRES6fSIieijJPyKih5L8IyJ6KMk/IqKHkvwjInooyT8i\nooeS/CMieijJPyKih5L8IyJ6KMk/IqKHkvwjInooyT8ioodWmuQvaU9JV0q6WtK/zPT+zGWSVugn\nIuaOlSL5S1oV+AiwJ7AN8DxJW09vqwumvV9zO46X8nPWUtrbWbBgQdPtJc7sipE4K0eMlSL5AzsD\nP7Z9re27gf8Bnjm9TS6Y/l4lzrQs7QziSU960kjOMObSB39UcebSe5lrceZq8t8IuG7o/vW1LWa9\nic4iDl9K+4qdYSztS+Ttb397urEilsb2jP8AzwE+OXT/BcCHxz3HEwEMnuDn8KW0l9dM1cobp2WM\nUcWZ6d/ZdOJM7WcUMRIn/zfLyo9eSt5VeXxmSXoMMN/2nvX+m4C/2H730HNmfkcjImYZ2xOe5q4s\nyX8ecBWwO/ALYCHwPNtXzOiORUTMUSvFGr6275H0SuAbwKrAMUn8ERHdWSmO/CMiYrRWltE+0QOS\n7j/T+xARxUrR7TNdKuP2Hmz7ukmfPP1Yq9r+c9dxBrGA9Rn6f7L988YxngJ8x/YfWm53XIzHAUcD\nawIPkbQD8DLbhzSOsxpl5NimLP6d2fYRrePY/uNkbdPY/qXLeNi2t28UZxXgMbbPa7G9KcRdh/J5\nvaSDbV8DnA+cC5xr+7LWMYZidfL5lHQnLHXcs20/YLoxYI50+9Tkf6ntvx1BrJ8ApwDH2r68wziv\nogyIvwVY9GVje7vGcU4AHgP8Bjin/nzb9m8axlgI7AN8yfaOte0y29u2ilG3+Q3gNuBClvydva9x\nnItsP3Kytmlsf9N6c/Dl+GlAwPMBbDcrfyLp+7Z3aLW9ZcQ5G3gGJVFeCPySctDx2sZxVgN2AXat\nP1tQcsPejeOM5PPZpTlx5G/bki6UtLPthR2H2wHYHzi6fvN/CjjR9h2N47wG2NL2rY23uwTbLwSQ\ntCElQX8U2JDGfxu2fz5uYtU9LbdfbWT77zvYLgCSNqD8bu4v6ZGUhGzgAUCzLi3b19Z4e4xLzJdI\nuhhoWfvqm5L2AU5xt0eCa9m+Q9JLgBNsHz7JGc6Kuge4m5KQ/0L5krm5gzidfT4lPaD+rtaZ6HHb\nv24RZ04k/+oxwAsk/Qz4XW1rdoo8UJP8J4BPSBoDPgN8QNLngH+z/eNGoX4OtP5CuRdJB1COkLan\nfFA+Any7cZifS3p8jXdf4FCgi9Fc50navovuhGoP4EWU2efDZxO/Bd7cQTxJ2tX2t+udx1O+cFp6\nOfA64M+SBt1WzboWhqxavzz3Bd46iNM4BpTPzKXA+4Gjbf+qgxjQ7efzROBpwEVM/DvarEWQOdHt\nA0ucKi9hcBTVMM48yn/Miyl9yycAn6Uk0CNtb9Eozqcop6xfBf5Um237/S22PxTnVuAa4Chgge2f\nttx+jbEe8CHgKZTkdTpwaOujJklXAA8HfgrcVZubHwBI2sf2/7bc5lLi7AQcC6xVm24DXmz7oq5j\ntybpucDbKF09B0t6GPAe289pHOeZwBOAR1POAM4DzrH9zcZxRvL57NKcSf4Dkh4IrDa438EF0p9Q\nqp8dPf5CmaQP235Vozjz683Bf5Aof1xvb7H9oTgCtqV8YJ5ASZ4/sv2ChjGeavvr49pebvu/WsWo\n29yECY6MOzgAeBDwTko3056StgEea/uYlnGG4q1F+aze1sG2V6FcS9jM9hGSNgYeNILu005J2gr4\nf5TumQfaXm2Sl0x1+/Prza4/n88Edqtxzrb95WbbnivJX9JelFPxDSkXYTYBrujgouIatu9suc2Z\nJOkBlLOW3erPusD5g2sBjWKcB7zN9pn1/huBJw/KeTSM8w7gbOA827+b7PnTiHMa5Yj8Lba3l3Qf\n4OLWAw5G8SUj6b8ofeNPtr1V7Wc+3fajGm3/w0N3zeIvZwPYPrRFnKF4p1Cuy11DGbxwLrCwy9Fs\nXZH0LsoZzGcov7f9gQtsv6nF9udSn/87gMcCZ9jeUdKTgANabXz4j7jWGRo+wnSrP2JJH7L9akkT\nfcPb9l4t4gz5NvAdyofkI7avb7x9gL2Ar0j6E2XNhq1qW2s/Af4B+M86XO4cynC/LzaOs67tkyQd\nBmD7bkldXMA+jvolU+9fDZwMtDzD2KV+Xi6GcjGxfpm1cmH993GUtTpOonx2ngt0MQzzXcBFXQ3H\nHvHn82nADoP3Iuk44PtAkv84d9v+laRV6lj8syR9qOH2R/VHfEL9t+nwxKUZ9IdLWpNuLsBR/1/2\nAs4ELgD26WJkie1PAZ+qR8z7Af8M/BOwRuNQd0r6m8GdWpjw9sYxYDRfMn+qo9aARddn/tJq47aP\nq9s9GNjVZb0OJB1F+4EFAD8AXilpt3p/AfBfg7gNjPLzaWBtYHBtbG0afkbnUvL/TU1g5wKfkXQL\n0Kx7ZlR/xLYvrP8uaLXNZZG0HeUP+m/q/V8CB9r+YYNtj5+scl/KSIV9JDUfUSLpGGBrytC+b1Mm\nfF3cMkb1euDLwENrl9Z6lGGyrY3iS+bDwBeAB0o6kvI+3rrsl6yQtSlDYgeJbM3a1tpRlLz2UcrB\n2QG17SUtNj7iz+e/AxdJGsR6InBYq43PpT7/NYA/UP7DX0D5Q/tMByNKrgIeN9hu7SP9P9tbNo6z\nBXAk5WLs4GKVbT+0cZz/A95s+6x6f4wyaulxLeOMgqQvUIZhXkbp8jnb9k86ijUP2JJSIuXKhkeW\nwzF2oiTnbSnvaT3KWdMPGsfZmlJRF+BMd1BUUdKLgfksXiruiZQy7sc1jnPJ+NFdE7U1iDOqz+eG\nlH5/U65d3NRq23PmyN/2nXW458NtH6dSR2bVZb9qhbyLe38bz+8gzrGUGYTvp/STv5hu3s/9B4kf\nyhGNpL9qHaTLUQsDtp9VY21N+Z2dVbsAH9wyTv39vA7Y2PZLJW0uaUvbX2kZh5Lwn0j5khGl7HkX\n9bh+RBmzPg+wpI1bj5KzfWy9UL4L5W/gX1omsiH3SHr4YL5NHVLaxfWYzj6fNY/dZvs227+Q9Ftg\nb+Dhkj5i+0/L3MDyxplDR/4vA14KrGP7YfWb+Sjbu0/y0hWJtQGL/4i/28UfsWq5AEmXuk4ZV8MS\nAkNxvki5njFcQmCnQSJtFKPTUQtDcZ7B4iGra1NrvNRrAS3jnEz5nb3Q9rb1y+A8249oHKfTMhJ1\neyMrUyBpIxbXXRqM9jmncYzdKYl5MF9lU8rciG81jtPZ51OlHMreNfHvQLlWdiTwCOBPtpt0Yc2Z\nI3/gFZSF4M8HsP2jOua/KUlPpPzhDmrfbCFpi9Z/xMAf64W4H6usdfALoPkROfCPwNuBz9f759a2\nljodtTBkT0p3zwdt/6Lxtoc9zPa+kvYHsP07NVwTWCMqI1GNpIyIpHdTLsJfztCXDOX/qxnbZ9YD\nvy0pv7OrbN81yctWRJefz9WG/n5fQFnf5H0qczKadfnNpeR/l+27Bh/C2ifbxWnNG4a2uxrlC+dC\n4MmN47yG8kE/FPg3ygf/wMYxBnVCmkxMW1YYOhy1sCiI/Yo60ufRNWkutH1L6zjAXZJWH9ypXQst\nE8woy0iMpIwI8CzKl0wXiXi8R1IGFswDdpCE7RMmec1Udfn5HD6S2J16kGT7Ly0PMuZS8j9b0lso\nR0t/R6mI2EW/8tOH70t6CKV0Qes4gxmWvwVeVI8y9qee2bQiaUvKkMhNWbIMcssvs05HLQxI2hd4\nL2Wil4CPSHqD7c81DjUfOA14sKTPAo+nJOsmbB8v6b+B/W1/ptV2h0l6fb35E2CBpK/QbZmCayij\nvTpN/vX39lDKmeXwGUbT5N/x5/MslVphN1IOlL4Fiy7+Nvv9zaU+/1WBgyhHTVCWhDy6i/Hk4+IK\nuNz21o22txalC2sj4EvAN+v91wOXtJ7kJekSylC4i1j8YfFgSFvDOJ2NWhiKcQnwlMHRfh2zfmbr\nkR512+tSiglCmRHdvICYpAtt79R6u3Xb8xlXmmD4cbcvU/B5Sp/1mSxZd6n1DN8rgG26+tyP4vNZ\nu3f2Ax4EnGz7htq+I6VUxTemGwPmUPIfFS05XX0VylTyn7pRLRxJpwK/Bv6Pcsq3PuXDeajt77eI\nMS5eZwlmXJxHcO+LfZ9f1mtWIMalwPaDD/6gj7TVxcs69HL4AzO+VEHTgmv1QvmvKBMKF5WrcKOS\nvjXGvrZPnqytQZwXTdBs28c3jvM54NVdXfMZ9eezS7M++WtEqx4NxXsRixPAPcC1tr/TcPvDowdW\npZz6beKOapPUI8BfUi74LjqlbJxgjgW2owxdXDR71PaLW8Wocd5LObr8LOUDuR/laOyNjba/gGVc\nq7D9pBZxhuJdO0G8pmPJJV3susDOstpmi/p/tAOwkCXPMJqcMY/68zku9jcplUo/2mJY8VxI/ptT\nvn3H16R5CHCj29XXH455P0o518FogmYTfMZ/8Lr+IC4lwWC7Sc3wGuNyYNsRdcE9m1KozpRhnl9o\nuP2dgets31jvH0iZRfwzyoSlTkfMtCTpqZSql/tRziwG1qR0m+zcKM7nbD93KQdpXRycjQ1vHxZV\n2zy70fZH+vkcF3sjYANKPaaPTnt7cyD5fxU4zPal49q3B95p+xmN440Bx1M+8AAbU8ohtPrj+jPw\n+6Gm1Skzl6Hl+p33TmQvoiSya2mcyFRqn7/PHa6nOgoqxc92dyl+thslab4S2BHYynbTEg8qC98c\nzNDkOBrVqandcDsC76aMVhHlTPZm4Cw3WsZT0oZ1vPqmEz3uxuW2h+KuS/m9/azl9atRfT5HYS4k\n/wu8lPKzkn7o9mV2LwKeZ/uqen8L4H/cePJV10aZyOrciFMpiaXLRVaeQ5mBPeiHHcRp9YX5A9eJ\nXJI+CvzS9vzxj7WiUqtoHuVgY1Cn5h43mOSjUrnznZSJkYMJURtTJki9ueXZbI33EsrM7qtbbndo\n+1+lzBr+YZ0ncTHwPeBhwCdtf6CLuF0YVVf2XBjquaziUE0XcKjmDRI/LJpMNht/j6sM9evvB3zc\n9inAKZKa1o6hrHN8APBDGlaMnMB7gKe7g9o01aqS7lMT41OAlw091sXfwKPHfdDPrCOaWngvpdrp\nxrZ/C6CytsP7gP8AXt0ozsDGwMclbUap7Doot93qIummXlyM8MWUNQleqFLs8Txg1iR/ykL3UIar\nw5Kz75uZjUlrvAskvcz2J4YbJb2UxWWYW7pQ0tHAf7P4P+SCDuJ0bZSJ7Bbbpzbe5kRu6jDxQ1lb\n9WxJv6Kc+p8Li647NV9li27r1Dwd2ML28AX4OyS9nFJDqGnyt/2vACqT414GvBH4IO3qVQ2fqTwF\n+GSN+1tJXR5wNDfoCpO0h+0dhh66pJ6x/0uLOHMh+b8G+IKk57M42e8E3I8yq7C1gynjegfjk88F\nPtZBnK6NMpFdXCdDfZklJxI1HepJORA4CfhiF3Fsv1PStyjjr08fSpyim1nSbwC+JWmJOjWNtv2X\n4cQ/YPvPXSRLSW+jrIWxBmUC1utpW8//epU6RTdQui5Pq3Hvz+zNc5K0q+1v1zuPZ8nZv9Pb+Gzv\n84dFozyeBPwt5cLYZW5cyGkukvRYFiey39W2LYA1Wo5ZV6nlM9GIotZDPUcSZ5QkrUYZWQYN69RI\n+hLw+fHj7CUdADy31dDIoe1eTDk6/yqly+e8Vu+lbn994AjK3/NHbZ9e259EKVT4H61ijUqdV3Is\nsFZtuo1SpK7JZ3NOJP9RGPV8gpg6SX8zfpSSpIe6o5r+XatdJIcwNHSVUqn2jw22/WDK3I4/sOQZ\n8/2BZ7mD5TzrNYXHU6quPhe42faureOMgkqZ5fFup1xkfn3LvzmVWcXYbrqQT5L/cpqJ+QSznUqZ\n7QX1orgoF34Hw0lf1PLsosY7D3jq4EOisuD552xv2zLOqKjMVr2DxdeX/gFYy/ZzG21flIKE21K+\nXC63fWaLbU8QaztK0t8NeBTlc3TO4FrAbCPpHcB1lO5TKHV9HkYZZfRy22MNYqxG+bxsypJ1t46Y\n7rYhyX+5jXo+wVwg6TJKKee7Jf0DpYDc31H6ZA+3/YTG8Z5GuRj2/yglfU8Ant9wRMlISbrc9jaT\ntc0G9fNzDuXs5QI3WpBkpmjiFcO+b3uHVsN+JX2D0tVzIUuutdBk/eDZeiFkJqw/PvED2L6kDl+L\ne7t7aLxRMV0rAAAPRUlEQVT404ETarfMN1VKMTRl+6t1YtQZlAuLzx4eljsLXSTpsbb/D0BlDd8u\nRrB1zvbTBrcl/bWkrWy3GrY6E34vaT9gUDF2H2DQHdfqiHoj23/faFv3kuS//EY9n2Au+ItKNc9f\nU4pgHTn02OoTv2TqtGSxPSi11a8BXqmyUHzTypEj9CjgO5KuoySUjYGr6vWnWXWdSaXmzl6UnHMh\n8EtJ37H92sZxPszisg7U23cA37P9pYahnk8p5T4os3A+8IJ6neaVjWKcJ2n7rr4kk/yX36jnE8wF\n/0q5ADYPOHUwCUelRMY1DeNcyJLliQf371WqeJbZc6Z3oKG16zyCl1DOAA+fZBDFilqN0uX3Ocr/\n/3MoM5i3l/Qk269pEcT2NZSz2Ym0GsL6BODFdahv85nx6fNfTiorRH2BMn78XvMJXGvkxJJqGYE1\nh2YTo7LmrWzfOXN7NjtIegLwcJcF0NejDMP96WSvW9nURL8HpVTFW20vnKjfvEGc7wKPt31PvT+P\nkox3BS51u3U3HkgpjbEpS16MbbYEqjquh5Qj/+Vk+yZJj2PJ+QRfyXyCZat9/uPLQ29DmYzTNPlL\n2pWyGPmmLPmBbFYCeZRUym3vRDmSPZayEtZ/U4ZLzjZHUBZY+k5N/A8Duqjzszbles9gouIawDq2\n75E07SGyQ75EuYB9BotLljQ9kh6a6ftAOuhazpF/jJykEyj1/X9ke7+G272KMuN7eFUy3MEqW6NQ\nayztCFzoWja4i6PluUTSQcBbKRVQoSwZeiRljYf5tt/QKM73x5VeaE7SXpRaSxsCtwCbAFe0Grq8\nSouNREyF7RfWZPbSxpu+zfbXbd9s+1eDn8YxRumu4RIMtbtsVpL0HkkPkHQfSWdK+lWdTdyU7WMo\nZ0ZfpHTTPt72J23/rlXir75ShxZ36R3AYykHSZtRBk18t9XGc+QfndO9lz+EMhvyZ4O+2UZx3kUp\nFPZ5lqzt03Qy2ahIegPwcEpf+b8D/wh81vZ/zuiOrYDB2HdJz6JcKH0dpapnF+srb8S9lww9p3GM\nOymzof/E4qJydsN6/qpLrNYzwEe61F1qduaXPv8YhY9S+q4HQ9YGSzquJelgN1qQGtil/jt+fYem\nyyuOiu33StoD+C2lvs/bbJ8xw7u1oga55unA/9q+XVLzI09J76aUKL+coa4/Sv98M7bXaLm9pfiN\nSknqc4HPSLqFhtfJcuQfnZP0eUriuqze34ayetQbKcXFpjUbUtLrxzWZsvD5t2dxXZ9Vgb8edFup\nLB16IPDaViNWRqmele1NmQi1M+XC7Jdt77LMF049zo+A7dywaNy47W9t+wpJEy7e1PIss3bz/ZHS\nPf98yvyVz4yvX7WicuQfo7Clh5ZwtH15neF5TaOjvzW5d7fSJsBbJM23feIEr1lpSdof+DhlFumP\nKCtufYqybkTTBT1GxfZhkt4D3F67L35H+TJo7RrKqKhOkj+lu+qlwPuZeHRPs7NML660+9eUiWqX\ntkr8kCP/GAFJJwO3Av9DmXizL7Ae8ALK0fmjO4q7DnCmR7TAdiu1JtIzbf+4HmGeDzzH9pdneNdW\nmJZcjxhgAY3WIx4X5/PAI4AzWXJi1KyZ5a0RLUmZ5B+dU1lQ4xAWj0//DmUBnD8Cf+W6jGBHsS+e\nhcl/iX1WB2tRj5o6XI94XJwXTdBsj1u3oEGcnwDvtX3UUNtXbC9t1u9Utn3ZYDinpDdT1tRetCSl\n7e2mGwPS7RMjYPv3tebK4MLuVV5c1bHLxP8k4Dddbb9D60l6HYvr06w9dN+23z9zu7bCulyPeBHb\nx7Xe5lLcDYxJ2plSwvkuYKOG2x7obEnKJP/oXK3lczzws9q0saQDbZ+99FdNafsT1Yj5a+BG4IUt\nYozY0ZTrGEu7Pxt1uR4xkj5n+7lL+Vvoogje723vJ+mNwDmS9m247ZEsSZlun+icpIuA57mWV1ZZ\nKvJ/bE84YmIFtr/puCYDt6Z20MpD0u6UEhVLrEfcqjyKpA1t/6LrejhD8RZ1zUl6CmU48zq212uw\n7ZEsSZnkH52baGJKyhT0Rx22+mrKdZ4ta/NVbrAc5QSxXgKcbbuLukHDcfayferQ/U2AA91ola1R\nSPKPzkk6ljLhZrAc4fOBVVpWQIyVm6TvdTWqa1ycIygVPDejDI09hzKTuOlqbpIGY+83s32EpI2B\nDWw3K7/QtST/6JzKWqSvYPFon3OBj3U1EWeu0ASLz0/UNhtI+gBwH+Ak4HcsvnjdSekNlUVVXkZZ\nOnRD26s23v5/Uap5Ptn2VnVY8em2x88uX2kl+UespCYapjqo9zJT+7SiVFbyuleysd209IaktwGP\no5Ry/j7lQOPbtn/ROM7Ftncc1/ffZO3eUclon+icpGdQLmBtypJ19psVwZpLJG1NWfNgLUnPZvGK\nZA9gli4ZantsRKGeTRkqOVgw/ryOzjD/VK9lAKCy0E6zYZh1m50uSZnkH6PwQeBZwA+HSxTHUm0B\nPANYq/478Fval8HuVK27dLvto8e1H0RZ4e2DLePVo/EHULoY/w74hKSbbe/aMg7wYUrJ6AdKOpKy\ngPtbG8fodEnKdPtE5+op/+62/zzZc2MxSY+zfd5M78d01GG+jxma1Ddovy9lkZoms1WHtrsdZe3b\n3SjVXa8HzrH9ry3j1FhbU2rsQykjckXj7Xe6JGWO/GMU3gh8rX4JDNfZn40zVUfpOklfoHzYoXRj\nvNr29TO4T1M1b3ziB7D9J0ma6AXT9O+Ufv7/pHSPNK0dNM66wO9c11eWtJnbrq/c6ZKUSf4xCu+k\ndFmsRqm4GMvnWOAzlEJ4UIYWHkvpzpgtJOlBtm8a17g+jde8BbD99Fr+egtgS0lXdfEFoNGsr/we\n4GJJSyxJWUs9f3O6G0+3T3RuLhQmmwkTjR6ZbSNKJL2QMsHr9cCFtflRwHuBj7SuxTNRKRHK5Ksm\npUSG4oxkfWVJG1LWPzDlTKbZqKUc+ccofE3S37vdil19cavKOrefpVzw25+ySM2sYfsESb+kjPYa\nLDx+GWVxn693EPL9wB7jS4kATUqJDLnL9l8GPVfqbn1lAb+k5OqH1/pITVYly5F/dE4jWO90Lqp1\naj4MPKY2nQe8yvbPZ2qfVnajKiWiEayvrKUsSWn7GUt90VS2n+QfEXPFKEqJ1AvVDwG2oiR/gG+4\n8frK6npJyiT/6Eo9cr3N9m31/pMpS/ddS+nvvdcokABJhy/lIQPMpuJhozaKUiI1+V/a9XUsSV8H\n9nVHix0l+UdnJC0E9q6ldnegLK13JGWZvT+58SpOc4Wkf+beI2H+CjgIWNd2V/3LsZwkHU8pt7yw\nwxidLkmZ5B+dGe5rlfQfwF9sv7FWRPxB6wk+c1GdrXooJfGfDLzP9i0zu1dTJ+lBlCG/G9neU9I2\nwGNtH9No+xMt4jLQfDEXSVdR+vx/RilU1zyOOl6SMqN9okvDk3h2B94EMDxKIiYm6W+A11L6rE8A\nHml7Ni5JOXAcZTz8W+r9qylfZk2SP0uWweiMpI3rBfe/Z8m6O821HgY7XpJ/dOksSZ+jLKe4NvAt\nWDR2OeWcl6KeJT0L+ASwfVd9viO2ru2TJB0GYPtuSc2WcfQEK3VJWpeyolvL7o0vATvavlbSKbaf\n03DbwOiWpEzyjy69hjJU7UHArkMXeNdn8RFg3NvrKMNi3wq8ddxZ0mwdIntnPZsBQNJjgNtbbVzS\nYylDLn8NvINytrQusKqkF3Y0p+ChHWwTyqQ46PhsJsk/OlMreJ44QfvFM7A7s4btVWZ6HzrweuDL\nwEMlnQesR6mE2cpHKN2Ka1HOMPe0fb6krSiTvLpI/p0YmsX7FDpckjLJP0ZO0jcpk70+avsrM70/\n0a1a9363+rMVpZ/8qsZDfVf14oXOj7B9PoDtKyW17PbZXtKgG271ods1XNOzso2Bj0vqZEnKjPaJ\nkZO0EbABsIvtj870/kT3ul7Dd9yKWkusgDbRimizSVdLUib5R0Tnul7DV9Kfgd/Xu6sDfxh6eHXb\ns66Xo+slKZP8o3N11ML4YXG3A98D3mH71hnZsRiZUa3hO5dIupgOl6RM8o/OSXovcA9LVqe8P3AT\nZaWikYzRjphthpakfALwXKDZkpSz7lQoZqWnjOtzvWTQDzvJzMyYI7qe4TsXLW1Jylbbn4tDymLl\ns6qkXQZ3JO3M4r+9ZhN9YqV2HHA6sGG9fzVlBnMs3b8Da1KWpNza9ljLtYjT7ROdk/RoytT+NWrT\nbym1ai4Dnmb75JnatxgNSRfYftS4UTnft73DTO/bymxoSUpThsc2W5Iy3T7ROdvfA/5W0lr1/vDM\nziT+fuh0hu9cNNGSlJKaLUmZI//onKS1gcMpfZcAC4Ajxn0JxBwmaSfKqmTbUs741gP2sf2DGd2x\nlZiki4DnjV+S0naTJSmT/KNztS75pZSjGAEHUAqWPXtGdyxGStI8upvhO+d0vSRlkn90TtIPbD9i\nsraYuyRdQqmxc5Lta2Z6f2aDrpekzGifGIU/SHrC4I6kXVk8GzP6YS9KIjtZ0gWS/lnSxjO9Uyu5\ng4ErKIv5vIrSXXZwq43nyD86V5dwPIFScRHgN8CB6e/tJ0mbA28Dnt+qTk1MXZJ/dGZo1aPB/YlG\n+0RPSNqUsr7DvpSzgJNsv28m92llNKolKZP8ozPjxnR3supRzA6SvgvclzK09yTbP5nhXVpp1S/J\npZpo1bIVkXH+MSpdrXoUs8OBtq+c6Z2YDUa1JGUu+EbEKNwm6RhJpwFI2kbSQTO9UysjSY+VtEDS\n5yU9UtIPgR8Ct0h6arM46faJrkxSY322rkUbK6Am/WOBt9jeXtJ9gItt/+0M79pKR9KFLF6S8pOM\nW5KyVUmMHPlHZ2yvanvN+jNv6PaaSfy9s67tkygXeqk1alLUb2Kr2j7d9ueAG4eXpGSCNRFWVJJ/\nRIxCavssv+EE/8eugqTbJyI6l9o+y29US1Im+UfESNR+/i3r3abliWPq0u0TEZ2RtLOkDWBRP/9O\nwJHA+yStM6M713NJ/hHRpY8DdwFI2g14F6W66x3AJ2Zwv3ovk7wiokur2P51vb0f8HHbpwCnSEp/\n/wzKkX9EdGnV2tcP8BTgrKHHcvA5g/LLj4gunQicLelXlBEs58Kiyp63zeSO9V1G+0REpyQ9FngQ\ncLrt39W2LYA1bF80ozvXY0n+ERE9lD7/iIgeSvKPiOihJP+IiB5K8o+I6KH/D7DzrU4AEuXgAAAA\nAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "USAccidents_Data.WeatherVal.value_counts().plot(kind='bar')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "** Weather has no significant impact on accidents. **" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3. Time of the day ###" ] }, { "cell_type": "code", "execution_count": 63, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 63, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEBCAYAAABysL6vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG/hJREFUeJzt3X+8HXV95/HXm4TwQ5GIaEwC3UQfN4WrKEIhdFubQDVN\nW0vSH5LQihFSXDZbRNdHa6I+zN32sZha20p3G6zFQKg1boqSR9KFmMjmWFwlwYoYuKQk1gA3bq7a\nRamufZjIZ/+Y7yXHw5lzcueee3Lu976fj8c87sx35jvf78w95zNzPmfOjCICMzPLz0knugNmZjY+\nHODNzDLlAG9mlikHeDOzTDnAm5llygHezCxTLQO8pA2ShiXtrSu7UNIDkh6S9KCkS+rmrZG0X9I+\nSYvqyi+WtDfNu2V8NsXMzOq1O4O/HVjcUPYhYG1EvA74QJpGUj+wDOhPddZLUqpzK7AyIvqAPkmN\n6zQzsw5rGeAj4n7g6YbiZ4Ez0/h04FAaXwJsiogjEXEQOADMlzQTOCMi9qTl7gSWdqDvZmbWwtQK\ndd4JfFbShykOED+bymcBD9QtNwTMBo6k8RGHUrmZmY2jKgF+FfDOiLhb0puBDcAbO9EZSb5vgplZ\nBRGhxrIqV9G8NSLuTuN3AZem8UPAuXXLnUNx5n4ojdeXH6JERDQd1q5dWzrvRNdx/7wvTnRbvd4/\n74vxbatMlQD/TUkL0vgVwONpfCuwXNI0SXOBPmBPRBwGnpE0P33peg2wpUK7ZmY2Ci1TNJI2AQuA\nsyU9RXHVzPXALZKmAj8E3g4QEYOSNgODwFFgVRw7tKwC7gBOA+6JiO3jsC1mZlanZYCPiKtLZv1M\nyfI3Azc3Kf9H4IJR967OwoULe7ZON9vq9f51s61e71832+r1/nWzrV7vXzfbUqv8TbdJil7qj5nZ\nRCCJ6NCXrGZmNgE4wJuZZcoB3swsUw7wZmaZqvJL1q44dp+y5vxlrJlZaz0b4AtlQbx18DczM6do\nzMyy5QBvZpapHk/RjF6r3L3z9mY2mWQX4AvNArnz9mY2uThFY2aWKQd4M7NMOcCbmWXKAd7MLFMO\n8GZmmWoZ4CVtkDQsaW9D+Y2SHpP0iKQ/ritfI2m/pH2SFtWVXyxpb5p3S+c3Y2wktRzMzCaidmfw\ntwOL6wskXQ5cCbwmIl4NfDiV9wPLgP5UZ72ORcdbgZUR0Qf0SfqJdfaGKBnMzCamlgE+Iu4Hnm4o\n/o/AByPiSFrm26l8CbApIo5ExEHgADBf0kzgjIjYk5a7E1jaof6fMD7rN7NeVyUH3wf8gqQHJNUk\njTyfdRYwVLfcEDC7SfmhVJ4Bn/WbWe+q8kvWqcCLI+IySZcAm4FXdKpDAwMDdVM1YGGnVm1mloVa\nrUatVmu7XNuHbkuaA2yLiAvS9L3Auoj4fJo+AFwG/C5ARKxL5duBtcATwK6IOD+VXw0siIgbmrT1\n3EO3izRH+e2Cy/pdXq9KnfJ6VftnZtZpnXzo9hbgirTSecC0iPgOsBVYLmmapLkUqZw9EXEYeEbS\n/PSl6zVpHWZmNo5apmgkbQIWAC+R9BTwAWADsCFdOvkj4K0AETEoaTMwCBwFVsWx09hVwB3AacA9\nEbF9HLZlQvDdLs2sW9qmaLppMqRoqvevXC/9D82s+8pSNJneLjhHfnyhmY2Ob1VgZpYpB3gzs0w5\nwJuZZcoB3swsUw7wZmaZcoA3M8uUA7yZWaYc4M3MMuUfOmXMt0Uwm9wc4LPX/LYIZpY/p2jMzDLl\nM3j7Cb6xmVk+HOCtCd/YzCwHTtGYmWXKAd7MLFMtA7ykDZKG09ObGue9W9Kzks6qK1sjab+kfZIW\n1ZVfLGlvmndLZzfBzMyaaXcGfzuwuLFQ0rnAGykeqD1S1g8sA/pTnfU69o3drcDKiOgD+iQ9b51m\nZtZZLQN8RNwPPN1k1p8Bf9BQtgTYFBFHIuIgcACYL2kmcEZE7EnL3QksHVOvzcysrVHn4CUtAYYi\n4msNs2YBQ3XTQ8DsJuWHUrmZmY2jUV0mKel04L0U6ZnnijvZoYGBgbqpGrCwk6s3M5vwarUatVqt\n7XJq98MVSXOAbRFxgaQLgM8B/y/NPofijHw+cC1ARKxL9bYDayny9Lsi4vxUfjWwICJuaNJWjPSn\nSN+XX49d1u/yelXqlNfr9f51uy0zO3EkERHPO9keVYomIvZGxIyImBsRcylSLxdFxDCwFVguaZqk\nuUAfsCciDgPPSJqfvnS9Btgy5i0yM7OW2l0muQn4IjBP0lOSrm1Y5LnTuYgYBDYDg8C9wKo4drq3\nCrgN2A8ciIjtHeq/mZmVaJui6SanaJyiMbPR60iKxszMJg4HeDOzTDnAm5llyrcLtjHzPeTNepMD\nvHWI7yFv1mucojEzy5QDvJlZphzgzcwy5Ry8nTCtvpz1F7NmY+cAbydY81/alvEVO2bHzwHeJiBf\nsWN2PJyDNzPLlAO8mVmmHODNzDLlAG9mlikHeDOzTLV7otMGScOS9taV/YmkxyQ9LOkzks6sm7dG\n0n5J+yQtqiu/WNLeNO+W8dkUMzOr1+4M/nZgcUPZDuBVEfFa4HFgDYCkfmAZ0J/qrNexi5ZvBVZG\nRB/QJ6lxnWZm1mEtA3xE3A883VC2MyKeTZO7gXPS+BJgU0QciYiDwAFgvqSZwBkRsSctdyewtEP9\nNzOzEmPNwV8H3JPGZwFDdfOGgNlNyg+lcjMzG0eVf8kq6X3AjyLikx3sDwMDA3VTNWBhJ1dvZjbh\n1Wo1arVa2+XU7t4dkuYA2yLigrqytwHXA78YEf+WylYDRMS6NL0dWAs8AeyKiPNT+dXAgoi4oUlb\nMdKfIn1f/pP0sn6X16tSp7xer/evm21NjP6Va93W6OuZdZskIuJ5L9hRp2jSF6S/DywZCe7JVmC5\npGmS5gJ9wJ6IOAw8I2l++tL1GmBLpa0wqyxKhvGqZ3bitUzRSNoELADOlvQUxRn5GmAasDOd4Xwp\nIlZFxKCkzcAgcBRYFcdOcVYBdwCnAfdExPbx2BgzMzumbYqmm5yicYpmYvSvXVvleun9ZvkoS9H4\ndsFmHefbGVtv8K0KzMwy5QBvZpYpB3gzs0w5wJuZZcoB3swsUw7wZmaZcoA3M8uUr4M36wGtfiDl\nH0dZVQ7wZj2j+a96y1T51ax/aTu5OMCbTWhVfjXrX9pOFs7Bm5llygHezCxTDvBmZplygDczy1TL\nAC9pg6RhSXvrys6StFPS45J2SJpeN2+NpP2S9klaVFd+saS9ad4t47MpZjZeJLUcrDe1O4O/HVjc\nULYa2BkR84D70jSS+oFlQH+qs17H/vO3Aisjog/oS4/9M7MJxY8vnGhaBviIuB94uqH4SmBjGt8I\nLE3jS4BNEXEkIg4CB4D5kmYCZ0TEnrTcnXV1zMxsnFTJwc+IiOE0PgzMSOOzgKG65YaA2U3KD6Vy\nMzMbR2P6oVNEhKSOfkYbGBiom6oBCzu5ejOzCa9Wq1Gr1dou1/ah25LmANsi4oI0vQ9YGBGHU/pl\nV0ScJ2k1QESsS8ttB9YCT6Rlzk/lVwMLIuKGJm35odsd6l8325p8/etmW73ev9b1rDvKHrpdJUWz\nFViRxlcAW+rKl0uaJmku0AfsiYjDwDOS5qcvXa+pq2NmmfKVNydeyxSNpE3AAuBsSU8BHwDWAZsl\nrQQOAlcBRMSgpM3AIHAUWBXHDuurgDuA04B7ImJ75zfFzHqP73tzIrVN0XSTUzRO0UyM/nWzrV7v\nX+fbstHrZIrGzMwmAAd4M7NMOcCbmWXKAd7MLFMO8GZmmfIj+8ysp/gB5J3jAG9mPWh0DyC35pyi\nMTPLlAO8mVmmnKIxswmv3b1tJmvu3gHezDLh+940corGzCxTPoM3s0lpMqR1HODNbBLLO63jFI2Z\nWaYc4M3MMlU5wEt6l6RHJO2V9ElJp0g6S9JOSY9L2iFpet3yayTtl7RP0qLOdN/MzMpUCvCSZgM3\nAhenh3FPAZYDq4GdETEPuC9NI6kfWAb0A4uB9ZL86cHMbByNJchOBU6XNBU4HfgmcCWwMc3fCCxN\n40uATRFxJCIOAgeAS8fQtpnZCTGRHiReKcBHxCHgT4EnKQL7dyNiJzAjIobTYsPAjDQ+CxiqW8UQ\nMLtSj83MTrhoMvSeSpdJSnoxxdn6HOB7wN9Jekv9MhERklptddN5AwMDdVM1YGGVLpqZZatWq1Gr\n1doupyoX80t6M/BLEfG7afoa4DLgCuDyiDgsaSawKyLOk7QaICLWpeW3A2sjYnfDemOkP5Pv6fSd\n7V8325p8/etmW73ev262NXH7N94kERHPyxFVzcE/AVwm6TQVW/sGYBDYBqxIy6wAtqTxrcBySdMk\nzQX6gD0V2zYzs+NQKUUTEXsk3QV8BTia/n4MOAPYLGklcBC4Ki0/KGkzxUHgKLAqTtShzsxskqiU\nohkvTtE4RTMx+tfNtnq9f91sa6L3r9xY43BZisb3ojEz64ru3/fGPzYyM8uUA7yZWaYc4M3MMuUA\nb2aWKQd4M7NMOcCbmWXKAd7MLFMO8GZmmXKANzPLlAO8mVmmHODNzDLlAG9mlikHeDOzTDnAm5ll\nygHezCxTlQO8pOmS7pL0mKRBSfMlnSVpp6THJe2QNL1u+TWS9kvaJ2lRZ7pvZmZlxnIGfwtwT0Sc\nD7wG2AesBnZGxDzgvjSNpH5gGdAPLAbWS/KnBzOzcVQpyEo6E3h9RGwAiIijEfE94EpgY1psI7A0\njS8BNkXEkYg4CBwALh1Lx83MrLWqZ9FzgW9Lul3SVyT9taQXADMiYjgtMwzMSOOzgKG6+kPA7Ipt\nm5nZcaj6TNapwEXA70XEg5I+QkrHjIiIkNTqSbJN5w0MDNRN1YCFFbtoZpanWq1GrVZru5yqPM1b\n0suBL0XE3DT988Aa4BXA5RFxWNJMYFdEnCdpNUBErEvLbwfWRsTuhvXGSH/8RPaJ09bk61832+r1\n/nWzrcnVv9GQREQ87+ndlVI0EXEYeErSvFT0BuBRYBuwIpWtALak8a3AcknTJM0F+oA9Vdo2M7Pj\nUzVFA3Aj8LeSpgFfB64FpgCbJa0EDgJXAUTEoKTNwCBwFFgVYz1kmZlZS5VSNOPFKRqnaCZG/7rZ\nVq/3r5ttTa7+jUZHUzRmZtb7HODNzDLlAG9mlikHeDOzTDnAm5llygHezCxTDvBmZplygDczy5QD\nvJlZphzgzcwy5QBvZpYpB3gzs0w5wJuZZcoB3swsUw7wZmaZcoA3M8vUmAK8pCmSHpK0LU2fJWmn\npMcl7ZA0vW7ZNZL2S9onadFYO25mZq2N9Qz+JorH8I08jmQ1sDMi5gH3pWkk9QPLgH5gMbBekj89\nmJmNo8pBVtI5wK8AtwEjj4q6EtiYxjcCS9P4EmBTRByJiIPAAeDSqm2bmVl7YzmL/nPg94Fn68pm\nRMRwGh8GZqTxWcBQ3XJDwOwxtG1mZm1MrVJJ0puAb0XEQ5IWNlsmIkJSqyfJNp03MDBQN1UDmq7e\nzGzSqtVq1Gq1tsupytO8Jd0MXAMcBU4FXgR8BrgEWBgRhyXNBHZFxHmSVgNExLpUfzuwNiJ2N6w3\nRvrT60887/X+dbOtyde/brbV6/3rZluTq3+jIYmIUGN5pRRNRLw3Is6NiLnAcuB/RcQ1wFZgRVps\nBbAljW8FlkuaJmku0AfsqdK2mZkdn0opmiZGDj/rgM2SVgIHgasAImJQ0maKK26OAqtirIcsMzNr\nqVKKZrw4ReMUzcToXzfb6vX+dbOtydW/0ehoisbMzHqfA7yZWaYc4M3MMuUAb2aWKQd4M7NMOcCb\nmWXKAd7MLFMO8GZmmXKANzPLlAO8mVmmHODNzDLlAG9mlikHeDOzTDnAm5llygHezCxTDvBmZpmq\nFOAlnStpl6RHJT0i6R2p/CxJOyU9LmmHpOl1ddZI2i9pn6RFndoAMzNrruoZ/BHgXRHxKuAy4D9J\nOh9YDeyMiHnAfWkaSf3AMqAfWAysl+RPD2ZmLUgqHY5H1YduH46Ir6bx7wOPAbOBK4GNabGNwNI0\nvgTYFBFHIuIgcAC4tErbZmaTSzQZjs+Yz6IlzQFeB+wGZkTEcJo1DMxI47OAobpqQxQHBDMzGydT\nx1JZ0guBTwM3RcS/1n9siIiQ1OpQ03TewMBA3VQNWDiWLpqZZeknY2Vzqvo0b0knA38P3BsRH0ll\n+4CFEXFY0kxgV0ScJ2k1QESsS8ttB9ZGxO6GdcZIf3r9iee93r9utjX5+tfNtnq9f91sa3L1bzRt\nSSIinpeYr3oVjYCPA4MjwT3ZCqxI4yuALXXlyyVNkzQX6AP2VGnbzMyOT9UUzc8BbwG+JumhVLYG\nWAdslrQSOAhcBRARg5I2A4PAUWBVlB2yzMysIyqnaMaDUzRO0UyM/nWzrV7vXzfbmlz9G01bHU3R\nmJlZ73OANzPLlAO8mVmmHODNzDLlAG9mlikHeDOzTDnAm5llygHezCxTDvBmZplygDczy5QDvJlZ\nphzgzcwy5QBvZpYpB3gzs0w5wJuZZaqrAV7SYkn7JO2X9J7R1a5VaLFbdbrZVpU6ubZVpU6ubVWp\nk2tbVerk2VbXArykKcB/BxYD/cDVks4//jXUKrTarTrdbKtKnVzbqlIn17aq1Mm1rSp18myrm2fw\nlwIHIuJgRBwBPgUs6WL7ZmaTSjcD/GzgqbrpoVRmZmbjoGvPZJX0m8DiiLg+Tb8FmB8RN9Yt0zsP\niDUzm0CaPZN1ahfbPwScWzd9LsVZ/HOaddDMzKrpZormy0CfpDmSpgHLgK1dbN/MbFLp2hl8RByV\n9HvAZ4EpwMcj4rFutW9mNtl0LQdvZmbd1c0cvNmEkH6fMQvYHRHfrytfHBHb29RbwrGrw4aArWWf\nVCVdBjwWEd+TdDqwGrgIeBS4OSK+15ENqkjSO4C7I+Kptgt3vu3XU1xavTcidpQscwqwHDgUEZ+T\n9DvAvwcGgY+ly7HL1j8feDYiHpT0Korf5zwWEfe0qPNK4DeAc4BngX8CPhkRz1TayC7wrQqakHRv\nF9p4yXi3cbwknSlpnaRPSPrthnnrR7mul3W2d92VgtoW4EbgUUlL62Z/sEW99wCb0uTuNJwEbJK0\npqTaBuAHafwW4EXAOuCHwO0V+n7taOu08UfAHklfkLRK0ks7vP7nSNpTN3498N+AFwJrW+y/24Ff\nAW6S9DfAbwEPUBwYbmvR1gDF/v6opA+mtk4HVkt6f0mdm4CPAqek9Z8C/BSwW9Llx7+lIOnO0Sw/\nJhHRUwNwCbAL+ATFlTY7ge8BDwKvK6nzEPB+4JWjaOeikuFi4HCLetMp3oT7gKeB/5vG1wHTS+r8\nMfDSNP4zwD8DB4AngYUldc4A/pDibO4Z4DsUQeNtbbZrcUNfPw7sBT4JzCip85nU/18HtgGfBk4d\n2bct2jqrYXgJcHBkuqTOTOBW4C/T8gOpf5uBmZ2qk+qdmbbrE8BvN8xbX1LnEeCFaXwOxcUB7zyO\nfbEfOLlJ+TSKH/g1q/NY3fhXGuY9XOG981TF99y9Ld5XJwGLKA5G3wa2AyuAMzr8HnmobvzLde+X\nFwCPlNTZm/5OBb4FTE3TGpnX4n88lSKo/ytwZio/raxeqjMljZ8OfD6N/xTw1RZtbaO4mGRb3fCD\nkfKSOpXe+82GXkzRrAc+QPEi+RLwLooX2BVp3s82qTM9DbskDVMEs/8REd9s0c6DwD+UzDuzRb3N\nwH3AQmA4IkLSTIoX/ebU10a/GhEj9975MLAsio+G8yjO+i5uUudvgbspPjq+meJs5lPA+yXNi4j3\nlvTvgxRvQoA/Bf4P8GsUwfuvgKVN6rwyIn4jjd8t6X3AfZLa/dL4O8ATDWWzgX8EAnhFkzp3AH+f\ntqdGsZ2/SpHa+CjNf91cpQ4UZ3iPUxywrku/xfidiPg3mr+OoPhe6vsAEXFQ0kLg05L+HUXgKPNj\nim0/2FA+K81r5lFJ10XEBuBhSZfUvS5+1LRz0t4WfSj99CTporJZwOvK6kXEs8AOYEe6+u2Xgasp\nXltnl1Sr8h6ZIums1J8pEfHt1P4PJB0t3yydQhFwT6N43/4LcCqt089HI+IocFTS1yOlwiLih5Ke\nLakTwMkU/8tTKQ48RMSTkk5u0dY5FCmj2yjSOqI4yftwizpV3/tNel3hiD+eAz95JH+yYV7TI+VI\nnbTzfoHibO8wxSeBt5fUeRSYVzKv9EwIeHy084DHSGd3wAMN88rOGL7WMP3l9Pck4J+Oc/89TPoi\nfWS6Rf9Oaih7W9pHT7Ro690UB5PX1JV9Ywz/37L+jbpOs3nA+4D/TRGYmp6Np9fMhQ1lJwN3UuRs\ny9paTPGpbDvw12nYDnwd+OWSOtOBjRSf6HYDR4BvUJx4vLakzjBFQJ7TZPhmi/79OG1bs+GH7fZ7\nk3kvGO37oNU8igPjN9Lwz6RPZhRns2Xv+3elZZ8EbqI4qNxGcbY90KIPu4HT0/hJdeXTafgkVTfv\nJopPjbdR5N6vS+UvA/6hRVtTgP8MfI6UgTiO90il937TdY1m4W4MFDm0XwKuojg7/PVUvmBkQ4/n\nhUhxBF8M3F5S583AeSXzlrbo307gD6hLdwAvB94DfK6kzo2p3hUU6YVb0vb8F+BvSup8CXh9Gl8C\nfLZuXqsAP5ReUO9Ob5r6AP+1kjp/AryxSfliYH+b/9e5wN8Bf06RQ2734n24bvy/NswrO9iNuk6a\nN+oDV9qelzcpF/DzbbZtCsUng98CfhO4jJQ2aFPvTOBCijO757XdsOyGkddFk3mbWtQb9QkN8NPt\n+l5Sb9TvkRbrOh2Y22L+bGB2Gn9xel9f2madp5aUnw1c0KLeq9P/tmncaNPmOel98pdl+7tu2Urv\n/abrqvIPHM8hvdB3UJz9nAf8BfBdio85P1dS51MV2zof+EVSzrWuvOkZV5p3FvAhjuUXn07jH6Ik\n75zqXU7x8fQhijOBe4H/QJO8bVr+tRRppO9SnHX+dCp/KfCOFu0MAGvrhpel8pnAnZ3cFw3LLaE4\nMxpus9wf0SR/C/QBd3WqTppf+cCV20DFE5qKbVV6j0yGAXgTxRVSrZapf+9/4Xjf+03XdaI3eJQ7\n59pO1QHeQfFRawvFJ4WldfNKP5qm+ecDb2gMOtR9wTnO23Rdxf3XtN5Y9kXDek4nnQF1ervGegDq\n1D7McejmvqjyupiMA/BKik9Bf0Hx6fgG4EWjXs+J3pBRbvSorxIoq0P1qyU6EgzHY5u6vS+6+L/q\n6D4fyz7McejmvvB+P659dBNFmuv9wBcpLi65mSLlePlo1tVzV9G0uUpgRoU6ZVcWVL1a4u3AxRHx\nfUlzUp05EfGRsgpV+ldlP4yhXqV90cXtatznd7Xb52NoK0vd3Bfe72N2PcWX7D+W9GcUl7EukPRR\niksuLzzeFfVcgKcIDIsp8naNvtjBOt+SdGFEfBUgBY83UVw3/poW/WsMhgtoHwy7tU1V61XdF93a\nrqoH46r7MEfd3Bfe72MTVLsk83l6McD/T4p0wUONMyR9voN13kpxWdpzIuKIpBXAx1r0r0ow7NY2\nVa1XdV90a7uqHoCq7sMcdXNfeL+PzW3Ag5J2A6+n+KHkyK/E/2U0K/LNxkZJ0rnAkYg43FAuiqt8\nvnBiepYv73ObbCS9muIqwkciYl/l9TjAm5nlyTcbMzPLlAO8mVmmHODNzDLlAG9mlqn/DxWAw4PP\nOIMZAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "USAccidents_Data.HOUR.value_counts().plot(kind='bar')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "** This chart shows most of the accidents happen during or after sunset. So sun light also has an impact on accidents. **" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 4. Day of the week ###" ] }, { "cell_type": "code", "execution_count": 52, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 52, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import HTML\n", "HTML('')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It is clear from the above chart that week end's are more prone to accidents than regular week hours. \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Conclusion:" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "From the above study I came up with the following conclusions.\n", "\n", "* There is a significant amount of reduction in crashes from 2006.\n", "* Alcohol related traffic accidents are still a huge issue \n", "* Most of the accidents happen in weekends\n", "* It is also interesting that most likely the accidents happen during or after sunset. \n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The study is not complete. There are a lot more variables to consider. I think including more aspects like hospital expenses, auto insurance expenses will be more interesting.\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.11" } }, "nbformat": 4, "nbformat_minor": 0 }