You need to sign in or sign up before continuing.

Exo 2 Final solution OK

parent 533d21cf
......@@ -25,7 +25,7 @@
},
{
"cell_type": "code",
"execution_count": 5,
"execution_count": 4,
"metadata": {},
"outputs": [
{
......@@ -990,7 +990,7 @@
"[1535 rows x 10 columns]"
]
},
"execution_count": 5,
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1002,7 +1002,7 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": 5,
"metadata": {},
"outputs": [
{
......@@ -1049,7 +1049,7 @@
"Index: []"
]
},
"execution_count": 6,
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1060,7 +1060,7 @@
},
{
"cell_type": "code",
"execution_count": 9,
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
......@@ -1076,7 +1076,7 @@
},
{
"cell_type": "code",
"execution_count": 11,
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
......@@ -1085,7 +1085,7 @@
},
{
"cell_type": "code",
"execution_count": 12,
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
......@@ -1098,16 +1098,16 @@
},
{
"cell_type": "code",
"execution_count": 13,
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fea89a6a898>"
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62e69e748>"
]
},
"execution_count": 13,
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
......@@ -1130,16 +1130,16 @@
},
{
"cell_type": "code",
"execution_count": 14,
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fea87933240>"
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62c59bcc0>"
]
},
"execution_count": 14,
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
......@@ -1160,6 +1160,210 @@
"sorted_data['inc'][-200:].plot()"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62c4dd898>"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEVCAYAAAALsCk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl8XGd97/9+Zp/Rvlqydi/xks2Jl+wkxEBMSptQAri0JJSUUJaWFmhL2t7b0l/TH3BvSeHXEhpugLAnhFwSlgRCQnZjx3Zix45jW95k2bK1SyONZn9+f5xlZjQzmtHMSJqRnvfrpZdHzzlnfGY0cz7nuwspJQqFQqFQxGNZ6BNQKBQKRfGhxEGhUCgUSShxUCgUCkUSShwUCoVCkYQSB4VCoVAkocRBoVAoFEkocVAoFApFEkocFAqFQpGEEgeFQqFQJGFb6BPIlfr6etnZ2bnQp6FQKBQlxZ49ewallA2Z9itZcejs7GT37t0LfRoKhUJRUgghTmWzn3IrKRQKhSIJJQ4KhUKhSEKJg0KhUCiSUOKgUCgUiiSUOCgUCoUiCSUOCoVCoUhCiYOiJIlGJdGommKoUMwVShwUJck//PQAH/3enoU+DYVi0VKyRXCKpU13v5cBb2ChT0OhWLQoy0FRknj9YcamQgt9GgrFokWJg6IkMcRBxR0UirlBiYOiJBn3h4hKmAiGF/pUFIpFSUZxEEK4hBC7hBD7hBAHhRCf19f/WQhxRgjxmv5zc9wxdwshuoUQh4UQN8WtbxRCvK5v+6oQQujrTiHEQ/r6TiFEZ+FfqmKxEI1KJgKaKIz5lGtJoZgLsrEcAsCNUspLgQ3ANiHElfq2e6WUG/SfXwIIIdYD24ELgW3A14QQVn3/+4C7gNX6zzZ9/U5gREq5CrgX+GL+L02xWJkIhpG6N0nFHRSKuSGjOEiNCf1Xu/4zk6P3FuBHUsqAlPIE0A1sEUI0A5VSyh1SSgl8B7g17pgH9cePAFsNq0KhmI7XH3MlKXFQKOaGrGIOQgirEOI1oB94Skq5U9/0SSHEfiHEN4UQNfpaC3A67vBefa1Ffzx9PeEYKWUYGAPqUpzHXUKI3UKI3QMDA1m9QMXiw+uPCYISB4VibshKHKSUESnlBqAVzQq4CM1FtBLN1dQH/Lu+e6o7fjnD+kzHTD+P+6WUm6SUmxoaMg4yUixS4i2HURVzUCjmhFllK0kpR4FngW1SyvO6aESBbwBb9N16gba4w1qBs/p6a4r1hGOEEDagChie1StRLBnGp5TloFDMNdlkKzUIIar1x27gbcCbegzB4N3AAf3x48B2PQOpCy3wvEtK2Qd4hRBX6vGE24HH4o65Q398G/CMHpdQKJJQMQeFYu7Jpn1GM/CgnnFkAR6WUv5cCPFdIcQGNPfPSeCjAFLKg0KIh4E3gDDwCSllRH+ujwHfBtzAE/oPwAPAd4UQ3WgWw/YCvDbFIsWIOVgtgrGp4AKfjUKxOMkoDlLK/cBlKdY/OMMx9wD3pFjfDVyUYt0PvDfTuSgUAOO65bC82lXSlkMkKglForjs1sw7KxTzjKqQVpQcXn8Yh9VCY4WrpAPS3/vdKd76v59d6NNQKFKixEFRcoz7Q1S6bVS57SVtOfSO+Ogb86PCa4piRImDouTw+sNUuOxUl7g4BMNRAEIRJQ6K4kOJg6Lk8PpDVLhsVLrtJd1bKRgxxCG6wGeiUCSjxEFRcmiWg+ZW8gbChEv04hrQLYewshwURYgSB0XJMT4VosJpp9pj1373l2bbbsOtFCxRcVMsbpQ4KEoOrz9sBqShdAvhTMshqsRBUXwocVCUHFrMwV7y4mAGpMPKraQoPpQ4KEqKSFQyGYxQ4bKZbqVRX2lWSZvioCwHRRGixEFRUkzo8YVFYTmobCVFEaPEQVFSjOt9lbRsJYe2VqrioLKVFEWMEgdFSWGIQ2Wc5VCqLTQCYa0fpcpWUhQjShwUJYXRrrvSZcNhs+C2W0vXraQsB0URo8RBUVJ442IOANUeO6MlLg4q5qAoRpQ4KEoKI75Q4dK6zZdy8z0VkFYUM0ocFCWFMeinUo83lLI4BFTjPUURo8RBUVLE3EpxlkPJBqSNmIOyHBTFhxIHRUnhDYRx2S3YrdpHt1QtByml6q2kKGqUOChKCqN1hoEWkC69Cul4V5LKVlIUIxnFQQjhEkLsEkLsE0IcFEJ8Xl+vFUI8JYQ4qv9bE3fM3UKIbiHEYSHETXHrG4UQr+vbviqEEPq6UwjxkL6+UwjRWfiXqlgMjE+FTZcSaJaDPxQ1awZKhXhrQQWkFcVINpZDALhRSnkpsAHYJoS4Evgc8LSUcjXwtP47Qoj1wHbgQmAb8DUhhDFB/T7gLmC1/rNNX78TGJFSrgLuBb5YgNemWISM+0NUxlkOpdpCw3ApAYSiynJQFB8ZxUFqTOi/2vUfCdwCPKivPwjcqj++BfiRlDIgpTwBdANbhBDNQKWUcofUhuZ+Z9oxxnM9Amw1rAqFIh5j0I9BlUdroVFqQel4SycUVpaDovjIKuYghLAKIV4D+oGnpJQ7gWVSyj4A/d9GffcW4HTc4b36Wov+ePp6wjFSyjAwBtSlOI+7hBC7hRC7BwYGsnuFikWFdxFaDmqeg6IYyUocpJQRKeUGoBXNCrhoht1T3fHLGdZnOmb6edwvpdwkpdzU0NCQ6bQVi5DxaZZD9SIQB1XnoChGZpWtJKUcBZ5FixWc111F6P/267v1Am1xh7UCZ/X11hTrCccIIWxAFTA8m3NTLA28/pBZAAeUbPO9QFgFpBXFTTbZSg1CiGr9sRt4G/Am8Dhwh77bHcBj+uPHge16BlIXWuB5l+568gohrtTjCbdPO8Z4rtuAZ/S4hEJhEopE8YeiVDgTs5WgBC0Hla2kKHJsmXehGXhQzziyAA9LKX8uhNgBPCyEuBPoAd4LIKU8KIR4GHgDCAOfkFIa0bePAd8G3MAT+g/AA8B3hRDdaBbD9kK8OMXiYnp1NMTaaBjN904P++jun+CtaxuTn6CICITiYg7KraQoQjKKg5RyP3BZivUhYGuaY+4B7kmxvhtIildIKf3o4qJQpMNrDvqJuZWsFkGFy8b4VIhAOMKHv/0KZ0aneONftqV7mqIg3nJQFdKKYiQby0GhKArGp5ItB4i10Pj/nu7maL+Wde0PRXDZrUnPUSwkZCspy0FRhKj2GYqSYXpHVoNqj529PSPc99wxasu0uocRX3G31AiqgLSiyFHioCgZxlPEHECzHE4N+agtc/B329YAMDxZ5OIQiSuCU5aDoghRbiVFyeCNmx8dT7VbsxbuufWiWIC6yFNbjYC0w2ZRloOiKFHioCgZ0lkO79vcxsWtVbzjwiYOn/MCJeBW0gWh3GlTFdKKokSJg6JkMCyHcmfix/b6Cxq4/gKtYr7Go1kOI0VuORgxB4/DSjCs3EqK4kPFHBQlg9cfpsxhxWZN/7Gt1hvxjRZ5zMGokFaWg6JYUeKgKBmmD/pJhcNmocxhLSnLQcUcFMWIEgdFyTC9XXc6qj0ORos85hAIR3FYtXGnKltJUYwocVCUDOP+UFbiUFNmnzEgPTYV4mvPdhNZwCE7wXAUh82ispUURYsSB0XJcGZkivpyZ8b9ajyOGd1KP331DF968jBvnhtP2jY0EcjrHLMlGIngsFmwWYSqkFYUJUocFCVBz5CPk0M+rl6ZNAMqiUxupdfPjAHJ0+NOD/vYfM9v+N3xofxONguC4ShOm+FWUpaDovhQ4qAoCZ47oo0LuX5N5m6rNR77jJbDAV0cRqeSxSEq4eTgZB5nmh2GW0mJQ+nz9eeO8V+/7V7o0yg4qs5BURI8d2SA9loPnXWejPtWexyM+0NEohKrJXHIoD8UMZvzTa+iNsRiaB7SYGMBaaEC0iXOj3b1cHpkit+7uJnO+rKFPp2CoSwHRdETCEd4+dgQ11/QgDYnamZqPHakTD0A6FDfuBmIHp1KFAEjiD0ffZkMy8FmtRBWlkPJEopEOT0yRSQq+crTRxf6dAqKEgdF0bPn5Ai+YIQb1mQ3N7zGk74zq+FSgmTxMCyJkfkQh0jMrRRUlkPJcnrYRyQqaa/18NPXznD0vHehT6lgKHFQFD3PHhnAYbVw5YrMwWiAGqNtd4qL/Otnxqjx2GmscCYFpMfm2a2kBaSFqpAuYU4OafGp//mu9XjsVv7jN4vHelDioCh6njs8wOauGsqc2YXIZuqv9PqZcS5qqaLaY0+OOWTpVgqGo+Q74lxzK1m1gHRYiUOpcnxAE4fLO2r402u6+MXrfbxxNjlFuhRR4qAoavrGpjh83ms21suGdG4lfyjC0fNeLm6potrtSBFz0MRiJnGQUrL1y8/yzZdOZn0+qTAC0jarILSAxXiK/Dg5NEmly0aNx85HrltBhcvGV54+stCnVRAyioMQok0I8VshxCEhxEEhxKf09X8WQpwRQrym/9wcd8zdQohuIcRhIcRNcesbhRCv69u+KvToohDCKYR4SF/fKYToLPxLVZQizx8ZAOD6CzKnsBpUe4yZDokX+cPnvISjkotbqqjy2BnTx44ajGUhDoMTQU4PT3FqKL9012A4gtNmwaFSWUuak4M+uhrKEUJQ5bHzgS3t/OZQ/7wVU84l2VgOYeAzUsp1wJXAJ4QQ6/Vt90opN+g/vwTQt20HLgS2AV8TQhjDfO8D7gJW6z/GFPg7gREp5SrgXuCL+b80xWLguSMDNFW6uGBZedbHlDtt2Cwiya1kFL9d1FKlzZ2eJh6GJTEVijAVjJCKnmFNFHxptmeLEZC2WSxIyYK28lDkzonBSbri0qvffXkLkajk5/v7FvCsCkNGcZBS9kkp9+qPvcAhoGWGQ24BfiSlDEgpTwDdwBYhRDNQKaXcITWH7XeAW+OOeVB//AiwVWSTs6hY9Lx+ZozNXbVZpbAaCCFSVkkfODNGtcdOa42barc9qQhuxBfCbtX+n+E0FdanhnwAacUjW8wKaZv2/ynrYW7oHfHxNz/elzCzu1D4QxHOjk0l1DasbapkbVMFP33tTMH/v/lmVjEH3d1zGbBTX/qkEGK/EOKbQogafa0FOB13WK++1qI/nr6ecIyUMgyMAUmpKUKIu4QQu4UQuwcGBmZz6ooSZXwqTK1n5jbdqajx2BmZTLYcLm6p0sXDji8YMS8aUkrGfCHaa7W7wOGJmcXBFwyn3J4tZoW0RfsKKnGYG3YcG+LHe3pNi6+Q9Az7kBK6phW+3bKhhVd7RvN2PS40WYuDEKIc+AnwV1LKcTQX0UpgA9AH/Luxa4rD5QzrMx2TuCDl/VLKTVLKTQ0N2QcoFaWJlJKJQJjyLDqxTkdrvhe7wAfCEY6c93Lh8ioAqvRZ00b66lQoQjASpatec1+lsxx6hg1xyM9yiK+QBlSV9BxhjGOdCuYmvr98vY/T+t98OkamUrI4LAfgsdfOZnz+7v4JPveT/UV5c5CVOAgh7GjC8H0p5aMAUsrzUsqIlDIKfAPYou/eC7TFHd4KnNXXW1OsJxwjhLABVcBwLi9IsXiYCkWIRGXGAT+pmJ6qeuTcBKGIFowGqNIzmsb0OIOx78oG7Ys+PJk6oGjcDU6F8ncrGRXSgKqSniMMyzCXv1c0KvmLH77KBx/YmbLa3qhxmN4yY3m1myu6avnpa2cypjx/86UT/OiV0xw9PzHr85trsslWEsADwCEp5Zfj1pvjdns3cEB//DiwXc9A6kILPO+SUvYBXiHElfpz3g48FnfMHfrj24BnZL6J5IqSx+vXXDfZzHCYznTLwQhGG+JQ7TYymvSqaH3fFaY4pG7cF3Mr5S4O0agkHJXaPAddHIJKHOYEQxxycQMaNycnh3x8+qHXiE5LGjg5OEldmYPKFDcvt17WwvGBSfNzl4pIVPLrg+cB5sTtlS/ZWA7XAB8EbpyWtvolPS11P/BW4K8BpJQHgYeBN4AngU9IKY1v0seA/4MWpD4GPKGvPwDUCSG6gU8DnyvIq1OUNF6/doHOyXIo0ywH4x5jf+8o1R47bbVubbsn0a1kpLG21XqwWkRKy2EiEDarp/MJSBtC4LRZseluJTXTYW4wxMGfg+UwqQvKpa1VPP1mP199JrH6+fjgZJJLyeDmi5pxWC389NX0rqU9p0YY1FNejZuOYiLjLZmU8kVSxwR+OcMx9wD3pFjfDVyUYt0PvDfTuSiWFuN5Wg7BSBRfMEKZ08a+3lgwGmIxB8NyMDKXassc1HgcKWsdDJdSXZkjr4B0QL9gGb2VQAWk5woz5pCDOPgC2jF3XN3Ji92D/MdvjnJpWzVv1dvGnxyc5C1pijOrPHZuXNvIQ6/08I4Ll6Vs/fLEgT4cNgtOq4VTaeIaC4mqkFYULaZbKcu2GfHEWmgEmQpqwehLW6vN7dVuLeZgiILhVqp2O6grSy0OPfrd3brmyrzcSoGwdqzDpgLSc40Zc8ghIG1YDmVOG//27ou5YFk5n3/8IKFIlMlAmH5vIK3lAPDPf3AhzdVubv/mLp5643zCNiklvzpwjresrmdFY7n52SomlDgoipYJ03LIJSCtX/x9Id7Q23Rf0lplbq9w2RACsxDOsCCqPXZqyuypLQf97m5NUwWBcDTnwjXjguW0KsthrgnkEZA2bgDKHDZcdit/t20tJ4d8PLKnlxODqTOV4mmqcvHwR69iXVMFf/69PfxkTyyTf3/vGGfH/Gy7qJmOWg+nSjTmoFAsCLGYQ25uJdBaYezvHQXgkjjLwWIRVLrssZjDVAiX3YLLbqWuzJnGreSjxmNnWaU2xzrXjKVgnFvJzFZSnVnnBMOtlFPMIaDdnHicWoOHG9c2cnl7NV99+iiHz2mtuTvrZh7uU1vm4PsfuZIrV9Ty2Uf28fBurQTsiQPnsFkEb1vXSEedh7Oj/qK7QVDioCha8stWirmV9veO0VjhpKnKlbBPtSdWJT3qC5qupnSWQ8/wJB11Zbgd2vnkGneIBaRjbqVgWLmV5oJ8spXiLQfQKu8/e9Ma+sb83PsbrbleZ33myYTlThsP3LGZa1fV83c/2c+je3t58kAfV62so9rjoL3WQyQqOTMyNetznEuUOCiKFq8/hBCxL+dsMGY6jPpC7OsdTbAaDKrd9rhU1pCZwVRb5mR0KpTkNjo15KOjzoPHrt1JGgHL2RIIJQekleUwN+QVczAsB4fVXLt6ZT3Xrqqnd2SKZZVOPFl+Nl12K/d/cBNXrajjMz/ex8khH9suagKgQ7c+ii0orcRBUbSM+8OUO2xYLLNvs2XUMZwe9nF8YJJL4+INBlUeh2k5jMWJQ12ZAykTu7oGw1HOjk7RUesxLxa5BqUNy0FlK809+RTBmZbDtISIz960Bpg53pAKt8PKA3ds5oquWpw2C+9Yb4iDZn30FFm7jdnfkikUs0BKybOHB7j+goZZX+QnAuGcXEoANquFCpeNF44OAnBxKnFw283WCKNTQVborTMMq2N4MkhduRZfODM6RVRCe10Zbl0cpkI5upWMmIPVgs2ispXmkrxiDsFkywFgQ1s1n9q62iyYnA1uh5XvfPgKBiYCNFRon62GcidOm6Xoah2U5aAoCP1eP199+mhSG4hnjwzwp99+heeOzL5RotcfyilTyaDG4+CwPtM3vVtJsw5GfCFqymKWAyTOdTBqHDrqPKYrIWfLIS4g7bApy2EuibmVcqtzsFoETlvyZfKv334Bt2yYqTl1ehw2Cy3VbvN3i0XQXutRbiXF4uTpQ/18+akj7DqR2BLr5W7tzr27f/a9Y7z+3C0HiAWl22rd1OoX/HiqPVq2UjSqdWSt0gPStSnEwWi4Vwi3kpFe6bRZTctBVUjPDfm4lSaDYTwO66zaxedKR52n6GodlDgoCoIRvHtBFwODHceHADiRgz81X3Ewah0uaUm2GkBzK0Ul9HsDBCPRuIC0dtxQguXgw2230lDhjLmVchaH+CI41VtpLglE8rMcckmGyIX22jK9BXjx3CQocVAUhEk9c+fFozFxGPOFOKgPWz8xkIs4hCjPy62kHXtJingDxFpoGC4jI4htzqCeJg7ttR6EEPkHpMPxqaxGV9biuSgsJvK2HJzWzDsWgI46D1OhCAPe4hkvqsRBURCMPPIDZ8dMd8zOE0NIqbl1TuZgOeQTkIY4yyFFvCF+u3Fuxu8Om4UKp22a5TBJu55V4rEXps4hsX2GshzmgqBupeUkDoHw/FkO+mermOIOShwUBcHI7JASXtJdSzuOD+G0WXj3hhb6xvyzNu3H83QrtVS7cdgsXNRSmXK74UY6qft6q+MmztWWx1p+R6OSnmEfHfqUuHzdSgnZSiqVdU4J5uFWmgxGkjKV5grjs1VMGUtKHBQFwReIsLzKRaXLZrqWdhwbYmNHDRc0VQDMynoIhLURnql65WfLB6/q4MlPXZc24ynJrRQvDnHN9/b2jBAIR7lQFxmHTUtB9eXZPsNpj81zUKmsc0M+Lbt9wXBSjcNc0VrjwSKKq9ZBiYOiIEwEwlS67Vy9sp4Xjg4wPBnkzXNerlpRZxYLGc3KsiGf1hkGLruVFQ3labcbMYaTg9rdmhFrAKj1OBjS50j/9LUzuOwW3q4XLYFmPeQekI63HIxsJWU5zAV5FcEFIvMmDg6bheYqt3IrKRYfPt0Ev+6Ces6O+fnhrh4ArlpZZzYnm404GB1Zy+fwy1k5zXIwLAnQLIcRX5BQJMov9vfxtnXLEs7F47DmHnMIR7EIrVAvVgSnxGEuiBeH2WYCTQbDlM2TWwm0oLRyKykWHZO6CX7dKm34ydefO4bbbuWS1mrKnDaWVTpztBxydytlwmW34rJbmAxGzI6sBrXlDoYmgzx/ZIARX4hbpxU8eRy2vNpnGMVvQgjsVkEox/bfipkJRjQhljJmsWWLLxDJundSIeio85j1NMWAEgdFQdC+SFba6zy013rw+sNs7qo1L4KddWWcnJU45N6uezaYnVg9iUVytR4HwXCUH+zsodpjT5r45bbn7lYKhqNmrAHAbrUQmuWFS5GZaFQSikjTQpzN30tKqd/wzJ/l0F5bxvBk0PzsLzRKHBQFQTPBtQv5davrAbgqbjTiioayWVkO+YwInQ1GEDrepQSxQrhnDvfzexc3myJnoLmVco85OOOsFJtFEFaWQ8ExMpWMv+1s4g6BcJSoZF4tB6NX0+u9Y/P2f85ERnEQQrQJIX4rhDgkhDgohPiUvl4rhHhKCHFU/7cm7pi7hRDdQojDQoib4tY3CiFe17d9Veh16UIIpxDiIX19pxCis/AvVTGX+IIRs2DobeuXIQRcH3e33VlXxtBk0ByukwnTcnDOnVsJYheO+EwliImDlHDrZck9dNwOa87ZSoFwJMFycNgsqkJ6DjDe0+ocxMGo+J9Py+Etqxuoctv53s5T8/Z/zkQ2lkMY+IyUch1wJfAJIcR64HPA01LK1cDT+u/o27YDFwLbgK8JIYx3+D7gLmC1/rNNX78TGJFSrgLuBb5YgNemmEcmArG0v7euaWTH57ayfnmsvqBTz1jK1rU0EZgfy6FqWlW0gSEOLdVuNrbXJB1X5rAxlUdAOr6Zm81iUdlKc4ARjM7FrWRYhfNpObgdVrZvbuNXB8/TN7bwg38yioOUsk9KuVd/7AUOAS3ALcCD+m4PArfqj28BfiSlDEgpTwDdwBYhRDNQKaXcIbW0ge9MO8Z4rkeArYZVoSh+QpEowXA0oZp0+tS1FYY4ZJnHbQSky+fJrTTdcqjXW3X/wYblKVuN5+NWCoajCW4qu02oOoc5wBCHXNxKRlHnfGYrAfzJlR1EpeT7v+sx10KRKD/ff9bsyTVfzCrmoLt7LgN2AsuklH2gCQjQqO/WApyOO6xXX2vRH09fTzhGShkGxoA6piGEuEsIsVsIsXtgYPYtoBVzQ+wuK/0Xqa3WgxBwPMseS15/CLfdavYemiuMlhlGR1aDtloP//7eS/nz61emPC6fOodgJNFysFssKpV1DkgSh1n8vYxeYZ55qnMwaKv1sHVtIz/c1UMgrKXf/v2jr/PJH7zKw6+czvwEBSTrb54Qohz4CfBXUsrxmXZNsSZnWJ/pmMQFKe+XUm6SUm5qaGhIcYhiITDy/WcqGHLZrbRUZ99jKd+OrNmSLuYA8J6NrUmBaoOCWg5WJQ5zQT4Bad8CWQ4Ad1zdydBkkF/s7+MrTx/lx3t6sVoEL3UPzet5ZPXtE0LY0YTh+1LKR/Xl80KIZilln+4y6tfXe4G2uMNbgbP6emuK9fhjeoUQNqAKSBwMoChazLusDF+krvrsM5a8/vCcu5QgPuYwu8C322FjKhQhGpWznnAXCEdx2eNiDlahurLOAdMth9m00Ih9pud/WOY1K+tZ0VDGPb84xNBkkPdc3opFwK/fOE8kKrHmMDY3F7LJVhLAA8AhKeWX4zY9DtyhP74DeCxufbuegdSFFnjepbuevEKIK/XnvH3aMcZz3QY8I4upsbliRmJ3WTN/kQxxyOZP6w2E57QAziCWypo8DGgmDCH05+AHTlXnoLKVCo9R9Gb8jWcXkJ7/bCUDi0Vwx1Wa9XDd6nq+8J6LuXZ1PWNTId44O5PTprBkI4vXAB8EXhdCvKav/T3wBeBhIcSdQA/wXgAp5UEhxMPAG2iZTp+QUhp/lY8B3wbcwBP6D2ji810hRDeaxbA9z9elmEcmApndSqCls3r9YYYmg2bANx1ef4jKebAcmqu0cY3Lq10Z9kwkfqbDbO8uk91KynKYC/ILSC+c5QCwfUsbbruVmy9pxm61mDVDLx0bTDkPfS7I+MqllC+SOiYAsDXNMfcA96RY3w1clGLdjy4uC0U0Knn19AjLq900VbrmZTTgYsGnm+CZ7rK6GmI9ljKLQ5jmqtldsHPh8vZqfvmX1yWk3WaDWy9i8wUikL63X0q0gHTsvbJbLeaFTFE4DGvMSGWdTYzItwB1DvE4bVbetznmnW+sdLG6sZyXjw2lTZIoNAsji0XIbw/3c+eDuwHNDL2ktZovv+/SjBcxRSztL9Nd1vrmSuxWwed/dpBvfWgLDRXp31uvPzTnBXCg9TaarTBA7LX6QrOvdZhuOdisFvNOVVE4ArqlUOG0YxGzjDkEIwgBLtvCiEMqrllVz0OvnE76/MwVqn12h/skAAAgAElEQVSGzpvnvAD8j3et59pV9Tx/ZIC9p0YW+KxKA+OOLNNd1rJKF/ffvonu/glu+/rLMw5Un6+AdK7kMyo0EI4kfLkdVqF6K80B8RP3ZtsLyxcI47FbZ51sMJdcvbKOqVCEV3vm57qkxEHn5OAkDRVO7ry2i7tvXgckDphXpMdoNZCNf/ataxr5wUeuZGwqxB/e9zLHBiaS9olEJb5gZF5SWXMln2lwgWkBaZvFQjiqxKHQmBP3bBatLmWWlsN81zhk4ooVdVgEvHRsflJalTjonByaNIfS1OmtE4aVOGSFaTlkmRN+eXsNj/z5VQRCEf7914eTtk/MQ7vufMnHcpjePsNus6gK6TkgXhxc9tmJg2+eZzlkQ5XbzsUtVbysj+Gda5Q46JwYnKRLH0rjslspc1jNSWCKmZkMhHHaYvOQs2FVYwV/clUHTxw4l1T7MD5P7brzISYOs4s5SClTVEgLVQQ3B5huJavmVpptncNCZSrNxNWr6nnt9Khprc8lShzQgp+DE0GzORwYw14CC3hWpcNkjrN2//SaTuwWC9944XjCutFXaT5SWXPFrV84ZutWCkclUjJjhXQoEmXvPPmVFzPxs7pnW9Hum+dZDtlyzcp6wlHJrpNzXyOsxIHYDOGuOHGoK3Mqt1KWGIN+ZktjhYv3bGzhkT299Hv95rrRrrt8HrKVcsVjz82tFIhzdRhMr5B+8sA5/vBrLxdFZ85SJn5Wt2uWAenJHOpX5oNNnTWsbChTlsN8cULv95MoDg7lVsqS+EE/s+Uj160gFIny4MsnzbX5atedD2ZAepYzHYJxFyyD6RXSA17NYh31FcdEsFIl/r12O2bnVvIFitNycNmtPP2ZG3jXJcvn/P9S4gCc0DuFdtR5zLXaMuVWypb4QT+zZUVDOTetb+K7O06ZouCdpylw+eC0WbCI2cccYkHS+CK4RMvBeP25NvZTaAQjUexWgcUitFTWWQWki9NymE+UOKBlKi2vciUMmK8r19xKqsVTZiYDuVsOAB+9fgXj/jA/3q21JI7Njy5et5IQAo/DNusLuOkHnyHmMBHQXn+uLcEVGvE9rGYrDpNFmK003yhxQM9UaihLWKsrcxCKSHOWsSI9k4FIXib4Ze01XNxSxaN7zwDzNz86X3KZ6RCMaPtPr5DWAtXajUjMclCfvXyIryR2OaxMBbPPCJsMhIuuzmG+UeKAJg6dddPEoVzVOmRLPjEHg1s2LOf1M2McG5hgIhDGbhUJd9fFSC4zHfyh5IC0w6pV4Rq1DoY4zDaeoUgkXhw8dmvWY12D4SihiFSWw0KfwEIzog+9jw9GQ2yG8LCKOyQQDEeTBDOfmIPB71+6HCHgsVfPaH2VXPaib36Yk1spkipbSXtsVEl7AyrmUAiCkZg4GBXS2biJfVn2ClvsLHlxSJWpBFoqK8CgylhK4JsvneAd9z6X8CXLN+YAWt+lq1fW8di+s/M2BS5fPA4rU7NsvGfGHKZlKwGEwobloMUclDjkR3zMwWW3EpVkNTdjMsteYYsdJQ56plLndHFQbqWUHB+YYHAiaKZZhiNRAuFoQe6ybrm0hVNDPnYcGyoZccg5IG2PFwfdrWRYDoZbScUc8iIQjppZYUaLdX8WcQffLHqFLWaWvDicHJrEIqCtxpOwXqv6K6XEsKT69Vx8X6hwd1nbLm7CYbPQ7w3MS7vufJltp0+Iz71PnOcAmBlLynIoDNPdSpBdHEdZDhpLXhxODE7SWuNJ6o/uslspd9oYnIjFHLz+EJf9y6956o3z832aRcOQ/n4YFc2+As7arXTZ2bq2EaCo23Ub5GI5pKyQ1ttCG7UOE6rOoSAEwxHTfWdYDtmIg7IcNJQ4DE4mxRsMasscCZbD0f4JRnwh9p0ena/TKzpMy2FcE4mJAk/MumWDVvlZCm4ld04B6eRUVuNxKBIlEpXmnauqc8iP+Gwl9ywaJZqWgxKHpYuUkpMZxCG+hcZJvXvo2SXa80ZKyYBpOehuJWMQe4G+SDesaaS+3EFLtbsgzzeXeBzZp0caBFNaDoY4SNNqgJjLTpEbCW4lI+aQjeUQLOwNT6mSURyEEN8UQvQLIQ7Erf2zEOKMEOI1/efmuG13CyG6hRCHhRA3xa1vFEK8rm/7qtDzFIUQTiHEQ/r6TiFEZ2FfYnoGJgJMBiN01nlSbq8vdyQM/DHFYXRpioM3EDYvbufHNbfSpOFWKtAXyWW38uu/vp5P3riqIM83l3gcVnxZpkcapK6QNuocoma7clAB6XxJqJA2hzNlka1kzkRXlkMmvg1sS7F+r5Ryg/7zSwAhxHpgO3ChfszXhBDGVeM+4C5gtf5jPOedwIiUchVwL/DFHF/LrDEylboaUk+I19xKsZjDCX2sZd+YP+X+i514K2pgjiwH0N53ZxHN7k2H22FFylgcIRtSxRziA9ITcd02VcwhPxLcSrOJOZh1DsX/GZxLMoqDlPJ5INvm4bcAP5JSBqSUJ4BuYIsQohmolFLukNpt1neAW+OOeVB//AiwVcxT9dMxXRxWpHUrJfZXMiyHvlE/0ejS67lkBOctIhaQXsqZHbm07Q6k6Mpq0y2HcFSaaaw2i1DikCcJ7TNmIQ6TBUyyKGXyiTl8UgixX3c71ehrLcDpuH169bUW/fH09YRjpJRhYAyoS/UfCiHuEkLsFkLsHhgYyOPUNfb2jFDjsdNak9q/XV8e669kxCecNq298lKcLz2oWwsrGspjMYclnNlhvObZ9NZP17IbIBSOmmmsDRVOFZDOk1SprP4s3lNfMIzLbsFqKe4K/bkmV3G4D1gJbAD6gH/X11O9m3KG9ZmOSV6U8n4p5SYp5aaGhobZnXEK9p4aYWNHTdo2DfG1DoMTQbyBMJs6NR1cinGHQV0QL1xeSf94ACllXLbS0hOHXGY6xLeRNogVwcUsh8YKJ75ZVl8rEglM68oK2WYr5V/xvxjISRyklOellBEpZRT4BrBF39QLtMXt2gqc1ddbU6wnHCOEsAFVZO/GypnhySDHByfZ2FGbdp+6cq2FxtBEgJN6m42rV9YDS1QcvAGEgDVNFUyFIkwEwqbrYyn6Z2NzpGchDuFoUjwlwXLQxbax0qUshzzR3mu98Z4p5NlUSOffK2wxkJM46DEEg3cDRibT48B2PQOpCy3wvEtK2Qd4hRBX6vGE24HH4o65Q398G/CMnIchCntPaTN6N3bUpN2nTrcchiaDnBg0xEHzeJ1dgkHpwYkANR4Hy6s0N1y/N8BkMIzDZjEvcEuJ2eTOGwxNBJK6zRqprOFozK20rNKpYg55IKVMcCsZ73l2FdLKcgDI+A4IIX4I3ADUCyF6gX8CbhBCbEBz/5wEPgogpTwohHgYeAMIA5+QUhp/jY+hZT65gSf0H4AHgO8KIbrRLIbthXhhmdjTM4LNIriktSrtPoZbaWgiSO+ID5tFcFFLFS67ZUlaDkMTQerKHDRWaBbV+XE/vkBkybY2NmIO2d7hHzw7xs/29/G+TW0J6w6b5lYKRjS3kt0qqPE4zC6ixd6dthjR5mPEYjtCaNPgsqtzyG0m+mIjozhIKf8oxfIDM+x/D3BPivXdwEUp1v3AezOdR6HZc3KEC1uqEqa/TSe+bffJoUnaaz3YrRaWV7uX5PD3wYkA9eVOGis1cRjQLYelGIyG2bmVIlHJ3//fA9R47PzdtjUJ20zLIRJlwh+m3Gkz02T9oahpoSiyJ1WxYbbDmSYD4SUZQ5vO0vMFoH1w9vWOsrE9vUsJYv2VNLeSz+zc2lLt5szo0nQr1Vc4aahwAVoLDV+eU+BKGTN3PosLzvd3nmLf6VH+x7vWU+1xJGyzx7XPMGZZeGYRQFUkk1IcshwVqiwHjSUpDm/0jRMIR83Mo5moK3cwOBHkZNy0uOYq15J2K1W6bLjsFvq9fs0/u0TvsjxZxhzOj/v5X08e5tpV9fzBpcuTttstsUlwxiwL02WlWmjkRKqhSi67JTvLQcUcgCUqDnuyCEYb1JY5eLNvnKlQhK56rc3G8mo3A94AgfDS+eL6QxG8gTANFU6EEDRWuLSAdAEG/ZQqxgU8Uw+kLz15mEAkyr/eelHK+EF8hbQhDrF2D0vnM1ZIUtWTeBy2LLuyqmwlWKLisPfUCC3VbpZVujLuW1fm4Gj/BABd9VqbDSNb5/zY0hkhalRH1+tDkBornJpbaQmb4C67BSFmvoD3e/08vu8MH9jSnjRQysCskI5IvIEw5U57Tmmyihip2pSkm78RCEf49EOv8cNdPYQjUWU56Cy5d0BKye5Tw1zRlbIIOwljXChAZ5zlAFp31vY0TfsWG0arbuP9aKx0cvicl2AkumTdSkIIPPaZZzr8YGcPoYjk9qs60u5jWA5BPeZQ6aqIS5NV4pALqRocuhxWxqdCSfse65/k0VfP8OirZ/jWSyfwhwoz2bDUWXKWw5nRKc6PB7JyKQHU6nfKDpvFtBiWV2sWx1KKOxhDfur1NFbDreQLLF3LAbSZDunaZwTDUb6/s4cb1jSwIk1zR4iJQzhlzEEFpHMhVczBbbekTGU1Zrb8+fUrTVGp9hT/JMK5ZsnJ42ziDRArhOuo9ZgtD5p1kVhK3Vmnu5UaKpx4/WGsFrFkLQeAzjoPh897U2574kAfA94AH7q6c8bnsFoEFqENApoIhCl32ZRbKU9SjWN1p7HyhvTOy++5vIXPvOMCnj8ywBUrsvMsLGaWnOXgD0Xoqi9jbVNFVvvX6RfDeH+x22GltszBmSVkORhupXq9pYgRr4lE5ZL2z27srOHAmbGUd6TfeukkXfVlvGV15j5gNqsFrz9MJCqpcNnjegEpcciFtHUOM1gOtWUO7FYLW9cto3wJ3/AYLDlxeP/mdn772RuwZdnuoVb3sU+fFrfU0lkHJwKUO21m0aBRJQ1Ls123weaOWkIRyf7esYT1106P8trpUe64qiOhyV46HFaLeZGqiLMcVLZSbqQax+q221J2ZR2eDCIESfUnS50lJw6zxXArGTUOBsur3fQt4kK4+549xmd/vM/8fXAiaLqUALNKGpZmu24Dwz35ysnEXpHfefkk5U4b79nYmuqwJGxWEScO9liarBKHnEiVyup2WFJaDkOTQWo8jiXfons6ShwysL65krvfuZZ3XdqcsN5S7V7UlsOL3QM8ureXEf2CNegNmC4l0ALSBkvZcqgpc7CyocyMZYHmuvzVwXO865JmKlzZBTbt8ZaD0xaXJqsC0rmQLpU1HJWEIomdWYcngmarHEUMJQ4ZsFgEH71+JZXTvuTNVS68gXDCzN/FxNBEkKiE545oQ5WGJgNm/AWgxmM35xAsZcsBYFNHLXtOjZjTAZ87MsBkMMLNFzdnODKG3SIY8cXcSkajOGU5JHJ8YIK9PSMZ90uZyppmGtzwpBKHVChxyBGj1mGxupaMu9in3+wHDLdSzHIQQtCg/75Uu7IabOysYWwqxLEBrVjyl6/3Ue2xc9XK7DNe7DYLI5PajYZhbXgc1ozV10uN//Wrw3z4268QyTCmN2Uqa5o4ztBkwHQfK2IoccgRs9YhQ3fW08M+zpVYyquU0hSH5w734w9FGPEligNAg56x5FnimR2b9LjD7lMj+EMRnj7Uz03rm2Y148JmEeYFrcKlvZ/ZdhFdSpwd8zPqC/FqBushleWQrlHiiC+kLIcUKHHIEbNKOkPc4ZM/fJV//Onr83FKBWN8Kkw4KtnUUcO4P8xvDp1HShIC0gDL9Iyl8iUccwAtk62uzMHukyM8f2SAiUCYmy/J3qUEJAhJuS4OHrtNdWWdxsC4dqP1jG7RpiNVzMGTYqxrJCoZ8QWV5ZACJQ450ljhotJl41cHz6fdR0rJ8f4JekdKK3BtFAXdclkLdqvg4d29AEmWg5GxtNRjDkIINnbUsPvUME8cOEe1x25ODMwWuzmUBsodMctBxRxiRKOSfq/22cwkDuka7wHmnG6AUV8QKVGWQwqUOOSI1SL4y62ref7IAL89nPqDOuoL4Q2EzeriheT0sC+rKVgQizd01HrY0lXLC0e1oHR9RaI4LNMzlpZyEZzBps4aTg35ePLAOd6xftmsx6Yawf1yh82si/Aot1ICw74g4aikpdrNm+e8Mw7cCoajWAQJ9UzNVdrnNf44swBu2o2PQolDXtx+VSdd9WX868/fSEqPA+gZ9gFaHnU4xfb5wh+K8M6vvMA9vziU1f5GNXRtmYO3rmnEmOg93fR+76Y27nn3RVSpPjRs7KgFNJfFbLKUDIyLmOFSAj0grcTBpH9cu8kyxqz+9s2BtPvGz482aDZdwbEY4JAhDqoALgklDnngsFn4h5vXcWxgku/97lTSdkMcpIzdoSwEu04MMxEI8+jeXibSNImLxzjXunIHW9ctM9enWw5NVS7++Ir03UaXEhe1VOKwWah02bh6Zf2sjzfcHxVx4uDOcv7AUuG8V7uoX7u6jpZq94yupWA4muBSAih32qhy2xPihPGtMxSJZBQHIcQ3hRD9QogDcWu1QoinhBBH9X9r4rbdLYToFkIcFkLcFLe+UQjxur7tq0KfeiKEcAohHtLXdwohOgv7EueWresauXZVPf/xm6NmwZiBIQ4AAwvoWnqpexCAyWCEx187m3H/YT3mUFvmoKu+jBX1ZThsFiqWeFbSTDhtVt67sZUPX9uVdMeaDcZMh/iiOa0duApIGwzolkNjhYsb1zbyUveg6Sod9QU5EtcAMRCO4rAlJ0osn1a8OhR3I6RIJJtP8beBbdPWPgc8LaVcDTyt/44QYj2wHbhQP+ZrQgjjL3QfcBewWv8xnvNOYERKuQq4F/hiri9mIRBC8I/vWofXH+JbL59M2HY6Xhy8CycOLxwd5IquWtYsq+AHu5ItnOkMTQapcNpw6l+u7VvauGZlXcopZooY97z7Yv7qbRfkdKw9peWg3ErxnNczlRorndy4tpGpUIRdJ4bZ2zPCtv94gVv+8yXTfRsMRxPSWA1aql0JDTOHdRdqjXIrJZFRHKSUzwPD05ZvAR7UHz8I3Bq3/iMpZUBKeQLoBrYIIZqBSinlDimlBL4z7RjjuR4BtooSuwqtbapkXXNlUu51z7DPTP9cKHEYmgjwRt84162u5wNXtHPgzDj7e0czHBM051gA3PWWlXzrT7fM9akuacyAtDMx5qAC0jHOe/1Ue+w4bVauWlmH02bhi0++yfv/eweDEwGmQhHTQk8Vc4Bky2F4MkCFy5aTtbfYyfUdWSal7APQ/23U11uA03H79eprLfrj6esJx0gpw8AYkDIPUAhxlxBitxBi98BA+mDUQrCuuZJDfYl9/XuGfVzWrnncFsqt9NKxIQCuWVXPrZe14LJb+MHOnhmPUe0E5p+Y5RDnVnJovYCMtMylTv94wMyQc9mtXL2yjoNnx7lmVT1fuu0SAM7rrqdgOJIUcwCtJ9q4P4xXb3szNKlqHNJRaLlMdccvZ1if6ZjkRSnvl1JuklJuamjI3CN/PlnbVMHgRMC0EEKRKGdHp1jbVEG507ZglsNLRwepdNm4pLWaKred379kOY/vO2t+OVKhfWFUat98YrNoX8XKaQFpUG27Dc57AwndgP/+5nV86T2X8M07NnPBMm0+i9GNIBhObzlAbFCXuhFKT67icF53FaH/a6QN9AJtcfu1Amf19dYU6wnHCCFsQBXJbqyiZ31zJQBvnhsHtMrpqIS2Wg/15Q4zPXQ+kVLyYvcgV6+sN9sRf+CKdnzBCJ/98T6+/twxHt3ba/pyDYYmVK+Z+cZhS3YrudM0iluq9I/7E7oBr15Wwfs2t2GxCHP4VL+e0TSTWwkw4w6aOKgboVTkKg6PA3foj+8AHotb365nIHWhBZ536a4nrxDiSj2ecPu0Y4znug14Ro9LlBRrdXE41KeJw+lh7cPXXuuhocLJgHf++yudHPJxZnSKa1bHUis3tFVz04XLeOHoIF944k0+/fA+/udjZiIaUmrtBGpV9sa8YlgOFdPqHACVsYRWHT3gDbCsMvWFvK7Mgc0iEi2HNG4liLW9GVZupbRkzE0UQvwQuAGoF0L0Av8EfAF4WAhxJ9ADvBdASnlQCPEw8AYQBj4hpTRuez6GlvnkBp7QfwAeAL4rhOhGsxi2F+SVzTO1ZQ6WVTp5U487GGmshjgcPpd6zvBc8qJe2Xzdqpg4CCH47w9uAmAiEOYvfrCX7v4Jc/u4P0woItUXZp5JFXNwT5sjHYlKfr7/LDdd2GS2n14qGNXRhoUwHYtF0FjhjIs5RPF4ki9vDRVObBbB2dEpdSOUgYziIKX8ozSbtqbZ/x7gnhTru4GLUqz70cWl1FnbVMkbuuXQM+zDYbWwrNJFQ7mTF72D834+L3YP0lLtpqPOk3J7udPGBU0VvNQ9RCQqsVpEQgGcYv6wm3UOyZaD4VZ65eQwn/rRa9x+VQf/ckvSV2lRY6axVqR3ATVWusz9AmliDlaLoKnKxdlRv7oRyoDK3yog65orOTYwQTAc5fSwj9YaN1aLoKHCybg/TCA8f77jQDjCy8eGuG51/Yz1CZ11ZQT14Dlo8QZA+WHnGXua9hkQsxxODU0C8J0dp3jx6PzfbCwkRsO9xjSWA0BTnDikizmAFnc4Mzpl3gipGofUKHEoIOuaKwhFJMcGJugZ9tFWq92xG91M5zMo/d0dp/D6w/z+pctn3M+YjX1Sv/CYFaPqbmpeMSqk4ycOuu1GtpIWczg15MNmEaxoKONvH9m3aKcQpqJfv+inizkY286Nx2IOzjTND40Rv2YnAGUlp0SJQwFZFxeU7hn20a6LQ4NuCs9XOuvYVIj//G03162u55pVM/f56arXxWFQEwfVa2ZhSFUhPd1y6Bn20VLj5svv28B5b4B/+dkb83+iC4QRS2iYwa20rMqF1x/GFwynTWUFbVDXuTE/A151IzQTShwKyIr6MhxWCzuPDzM2FVowcbjv2WOMTYX43DvXZtx3WaUTt93KiUEtgK7EYWEwWj1ML4KDmDic1m84NrRV8/EbVvLInl52Hh+a/5NdAPq9fmr06uh0GAVy58cDGd1K4ag0087VZz01ShwKiM1qYfWycn5zSBsA1LYA4nB2dIpvvXSCd29o4cLlVRn3F0LQUecx3UqDEwHKnbYllw2z0LzrkuX8660XJVyops88PhVnjX7shpUIATtPlFxJUE6cHw+kzVQyaKoyxMGfNpUVYrUOB86MAaiCzzQocSgw65orTb+98UU2PnzzMfTny08dQQKffkf2DeA668pMcVAVowtDU5WLP7kysf25MbnMF4wwNhVi1BezRj0OGx21HvPud7HTP+6f0aUEsXiEKQ5pLAej1uH1M2O47VZThBWJKHEoMGubKszHbbXah9Bhs1Dtsc+55fDUG+f5yd5ePnR1J601qdNXU9FZX8bpYR/hSFSJQxFhtQgcNgu+UNjs8GuIA2ip028uQP3MQtDvzWw5GNv7xvyEozKtODRXxdxP6rOeHiUOBcZoo1HjsSf4jxvKnXMqDr87PsQnfrCXS1qr+dTW1bM6tqveQygiOTvqZ2giaHaSVSw8RmdWs6gyrmZlTVMFJwcnsx7/WqoYs6NnylQCrW7H47CaQppOHCpcdrOHlarnSY8ShwJjtNGIv8MDLe4wV51ZD5wZ488e3E17rYdvf2gzZbMcymOks54YmmRoUt1NFRPawJ+YOLTFfa7WNVcQlXD0/ES6wxP47+eO8Yv9fXNynnPJ0GSQyAzV0QZCCJoqXeZ7lS7mALG4g/qsp0eJQ4GpLXPQWuNmZWN5wnpDhXNOYg7nxvx86Fu7qHLb+e6dW6jJ4cMen86qGpEVF+44y6HGY0+og1jTpKdOZxl3+O/nj/Pd350s+DnOdb2F0Uxvpupog8ZKJ70jWkFnqmE/Bi1KHDKixGEO+P6fXcE/3LwuYa1+jtxKn//ZQbz+MA9+eDPNVe6cnqOhwonHYeXAmTHVTqDI8Dhs+IJheoZ8tOsWnkF7rQe33ZpV367JQJjhySDd/ZMFPb/fHR/i8n95KmHqYaHpH89cHW3QVOmid2RmtxLELAf1WU+PEoc5oKOujLryxLuchgonvmCEyUDhOmw+feg8Txw4x19uXc2qxorMB6RBS2ctY88pbZKd8sMWD8ao0PiiSgOrRXDBsvKsMpaMu+nBiQBjvsLd6R897yUclRw8O3dZU+fN6ujM4rCs0kUoojV1zkYclJWcHiUO80RDeWFrHXzBMP/zsYOsbiznI9etyPv5uuo9HNerpJWpXTx4HFa8/jBnRqdor022DNc2VZqdgGfCuJsG6B4oXIbTgN4SxkiFnguMvkoN5dm4lWIC4rCmT1FdXq3tV1tmT7vPUkeJwzxhFsIVKO7wH785ypnRKf7tDy8uyPzbzjiXhSoKKh48DisnBieJRCUdtWVJ29c0VTA0Gcx402FYDkBCi/Z8MeJoRvuVueD8uJ/aMkdWn/OmeHGYYX8jzparK3YpoMRhnjDEYbAAlsPxgQkeePEEf7Sljc2dtXk/HySKg2pEVjy47TazZXdbbXLtytpmzZ2YybXUO+LDZbfgtFkKKg6GKJ2YQ3E4N+bPKhgNiY35ZhKHS1qreeiuK7k2Q++xpYwSh3nC6MxaCMvhyYPniEQlf/W27KugM9FZH285KHEoFjxx1bvtKeZyrNUzljIFpU8PT9Fa46GrvmxuLIc5ciuN+UK8dGyQy9qrs9p/WYJbaebL2xUr6rBY0rezX+rMLiFekTO1ZQ4sojAxh+ePDLCuuTKrAF22dNZrF54yh1X1VSoiDHFwWC0JLhOD2jIHjRVODmWIO/SOavNFyp029vWOFuz8jM/z+fEAvmDYbPlRKH685zT+UDSptUg6GrO0HBSZUe/ePGG1COoKkM46GQiz59QI11/QUKAz02god1LmsCqXUpFh9P0xBkelYk1TBYfPZ3IrTdFa42ZVYzm9I1MFqaqWUjI4EWC53o7i5GBh01mjUV5ROFMAABRTSURBVMl3dpxic2dNVk0kAZw2q5lQMVOdgyIz6t2bRxrKnRw57+XZw/08e7g/p9zwHceGCEUkb7mgsL5SIQSd9WUqGF1kGJZDqniDwbrmSo6cnyAciabc7vVrTfvaajysbqxASjg2kL9raTIYwR+KsrlLi3udKrBr6bkjA/QM+7j9qs5ZHWfEJ5TlkB95vXtCiJNCiNeFEK8JIXbra7VCiKeEEEf1f2vi9r9bCNEthDgshLgpbn2j/jzdQoivipnmWpYwHXUe9vaM8qFvvcKHvvUKd3xr16yf47kjA3gcVjZ1FCYQHc/d71zHZ9+xpuDPq8gdt+6mmV7jEM+aZRUEw9G0fn8jU6m1xsMqvXK/EHEHwwre1KF9xU8UWBwe3HGSxgonN13YNKvjjNbdmWIOipkpxLv3VinlBinlJv33zwFPSylXA0/rvyOEWA9sBy4EtgFfE0IYzu37gLuA1frPtgKcV9Hxpdsu4dGPX82jH7+a929q49SQj1Cau710PH90gKtW1M3JXdG1q+u5drXK3igmPHr8pyNFMNrAyFhKF3eIiYObznoPFgHHCiAORjC6o66MhgpnQdNZTw5O8uzhAT5wRfusP+vG0B9lOeTHXLx7twAP6o8fBG6NW/+RlDIgpTwBdANbhBDNQKWUcoeUUgLfiTtmUVHhsnN5ew2Xt9ewsaOGSFTSN+rP+vhTQ5OcGvLxlgLHGxTFSzZupVWN5VgtIm3GklEA11rjxmmz0lFXRncB3EpGWnZ9uZOuurKCxRyiUcl/P38Mm0XwgS3tsz5+WZUSh0KQ77sngV8LIfYIIe7S15ZJKfsA9H8b9fUW4HTcsb36Wov+ePp6EkKIu4QQu4UQuwcGBvI89YXF+LKfHsn+C/X8Ee01K3FYOqxpqqCzThsNmg6nzcrKhrK0tQ69I1N4HLFA7cqG8qw7uc6EkZbdUOGks94zo1up3+tHu/ebmRePDvL7//kiP9x1mts2tmbVT2k6119Qz41rG6l2q+rnfMhXHK6RUl4OvBP4hBDiLTPsmyqOIGdYT16U8n4p5SYp5aaGhtK+QBqDgHpmEZR+7sggbbVuOmdwMSgWFysaynn2b96aMW15bVNlWrfS6WEtjdUI5a1qLOfk0GTaAHa2DHoDWISWTttZX8aAN8BEit5hLx8b5Mp/e5oHXjwx4/N9+qHX+JMHdjLqC/GV7Rv4t3dfnNN5beyo5Zsf2oxNxRzyIq93T0p5Vv+3H/i/wBbgvO4qQv+3X9+9F2iLO7wVOKuvt6ZYX9Q0V7mxWUTWGUvBcJQdxwZ5y+oGFmm8XpEHa5srODM6lbJ9tpbGGruhWNVYTigiOZVnJ9WBCW1qoNUi6KqLtX2PZyIQ5m8f2U9Uwv3PHycQTp1Cu+PYEI++eoY7r+3imc9ezy0bWlSB2gKTszgIIcqEEBXGY+AdwAHgceAOfbc7gMf0x48D24UQTiFEF1rgeZfuevIKIa7Us5Rujztm0WK1CFpq3FlbDq/2jDAZjCiXkiIl62aolO4d0SwHg0JlLA14A2blv1FhPz1j6t9+eYgzo1N8autq+r0BHnst+b5PSsm9vzlCY4WTv7lpDU6bKsIsBvKxHJYBLwoh9gG7gF9IKZ8EvgC8XQhxFHi7/jtSyoPAw8AbwJPAJ6SUxm3Ex4D/gxakPgY8kcd5lQxtNR5OxzVEm4nXTmtVrVsK1EtJsbgweyz1JcYdxqZCjPvDtMVZDisbtAt5vuIwOBEwe4Z1prAcnj8ywA929vCR61bwV29bzbrmSr7x/HGi0USv8Y5jQ+w6MczHb1ipqvOLiJxr3aWUx4FLU6wPAVvTHHMPcE+K9d3ARbmeS6nSVuvhVwfPZbXvm+e8NFe5cpr0plj8NFW6qHLbeWNa3CE+U8mgwmWnqdKVdzrr4ETA7G7qdlhpqnRxQs9YGvOF+NxP9rOyoYxPv/0ChBDc9ZYu/vqhfTx7pJ8b1y4DYlbDskon23PITFLMHSpis4C01boZngymDOJN51DfOGubch/oo1jcCCFY21SRlLEUXwAXz5qmipQDevyhCN4sxn5KKRnwxiwH0PpzGYHuT/5wLwMTAb78vg2mNfCuS5azvMrF1587bh7z8rEhXjk5wsdvWKWshiJDicMCYlS9ZgpKB8NRuvsnWNdcOR+npShR1jVXcvicN8FtE18AF89l7dUc6fcm3Zj8Pz9/g6u/8AyH+mbu1TQRCBMIR6mP68XVVV/GycFJvvjkm7xwdJB/vfUiLo1LwbVbLXz42i52nRjm7kf389kf7+PuR1+nqdLF+ze3pfpvFAuIEocFxPADZxKH7v4JwlHJWiUOihlY11yBLxhJqJ3pHfFR7rRR7UnM+b+svQYpYf/pWIdWKSXPvNmP1x/m9m/umvFzORBXAGfQWVfG0GSQb7xwgjuu6uD9m5PdRNu3tNNW6+Zn+/p4qXsQj8PKv956kbIaihDVsnsBMSyHTBlLhqtgfbNyKynSY8x2ONTnpUMPEGtzHNxJ6c9GUd3enhGu1gfe9I5M0Tfm54NXdvD4vrP8yQM7eeTPr05wHRkM6uNBE91K2v951Yo6/vFd61OeY7nTxgt/e2M+L1MxTyjLYQGp9tgpd9oSRjim4lDfOA6bJWFam0IxnQuWVSBE4lS46WmsBlVuO6say3m1J2Y57DwxDMAfX9nOt/50M/3jAf7swVdSVjYbfZXiLYdrV9XzFzeu4r/++HLsqgCt5FF/wQVECEFbrSejW+nNc14uWFauKj4VM+J2WOmqK+NNPWPpyHkv3f0TrGpMbXFe1lbNq6dHzYv/rhNDVHvsXNBYweXtNfzzH6xnX+8YO44NJR2byq1U5rTxmXesMdt0KEobdbVZYNqmFcJ993enuOr/fRpfMBYoPNQ3bhY5KRQzsba5gkPnxolEJX/7yH4qXDY+cl1Xyn0va69heDJofv52nhhmc2etWZl8y4YWqtx2frCrJ+nYwYlY6wzF4kSJwwLTXuvh9IgPKSVSSr790gn6xvxm/cOAN8DgRFAFoxVZsa6pklNDPr7+3DFeOz3KP/3+hdSVpx7gdHlHLO5wbszPqSEfV3TFiixddivvubyVXx08Z7qRDAYnAtSWOdNOp1OUPkocFpi2Wg/+UJSBiQD7esc4NjCJEPDj3VqjWiOlcJ0KRiuywLiJ+N+/PswNaxq4ZcPytPuubqygzGHl1Z5Rdp7QXEdXdNUl7POBK9oIRSSP7OlNWNdaZyirYTGjxGGBidU6TPHIntM4bRb+7NouXj42RO+IzwwuKreSIhuMQkm3XUsRnalJo9UiuLStmld7Rtl1Yphypy3pJmRVYwVbumr54a6ehPqJgYlgyiwmxeJBicMCY7TuPtY/wc/29XHThU3ccXUnAI/uPcOhPi9NlapthiI7WmvcXLOqjs//wYVJVdGpuKy9mkN947xwdJBNnTUpkx7++Ip2Tg35eDkuMD3oDdCQxl2lWByoOocFxvgCf/vlk4xNhXjPxlZaazxcvbKOR/b04rZbzaZqCkUmhBB8/8+uzHr/y9pqCEclPcM+tm9JXaV804VN1Hjs/GDXKa5dXa+1zpgIUK8sh0WNshwWGJfdSmOFkzf6xllW6eRavSDpto2t9Az7OHzeq9pmKOaMy9pj7S3ig9HxuOxWbtvYyq8OnmfXiWG8gTDBcFRZDoscJQ5FgBF3uPWyFjP7Y9tFTZQ7NcNONdxTzBV15U466jy47BYubkk/ivSTN66mo9bDx763h/2nxwCor1CuzsWMEociwJgnfdvlsYF4HoeN37u4GYD1ynJQzCHv39zGB7Z04LClvxxUue3cf/smgmGt4yokFsApFh8q5lAEvG9TG221HlYvS7QQ/mLrKlpr3KxsKF+gM1MsBT5+w6qs9lvVWM5X/mgDdz64G0BlKy1ylOVQBFy1so5Pv/2CpPXWGg9/sXW1mqWrKBpuXLuMu9+5liq3nZbq5J5NisWDshwUCsWsuOstK7nz2hWqOnqRoywHhUIxa5QwLH6KRhyEENuEEIeFEN1CiM8t9PkoFArFUqYoxEEIYQX+C3gnsB74IyFE6mkhCoVCoZhzikIcgC1At5TyuJQyCPwIuGWBz0mhUCiWLMUiDi3A6bjfe/W1BIQQdwkhdgshdg8MDMzbySkUCsVSo1jEIVV0K2k2oZTyfinlJinlpoaGhnk4LYVC8f+3d/+xXtV1HMefL71KA4yyzGEglxWWQhriKKYV1B+WazNMCnKKP/qDjWWUbZVra6ts/SHO7NcNC8NiFYXmzRVKhDb6oaElcKVApYaLZU5FyGDB3v3x+Xzp6/3eH9zv93vP9xx4PTbG937Or/c5fA7vc873ft7Hjk9lSQ5PA/VVvyYB/+hQLGZmx72yJIc/AtMkTZV0MrAQ6O1wTGZmxy3VXi7eaZIuAW4FTgRWRsRNw8z/H6BvBJuYAOxtPkLOBBpfptu+7Tm+1rbXSnwjja2Z7ZU9vlaWdXytLVtEfK8Fns2fp0TE8M/la+8urtof4F8jnH9Fmbfn+DoX30hjOxbja3HfHF/J4wM2j3QbZXms1IwXRjj/z0u+PcfX2vZaiW+ksTWzvbLH18qyjq+1ZYuO76iU5rHSSEnaHBEXHKvbGynH17wyxwaOr1WOr7ltVPnOYcUxvr2RcnzNK3Ns4Pha5fia2EZl7xzMzGz0VPnOwczMRomTg5mZNXByGICk/cNMf0BSoV9wSZok6R5JOyU9KemrecDgYPMvkzS24BiHPG6dJmm+pJD05k7HMhT3v6ZjdP9rIyeHCpAk4C7gZxExDTgLGA8MNVBwGVDoyVkBi4BNpBH4Ry2XlD9uuf+1TaX6n5PDICTNlXRv3c9fl3R1h8J5N3AgIu4AiIjDwCeAayWNk3SzpK2Stkj6mKTrgTOAjZI2FhmopPGSNkh6NMd0aW7vlrRd0u2S+iTdL6mwlxBLGg9cCFxHPjnzv/FvJN0t6XFJPZJOyNP2S/qCpIeAOUXFWRev+18T3P/ax8mhGqYDj9Q3RMSLpCH3HwWmAjMj4lxgdUTcRipcOC8i5hUc6wFgfkScD8wDlucrT4BpwDciYjpp4M8HC4zrA8C6iNgBPCfp/Nw+G7gBeAvwBuCy3D4O2BYRb4uITQXGWUbuf62rXP9zcqgGMUAJ89z+TqAnIg4BRMRzRQY2AAFflrQF+BXpvRyn52m7IuLP+fMjQHeBcS0ivUSK/Pei/PnhSC+ZOgz8ELgotx8G1hYYX5m5/7Wucv2vq5MbL7lDvDx5vqJTgZAKDL7sKkfSK0llzp9i4BO3U64ATgNmRcR/Jf2N/x+7g3XzHQYKua2X9BrSo5EZkoJU3DGAX9B47Go/H8gnbKe4/zXH/a9NfOcwuL8D50gaI2kC8J4OxrIBGCvpKjjyBdVy4HvA/cASSV152ql5mX3AKcWHygTgmXxizgOmdCCG/i4H7oyIKRHRHRGTgV2kq7TZSqXiTwA+TPrCsAzc/5rj/tcmTg795E5+MCJ2A2uALcBq4E+diinSMPb5wAJJO4EdpGerNwLfIT373SLpMeAjebEVwC+L+kKwdtxIx+oCSZtJV3F/KWL7w1gE3N2vbS3pWP0e+AqwjXTC9p+vUO5/zXH/az+Xz+hH0nnA7RExu9OxVEkVj5ukucCnIuL9nY6lporHsQyqeNzK2P/q+c6hjqQlpC+FPtfpWKrEx609fByb4+M2OnznYGZmDXznYGZmDZwcrCmSJkvamEed9kn6eG4/VdJ6pRo86yW9um6Zz0p6QtJfJV1c175O0mN5PT06zstV2NDa2ffqpvdK2lbkfpSdk4M16xBwQ0ScDbwdWCrpHOAzwIZcg2dD/pk8bSFptO17gW/WJYEPRcR5wAzS76gvKHRPrGra2feQdBlQ6qJ9neDkYE2JiD0R8Wj+vA/YThqNeimwKs+2ilQ2gNz+o4g4GBG7gCdIpQNqpRggDco8mXINqrKSaWffyzWPPgl8qbg9qAYnB2uZpG5gJvAQcHpE7IF0EgOvy7O9Hthdt9jTua22jvuAZ0iDp3466kHbMaENfe+LpAF9LxUQbqU4OVhL8pXXWmBZ3R3AgLMO0HbkDiEiLgYmAmNIpQbMhtRq35P0VuCNEVGagWdl4uRgTZN0EunkXB0Rd+Xmf0qamKdPJN0NQLpam1y3+CRS5c4jIuIA0Et6DGA2qDb1vTnArFx/aRNwlqQHRj/6anBysKbkMsjfBbZHxC11k3qBxfnzYuCeuvaFuVbQVFL55IeV6u/XTugu4BLKUfLASqpdfS8ivhURZ0REN6nO0Y6ImFvEPlSBq7Jasy4ErgS2SqqVQb6RVCdmjaTrSDV3FgBERJ+kNcDjpN82WRoRhyWNA3oljSFVq/w10FPsrljFtKXvFR92tXiEtJmZNfBjJTMza+DkYGZmDZwczMysgZODmZk1cHIwM7MGTg5mo0DSkto7l49y/m5XBbUy8TgHszaT1BURHqthlebkYDaAXNBtHamg20xgB3AVcDZwCzAeeBa4OiL25LILvyMN0OqVdAqwPyJuzjV8eoCxwJPAtRHxvKRZwEpS0bdNxe2d2fD8WMlscG8CVkTEucCLwFLga8DlEVH7j/2muvlfFRHviojl/dZzJ/DpvJ6twOdz+x3A9RExZzR3wqwZvnMwG9zuiPht/vwDUomGGcD6VN6HE4E9dfP/uP8KJE0gJY0Hc9Mq4CcDtH8feF/7d8GsOU4OZoPrX1tmH9A3xJX+v0ewbg2wfrPS8GMls8GdKamWCBYBfwBOq7VJOknS9KFWEBF7geclvSM3XQk8GBEvAHslXZTbr2h/+GbN852D2eC2A4slfRvYSfq+4T7gtvxYqAu4FegbZj2LgR5JY4GngGty+zXASkkv5fWalYarspoNIP+20r0RMaPDoZh1hB8rmZlZA985mJlZA985mJlZAycHMzNr4ORgZmYNnBzMzKyBk4OZmTX4Hw0n+bp1HRr0AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sorted_data['inc'][600:700].plot()"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62c488d30>"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEVCAYAAAALsCk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXl8XOV577/PjGYkjfZdtrzIxgvYBgwYY8ISCGshBbKQQtuEtqTkpmSlzW3ocpv0lra5N0uTtkkvWRpIkwDNBkmAsAdCwMZsXvGCVy3Wvkuj2d77xzlnNCONpJE0I82Mnu/no49H75wzeiWfmd95djHGoCiKoiixuBZ6A4qiKErmoeKgKIqiTEDFQVEURZmAioOiKIoyARUHRVEUZQIqDoqiKMoEVBwURVGUCag4KIqiKBNQcVAURVEmkLfQG5gt1dXVprGxcaG3oSiKklW8+uqrncaYmumOy1pxaGxsZOfOnQu9DUVRlKxCRI4nc5y6lRRFUZQJqDgoiqIoE1BxUBRFUSag4qAoiqJMQMVBURRFmYCKg6IoijIBFQdFSRGRiEEnKyq5goqDoqSID9+/k889sneht6EoKSFri+AUJdM42jlEMBxZ6G0oSkpQy0FRUsRwIIQ/GF7obShKSlBxUJQUMRwIM6LioOQIKg6KkiL8wTDDARUHJTdQcVCUFBAMRwiGDX4VByVHUHFQlBTgWAzqVlJyBRUHRUkBTiBaxUHJFVQcFCUFOJaDPxghEtFCOCX7UXFQlBQwHAhFH/tDaj0o2Y+Kg6KkgNj6hhENSis5gIqDoqSA2BRWjTsouYCKg6KkgDhxUMtByQFUHBQlBcS5ldRyUHIAFQdFSQFqOSi5hoqDoqSAWHEYVstByQGmFQcRKRCRHSLypojsFZHP2+uVIvKkiByy/62IOeduETksIgdE5JqY9fNEZLf93NdEROz1fBF50F7fLiKNqf9VFSV9xLqVtIWGkgskYzmMAu8yxpwNbAauFZFtwGeBp40xa4Gn7e8RkQ3ALcBG4Frg6yLitl/rG8AdwFr761p7/XagxxizBvgK8IUU/G6KMm/E1jlozEHJBaYVB2MxaH/rsb8McCNwn71+H3CT/fhG4AFjzKgx5ihwGNgqIkuAUmPMS8aapXj/uHOc1/oRcIVjVShKNqCprEqukVTMQUTcIvIG0A48aYzZDtQZY1oB7H9r7cMbgJMxpzfZaw324/HrcecYY0JAH1CVYB93iMhOEdnZ0dGR3G+oKPOAPximwGO9nTQgreQCSYmDMSZsjNkMLMOyAjZNcXiiO34zxfpU54zfx73GmC3GmC01NTXTbVtR5o3hQJiqonxAxUHJDWaUrWSM6QWew4oVtNmuIux/2+3DmoDlMactA1rs9WUJ1uPOEZE8oAzonsneFGUhGQ6EKSnII88l6lZScoJkspVqRKTcflwIXAm8BTwC3GYfdhvwsP34EeAWOwNpFVbgeYftehoQkW12POFD485xXuv9wDN2XEJRsgJ/MEyh102hx63ioOQEeUkcswS4z844cgEPGWN+ISIvAQ+JyO3ACeBmAGPMXhF5CNgHhIA7jTHOu+WjwHeBQuAx+wvg28D3ROQwlsVwSyp+OUWZL4YDYXxeN4Vet7qVlJxgWnEwxuwCzkmw3gVcMck59wD3JFjfCUyIVxhj/NjioijZyHAgTIXPa4mDWg5KDqAV0oqSAuLcSmo5KDmAioOipIDhQAifx02BxhyUHEHFQVFSwEjAshx8GnNQcgQVB0VJASOaraTkGCoOijJHguEIwbCx3EoakFZyBBUHRZkjjhgUet34NCCt5AgqDooyRxwxKHTqHNRyUHIAFQdFmSNOR1afprIqOYSKg6LMkajl4MmjwONmNBQhEtHuL0p2o+KgKHNkJGgN+nFSWa01tR6U7EbFQVHmSJxbScVByRFUHBRljoy5lawK6dg1RclWVBwUZY7EprIWetRyUHIDFQdFmSOxbqVozEEtByXLUXFQlDniCIHPk5fRlsMLhzrY39q/0NtQsgQVB0WZI44QFHhdFGRwQPrTD77JV586tNDbULIEFQdFmSPDgRBul+B1uzLWrdQ9FKBzcJTOwdGF3oqSJag4KMocGQlE8HnciMiYWynDxOFg2wCAioOSNCoOijJHRoKhqDspU2MOh2xx6BoMLPBOlGxBxUFR5shwIBx1JxXMwK307Fvt/Nn3X8WY9LfaOGCLw8BoCH+GCZeSmag4KMocGQmEoxbDTCyH3xzu5NHdpzjeNZzW/QEcbBuMPu4eUutBmZ5pxUFElovIsyKyX0T2isgn7fXPiUiziLxhf10Xc87dInJYRA6IyDUx6+eJyG77ua+JiNjr+SLyoL2+XUQaU/+rKkp6cKbAAXjcLjxuSUocBvxBAF470ZPW/RljONQ2QHVxPqCuJSU5krEcQsCfG2POALYBd4rIBvu5rxhjNttfjwLYz90CbASuBb4uIm77+G8AdwBr7a9r7fXbgR5jzBrgK8AX5v6rKcr8EOtWAihIsm33gN9q2JducegcDNAzHOTC06rs7zUorUzPtOJgjGk1xrxmPx4A9gMNU5xyI/CAMWbUGHMUOAxsFZElQKkx5iVjOVnvB26KOec++/GPgCscq0JRMp1YtxKQ9EyHqDgc703b3mAsU+nC1SoOSvLMKOZgu3vOAbbbSx8TkV0i8h0RqbDXGoCTMac12WsN9uPx63HnGGNCQB9QleDn3yEiO0VkZ0dHx0y2rihpw3Ir5UW/9yU5Dc5xK711qp+h0VDa9hcVB9ty6NKYg5IESYuDiBQDPwY+ZYzpx3IRnQZsBlqBLzmHJjjdTLE+1TnxC8bca4zZYozZUlNTk+zWFSWtDAdC+Dzj3Eox4nCsc4jdTX0TzhvwhygtyCNiYFeC51PFwbZByn0eGqt8FHhcdKnloCRBUuIgIh4sYfi+MeYnAMaYNmNM2BgTAb4JbLUPbwKWx5y+DGix15clWI87R0TygDKgeza/kKLMNyOBsYA0WN1ZY9NF//HR/Xz6oTcmnNfvD3HRmmogvXGHQ20DrKstQUSoKsrXgLSSFMlkKwnwbWC/MebLMetLYg57D7DHfvwIcIudgbQKK/C8wxjTCgyIyDb7NT8EPBxzzm324/cDz5j5SP5WlBQQm60ElltpOCbmcLJnJGH66IA/yIpKH6uri3g9TeJgjOFA2wBr64oBqC720qGWg5IEedMfwkXAB4HdIuLc/vwVcKuIbMZy/xwDPgJgjNkrIg8B+7Ayne40xjjvlI8C3wUKgcfsL7DE53sichjLYrhlbr+WoswPwXCEYNjEuZUKPW56hoLR75t7hhkKhDHG4ORZBEIRRkMRSgryOGdFBc8eaI97PlW09Y8y4A+xrq4EgOrifFr7/Cn9GUpuMq04GGN+Q+KYwKNTnHMPcE+C9Z3ApgTrfuDm6faiKJlG7KAfhwLPmFtpwB+k385KGgqEKc7Pi64DlBR4OHdlOT9+rYkT3cOsrCpK6f6cYLQjDlXFXva0pC++oeQOWiGtKHMgOiLUG285OG6l2Lv03uEx15KTxlpSkMe5K6xEPyfuEI4YOgZS4/oZEwfLrVRVbMUc1GurTIeKg6LMgdgpcA6xqazNvSPR9b6RMVfTmDh4WFdXQpHXzWvHe2npHeH3v/ky7/jnp2mJOXe2HGwboKrIS5VdHV1V5CUUMfSPpC91VskNkok5KIoyCVHLITaVNVYceiYTB8etlIfbJZy9vJxf7T3FI2+2MDgaIhwx7GvpZ2l54Zz2d7BtMBqMBqItNDqHRinzeeb02kpuo5aDosyBkaB1Bx5bBFfocRMIRQhHTNzdf3+MOPTHuJUAzltZQfvAKCurfPzszy4C4FD7WLO82XKkY5A1tQnEIUVuKyV3UctBUeZAIrdSbGfWlt4RvG4XgXAkoeVQWmDdvX/44tWsqi7i3WctxZvnor60gEPtA3Pa2+BoiH5/iIZyX3StqtgLaJW0Mj1qOSjKHEjkVoodFdrS64+6dRLHHKz7szKfh/eeuwxvnvWWXFtXzOE5Wg6n+iyrZUlZQXQtKg5a66BMg4qDosyByVJZAfzBMM29I6yrK8HtkoTi4KS2jmdNrSUOkcjss4qcTKlYcaj0WeLQqVXSyjSoOCjKHEjoVrIfD/hDnOr301BeSGlB3gS3ks/rJs+d+C24praY4UCYlr6xmEXX4Ch/9J87OJVkEduYOIwFtfPcLip8HrqG1HJQpkbFQVHmwFRupeNdQ4QjhqXlhZT7vPQOx1sOjkspEWtrraK1WNfSU/vbeO5AB68cS67tmCMidWX5cetVxfl0DqjloEyNioOizIGp3Epvd1gf7A0VhZQWeuIth9EgJQWTp5KutTOMYsVh+xFLFNr6k7UcRqgu9pKf545bry72quWgTIuKg6LMgeFACLdL8Ma4hxwrwvlgbygvoKzQE5fKOp3lUFHkpbrYyyF79rMxhpePdAEkXT3d2uenPibe4OBUSSvKVKg4KMocGAlEKPS44xrmOVbEYdtyWFJWSNk4y6HfH5rScgAr7uCkszb1jNBiu4nakxSHU31+6ksnFtFVF3l1GpwyLSoOijIHRoKhOJcSgM9jWQRHOoYo93koys+jrHBiQHoqywGsuMOh9kGMMbxkWw0VPg/tA8m5lVp6R+IylRyqivPp94cIhCJJvY6yOFFxUJQ5MBwIx2UqARR4XdHnGuz2F2WFHvr9oWjDO2cK3FSsrStmwB+ifWCU7Ue6qSzycsGqKtr6p7/rH7IL4JaUJxIHpxBOrQdlclQcFGUOjATCcZlKEJ+5tDRGHMIRw6A9K9qyHKZ3KwEcahvk5SNdbG2spK40n/YkAtKn+ifWODg4LTQ07pA97G7qm1PNy2xQcVCUOTB+ChzEi0Os5QBWlXQgFMEfjFAySQGcg5PO+uyBdpp7R9i2upLa0gL6/aG4MaSJcNJYE8Ycip1COLUcsoFDbQP87r/9hntfODKvP1fFQVHmQCK3Up7bFc1eGhMH6wO5byQY15F1KqqLvZT7PPzo1SYALlhdRU2Jddc/XcaS0/BvaSK3UpFaDtnEQTtj7RvPvR0Xt0o3Kg6KkgBjDA/sOBE3oCcRidxKAAUe6621dLzlMByMm+UwFSLC2tpi+kaClPs8rK8roa7U+rCfrtYhWgBXqjGHTOPtjkH+/dnDSQ9cOtY1BFg3Fvc+/3Y6txaHioOiJKCpZ4TP/mQ39z4/tSlvuZUmWgCOq8m5c491K41vujcVTtxha2MlLpdQa1sO06Wztvb7qSzyRgvyYinOz8Ob59L+Smmgvd/Pobapu+n+v1+/zf/91YGk//5HO4eoLcnnhrOX8p3fHEs6W22uqDgoSgKclta/3N066R1eJGLoGhxNmHXkswWjocK2HHyx4jA2P3o61thxh22rqwDGxCEJy6E+gdUAlkWitQ7p4e6f7OaqrzzPnT94jeP2HX8skYjhmbc6AJKe9Hesc4jG6iLuumodwXCEf3vmcEr3PBkqDkpW8k+P7ufszz/B+77xW+7+yS6eeastpa/fbbtcjncNs6e5P+ExB9oG6PeHojOgYynwuPG6XVTb/v1Yy2H8oJ+puGBVJUVeN5efXgtAhc9LnktoSyLmkCje4FBdkq+WQxp4s6mPxiofz+xv54ov/ZqvPX0o7vk9LX1RUU5aHLqGWFVVRGN1ER84fzk/2H6CE13DKd/7eKYVBxFZLiLPish+EdkrIp+01ytF5EkROWT/WxFzzt0iclhEDojINTHr54nIbvu5r4ldVioi+SLyoL2+XUQaU/+rKrnEGyd78bgFt0v4+Zut/OWPd6f09buHxgJ/v9jdkvCYl962CtO2nVY14blCj4sl5QW4XFbldJHXHW3bPX7Qz1Rsaihj799fy6rqIoCoa6l9mlqHU/2JW2c41BTnJ92GQ0mOjoFROgdH+eCFjfz6M5dx+em1fOWpg5zsHvsgf3p/e/RxcxLiMOAP0jkYoNH+///kFWvxuF38+mD7NGfOnWQshxDw58aYM4BtwJ0isgH4LPC0MWYt8LT9PfZztwAbgWuBr4uI4/j8BnAHsNb+utZevx3oMcasAb4CfCEFv5uSwwyOhti8vJyHPnIhH7l0NR0Do9EOqamgx3YrbVlZwS93JXYtvXykixWVvmhGUixnNpSxbdWYaIhItIXGTGIOiagpLZjS7zwSCNM7HIxr1T2e2tKCpOollOTZ32pZmGcsKaG2tIDP37ARAb6//UT0mGfeaue8lRX4vO6kxOFYpyUsq6qtaX51pQW88JeX88ELG1O+//FMKw7GmFZjzGv24wFgP9AA3AjcZx92H3CT/fhG4AFjzKgx5ihwGNgqIkuAUmPMS8Z6p90/7hzntX4EXOFYFYqSiMHRUHRQzvJK643T3Js6U7t7OIDHLXzg/OU09Yywq6kv7vlIxLD9aDcXrp5oNQB8/sZNfOH9Z8WtlY8Th+JZikNtydR3/VMVwDnUlebTNRQgGNYWGjPFGMOH73uFf3nqYNz6PlscNiwpBaxMtSvPqOOhnScZDYVp7/ezu7mPd51ey9LywqTcSkftuIVjOcBYEWO6mVHMwXb3nANsB+qMMa1gCQhQax/WAJyMOa3JXmuwH49fjzvHGBMC+oAJ7zoRuUNEdorIzo6OjplsXckxBvyh6Ifr8krrDvlkd3I+3GToGQpQ4fNyzYZ6PG7hl7tb457f19pP30iQbadVJv2apVFxCFLoceOZZNDPdNSW5E+Zytpqf+hM5VaqLbGeU9fSzHl6fztP7W/nJ681x63vb+1naVkB5fa0PYAPXriS7qEAj+0+xbMHLFfQu06vpaG8kJbe6S23Y52WOKysLJrmyNST9NUpIsXAj4FPGWMSR+jsQxOsmSnWpzonfsGYe40xW4wxW2pqaqbbspLDDMZ0NV1eYVkOJ3uSsxz6hoN87elDhKdoR9A9FKCyyEuZz8Mla2smuJac9tkXrq5Oes+xbqXZupTAci30DAcnbZyXaALcxNew7j6TnQ2hWEQihi89aVkMJ7qH41xD+1v7OcO2GhwuOq2aVdVFfO/l4zzzVjtLywo4vb6EpeWFSbqVhqgvLZhQhT8fJCUOIuLBEobvG2N+Yi+32a4i7H+dCEkTsDzm9GVAi72+LMF63DkikgeUAcmNu1IWHaOhMIFwJOpWqinJJz/PFRf4m4pH97Ty5ScPsqupd9JjeoYtywHg+jOX0Nw7wusnx45/+UgXq6qLprw7H09UHEan78g6FU46a8ckqaiOW2myVFbrNaznkm3/rVg8vvcU+1v7+cilqwHYbt8k+INh3u4YmiAOLpfwBxes4NXjPTzzVjuXn16LiNBQXkD3UGDaONmxriEa7XjDfJNMtpIA3wb2G2O+HPPUI8Bt9uPbgIdj1m+xM5BWYQWed9iupwER2Wa/5ofGneO81vuBZ0yy5YOzwMkWUbKTwXEBXRFhWUVh0m6lo7apPpVZ32VbDgBXbazD53XzxV8dIBwxhO14w7bVybuUYLzlMH2m0mTUlk5d69DSO0KFzzPl3WbdNK+hTCQcMXz5yYOsrS3mL65ZT1mhJ2pBHmobJBwxbFhaOuG895+3jPw8F8Gw4V12SrJT/xI7IzwRx7qGo5lq800ylsNFwAeBd4nIG/bXdcA/A1eJyCHgKvt7jDF7gYeAfcDjwJ3GGEcePwp8CytI/TbwmL3+baBKRA4Dd2FnPqWDb71whHf+3+dUILKYaEA3pnHdsgpf0m6lIx2OOEz+xuwZClBRZH2AlxZ4+F/v3sBv3+7i/z3/Nvta+hnwh6KFacniTIPrH5mr5eC00JjEcujzUz+FSwmsmQ4uUcthJjzyZjOH2we566p1eNwuLlhVycv26NaxTKWJ4lDu83LT5gaK8/N4x2mWG3Kp/f/T3DP5Ndg3EqR7KEBj1cKIw7RXqDHmNySOCQBcMck59wD3JFjfCWxKsO4Hbp5uL6lg66pK/uGX+/ney8f5s8vWzMePVFKM0/Y6VhyWVxbyxsnJ3USxHO20GplN5vMNRwy9I0Eqi8ayQn7v/OU8f6iDLz9xkP2tVnuEyTKVJqOs0EPEWDGBZRWzdxVE3UqTpLO29vmnzFQCcLuE6uKpA9vKGP3+IF964iAbl5ZyzcZ6AC48rYon9rXR3DvCvtZ+fF43KysT/7/+3Q0buPPyNTFtVWzLYYobFCcY3ZjBlkNOcdayct65roZvvXCU4UBoobejzIJEjeuWV/js6uOpLcJQOMIJOzYxmTj0jQQxBip9Y68vIvzTe86irrSAn7/ZwuqaImqn8OknwqmSbh8YnZPlMN1d/6n+6cUBrMB2MoODFjvGGP7ioTc51efn72/cGC1sdCzH7Ue62Nfaz/r6kuhz4/F581hRNSYc9WUFuGTqQjin4V4mu5Vyjk9csZbuoQA/iClOUbIHx3KI/YB1ah2mC0o3944QDFvhrMnu2rrtAriKIm/cepnPw7/cshmXWFkoM6UsRmzmIg7OXX+iKumh0RDdQ4HonelU1JXmq1spCb75whGe2NfGZ3/ndM5bORZnWl9XQrnPw0tvdyXMVJoKj9tFXWnBlOJwtHMIEVgxiTWSbhalOJy3soKL1lTxH78+Mu3QFCXzcOJFcW4lJ511mqD0EdtUP72+ZFJx6LHbdFeOEweA8xsr+emfXcRdV62b8b4dywGSa7o3FbWl+bQlcCvN5G6zpkSrpKdj+5EuvvD4Aa47s57bL14V95zLJVywqpLH95xiwB+KFr8lS8M0hXDHOodYWlaYsLPufLAoxQHg4+9aS+fgKA/syCzrIRiO0KXdMqckGnMoiI85ADRNE5Q+agejL15TTc9wMKFrMWo5+CaKA8DZy8snWBXJEC8Os7ccwApKJ7IcnEysZMRBq6SnxhjDpx98g5WVPr7wvrNI1LRh2+oqBuzrcSaWAzBtrcPRruEFS2OFRSwO21ZXsbWxkv/49ZGMylz6r5ePc/kXn9M37BQk6k1UVuihJD9vWrfS0c4hSgvy2NRQBiR2LTnikMhymAuptBwmcwk54pdMhoszCEirpBMzEgzT0ufn5i3LJ/3/utBuuihiWaMzYWl5Iaf6/JMWYx7rHFqwTCVYxOIA8Jlr19M5OMqf3r8zY9xLRzqG6PeHtNf+FAyOhvC6XeTnjZnbIkJDRSEnp0gNBEscVtUUR/PMmxPUOkxnOcyWWHFINANiJtSUFNA1NEpo3E3E0a4hlpQlV1HrZD1pxlJinJGcsf9v41lXW0KFz8PKSh9F08wEH09DRSHBsEn4Xu8ZCtA3ElywYDQscnE4v7GSL33gbLYf7eZjP3gtI+7WnQtlupbMi5kBfzBh07rllb7p3UqdQ6yuLooGbBPlmfcMBSj0uFPessDndZNnZ7PMOeZQko8xTJjJcHQGd5uO5aBB6cQ44lDum/z/yuUS/uyyNdz2jsYZv36DPW+jKcE1eLjDSrdWy2EBuXFzA39/4yae2t/OZ/77zaTnuqYLZ+i7vmEnZ9AfigtGOyyv8HGye2TS/0N/MExz7wirqouoK8nH7ZLEbqXhQMpdSjDWthtSEXNwxoXG3/Uf6xxiVU1yHyjTVVovdvqGp7ccAP700tX88UWrpjwmEVPVOjz7Vjtul3DeyomDpOaLRS8OAB/ctpK7rlrHz95o4dcHF7bba9RymKc5sdnI4GjixnXLKwsZCYajIz7HE5vJk+d2UV9akPCN2TOUHnGAsQ+aZAb9TIVz1x9bp9A7HKBnOMiqJO82q4q0SnoqknErzYWGScTBGMPje06xbXXlrBIfUoWKg81H3rmacp9nQhve+UbdStMzMIXlAJPXOjjBWsePu7S8gKaElkMwbW9Kp9ZhrpbDSrug6mDMMPuZZCqBVS9RM03778VMusWhpMBDSUHehIylg22DHOkc4tpNS9Lyc5NFxcEmP8/N7561lF/tPbVg2UujoXB0vrDezU3OZC2vo4VwkwSlj4z78Jwsz7xnKBBXHZ1KnA+a2Q76cSj3eVlXV8zOY2PNi6PikKRbCayU2GSqpB/b3cr1X3thyjbnuYYjDqVpEgdIfA0+vucUInDNhrq0/dxkUHGI4b3nNjAaivDY7lML8vO7Y9whk/XNUeKnwMWyrMIZ+jOJ5dA5RF1pfjSrZLJUQqvpXvrcSnMZ9BPLlsZKdh7vIWLv/2jnEC4Zs6CSIdkq6Uf3nGJvS3/0A3Mx0DcSRARKZpiFNBMaygsnBKQf29PKeSsqZtyeJdWoOMSweXk5q6uL+PFrTdMfnAY6ByxxcMnkHTcVJ+Yw8W6uKD+PyiJvwuwPsNNYY1wuS8sLCUVMXHxnNBRmYDREZYrTWB3OWlbOOSvKU/Ja5zdWMOAPcbDdci0d7RxieaUPb17yb+tkZ0m/fqIHYNGJQ1mhZ9J+Salg/LjQY51DvHVqgGs31aftZyaLikMMIsJ7z21g+9HupAfHpBIn3rC6plgD0lMwGDMidDzLKwonTWe1xKE4+n20p37Mm7PXzlBJl+Vw+8Wr+MGfbkvJa22x+/y8csz64J5JGqtDbcn0VdIdA6NRwe0dThzsz0UccUgnDRWF9PtDHLFTVx/fa3ktVBwykJvOscZa/+z1+Q9MO+JwxpJSOgcDi8q/myz+YPwUuPEsq/QlFPbe4QDdQwFWxbQjcLJFYgvh0lUdnQ6WVRRSX1rAzmPdGGOsNNYZFk0lUyUd2wq9dxFaDunk8vW1lBTkceO/vcjP32zhsT2nOLOhbE4t3VOFisM4llX42La6kp+83jzvNQ9OQdOGJaWEIyYuBqFYJOrIGsumpWUc6xrmlWPxU2bHMnnGLIdEhXA9WSQOIsKWxgp2HuuhY2CUoUCY1TMIRkNys6QdlxKM5f4vBuZDHNbXl/DoJy5hbV0xH//h67x5sjcjrAZQcUjIe89dxtHOIfY098/rz+0aHKXQ46bRTlNU19JExo8IHc9t71jJ0rIC/vZne+JaS7x63PqAi7UcivPzKCv0xLmVuqfoyJqJnN9YSXPvCC++3QnMvKI2mVnSr5/ojVpZi8qtNBxMa6aSw/JKHw9+5ELuvPw0akvyueHspWn/mcmg4pCAc+2A4RF7Yth80Tk4SlWxd6xyVdNZJzA2BS7xm9bnzeNv372Bt04N8F8/HJMJAAAgAElEQVQvHwfg5SNd/J/HD7B1VSWrYywHmBgQ7ElTX6V0saXRqqD9751WEsVM3UrTVUmHI4ZdTb28c30NoG6ldOFxu/jMNaez46+vjKZkLzQqDglw5u+e6pvfO/fOwQDVxfnRu7kOzViaQH+CWQ7juXZTPZesreZLTxzkt4c7ueP+nayo8nHvB8+bkHnSUB4/cKV7aPp+OpnE6fWlFOfn8du3u/C6XUkN+YmlqshqIzLZjcih9gGGAmHOb6ygpCAvGrDPdYwx8yoOmYiKQwKK8/Mozs+jdd7FYZTq4nxqJumbo0zvVgLLF//5GzbiD4X5/W9tJ9/j5rt/fD7lCayBhnE99buHRiktyEtJHcJ84HYJ59r9d1ZW+XDPMO3SmirnnTTm8PoJKxi9eXkF5T4P/YvEchgOhAlFDOUqDsp46ssKZmQ5DI2GJrRPnimW5eClwOOmtCBP3UoJmC4g7bC6ppiPv2stpQV5/OcfnT9p9sfS8kIG/KGoRdI9HMyaeIPD+bY4zLa981SzpN840Uu5z0NjlY/yQu+icSulu3VGNjCtOIjId0SkXUT2xKx9TkSaReQN++u6mOfuFpHDInJARK6JWT9PRHbbz31N7LFKIpIvIg/a69tFpDG1v+LsWFJWwKkZ9Jz5na++wL8+c3jWPy8SMXQPWZYDOMVJKg7jGYs5TF+1+okr1rLzb66KDvZJxPiMpXRWR6eLLY1WvcNsxWFpWSHH7aaE43n9ZA/nLC+PdpRdLAFpFYfkLIfvAtcmWP+KMWaz/fUogIhsAG4BNtrnfF1EnKb43wDuANbaX85r3g70GGPWAF8BvjDL3yWl1JUmbzmMBMKc6B5m+9GuWf+8nuEAEQNVxV775+erWykBzhS4ZHsTTVctfPayctwuiQavu4cCVGWZOJyzopyzl5dz6bqaWZ2/pbGCY13DE3r8DPiDHGof5JwVlmVS5vMsGsuhN8l23bnMtOJgjHke6J7uOJsbgQeMMaPGmKPAYWCriCwBSo0xLxmreOB+4KaYc+6zH/8IuMKxKhaSJWUFtA/4k3IVOR/i+1sHZl0b4bSZjloOJQXqVkrAgH/iFLi5sKLKxwe3reQHO06wp7mPnuFA1mQqORR43Dx850VctKZ6Vuc75714uDNufVdTH8ZYbWUAygs9i6bOYT6a7mU6c4k5fExEdtluJ2ciRQNwMuaYJnutwX48fj3uHGNMCOgDqhL9QBG5Q0R2isjOjo70zl2oLysgkmDSViIcf23fSJCWKayNvuHgpIVtnbYQjIlDPu39ows+fCjTGBwNzrnd9Xg+feU6KnxePv/zvXSncZZDprK+roSqIi+/fTve8nVqQ852xMG2HBbDNdmvbqVZi8M3gNOAzUAr8CV7PdEdv5lifapzJi4ac68xZosxZktNzexM6GSpt9sKtPZNPZMY4qtL97dMXjh310NvcN1XX0g4M7Zj0BEH64OppiSfQDiyqBqdJcNUfZVmS5nPw2euWc8rx3oYDUWyLuYwV1wu4cLTqnjxcGf0g98Ywy93tbJ5eXn0A7K80Es4YhgKZMa89XSSzIjQXGdW4mCMaTPGhI0xEeCbwFb7qSZgecyhy4AWe31ZgvW4c0QkDygjeTdW2qgvs8QhmbhDrDjsa00sDpGI4ZVj3Zzq9/PxH7w+wV3ljAeNDUiDFsKNZ7JBP3PlA1uWs6mhFCBtHVkzmYvWVNM+MMrbdgO4vS39HGgb4H3njb1tHZFYDEHpvpEgbpek5VrLFmYlDnYMweE9gJPJ9Ahwi52BtAor8LzDGNMKDIjINjue8CHg4ZhzbrMfvx94xmSA3brEKYRLImOpfWCU/DwXK6t87J9EHI51DdHvD3HJ2mpeOtLFl548GPd85+AobtfYjOHojGDNWIpjYJJZDnPF7RI+f8MmvHkuTqtduKHuC8VFpzlxB8u19KNXm/C6Xdxw1lgrB2eK3WIohOsbCVJakEcGhD8XjGnfZSLyQ+AyoFpEmoC/Ay4Tkc1Y7p9jwEcAjDF7ReQhYB8QAu40xjg26EexMp8KgcfsL4BvA98TkcNYFsMtqfjF5kqFz4M3z5WU5dDe76e2NJ+NS0vZO4lb6c0mq5job67fwHd/e4xvPPc25ywv5+qNVpOtrkErS8ap4J1sgPxiZ9AfmnEVcLKct7KC3Z+7OmXB7mxiRZWPZRWFvHi4k1u3ruDhN5q5akNdVBCAaEHYYnB1LvbqaEhCHIwxtyZY/vYUx98D3JNgfSewKcG6H7h5un3MNyJCfWlBUlXSbf2j1JUUcEZ9KY/uPpVwUtmbJ/vwed2sqS3m7353A7ube/nbh/dw5Rl1uFxi91XKjx6vbqXEWIN+0mfqL0ZhcLjotGoe29PK0/vb6BkO8r7zGuKedyrMF4Pl0KvioBXSU1E/rhAuGI7wmf9+kwOnBuKOaxvwU1dawIalls/6rQSupTebetnUUIbbJRR43Nx+8Sra+kejFkXnUCAajAaryMvndatbaRwD/uCi9gOnk3esqaLfH+ILj79FdXE+l66NT/pwgrO9I4sj5rCY01hBxWFK6scVwu1p7uO/X23isT2tcce1949SW5rPGUsscRgfdwiEIuxt6Y/mi4M15MPtEp7Y1wZYqaw1MZYD2Oms6laKYoyxrLI0Wg6LmXfYcYdjXcO855yl5I3rL1W2iNxK/SPBhL24FhMqDlPgtNBw4uNOEzJncAxYbo7B0RB1pQUsKSug3OeZkLF04NQAgVCEs5eNiUO5z8sFqyp5cl8bxphou+5YtBAuntFQhGDYpNWttJipKclnfV0JQFyWkkOBx01+nmtRFMJZMYfFfZ2pOExBXWkBgVCEHvvN4IxLPBYjDk4f/LrSfESEM+pL2dca73Z6w3Ydnb08vsfPVRvqONw+yJ7mfkZDkWgaq0NNaf6U4xsXG9Gme+pWShu/f8EKfvfspZxeX5rw+XKfJ+djDtqu20LFYQqWlMUXwjnicKRzKGpNONXRdfYMhg1LSzlwqj9u/vObJ3upKvJGp2k5XLWhDoAfvnICIC4gDU6VtLqVHGbaV0mZObe9o5F/vfWcSZ+3OrPmdsxhKBAmHDEqDgu9gUwmthCua3CUE93DLC0rYMAfirbBcGICTnbRhiWl+IORONfTmyd7OdvubBnLsgofG5aU8vDrzQBxAWmw3EpDgXD0jnmx48xymGwKnJJ+yhaB5eAU+ak4KJMSFYd+f9RquPEcK73vmN3iuC3GrQREg9JO3GFwNMThjsG4eEMsV2+si7YjGO9WcmIQ3Un0d1oMDIxaH0oac1g4ygo9OR+Q1nbdFioOU1BTnI9LLMvhjZO9uF0SHf59pMMRh1EKPe5oeuWa2mI8buHlI10YY9htd7YcH29wcFxLMFEcnDYOPYugXUEyjFkOKg4LRXlh7lsOY+KwuLOV9F02BXluF7UlViFcW7+f9XUlrK0tJs8lUcuhfWA0GowGa37AVRvq+MH2E7T1+aPDwiezHDYsKY2OqhzfDbSiyLpzWczi8L2Xj7O0rIArzqiLxhzUclg4yn25bzloR1YLfZdNQ11ZAS29I+xu6uOGzVbu94pKXzSm0Nbvj8YbHP711nM5d8VRvvjEAfxvRVhZ5Zu006eI8J5zGnh0d+uEwTSLqSJ1Mr761EHyXC4uXVczoylwSnoo93kZCYbxB8MUeLKrmrxrcJT//Yt9uOweZjUl+fzxO1ZR6I3/PaKWwyLuyAoqDtOypLSAXx/sYCQYjhaxNVYXcbRzGLBSWc8cZxW4XcKHL1nN1Rvq+Ydf7mPzisRWg8Onr1rHx69YM2G9YpG7lfzBcHSexqO7W8fEQS2HBcO5m+4fCWadOPz27S5+9kYLtSX5DNuJHqUFHv5w28q44zTmYKExh2moLytgJGgFjM+xP+Qbq4o4ZqezWn2V8hOeu6LKx70f2sKfXTbxgz8Wt0sS9vQpK/QgYs01XozEVqd/64Wj9PuDePNSNwVOmTljLTSyz5p1xqA+/efvZPfnrqaxysev9p6acFzvsNWuu8i7uK8zFYdpcDKWSgryWF1dDMCqmiJGgmEOtw8yEgxTN86tlCqcFt49i9St1GLXl1x/1hJ2N/fx6wMdWgC3wIzNdMi+a7K1z09JQR4lBR5EhGs21vPS210TYihOAdxibtcNKg7T4hTCbV5eHm2nvarK6vf/8hGr931taWLLIRVU+LyL1q3U2mtZDh+7fA3lPg9vnRpQl9ICU17oxMGy75ps7h1hadlYIerVG+sJRQzPvtUed1zfSDDannwxo+IwDc640NimeatqbHE4ag2sS5flAIujXcFkOJXpq6qL+MMLLL+wZiotLI5bKRszllr7RlhSPvZePWd5ObUl+RNcS9qR1ULFYRrW15ewpraYqzfUR9eWlBaQn+di+xFLHGoniTmkgsVsObT0+aks8lLgcfOhC1ficS/usY2ZQFk2i0OvP25QlMslXLWhjucOdOAPjs3F7te+SoCKw7SU+7w8ddc7OXPZWBGbyyU0VhXROWj1VRqfyppKKnzeRRuQbukdYal9p1dbWsBfXXcGH9iyfJqzlHRSkp+H2yVZZ836g2G6hgIsLYt/r16zsZ6RYJgXDnVG17TpnoWKwyxprLaK24rz89J6N1vhW7wB6dZef3SWN8AfX7SK9547sZW0Mn+ICKUFeVnXfM+Z6Dh+xOy21VWUFOTFuZZUHCxUHGZJY7UVd0hnMBqgomis6Gix0dI3MuFOT1l4yn3erLMcnDTW2JsNsDoaXHlGHU/vbyMUjjAaCqs42KgDd5astsXBadWdLqJ55cNB6ssWT9714GiIAX+IJePu9JSFJxub7znisLR84vv1mo11/PT1Zrb+49P0DAcwZmKH5MWIisMsabTTWevSbTnEVEnXL6K76Nbond7i+Z2zhXKfJ9qyPltw3EqJ3kOXra/lfecuw+2y3E7LK3z8zpn1E45bbEwrDiLyHeDdQLsxZpO9Vgk8CDQCx4APGGN67OfuBm4HwsAnjDG/stfPA74LFAKPAp80xhgRyQfuB84DuoDfM8YcS9lvmCacdNZ0prFCjDhk2ZtxrjRH7/TUcsg0ygs90a7E2UJL7wjVxfkJq+sLPG6+9IGzF2BXmU0yMYfvAteOW/ss8LQxZi3wtP09IrIBuAXYaJ/zdRFx/je+AdwBrLW/nNe8HegxxqwBvgJ8Yba/zHxSU5zPn16yiuvPWpLWnzPWmTW7zPi5MlkAUVl4ygo9WVcE19LnT+hSUiZnWnEwxjwPdI9bvhG4z358H3BTzPoDxphRY8xR4DCwVUSWAKXGmJeMNV/z/nHnOK/1I+AKyYK6dRHhr6/fwFmTtOJOFYu1+V5r7wguYdK+VcrCUebz0u8PEY4Y+v3BaEp3JtMyrjpamZ7ZxhzqjDGtAMaYVhGptdcbgJdjjmuy14L24/Hrzjkn7dcKiUgfUAV0Mg4RuQPL+mDFihWz3Hp2MRaQXlzi0NLnp7akgDy3JtRlGk5riY1/9zj+YASXwLN/cRkr7ThcpmGMobV3hEvWVi/0VrKKVAekE93xmynWpzpn4qIx9wL3AmzZsiXhMblGfp4bn9e9CN1K8a0OlMzhyjPq2NfaT3mhB2+ei68/9zavHu/JWHHo94cYCoRpUBfljJitOLSJyBLbalgCOJ2rmoDYEtZlQIu9vizBeuw5TSKSB5Qx0Y21qFmMVdItvX422PO4lcxiRZWPL95sBXDDEcN3XjzKnuZ+3nvuAm9sEiarcVCmZrY2+yPAbfbj24CHY9ZvEZF8EVmFFXjeYbugBkRkmx1P+NC4c5zXej/wjB2XUGwqijyLKuZgjIlrnaFkLm6XsGFJKXta+hZ6K5PiNHBUS3RmJJPK+kPgMqBaRJqAvwP+GXhIRG4HTgA3Axhj9orIQ8A+IATcaYxxSns/ylgq62P2F8C3ge+JyGEsi+GWlPxmOYTVfG/xuJV6hoOMhiJ6p5clbGoo4yevNROJmGhb+0yi2W79rm6lmTGtOBhjbp3kqSsmOf4e4J4E6zuBTQnW/djioiSm3OflZPfwQm9j3piqmlXJPDYtLeP+l45zrGuI1TXFC72dCbT2jpDnEqqLNfNtJmgqSBZQucia7zk1Dmo5ZAcbG6zY0J6W/gXeSWJa+/zUlxXgzkCrJpNRccgCyn1e+kaChMKRhd7KvBANIKrlkBWsrS3B63axtzkz4w7jJ8ApyaHikAVUZPGAldnQ0jeCxy1UF6kbIBvw5rk4fUlJxgalNS16dqg4ZAEVRU6V9OIQB2eOQyYGN5XEbFxaxp7mfjIt0TASMZzq82sbllmg4pAFlPuyd6j7bGjtG9FurFnGpoZS+kaCNPWMpPXnnOweZldTb9LHdw6OEgwbnQsyC1QcsoBK3+KyHFp69U4v29i01BqjuyfNcYd/emw/t977ctJFodrdd/aoOGQBTn+lxVAl3d7vp7l3hLV1mZcSqUzO+voS3C5Je9zheNcwQ4Ew33zhSFLHt/Rq5ttsUXHIAsZiDrkvDs8d7ADgsnW10xypZBIFHjdra4vZ05zedFbHErjvt8eSGjj0ZlMvXreL1TWZ2fcpk1FxyAKKvG48blkUbqXnDrRTV5rPGUtKFnorygw5s6GMPc19aQtKD42G6B0O8t5zGxgOhrn3+emth+1Hu9m8vJwCz+IZsZsqVByyABGxh7rntuUQCkd44VAn71xXQxaM9FDGsamhjK6hAKf6/Wl5fcdqeOe6Gn73rKXc/9IxuqaYJTE4GmJPcx9bV1WmZT+5jopDllDp80bdSpGIYeexbiKRzEobnCuvnehlwB/isvXqUspGzrC76B44NZCW12+2M6GWVRTyiSvWMhIMc+8UsYfXjvcQjhgVh1mi4pAllPs89AxZbqVv/+Yo7/+Pl/jey8cXeFep5bkD7bhdwsU6lCUrWV5pBX3Tlc7aZFsODeU+1tQWc8PZS/neS8fpm8TduuNoN26XcO7KirTsJ9dRccgSKmzLoXsowNeeOQTAl588mFMZTM8e6OC8lRWUFngWeivKLKgrKcDjlrSJQ3OPVTlfa4+O/cilpzEcCPP9HYlvknYc7WZTQxnF+ameabY4UHHIEqyZDkG+9vQhhkZD/Out5zDgD/IvTx1c6K2lhLZ+P/tb+7lsfc1Cb0WZJS6X0FBeSFNPejoIN/eOxFXOb1haysVrqvnui8cIhOL7jvmDYd442csF6lKaNSoOWUKFz0v30Cj/9fJxbt26gt89eyl/cMFK/mv7ibT5eOeTXx+wUlgv13hDVrOswpdGy2F4wkyGD1+yivaBUR55syVu/c2TvQTCEbY2qjjMFhWHLKHC5yVirHzyT1+1DoC7rlpHcX4ef/+LvRnX02amPHewnfrSAk6v1xTWbGZZRWH6xKF3hIaKeHF457oa1teV8K0XjsS9B3Yc7UYEzldxmDUqDlmCUwj30ctOiw4tqSjy8qkr1/Li4S52Hu9ZyO3NCX8wzAsHNYU1F1hWUUjn4CgjgfD0B8+AQChC+8DoBMtBRLj9klW8dWqA3xzujK7vONbN+roSynwav5otKg5ZwuXra/j0leu4/eJVcetXnlEHwNGOoYXYVkr41d5TDIyGuP6sJQu9FWWOLKvwAdDcm9q4Q2vfCMYwwXIAuHHzUmpK8vnXZw7TOxwgGI7w6vEejTfMERWHLKGqOJ9PXrl2QqVnVbFlUXQOTV4MlOn8cMcJllcWcvEaTWHNdpbZH94nU+xaitY4JGigl5/n5s7LTmPH0W62/uPTfPi+nQwHwmxdVZXSPSw2NMcry/F58/B53XQOZGdK65GOQV4+0s1nrlmv8xtyAMdySHXcIVrjkMByAPiji1Zx/qpKHthxkp+93ozHLVr8NkdUHHKA6uJ8OqdoI5DJPPjKSdwu4ebzli30VpQUUFuSb9c6pNat1NI7gsjU3VU3Li3jf99Uxt3XnU7nQICaEp0kOBfm5FYSkWMisltE3hCRnfZapYg8KSKH7H8rYo6/W0QOi8gBEbkmZv08+3UOi8jXRKOSM6K62EtXhrqVjDG8cbI3YTbVaCjMf7/axJVn1FJbqsNYcoGxWofUu5VqS/Lx5k3/keXz5rGiypfSn78YSUXM4XJjzGZjzBb7+88CTxtj1gJP298jIhuAW4CNwLXA10XEcaB/A7gDWGt/XZuCfS0aqovzM9at9Nu3u7jp31/kZ280T3juyX1tdA8FuHXrigXYmZIu0lHr0Nw7MiFTSUkv6QhI3wjcZz++D7gpZv0BY8yoMeYocBjYKiJLgFJjzEvGur28P+YcJQmqMtit9NT+NgDuf2lii4MHdpykobyQS9ZqVXQusbyykOYUu5WsGge1BuaTuYqDAZ4QkVdF5A57rc4Y0wpg/+uUvDYAJ2PObbLXGuzH49cnICJ3iMhOEdnZ0dExx63nDjXFXrqHA4QzsEvrcwc68LpdvH6iN26E5J7mPn5zuJPfO385bg1E5xTLKnx0DgZSVusQiRhae/1qOcwzcxWHi4wx5wK/A9wpIpdOcWyiTwAzxfrERWPuNcZsMcZsqanRu02H6pJ8jCGpyVjzybHOIY52DvGxd62hwOPiv+wussYY/v7n+6gs8nLbOxoXdpNKynHSWVNV69AxOEogHJk0U0lJD3MSB2NMi/1vO/BTYCvQZruKsP9ttw9vApbHnL4MaLHXlyVYV5KkqsjKysg019JzB6z/+hvOXspNmxv42RvN9I0E+eXuVnYc6+Yvrl5PWaFWsOYaqa51cOIXDeWatDCfzFocRKRIREqcx8DVwB7gEeA2+7DbgIftx48At4hIvoiswgo877BdTwMiss3OUvpQzDlKElTbhXBdg/NrOTz7Vjs3/vuL+IOJ3QfPHuhgVXURjdVF/OG2lfiDEb6//Tj/9OhbnLGklN87f3nC85TsJtW1Ds0xcxyU+WMudQ51wE/trNM84AfGmMdF5BXgIRG5HTgB3AxgjNkrIg8B+4AQcKcxxvlU+SjwXaAQeMz+UpKkumRhLIcfvdrEmyd7ee5AO9duim99MRII8/KRLn7/AisTaVNDGeesKOeLvzpAxMAXbz5bYw05Sk1xPl63a8a1DqFwBBGZcF041dHqVppfZi0OxpgjwNkJ1ruAKyY55x7gngTrO4FNs93LYqd6AdxK4YjhxbetRmc/39U6QRxePtLFaCgSN/Lzg9tW8vqJXq7dWM+Fp2lrg1zF5RIaZtGd9aavv8iWlZV87oaNcevNvcOUFXp0aM88o3/tHKC0MA+v20XnPLqV9jT30TscpK40n2f2tzMcCOHzjl1Ozx1op8Djimt+9u6zlnKsc4g/2LZy3vapLAzLKgpp6k7ecugZCrCnuZ9jncP85bWnU+gd6yF2tHNIM5UWAG28lwOICFXF3nm1HJz2yH9z/QZGgmGefWsstdgYw7MHOnjHadVxjQK9eS7uuno9dVoNnfPMdK7DbjvNeXA0xJN2bQxAU88wL73dxeWna3bifKPikCPMtzi8cKiDDUtKue7MJVQX5/PL3WMJZm93DHKie5jLdeTnomVZhY+uoQDDgVDc+tBoiJePdE1Yd8Shujifn742Vvb0g+0nAPj9C9TanG/UrZQjVBfnz1u20nAgxKvHe/iTi1bhdgnXnVnPQztPMjQawu0S/uK/d+HzurlqQ/287EfJPJx01leO9RAMRdjb0s+Lb3fy+okegmHDJ961hruuXh89fldTL41VPn7nzCXc+/wROgZGKS3M48FXTnLlGXXqVloA1HLIEeazM+v2I90Ew4aL11rzF64/cwn+YISn9rfxvx7ewxsne/nSzWdTX6buo8WKk85623d28OH7d/KVpw4yHAhx+8WrOb2+hF8fjO9wsLupjzOXlfPecxoIRww/f7OFR3e30jUU4EMXNi7Ab6Co5ZAjVBV76RoMYIxJ+6jNFw514s1zRefzbmmspLYkn3/45X46Bkb5+LvW8Dtn6lS3xczm5eV85pr1lPs8nF5fyvr6kmi2UYHHxVefPkTvcIByn5eOgVFa+vz8cUMZa+tK2NRQyk9fbybPLayuKeKiNZrZthCo5ZAj1BTnEwhH6PeHpj94jvzmcAdbGyujwWbLtbSEjoFRrjyjlk9fuS7te1AyG7dLuPPyNfzBBSs5b2VFXBrqJWurMcbq2AtEe26duawMgJs2N7C7uY/XT/TywW0rda74AqHikCNUF89PrUNbv5+DbYNcsjZ+pOefXrqa2y9exVd+b7NOdFOm5Oxl5ZTk5/HCIcu1tKupDxGrUBLghs1LcQkUety8T4dALRjqVsoRouIwMMppNcVp+zkvHLJSWC8eJw4N5YX87bs3pO3nKrlDntvFhadV8cKhTowx7G7u5bSa4qh1UVtSwO0Xr6KiyEtpgfbeWihUHHKEKqe/Upo7sz61r42aknzOqC9N689RcptL1lbzxL42jncNs6upj4vXxN9s/PX1eqOx0KhbKUeYD7fS4GiIZw+0c92menUdKXPCGfD049eaaB8YjcYblMxBxSFHqCzyImK5ldLF0/vbGA1FuP6spWn7GcriYGWVj2UVhdz322MAnKXikHGoOOQIbpdQ6fPSkcZCuF/uaqWuNJ8tKyvS9jOUxYGIcMnaavr9IVwCG5aoOGQaKg45hFUlnR7LYcAf5LmDHVx35hJ1KSkp4eI1lmtpXV1JXKM9JTNQccghqkvS11/p6f3tBEIR3n2WFrcpqeGiNVVxKaxKZqHZSjlEVVE+b3T3zvl1jDE8+MpJuoYC3HHpajxuF7/Y1cqSsgLOWa4uJSU1lPu8/Out57BxqYpDJqLikENM5Vba09xHQ3khFUXeKV+jeyjA//zRLp6y2yY/sa+Ne27axPMHO/jghSvVpaSklHdrckPGouKQQ1SXeBkKhBkJhON8uEc6Brnx319kZaWPBz9yITX2WNFYRkNhntnfzud+vpfuoQB/++4N1JcWcPdPdnHDv/2GiIHr1aWkKIsGFYccIoYUQDkAAAq5SURBVHZc6PLKsWHs//LUITxuobXPzx9+azsP3LGNiiIv4YjhhUMdPPJmC0/ubWNgNMTqmiK+fdv5UT/w2cvLuOuhN+kfCXLO8vIF+b0URZl/VBxyiOoSy2UUKw5vnern57ta+B/vPI2L11Tzx999hdv+cwdXb6jjhztO0tw7QmlBHtduquf6s5Zw0ZpqPO6xPIVlFT4e+siFRCLp7/aqKErmoOKQQ4xVSY/VOnz5iYMUe/P4yKWrKfd5+cYfnMtHvvcqu5r6eMdpVdx93elcvaEeb97UiWsaa1CUxUXGiIOIXAt8FXAD3zLG/PMCbynrWFbhw5vn4u8e3oMxhrrSAp7Y18ZdV62j3GdZFVecUcfDH7sInzePVdVFC7xjRVEyFTHGLPQeEBE3cBC4CmgCXgFuNcbsm+ycLVu2mJ07d87TDrOHnce6+euf7uFA2wClBXm4XcLz//NySrS7paIogIi8aozZMt1xmVIEtxU4bIw5YowJAA8ANy7wnrKSLY2V/OITF/NX151OOGK466p1KgyKosyYTHErNQAnY75vAi4Yf5CI3AHcAbBixYr52VkW4nG7uOPS0/jwxas1VqAoyqzIFMsh0SfYBH+XMeZeY8wWY8yWmpqaedhWdqPCoCjKbMkUcWgClsd8vwxoWaC9KIqiLHoyRRxeAdaKyCoR8QK3AI8s8J4URVEWLRkRczDGhETkY8CvsFJZv2OM2bvA21IURVm0ZIQ4ABhjHgUeXeh9KIqiKJnjVlIURVEyCBUHRVEUZQIqDoqiKMoEMqJ9xmwQkRFgJkHrFcCJOfzIMqAvjcfr/uZ2fCwz3etsftZ8/v10f/Fk4v9vLJm4v2qg03680hgzfaGYMSYrv4COdB6f4Px703y87m8Ox8/x2pjxz5rPv5/uL/P/fzN9f8DOmf6MbHYrzXRY8lyHK/88zcfr/uZ2fCwz3etsftZ8/v10f/Fk4v9vLJm+v6TIZrfSTpNEZ8HZHj/f6P5SRybuNRP3FEum7y+WTN9rJu5vNnvKZsvh3jQfP9/o/lJHJu41E/cUS6bvL5ZM32sm7m/Ge8pay0FRFEVJH9lsOSiKoihpQsVBURRFmUDOiYOIvEdEjIicvtB7mQoRGZzm+edEZF6DWiKyTEQeFpFDIvK2iHzV7pI72fGfEhHffO5x3M+f8m+4EGTD9ZeJ116CPWT8tZiJ118qyTlxAG4FfoPV9jtp7DnWixYREeAnwM+MMWuBdUAxcM8Up30KWDBxyFD0+psjei1mBjklDiJSDFwE3I795hSRy0TkeRH5qYjsE5H/EBGX/dygiPy9iGwHLlyA/V4mIr+I+f7fROSP5nsfNu8C/MaY/wQwxoSBTwN/IiJFIvJFEdktIrtE5OMi8glgKfCsiDy7QHtGRIpF5GkRec3e3432eqOI7BeRb4rIXhF5QkQK070XsuT6y7BrbzxZcy1m0vWXanJKHICbgMeNMQeBbhE5117fCvw5cCZwGvBee70I2GOMucAY85t5321msRF4NXbBGNOP1Qbgw8Aq4BxjzFnA940xX8Oa1ne5Meby+d5sDH7gPcaYc4HLgS/Zd54Aa4F/N8ZsxCpMel+a96LXX2rIpmsxk66/lJJr4nAr8ID9+AH7e4Adxpgj9h3ID4GL7fUw8OP53WLGIiSY222vXwr8hzEmBGCM6Z7PjU2DAP8oIruAp4AGoM5+7qgx5g378atAY5r3otdfasimazGTrr+UkjHDfuaKiFRhmaObRMRgTZQzWAOExl9ozvd++w27UISIF+iChdoIVhPDuDsbESnFmu19hMRv1kzgD4Aa4DxjTFBEjjH2dxyNOS4MpM2sz8LrL5OuvfFk07WYEddfOsgly+H9wP3GmJXGmEZjzHLgKNZd2lax5lO7gN/DChhmAseBDSKSLyJlwBULuJenAZ+IfAiiAdIvAd8FngD+h4jk2c9V2ucMACXzv9U4yoB2+415ObBygfaRbddfJl1748mmazFTrr+Uk0vicCvw03FrPwZ+H3gJ+GdgD9Ybdvxx84p9YY8aY04CDwG7gO8Dry/UnoxVKv8e4GYROQQcxPKn/hXwLSx/7y4ReRPrbwpWSf5jCxGQdv6GWH+3LSKyE+su7q353otNVlx/mXjtjScbrsUMvP5STs63zxCRy4C/MMa8e6H34iAiZwPfNMZsXei9ZCvZ8jfMtOsvW/5umc5i+DvmkuWQFYjI/8AKSv7NQu8lW9G/4ezQv1tqWCx/x5y3HBRFUZSZo5aDoiiKMgEVByXjEZHlIvKsXXG6V0Q+aa9XisiTYvXfeVJEKmLOuVtEDovIARG5xl4rEZE3Yr46ReRfFur3UrKDVF1/9vqtMdXdj4tI9UL8TsmgbiUl4xGRJcASY8xrIlKCVVB0E/BHQLcx5p9F5LNAhTHmL0VkA5ZPeCtWW4WngHXjawpE5FXg08aY5+fx11GyjFRdf1gFcy3ABmNMp4j8H2DYGPO5ef+lkkAtByXjMca0GmNesx8PAPuxKlFvBO6zD7sP6w2Lvf6AMWbUGHMUOIz1Ro0iImuBWuCF9P8GSjaTwutP7K8iu8VGKZZYZCQqDkpWISKNwDnAdqDOGNMK1hsY68MerDfuyZjTmuy1WG4FHjRqOiszYC7XnzEmCHwU2I1tQQDfnpeNzwIVByVrEKvr6Y+BT9mN2CY9NMHaeBG4Bcv0V5SkmOv1JyIeLHE4B8vdtAu4O+UbTREqDkpWYL+xfozVhfMn9nKb7Q92/MLt9noTVh8eh2XEmO92AVOeMSau86eiTEaKrr/NAMaYt22L9SHgHfOw/Vmh4qBkPLZ/9tvAfmPMl2OeegS4zX58G/BwzPotdt+gVVitk3fEnHcrajUoSZLC668Zq59VjX3cVVjxi4xEs5WUjEdELsYKHO8GIvbyX2H5fR8CVmD127nZaeEsIn8N/AlW99FPGWMei3m9I8B1xpic6YOjpI9UXn92dfUngSBW88M/MsZ0zd9vkzwqDoqiKMoE1K2kKIqiTEDFQVEURZmAioOiKIoyARUHRVEUZQIqDoqiKMoEVBz+f3t3DxpFFIVh+P0wNkHQxlbsRAyiWImItb2NiIZYCYKtpZWdimgTLSL+NGInFoqNgSj2GixEG4s0Qow/aY/FzELI7CYmJm6Q96mGw+VwbzNn7t2dM9ImSHKh9w3kPxy/N8n7zZyTtBYjw56A9L9JMlJVk8Oeh/Q3LA5SH22Dtec0LzodpvnI/TlgP3AD2AF8pXmJaS7JK+ANcAx42rZ2/llV15IcAiaBUeATcL6q5pMcAaaARWDm361OWp3HStJg+4C7VXUQ+A5cBG4Dp6qqd2O/umT8rqo6UVXXl+V5AFxu87wDrrTxe8Clqjq6mYuQ1sOdgzTYl6p63V4/ommZMAa8bNrtsA2YWzL+8fIESXbSFI3pNnQfeNIn/hA4ufFLkNbH4iANtry3zA9gdoUn/V9ryJ0++aUtw2MlabA9SXqF4DTwFtjdiyXZnuTASgmqagGYT3K8DZ0FpqvqG7DQNnUDOLPx05fWz52DNNgHYDzJHeAjze8NL4Bb7bHQCHATmF0lzzgwmWQU+AxMtPEJYCrJYptX2jLsyir10f5b6VlVjQ15KtJQeKwkSepw5yBJ6nDnIEnqsDhIkjosDpKkDouDJKnD4iBJ6vgNTDKVCoBvNT8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sorted_data['inc'][800:900].plot()"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"first_dec_week = [pd.Period(pd.Timestamp(y, 12, 3), 'W')\n",
" for y in range(1991,\n",
" sorted_data.index[-1].year)]"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [],
"source": [
"year = []\n",
"yearly_incidence = []\n",
"for week1, week2 in zip(first_dec_week[:-1],\n",
" first_dec_week[1:]):\n",
" one_year = sorted_data['inc'][week1:week2-1]\n",
" assert abs(len(one_year)-52) < 2\n",
" yearly_incidence.append(one_year.sum())\n",
" year.append(week2.year)\n",
"yearly_incidence = pd.Series(data=yearly_incidence, index=year)"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62c3b9e48>"
]
},
"execution_count": 22,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAD8CAYAAACyyUlaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHJBJREFUeJzt3X+QVeWd5/H3B1tBHTGAYATUdiJhBskMhi4kk51sXEdQsyVYYyYdHWUTqzD+mE2yfyRQcUtLmdk4ZeIMY8XV+As1ooaNJbtKFHWscXYRaNQoSBjakSCK0tkmgpmF2PDdP87Tcrptuu+9p+H+6M+r6tY9/b3nOZzHc+1vn+d5zvMoIjAzMytiWLVPwMzM6p+TiZmZFeZkYmZmhTmZmJlZYU4mZmZWmJOJmZkV5mRiZmaFOZmYmVlhJSUTSd+WtEHSeklLJY2QdIOktyW9kl4X5PZfKKld0iZJs3Px6ZJeS58tlqQUHy7pkRRfLak5V2aepM3pNW/wqm5mZoNFAz0BL2kC8M/AlIj4f5IeBZ4EmoEPIuKWXvtPAZYCM4DxwDPApyNin6Q1wDeBF9MxFkfECklXA38UEd+Q1ApcFBFfkTQaaANagADWAdMjYufBzveEE06I5ubmcv87mJkNaevWrft1RIyttHxTGfsdLelD4BjgHbJk0pc5wMMRsRd4U1I7MEPSFmBkRKwCkHQ/MBdYkcrckMovA25Ldy2zgZUR0ZnKrATOI0tWfWpubqatra3EapmZGYCkXxUpP2AzV0S8DdwCbAW2A+9HxNPp42slvSrpHkmjUmwC8FbuENtSbELa7h3vUSYiuoD3gTH9HMvMzGrIgMkkJYk5wGlkzVbHSvpL4HbgU8A0siTzg+4ifRwm+olXWiZ/jvMltUlq6+jo6Kc2ZmZ2KJTSAf9nwJsR0RERHwI/A/4kIt6LiH0RsR/4MVkfCWR3Dyfnyk8kaxbblrZ7x3uUkdQEHA909nOsHiLizohoiYiWsWMrbvIzM7MKlZJMtgIzJR2T+jHOATZKOim3z0XA+rS9HGhNI7ROAyYBayJiO7Bb0sx0nMuBx3NlukdqXQw8F9nIgKeAWZJGpTukWSlmZmY1ZMAO+IhYLWkZ8BLQBbwM3AncJWkaWbPTFuDKtP+GNOLr9bT/NRGxLx3uKuA+4GiyjvcVKX438EDqrO8EWtOxOiXdBKxN+93Y3RlvZma1Y8ChwfWmpaUlPJrLrDw7du3h2qUvc9slZzLuuBHVPh2rAknrIqKl0vJ+At7MWPzsZtZu6WTxM5urfSpWp0p9zsTMGtDk61awt2v/Rz8/uHorD67eyvCmYWxadH4Vz8zqje9MzIawF75zNhdOG8+II7NfBSOOHMacaeN54btnV/nMrN44mZgNYeNGjuC44U3s7drP8KZh7O3az3HDm9xvYmVzM5fZEPfrD/Zy6VmncsmMU3hozVY6du+p9ilZHfJoLjMz82guMzOrPicTMzMrzMnEzMwKczLJ2bFrD39xxyp2uAPSzKwsTiY5fgrYzKwyHhqMnwI2MyvKdyb4KWAzs6KcTPBTwGZmRbmZK/FTwGZmlfMT8GZm5ifgzcys+pxMzMysMCcTMzMrzMnEzMwKczIxM7PCnEzMzKwwJxMzMyvMycSsBJ5R2qx/TiZmJfCM0mb9KymZSPq2pA2S1ktaKmmEpNGSVkranN5H5fZfKKld0iZJs3Px6ZJeS58tlqQUHy7pkRRfLak5V2Ze+jc2S5o3eFU3G9jk61bQvOAJHly9lYhsRunmBU8w+boV1T41s5oyYDKRNAH4z0BLREwFjgBagQXAsxExCXg2/YykKenzM4DzgB9JOiId7nZgPjApvc5L8SuAnRFxOnArcHM61mjgeuAsYAZwfT5pmR1qnlHarDSlNnM1AUdLagKOAd4B5gBL0udLgLlpew7wcETsjYg3gXZghqSTgJERsSqyCcHu71Wm+1jLgHPSXctsYGVEdEbETmAlBxKQ2SHnGaXNSjNgMomIt4FbgK3AduD9iHgaODEitqd9tgPjUpEJwFu5Q2xLsQlpu3e8R5mI6ALeB8b0cyyzw6Z7RunHrv48l551Kh0f7K32KVWdByRYbwNOQZ+aleYApwG/AX4q6S/7K9JHLPqJV1omf47zyZrPOOWUU/o5NbPy3XHZgYlUF82dWsUzqR35AQmLLvpMtU/HakAp65n8GfBmRHQASPoZ8CfAe5JOiojtqQlrR9p/G3ByrvxEsmaxbWm7dzxfZltqSjse6EzxL/Yq83zvE4yIO4E7IZuCvoQ6mVkFvMS1HUwpfSZbgZmSjkn9GOcAG4HlQPfoqnnA42l7OdCaRmidRtbRviY1he2WNDMd5/JeZbqPdTHwXOpXeQqYJWlUukOalWJmVgUekGAHM+CdSUSslrQMeAnoAl4muwv4PeBRSVeQJZwvp/03SHoUeD3tf01E7EuHuwq4DzgaWJFeAHcDD0hqJ7sjaU3H6pR0E7A27XdjRHQWqrGZVcwDEuxgvNKimZXlygfaGHvciB5LXOf7law+FV1p0cnEzMy8bK+ZmVWfk4mZmRXmZGJmZoU5mZiZWWFOJmZmVpiTiZmZFeZkYmZmhTmZmNURz9ZrtcrJxKyOePlgq1WlzBpsZlXm2Xqt1vnOxKwOeLZeq3VOJmZ1wLP1Wq1zM5dZnehePjg/W69ZrfCswWZm5lmDzaxvHkZsh5OTiVmD8jBiO5zcZ2LWYDyM2KrBdyZmDcbDiK0anEzMGoyHEVs1uJnLrAF5GLEdbh4abGZWwI5de7h26cvcdsmZdX3356HBZmZV5FFzGTdzmZlVwKPmevKdiZlZBTxqrqcBk4mkyZJeyb12SfqWpBskvZ2LX5Ars1BSu6RNkmbn4tMlvZY+WyxJKT5c0iMpvlpSc67MPEmb02ve4FbfzKwyHjXX04DNXBGxCZgGIOkI4G3gMeBrwK0RcUt+f0lTgFbgDGA88IykT0fEPuB2YD7wIvAkcB6wArgC2BkRp0tqBW4GviJpNHA90AIEsE7S8ojYWbjmZmYFedTcAeX2mZwDvBERv0o3FX2ZAzwcEXuBNyW1AzMkbQFGRsQqAEn3A3PJkskc4IZUfhlwW7prmQ2sjIjOVGYlWQJaWuZ5m5kNujsuOzD4adHcqVU8k+ort8+klZ6/yK+V9KqkeySNSrEJwFu5fbal2IS03Tveo0xEdAHvA2P6OZaZmdWQkpOJpKOAC4GfptDtwKfImsC2Az/o3rWP4tFPvNIy+XObL6lNUltHR8dB62BmZodGOXcm5wMvRcR7ABHxXkTsi4j9wI+BGWm/bcDJuXITgXdSfGIf8R5lJDUBxwOd/Ryrh4i4MyJaIqJl7NixZVTJzMwGQznJ5KvkmrgknZT77CJgfdpeDrSmEVqnAZOANRGxHdgtaWbqD7kceDxXpnuk1sXAc5E9mv8UMEvSqNSMNivFzMyshpTUAS/pGOBc4Mpc+G8lTSNrdtrS/VlEbJD0KPA60AVck0ZyAVwF3AccTdbxviLF7wYeSJ31nWR9M0REp6SbgLVpvxu7O+PNzKx2eG4uMzPz3Fx2aHjJVzMrh5OJ9cmT15lZOTzRo/XgyevMrBK+M7EePHmdmVXCycR68OR1ZlYJN3PZx3jyOjMrl4cGm5mZhwabmVn1OZmYmVlhTiZmZlaYk4mZmRXmZGJmZoU5mZiZWWFOJmZmVpiTiZmZFeZkYmZmhTmZmJlZYU4mNmR5ATCzweNkYkOWFwAzGzyeNdiGHC8AZjb4fGdiQ44XADMbfE4mNuR4ATCzwedmLhuSvACY2eDy4lhmZubFsczMrPoGTCaSJkt6JffaJelbkkZLWilpc3oflSuzUFK7pE2SZufi0yW9lj5bLEkpPlzSIym+WlJzrsy89G9sljRvcKtvZmaDYcBkEhGbImJaREwDpgP/BjwGLACejYhJwLPpZyRNAVqBM4DzgB9JOiId7nZgPjApvc5L8SuAnRFxOnArcHM61mjgeuAsYAZwfT5pmZlZbSi3mesc4I2I+BUwB1iS4kuAuWl7DvBwROyNiDeBdmCGpJOAkRGxKrKOmvt7lek+1jLgnHTXMhtYGRGdEbETWMmBBGRmZjWi3GTSCixN2ydGxHaA9D4uxScAb+XKbEuxCWm7d7xHmYjoAt4HxvRzLDMzqyElJxNJRwEXAj8daNc+YtFPvNIy+XObL6lNUltHR8cAp2eHgue5MhvayrkzOR94KSLeSz+/l5quSO87UnwbcHKu3ETgnRSf2Ee8RxlJTcDxQGc/x+ohIu6MiJaIaBk7dmwZVbLB4nmuzIa2ch5a/CoHmrgAlgPzgO+n98dz8Yck/RAYT9bRviYi9knaLWkmsBq4HPiHXsdaBVwMPBcRIekp4G9yne6zgIVl1tEOIc9zZWZQ4p2JpGOAc4Gf5cLfB86VtDl99n2AiNgAPAq8DvwcuCYi9qUyVwF3kXXKvwGsSPG7gTGS2oH/QhoZFhGdwE3A2vS6McWsRnieKzODEu9MIuLfyDrE87H/Sza6q6/9/xr46z7ibcDUPuJ7gC8f5Fj3APeUcp52+HmeKzMDz81lg8DzXJmZ5+YyMzPPzWVmZtXnZGJmZoU5mZiZWWFOJmZWMzyTQv1yMjGzmuGZFOqXhwabWdV5JoX65zsTM6s6z6RQ/5xMzKzqPJNC/XMzl5nVBM+kUN/8BLyZ1aUdu/Zw7dKXue2SM30HMwj8BLyZDUke+VVb3MxlZnXFI79qk+9MzKyueORXbXIyMbO64pFftcnNXGZWdzzyq/Z4NJeZmXk0l5mZVZ+TiZmZFeZkYmZmhTmZHAZeo8HMGp2TyWHgJ3XNrNF5aPAh5Cd1zWyo8J3JIVTJk7puEqt/voa1ydfl0CopmUj6hKRlkn4paaOkz0m6QdLbkl5Jrwty+y+U1C5pk6TZufh0Sa+lzxZLUooPl/RIiq+W1JwrM0/S5vSaN3hVP/QqeVLXTWL1z9ewNvm6HFolPbQoaQnwQkTcJeko4BjgW8AHEXFLr32nAEuBGcB44Bng0xGxT9Ia4JvAi8CTwOKIWCHpauCPIuIbklqBiyLiK5JGA21ACxDAOmB6ROw82LnW2kOLVz7QxtjjRvR4UveOyz7+XFDvJrFubhKrH76GtcnXpTSH/KFFSSOBLwB3A0TE7yLiN/0UmQM8HBF7I+JNoB2YIekkYGRErIosg90PzM2VWZK2lwHnpLuW2cDKiOhMCWQlcF7ZtayiOy5rYdHcqUwZP5JFc6f2mUjAk9c1Al/D2uTrcniU0sz1+0AHcK+klyXdJenY9Nm1kl6VdI+kUSk2AXgrV35bik1I273jPcpERBfwPjCmn2M1HE9eV/98DWuTr8vhUUoyaQI+C9weEWcCvwUWALcDnwKmAduBH6T91ccxop94pWU+Imm+pDZJbR0dHf1UpbZ1T1732NWf59KzTqXjg73VPiUrk69hbfJ1OfQG7DOR9EngxYhoTj//KbAgIr6U26cZ+F8RMVXSQoCI+G/ps6eAG4AtwD9GxB+k+FeBL0bEld37RMQqSU3Au8BYoLV7n1TmDuD5iFh6sPOttT4TM7N6cMj7TCLiXeAtSZNT6Bzg9dQH0u0iYH3aXg60phFapwGTgDURsR3YLWlm6g+5HHg8V6Z7pNbFwHOpX+UpYJakUakZbVaKmZlZDSn1ocW/An6SRnL9K/A1YLGkaWTNTluAKwEiYoOkR4HXgS7gmojYl45zFXAfcDSwIr0g69x/QFI70El2R0JEdEq6CVib9rsxIjorq6qZmR0qXs/EzMy8nomZmVWfk4mZmRXmZGJmZoU5mZiZWWFOJmZ2SHm23qHBycTMDinP1js0eHEsO+x27NrDtUtf5rZLzvT8SA3Mi8MNLb4zqVA93rrXyjn7L9WhwbP1Di1OJhWqx1+I1T7nydetoHnBEzy4eisR2V+qzQueYPJ1KwYubHWnnmfrrZU/vOqJn4AvU60ttFNKk1GtnPOOXXtY9ORGnt7wLns+3M+II4cx+4xP8r0v/WFd/IKx8pW6OFytue6x1/jJmq1cOuMUFl30mWqfzmFR9Al495mU6YXvnH3QX4jVkL/bONiXvlbOuZ7/UrXK5BPHorlTq3gmpXE/T+WcTMpUK78Qy/nS18o5w4F1JfJ/qZrVilr5w6seOZlUoBZ+IZb7pa+Fc4b6+0vVhpZa+sOr3jiZVKAWfiGW+6WvhXM2qwe18odXvXEyqWP+0psNPv/hVRmP5jIzM69nYmZm1edkYmZmhTmZmJlZYU4mZmZWmJOJmZkV5mRiZmaFOZmYmVlhTiZmZlaYk4mZmRVWUjKR9AlJyyT9UtJGSZ+TNFrSSkmb0/uo3P4LJbVL2iRpdi4+XdJr6bPFkpTiwyU9kuKrJTXnysxL/8ZmSfMGr+pWD7xIkVl9KPXO5O+Bn0fEHwB/DGwEFgDPRsQk4Nn0M5KmAK3AGcB5wI8kHZGOczswH5iUXuel+BXAzog4HbgVuDkdazRwPXAWMAO4Pp+0rPFVe3VIMyvNgBM9ShoJfAH4TwAR8Tvgd5LmAF9Muy0Bnge+C8wBHo6IvcCbktqBGZK2ACMjYlU67v3AXGBFKnNDOtYy4LZ01zIbWBkRnanMSrIEtLRAna0OeJEis/pSyp3J7wMdwL2SXpZ0l6RjgRMjYjtAeh+X9p8AvJUrvy3FJqTt3vEeZSKiC3gfGNPPsazBvfCds7lw2nhGHJl9RUccOYw508bzwnfPrvKZmVlfSkkmTcBngdsj4kzgt6QmrYNQH7HoJ15pmQP/oDRfUpukto6Ojn5OzeqFFykyqy+lJJNtwLaIWJ1+XkaWXN6TdBJAet+R2//kXPmJwDspPrGPeI8ykpqA44HOfo7VQ0TcGREtEdEyduzYEqpk9aB7vZbHrv48l551Kh0f7K32KZnZQQzYZxIR70p6S9LkiNgEnAO8nl7zgO+n98dTkeXAQ5J+CIwn62hfExH7JO2WNBNYDVwO/EOuzDxgFXAx8FxEhKSngL/JdbrPAhYWrrXVBS9SZFY/Sl1p8a+An0g6CvhX4GtkdzWPSroC2Ap8GSAiNkh6lCzZdAHXRMS+dJyrgPuAo8k63lek+N3AA6mzvpNsNBgR0SnpJmBt2u/G7s54MzOrHV5p0czMvNKimZlVn5OJmZkV5mRiZmaFOZmYmVlhTiZmZlaYk4mZmRXmZGJmZoU5mZiZWWFOJmZmVpiTiZmZFeZkYmZmhTmZmJk1gB279vAXd6xix+49Vfn3nUzMqqjavwCscSx+djNrt3Sy+JnNVfn3S52C3swOgfwvgEUXfabap2N1aPJ1K9jbtf+jnx9cvZUHV29leNMwNi06/7Cdh6egN6uC3r8Auh3uXwBW/3bs2sOiJzfy9IZ32fPhfkYcOYzZZ3yS733pD8ta5tpT0JvVoRe+czYXThvPiCOz/wVHHDmMOdPG88J3z67ymVm9GTdyBMcNb2Jv136GNw1jb9d+jhveVFYiGQxu5jKrglr5BWCN4dcf7OXSs07lkhmn8NCarXRUoQ/OycSsSmrhF4A1hjsuO9A6tWju1Kqcg/tMzMzMfSZmZlZ9TiZmZlaYk4mZmRXmZGJmZoU5mZiZWWFOJmZmNaqe5m4rKZlI2iLpNUmvSGpLsRskvZ1ir0i6ILf/QkntkjZJmp2LT0/HaZe0WJJSfLikR1J8taTmXJl5kjan17zBqriZWa2r9uSN5SjnocWzI+LXvWK3RsQt+YCkKUArcAYwHnhG0qcjYh9wOzAfeBF4EjgPWAFcAeyMiNMltQI3A1+RNBq4HmgBAlgnaXlE7Cy3omZm9aJWJm8sx6Fo5poDPBwReyPiTaAdmCHpJGBkRKyK7EnJ+4G5uTJL0vYy4Jx01zIbWBkRnSmBrCRLQGYfU09NAjY0lfodrce520pNJgE8LWmdpPm5+LWSXpV0j6RRKTYBeCu3z7YUm5C2e8d7lImILuB9YEw/xzL7mHpqErChqdTvaD3O3VZqM9fnI+IdSeOAlZJ+SdZkdRNZorkJ+AHwdUB9lI9+4lRY5iMpwc0HOOWUU/qviTWcemwSsKGlku9ovc3dVtKdSUS8k953AI8BMyLivYjYFxH7gR8DM9Lu24CTc8UnAu+k+MQ+4j3KSGoCjgc6+zlW7/O7MyJaIqJl7NixpVTJGkg9NgnY0FLJd/SOy1pYNHcqU8aPZNHcqT0mc6xFAyYTScdKOq57G5gFrE99IN0uAtan7eVAaxqhdRowCVgTEduB3ZJmpv6Qy4HHc2W6R2pdDDyX+lWeAmZJGpWa0WalmNlH6rFJwIaWofAdLaWZ60TgsTSKtwl4KCJ+LukBSdPImp22AFcCRMQGSY8CrwNdwDVpJBfAVcB9wNFko7hWpPjdwAOS2snuSFrTsTol3QSsTfvdGBGdlVfXGlW9NQnY0NPo31FPQW9mZp6C3szMqs/JxMzMCnMyMTOzwpxMzMysMCcTMzMrzMnEzMwKa7ihwZI6gF9V+zwKOgHoPUNzo2n0OjZ6/aDx6zjU6ndqRFQ8hUjDJZNGIKmtyHjvetDodWz0+kHj19H1K4+buczMrDAnEzMzK8zJpDbdWe0TOAwavY6NXj9o/Dq6fmVwn4mZmRXmOxMzMyvMyeQwSUsb75C0Phf7Y0mrJL0m6X9KGpniR0m6N8V/IemLuTLPS9ok6ZX0GleF6nyMpJMl/aOkjZI2SPpmio+WtFLS5vQ+KldmoaT2VJ/Zufj0VPd2SYvT+jdVNcj1a4hrKGlM2v8DSbf1OlbdX8MB6ldz17CC+p2rbCn219L7f8gdq/zrFxF+HYYX8AXgs8D6XGwt8O/T9teBm9L2NcC9aXscsA4Yln5+Hmipdn36qN9JwGfT9nHAvwBTgL8FFqT4AuDmtD0F+AUwHDgNeAM4In22Bvgc2bLNK4DzG6x+jXINjwX+HfAN4LZex2qEa9hf/WruGlZQvzOB8Wl7KvB2kevnO5PDJCL+iWzhr7zJwD+l7ZXAn6ftKcCzqdwO4DdATY93j4jtEfFS2t4NbAQmAHOAJWm3JcDctD0HeDgi9kbEm0A7MEPZCp4jI2JVZN/q+3Nlqmaw6nd4z7o85dYxIn4bEf8M9FjlqVGu4cHqV6sqqN/LkZZkBzYAI5StkFvR9XMyqa71wIVp+8scWO/+F8AcSU3Klj6envsM4N50a/1fa6H5oDdJzWR/9awGToxsyWbSe3dzwATgrVyxbSk2IW33jteMgvXr1gjX8GAa5RoOpGavYQX1+3Pg5YjYS4XXz8mkur4OXCNpHdlt6e9S/B6yC9gG/B3wf8iWQAa4NCI+A/xpel12WM94AJJ+D/gfwLciYld/u/YRi37iNWEQ6geNcw0Peog+YvV4DftTs9ew3PpJOgO4mbT0OhVePyeTKoqIX0bErIiYDiwla1cnIroi4tsRMS0i5gCfADanz95O77uBh6ihphNJR5J9iX8SET9L4ffSbXN388eOFN9Gz7uticA7KT6xj3jVDVL9GukaHkyjXMODqtVrWG79JE0EHgMuj4g3Urii6+dkUkXdI0AkDQOuA/57+vkYScem7XOBroh4PTV7nZDiRwL/kayprOrSbf7dwMaI+GHuo+XAvLQ9D3g8F29NbbSnAZOANek2fLekmemYl+fKVM1g1a/BrmGfGugaHuw4NXkNy62fpE8ATwALI+J/d+9c8fWr1siDofYiu/PYDnxIlvmvAL5JNuLiX4Dvc+Ah0mZgE1kH2jNks3lCNrpkHfAqWYfZ35NGCFX7RTbqJdK5vZJeFwBjyAYTbE7vo3Nlvkd2N7aJ3GgRssEG69Nnt3X/d2mE+jXgNdxCNrDkg/S9ntJg1/Bj9avVa1hu/cj+gP1tbt9XgHGVXj8/AW9mZoW5mcvMzApzMjEzs8KcTMzMrDAnEzMzK8zJxMzMCnMyMTOzwpxMzMysMCcTMzMr7P8DIQIRMMTgjDAAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"yearly_incidence.plot(style='*')"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"2002 548385\n",
"2017 548682\n",
"2018 553291\n",
"2006 570156\n",
"2019 571344\n",
"2001 601235\n",
"1996 625186\n",
"2015 625728\n",
"1997 627588\n",
"1993 627750\n",
"2000 629803\n",
"2011 631660\n",
"2014 652534\n",
"1995 652845\n",
"2005 652865\n",
"1994 656733\n",
"2013 672421\n",
"2012 680066\n",
"1999 695096\n",
"2003 712485\n",
"1998 737255\n",
"2016 757015\n",
"2008 765806\n",
"2007 775543\n",
"2010 797113\n",
"2004 803801\n",
"1992 824689\n",
"2009 838983\n",
"dtype: int64"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"yearly_incidence.sort_values()"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fc62c4cd5f8>"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEp1JREFUeJzt3X+QJHV5x/H3w3HowSISDlY9NGfEENFTkVWjqTK7hLJQSCwtYySIMdGcVQhqQkrxR8UYIV5UoqlEU16pgfxyIREq/DAmGF1RY9Q7STyQoJScCBgIgsgigidP/uhe3bsb2Lm9mdnpp96vqq3b6e3r+X72u/uZnp7u2chMJEk17bPSA5AkDY8lL0mFWfKSVJglL0mFWfKSVJglL0mFWfKSVJglL0mFWfKSVNi+w9jo2rVrc/369QPb3t13380BBxwwsO2tpEpZwDzjrFIWqJWnV5atW7felpmHDvq+hlLy69evZ8uWLQPb3tzcHNPT0wPb3kqqlAXMM84qZYFaeXpliYhvDeO+PFwjSYVZ8pJUmCUvSYVZ8pJUmCUvSYX1dXZNRGwH7gJ+DOzIzKlhDkqSNBh7cgrlTGbeNrSRSJIGzsM1klRY9PM3XiPieuAOIIEPZubmHutsBDYCTE5OHjM7O7usAW276c7dlk2ugVvuWdbm+rZh3UHDvYPW/Pw8ExMTOy3rlbkrljs3o/p+76le89NVlbJArTy9sszMzGwdxqHwfkv+UZl5c0QcBlwOnJ6ZVzzQ+lNTU7ncK17Xn3nZbsvO2LCDc7YN5eLcn9i+6YShbn9BryvdemXuiuXOzai+33uq+lWVXVYpzwNc8TqUku/rcE1m3tz+eytwEfCMQQ9EkjR4S5Z8RBwQEQcufA48F7hq2AOTJO29fp5nTwIXRcTC+v+QmZ8Y6qgkSQOxZMln5jeBp4xgLJKkAfMUSkkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpMIseUkqzJKXpML6LvmIWBURV0bEpcMckCRpcPZkT/51wDXDGogkafD6KvmIOBw4AfjQcIcjSRqkfvfk3we8Abh/iGORJA1YZOaDrxBxIvD8zDw1IqaBP8jME3ustxHYCDA5OXnM7Ozssga07aY7d1s2uQZuuWdZm+vbhnUHDfcOWvPz80xMTOy0rFfmrlju3Izq+72nes1PV1XKArXy9MoyMzOzNTOnBn1f/ZT8O4FTgB3AQ4GHARdm5sse6P9MTU3lli1bljWg9WdettuyMzbs4Jxt+y5re/3avumEoW5/wdzcHNPT0zst65W5K5Y7N6P6fu+pXvPTVZWyQK08vbJExFBKfsnDNZn5psw8PDPXAy8FPvVgBS9JGh+eJy9Jhe3R8+zMnAPmhjISSdLAuScvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYVZ8pJUmCUvSYUtWfIR8dCI+FJE/HdEXB0Rbx/FwCRJe2/fPta5Fzg2M+cjYjXwuYj4l8z8zyGPTZK0l5Ys+cxMYL69ubr9yGEOSpI0GNF0+BIrRawCtgJHAO/PzDf2WGcjsBFgcnLymNnZ2WUNaNtNd+62bHIN3HLPsjbXtw3rDhruHbTm5+eZmJjYaVmvzF2x3LkZ1fd7T/Wan65ayLJSP1+DnuOKc7PYzMzM1sycGvR99VXyP1k54uHARcDpmXnVA603NTWVW7ZsWdaA1p952W7Lztiwg3O29XNkafm2bzphqNtfMDc3x/T09E7LemXuiuXOzai+33uq1/x01UKWlfr5GvQcV5ybxSJiKCW/R2fXZOb3gDng+EEPRJI0eP2cXXNouwdPRKwBjgP+Z9gDkyTtvX6eZz8SOK89Lr8PcEFmXjrcYUmSBqGfs2u+Chw9grFIkgbMK14lqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKs+QlqTBLXpIKW7LkI+LREfHpiLgmIq6OiNeNYmCSpL23bx/r7ADOyMyvRMSBwNaIuDwzvzbksUmS9tKSe/KZ+Z3M/Er7+V3ANcC6YQ9MkrT3IjP7XzliPXAF8KTM/P4uX9sIbASYnJw8ZnZ2dlkD2nbTnbstm1wDt9yzrM31bcO6g4Z7B635+XkmJiZ2WtYrc1eMYm4Gaal57jU/g7ASc9y1uVlKF/L02yO9fs5mZma2ZubUoMfUd8lHxATwGeDszLzwwdadmprKLVu2LGtA68+8bLdlZ2zYwTnb+jmytHzbN50w1O0vmJubY3p6eqdlvTJ3xSjmZpCWmude8zMIKzHHXZubpXQhT7890uvnLCKGUvJ9nV0TEauBjwF/v1TBS5LGRz9n1wTwYeCazPyz4Q9JkjQo/ezJ/xJwCnBsRPxX+/H8IY9LkjQASx7gyszPATGCsUiSBswrXiWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpMEtekgqz5CWpsCVLPiI+EhG3RsRVoxiQJGlw+tmTPxc4fsjjkCQNwZIln5lXALePYCySpAGLzFx6pYj1wKWZ+aQHWWcjsBFgcnLymNnZ2WUNaNtNd+62bHIN3HLPsjY3diplAfOMs0pZoBt5Nqw7qK/15ufnmZiY2GnZzMzM1sycGvSY9h3UhjJzM7AZYGpqKqenp5e1nVecedluy87YsINztg1sqCuqUhYwzzirlAW6kWf7ydN9rTc3N8dyO3JPeXaNJBVmyUtSYf2cQvlR4AvAkRFxY0S8cvjDkiQNwpIHuDLzpFEMRJI0eB6ukaTCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKsySl6TCLHlJKqyvko+I4yPi2oi4LiLOHPagJEmDsWTJR8Qq4P3A84CjgJMi4qhhD0yStPf62ZN/BnBdZn4zM+8DZoEXDHdYkqRBiMx88BUiXgwcn5mvam+fAjwzM0/bZb2NwMb25pHAtQMc51rgtgFubyVVygLmGWeVskCtPL2y/GxmHjroO9q3j3Wix7LdHhkyczOwea9H1GsAEVsyc2oY2x61SlnAPOOsUhaolWeUWfo5XHMj8OhFtw8Hbh7OcCRJg9RPyX8ZeHxEPDYi9gNeClw83GFJkgZhycM1mbkjIk4D/hVYBXwkM68e+sh2NpTDQCukUhYwzzirlAVq5RlZliVfeJUkdZdXvEpSYZa8JBVmyUtSYWVLPiKOjYjHrvQ4BqFSFqiVp1IWqJWnUhZYfp5yL7y276szC3wPuB94W2Z+ZmVHtTyVskCtPJWyQK08lbLA3ufp/J58RBweEQ9btOg3gI9l5nNovjEnRcSzVmZ0e6ZSFqiVp1IWqJWnUhYYfJ7OlnxEPCEiPg58DvjjiFh407QfAvu3n19A8/4QvxgRvd6eYSxUygK18lTKArXyVMoCw8vTqZKPiAMW3XwqcGNmrgc+BbynXX47cG9EHJiZtwNfByaB9SMc6pIqZYFaeSplgVp5KmWB0eQZ+5KPiIMj4tyI+DKwKSIObR/Bngx8PiIiMy8GvhcRJ9B8Aw4ENrSb+AbNO77dtxLjX6xSFqiVp1IWqJWnUhYYfZ6xL3ngOcAO4Pk074j5ZuBhNGN/RP70lePzgN8EvgTcRfNHTsjMLwDHAt8f7bB7qpQFauWplAVq5amUBUadJzPH4oPmfXFeDXyG5n3p17bLLwBe237+WGBT+/Wn0xy7WtV+bQL4v3Y764BrgNOAvwY+AOxvFvNUylItT6Us45RnxffkF7148KvArwFvB54FvKtdfjnw7PbzbwOfBZ6XmV+meRScAcjMeeCLwNMz8ybgFJrjVv8LvDUzfzD8ND9xIgWyODfjm8W5Gd8s4zY3/fzRkIGLiCngVTRPQd4N3Ar8HPDNzPxURFwPvDsingtsBV4YEWsz87aI+AZwd0Q8BvgL4GURcRjN+9x/l+apDZm5BdgygixPp3lK9Vngssy8F/j5LmZp8zg3Y5ilzdP5uWmPN6dzM7o8I9uTj4iJiHhoRJwHfAy4HvjzzLw1IvaheQT7ekSsyczraZ6aPJnmm3YzzbmiAD+mefqyT7udzcDJwDHA5sy8fwRZDmn/fTZwLs0pTscB72xXub8rWdocB0dzbu5m4EK6PTeHRcQh7S/deXR/btZGxNMi4i+Bi+j23OzXFvwv0xxy6PrcrI6INRFxLuPcaUM+JrU/8Fs0pwOd3y77U+DVi9bZt/33NOAsYH17+8Q27Nr2823AQTQvOHwc2G/RNvYZco6FK4N/neaY2RzN37F9PfCa9msHA18Fjm4nb9M4ZumR5/PAJ2l++Do3N+19HAC8guZp8J3ACV2dm13yfJJmT+5FHZ6b1cCpwCXAXwFHAK/r8Nws5LmUZq/7ie2Yx3ZuhrYnHxGrgauBFwPvysyFR61twFER8c72EfB3IuJngE8AjwCe0K53BfDLwH2ZeSnwYeCfgPfT7KH9aOG+csiP3JmZEXEQ8BLgfZk5nZnX0jzSLqxzB/DPwGtpyuawcczSI897M/O4zPwqHZybiDiC5in/rwBvBW4CbqB5QWtVO4bOzM0ued5Cc/z1W+04N3RpblqvoSmx99H8begXtbd/3I6hM3PTWsjzXpqfrxfS/LwdGRGbxnJuhvyodyFw8i7LHkNzae75wEk0k/+J9muvab92MLCG5tH/MYv+79pRPFo/QJZTgXe0ny/sCb8I+OKidR4F3DzuWXbN095eAzyyHfMFXZkbml+0hyy6/RGaHYsXdHFueuT5EM1e/WEd/b25BHh5+/krgdNp9ti/1LW56ZHnt4E3tZ12/rjOzbC/ISfSnMh/DvBp4I9oXkxYt2id1TSHP45tb58N/BvN09Q3ruSE7pLlJe04Twa+QvPI+wKawwNrF613OfDM9vOzxjFLjzxbgb+jORtg8Q9gJ+Zm0XgPoTnWexzN8c3bac477tTc9Mjz3PZ2F39vTgX+HfhHmjNJPkRziOMO4NCuzc0ueW6g2al4AvDocZ2boZ5dk5mXRsTpbbiX0xy7ej3wwUWrHQJsp3k/BmgeCI4ErsvMHw5zfHvoSuAhwDNpnko/j/apF83Ts/cAD6d5av2d9v+8nfHMAjvnOa79OIXmIo0b2nW6MjcAZOZ3I+KpwCWZeX807wPye8Ab26fPXZkbYKc8F7WLvrPoy52Ym8z8QETcR3PK4NE0x6lfTHP2ye9GxCY69HuzS56n0ZwH//s0h16+3a42VnMzirNrXpiZf5KZ36Z5Ff0RwMMj4ikR8RaaQzo/yOaYMJn5o8y8atwml6b47qN5UeUOmhdebmj/3R+4mOYRflVm3gBjnQV2z3MZ8DXg6Ij4hY7NDe3ZDNC8MH5E+/nZzZfiEro1Nw+UZ5/2TJs306G5AY4C5jLzNppnJquAj9Icvuja7w3snOc84B7gcRHxxHH8vRn6efK58wn7d9GczH8NzYtKq2leZb9y2OPYW5l5b7vXcVa7aB+aCxw2ZeaVEfEy4OouZIGeefajedr5DpozovajI3MDzQtVEbHwTn3XtcuuAd4QEScDX+tKFtgtz/Z22Y6IeCnNz14n5iYiVgE30lzK/zfA3TRXdp6amdd37femR54fAo+nudDpVJpnx2M1N0P/oyER8RDgeJpDAU+kOe/zvUO90yGKiLNpLnI4muY429vaPeFO2iXPJ4E35GivchyoiLgW+MPMPH/hwpuVHtPeWJxnpceyXBHxOJpTB++j+Vm7CDgrmys6O6dHnouBt2TmWLwB2q5G8pehIuLVNBc6/O0YPwXrW0QcCXyrQhaokWfRlZRPpTl1d0eXC75gnsfQ7PH+R2bes9Lj2VtdylPuz/9Jkn5qxd+gTJI0PJa8JBVmyUtSYZa8JBVmyUtSYZa8JBVmyUtSYf8PAj4ci+HvWWMAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"yearly_incidence.hist(xrot=20)"
]
},
{
"cell_type": "code",
"execution_count": null,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment