{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import isoweek" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "data_file = \"syndrome-varicelle.csv\"\n", "\n", "import os\n", "import urllib.request\n", "if not os.path.exists(data_file):\n", " urllib.request.urlretrieve(data_url, data_file)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "raw_data = pd.read_csv(data_file, skiprows=1)\n", "raw_data[raw_data.isnull().any(axis=1)]\n", "data = raw_data.dropna().copy()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "def convert_week(year_and_week_int):\n", " year_and_week_str = str(year_and_week_int)\n", " year = int(year_and_week_str[:4])\n", " week = int(year_and_week_str[4:])\n", " w = isoweek.Week(year, week)\n", " return pd.Period(w.day(0), 'W')\n", "\n", "data['period'] = [convert_week(yw) for yw in data['week']]" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXm8HUWZPv68fc69N/u+khASIOxLIDFsKktAAuiAChp1BL+DRhF/7iKIjjqYUXSUkRFQBEdAFBhQQVaRgCwJhJsA2QgkISvZ9+Qmdz31+6O7uqura+u+59xzl34+n3vPOd3VVdXdVfXWuxNjDDly5MiRI4cIr9odyJEjR44cnQ85cciRI0eOHAnkxCFHjhw5ciSQE4ccOXLkyJFAThxy5MiRI0cCOXHIkSNHjhwJ5MQhR44cOXIkkBOHHDly5MiRQE4ccuTIkSNHAsVqdyArhg0bxsaPH1/tbuTIkSNHl8L8+fO3McaG28p1WeIwfvx41NfXV7sbOXLkyNGlQERrXMrlYqUcOXLkyJFAThxy5MiRI0cCOXHIkSNHjhwJ5MQhR44cOXIkkBOHHDly5MiRQE4ccuTIkSNHAjlxyJEjR44cCeTEIUe3xOOLNmJHQ3O1u5EjR5dFThxydDts3duEL967ADPvzp0kc+TIipw45Oh2aGkrAQDe3XWgyj3JkaPrIicOObotGKt2D3Lk6LrIiUOObgeiavcgR46uDytxIKJeRDSPiN4goiVE9MPg+A+I6F0iej34u1C45joiWkFEbxHR+cLxyUS0KDh3M5E/jYmojojuD46/QkTjy3+rOXLkyJHDFS6cQxOAcxhjJwKYBGA6EZ0anLuJMTYp+HscAIjoGAAzABwLYDqAW4moEJS/DcBMABODv+nB8SsB7GSMHQ7gJgA3tv/WcuTIkSNHVliJA/OxL/hZE/yZpLkXA7iPMdbEGFsFYAWAqUQ0GsAAxthcxhgDcDeAS4Rr7gq+PwhgGucqcuTIkSNHx8NJ50BEBSJ6HcAWAE8zxl4JTn2JiBYS0e+IaHBwbAyAdcLl64NjY4Lv8vHYNYyxVgC7AQxV9GMmEdUTUf3WrVudbjBHjhw5cqSHE3FgjLUxxiYBGAufCzgOvojoMPiipo0Afh4UV+34meG46Rq5H7czxqYwxqYMH25NZJQjR2o8uXgTGlvaqt2NHDmqjlTWSoyxXQCeAzCdMbY5IBolAL8FMDUoth7AwcJlYwFsCI6PVRyPXUNERQADAexIdSc5crQT89fsxBf+MB//8ejSanclR46qw8VaaTgRDQq+9wZwLoBlgQ6B48MAFgffHwEwI7BAmgBf8TyPMbYRwF4iOjXQJ1wO4GHhmiuC75cCmB3oJXLk6DDsaWwBALy7M3eey5HDJYf0aAB3BRZHHoAHGGOPEtE9RDQJvvhnNYDPAwBjbAkRPQBgKYBWAFczxjiffhWA3wPoDeCJ4A8A7gRwDxGtgM8xzCjDveXo4WBGuwnTdTly5LASB8bYQgAnKY5/2nDNLACzFMfrARynON4I4DJbX3LkqCRy87gcOSLkHtI5ujUYY/j4b+bi70s2pbomR46ejpw45Oh2ENf2EgNeWbUDX/jDfOt1uWtNjhwRcuKQo9sh676fk4accciRIycOObohsoqFcsYhR44IOXHI0e3QFXb+jy3ciNueW1ntbuTIoYWLKWuOHF0WWbiIrCawaXD1HxcAAK4667CKt5UjRxbknEOOboesnAPlxqw5coTIiUOObgdx55+FTnQFsVSOHJVGThxydDvwxV1c5F3MVHOFdI4cEXLikKPbQdz451xAjhzZkBOHHN0O7fVwzglKjhw5cegyWLphD7bubap2N7oUGNJZHuVSpRw5IuTEoYvgwptfwDk/f67a3egSaO/GvyNMWXPk6OzIiUMXwt7G1mp3oUtAFAulEhHlrEOOHCFy4pCjGyJJEfJ1P0eOdMiJQ45uBxW3kIaByBXS3Qebdjfi/+rXoa2Uv9S0yIlDjm6H7FFZc/6iu+GnTy3Dtx5ciOVb9la7K10OOXHI0e2gdIJLc31Ze5Ojmli6YQ8AoLUtf6tpYSUORNSLiOYR0RtEtISIfhgcH0JETxPR8uBzsHDNdUS0gojeIqLzheOTiWhRcO5mCtxWiaiOiO4Pjr9CROPLf6s5egpi4TNSrAnd3UP6xeXbMP7ax7CiB+6ic1FherhwDk0AzmGMnQhgEoDpRHQqgGsBPMMYmwjgmeA3iOgYADMAHAtgOoBbiagQ1HUbgJkAJgZ/04PjVwLYyRg7HMBNAG4sw72VFWu2N+DJxRur3Y0cDsgXAjX+9sYGAMCrq3dWuScdBx42pZQPitSwEgfmY1/wsyb4YwAuBnBXcPwuAJcE3y8GcB9jrIkxtgrACgBTiWg0gAGMsbnMd2G9W7qG1/UggGnkEgynA3HeTc/jC39YUO1u5HBAzJQ1ixNcN19HOtXEqjD4vebEIT2cdA5EVCCi1wFsAfA0Y+wVACMZYxsBIPgcERQfA2CdcPn64NiY4Lt8PHYNY6wVwG4AQ7PcUKXQ3FqqdhdyOCKrExvfj3RXJ7juel8m8C1mz7vz9sOJODDG2hhjkwCMhc8FHGcortqYMMNx0zXxiolmElE9EdVv3brV1u0cPRSZneDacU1XQufiySuLkDh095daAaSyVmKM7QLwHHxdweZAVITgc0tQbD2Ag4XLxgLYEBwfqzgeu4aIigAGAtihaP92xtgUxtiU4cOHp+l6jk6MdTv2V7wNlwWxqy+aza0lPPz6u9qFsCeuj9w8OXdzSA8Xa6XhRDQo+N4bwLkAlgF4BMAVQbErADwcfH8EwIzAAmkCfMXzvED0tJeITg30CZdL1/C6LgUwm3VSUn//q2sxe9nmanej22Duyu1430+fxZ8XrLcXhr8D/Lffv4p/vq3nHOM6h/TolAPPAbc8uwJfue91PLVkk/J8xKZ3cSqYAZ1zNenccOEcRgN4logWAngVvs7hUQA/AXAeES0HcF7wG4yxJQAeALAUwJMArmaMtQV1XQXgDvhK6pUAngiO3wlgKBGtAPB1BJZPnRHffmgR/u339WWpq7m1hEXrd5elrq6Ktzb5duivr9vlVL65rYTZy7bgc3fp30FmnUOmq6qPVdsasG7Hfmzd50ft3bavWVmOmYS43RScG8wV0ulRtBVgjC0EcJLi+HYA0zTXzAIwS3G8HkBCX8EYawRwmUN/uxV+9NhS3D13DZ795lmYMKxvtbtTFYRK4DLOXV7Xtn09I8T52f/1HADg8tMOAQBtqAhONHsQbchNWduB3EO6ilj0rs817GhQ7/R6AiJrkvJN3ngmuPT1dlKJphUFz3+YLW0ay7rgtjqZlXhFwe+0i77SqiInDlVENHB77sitxOTN+jy7+ppZU/Cnsy3IXBe/zUzowVMsM3LiUEVEdvU9GBmfgYnTYJrv7nV3TXDOoVUrVopj4fpd+NGjS7v15iTXOWRHThyqCC+0wa5uP6qJynAOWa/s2nvqYjCgVJzDlr2N2H2gBUC0YF5yy0u448VVWmLSHZB7SGeHVSGdo3KIbLB77sCtjCgnW+C9rssz+Ag5B4XOYeqsZxLHOE3o2iTRDV37zVYHOedQRbiyvFv2NHZAb6qDyObebfq6LPbtpbVdlVZznUOLTecgUYMuertO8EJruO58l5VBThyqCM/RjPPm2cs7oDfVQUggyxi6iml/WK7r4uuHiXMQITvBdfX7NqES46unICcOVYQXPH3b5OzWk7cCdaqeVxqv4K76uAuhTX+663pCQL7uf4flR04cqghXnUNPGNhpFyjTIxNFCGnq7erPmRwNHBJipa5+4wbker3syIlDFeEaTrg7j2vXBY3DSeeQsnzast/96yK8Z9Y/YsdeW7sTq7c1uDdWAXi5N3ASeVTWzMiJQxXh7trffQc239mV8w7Fx/m/c1a3rwIF/vDyWmzdGw/N8eFb5+CsIIxFufH6ul34v/p1ynMX/+rF8HvW8NQ9Yd3sCfdYbuSmrFWEa+axbj2wUyodXMREIrG9+Rl3ZX5n3V1ecstLAIDLphycOPeGELjxsYV+GlvbXfSk8BleaBFY3X50ReScQxXhOZqydtI1qyxQOcE1trThZ08tQ2NLm/IaG7SxhRR4be1ObA8C9M19Z7vfl0ytVh/z1/q5odOKlbqzQjrXOWRHThyqCM/RuqRbT16FHfodL7yDW55dif99aXWivMsct8UWEvHhW+fgQ//ji2b++x/LndsoF5pa23DNg29gy972+7K4jqemlrbYM+oJ62ZOHNIjJw5VhGphVKE7jOtSieGnTy7DZsmhTyXg4Pm6TXm7TY+kpS3dA9uwuzEVQSknnly8CQ/Ur8cNj77Z7rpcQ5F868GF+Nr9r4e/u8Hw0qIHSdDKjpw4VBHkKA/tDpN3wdqduPW5lbFFSYTrPbqUy7LQ7wniDvltdPwTL8calkYh/cgbG8LvnVXXUg7wZ9LUUsIvnn4b+5tbq9uhLoRcIV1FeI4a6e4wd/l6LesDlAtaGKk1eeO8nGkxbVW5w1pW3yYDl1JJ8Nsuxw7X50RZBp1D9wXXOdxfvw7z1+xES1sJ355+VJV71TWQcw5VhGvy8+6gc9DtTlW+Hi7iEdMTaU0pVgLiBKUjifE7gW9EWTiH4DO1h3TXH15ahJxDq2/csL8p5xxcYSUORHQwET1LRG8S0RIi+kpw/AdE9C4RvR78XShccx0RrSCit4jofOH4ZCJaFJy7mQKhOxHVEdH9wfFXiGh8+W+188HZAawbTV45jEXo5+DqBOdQRsk5qOoSGq1W7J00prY2pHUoDFGm8bVw/a5OK6Li46ytk/avM8KFc2gF8A3G2NEATgVwNREdE5y7iTE2Kfh7HACCczMAHAtgOoBbiagQlL8NwEwAE4O/6cHxKwHsZIwdDuAmADe2/9Y6P1w9WruzpYWSc3D0HNfBVSEt7rCrxTlwlMP3IHIo7HhT1scXbcS//Ool/OW1d9tdVznBnysX4aawcu7xsBIHxthGxtiC4PteAG8CGGO45GIA9zHGmhhjqwCsADCViEYDGMAYm8v87cXdAC4Rrrkr+P4ggGlUjtnSyeEasrs7kAZ7iJCoBEG/BS6nKatYrlrWShzlHOzVIG6rAvHY8i37tGW+et9ruOOFdzqqS3HkobtTI5XOIRD3nATgleDQl4hoIRH9jogGB8fGABB9/dcHx8YE3+XjsWsYY60AdgMYmqZvXRGeI/3rzuO5EnsAl8xmm3Y34q+vR7vcqosbNI8hzWLGOYDUCukOuvW/vr4BP3qs/Sa7aSA/1mpvAroSnIkDEfUD8BCArzLG9sAXER0GYBKAjQB+zosqLmeG46Zr5D7MJKJ6IqrfunWra9c7LSJWt+twDi1tJdz63IrM3ssyTKl+lPftEnhPsdrJA+yTv30Z1zy4MPwtKrGr8bx1IcXTrGX8th9+fQMWCWE1rNe5N6FFZn1HhcH7lVVZ35PhRByIqAY+YbiXMfZnAGCMbWaMtTHGSgB+C2BqUHw9ADEIzFgAG4LjYxXHY9cQURHAQAA75H4wxm5njE1hjE0ZPny42x12YoSZuyyC0M6kc7jv1XX46ZNv4df/XJmtAmkNJAV1aO9Co7pOPrRxd9wZb+2O/dkaqzDS7HTF+/6QEJDPfl37x1dkWGCv63cvrupw8U7WoIQ9GS7WSgTgTgBvMsZ+IRwfLRT7MIDFwfdHAMwILJAmwFc8z2OMbQSwl4hODeq8HMDDwjVXBN8vBTCb9YC3WFv0H39zBtPLjsA/lm7GFb+bF+MSdu9vBpAufpELRKWoSdDkojzNomD94r0LouurMPR00rWO2BiUlXNwKPsfjy7FP9+uLOffVmL4wE3/xHNvxdupuviwC8HFCe4MAJ8GsIiIuHvrdwB8gogmwR8PqwF8HgAYY0uI6AEAS+FbOl3NGOOry1UAfg+gN4Angj/AJz73ENEK+BzDjPbdVnps3tOI1dsacMqhHafqCDkHmwNWlcbzNQ8txI6GZmze04hDhvYFEBEy3ndX6OZkWlPWqD79BepMcJ0buv6lIQ5ZrY7KsV56KXfmlbY3aWhuxdubI+U4by3XObjDShwYYy9CPXYfN1wzC8AsxfF6AMcpjjcCuMzWl0rio7fNwfqdB7D6Jxd1WJsR52AmDtVyguMLkziheLwj3ve0kAeSSoQU7UKzWSu5PK3OZgun609WsVIalGN82Rw6S9KJvrUFdcEAG3cfQE3Bw7B+dRn7E0eeCCk9cg/pAOt3HgAAHGguj6LVBTUFf8DaOIdqjWeVEi8kDgLnsKOhGU8v3dyuNmJiJcrGTXB0xfmvVUinkN65WGkpUYbnZSO28gaoV42ZOJz249mY8qN/GMuYIN9SaDae+zk4IycOAfhOeGcgU+8IuCqkq0YcFLbhzW0+8RQ5h8/e9So+d3c9du9vgQ663WnaHXwXXPfbhY6QkZezBV13ZcJV6dtimimV6xzckROHAHXBQi0GYJt8w9OY/t/PV6xNThxsCulxQ/tUrA8uECdUS2tS57B6u2/p4xq2Ig49l5B1GrdXTFIdD2n18a4iIw83EppnL+siKi0q1RGBnDa4I4/KGoDLIsVBvL2hGdsbKsdJFAItnilvAQD0q/Nf0/FjBlasLyqEYiWhe3xSeyl2/Bt3H8Cf5qlzIKvClptMWV0Unl1xAdARh46wnCpHE7ZgifLhStM8HVHtAUaQZUPOOQRg0mel8fclm9AQRIjsjH4Otz+/MiSMYvtZuvKZ372KvwX5A+RFMKteOG03ytVu5aDuUceIlSxOmIzhgVfXGTcxNj8C+XClx3SCUwl+5qTBHTnnECDiHJLnNuw6UNa2lm/ei5n3zA9/2ziHakgW/vPxZUL76g40trTFkuToYOK+VH7JpgByLo+ivbvDqiT7qaJYyfa4Hl+0Cdc8tBBrdjTgW+ercyF4oVhJ10i6NtsLmajyMZxbK7kj5xwCRHMwOXhO/8nssra1V4opb+Ucgs5Vz6Q1eYxAuPzOeZj6n8841JBOFKTymk6D9s7/zrR+dIR1je129zT6G4Dt+/RE3ubVLo9dEwHfuNu+GduytxHn/NdzWLO9QXleJqp7Gv0510VUOJ0COXHg4GxnhQfPhl0H8JFb58SO2fwcTFxNR0C1e2VgmLc6EeHEipffiV+Tlt23PYOGplY8sXhT6n7F2mjX1dmgE3N1iFipjG3IRGBHQzMaW9oUYiV9Haf92L4Ze2zhRryzrQH/+9JqdT+k+lcE0WJznYM7cuIQIFyAK9zOPS+vSRyz6xz8z2qN61g4bYOwXt89c0CMNGds3NOx338KSzfuMZbpjNA91yfbSehcUEmF9Mk3PI2P3Don8dbkRXr8tY/hv556y7m9uqLvJ6ELANlVrLw6M3LiEIAPpUrLJFVOdrbkNKyDCJcO4jwzPZ72PLt4Pofkse4OnRPcjU8uUx7vCGzd24RWxxhaZNA5LN24J/EuVWv3r55d4dy3usDP5r5X1ymz6ek4rlzn4I6cOAToKNGNiji4ZoKr1mKp2oUpF7MM3eO3pDJlVV+Qvg0VxNg+cns9iSgB6jG/t7EF75n1D9zw6FJjOY51QVRbrbVSok2m/O6Kuppo6frF028n29MRh9xD2hk5cQgQyr4rvC60KEanrc1qc8i2ycvP6/qZTRTVcXBNulRJdIIuxLAvMJp4aolbWJTfPG/O8CYPIWY454KiZ166dAxPtYw6uiJy4iDBdfC8tnZnpvoLilXA1mK1FdI24tRaBmsq8UpTpNZKPALZoa8zKaRF7Ghorogs3fW9uRAwV2ulmO+MU+vJGk3QO8Gpyz++aCOefWtLpp50V+TEQYLrAjxn5fZM9RcLCuJg3ZkHn1Xa9dgsZripbZZ1K+LYkkrvct6tTqZvO9dRUIWwlsfFyTc8jZ8+VX4dhOn1btrTiK17m6zlOLQiUoO1Uhaxkp3bTkccvnjvAvy//301dT902NfU6mSS25nR44nDows34OQbnu6w9rKIMKrPOZh3eW1h/9QdbG/ynlj5Mj8DIiQ72EkkD/9Xvz5x7O+OYp40sN2uSqb/j6WbsWVPY+K4a/gM1k7OwZXbdj1ebnzk1pecTHI7M3o8cfjBI0uwQ/DgdR07WWXEBUVQItddULUsLeRY/ABiCypXozQ0uYU7jysj45/xcqY6nJrSQnwLaeJExftQ2fexIKPoMi3cE/T4n20lhs/eXY8Zt7+crEvbhv53lsdomwtasZKl3u8/vBj7m1stpewQEw11VfR44iBvGyu9AKs4B9vuOfRzqESHFJi/Jr4o2cRFnHO44JduEWzVugRRrGQKn1H+pyDrgVxbyDpU1m5P5quWh0VrWwn3vZoMVlgJAVja2+DEZLXCO1krzjHqHNJbLtl0L9qkQ5b675q7Br+fs9qpD90dPZ44JMwYK9xeFs4hnDCWcm9u3IOLf/Uixl/7GPY26mMe7drfbJSHfvS2uAf3+p3JxUxEm0XnID9jlZiqmtajsrzfNT921i5/6U8LEsdkvcdDC5IipUrB9dm7DMMte5qc2tBxDq59sROHdDqH2LXVNg/sJOjxxEFGpUUFKuJg9XMI1ipbzy745Qt4Y/1uAMA2QxycU/7zmVTy0NnLIiuOd7a2n12OmzHq7yqtqKlcaLU4JXJkGSsL1u5UxiiSCei3H1qkrqAiuvPkfZjWRz5eVUr0F1dsUwaSlKtrr05AJg7ygp5VrARUPr91V4GVOBDRwUT0LBG9SURLiOgrwfEhRPQ0ES0PPgcL11xHRCuI6C0iOl84PpmIFgXnbqbgLRBRHRHdHxx/hYjGl/9WNfcn/a702qMadq46hzSLkWl4N1miwMp4Yfk2vBrEUVqwdleqa1VQLQDioY6emnJ7pqRF+4SgiVnGykdunYN3FVF+Xe+5ImIlxY2ods+24HocbSWGthLD9n0RF5FM9qNu3/WZyov/44s3xn7rOQd7Czlt8OHCObQC+AZj7GgApwK4moiOAXAtgGcYYxMBPBP8RnBuBoBjAUwHcCsR8YSxtwGYCWBi8Dc9OH4lgJ2MscMB3ATgxjLcmxOS3rGO12Wdphkci0NroDTNlHmAP9WOGD/ys7LrHPR1lZt4E4CaYnwamEKo//jxN6O+9DDpg2w8oHtNRH7Yj8lCDuikWEmtc3BN0yubV7+6Kh7QUUffXTiTzmDa3BlgJQ6MsY2MsQXB970A3gQwBsDFAO4Kit0F4JLg+8UA7mOMNTHGVgFYAWAqEY0GMIAxNpf5I+Nu6Rpe14MAplEH8XbJgVDZGZ8l7EQW7+1yD/Byvg1V8iC1CKmC70K4n6Ik6pPzHYs4IAR6K6dy3PX5VmJaqO7CtIi6LLC2gIG6eF03PPpmsrACLdYcKNl1Djl8pNI5BOKekwC8AmAkY2wj4BMQACOCYmMAiGYW64NjY4Lv8vHYNYyxVgC7AQxVtD+TiOqJqH7r1q1puu6MSg8e1dy2Wysxp3K2djoL4mKEJFdkSudQCYKxZW9ciWrSOfSvi/Jjqbqy2yH5kQpiXbv3Z6sjK5RiJcUxV+dEF12RTpTEMwbasHJr3FJK5iTak0O6M8+djoQzcSCifgAeAvBVxpgpJrJSrG44bromfoCx2xljUxhjU4YPH27rshOSljTZrnNuT3HM1ma1Q3YD7duxmqyVQog6hyrPTlN+DVvfsnZdfCJNrXp/kcqYsqoU0nq9UKSQ1teXtAKUFMgx7jH9wJafkTyHdBZHbmKlHIAjcSCiGviE4V7G2J+Dw5sDURGCT27Ssh7AwcLlYwFsCI6PVRyPXUNERQADAaTPJJMBCaGSq0OQ4dz2fU0Yf+1jeHC+mzmie1RWp+o6HeR+q0QKrlxR1mdQCXqj6kvWIH7iGGhssedqLge+e9HRANKL9Oym1/Zj7Y2tZCMGug2Xy/gxiRV7ElyslQjAnQDeZIz9Qjj1CIArgu9XAHhYOD4jsECaAF/xPC8QPe0lolODOi+XruF1XQpgNqugwPkjt76EW5/zY8fLO8FyNLpqm8/y/vGVZGIfpVjJOtnSWytlQbmsoewNiW0q6na0iqkUjj1oQOZrsz4X8V4PaBLY+PWXjzqEzoaOYiUOaywwxbHtDWr/B137NsjXyNZL4u/3HxFJGVw2ISpz86x4c+MebUKizg4XzuEMAJ8GcA4RvR78XQjgJwDOI6LlAM4LfoMxtgTAAwCWAngSwNWMMf50rgJwB3wl9UoATwTH7wQwlIhWAPg6AsunSqBUYliwdhd++qQ661Q5FqRQVua4zbNNNj7Q29u1xpY2zF6mjs3z7q4DmHDd4/izxvmqPdPFRaykMmWtZKBB0/2M6F+nv064UO3B7Wc1m/XY0sQ5E8QxUI7wDS7ga6CrWCk653/qCJVqPH/+nvn6+jMRB7OOQay/RljsXZiCvrUFeyFHXPDLF3DtQwvLVl9HomgrwBh7Efq5NE1zzSwAsxTH6wEcpzjeCOAyW1/KgT2S57CsQGzvgtTY0oZlQZpKV/2CrcWQOKTommpyf//hJbi/PhmSAQAWv+s7zz2+SGNlUkZPcpVIoRJhSy6dPNZJtPe5903Ab19YleiTCuKCqOpyfeAP8tsXVuH6i45x7qs4LlQJocL2yyhW4jtktZ+D/rosHKzslCnWn2XO2a4Qx5MYCdnl8dkyM6ZF/ZqOiZFVbvQ4D2mZ/RSdmgAYR90pE4bgs++dAAAY2LsmViev95oHF+J7Dy8BoImjpGThzYORD9Y0k0hFhN7esldbviF4Dv3q3HZN//m4m8khkJyQruGabbLr8dc+hqvvTYai4PjBvxxr7xtR4j2ZXkecc0jiMxnDPotjwCRWSgN5IyTDxNm6cA66ldZllIplsuwL5P7JXIw4z4uFdMuc7t6bW0s457+ew7PL0uV9WL/zQJcMydHjiIN9x6E/V/AIl582HgBQIwy4C375PA6//nEAwDzRGcfRbNWWujAL56BadE3xaDhx6FOnZiblyberHeaWOnFM1Bglj2mufWxR3DNWFBe77BIZY/AkGbPc7rod+3HMvz+JlVv3xeospw5IrGm/kXNwZx0+e1e98XxdykWTw7ZJcYpfJBTKwjXKl5hEl6JYqT0JizbvacQ72xqOLQTRAAAgAElEQVTw3b8udu1mCFUGyM4Oq1ipu8Gq/DUMfI8oHFzi4BPD84rXq/RaFitOJXg4hzRTSEUHTOwyD6lRV6z8fkElbn5t7S68/M529K4ptEu/MaB3jZJwyXWKi2wiE5z0kh55YwP2N7fhwfnrrZxDVohtmsRKafDGOnOok9rgXaflZhOWR/JgY3bC3N58DknOQe5T9N3GOcgKY52PBN9EqMKf2FCJDH6VRg/kHLLvejyPwgGiKxdXrKrESiplrKNCut2cg333UonQAfJuV8yfIb6OGbe/jItveSk6ZRErqSCK+1x2iUSkTN3qgpVbzEEIn1y8CZ++8xWnusT7ajL5WTjVFtRpOR8Sh6yB94Lf8mLqIv6MbRAMxXVzw6qnEzmHmM4h+QTlqMO6hVzcRMxfk87Sviuax/Y44mANVWE451E0QPRRJcW6VIRA0abi2Cdufxl3vugrSaOBpe/dGim2vqqkS7RR3TpZTkXorMeE+ESKnkZtpZ9QohmiK6FLmDMbdQ5R2b++9q6x3i/8YT5eWL7NqQ8l19UyBWybjtpC+ziHptYSzvjJbNz78hrleRNs88RWl9VaSWjALrY1180hjiddaHIdcp1DF4DtFZkmVEFQXrq8a9WgVF2nmohz39mOGx71zSFdOIczf/actc727F5cltkjRvZTHpfZcBuLbVrUrXeQ8hYJSbt2cbGav2ZnSKQJ8efQq4wmj+IjkR/P9GNHhd/TEGnb6/aC2a/U7RgV0tG5d3cdwA/+FjfbdXkFrqasulM20Zb405afQyZOuuJin0UO1QU559AFYNc56EEanYO2BkeF9PIt+4yOMq2htZI7VISp1TGJTVa4egfXCnoNs0ghfR/ES8TuNGjk+AzAqYfGw3iF+TMYw0dvmxMXgwnoXeNGHFwU13HGIV7+GMEpT/eIL508NnHM1i4nwqpyZic4Y7VgjFkV50zzXVWXCvL8k5libWwlB25eN7fFTY1sxGBDrnPoArDKQw2nC160AOoHbfRdGTBKU/+9r6zVthtxDu4DTHWfO9thYdQeKw8ZMeJgaEstgiv/JJs6YUjsN3/er6xKypVfFo7101h2yeBdNpvtRufkUuLCouOqCMDtn54cO8Yve3qp2vGRV6XqlUkM0p4NVlSHcL9G/YZbH+T+HhAcCcXFXrVIy8TAJRFR2sU+5xy6AGwD2yRr9VKKldShMtQX1hqshLi1UprxxRiwaP1u3CuE8DDZz5djzXU1SRzSp9Z4Pq0deXyhSXcjqqWWP+99jXEfGKK4BdCEYX2d2nCJjSXmMXDdyYpgAD5w7CilaO/5t5MRjL/5gSMiT3Sp+l37m42bFVt/nHQOwmA26hw05+S5IC/WextblWWVTqiOz1u8Nu1i31Zmx7qOQM8jDrbzhgI+cfC/a9MQChUorZU0dQ/uo5dhtmbgHEqM4UO/ehHX/2Ux2koM336w8i78rr17j7BTV93S34OdrtUfQlGHTqwkw3SOP2/bAuD6Opj0qcKzb23FU0s2KcuJXbVxcErHS0XLXzpnoiD6iZ//xgNv4BFD6GzrHHIYCX8RlPnZRIuyWCn+e48Q+SDtjl8nfRWvNVn+vf+nzyaOmbILdlb0POJg2/VYruc206pB9o0H3rCm4NQ137+XmjjsPtCSSecgtrNux35t2AwO26KT1rHJXJfbrtEVMfl1TKyXzcQqfN7S/cgWKq7365rmlWczk8t5MQssNfglKuKgM9XV1bVNo2PhsN43sxsw8FznQXF9VZqTNs5BjGwr1qGat3IbLnoOleVfc2sJ+5tbsXbH/sS5XOfQBWAb149LHrcyuJOYSoH80IL1Me9W1ZzUTSydfuvEH/49GlgpxpfYTi9Hxam5Pl9MZYTrTpqpv5vKmY+pK8lifnvmEcPDXZ68G/0/KU5TyBE4ilnkUl8867DYb12so9jibrkpTzGjbcrTRPeFA//+wWR8qHLoHOL1Ja/4yrSJANxzQcuLr7jpEHftSj+jhLWSnTgsV/i4fPz2uTjm359SXpvrHLoYJgYhL0T8xWK7XlPwUPDIyiEAuqxvapgmXChWsrYYwcWO3LrYC3h781586FcvWtp06+Hbm6MYT+VQMLvazLugpuCFu0Lbbi9SNJvr1JW7ZvpRsd+hg6V0D2ksY1JxDgaFNMfphycSMmrf2SmBuDDtK1WV7xvE+NJVlcjnIFUi/hzRv1f4XWXFJB/SWTqJ4+FnT70VBqvkeG2t3iPdxceos6HHEQfxvaeNvsgnbV3Rc4rRrvaQTtUkgEi+mW4hjcrq1jjbYi9i6167049r7+4I/AYAi8mk41H+XjbuPoB1OyKfCtOSqjtXU6CQGNsVr25EO0qFai4ZhtCWOYcUxEFlQqq7PjJl1ddXVLAiuuI8+mk5RIV2q8Do+LghfRScA3DQwF7442dPwbcviIiwSp8gt/C/L61Wtyld+8H/cZ8/lYg6XGn0POKQYuAmUh0Gl9YVveycg9b+Wg++WDU0t2Hjbre4LjELDQeW1i4qsC9s2QKo6a9xVvgG5WROKEu60YJHoT+ITYfoereuaV7DBVHuk3AftjtSnddxHpHJsL5jYugJDt175oSkHJxD1Jb9+j61hRhx2NHQjAfnr8eG3Y04/fBh6FdXxAvXnI0jRvbT5BJxO9aeBb7rkYaeSBza8Zb4tb1qCpmzO+naN01Q0Sxv6YZk+u45K5MhGmKmgka5fvxk1igZRc/L9GzTimJd9RBZUVv00OLYKVdFsyvHF+2W/d8vXHM2Xrr2HIweFIlFTHmb/TpU9aqvUdsqxaEKWqcjmkVPTdxsUG08PIvMS1yoCx7FxsCHFDv6g4f0wWmHDlUrpFV9UhzUiZtckHMOXQBpXlGNxFJHE5AyDxTdDlw8alpMZOXyrv3N+ORvk8HdxBpMA7NcVhTFAmUkDgbOwfFtvbRCHb/IhdDJi23f2mL4TFTKXRE6RbMM10fMxT/8mQzvX4cxg3rjouNHR/0Vys9TOOmpREhaz3WDsyFHjUdJDlpzx6FYKeVAUBW3RSIw+SbooqYSkTULoaq+8JjlRZpCalQ6xW8l0POIg+UlnSLY4OvESp4HJyqjEmtox5dw/KUV27V1ysRBJ95yjZcv90fUB8S6Z7lff/fmPgFeXZ0uqmXYD8WxOzV9zmKtJIqVbOFAnHeDjmIleWHn7evEYx/7zdzEMblsqcSwfLM6emyoc5CeqvirWPCckyEVDYH8TFAV5y26eCv3rS06bdYKHikXeNfwIbc+t9JYf9GgG+qCtMFOHIjod0S0hYgWC8d+QETvSjml+bnriGgFEb1FROcLxycT0aLg3M0UjGIiqiOi+4PjrxDR+PLeYhymd9S/V9EY2oHvjAjqHYhT+zqxktDaY4v0Dki1EpuvF1NF38/5+T+19bnehymvMuBb+aRhQh5buNHefjsnVBqdw/0zT8U/vn4miKJm3XOAx3/LXsqh+MlyQ7IOQNW8tk984yKd/p/ZK/DkEnXqV40PXAy1RU+R70JdNrY4piDMqsW5wP2JpHMrt+5Dc2sJjAFTxw/B4h+ej961BaexV/DUHL/qUnlc+vnX03nux+vLfGnV4MI5/B7AdMXxmxhjk4K/xwGAiI4BMAPAscE1txIR3+reBmAmgInBH6/zSgA7GWOHA7gJwI0Z78UJprVIlF1eetscNEu78sOH+5PeExaQvYZUjKqwBYveVZu7if0yKUITJns6MZXjYLSJlY4a1R+APQplwSNlX2yseFprJaXOAXwxjVakV76jTG+uxSmHDsXhI/rBo2gM2Dgh3aIvx1ziZ50V0kG5LPofeZe/YK0+f7FL/X7yJYkb0dxIGqsqEaraePY20QR0855GTPv5P/Efjy4BYz4H36+u6M9HhwHvESnnlurS+1+NO43+URNOJBb3ynD73VLnwBh7HoCrDOBiAPcxxpoYY6sArAAwlYhGAxjAGJvL/Ld4N4BLhGvuCr4/CGAapdnypYb+JRUoWuBUScHrApGOL7v0j33jgTecW24rMaXzDCARB8VAGj+0j/Kcbsy5DsY2y67WdcLXeKRc6LVZtaSF0BWuOTJGDuiVPChAN8TSWF3xYrINe1IMw59xEjd/4qTwO7dKMnEuWg/p8Jr4cde4TLHfwoGCoHPgz15LHLK+U0V5LqISny2PPzVv1Q4wsPA5expdQqJ/ns7PIXns+48sif3etV/tNb5wvbjZ62FiJQO+REQLA7HT4ODYGAAiyV0fHBsTfJePx65hjLUC2A0g6XlTJphekucllao//egJ+MTUcQAisz4SdiqyI4wJbSUGxoB/O2NCIhKoCNUi64XKyvhx3e04h3aw7Oxdw3AXNAppHWdiD31eHSWeuNDYmue3ds1D8bhVsuloZMqarPBfTjwocZ2Jc9C9jlCJLuscDDdBjot5gthp6+Pn1SW+HHg9J5Esz+daSykZBsMX60btEZE2HpIIHRFxGWV7pCCMHFlCxnQVZCUOtwE4DMAkABsB/Dw4rhq6zHDcdE0CRDSTiOqJqH7r1qTIxgWmV1Sg5AI3tF9tGDKD23ETooFaY8m5LA4KPpCG9a/F6IHxna3YrGrAFSVLFlX9uvpMsImVosx35npqPE/ZFy1xCD6Nfg6qY0qxUnbIIhORK7SnlPXPc/0Jh+yRHDnBmSFnGVTqHDS7Ux4oT17ITe9Xt5jL7YaLsCH/Q/y8xt9Cpy5RVFej4BxC4hBszvi9Fjx3sRJjybJvbkyah8vY16QmDmL/zGIlaxOdDpmIA2NsM2OsjTFWAvBbAFODU+sBHCwUHQtgQ3B8rOJ47BoiKgIYCI0YizF2O2NsCmNsyvDhw7N03a5zkCZK79pCmEmKEwISxE81luTl4uQMrZ0UREhFRESEocJlT1DN/bQ4OOn5bZnP892szRrEf3ZJ6K4jB6LjPKHKOPFE+bVtvdGdl01gIx2GuT4+lkxiJRtOPHhQ7LfpGUYEOn5c/u1qrWTKLOe3p74fVXm+GRKzuPE5t2zTXqzfGZmruouVgrEsPZQfSpnsVBB1iwN6RTolcXyb3lY5vMY7GpmIQ6BD4PgwAG7J9AiAGYEF0gT4iud5jLGNAPYS0amBPuFyAA8L11wRfL8UwGxWQR7M9JK4Qnq/kCikd00BwwNLneH9/DwEHkVKYxtxEANu8QEsKrSjfkEop+6b6pxuZ2jK3SDCNqn4wmCLR18seKlsw10WTLXZoaIurpA29tANnsg52IiDZizJi2BYj5U4RI4HWrpgucmvTJuID54QTU/TVFL5mc1ZsQ2L3pU9zePX6QhOJKZSF8jCOfx/f3pNWW7L3qaYzsHFX6fguNFRQTQZF8WGrpG4uyLnYE1lRUR/AnAWgGFEtB7A9wGcRUST4I+r1QA+DwCMsSVE9ACApQBaAVzNGOOr1FXwLZ96A3gi+AOAOwHcQ0Qr4HMMM8pxYzrYOQfgtB/PDo/1ri3gi2cdjgnD+uL8IJcvIeIcRg6ow5uGQK5tSuKQ9AmIWyu5i5V0A72hyY04uIqVbBOqqNDXmOrnh//82nrleZe+yXWlgXbtFRTStmp17SZ1Drw+//PzZx6KD51wUOK6MPyEoX+6IHrheY9wzEED8GhoKmwqnVzMP3lH0qFSblG3obBZxurCeKiILDcbX7WtQVNbRLRUukJl+xYdy6//dTJWbWvAjU8uS5wTn7toshvjHLqZtZKVODDGPqE4fKeh/CwAsxTH6wEcpzjeCOAyWz/KBaNCmgAwP4cCR58a3/fh4kljwmO+zNP/rsvDwCHKJPlE9cVSiZ4J5RRiJQ1x0O3MRe7HBDtxUIuzZBQLatbeFv54zfZk7PvwWqXysLKTjFKYsupOy+GIQsYh+DJqQC8cN2Zg4rpQLMP0IiWdjotbs8mYr7C643CJyur3y1GsJCy+qv7rjBvk+moLnpIjl8vx2lwtzDhn1txWUoaxH9K3Ftv2qQNMis9AvA8x6Y8xf0jXow090EPaJlaSzveqTT4iUWlpG5SqWPKcCMX6JfxWLah855IQK2nab3Yx34C7WMkWj76okpUBWLFVbbrrAlfOoZxSyCgyqp0MacVKCi9lv3xwXlcfJ0pghjLqNq8++/CofUcBW1jKcqNRQEBm7ANZKlQF8fPri76/b+IwPPets2K7c96ePFbrioWwfy6iov6BrmCvxvLIIz0BEzmHOHGIynQ3zqHnEQcj55BkT3srdhj+uDVPFI64WClq5zFDUqG/KxLChzoHacHUyTxd48fb1l+dEk9G0VPrHFRxnwD/udm4G9kjtbWtpLyvcspzYznCLfXq2hVTVIqIPJ/Vq0i04TAE2NNxK2K2uHIoXwQkPKS15cxiGzm8BDcDF8fWCWMH4qBBvWM7da4/k6vtVcMNRNxk/9yRc/d+9fsh0utFxOcrfo9zDnqoxsqrq3c4WUo1trRVx6y7w1vsxFBZ3KiIgy+X9r/bBqVOIS3D9uq5DDYhVtIMmhZHzsG26POFRiyn2gEWC2prJR0YgGeXpTNHPu0ns3HxLS9pz5djURTNSa28g1bxqtM5+NAv/HzDkdz9y0H5ZGTxTiaJIzCUDP6b+zB2cO+gPjXkCK9Nrf6i39wW6cd41YXY7ly9EeOioYJCh6dC38BzvUG7KdE/Q/H5itZocc4hfv0dl08Jv6v6d9mv5+KCX75g6DGwYdcBHPW9J3Hfq+Y0v5VAjyMOVoW0VEAVslhUKFvFSm1J4iAOosF9aqz94m36dcSP6xZ3XSKjn116Quy3rf+cFRfbmXbUyES5Go21kg5Z2GyXhEOpIa0HRNECqOsij7+lpavyTpvFP+Ul6OeXneif5+XBEoXmXnsOhvWrdSIOrmRCZ8oqIxS1hR7SyTJfnjYRowf2NtYnbyoKwSqrCh4p3g/fgIkpeIGIOIgWZibYkhuJcbVkFHU6B8PDE89k5W555N0Xl6sjD1cSPY84GHZJnlJRnAQhetlzVuojqAJqD09xcI3iE8rSclEj3tENzlYN51ArKTQ5h6Eb42MG+f0T70O1QBU04TN0kKuQg9WlQZZ5x1+BHD9LFI3YFgpXVj8Rg0naYZ5yaOAtH1GHBHc5YkAvHDa8n0EJ7kYSfvChY/DG9z8Q64brxoRDdd/964pWD2k5qxx/jiJx4FeKbXK93dV/XBC7PiQOmrAYMuTAhonzsD8LIP6sbYYav/7Xk41t2sDHZ5/a9ueBT4ueRxysnEP0+1OnjFOW49ZGza0lreckx37BpFQUK40KYv9EKUDN/S5oFiTd4HRNaM4npq40F6uJ+XFVVdcUktsu04Q4RLKuaZfeoAIKaRPnoPM5ifqj+anhHGTxDoNaqazSiYXnHHUOvWsLoew9Ctlthlyfqg9Edk6kmOAcAushC+fAF/5dkq6A6xxUpuEqhP3TnTc40zEAQ/vWYv53z409a1F8m3xODBOG9TO2aQOfx/Kz6wj0POJgOFeQOAdt7l3yX7yLaGSekLdAVEj/aeapAIDzjhlp7RcAHBZEhHUXK6k5B1kuyjPa6cUVKpNCFeeQFCuZCFTv2mJsMnWWFIyi+E6/AzYvqon8CBadg5w7ulRSO8F5nj3oHWC2VhLfiW0nHfWPc61BecWdE1GME1H1QDZPLRqJQ/Rdp9fju2lXsZItuZHCiDDqA2Po16uIof3qYlydbRMmh0RJC755zBrxtj2w+jl0N5gmgucBaI3O68zaTDs4GTc8uhRXvncCgGiXTwRMGNYXy26YjvU7D+CWZ1cq+/XCNWdjaL9aLNmwB4P71OA3z7+TYJ91Y9PVWinkHBx3e36bycI1CmX+z//+trZd+X5VRE7mLvR1+Z9ZFNInjo37G5AwmfWcA09qo9llJjgo/5MT7EQ8J0QECeCcQxKmMBGjBpqj0HKY0mTqdDoyEVEt1iT8txFVGSJxUIleWzXUoXeoc3Aze9YlNwrPG+RKIsET+9YS5htnsZAeADB6YG+nMDEmcOLjKjYsJ3LOQYBsraR7H77Owfy2Vclx5IHfq6ZgXNAOHtIHfWqLeM/4IYIsXCYO8d8HBYuEq7VSkyXMhopAqga6ygnulVV6fYwsDlPtHp0jXraDd5B1MKGIp6QfKzwQI2PqOP9JizL/81ezVwAAXlget9ISZfWlEsM9L69BQ7P6vegWGdGpzjSmVIHieHcv+/Uc5TWys5w69pddhyFzDiUGvL5uF/YocqKoFNIyesVC6LvrHHQv1o98oAZDNBfEvvFxK3MQw/rV4fixA60hRWzgRE/FwVcaPY84GN6RzBHo2HPuIS0OyGMPGhArc+So/hg9sBdOGhcFQgt1DsJTl1vQ7YAicYd5x33/508DALRo6nnP+MGx35xzECOLThKCtxW8pHJUrZD2Es/2gGaB8+uI/1YRM9f5xMvJ1ixm6LjCoE4wbQdqBLPi7/xlUbI/id/+Eb6z3LynUdkTxnynQRWh9PumXrwOHdZXWV6FWHIaSQq/WuOtLo89VR+I7K53Mhe6fud+XHLLS7HAd2Ke9rDPmvfQOxArybpCmRsM+xi2oYbvL6E+WxJYB5E4cKIgb1COHu0nycqat4SDWx1aQrhVBD2OOJh4h6LkIa0T8/GorOILf+zL70uUGzu4N3oVIysDMbZSoldBXToWWpb7csgThyvLdNZK3NyQg+sclgrOOOLi7xEl9BSqgV6rSC5vWqxLLO4FrFoQXSfUuCG++OlLf3zNUjIJeUmL6xzUKNh0DtIJ/kr54pjYAAg7c5M4UJfxLI2y0tVQId6uJPbScg7mfshyc1nBDACD+wTBLTWOZiJEsZI4ZvWRgM0LtcmUVbQujomVLGJZXjKLzmHOym1hnCduUv8/zyzHowv1aYTLiR5HHFTviIfgTXAORrGS+YU3t5Ziob0BoH71zqBe0bIkLgfVLQ6c25DbnCuZ0vI5Jfo5DAuiyaqgsjEXF6+Cl9wRfvZ9ExLX8Cx5YvdMxEF+dCq/DNcJNXFkf6dySmiUwyadAzfJ1PWvThJV8XfLr5Pv1RPkMaZ71ukcZBNR0yId8+i1iIHuudKPxB9xNkxfXmiTsSQneNK4QQmxklzNDRcfG+rn4k5w/ueFx4+Kle8t+DmIY1YnUbWZ2hIIQ/qq50o881x0XBdGRxYhZ+EcvvPniCvlHM3Pn3470yYoC3oecVAcu+njk/DGv38gFlAP0E8yL9himN53Y0tbzB8CiDKGqfKwh5yDjjhoZJe3P/9O7DdfykUOZNu+eIrDvoLNdKNC5/Dhk6Iggx5RgtM5cewg+RLUFpILpooLeuIr70uUA9SxoMQySzaoM+655g8WsXH3AW2AtbgTnF/v9RceHW/TM0/4aUePiP3m5bhCVn4u0e7SvIjoMp7JzmU60jBmUG98+OQorYqsS5DxvonD4+XChVDRN4hiG4bGlnhHix4lFNLyGPjXUw8JCUjcC1lNlHpxayVJrKQTDclzTUaJMVx0/GgM71+HCZKorlSKnoPKzFZuUk7YlIVzEMV8WTi+9qLnEQfFM64rFjCwTw0AR4V0wMYyg873QEtbYP2QPOfFOId4v9KKlZLlgnoM4okzj4wSJTW1JCu88r0TcHKgK/E8UngRJ+vkil3e6vw1O9Ag+YBMGNY3nHQuJrnioYtuflF1K86JXkQ8ND8KEy7fisoJjpsbc9ic4BI7ZE4cClzkF79OVFqa7qXokVLE4mrm+NK154ROjYDdY5gjqXNIXkCSQnqTpFcpeKR4LlE9dUUvthmLB7fjC3C83fRipaBd5Vm/jmLBw2mHDo31rbm1hCeXbML2YJMl9k0XKUEmDll1Dhx3vriqw72keyBxSL4lPrm4/wKHTsXGFYMmK5mRA3rF8j7E+xB9lx2R+A5BzjEt20sv2bA7Iac/edygcIKZrJV+8bFJePabZ6HgERpbk5wDEeGIQFRTUCgaVRyVGPvnycUb8dHb5ipFReJO6qp7FyTOi3DhCFpLDLc8u1KrxFXB5DCmEivJnFPBwjnwnOMc/J2N6O9bkh02Iu4NLipKjcShQModZFKspK1CWY6B4d1dB7TlZNGIPtkPP2/vo1gfkLQaE8EX+wTnUKOOyqr3WjYTdU6EZN3DXXNWAwC2N/jEIR54L+ibNPwSYqUyeON8/5HF9kJlRM8jDopjfEeX3EWq6+Ccg24M3vapk3HzjJMCp6XkeZVXJR+w/NylJ4+NXeMJu9V3dx3ARTe/mLCUaWNRfSY2tFdNAROG9UVd0VNyDoBoQpdcHFWLj7hD+sIf9Iu+TjymQhpOWid2UqEQ2/nFz6kU0vL9mpzgvnHeEYlcAbzcUYEFy3UXHBU7Lz470z0XPXLyX3FVT4tilit//6q+nLToq94dYxA2JsnzsjWfX1/0XdbTiPcQcQ7x6zlBkXNDa5MROXAOvG2xilmPvxkrJ84Hrc4h+ORrS5rNiw4rt+oTH1UCPY84KEZGnHOIjuvFSupE5RwXHD8ag/vW+pyDokyTYrfOwSe/bIEiLlo8JLSsjBaTruuslUQUPH0cfH7ct1aKn1NZW7mIKMT+uSz8acRFafItm8JbhwthKdI5yGVMEVLrapJTSl5U6yTiERklmIlmseCF73XFlr3h8bFDeusuMUJcLHdrwozL/QP079g09r55/pGJ5yg+v1qjraaaKHEiLXM2WrFSvLoEeLc9UnP8HLE0oRqCxPvav84PVbLXEmanM6LnEQfFS+eDTBYDaf0c4L/8xxdvih3nEVbDcqQeh0eNSu6iZLGSLEcO03WWWNhfOR/C6IG9wnLNDjtMk6d3SehHgnMAElYdcuROfZtB/VLD//j6mXjoqtOVfdBBVJyL3r02u3+d57t4Tuye/IyKBT0hVI0ZXi4SNUjXCNyj6ZZrBLHSub94HgBw1Kj++NEliQSLjnDTdXgCwQTsITxUxgXKzG7Cd5lgDulbG3Iboue4iIg4IOxXc2sJLa0a4mAR8YR6L9I73gHAtz5wZPi9pOFq+CPqVeOhpkDaBEOdGVbiQB69Q88AACAASURBVES/I6ItRLRYODaEiJ4mouXB52Dh3HVEtIKI3iKi84Xjk4loUXDuZgreFBHVEdH9wfFXiGh8eW9RQgrOQSdW4knPfiblmp39jbPw9NfeH/4mafE9alR/nHPUCJwoOJmF61RQjiuk5ckkpgnl/d0p2Il/58Kj8LPLToyslRw4h1CxLnTyb196LwDgmulH4dyjR+L8Y0cplbbPfessLPjeebjlkyfjM6ePdwoTQILPhFzu8BH9MPmQuIOejXH4+nlHhN//Z/by8Pv3PniM8TqRK5PbEEUo8qLOMXFE/+C4G2cTiTv83ypiy9vhC68sZgH8cSqLCw8d3hd9auNRcFy5KLGY6b0lFNJazsH/1Im+eMA/rhQXn5/MORBRaCWmI0rF0LLJb7iNMRzx3ScSyvCwzuDTZK3klzM/P1FEFulD1GIlIkK/uiL2dUfiAOD3AKZLx64F8AxjbCKAZ4LfIKJjAMwAcGxwza1ExLcEtwGYCWBi8MfrvBLATsbY4QBuAnBj1ptxgVLn4EUZpWLnNZOMu+vLEVkH962N2dxzDoOjqbUUJhwR6/L75ZcLxUoJziHaraqUexccNxoDetWAglMupm+cc/jLa++Gx44PvEsPGtQbd1wxxe+vYqc7oFcNhvStxUUnjMYP/uVYJ12CmCbVTedg4UJi4qGkklB7nQvngOidyDvN71x4dLhBkKGqmpeTLVjkvjNEO/i7/m1qop6i5yWI/s6GpDioX136kGmmRx06Vho4BwZR5xD18SAh5tPYwX3w16vPCDkdm0LaJs7iJryuvgQqyyFxHIrvxzQ+xfenU9LHRGZFL6FzqEZmt7SwEgfG2PMAdkiHLwZwV/D9LgCXCMfvY4w1McZWAVgBYCoRjQYwgDE2l/lP5W7pGl7XgwCmkevWJwOMOgdJR6BVSMOuPATixKa1rYTt+5rQS6N4483qQvSGYiWmjtgZlfNPrtpmV15xE8BF75qVuS4KaQ4XXYKr+amtrn7Cjll8XIP7RuK9F645W9k+h9yEinMY3LcWR4+Odou1gdml6h449/Ovp47DBwIT2DCHtKDHibXJ+8KSVi4iWtpK2NPYiqeWROJM1bu4RBC3mRAfe4bFMPjkGxflqxM4WlGs9Nerz8Afrjwl/D3p4EGhXkYk4mri4H/q5Pp8k+Qa+VQVolzsKz9PMD0N9UZEblt8LTUFL2E9KCrtOyuhyKpzGMkY2wgAwSf3+hkDQMxntz44Nib4Lh+PXcMYawWwG8DQjP2yQmej7X+JDwp9bCW3EMGiTP/L972GPY2tCUsW2YKC7wxl7kBk7VVjiQ9Os2IvDn4frjuu8Lcm1wAA86wSyopz5YMnjA6/3/apk8Pv4aKgeNjPfvOswDfFh6ijmXxIZAZ88JBkZNcYVyZVLXJATCjPnfeickni9daPpuOkcT5x+NElx+MzZ4yPNaETK/Hfuw60CCFWEt3Gsk2+Evq//xGJ0HSJl1ziLcU5Fn05HkSSb1i4M+e3zj8SM95zMK9MKVYaMaAX3jtxWLzdYPzEPPEVVE4eUzbOwTXlLV+M20os1FX1ring9MOGhuVCjsBSp45w/XLGSeH32oKXiHUmEgs5mqsLymH9ZEO5FdKq1ZQZjpuuSVZONJOI6omofuvWdPmHw4pVCyu3b5Za1u2QZbHIaYeqaRkhGjSPL/J3e71qZM4hzhK3WKyVGFMTOH59bdELU3vawO/DOqkU1yXKWHZvkw8ZjJ989ISwrPj8vjxtYvj9yFGRWI4X+eO8ZORT2YPVJCpKdlZ/itfT2FLCH15eExRPXkBIcg51xUKiDJBcQOTnx/cBtz23MlykVczzuh2+xyzP1QwYFKcOjyPOsejHwM+CNKZHjYqHKfniWYeF/fQoEj/ZIgKLXDCHmPck7J+gx2ppKyXGadi2JvmSLgAfL3bDo0vx3hufBQB8e/qRYT2iYYpNPKvTw4ibkpqCF8ZgAoBNuxvxtftfD39nydVwf33lc0pnJQ6bA1ERgs8twfH1AA4Wyo0FsCE4PlZxPHYNERUBDERSjAUAYIzdzhibwhibMnz4cFURK1SvWpyQLHZcPTCKhXgsF1HBLILI351c8+Ab4bGh/eoSZfx+RbsZQKGQ5hOqpN7liQSBZ5mzgS9wVtl+Qqyk3+XpanroqtNxakBEZbHSEYKeJhaqOSizRaNgNPXRBFEmLxNa3vxvnn8n3NGpqm5uK2HDLnO/wluR5NLy8+NE5aRxg4ycw68+6XNVYjh4q9mmASb/io9PiabxwN416F9XTCyUvsEF31hFVm26iMAcoQI5RaKcidc/gbnvqEPAR8mSovr+7YwJeDgwroj6y8v5n39eEAkzagXCHuMcLHODE2dTuZoixQjm9x5ejL8v3ZzoVxrUdUCY1qwtPALgiuD7FQAeFo7PCCyQJsBXPM8LRE97iejUQJ9wuXQNr+tSALNZBYVwqqq5jJqk8zqWUg70pTffJGzc3YgH6qNBOHKARBykK1o0mZ9E23rVPYhER7VQ/ubTk/GLj50o3UfcQkaHpJ+DvoyLLkElkonOJYnDHxQ5E2Sk2X3Fwx/Ez/GFe58ix4CMv71hjo4pW2Yxw8J/1Kj+GNG/TquXAHyv+WH96pw2MC5qO1EGL4+pGy89Ifbb8wilEkvMCX4ZUXRfLRaRh6zDsJW05YdW+c6o1k45RLk4h0XrMHGTaOMcdLGVRNQUvJhuQ37Wrom5RJg8yssFF1PWPwGYC+BIIlpPRFcC+AmA84hoOYDzgt9gjC0B8ACApQCeBHA1Y4x7fF0F4A74SuqVAJ4Ijt8JYCgRrQDwdQSWT5WC+BoOH9EP9d89NwxrIFsr6QZG0eA8JkJMgMKhsvcG/Ek2f80O3PGCH0ivRhNps6TY5SXaVaw+5x87Ch+RvK65zoGLsg7SZhOzcw62cMix/hkU0nHOwf/c0dCsLCsiTQJ2YzbAoPk6aSeZBbLvhylkOw+qVxIcsXR1umxg0nEOyUB5MrjDpDzuo/wLUZ9t74uPFfui73+qiM3Fkw5KlBMXe9UcEDmHe1+JJ1QSrQhFsaeOu/nt5VMAAA8GcbqMnIMXV0jLRV2y2CXq7ADOwSqcZox9QnNqmqb8LACzFMfrASS8dRhjjQAus/WjbBDeQ9EjDBN23D7nEJ3XvbSC58WpvebdEiUneUJ+LejcPnrb3Khvihj9HvHFIN7gdy86OlHOBZ6HmONVUTPgXOqL5Nf2gS77k4iQOQBVfZ86ZVziWD9HPQsQJ65y7fx9iTszexobNWTfD51CGvB3uqKIT0eQiOJ6Bt0Cm4aTYohb7dygcKrjRgTynBBFZfy+fvXsCmN7keLappvgOozkPU47OgqGqBJTqRTcon77+r/E4xSJkW3FdUC3BsjBGE3jvqZIMeIrE5IsEVc7gnPoeTmkheVAJUsXzxs5B+GclrVXhM/QWv4kzPRUCwgprZXOPzYe595V/s51DqH5rGZBcanOpHPoL9nde9LzU9XDoSqmTF3qMMGeeXMzph09MvaO5ffDGTaROPPm7rh8Chqa3Z2Z5BSRpoW/QFz/A20ZIMl16RTSLlnxtARIcazg+c9YnhNiPCJXghRyDo7WReosgcl5LJYzhXhRtStyEaIEwXVXbxMriR7SclldJGYTcuJQAYhrgexL5so5eApPVXW5FOZ10nHVLp6LgeQqE5ZNjpOUO3LxPuoyirnsnHU6h5evm5bQsxRIL5aTFxidqaYM1e5SxpyV2zHt6JFmL+7Q1j2+kwSAc6Xdog0RNxX/1IuVIq5QL1aSIpBqnuPaHeqUn0o4rH/8nfHd/jc/cETsWhJMWW1QiYHU5ZKLvqm+GHEwiJVU9Z0qREEWN3WmaLUcjDHc9txK7Xnfz0H/zlx0Dl+eNhE3PxOZMKcxWc+KHk0cEqwnuRGHohffvek4SoJ+hyw0qaxDtYvnTmuyAly+D/dJ6hObKJ5TdrFSZDESPz5KoccgIu1OX74X1eKnWjdNwQw5eGIjF+/XWIytjEoHOVwzv2fV8+RcoUn0xPtnWmjSQBdrSNU0V0jzscJDYfArff2a66bEL+cqTlHFalLVJ4pulGKl4FMmDsP61WKEYOEncg5vb94LG15asT0WZUBGbaF9OgdfnyPV2RkU0t0N4muQB7O8Q9YN3oJHMXmp9tVSUiYsD1lZ9MCh2sUXgkWV+0xwyLskPll0YqKobX9xaTMsWmIfjXVJNv0meKR/tjJ9UlWnmvgNTXbicCAkDkL9ib4lWblspEHgpkJzR6kNAf671ftCiP0TI6hmUWaKbfp1OJTlCmlpIyEq2dOa6ztzDhbrJ17u+r9GIeyV1kpB/+SglHJmQ1GCIItEVRDFQlMnDMGca8+Jna8pxE1ZkzoH8/0VC15ibVLF3io3eh5xiMkq4+dEKwVAn9i8IMnM9ZyDg1iJ1yEdV8VP4jt9kb3kx2P948TBknjeC+zU+ULdHgNi2YbchMaWttDKQ4aLWEklMpAj1Ora1dUZ1q3QnbhsiHm+ZVNdJp0DBRsJnqVOn6IW2CZEoG0HbQgXGJnr4oEFRRQCsZesnxJ1JHJX/vmts5Tthh7NjibUsshw4oh+uOC4yKueT5XX1u4Sjqmen5rYyGVF/w1THz9z+ngM6FWM7eIPG94XBw2Kh1CXneDkKm1ipRovSXhznUMFIL4GVYwb8bxJIS2e0/k5yP4QgLtZpJwXmF+rlMEnFOu8Ds9oosgtXzgXpBuiGmmTVBfngOxl9xgiVLoopFU6BzkIogotithAcn/VPhzml3b9hUeH+ZZVkPM56DLpLVi7M3xfJi5ub1PEOXxF8C5PCx7jqElaLOUshEAgVmIMbW2cc0hyvLKosH+veAh7Dn77rtZKoljpnKNG4HefeY+ynO1YxDnE21XNHxeFNN+sibv4Yw5KemXXFL0YtyLPYZu+rOAlc6p0hClrD+Qcou+nHy7FfJF0Dlq5uKvOgdyC86nqGNSnNlFWbpdDXrz5xLUNIJ7UxLZzcVFIuwY/s9dj5xwOKCxxTNY5PGieyutVrl0n5jPBZHYqtttaYlpRX8GLmzuadA78fd1z5VTnIHsqcF8OXTbAWP+CjQ530gxzWoT9ShppmERjgF3nwB/B6+sijkBVpU5Ml6gv+JQ5h4K0EfMV0v53cax8bErcT6gQGJyI8yyMNSXApnOwcb3FgpfYUOScQ0Xgv5nffHoyzjs6bn0iJ/vRhQGQOQcdXJgEVaTIMYPUmb10zmPJxECcONh0DlwhHXAO7bB84WNXlZPahIuEoHuAP+gf+dIZeO6trfjF028ncvMCwLJNexLHGgycg/wUTG8utHxpj7wmrCtOaFpLTG8RRhQTV+rMQj2K5NdZYvKIKHiEmgI5KfN982Mk+id6fetiHyXqCnUxdkUsADwdCzXhxiWorZUCsZLEOcjzTRQv86LPf+tsjBsaD+LoKTZrKuJf9KJ31lZi2CN533NdmAw+Hn78keOxcuu+2LmOsFbqsZzDuCF9FLJGyVpJs6P2PH0GtVg5FxlSuLt0IDZESuWhLqQ211sM7ZvkQvzrENM5uOCa6Ucqj6ssRsRQzXK7HKrgaCeMHYQBgVObihhySxkRJs6Bc1ZMEu8EP6SygTw8RUgD3SIo67Zb2koJz3eOAsUXex0R8SgSQ6h2x2lRVywkxErK/gVOelG+kUAhLXh0y4usXjTmf9rHnV/w1EMjMZeqTqX1l0o8GHzKCumvnXtEolyoJ+LGGorXxjdrIlFUjQU/vatf5nsPLw6j63Loxq5HfnBDP+GWxDnkxKH8iNjg5DlZqaZTRsm7A+3C7kIbNP1QoeCp29LFYSoWCD+/7ET85YtnKOvjMlO+2L7/CHMww+PGDMAXzzrcWGaDYBc+oLeaMY3lcNY8JDHznQyVst600PA2nlm2BY0tbRYPaf+zfo0y9qOmfjXkJDStbXrOQX6HqnvkdbaWiXMAfO6yUbNzjffPQ2spWgj5HPj8mYeid00Bp0wYgsF9a3HR8REnqCea+ncrIvKkFhdfRTmNw2iyXf9TDHf9mdPHo7cUeqUYiIFKpShciKo+LmqzhrMRdIUPKQwxdMShrRTlXJfv20UP2F70POLArSsUU1q2QjKFz4jVqWlLvfBJC3loThjVoovvrhUrSSOHD6gCET46eWyCHQ57Egzagb1r0LumgG9PP0pzJ0E7DiPyi/cuCL/r5KLioqGX1/MFJHkurfmmKCr8w8trjISEL0Tb9tnjOXFo7yH4/OlTfjrZ1lJJG6IkcoLzofdWp4hzKANxKHhxM8sLjhulLkfcQzogTAGRO2ncYLx5w/Qw8KPYJ6vOwcKduXo+q449+1YypD+fj2J9cn4VABgzuDda2hi27G2KxGga0ZXvlGpXKJcMREZnFVliEeETLztoYK9M2f7SoucRB3CLkeS5ohQgy2StFKtTMzZc5i4fLD9+YpmlJLd+UhxP6BzidevAc1y3tjGMG9JHW/7/BUlrTCqMrYJ5JYdOIV6IEQezXPq1tTsT5044OBJF3fLJkxPnZYjv50ePvYnv/XWx8hwANGZIoqJ7LPze1mz3dTYtbQw1OoW0RBxkJanYVqsmcq8J4xRJjwB/TIk7ad1Ggm+cZM5BRjFGHPT3AKQJn2HhHBTHtu9LjseoPpE4JMcoF2k2NLfiV7P9OFFK7sSRyInJwcRnEnJGiucgO0yKG82vnXdEZsfMNOhxCumIc0iipkCxF62j6PJAufpstajF5f3ZHNXi7brpJkLOwVI3VyLu2N9s9Ing8k1TfXsOJENc6+SiA3oXQyWcbXf55OK4w9+sDx+HT7wnCrw3RKNPEWFag2Qz5CZJxDKoj9ocU4R+osbrbm3Tcw5JsZJO3GZ2ptNBzl0utivqHLRiPvKd4H4/Z7Wyv2J9NrhaK6k4B1eFtKlu0ZR17OAkMeT3sG7HfmwKcomoOAf+Kjmx/uPnzDo2xlhsvBc9P5S3SjG/K5hPfIMlNp8qsVU70POIQ/Cp5BwKXsxbUbcj4BO3pkC4dPJYZYgIQD3R5HbT5iFw80B2JQ6EF1dsA6DfWQKReMg0KPcpPJR1nMPw/nXYvMff2enl9f6nbNlx9OgBMeLs8vjSmKXKnMNTX32/9RpdH2QfkxajtVL8t0nnwJFm7Oj0lx5JxEFTZcHzOYxHF2409k+8P5M5LhDnHMQQ3DJEj3BVjapd/Tc+cETimErnMFoxdzkhEEeNqg1+jK8TpnhYgE/URR+fYoHQ3Kb2UP/W//kJwt4KlNcuothyo+eJlcKFIvmEazyKsbCnHaZO/8knZUsbM05QJ2OlFG/a3705lAveqlWsJHzftV8vY3dxuJGLHD6iXyLgnqo+XZBAnc4h6bBkf34m2jDlkLjD1xDJv8SFc9C9aNlEtNVkreTIOYj3a+I6v3xOnJvVWTbxRd8GHj4jbNtBsW63Vora/dAJSeLAx8YWQWSpGtOqds48YkTimMqUVVmfwsFPp5AGIk5Eb37sf8qK/9YS8z3jFRIKbj7O+yrW3J6QKWnQ44gDh45z4PjISWPw5XPU3qfipNTtoPw2ko24iJFUjjR+fa6xi9w5Bw5TWb6Ym9r+amAOyM0Zb7j4OMeYTOa+yRNBlz7VBH2mPuB6KRfGeyfGHSNdWHhdidpCXNlptFaSjQosC42pDAAcNqKfU32+WEkIV62pz6N4sES9zkEg/LadtLAmKgMsavqrq4/juxcdrbZWCj7FDaAyND5f9AWiqRYryZyDosOInv11f14UO97cWtJGKOZiMT7exbpz4lAhmHQO4sQ9YexA/QQVjqfhHEYOqMNZRyZ3NK4oUDI/hLJ/grWSCeJpkyUSd6YzyXH71hVx0MBe4U7H5MEp9srqKCXdb9+6+ILbXs5BxRWJOZp1z3Dm+w8V+qCu+/jAh+P4Mf5nU2tJGzDNlYMUS6XxczBllnMRKxEBb6zfHf7WPdOYmbJusQx1Dubgla6ez2K5I0f2x2ffd2iijNgfG+fAj4nPRefnAET3oR/L/vE5K7clz3lqIxPeR74miXWXw0HTBT2POITWSiqxUvQ4CgZRSpxz0E9Q+dTXzzuiXeaHHhGWbkh6B+vaTcc56MvxBc3m0UoU2cybokaKtdgWEJk4yDbpLutj2pAeLmIvMaaRKbzIaYcODS1imlrbYulHRbj6NLlye4lgjJqyuw+0YKGw6Ovu5TnJNFTWBXEUY8RBJxrzP23BK1WX28RKdQrro7A+JBd9ZYBLBXFw4xzs3J6In116QpA4LEkdeJ18TYpxDi5hdMuAnkccDJxDwXHRF8uZPHPlgddeKwMiYMPuxvD3PVdOxbzrk9la+eC2Egehe6a+8cXSmvPXixS6rrFfbEpLed70rZWyyrWTc1CBPzfT4xPbNXXBtzDzvze1lpSmk2KbAIw27JmJg6aTsj+H6xA9YmQycqvYJ9Ozi3QOInFIviSV+MSkGAbMmxLerij7N+kSYsTBQJQ452CzvNu1PyKoJ48bhMumHByYMCevCevktyO8mCxpRbOgXcSBiFYT0SIiep2I6oNjQ4joaSJaHnwOFspfR0QriOgtIjpfOD45qGcFEd1MFTTiDYmDogUxFpGJZRcHihgUzFTfZZPH4uJJ2YOk8XZFlrhPbQEj+ietLVx1DuIu0YU4uDgtcTmtq3u/rlUt5yA5LbkMlLRTqehAXEXCauoDT8UK+AHudJyDOORv//RkbX0li5I0qi/+2zmFp0OZr0ybiJED1BZ6/NmZprAqXafqHQ1QhElRi5Wi77rn67frQ9zQKXO1B++2OSZu0xOlZgvnwK8VF/VrLzg6rEPF2coRhMWau5LO4WzG2CTG2JTg97UAnmGMTQTwTPAbRHQMgBkAjgUwHcCtRMTf5G0AZgKYGPxNL0O/lIhslcw7ENNkErkKU3A7UY7//X85tt2RFP38Cy7lePsW4uDIjvN+28MdRBW63msasdKqH1+Y2DkOF/QDV599GJ7+WtL0NGukWNMcdOUciPwx9/bmvXhr815s3tuoLOfiGAgA9Wsip8A0Is1ybrfMXLX9vfPLxRhDqnfUr66YMHG1KaR1nBmAcIXdLwRpNFohWay4Iic4zjmon4tqieBh0QseKZP98E0gfyxeV+McNLgYwF3B97sAXCIcv48x1sQYWwVgBYCpRDQawADG2Fzm85Z3C9eUHVFM/eQ5V5ZdnAAmM09XTsQVi97dHfttUwraLKPE++1Tq99xhWIlR49WAOjXSy8aEfttk0vzHVSf2oKyrGinftjwfpioEHekpQ3vbGsAYL5f18Wce8fymDpiQppYfY6+G5MOjrKW6RzqVH1y1nU5jFOTlRTfibskVBKh4zR1oWF0x8ycg1+uQeQcDAppG3FI+DnofEksImqlWIlzDoqIDl2Fc2AA/k5E84loZnBsJGNsIwAEn9w8ZwyAdcK164NjY4Lv8vEEiGgmEdUTUf3WrcnYKa4d1kF8h66cg2mHHDPrS/Gkp2vi28g4XDJX5OALg00ev0nQX8jiGhG1RT7hzf3h7Q3pW4sBmkQvruCPmDepC2NODgS9HTYAWoiP1lS9RwAYs1ojxcIqGDp83jFRmHkT1yqfcfXEdynloo8zEWT5UXzhzMNwrhQ+PyorE7lkGfHZiZykrl0xRLaKwPLn/9iiDdq6gIhwtZTMnIPp3fPUvzImjfM3AR868aCg7uicLUlSudBe4nAGY+xkABcAuJqITO6kqifEDMeTBxm7nTE2hTE2ZfhwcwRRLQw6B5eFBohPXpNs3cVbVMa870xzMnedc+05yoRAflvJ9lV4S0iefvyYQdpyoc7BkuuWt9srjfhMs4jwd8Hv4N7PqkMTiNC9szuueI/y+K2fssdl0oEcOQce+da2Nsf8FwxlRSJu8rGRx5vr+DtqlFrRLMJ146SD/LyuOuswgx+G/Fu1mEffjcQh+LQFOOSL/tub/RwKfIFO9k3iHLQiUv9T6Y3tqXPDjB3cG2MH9w4zDIr3fdK4wYnylUC7iANjbEPwuQXAXwBMBbA5EBUh+NwSFF8PQPTuGgtgQ3B8rOJ4RWAyZXUVK8VMWQ0LsKtY6TeCAtI0uEXRj0lssyQwd3VJ0sNx3YX6iKycAFpoQ/j8TOIOGaYUq4BvMTJqQC+M0ChARegWJhWHNXJAHS48frSidHqYlkOCL2Kxrc1x/wB9YdGU18Q5yGvox6aoHStF/OQjx+MCh2diFrk6iKXI/Dt+Lh3n4MLVWImDdOzy0w5R9y0o12LROfDjxypSiOqyO7a0sdjGkxPtPrUFnH1Udl+pNMhMHIioLxH1598BfADAYgCPALgiKHYFgIeD748AmEFEdUQ0Ab7ieV4getpLRKcGVkqXC9eUHSZT1phYydFayUxE7Pbyfjm3heHy08YD8BeF/gZzRy7XXvyu3SeCw6g7cVRIU0gczAuEWIuOGPLuNLe2WeuLrnEfzuVMs2hWSPtGBLc8uxKAfkGPx4vSVygqXG26Do4Vsy7AuceoxTYieDpVG8rNObhy6YB6XsY2YYa6VO+8T21yHrn6iMgmr7YoxHLual63rEN4e/Ne/O2NDaHuC0AoplX1t1JoT0sjAfwleNFFAH9kjD1JRK8CeICIrgSwFsBlAMAYW0JEDwBYCqAVwNWMMS78uwrA7wH0BvBE8FcRmALvicd0IZOB+AQwe9+6LWomwqFqd2jfug4J2ctRG7rw2+Tm/meaSLNna0RoJEwo12TqaZJjyf4S7YFZMRwnqrp7cd2YmERJsXZjfXAlro7j1bhxcrdW4jA5EcrPYonCCTTmtOhgki1C5RfhGueKV8cj+dqCKj7/dlJP6lHSf+hvbyQFJ5xjdEnpWi5kniGMsXcAnKg4vh1A0jPLPzcLwCzF8XoAx2XtSxqYkv3ExEomObLrYu64WrkuprxdlzDV5QTfcdnWmFCslGIHbws50NLK0KtX+TmHciTK4UhjUnr6YcOU5eLWT9naircr1Od0hVuApYEQxAAAElNJREFURcBmhOEiVpIWX5N4TDr1zLItiTLiPLNxNUTxDZ1qLsv3Z8rnDUScg0skXRkFjxIpabcpclHwsDEuKV3LhZ7nIW1I9uMa8dJ1grqWc12oeH1OkULLCN6uK+fgyjH97NITrHU1G3IgyEhjLlxOxsvojIa4THnWh9V7IFexUpbdvWkzc9PHo/2daZEW2zV5Ic9fk0zOZMJvPj3ZSJRcNmLieDtgSHlKRE7OmfL92Yjmfa/6Rpi6uk17loLnJTgHla6wV2CiqwvKWQn0vHwOjjoHm22yXJ8KzuKQlGy/yey0krDNU251sVuR+CeG4KEdpDFPBQR5bktb2YksYBZleGQ32xVh5By8+BiRY0NFbUZ1GHfmjoTXlfidc1Ski7BtiLhsXJVak2PjHrWTnwhXBTKQnBvXXZA0nBBjou1oMKd3rS164e77+guPVpaRiYNuXInhMAD9uzET+2TMMm74IYa89zzCshumO0ceKAd6IOcQQPG+xJfougDPPFMdARJwn8juCle/XLllyDa4Zh7j8uDV292spEy1iTqHtM/HBbpcE0A6aytbuzLnoMvnIIZ0MO3MuejMZnLqaroqdt103+KiZOqfiyms2DUrNyo92zMOT4rlaoT+2Myoed8PG94Xn3u/eu7KxE9HwPYJntaAfh6bYmUVKGnK2loqYWjfWvzj62cm+uUq0i4Hehxx4Ns4tc4h+m4yFRVNOk822By7yt5dZeWReMepuHESp8EhQ/vgfROH4b8+llAxtQtm/wD/08+77Poc3SfOjz+iF2mZYhupYM7pEbfO0oV3qF+9I/xuem+6HNSJcq7GEKKXvGFDFFuADeW+cOZhqdo8/XB1Qi1VWW3fHO8ViDg301hx5Rz2N8eJg26c6uJQAerYSi1tDMccNAD92+lI2l70OOJgtlaKDppMRV1j9bimp3Rd0/hEcd0Vlstcs6bg4Z4rTzESwiww3Ya4Q6oE52DyJ0mbc8McAI+cwoV85dwoBLjpvblyj67PQixn0mWJfhKm/hljGwUQu2YKd+H3T742eV/iomybcQcN9EWZxmjKBTfiIAfS1O3q+xs2mqK4jqO15G6hV0lUvwcdDJPOQQzla3K6soWu5nBNyuEa+4cvkq4TvzMMMBX4YzFb+YjEobxWX9+9SC1rFvG5903AWUe6eeHbrG1cYuG8Z3yUrtRF52C7U9d3HwsDYiA43xRyMpvEJC7t8nY+eILd4U6261cm3XHUAQLAht0HAADrdx4wlpsqvA/dZuzgIX1wpCZ0uQjT+/RUYqU2lsocvFLogQppvYe0K0xsoohyx0Dhg9S16wcP1it8qwnOeZnzEUTfXUUprhyVywJ2/UXHWMsM6lODXftbrHkVeAA3085c7JNJ6RjGzbLcQjl8FnT12XbCLljwvfOM9XDIMbV0/f36eUfgF0+/beXqXUWUUycMwTxB1KcDJ0zHjdE7EJq4IzGH957GFixYsxOtJX062Y5EzyMOwWd7Hv2YQb0x7zvTMKyfXjQB2PMfpIWrSemnThmHe19Zi187ys4XfO+8dvctDSLOwS02kLOpr+OEskXbdAVvzbQgEqLwClcZ5PGuimHO3do8ZbPoHEwQN1MmzsF10+XqqyNHCx6pyF8CCAH/LPW5ilrTOl6edLBe5Go2MIhySH/1vtcxe9kW9K0t4BhHb/VKoucRB0PgvTRx/11i/bTYghFJ4LmGdRgaECPb4vYfFx+Ha84/CgMd/SE6epfCTfdcM5m5T1RH4lAmjo4vhDadQ7MlvIJYlw37m3ziYFqg/T65PbMs0gubd/nHpxyM9x+RMTCmBP7Mxg3pg9nfONMqYsyYuiOBBkHZPNywCeTmrIcM7aMtIxKk+2aeGhNZieEzlm/ZG7TtHjKmkuicQukKwpTsp9yRcD98UrrMb7b1YcohgzF2cO9EAhQZBY+cCQOQLtxFOeAkVhJGprNC2nGBbTI4SqUBF1GaF/2IGKWxqtHhpHGDUFvw8MWzzFZBziG6M4hXbeaUN156Ai5y0Ce4gL97BmYNU8LLueA2S0Te2YIntul+OXEYpYi4yhFPDubFk4pRRBx2NUR+E2mdCSuBHsg56DXSWTOG6eAaJIsvkiYTQQAY3LcWL377nHb3S0aacBcusC0M/Dmb1pi4o1T7OYeiFyn+msq0C2hs8esxyc49gTik9Z9QYWi/Orw96wJruc5qjJAWru+ejxfXKTzYItYaGKQo/dq5RxjLcV+HQb319ZkiL3gC5yByK2kiKlcK3WMEZYBqw+RqelpuHHvQgP+/vTOPsao6A/jvYxgYGbYygIjDwGBQS0FApijFLUYUqy3EJXWpUDRRW02laRqX1jSNNWIX60ITSyxU7WK1te3Y2hq7qNVaF1QURMSFCJZIiSIjRhT79Y977szlbfe+N/fdZeb7JTfvzpnzzjvf/e573/nO/c53WHbCFG4+a1Yqnx/HiDbItYsqp8nq3vow8rRS70NZn7iqJ91XXM8c/FQNleLRhZ5Q1iQ9NP9axK3bpPH7H/bV9C9tqY1zgvjthBnPU1zq8rOPiJauYvh+1Q0EfQY1DOh+JhXsem83y4qDfug5eK+lvjL12H1v0czxHFAhTQR4I4tlISOUehJ3htew2PWPuz2HaNFKcaxzaAnMG8dlHGa1jeTZN3ZWXDwWHPjGmSY8DD/iqdwOa3nB97bCjMOMVm+zqk+3j6pYz38OGJaGYum8SZxzRFuoN+8Tds/7FN7LTY0N3R5okCiRXPUm/R4kTKXNfuKeVgK4MSVvIE3CfgQjTSsNqH5aKar64jIOP1s6hy1vvx8yBx9djqOnjC6ZHqIWRgxppPPSeRwcIQ4/y0T1fI6Y3MIzV88PjYL68rEH8Y3fPE/7mOaK9UQksmGA6Ia/8B4YMqihaKU1xJtSvlbS70HCRPEcwqKGjMqERQ35QVyVPYfqvZmoixPjilYasV8jI0LuleCleCskKd2dF4RvhVoNh7WW3/o1L/j3UpQf3yjhsWd2TOCM2a2xe8tRjVjh1KJnHIoDJJLMoVSOfvfMoVL6DP+Zw8wJ+f9SReHui+bus/K1txTGpJfDv86VvgBBL+6fm4o3SSnF6KGVfxyuO206AHtKuPH14sVtPZvTdH1QPELMEyJwfEJbVPrsdg98/QfEcVCPjbJqnTJsamxgz97/7ZOdAWoLMY6bfuc5nDxtHAfvP7TkHKH/MCsLikmCOe2jmBMyR1sNDyw7pjtWuxIfR5hW+igwug+L+tq8/JRI/Zvm9vA9OELm0LjY8nZPmoY4r3WcfHfRtEhTKK9fF+06x8mh44YzeUwzV58avmI9TQY3RBsYFaZs9wdUJ/7okX3Kv39GvEkua6HfGYeJLc1MbCk93zjTJZaLawFPf2PCqCFMGFV+MZDP/sObeGvXnooRI8HnAstPnx5L/6a3jqDz0nmJrj4NJqKbe1DlDKRp8cUjJ6bdhbI0Dx7I379+XNrdCCXMc7j3K59hR9eeotQ7fshqMHT1czPGM3W8rZDuRkQWADcBDcBtqro86T7MnDCSl65ZUNWDKKN6blvSwaObdlRMPzJpdI8BP3RcfF+UpOfhgyt8jb5H++hmXt+xO/SZQ7mMxoUDpHXfOSl09XtSZKIXItIA/BiYD2wFnhKRTlV9Mem+mGGoP2OHNXHa4a0V6wxvamTy6GZe27E71v2ek8bv+Ypz+l/UWn/g7ovmsnbLzlgWOEJ4WpQkycoD6TnAK6r6mqp+CNwFLEy5T0bK3HPxXH5/yby0u9Er/K1QC9MyG32DMcMGc8LU2teSLJg2LsbexEtWzNSBwJbA31uBeOP6jNzRMnTwPovX8sgNX5jB6sc2dy/SMowgR05uYfPyU9i+64O6RFH1hqwYh1JXpWioJSIXAhcCtLW11btPhtFrxg5r4vIFh6bdDSPjRMnynDRZmVbaCgSTmLQC/ymspKorVbVDVTvGjLGIIsMwjHqRFePwFDBFRNpFZBBwFtCZcp8MwzD6LZmYVlLVvSJyKfAAXijrKlVdn3K3DMMw+i2ZMA4Aqno/cH/a/TAMwzCyM61kGIZhZAgzDoZhGEYRZhwMwzCMIsw4GIZhGEVIWvsm9xYR6QI2lvhXG/BGhCZGAO9muF5fkCOqDFHbM130rp7JUftnxt1emroYp6rheetVNZcH8HSZ8v9GfP/KjNfLvRxRZaiiPdOFyRGrHFE/M+720tRFud/OwqMvTivtjFjvvozX6wtyRJUhanumi97VMzlq/8y428u6LnI9rfS0qnZELc8bfUGOviADmBxZoy/IkaYMUT87z57DyirL80ZfkKMvyAAmR9boC3KkKUOkz86t52AYhmHUjzx7DoZhGEadyIVxEJFVIrJdRNYFymaIyOMi8oKI3Cciw135IBFZ7crXishxgffMduWviMjNkuDuGjHK8JCIbBSR59wxNkEZJojIP0Rkg4isF5HLXPkoEXlQRDa5108E3nOlu94bReSkQHmauohTjtzoQ0RaXP33RGRFQVu50UeIHKnoowYZ5ovIGnfN14jI8YG2UtPFPkQN40rzAI4BDgfWBcqeAo515+cD17jzS4DV7nwssAYY4P5+EpiLt7nQn4GTcyjDQ0BHSno4ADjcnQ8DXgamAt8DrnDlVwDXu/OpwFpgMNAOvAo0ZEAXccqRJ300A0cBFwMrCtrKkz4qyZGKPmqQYRYw3p1PA97Mgi6CRy48B1V9BHi7oPgQ4BF3/iBwujufCvzNvW87XshYh4gcAAxX1cfV08AdwKJ6990nDhkS6GZFVHWbqj7jzruADXhbvC4EbnfVbqfnui4E7lLVPar6OvAKMCcDuohFjqT6W45q5VDV3ar6KPBBsJ286aOcHGlSgwzPqqq/odl6oElEBqetiyC5MA5lWAd83p2fSc9OcmuBhSIyUETagdnufwfi7Tjns9WVpUm1Mvisdi7z1Wm5nCIyCW/08wSwv6puA+9LguftQOm9wQ8kQ7ropRw+edFHOfKmjzBS1UcNMpwOPKuqe8iQLvJsHM4HLhGRNXhu3IeufBXeBX0auBH4F7CXiPtUJ0y1MgCcq6rTgaPdcV6iPQZEZCjwW2CZqu6qVLVEmVYoT5QY5IB86aNsEyXKsqyPSqSqj2plEJFPAdcDF/lFJaql8juVW+Ogqi+p6omqOhv4Fd48MKq6V1W/pqozVXUhMBLYhPdj2xpoouQ+1UlSgwyo6pvutQv4JQlPb4hII97N/wtVvdcVv+XcYX+KYrsrL7c3eOq6iEmOvOmjHHnTR1nS1Ee1MohIK/A7YLGqvuqKU9eFT26Ngx+FICIDgG8Bt7q/h4hIszufD+xV1RedS9clIkc6V3Mx8Id0eu9RrQxummm0K28ETsWbmkqqvwL8FNigqjcE/tUJLHHnS+i5rp3AWW4utR2YAjyZti7ikiOH+ihJDvVRrp3U9FGtDCIyEvgTcKWqPuZXTlsX+5DGU/BqD7xR9TbgIzzLegFwGV5EwMvAcnoW9E3Cy9a6AfgrMDHQTgfezfIqsMJ/T15kwIvSWAM8j/cQ6yZc1ExCMhyF5+I+Dzznjs8CLXgP0De511GB93zTXe+NBKIuUtZFLHLkVB+b8QIj3nP34dSc6qNIjjT1Ua0MeIPB3YG6zwFj09ZF8LAV0oZhGEYRuZ1WMgzDMOqHGQfDMAyjCDMOhmEYRhFmHAzDMIwizDgYhmEYRZhxMIw6ICIXi8jiKupPkkDGXsNIm4Fpd8Aw+hoiMlBVb027H4bRG8w4GEYJXPK0v+AlT5uFt1BxMfBJ4AZgKLAD+JKqbhORh/ByYM0DOkVkGPCeqv5ARGbirX4fgrew6XxVfUdEZuPl0XofeDQ56QwjHJtWMozyHAKsVNXDgF14+2zcApyhXj6sVcC1gfojVfVYVf1hQTt3AJe7dl4Avu3KVwNfVdW59RTCMGrBPAfDKM8W7cl783PgKryNWR50maAb8FKi+Py6sAERGYFnNB52RbcD95QovxM4OX4RDKM2zDgYRnkKc8t0AesrjPR3V9G2lGjfMDKDTSsZRnnaRMQ3BGcD/wbG+GUi0ujy8ZdFVd8F3hGRo13RecDDqroTeFdEjnLl58bffcOoHfMcDKM8G4AlIvITvKyatwAPADe7aaGBeJsxrQ9pZwlwq4gMAV4DlrrypcAqEXnftWsYmcGyshpGCVy00h9VdVrKXTGMVLBpJcMwDKMI8xwMwzCMIsxzMAzDMIow42AYhmEUYcbBMAzDKMKMg2EYhlGEGQfDMAyjCDMOhmEYRhH/B0sdApIXoyULAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'].plot()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmcI2d95/956tLdrVYfM33N6fExY3zgwTb4wMGAnbDEsDHEJCzsj8MsSxLYTUg49pdr8W4SsuG1yQ9ISCAcCXE4Y05zmMPB2MZjezyXPfbcfU5fUusuVame3x9VT6kkldRSd0uqlp7369Wv7i6p1CW19Hyfz/cklFJwOBwOh+NE6PQFcDgcDsd7cOPA4XA4nCq4ceBwOBxOFdw4cDgcDqcKbhw4HA6HUwU3DhwOh8OpghsHDofD4VTBjQOHw+FwquDGgcPhcDhVSJ2+gPUyNDREd+3a1enL4HA4nC3Fk08+uUQpHV7rflvWOOzatQuHDh3q9GVwOBzOloIQcr6R+3G3EofD4XCq4MaBw+FwOFVw48DhcDicKrhx4HA4HE4V3DhwOBwOpwpuHDgcDodTBTcOHA6Hw6mCGwdOV3NhOYufPr/Y6cvgcLYc3DhwuprPPHIWv/MvT3f6MjicLQc3DpyuJlvQkS3onb4MDmfLwY0Dp6tRdQNakUIrGp2+FA5nS8GNA6erKeimUchrxQ5fCYezteDGgdPVqLZx4MqBw2kGbhw4XY2qm4qBKwcOpzm4ceB0NaqlGHLcOHA4TcGNA6erUXnMgcNZF9w4cLoa5lbKFbhx4HCagRsHTlfDlAN3K3E4zbGmcSCETBJCfkwIeZYQcpwQ8l7reIwQ8gNCyAvW9wHHOR8khJwihJwkhNzhOH4dIeSoddtfE0KIddxHCPlX6/jjhJBdm/9UOb0IizlwtxKH0xyNKAcdwO9SSq8AcCOA9xBC9gP4AICHKKX7ADxk/Q7rtnsAHABwJ4BPEEJE67E+CeBeAPusrzut428HEKeUXgLgYwD+fBOeG4fjyFbiqawcTjOsaRwopXOU0qesn1MAngUwDuAuAJ+z7vY5AK+zfr4LwP2UUpVSehbAKQDXE0JGAfRRSh+llFIAn684hz3WVwDczlQFh7MRuFuJw1kfTcUcLHfPtQAeB7CNUjoHmAYEwIh1t3EAU47Tpq1j49bPlcfLzqGU6gBWAQw2c20cjhu2ceABaQ6nKRo2DoSQMICvAngfpTRZ764ux2id4/XOqbyGewkhhwghhxYXeRtmTn30ooGiYb6N8jo3DhxOMzRkHAghMkzD8M+U0q9Zhy9ariJY3xes49MAJh2nTwCYtY5PuBwvO4cQIgHoB7BSeR2U0k9RSg9SSg8ODw83cumcHqbgaLaX58qBw2mKRrKVCIBPA3iWUvpXjpu+AeCt1s9vBfCA4/g9VgbSbpiB519YrqcUIeRG6zHfUnEOe6y7AfzIiktwOOtGdQShecyBw2kOqYH73ATgPwE4Sgg5bB37EIA/A/AlQsjbAVwA8AYAoJQeJ4R8CcAJmJlO76GUsk/muwF8FkAAwHetL8A0Pl8ghJyCqRju2eDz4nDseAPAjQOH0yxrGgdK6c/gHhMAgNtrnHMfgPtcjh8CcKXL8Tws48LhbBaqI87AU1k5nObgFdKcrmW9yiGj8slxHA43DpyuxRlzaDQg/diZZVzzp9/HQjLfqsvicLYE3DhwuhanW6lR5XBmMQOtSLGQUlt1WRzOloAbB07XwtxKiiQ03FspkSsAKE+D5XB6EW4cOF0LUw7RgIxcgwHp1ZwGoDR7GgAWknnuZuL0HNw4eIy//N5J/O6Xnun0ZXQFLObQH5AbVg6r2Wrj8M4vPInf/TL/n3B6C24cPMbTU3E8ca6qOJyzDphrKBps3DgkKozD3GoOz0wlMLfKlQOnt+DGwWNk1KLt2uBsjJJyUBoOSNtuJcuwPPSs2RWGGQ0Op1fgxsFjZFQdybwGw+DdQzaKHXMIyg13ZU1YxoGd+9CzF83j2QJ4RxdOL8GNg8fIFoqgFEhVFGKdnE/hW0dma5zFcYNlK/UHZKi60ZDBTToC0tmCjkdOL8MvC9ANigxv3sfpIbhx8BhpyyisVrgxPvz1o/iDrxzpxCVtWZhxiAbkst/rkchaqay6gcMXEijoBm6/YhsAIJ4ptOhKORzvwY2Dh6CU2q0bnHGHE7NJHDofR6ZQLCvs4tRH1YogBIj4zRZia8UdtKJhqwNVN2z1dslwGACPO3B6C24cPEShaEC3XB9O4/BPj5+3f+YLVOOougFFFBBQzBHmaxkH52teKBq20tje7wdQKpDjcHoBbhw8REYtLV5socoVivi3p2cwFFYAACvctdEwqm7AJwnwy6ZxWCudtcw46IZ9/+19pnGIc8PM6SG4cfAQzm6gbKG6mMwjWyjitsvMEd3xLDcOjaLqRfhkEQHLOKyVseRUZQXdRTnw157TQ3Dj4CEyhWrjkMqbx3bGggC4W6kZVM1UDsyttJZySDqUg6obUK37b+tjxoG/9pzeoZFJcJw24XQrMf92SjUXpEnLOHDl0DhqsdytVCvmkNeKePpCoiym4FQOIZ+IsE/irz2np+DKwUM43UrJCuUwGQsA4OmUzWAqh5JbqdY0uK89NYM3/f1jePpCAgCgiIIdcyDE/D0alLly4PQU3Dh4CLeYQ9oyDkNhHwKyyIOiTWDGHNZWDudXMgCAH54wq6GHwoqdreSXRBBCMBBUeMyB01Nwt5KHYDn2Q2GfI+Zgfg/7JMRCCndtNMALF1OYT+btbCU75lAjID2bMJvqza7mEfFJCCiirRx8srl/igZlbpg5PQVXDh6CKYfxqL+kHKxjYb/EXRsWS2kV7/6nJ5HMu78Wn/jJabzv/sOWcRDhl8y3eb5GAeFsImf/3BeQ4ZNEKyBtKgcAiHLlwOkxuHHwECxbaXwgUJatpEgCfJKIgSBXDgBwZDqB7x6bx/GZpOvtyZyG5UwBSykVPklAyGcKZBa/qWQmnoNAzJ+jQRmKJEDVi8jrJeUwEJTtpnwcTi/AjYOHyKg6RIFgJOK3FUJK1dFntX+IBmUekAZQ0M0q8mzBfbFnbS9mEjn4ZBF+WcRAUMaMQyEwtKKBi6k8XrZ3CIDZpE+RSgFpp3JYzWko8m65nB6BGwcPkVGLCCoi+gIyUnkdRYMildcRtna+pnLont0rpRR/+MAxHJtZbeo8zZq1UKtLqjOw77NcSjtiQUytZKvuO7+aB6XAHQe2QZHMrCSfJNgBaTvmEJBBaXktBIfTzXDj4CEyqmkIWBfRVF5DOq8hbCmHgZCCZL60e/356SUspdWOXe9GSeZ1fP7R8/j8o+eaOo8Zh6zqrhzSLsZhooZxYPGGXUMhvP/Vl+Hu6ybKUlmZchgImf8T7lri9ArcOHiITEFHUBHRbxmH1ZyGtKoj4jN/Hwiau9fVnIaVTAFv/ofH8fmfn+vgFW8M1mH2kVPLTQ3SaUY5KJZxmBwIYiaRsw3rsZlVfOCrR3DBMhjj0QDeeesevOLybbZbqUw5BM3eVjzmw+kVuHHwEBm1iLBPKjMOqbxeUg7BUvO9R04twaCoOVL01EIa//ETj9ipsF6EjfGcSeRwfrl6V1+LQtGKOTSkHMyd/45YEFqRYj5ppq3++LkF3P/EFL785DQAYCwasM9RLLdS3iqiA1D2P+FwegFuHDxERtUR8knoD5Ybh4gjIA2YDeD+/YVFAEBadd89H5tZxVMXEphaqQ7CegU2pxkAHjm91PB5mtXWIu0SkNatRZ3B3Eqswpy5lljQ+hdnVzAYUuxCOXaOqhl2ER0A272k1qiy5nC6DW4cPESmUERQkTBgGYGVTAGpvIaIr1w5xLMaHn7eXEwzNXbPbOH18nAg50L781PLDZ9XijlUPzfWn2rcUgJscd9h9aayjYNDUTlVA1BSDs46B/Y4ToPG4XQz3Dh4CDMgLWJ7v7lYzSRyZszBbxqLWMg0Do+fWbbdI5ka6Zy65XppZDRmp2CGayjsw89OLa3ZUptRijlUP3fWqPDK8T4AJbfSWDQAgZSMQzKv27UN45XGQRStmENJOSii+V1do7Mrh9MtcOPgIbIFHUGfZGYsBWWcWkjDoLBjDrGQAlEg+IefnQUA7BkK1VQOmq0cvGwczGt7y0t3YjWn4W9/erqh80oxh9rK4UXj/QBgN92TRQGj/QFMxU03Wyqv48BYPy4ZCePqyWjZY5TqHKqVg5dfTw5nM+G9lTxEWi3VNIxHA3huLgWgNAM55JPwr/feiJ+dWoJPEnF4Ko5zS6VAbipvuptec9VoyTh4eKfLFtqbLhnC8xdT+NufnsYbXzJZtZOvpJ5ySNvKoR8fvfsqvGr/Nvu2yVjAzk5K5zX0B2Q88J6bIDAJYcHcSgal8FtGgSkQbhw4vQJXDh6BBVJDimkIJgYCOLWQBgDbYADAwV0xvO+Vl+Ldt+1FyCeVZeZ8+8gc3vPFp7CQzDtiDt5dzJjh8kkCPvQrV6BQNPDlQ1NrnscC0lkXNxQL0Id9Et5wcNJOQQXMdNZSzMEM9FcaBnY9AKAb1DYK7JiXYzgczmbClYMHODGbtFMkQz5zMRqPBu0Fvs+KOVQSUqSyFhIs7z+nFaHpWyHmYF6bXxYwFg0grEhI5tzdZE5s5eDiUss4GhVWMj4QwEJKhVY0yrLAKmHxBXZtgMM48GwlTo/AjYMH+ONvHseRaXPQDGsSNzFQcq24LXTsvs7pcWxXqxUNR8zBuztdZhwU0TSIAUVETlvbOLCYg6tbyWquxxSYk0EroJ/IakjlNYR97kaXxReAklEghFgN+bhx4PQG3K3kAeKZgp2bz4zDuNM4+NyNQ9gnolA0ULAWLPYYqu4wDh7e6TLDxRbjgCK6uooqqZfKytxsbqpgwDIOS2kVmUKxQeVQXv9Q4MaB0yNw49AmKKX4jb9/DB//8amq21ZzGiTL9x2yBtM4lUOtRYwZEuZKKSkHukViDua1sd15QBYbSmetH5C2lIOLQY1Z8QcWlK5pHCSHcqhQEV5WYhzOZsLdSm3iqQtx/Pz0MoKKWHXbak7DG18yiaGwDzfsGQQATESD9u2RGu4P5jrJFHQMhBR7sS3o5W4lw6D4+ell3HTJIAipDsB2Cma4WNA3qIg1R3k6Yc8trxkoGhSiI6icUXX4JAGyWL3viYUt42C16qgVy3EaB5bKyq7Ty8aWw9lM1lQOhJDPEEIWCCHHHMf+mBAyQwg5bH39iuO2DxJCThFCThJC7nAcv44QctS67a+JtUoRQnyEkH+1jj9OCNm1uU/RG3zxcTMLZ6ViHkNeK0LVDYz1+/HfX3Wp7ULqC0j2z/ViDkApt78s5uAISD92Zhlv/vTjeOpCYpOf1cZgLhq2GDfqVmLzHIDqmQ7OdOBKmHJgM6MbcStVKwduHDi9QSNupc8CuNPl+McopddYX98BAELIfgD3ADhgnfMJQgjben0SwL0A9llf7DHfDiBOKb0EwMcA/Pk6n4tnWc1q+NaRWQDVxoGNumSN3RiEEEwMBBBUxLKdsROW2cRcKa7KQTPsNtPnlzOb8XQ2DVUvQhaJ/fwCstSUWwmoTmdNW/2p3GBprazJXy2jW0s5KJLg6boRDmczWdM4UEofBrDS4OPdBeB+SqlKKT0L4BSA6wkhowD6KKWPUrM38+cBvM5xzuesn78C4HbiJd/HJvDNI7NQdQMHdw5gudI4WKmbfYFqF8d4NFBzdwuUlAPbPbNdbaFolPVWYgvudNxbTfjYjGdGoEm3ElCdzpqpoxwUSUDEL9nGIdKAW6lMOcjcrcTpHTYSkP4tQsgRy+00YB0bB+CsYpq2jo1bP1ceLzuHUqoDWAUw6PYHCSH3EkIOEUIOLS4ubuDS28uDx+axZyiEWy8dRiqvly1urL7BzTi8+caduPfWvTUf1445WAtk3lpYy2MOBrLW8RnPGYeiHYwGgKAs1hz96aTcOJQbE+fkPDdiIcUeF1rL8Poq4gyln3lAmtM7rNc4fBLAXgDXAJgD8H+s4247flrneL1zqg9S+ilK6UFK6cHh4eHmrrhDrGY1PHZmGa8+sN1unOecA83cSm7B0V+6fARvv3l3zcdmi2DajjmYi6ZZ52C+hAXdQN5SDm4zlDuJqhllxiGgNJatVChSOy5QmbGUKeg13UWA2dmWDfypbRyqi+DYcZ7KyukV1mUcKKUXKaVFSqkB4O8BXG/dNA1g0nHXCQCz1vEJl+Nl5xBCJAD9aNyN5Xkeeu4idIPizitLxmHFMU2MzSSujDk0Aos5VKayVmYrMVfNdLzxgTrtwJy0tg63km7YMy+qAtL52jEHoFQIBzSWrVStHLhx4PQG6zIOVgyB8XoALJPpGwDusTKQdsMMPP+CUjoHIEUIudGKJ7wFwAOOc95q/Xw3gB/RZmZGepwHj81je58fV433l4xDumQcVjdkHEqprECpCE5zFMapumEHbWcTeRiGd15aN7eSVqRlbiM3tKJhz9mudCulrWl6tWCFcJJAyv62k9pFcDzmwOkd1qxzIIT8C4DbAAwRQqYB/BGA2wgh18B0/5wD8C4AoJQeJ4R8CcAJADqA91BK2af33TAznwIAvmt9AcCnAXyBEHIKpmK4ZzOemFc4dD6OV12xDYJAbOPgDEon7ZhD8yUnPkmAKJAq5VBZIW3HIooGltIqRvr8639Cm4iqG2W79IBVA5LTiq51CgytaGAkYj6HSuXAZmLUgv0PIn6pZs0HL4LjcBowDpTSN7kc/nSd+98H4D6X44cAXOlyPA/gDWtdx1YlVyjaLhA75pAtVw5+WShzXzQKIQQhRXTUOTDlQO2Yg6oXyxbQqXjOO8bBJeYAmK9ZLZcPYD4/9po6lYNeNJDTijV7JgFO41D7PjWL4GTB0+1IOJzNhLfPaDFa0YAsmjvUaEAGIcBy2qkc9LoL4VqEHW27a2Ur5TTDnnrmpaC06VYqLb6senytQrhC0bDdcG5daUP1lEOwpBxqodjN9mD/7wDuVuL0Ftw4tBDDoNANartIJFFAf0AuK4RbzWnrijcwgj6pqs5BK6tzMJAr6PYMZS+ls5p1Dg7lIJsL9loZS1rRQEgRoYiCnakFlIoBG4k51DMO7Jr8kljmeuJuJU4vwY1DC2ELtNNNEQsp5dlK+Y0ZB3Pgj+VW0kpFcJXZSoNhH6JB2VMZSwXdKPPpl2IO9WsdNN2ALAoI+kTMrebwP/7tKOKZAlaztWtGGA25lSxj7rw2oDQ+tIvyJTicmvDGey2ELdDO7JdYUKnKVtq2gRhA2Ccio+qglCLvTGVlvZU0A7lCESGfhPFowGNuJWNdbiWtSCFLAkKKhG8+MwuDAjfsHrQVQ73X0zYOddQFIQSKKJTFGwBTORgUlhrsqiJ+DqcKrhxaCEsndWbexELKprqVQoqEjKpDK1KwDW25cjBTWf2yiO19fiwk1XX/rc2mMpU1IJcC0rWg1GxHLosCgooIlpk7t5rD3GoeADDaX8c4NBBzAEyVUKkc+BxpTi/BjUMLYRlDTrfSYFipKoLrW2OhqkfIJyFT0Mt84Zpe3lsprxURkEUMhhUspb1kHGpkK9UphNMta6CIxK7zkEWC2UQec6s5CAQYifhqnt8XkBALKZgYCNa8D2D+z6qUg8xGhfK4A6f74W6lFsJ2707lMBBUEM8UQKm500+p+gZjDmYqa96RYummHIKKaKsWw6AQanR6bSeqVl4h3Yhbyfma3rhnEFeO9+HQuThmEjn0B2SMRPyQ6tRIEELwg/92a92YA2C6AquVg2UcukA55LUi/uLBk3jfq/ZtKFuuF6GUemouSqvgyqGF2DOSKwLSukGRzOtIqToorR9AXYuQlcpaphys3kqEAJSaQW+/LGIo7INuULsqu5NQSptyK+UKRXz50FSZq+4Dv3w5PvK6F2EsGsBsIoe51RxGo2vHbwbDvrL/iRuuyqGL3ErPTCXwmUfO4hdnuqZTTVv41pFZXPOnP8CJ2WSnL6XlcOPQQkoB6dIuYyhsujx+/yvP4NHTSwA2aBwUCQXdKCsGYxPSwlbX1rxmIKiYbiUAWM503rWkGxQGRcNupa88OYX3f+WI/aGUHeeN9vst45CvG29oBp9rzIEph63vVrJbrnTBc2kX3zk6h9/5l6exmtNwajHd6ctpOdw4tBC3gPTtV4zgbTftxuNnV/Bf//kpAOvrq8RggVVnLIG103B2Jw3IIoYtw7SYKp8p0QncVJUimu1A3Np2P3EuDqDUesRpcMeiAcSzGqZXchjtD1Sdux4mBgKYjJXHJdi1dkNnVpb+nOcV3w3zdz89bb+/EtnOf4ZaDTcOLURzqXOI+GX84Wv348H33orxAfONthGfLzt3IZW3j7FdobMYLKCIGLSMgxeUAwvqOlNZCSEIyCJyheoF69A50/3BPpROgzseNV/HQtHYNOXwqbccxJ/+6oGyY93kVqqcAcJZm3PLWdyybwgAkMh23jXbarhxaCEFl4A0Y3u/H198x41460t34urJ/nX/DaY6FlPmgi8KBFlrV1imHBQRQ5ZbaSnlAeNgLbCVnVHNtt3lymEmkcOslaYatz6Uztd0LFpSC9s3yTjIolAV2C5lK3Hj0GsksgWs5jTsHQ4j7JN6wjjwbKUWUnBxnTiZjAXxJ3dV9SJsCtaAjtUvRPySaxuJgCxiIKhAIMBSuvOS2DYOFX79oCJWZSsx1QCUdmzlxqFkEDbLreRGN8UcUnluHJqBjZbdORhENChztxJnY9h1DnVSKzcKUw4XUyXjwBZXZ6FXUBGttuE+b7iV9Gq3EgDLrVRpHOIQrdRb9qFUpFLMYVuf324suFluJTe606209Z9LOzi/woxDyDQOHsj4azXcOLQQtzqHzcaOOSRNt0vEJ7vGHNjQmqGw4omAdKGuW6nCOJyP47qd5pjyuEvMQRYFjET8axbAbZRuUg6lAVFb/7m0gwvLGQDAjljQrFXiyoGzEdZyK20GdszBylYK+yW7jYZzrkHQSmsdCndeOTxxbsV2a7gpB6dbqaAbeOFiCtftHIAsEteYA2C6ltYqgNso7P/YDTEHO1upCwxdOzi3nMVIxIeAIqI/INtNHrsZHnNoIaWAdOuqKf2yAEUUsGjFHJytOCpTWQFTOZy/kGnZ9azFUlrFG//uUdx5YDsA95iDM9h3bjkD3aC4bFsEIZ/kmq0EAHdeuR0XW9w3iimHwhpjTLcC3K3UHBeWs9g5aKY294py4MahhdjKoYW7WUII+gKyXefgbAsRqUhlBUzlsNRBt1Iyp4FS4Mj0KgA3t5JU5lY6OZ8CAFy6LYKQItnKofI1vffWva28bACwW310h3LgbqVmOL+SwS37hgEA0aCM1ZzmmTY0rYK7lVqIW51DK3DOn47UUg6WcRgM+5DTiq6FZu2A7VRZ6/Bqt5JQdm3PX0xBFAj2DIcQ8olI5i23ktT+D2VXxRy4cmiYXKGIi0kVO62iyP6ADIOWMr66FW4cWkg7AtJAKe5ASCm2AFSnsgJw1Dp0Rj1U+rgrlUNQkcqylU7Op7BrMAi/LCLkK8VTWv2auiEJBALprmylbjB0reaClam0w+FWAoBErrtdS9w4tJB2BKSBknHwSULZ3wpXpLICwJCVzbPUoaB0pRuj8rWpzFZ6YSGNy7ZHAJQbu1a66mpBCOmaOdKl9hncONRDLxr46PdOghDgwJhZrBq1aoviXR6U5sahhRSsOgepxX5JZhz8sljWc4jFHAgp7dDZriee6ZBy0Oorh4AsQitSaEUDea2Ic8sZXLrNNA7MwAGdUQ6AGUDvhnkOTDnUm53RC0ytZPG/v/ssvnxoynZZOvno907ih89exJ/86gFcMhIGAESZcmhhUDqZ13Dv5w/ZKeqdgBuHFqIVDSii0PLe76zWoVI5sGE4AVm0ryHU4CjOVlHp43bOcwDKZzqcWkiDUuAyyziEHMqhU2M6FVHY8spBLxq2Uej1mMPXnprB3/30DN7/lSP4q+8/X3X7j55bwK2XDuMtL91lH2PKoZUtNE7MJvH9Exfx5Pl4y/7GWvBspRZS0I2Wu5QAp1tJLNtRB2QRkkDseAMABK0Ftt4ozlaylnJgBiBb0PHCgpmptG+buWNzupXkNryubvhkYct3Zc04/ve97laaTeQwFPZhKKzYsQUnGVWvKqwcaINyYMFuNzXTLrhxaCFa0WjLDrfkVipXDpJI4JMEO1MJAIKWoch0OFvJJwnQikaVy40ph4xaxErG/GAMh/3WbZ2NOQDoipgDcymJAul55TC7msN41I/hiA/T8VzV7SlVL9uUAKVaola20Eir5mN3MiOKu5VaiGkcWv8Ss1TWSuVgjroUK5RDp91K5t+9fHsEPkmscrmxD2JG1e1FLGRdc9jngZiDJGz5DB/2usZCSlfETzaCOSAqgNH+AOZWy/37lFJkXIyDJAro87e2M2uaKYcO9nDixmETOb2Yxr8+ccH+XW27W0koc9PIolClHBRRgFRjoE47YL7uG/cOuo70ZOogU9CRVnX4pFLrbOZyEgjsRnztxjQOW3u3zQrgBkNKT7fPoJRiLmGOlh2N+rGa08o+FzmtCIOWZ/0xokHF1a30hcfO441/++iGry1pu5W4cugK/sfXj+EPvnoUU5bvUivStrg/2JhRnyyUN6SzjIVTORBCEFTEsrGi7UTViiAE+L1XX4Zv/fbNVbeXlEMR6YpdGzMOnVINgOVW2uKuGPa/Hwr7oBUpigbt8BV1hmReR6ZQxFh/wO7mO5soqYe0rVzdjIPsmsr69IU4nji/Ytc4rRf2tzsZc+DGYR0UdKNq5314KoFHzywDAL59dM66X7GtysEviWXGSBbNvHyncgDM3XnHKqR1Az7JNGLOGAKj5PYy3UrOD2bIun+n4g2AWZex1d1KbOFhBZG9GpSeWzVjDKNRvz0HhB0DSq6diKtxUMpG8zLimQIoBVasVPFsQceF5WzTBjhlGYVkjiuHLcWfP/gcfv3vHis79rc/OY0+v4TLt0fwHcs4aEXanpiD36EcnG4lQcA7b92D/3TjzrL7B33VA3XaRV4r2u3D3WAGIKMWq42DZTg6lakEAINhBTOJHIwtvNtmMQc2NrZnjYOlEkb7AxhjxsGhHJhFYiV3AAAgAElEQVTCclMON+yO4fhsEk84BlEBpcK4xZSKLx+awv4//B5u/eiP8emfnWnq2tIeyFbixmEdnFvK4Pjsqp3SmCsU8b0T83jT9TvwumvHcWR6FVMr2fZlKwVLqaxsVy0JBIJAcPd1E7j9im1l9w8pUseMQ65QLHNzVcIMQEbVLbdS6b5h263UuWZnL9s7hKV0ASfmkh27ho3CMtUGmXLY4jGU9TJrqYSxqB/b+n1lxwAgZWUMVQakAeD/uWkXtvf58ZFvP1u2UWDdWpfSKp66kEDEL5ldk5sczcvUHc9W2mIk8xoMWuq5ksybnUZ3DobwmheNAgC+f+Ji2wLSYUWCQFgRnLlw1lMsAUW0d4/tJq8bdZWDMyCdUYsVyqHzMYdbrQHzD7+w2LFr2Ci2WylkLoidqnnpNHOJPESBYCTih08SMRT2uSoHN+MQVCT83h2X4ZmpBP791JJ9nLmTltIFzK/msHMwCL8s2FMhGyXJs5W2JqvWP+zckjkXgVn3sF/CZCyIoCJiLpFrWyqrIBDsGgphLBqAIlqulzq765DLtLV2kdeKVYVvTkSraC9bcHEreSDmMNLnx+XbI3j4+S1sHPI6RIHYiQy96laaXc1hW8RnZ76NRf1lyoHVGrhlKwHAyy81W3ift6bEaUXDXgsWUyrmVvPY3heAso7aGO5W2qIw43DWMg5sJxa28/ElpFXdbp/RDr75Wzfj3bfttVtZ11MsQZ/UOeWwRswBMF1LaeZWUlxiDh00DoC5KDx5Pt6x13CjZFQdIUWEX+6eFuTrYTaRw2g0YP8+2u/H/KozW4nFHNzfrwNBGYQAy2lTLTjrHpbSqlVD4Ycikqar6tmaklb1jsW3uHFYByyD4Ky1Y2BWno3lDPslpFS9be0zANPlIouCbYzqLaBBubMB6XoxB8DKplJdspWYW6kDsxyc3HrpMLQixS/Orqx9Zw+SVosI+yTbSPdqlTRbvBmVhXClbCW56lzALIYbCCr22F3ndLiplSxWcxpGo34oVjeAZmDZSpQC6Q5lFnLj0CQFvdS07OxipXIwF6+IT0I6r7ctW8mJ3IBxCPk6F5DOa4a9Y61FyFJemUKxLCDtkwSIAumoWwkA9lndOafj1b14tgLM6AZs49CbymG+wjiMRf1Iq7rtysmoOgSCuu/XWEixlcOKo9PxsRlz0uFov2kcmlEOlFKkVd1ONe5U3IEbhyZhLiVCzPnGQMk4sClsYb+5uLVTOTCYP79ezCGoiB2cBNeAW0kR7ewOp3IghCCkiB13K8VCCggBFtNbc9hLpmAah15WDgXdgKobdo0QYBYFAsCK9X9lRZj1uioPOowDq5iejAUwaymQ7X0ByGJzykHVDWhFijHL5dWpjKU1P2WEkM8QQhYIIcccx2KEkB8QQl6wvg84bvsgIeQUIeQkIeQOx/HrCCFHrdv+mlivOCHERwj5V+v444SQXZv7FDcXtqvYOxzG3GoeuUIRaesYW8jClnIotCkg7aQR5RBUzJkJnegumtcbiTlIWHAxDoD52rbb4FYiiQIGQ0rT6YleYTGlYjCk2DviXlQObHPkLMRk89fZZ7yyQt+NobDPHpzFGkVeOhKxbx+z3EqFJowDMwas9sLLyuGzAO6sOPYBAA9RSvcBeMj6HYSQ/QDuAXDAOucThBC2EnwSwL0A9llf7DHfDiBOKb0EwMcA/Pl6n0w7YMrhmskoAFM9sBbIpQZxMlJ5zQpIt9c/zhbOugFppXNtu3OF+qmsgPk6soW38sMZ8kktH57UCENhn2uF7FZgPpnHaNRv/x9SeQ1//I3jHR0s024qP7NAqdsqW5zTeb1mphJjMKzY7iQWc9i3rWQctvX5IYvNuZWYJ8LzyoFS+jCAysjbXQA+Z/38OQCvcxy/n1KqUkrPAjgF4HpCyCiAPkrpo5RSCuDzFeewx/oKgNtJPR3XYZgVv9oyDmeXMkjldSiSAJ9kvtEiHQhIMxqLOXSubbeqFdeOOSgSdCtDo1I5vOOW3fj1l+xo2fU1ynDEtyWVQ65QRCKrYbQ/AL/1fn3iXByf/fk5/OTk1k3PbZasWkc5WJ9x5n6rRyykIJE1N4LxTAFBRcT4QMC+zS+L8DWtHMy/P2Y1puxUOut65zlso5TOAQCldI4QMmIdHwfg7CsxbR3TrJ8rj7NzpqzH0gkhqwAGASzBgzDlsH/U3B3MreaRVrWyHW7EijmIhHTArUTKvrsRUEoDddpNo26l0s/l9/WCYQBM5XDGSkjYSrA8/tF+P3yWkX5u3qz2Xmnh8Bqv4aYcIhXKIZXX7WO1YC1I4pkCVrIFDAQVDFuBZBbsbjbmwLKkxqPedys1g9uKROscr3dO9YMTci8h5BAh5NDiYmd2OewfNTEQhCgQrGRUZKzUQEbYJ4FSQDdo25UDIWY2T13l0KFRoXrRDLSxHWstnLOi1/L5dorhiOlWMoVw58lrRZycT615P5bHP9ofgE8SrMQKM+uqU3PFO4FbzIEVBTqzldaMOYRMQ7CULiCeKSAWUjBsTY5jxkFp0q3EqqPd3EoZVcf19/0QX3962vXczWS9K9dFy1UE6/uCdXwawKTjfhMAZq3jEy7Hy84hhEgA+lHtxgIAUEo/RSk9SCk9ODw8vM5L3xjsH9cfkDEQVLCS0ZDKl7+JnH7KTmTWKJJQN90z6Ghu105YD5+AsnYqq9vPXmIorEDVDaQ8Ugh3/y8u4LV/87M140iziVI/IULMSYGsY+hyLxkHVuCmlG/qgNJnvJGANFMOyxkVK1kN0aBsZz1tZ8pBaq59Bos5DAQVBGSxzK303HwSCynVrqlqJetdub4B4K3Wz28F8IDj+D1WBtJumIHnX1guqBQh5EYrnvCWinPYY90N4EfUK9sxF1ZzGnySAL8sIhaSsZJRq9xKzp87kZMvi/XdWSFHW+x2wrJiGkllZbi1S/YCbHe45JG4w3Q8h0LRWNM/zYq8tvWZC5ezILGXlAOLtwUdbiVRIIj4JNvnn1YbizkAZo1DImsqB7NXk4BdgyEAzSsHlv0Y8UvoC0hlbbtPzJouwANjfQ0/3npZ85NHCPkXALcBGCKETAP4IwB/BuBLhJC3A7gA4A0AQCk9Tgj5EoATAHQA76GUsq3Mu2FmPgUAfNf6AoBPA/gCIeQUTMVwz6Y8sxaRzGl2bnQsZGYq5LQiRiKlYhqnn7ITaZeKJNRtax3skFvJNg5ruJW2hnIwjcNiSsWe4XCHr6a060+rOrbVud/cat5KYzX/B+Z3czHqpZgDe++HKmaKRPzmYsxGhK4Vc2CFakvpAlYyZswhoIj4zntvsWMGikTWlcoa9kuI+GW7OywAHJ9NIhqUy4r3WsWanzxK6Ztq3HR7jfvfB+A+l+OHAFzpcjwPy7hsBVZzmu2bHAz58Nx80hwlOORUDiXJ1wm30vb+ALZZO1s3gh0KSLNiK18DFdKAWWgYVOobkk5hKwePFMKxtNq1+j3NrebKxrP6e1U5qNXKATAzllJ5zR4RutbmpM8vQxIILibzSOV1DARNY7HXsWFoWjmoOvzWVMc+f4VymEviwFhf3cK8zYJXSDfJqkM5DIRkrGQKZszB7+5W6sTsgX9+xw34vTsuq3k72y3Fsxr+5JvHsZBqT347Uw5r91Yybw8p9atTO8mwrRy8URvAqnTXiiOZLSNKzeZYRf32Pn9vxRws5RCseC/2BSSk8rqjX1p94yAIBLGQgkPW0B82I8NJs+0zUqpubzD7ArLtKtSKBp6bT2H/aOtdSgA3Dk2TzDvdSj4kchqS+epUVkYn3ErOpmpusLGhPz25iH985Bwefr49WcONxhzYa1mrG6YXGAgqEAXiGeXAmr/VUg6//S9P4wNfPWJ2Iu2vVg7X7ogildc3PPt4q5Ap6PBJAqQKZR/xm4txZb+0egyGfXjqQgJ+WcCrD1Q79ZpNZXWm0EYDst3t9cxiBgXdwIGx/oYfayNw49AkqznNrqSMBWVQavZp8VJAei0USYAsEjx1IQ6gVHTTaphbaS3jwNxeXo03AOaO0SstNAyDlpSDi6uwoBv4/vF53P/EFJJ5vUw5sIJEVtQZ75G4Q7ZikBSjz28ph2aMgxWU/s0bdpbFHhmKJEA3aMOtt9N5zTYOAyHFdvcdnzWb+e1vQzAa4MahaZI5vaQcwiW/frhGELXTfYBqEVQkewBJu8rzS8qh/mvCXkuv1jgwhsI+LHqghUYyr9kV5WkX5fDcfLJs2MxYRcxhKOzDhFXVG890brhMO8kUdNd4Fos5sNexkQ3KSJ8PPknAu16+x/V2FndsNCidUYv2tcWCClKqDlUv4sRsEj5JwJ6hUEOPs1G8uXJ5FMOgZW4ltmMAKtSCJDi6o3rzJXZ+MNqlHHKNxhx8pZiDl2GFcJ3G6drKusQcDk8lAADvutVcvCZjQfu2V+3fhjddP1mWktkLZNWi6/sr4peQdMQc1spWAoD/9spL8cV33uiqGoBSXKdR45DXSzNPYlYMI5HVMBXPYkcsWOUKaxXe/vR5iKJBsZxWQWmpkpJlJgDVowQjfglqurAljIMzG6KVNF7n4H23EgCMRHx4di7Z6cvAssNAuSmHpy8kMBzx4Q/uvByv2r8N11ouJMB0hQCOFho9YhwyBb0qUwkwP9tFg9pTHt0CzJVMxoJlBrcStgZoDQalc4UiAgMl5QCY/5fFlIqRvtpZiJuNN1cuD/Lhrx/FS//sRwBKxsH5xql0gbDfvepWci68zjzqjbCYUvFfvvCk3deeQSnFA4dn7L5Ua6Wy+mUBAimNXfUqEwNBLKTUjre8dmYZuQWkD08lcM1kFIJAcHBXzDUDzF6EeiXmUKitHADgyPQqArKIbTXUQDMoTSqHnGPmyYBD0S2mVTtLrh14c+XyIOeXs+jzS7hqoh8v3mHuvJzKoWrugPUm82JAGii5dnYPhTZNORyeSuDB4/N49PRy2fHn5lN47/2H8Y+PnAOwtnIghCDil+0umV5lMmb66WcSuTXu2VqYcpAEUhWQjmcKOLuUwbU7om6n2kSt93Kv1Dpk1NoxBwA4MpPArqEQhE1oD19SDo0FpJ2jdJnrejlTwEJStetr2oE3Vy4PklZ1XDMZxTd+62ZcYg3zUCTBbu9Q6ZtkyqHT845rEfHLmBgIYEcsuGkxBxZTOHmxvAEca/bGFtG1KqQB4K/eeDXedvPuTbmuVsFcCRdWOjsudDFdACHAaNRfVefw9JSZkXbNZH3joEgCIn6pZ9xK2ULtbCUAmFrJYfdQbVdRM5SUQ2MKM1coGQemHKZWslB1o63GwdtOXQ+RVnXscskSiIUVq2il0jiYOxCvKoffu+NSZFQd//jIOUxt0uKWtwqLnq8wDhcdQ2REgTRUGHj7FfWaQHiDHZZxmO6wcVhOqxgIKujzy1VupR8+u4CgIuLFOwZqnF2CtYPpBbJ1spUYuzcpK4gN/Co0oBwopcjrpYFY0YAMQkz1DaBm0LsVeHPl8iCVnVcZzLVUuQthSsKrAenLt/fhup2xsgrMjcLacVS2jr6YLAVM/ZLg2arnZhkO+6BIAqbinXYrFTAYUhDySWUB6aJB8f3jF/FLl42s6coDzPdyr9Q5ZGrUOfQHSsd2D21Oz6xmYg5akaJoULtQVRIF9AdkPG99prhbyYOkVc01rY35BGsFpH0eDUgzWOreZpCzitzOLWfLgrQLqTwifgl+WWhokdoqCALBxEBg05TXelnOqBgMKwj7pLKYw9MX4lhKq65Vu270inIoGhQ5rdiActgkt5Jo/p1GqqRzLhl9saCC04tpANw4eA6taCCvGa7KIRZSEFREiBWBq7DHlQOjzy+joBubknHD3thFg5ZNSbuYVDEeDeCmvUN2ple3sCMW7HjMYTldwGDYh6AiltU5PHhsHooo4BWXj9Q5u8RIxIe51bxnBhi1CvY+rZetBGyecpBtt9LaxsGt/1gspNhFju3MVuIxhwbI1Cmlv+uacXtik5M+awfi1VRWhnOo+kZ39TnHrvX5iym7zH8hlcdInx8fef2VXbcznRwI4qnz8Y5ew1JaxVDIHD7kdCv99PlF3Lh3sOGsrwNjfbj/iSnMJHKYGNicXbMXydboyAqYi7IkEIR8EgaCm7ORacatxIY1ObsIsKC0LBK7ALcdcOPQAHafFRe30s37hnDzvqGq43dfN4GxqN/zhVxs4UjltQ1L1pxWRMQvIa8V7QAaYAakL9sWwWh/oKyvTzcwGQsgmdfLuvW2E61oIJnXEQv5kMxrZQHp+dU8brqk+r1Zi6smzIymo9OrXW0cMjVmOQAsjVrCjsHQpsXG7PYZjSgH3UU5WHHNobBvU1JrG8Xb21qP0EwTLsZwxIe7rhlv1SVtGn2B8tGIGyFXMNDnl7FnKIxnrJYNRYNiKV2wJ491G5PWItqpuAMLIMfCCkKKiEyhCMOgUPUiUqput8VohMtHI5BFgmemV1t1uZ7AnuVQY1bIjsEQrp7YvM6ndvuMBoyDrRwc18ZaaLQz3gBw49AQjfZ234o4lcNGyWk6AoqI1149ikfPLOPh5xexnFFRNGhby/7bCat16JhxsBrlDQRlW6XmtKLd5rkZ4+CTRFy+vQ9HphObf6Eewp4CV+Pz/MV33IAPv+aKTft7dhFcEwFpN+Uwwo2D90jVcSttdSKOmMNGyRXMDJB33LIHu4dC+KNvHMfUipnm2c787HbCFFEru7Pe/4sL+ODXjrrexmI4saBiL3YZVbdbeA82YRwA4KqJfhydXm24vfRWxJ4fXUM5hHwSfA0UajaKUkc5sHbqLAnArf8Yizlw5eBB7A6NXagcWOA8mdsM5WD2hPHLIv7otftxdimDv3/4DABgW5cqB7bA5Fo4j/uR08v4ztE519uYW2kgpNjDkdKqXjIaTRqHqyeiSKk6zi1n1r7zFoVldLUrHlhPOTxweAb3fuFJPDtnxujYzBOncmAGvp2ZSgA3Dg1RLyC91dls5cDe1C+/dBi7h0L43ol5AOjamAPb4eVa2HxP1YpYzWmui4vTCLAAa0Yt2pPhGukq6uRFlq/96Ez3xh3WUg6bDVMOqotyOGLFd6bipluSbTICXDlsDbo55hBSJAgEm1Il7SwsIoTg7usmwFLm2/3GbheiQOCThJYqh7y1qLhVL7MOuNGgbL8/MwXdbqDnbA7ZCGNWNplXxp+2ApbK2q55IYqtHKpddces6W6zVt8xuwhOKS3Ne4ZDODBmdjRoJ9w4NEC6zW+mdiIIBGGftDnKwdFNEgBef+04CDFlsdeLATdCQBFbqhyYH9qtRmQlY84v90liWcxhJWM244s2aRxY24ZsjVnU3QBLZXWrc2gFtWIOetGw54Ew4+AWc+jzy/j279zStvGgjO5b7VpAWtURUsS25hi3EzZUfaPkCsWyFLyxaAC3XTrsOoCmmwjKop0B0wqYO2LFZTcfzxYwEDLjRs6Yw3KmgIGgUlW5vxaKJEARBXsB7UZSeR2ySNrWFFMUCESBVLkFTy9m7BgD61js5lbqFNw4VHDft09gKOzDu16+1z6WzutdGW9g9AXkTZnpkCsUEax4U//Nb7wYRRc53U0EFLGlbiXV2k0uuSqHgp3qWFIORfN4k8FoRtAn2k0Uu5GVjIpYSGlrA0hZJFUV0iyuM9bvx0zC7Fyc14uQBOIJpd35K/AY3z9xEd89Nl92LO3SkrubiPilDdc5UGo2MwtUBPnCPgn9m9SGwKu02q1UUg7V6bLxbMF2HTHjkC2YymG9xiGkSFVzIbqJlYyGWKi9MTBFFKrcSsdmVhFURNx0yVAp5lAwPKEaAG4cqkhktaqCppSqI+zxqWQboc8vb7hCWtUNGHTtKW/dSFCWWrrTdsYcFlMqnrpQ6uUUz5aMAFNtadUMSMeajDcwQl2oHL70xBTe9tknADDl0N7PsyIJVcrh2Mwq9o/2YTIWxGJKhaoXzXTwNmVRrQU3Dg4MgyKZ17CcKZT1qEnnta6scWD0bYJyYAtYu9IDvYRfEe125a2AvbbLmQL+70PP4+5P/hw/PrkAwKyQZhlJkijALwt2QDrWZBorI6hIXRcneuT0En58cgF60bBcbu1XDlqFcphJ5LB7KGQ37pxfzSOvFcua7nUSb1yFR0jldTv1kuUdA93vVgr5pA0HVN3K/nuFoCyWdaTdbGy3UqaA5y+mYVDgt7/4NJ6dSyKt6mW74LH+AJ6/mEY8W2i6OpphKofucistplRQaqbormTW/9qsF9lFOSRzGvoCMsaiZg3QTDxXVivUabhxcLDqqBK+sOwwDl0ekA76xKrxks3CFpPKmEMvsBkxh6Ua7TcopWXK4cxiBrdeOgzdMPDXD70AoFQkBQC3XjqMn51agkGbr45mBBVpw+8Hr7GQMl/fmUQWybzedP3HRqmMORQNikzB7GI8bimHmUQOeZ0bB0+SyJWyQZwDXLpeOSgSVN2A3kBjsFp4KQWv3Ww0W+mZqQRect8Pq2ZvA4BuULA2R+eXM1hKq3jZ3kHcsHsQDz1nupacsYXbLhtG0TphvcYhvAlK0mssWHPMT86bE9XW63JbL7IolKWyssLaPr+M7f2mcphN5M10cI98hrhxcOBUDiwoTSntfuNgV9auf0GwJ1j1onLYYJ3DueUMKIWrcXBO6GOzuHcPhXDLviF7J+pUDjfuGbR91utXDt0VkM5rRTvh4rl5s+is3W4lRRLK2mewuqKI3yxgHIn4MJvIWTEHb3yGuHFwwNoc+xxD43NaEQbtzr5KjJC1oG/ElZDtYeUQtNxK6x2vySqf56xcdyesSMo5lWzvcAi3Xjps/+40An5ZxE17h6qON0PI110B6cVUyWX3nNXgrhNuJadyKBkH8/86Fg1gOpGt6jLQSbhxcMCUw/6xPtut1M19lRjO/Pj1kuth5eCXRVDq3litEeLWpmR2NVd1m2pNBmMZLaJAsCMWwr6RsN3pNlpRR3LnlduhSILdJ6lZgoqIvGbY7qmtzkKqZHSfZcqhzW4lRSqPObB2NWzY1mQsiKmVnGutUKfgxsEBMw4vGu/H1EoWhkHtWQ6RblYOdtuFTXAreWTX00422rY73oByYMZhciAARRJACMEt+4ZBSPUu+O7rJvDIH7yizN3UDKyHWLe4lhYsd5wiCvai3HblIAlljfdYi3zWMn9yIIDZRA5ZlbuVPMlqToNfFrBvJAxVN7CYVu03Uzc23WME2WKwGW4lj+x62gkzDtl1ZiytWJ1VZ1dzOLeUwYE/fBAnZs0drq0crKDlnuGwfd57b9+Hj73xmqpWC4SQDXXBLSnJ7ghKs0yly7ZH7GMDba7al0XiqhzYpnMyFoRuUCxnCrzOwYusZs0h8TsGQwCA04tpnF4wsxt2DHbvwPXwJgSk2a45KHevEa2FPdNhna8fa7s9m8jjF2dXkCkUcdiawV2pHPYMhezzJmNBvO7azZ9TzpRkt6SzLqTyEAWCyy3jEA3KkNrcu0iRxLKYAys6ZcphYqDkAvSK+u69T3IdErkCogEFV42bA08OTyVwcTWPoCJir2PH1m0ENyEg7daHvldgymu9xmHFmgO9lFbxjDW/mcW8WNO93UMhBGQR1+yIbvRy18RWkt2iHJIqhsKKnTK63rYiG0EWSUW2UvkAscmB0uaTGwcPspozlcNASMHuoRCevpDAUlrFlWP9Tbc+3ko4h8Ssl1yhCIGgbW2QvQT7MK/XR5/IFiCLBFqR4icnFwGUUqnZgjLS58djH7odfW2IfbHstW7JWFpIqRiJ+DFiTSNcbxbXRvBJQpVyCMii7RIciwZACECpd1yzG/okE0LOEUKOEkIOE0IOWcdihJAfEEJesL4POO7/QULIKULISULIHY7j11mPc4oQ8teknb10HSSymt1B9NodUTx1Po4Ts0l7dGK3EnQMiVkv5hQ4qa1tkL0C+zCvt0p6JVPApdtMlwfr68+UAwv0+yQB/QG5La9vcBOy17yEaRx8GLHiMJ0wDrJY3j4jldftTCXADFiPWsbL5xHlsBnbvF+ilF5DKT1o/f4BAA9RSvcBeMj6HYSQ/QDuAXAAwJ0APkEIYa/CJwHcC2Cf9XXnJlxX0yQt5QAA1+4YwHKmAFU3cFW3GweZuZU2EHPwUPFOuwlsIOaQKxSh6gb2j5amfPlloeRWspRDO1/bUt1Ld7iVFlN5jPT57DnmnTAOlY33knnNrnFgTMRM15JX3Eqt8AHcBeBz1s+fA/A6x/H7KaUqpfQsgFMArieEjALoo5Q+Ss0qos87zmkrCYdxeLHDt3vVROv9vJ1EEAiCysb6K+UKRQR6MN4AOFJZ16EcWKaScwTkyy8dxmpOw2pWK1MO7WIz6l68gl40sJwpmG6lTioHqVo5VKbHs7hDtxgHCuD7hJAnCSH3Wse2UUrnAMD6PmIdHwcw5Th32jo2bv1cebwKQsi9hJBDhJBDi4uLG7x0U7I/eX4FgDnfNVsoImoZh8u2RRBURET8EnbGujdTiRFUpA1nK/ViphLgmLu8jteP1TiMRQN2Mdur9m8HYHYGdpsp3GpY2nY3KIfFtNmNdaTPdCtdPRnFdTsH1j5xkzErpKldRZ/M63amEmMyZmYseWWTtdFP802U0llCyAiAHxBCnqtzXzdnKa1zvPogpZ8C8CkAOHjw4IbLN79xeBZ/8LUjePxDt4NYl8FiDpIo4OZLhiAQ0rWzo52ENzjgxUtDStoNMw75dSiHuKUcBoIKRvsD6A/IuGLUjD9cWMnabqV2KofAJmSveYWpFTOGMzkQhCQKeOA9N3XkOhTr/1coGvBJIlI5DZMD5RXsTDl4xT27IeNAKZ21vi8QQr4O4HoAFwkho5TSOctltGDdfRrApOP0CQCz1vEJl+MtZ3Y1B0qB5XTBzhpgbiUA+Phvvrgdl+EJNtqm2exD740dT7spZSutw61kKYdYSMZbX7oTRUoxaSnVCytZqFr7jYMiCVBEYUNK0iuw2M2ODqv/Iatdxz8t+ngAABdYSURBVPvuP4w/uesAknm9KuZw/e4YrhzvwyUj3kibX7dxIISEAAiU0pT186sB/CmAbwB4K4A/s74/YJ3yDQBfJIT8FYAxmIHnX1BKi4SQFCHkRgCPA3gLgL9Z73U1A/tgJrIaFMlSDg7j4IUh3+0i5BM35EbIajqGw+2druUVZFGALJJ1GQfW7DEaVHDP9Tvs4wNBGRdWsugPyJAE0vairWCXjAq9sJKFQEpFhJ3i9ddOYH5Vxcd/fMoay6uVZSsBZlHjt377lg5dYTUbUQ7bAHzdSq2TAHyRUvogIeQJAF8ihLwdwAUAbwAASulxQsiXAJwAoAN4D6WUfZreDeCzAAIAvmt9tZxlyzis5gq27It2oEDGC4R8km0s10M8o9npmL1IQBbX5VZir3k0UL6L3BELYmolC/9IpCNuhpAidUXMYXoli9H+gP357hSKJOC9r9yHZ6YTeOzsMgq6URVz8BrrNg6U0jMArnY5vgzg9hrn3AfgPpfjhwBcud5rWS8raWYcNEiCZRwC3v6HtYqQItmFV81CKcWiVWjUqwTWOQMhkS2gzy9VKYPxgQCem09hMhbsSK+d0BZTDnrRwP1PTOHUQhrX7ojirmvMnJYLK1k70OsFrpmM4kfWkCavN/P09tW1GKdbSbCKi9o9IcorbGRucCKroVA07FTBXiSoSMhpzbfsXslqrqmVQ2EfltPLUDUzgNluNpq91m6+9vQM/se/HYMoEDxweAa/evUYCCG4sJLFbZcNr/0AbeLqyVJavNeVQ+841V0ouZU0LGVUKKKASBfPbahHUFn/gBfW9XKkr3eNg18WkWtip/3sXBK3/5+f4LEzy66uzKGwD6s5Dam8Bl+HlMNWyVailOIzPzuLy7dH8L9efyXiWQ2nF9PIa0UspNSyvkWd5mpHQa3XlUPPGgfDoHYaYSKnYSVdQCyk9GT7B6CkHNYzzYwNU+lltxKbBtcoP3thCacXM1hMqRiLVr9ubBjN7Gquc8phixiHR88s47n5FN520268ZFcMAPDEuTim41amkoc6KkeDZt82AFXZSl7D26arhSTzmj3pilWidqJy0iuEfBKKBoWqG00HQNkwlV52KwVksanGhWeW0oiFFHz5v7zU1b0wZGV+zcRz2OVo090uwj7J811ZV3MaPvKtE/jxyQXEQgp+9Zox+CQBQ2EFT5xdsSflTXqsiPXqiX6cXcpUZSt5jZ5VDsuOzBzTrVRo++hAL8GqYpN5reFA5JcPTeGxM8u2W2kjA2a2OgFFbKq30unFDPYMhbB3OOz6ujHjEM9qba1xYGy0nUo7ePzMMr785DSuHO/Hx3/jxfDLIgghOLgzhifOr+DCsqkcvORWAoDrdsVASGfaeDSDt01XC2HBaEUSkMgVsJrTsNtD8rPdsH46f/rNE/je8XnccWA7PvyaKzBaZw7xfd95FldNRLF3OISQItqP0YsE5ObcSmcWM3jF5bUDpUOOjUonUln7AzJWcxoopZ51tbIN3v96/YvK6hhesjuGB4/P4ztH5xGQxbLX0gvc85JJHBjr87wbtneVg5XGumswiFU75tC7O1/WifN7x+cxHPbhW0fm8PWnZ2reP5XXkMhqeG4uabZE7vP2G73VBJXGs71WcxqW0mrZyM9KhhwFhf4OxBxiIQW6QZHMeVc9LFmKtXIHfsNuK+5wfgWvvXrUc8ZNFgW8eEf7+zs1S88aB6Yc9gyFsZBUkSkUe9qtxHr4a0WKd718L2SR1F0YWM+ahZSK5+aSPe1SAswK3KW0inNLmTXve2bRHD1bb7pgyCfZbTk6ka3EjNNyRm37326U5UwBEb9UpayuHO/HF995Ax79wO34i7urSrE4DdLDxsF80+8eDtnNzQY97gNsJWFf6QN222XDiPhle86tG1PxUsHc6cVMTwejAdNVIIsC/u7hM2ve98yiaUD2DNcPNA9FzPdjp5QDUB6b8xpLabVMYTl52d4heywoZ330rHFYzhQQ9klli5rXA0SthM0N3jUYxM7BEMK++nUPldXUXveftpqRPj/uvm4CX31yGgvJfN37nllKQxLIms3gBi03ZyeUA1PRzP3qRZbSak9v6FpNzxqHlYxZ1+BstNfLbiU2R/q2y8zxGxG/hFS+tnGYjufKjGsvF8Ax3nXrHhSKBh44XL+p8JnFDHbEgms2dmS74k4EpJlh8rRbKV2oqRw4G6fnjQMbsAKUPhC9yGi/H79+cBJvvnEnAGYc6riVVrKYGAjgCmu8Za+7lQBg52AIA0EZZ5fd4w7Zgo7/70cv4JFTS2u6lABg2HIrdSKVdSBkfi5WPKwclns8/bzV9KxxWE4XMFihHHq1rxJgDjf687uvsnvJmzGHOm6leBaTsSAutwbT9LpbiTFpdVN144HDs/jL7z+Py7f34X2vvHTNx2KblU4oB59kTkH0asxBLxqIZwsY5MqhZfSscSi5lUyDIIukZ/squRHx1XYrUUoxtZLD5EAQL94xAEKAnT1cI+JkciCI6XjO9TaWevlP77gBV473u97HCcvP74RyAMwEDa8ah5VsAZQCwz28oWs1PWkcVnMa5pN57IgFbbfSYMjnuXzoTlLPrbScKSCnFTEZC+DV+7fhJ793m+daFHSKiYEAZuI5GEZ1j6qVbAERn9TwbIGhCAtId2Zs5GDYh+W0N2MOLFDOlUPr6EnjcGQ6AQC4ZkfUdiv1cqaSGxG/jLSquzbiY26TyYEgCCHYOdj+3j9eZSIWRKFo2C1FnMQzBQw08T6zA9IdUg6xkLKhAVDNcOjcSt0YVyVLltHi2UqtoyeNw+ELpnG4aiIKWRQQUkQe2Kog7JdgUPe5yCfnUwC819DMC7Ch8c46EMZKVmvKOOyIBSEKBNs6VH3eLrdSMq/h1z/1GP7psQsNn8OUwxBPhGgZvWkcphLYOxyyVcO2Pj9GecFMGazXfGXc4fjsKj7y7Wdx2bZIQxk3vQYzmG5B6XimgFiw8TbNY9EAHv3AK3DLvqFNu75mGAybysHNRbaZTK/kUDQoztfI8nKDKYehHs4wbDU9F4GllOLwVMLO5weAf3jrQfT16HjQWrBe82lVA2AaTsOgeNcXnkTEL+Gzb3vJmnn6vci41QCOtRdxspIpYN+22i0z3Ohkz6pYyIeiQZHMay2drc7mLswk3AP5biylC5BF4vm211uZnvt0T8dzWM4UcM2O0ri+PcNhXkxTAcvcSjqUw1MX4piO5/CBX768brfWXsYvixiJ+FzdSvFsAbEWLrKbDcuWWmpxrQMzCjPxHCileM8Xn8JDz16sut8/PXYe//7CIgBgOa3yJJIW03PG4ekpM95wrWOWK6caN7fSg8fmoYgCXnH5SK3TOHCvdchrRWQLxaZiDp2GJWm0Oig9Y6X+ziRyWEyp+PaROfzw2YWy+6h6Ef/zWyfwmZ+dBcAL4NpBzxmHZE7D9j4/Lt8e6fSleBrbrWQZB0opvntsHjfvG/L8eMNOMzkQqKp1SGTNTJytlBVnt9BocTorUw6qbuCxsysAgLnV8tfvmalVqLqBC5bRXUypPI21xfSccXjzjTvx6AdfAYn7y+sStpWDuagdn01iJpHDnQe2d/KytgSTsSDmVnNQ9VKmF9t9D2wht5LdfK/FymE6noMkmO6hnzxnKoa5RHnzwsfPLAMApuJm8PrsUga7eOFlS+nJFZL7Kdem0q307y8sAQBeuX9bx65pq7B7KASDwh5TCZjxBmBrKYdYSIEoEMyv1u8yu1FmEjm7Yvynz5sxhdkK5fDYWdM4FHQDR6YTSKu63eqF0xp60jhw1iasSCAESFltuy+sZBELKVtqcesUbIjPGcfgH6YcYqGt45KTRQHj0QDONZFi2izZgo6VTAHXW9PbmEpJ5XW7ZXxBN/Dk+Tj2WcbgR5a6qDcsibNxuHHguCIIBGGl1EJjJpHDxADPUGoEVv/BhvoAJeWwldxKALBrKITzy+6NBDeDWSvesH+0zx5Vy+qP5qzbjs4kkNcMvPHgJADgIStYzZVDa+HGgVOTsGOmw3Q8a+fwc+oT8csYjvjscaCAqRwIQVkX4K3ArsEgzi1nXNuobAZTVuB+YiCAcWvzcfMlZtHfrOXOeuyMGaS+69oxSALBibkkIhWDujibDzcOnJpE/BLSebO/0kycK4dm2DMUKnMrxTMF9PnlLZcIsXMwhFReRzzbeN+jRigaFB//8Sl8+8gcAGB8IGBvPl5+6TCAknJ4/OwK9o2EMRLx2+/BvSNhHjtsMby8kFOTiF9GStWwlC5A1Q2uHJpgz3AYDx6bw9mlDL5zdA5LVov4rcbuITMj6OxSZlOv/9C5FXz0eycBAJJAMBLxl5TDviEQYioHvWjgyXMreP2LxwEAOwZDOLec5S6lNsCNA6cmYZ+EeLZg56FPDPDUwUbZOxxCPKvh/V9+BofOxxGQRVwxuvVqa1jH3fPLGVy3c2DTHveh5xagiAI+/JorIBBAFAjuPDCKXMHAaL8fIxEf5hI5HJtNIlMo4obdg+b1WL2ruHFoPdw4cGoS8UuYWsnavW/GuVupYVhQ+tD5OAQC5LTillQOEwMBCAQ4t8lB6R8+exE37InhrS/bZR+7ed8QbraaDI72BzC7mrPrG27YY2YzsaFSl/BMpZaztRygnLYS8ctI5nW7vQE3Do2zZ8hcvPyygL+4+2oAWy9TCTDHhY5FA011THXjbx56AR/46hFQahawnVnM4JVX1K6ZGYv6MZfI47Ezy9gzHLLH0L5kVwzb+ny4anLtSXqcjcGVA6cmEb+EZF7DueUs+vwS+njbjIaZGAggGpTxay+ewK+9eBzPX0zhpks603p7o+weCtnK4fxyBt8+Ood3v3wv8pqBJ86t4JZ9Q3WDw8/NJ/GxHz4PgwLX747Zm43br6jdo2u0P4DvHpvHmaUM/rNDXVw9GcXjH3rl5jwxTl24ceDU5IbdMXzq4TP4+tPT9k6Y0xiSKOCH//3liAZkEELwoV+5otOXtG52DgbxjcOzoJTiUw+fwT8/fgG37hvGwy8s4i8ePIk/+48vwj3X73A9l1KK+779LCJ+GTsHg/jg145C1Q28bO9g3RjWNZNR+CURb7t5F/7rbZe06qlx6sDdSpyavOLyEVwzGUVeM7hLaR0MhX1bLnXVjasnokjmdRyZXrXbW/zgxEU8eGweAPDH3zxuTwes5KkLcfz7C0v47Vdcgr98w9WIBmW885bd+Mx/fkndv/naq8fw7P+8E++/43KEfHwP2wm2/juX0zIIIfj9Oy4DAF7j0MO8ev92yCLB/33oBUzHcxAI8JUnp3FkehVvv3k3wj4JH/76UVBK8ftfeQYf+vpR6EUDAPD1p2fglwW86foduHRbBI9/6JX48Gv2wy+LHX5WnLXgJplTl5ddMoSPvO5KvHTvYKcvhdMh+oMybtk3bPc0+s0bduILj523ft6BXYNB/L8PHMeffPMEvnRoGoDZG+mjd1+Fbx+Zw6v2b+e7/y2IZ5QDIeROQshJQsgpQsgHOn09nBJvvnEnb3LW4/yHq0YBmPUFLP3UnCMexhtfMomxfj8++/Nz2DcSxvvvuAzffGYWv/bJnyOe1XDX1WMdvHLOevGEcSCEiAA+DuCXAewH8CZCyP7OXhWHw2G8cv82BBURr7xiGy4ZCeNXXrQdb7t5FwAz3fV9r7wUAgH+6LUH8J5fugS/f+dlOD6bRDQo41arHQZna0Fa1VCrqYsg5KUA/phSeof1+wcBgFL6v2udc/DgQXro0KE2XSGHw5mOZzEU9tWMF6xUtAj58qEpBBQR/+Eqrhy8BCHkSUrpwbXu5xVH4DiAKcfv0wBu6NC1cDgcF9Zqn1JZAf4Gq8U2Z2viCbcSALcKmipJQwi5lxBy6P9v795CpariOI5/f3S6UJnlLQorC6IyCS0hKyMofKiXAoOSSKuX7EL1VkZQLz0kFVE+mKTRjbCwSLtSUpLdLUU9SakhZUkpmalRFP172GtoOHOOx5mz99l7Zn4f2MycNXsW//U/c+a/9zoza0tas3PnzmEIy8ysO1WlOGwH6g8zxgM/9d0pIhZFxNSImDp2rOcxzcyKUpXi8AVwuqRTJR0GXAssLzkmM7OuVYn/OUTEP5JuB94BDgGWRERvyWGZmXWtShQHgIh4E3iz7DjMzKw600pmZlYhLg5mZtbAxcHMzBpU4hvSrZC0F/hmiN2MBPbkEE69McCunPrKM74ixlpEv3nmD/KPr8q/k27KXRH9dXr+auM7JSIG/y5ARLTlBqzJoY9FVYyriPiKGGsR/eaZv4Liq+zvpJty5/wVP75un1ZaUXYAg8gzvqLG2k05zLs/565a/eWtrcfbztNKa+IgFo8ablWNq104f61z7oam0/PX7Pja+cxhUdkBDKCqcbUL5691zt3QdHr+mhpf2545mJlZcdr5zMHMzAri4jAISSdJel/SJkm9ku5M7aMkvStpc7o9LrWPTvvvk7SgT1+zJG2QtF7S25LGlDGm4ZRz/q5JueuVNL+M8QynFnI3Q9KX6TX2paRL6/o6L7VvkfS4pP6Wye8oOefvQUk/SNpX1niGXZ4fterEDTgBODfdHwF8S3Yp0/nAPan9HuChdP8oYDowF1hQ108P8AswJv08n+zqd6WPsU3yNxr4Hhibfn4GuKzs8VUsd1OAE9P9ScCPdX19DlxAdu2Ut4DLyx5fm+VvWupvX9njGq7NZw6DiIgdEfFVur8X2ER25boryd6gSLdXpX32R8Rq4M8+XSltR6WjtmPo55oVnSbH/J0GfBsRtas8vQfMLDj8UrWQu7URUXtN9QJHSDpc0gnAMRHxSWTvdM/WntPJ8spfeuzTiNgxnPGXzcWhCZImkB1dfAYcX3uxpNtxB3puRPwN3AJsICsKE4HFBYZbOUPJH7AFOFPSBEk9ZH/QXXMdyhZyNxNYGxF/kb0hbq97bHtq6xpDzF9XcnE4SJKOBpYBd0XE7y08/1Cy4jAFOBFYD8zLNcgKG2r+ImI3Wf6WAh8C24B/8oyxqprNnaSzgYeAm2tN/ezWNR9TzCF/XcnF4SCkN/ZlwAsR8Upq/jmdrpNufxmkm8kAEbE1ndq/BFxYUMiVklP+iIgVEXF+RFxAtq7W5qJiropmcydpPPAqMDsitqbm7WSX3q3p9zK8nSin/HUlF4dBpP8PLAY2RcSjdQ8tB+ak+3OA1wbp6kdgoqTaglczyOZAO1qO+UPSuHR7HHAr8FS+0VZLs7mTdCzwBjAvIj6q7ZymTvZKmpb6nM1B5Lvd5ZW/rlX2f8SrvpF9cibIpoHWpe0Ksk/PrCQ7el0JjKp7zjbgV2Af2VHbxNQ+l6wgrCdbJ2V02eNrs/y9CHydtmvLHlvVcgfcB+yv23cdMC49NhXYCGwFFpC+ANvJW875m59ei/+m2wfKHl/Rm78hbWZmDTytZGZmDVwczMysgYuDmZk1cHEwM7MGLg5mZtbAxcGsAJLmSprdxP4TJG0sMiazZvSUHYBZp5HUExELy47DbChcHMz6kRZqe5tsobYpZMs9zwbOAh4FjgZ2ATdExA5JHwAfAxcByyWNIFve+WFJk4GFwJFkX0K7KSJ2SzoPWAL8AawevtGZDc7TSmYDOwNYFBHnAL8DtwFPAFdHRO2N/cG6/Y+NiEsi4pE+/TwL3J362QDcn9qfBu6IbK0os0rxmYPZwH6I/9fYeR64l+wiMO+mC6kdAtSv8b+0bweSRpIVjVWp6Rng5X7anwMuz38IZq1xcTAbWN+1ZfYCvQc40t/fRN/qp3+zyvC0ktnATpZUKwSzgE+BsbU2SYemtf8HFBF7gN2SLk5N1wOrIuI3YI+k6an9uvzDN2udzxzMBrYJmCPpSbIVPJ8A3gEeT9NCPcBjZJeUPJA5wEJJRwLfATem9huBJZL+SP2aVYZXZTXrR/q00usRMankUMxK4WklMzNr4DMHMzNr4DMHMzNr4OJgZmYNXBzMzKyBi4OZmTVwcTAzswYuDmZm1uA/pLqmxS2EdG8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'][-200:].plot()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "first_september_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", " for y in range(1985,\n", " sorted_data.index[-1].year)]" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "year = []\n", "yearly_incidence = []\n", "for week1, week2 in zip(first_september_week[:-1],\n", " first_september_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": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGDpJREFUeJzt3XGQndV53/HvI60sYVuyEZYckDAiYw21IGNsdgSp29TYDYI4Y5GJ7VFNQG3o0DqkdTKZYlHokLHVGdPJxImGmsIU29gYMCFmTGsrtizHE5pSiVVsF4OsSAlEyMjWepaC8IwEQk//uGfNlbzaPbt79+57734/Mzt777nvee/ZI+3+3vec8743MhNJkmrMm+0GSJJ6h6EhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKnawGw3oNPe9KY35apVq2a7GZLUU3bt2vWTzFw20XZ9FxqrVq1iaGhotpshST0lIv6hZjuHpyRJ1QwNSVI1Q0OSVM3QkCRVMzQkSdUMDUnTcuiFI3zojkc5dPjIbDdFXWBoSJqWLdv38tjTI2z55t7Zboq6oO+u09CJDr1whN+97zvc9uF3sHzxotlujvrIeTdv5eix4z97fs+O/dyzYz8LB+axZ/MVs9gyzSTPNPqcR4GaKY/ccCnvv/AsFi1o/RlZtGAe6y88i0c+dunPbesQVv/wTKNPeRSombZ8ySIWLxzg6LHjLByYx9Fjx1m8cGDMM9r2g5fNv/FLs9BadYqh0aceueFSNn9tN9944kccefk4ixbMY935v8BN73vbbDdNfeQnLx7lqovP4cNr38K9O/czfNKZhAcv/cfQ6FOTOQqUpuqOqwd/9njzlRf83OsevPQfQ6OPTXQUKM00D176j6HRxyY6CpS6wYOX/hKZOdtt6KjBwcH01uiSNDkRsSszByfarmrJbUT8fkQ8ERHfj4j7ImJRRCyNiG0Rsbd8P71t+xsjYl9E7ImIdW3lF0XE4+W1LRERpXxhRHyplO+IiFVtdTaW99gbERsn0wmSpM6aMDQiYgXw74HBzLwAmA9sADYB2zNzNbC9PCci1pTXzwcuBz4dEfPL7m4HrgNWl6/LS/m1wHOZ+VbgU8CtZV9LgVuAi4G1wC3t4SRJ6q7ai/sGgNMiYgB4LfAssB64u7x+N3BlebweuD8zj2bmU8A+YG1EnAksycxHszUm9vmT6ozu60HgveUsZB2wLTNHMvM5YBuvBo0kqcsmDI3M/CHwR8B+4CDwfGZ+A3hzZh4s2xwElpcqK4Bn2nZxoJStKI9PLj+hTmYeA54HzhhnXyeIiOsiYigihoaHhyf6kSRJU1QzPHU6rTOBc4GzgNdFxG+NV2WMshynfKp1Xi3IvDMzBzNzcNmyCT8XXZKmba7eGqVmeOqfA09l5nBmvgx8GfjHwI/LkBPl+6Gy/QHg7Lb6K2kNZx0oj08uP6FOGQJ7AzAyzr6kvjRX/xD1orl6X7ea0NgPXBIRry3zDO8FdgMPA6OrmTYCXymPHwY2lBVR59Ka8N5ZhrAOR8QlZT/XnFRndF8fAL5V5j2+DlwWEaeXM57LSpnUl+bqH6Ject7NW1m16avcs2M/ma1bo6za9FXOu3nrbDetKya8uC8zd0TEg8DfAMeA7wB3Aq8HHoiIa2kFywfL9k9ExAPAk2X76zPzlbK7jwCfA04DtpYvgLuAL0TEPlpnGBvKvkYi4hPAY2W7j2fmyLR+YqmBvEdT75jrt0bx4j6pAQ69cOSUf4i85Ubz3PTQ49y7cz+vmT+Pl145zlVr39Lzd++tvbjP24hIDeA9mnrLXL41iqGhvtW0Ty2cqD1z+Q9Rr5nL93UzNNS3mvbBPxO1Zy7/IVLvcE5DfefkSeVRszWp3LT2TEbTztY0czp6w0Kpl0zms6vnYnsmwyXAOpnDU+o7TZtUblp7argEWKfimYb60uik8kO/8y6uuvgchl88ansmoZfPjjSzPNNQX2rapHLT2jORXjw7UncYGpLG5BJgjcXVU5IkV09JkjrP0JDmKG/DrqkwNKQ5ymswNBVOhEtzjNdgaDo805DmGK/B0HQYGtIc08/XYDhPM/MMDWkO6rUr1Gs5TzPzvE5DUs/r5TsJN4XXaUiaMyYzT+MQ1vQYGmocf6k1WZOZp3EIa3pccqvGqfnEPT8cSCeb6F5ZLjXuDOc01BiTGZe++aHH+eLO/Vy19i2N+ChXNd+hF46w+Wu7+cYTP+LIy8dZtGAe687/BW5639s88KB+TsMzDTXGIzdcespf6lEeLWqq+nmpcTc5p6HGqPml9sI0TUe/LjXuJs80VKVbcwgTjUt7tKjp6LUPw2oiQ0NVaianO6Hml9oPB5JmjxPhGtdkJqdd0ST1Li/uU0dMZg7B9e9S/zM0NO7FdDVzCOfdvJVVm77KPTv2k9la0bRq01c57+at3fwx1MO8oLN3GBqa8AxhohUnrmjSdHmW2juc05jDOnmTt5seepx7d+7nNfPn8dIrx73obpb1yvySNxpsDuc0NKFOniG4/r17aoZyeuXI3bPU3uOS2zmsk9c8uP69e8Zb/txrV8xP5v9gr5w99TtDY47zmofeURMINbdiaZra/4PdulZI43NOQ3NazdFrU45wa2+412/zS857dIdzGlKFmrH/pswP1A7l9Nv8kvMezeLwlOakmqGeJs4P1Azl9Nv8kvcba5aqM42IeGNEPBgRP4iI3RHxyxGxNCK2RcTe8v30tu1vjIh9EbEnIta1lV8UEY+X17ZERJTyhRHxpVK+IyJWtdXZWN5jb0Rs7NyPrrms5ui1iUe4d1w9yOYrL2DNWUvYfOUFJwREP+v22ZMXG55a7fDUnwJ/kZn/CHg7sBvYBGzPzNXA9vKciFgDbADOBy4HPh0R88t+bgeuA1aXr8tL+bXAc5n5VuBTwK1lX0uBW4CLgbXALe3hJE1VzdGrR7jN0e2wbMqQZBNNODwVEUuAXwH+JUBmvgS8FBHrgXeXze4Gvg18DFgP3J+ZR4GnImIfsDYingaWZOajZb+fB64EtpY6f1j29SBwWzkLWQdsy8yRUmcbraC5bxo/swTUDfW4umxuaeKQZNPUzGn8IjAMfDYi3g7sAj4KvDkzDwJk5sGIWF62XwH8n7b6B0rZy+XxyeWjdZ4p+zoWEc8DZ7SXj1FHmpaasf9+mx/Q+HpxyXK31QxPDQDvBG7PzHcAP6UMRZ1CjFGW45RPtc6rbxhxXUQMRcTQ8PDwOE2TpFNzSHJiNaFxADiQmTvK8wdphciPI+JMgPL9UNv2Z7fVXwk8W8pXjlF+Qp2IGADeAIyMs68TZOadmTmYmYPLli2r+JEkaWxNXLLcpIn5CUMjM38EPBMR55Wi9wJPAg8Do6uZNgJfKY8fBjaUFVHn0prw3lmGsg5HxCVlvuKak+qM7usDwLeyddXh14HLIuL0MgF+WSmTpBnRxBVqTZqYr71O498BX4yI1wB/D/wrWoHzQERcC+wHPgiQmU9ExAO0guUYcH1mvlL28xHgc8BptCbARz9w4S7gC2XSfITW6isycyQiPgE8Vrb7+OikuCT1uyZOzHsbEUlqqNpbx3SCtxFRIzVpbFZquiZOzBsa6qomjc1KvaBpE/MOT6krvFOp1GwOT80BvTTU08T7OEmaPEOjh/XSUE8Tx2almdZLB3a1vDV6D2riMrwa3sdJc00/ftqgcxo9qJvL8CRN3mTm8JryyZDOafQxh3qkZpvMHF4vDTODw1M9y6EeqblqDux6dZjZ0OhR3rJbaraJDux69TbshoYkzYCJDux6dZjZ0JCkWdKLw8yunpI6oCkrYKSpcvWU1EW9tgJGmiqHp6Rp6NUVMNJUeaYhTYP31NJcY2hI09CrK2CkqXJ4SpqmXlwBI02Vq6ckSa6ekiR1nqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqpVh0ZEzI+I70TE/yzPl0bEtojYW76f3rbtjRGxLyL2RMS6tvKLIuLx8tqWiIhSvjAivlTKd0TEqrY6G8t77I2IjZ34oSVJUzOZM42PArvbnm8CtmfmamB7eU5ErAE2AOcDlwOfjoj5pc7twHXA6vJ1eSm/FnguM98KfAq4texrKXALcDGwFrilPZwkSd1VFRoRsRJ4H/Df24rXA3eXx3cDV7aV35+ZRzPzKWAfsDYizgSWZOajmZnA50+qM7qvB4H3lrOQdcC2zBzJzOeAbbwaNJKkLqs90/gT4AbgeFvZmzPzIED5vryUrwCeadvuQClbUR6fXH5Cncw8BjwPnDHOviRJs2DC0IiIXwcOZeauyn3GGGU5TvlU67S38bqIGIqIoeHh4cpmSpImq+ZM413A+yPiaeB+4D0RcQ/w4zLkRPl+qGx/ADi7rf5K4NlSvnKM8hPqRMQA8AZgZJx9nSAz78zMwcwcXLZsWcWPJEmaiglDIzNvzMyVmbmK1gT3tzLzt4CHgdHVTBuBr5THDwMbyoqoc2lNeO8sQ1iHI+KSMl9xzUl1Rvf1gfIeCXwduCwiTi8T4JeVMknSLBiYRt1PAg9ExLXAfuCDAJn5REQ8ADwJHAOuz8xXSp2PAJ8DTgO2li+Au4AvRMQ+WmcYG8q+RiLiE8BjZbuPZ+bINNosSZqGaB3Q94/BwcEcGhqa7WZIUk+JiF2ZOTjRdl4RLkmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkaktQHDr1whA/d8SiHDh+Z0fcxNCSpD2zZvpfHnh5hyzf3zuj7DMzo3iVJM+q8m7dy9Njxnz2/Z8d+7tmxn4UD89iz+YqOv59nGpLUwx654VLef+FZLFrQ+nO+aME81l94Fo987NIZeb8JQyMizo6Iv4yI3RHxRER8tJQvjYhtEbG3fD+9rc6NEbEvIvZExLq28osi4vHy2paIiFK+MCK+VMp3RMSqtjoby3vsjYiNnfzhJanXLV+yiMULBzh67DgLB+Zx9NhxFi8cYPniRTPyfjVnGseAP8jMtwGXANdHxBpgE7A9M1cD28tzymsbgPOBy4FPR8T8sq/bgeuA1eXr8lJ+LfBcZr4V+BRwa9nXUuAW4GJgLXBLezhJkuAnLx7lqovP4aHfeRdXXXwOwy8enbH3mnBOIzMPAgfL48MRsRtYAawH3l02uxv4NvCxUn5/Zh4FnoqIfcDaiHgaWJKZjwJExOeBK4Gtpc4fln09CNxWzkLWAdsyc6TU2UYraO6bzg8tSf3kjqsHf/Z485UXzOh7TWpOowwbvQPYAby5BMposCwvm60AnmmrdqCUrSiPTy4/oU5mHgOeB84YZ1+SpFlQHRoR8Xrgz4Hfy8wXxtt0jLIcp3yqddrbdl1EDEXE0PDw8DhNkyRNR1VoRMQCWoHxxcz8cin+cUScWV4/EzhUyg8AZ7dVXwk8W8pXjlF+Qp2IGADeAIyMs68TZOadmTmYmYPLli2r+ZEkSVNQs3oqgLuA3Zn5x20vPQyMrmbaCHylrXxDWRF1Lq0J751lCOtwRFxS9nnNSXVG9/UB4FuZmcDXgcsi4vQyAX5ZKZMkzYKai/veBVwNPB4R3y1l/xH4JPBARFwL7Ac+CJCZT0TEA8CTtFZeXZ+Zr5R6HwE+B5xGawJ8aym/C/hCmTQfobX6iswciYhPAI+V7T4+OikuSeq+aB3Q94/BwcEcGhqa7WZIUk+JiF2ZOTjRdl4RLkmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqr1RGhExOURsSci9kXEptlujyTNVY0PjYiYD/xX4ApgDfAvImLNTLzXoReO8KE7HuXQ4SOzvk3NPiSp2xofGsBaYF9m/n1mvgTcD6yfiTfasn0vjz09wpZv7p31bWr2IUndFpk5220YV0R8ALg8M/91eX41cHFm/u5Y2w8ODubQ0NCk3uO8m7dy9NjxnytfODCPPZuv6Oo2NfuQpE6LiF2ZOTjRdr1wphFjlJ2QdBFxXUQMRcTQ8PDwpN/gkRsu5f0XnsWiBa3uWLRgHusvPItHPnZp17ep2YckzZZeCI0DwNltz1cCz7ZvkJl3ZuZgZg4uW7Zs0m+wfMkiFi8c4Oix4ywcmMfRY8dZvHCA5YsXdX2bmn1I0mwZmO0GVHgMWB0R5wI/BDYAH+70m/zkxaNcdfE5fHjtW7h3536Gx5iA7tY2NfuQpNnQ+DkNgIj4NeBPgPnAZzLzP59q26nMaUjSXFc7p9ELZxpk5teAr812OyRpruuFOQ1JUkMYGpKkaoaGJKmaoSFJqmZoSJKq9cSS28mIiGHgH07x8puAn3SxOdPVa+0F29wttnnm9Vp7YXptPiczJ7w6uu9CYzwRMVSzDrkpeq29YJu7xTbPvF5rL3SnzQ5PSZKqGRqSpGpzLTTunO0GTFKvtRdsc7fY5pnXa+2FLrR5Ts1pSJKmZ66daUiSpqGnQyMiPhMRhyLi+21lb4+IRyPi8Yj4HxGxpJQviIi7S/nuiLixrc63I2JPRHy3fC1vSJtfExGfLeXfi4h3t9W5qJTvi4gtETHWh1U1rc1d6eeIODsi/rL8Oz8RER8t5UsjYltE7C3fT2+rc2Ppyz0Rsa6tvCv93OE2N7KfI+KMsv2LEXHbSfua8X7ucHub2se/GhG7Sl/uioj3tO2rM32cmT37BfwK8E7g+21ljwH/rDz+beAT5fGHgfvL49cCTwOryvNvA4MNbPP1wGfL4+XALmBeeb4T+GVan2y4FbiiB9rclX4GzgTeWR4vBv4WWAP8F2BTKd8E3FoerwG+BywEzgX+DpjfzX7ucJub2s+vA/4J8G+B207a14z3c4fb29Q+fgdwVnl8AfDDTvdxT59pZOZfASMnFZ8H/FV5vA34zdHNgddFxABwGvAS8EI32tlukm1eA2wv9Q4B/w8YjIgzgSWZ+Wi2/jd8HriyyW2eqbaNJTMPZubflMeHgd3ACmA9cHfZ7G5e7bP1tA4ojmbmU8A+YG03+7lTbZ6JtnWqzZn508z8X8AJnyrWrX7uVHu7aQpt/k5mjn6y6RPAoohY2Mk+7unQOIXvA+8vjz/Iqx8V+yDwU+AgsB/4o8xs/0P42XKa+Z9mcqjnFE7V5u8B6yNiIFqfXHhReW0FrY/BHXWglHXTZNs8qqv9HBGraB197QDenJkHofXLSOtMCFp990xbtdH+nJV+nmabRzWxn0+l6/08zfaOanof/ybwncw8Sgf7uB9D47eB6yNiF63TuZdK+VrgFeAsWqfzfxARv1heuyozfwn4p+Xr6u42+ZRt/gytf9whWp9c+L+BY7ROL0/W7WVwk20zdLmfI+L1wJ8Dv5eZ451Vnqo/u97PHWgzNLefT7mLMcpmrJ870F5oeB9HxPnArcC/GS0aY7Mp9XHfhUZm/iAzL8vMi4D7aI31QmtO4y8y8+UybPLXlGGTzPxh+X4YuJfun+aP2ebMPJaZv5+ZF2bmeuCNwF5af5RXtu1iJfDsyfttWJu72s8RsYDWL9kXM/PLpfjH5TR9dEjkUCk/wIlnQ6P92dV+7lCbm9zPp9K1fu5QexvdxxGxEngIuCYzR//+dayP+y40RlcxRMQ84Gbgv5WX9gPviZbXAZcAPyjDKG8qdRYAv05r6GXW2xwRry1tJSJ+FTiWmU+W09HDEXFJOS2+BvhKk9vczX4ufXIXsDsz/7jtpYeBjeXxRl7ts4eBDWXs91xgNbCzm/3cqTY3vJ/H1K1+7lR7m9zHEfFG4KvAjZn516Mbd7SPpzJ73pQvWke4B4GXaSXptcBHaa0w+Fvgk7x6AePrgT+jNTn0JPAf8tUVEruA/1te+1PKKpQGtHkVsIfW5Nc3ad2FcnQ/g7T+o/4dcNtonaa2uZv9TGvFS5b3+m75+jXgDFqT9HvL96VtdW4qfbmHtlUl3ernTrW5B/r5aVqLKl4s/5fWdKufO9XeJvcxrQO4n7Zt+11geSf72CvCJUnV+m54SpI0cwwNSVI1Q0OSVM3QkCRVMzQkSdUMDUlSNUNDklTN0JAkVfv/ZquYoDCrbPcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "yearly_incidence.plot(style='*')" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1986 0\n", "1987 0\n", "1988 0\n", "1989 0\n", "1990 0\n", "2020 221186\n", "2002 516689\n", "2018 542312\n", "2017 551041\n", "1991 553090\n", "1996 564901\n", "2019 584066\n", "2015 604382\n", "2000 617597\n", "2001 619041\n", "2012 624573\n", "2005 628464\n", "2006 632833\n", "2011 642368\n", "1993 643387\n", "1995 652478\n", "1994 661409\n", "1998 677775\n", "1997 683434\n", "2014 685769\n", "2013 698332\n", "2007 717352\n", "2008 749478\n", "1999 756456\n", "2003 758363\n", "2004 777388\n", "2016 782114\n", "2010 829911\n", "1992 832939\n", "2009 842373\n", "dtype: int64" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "yearly_incidence.sort_values()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }