diff --git a/module3/exo2/analyse_varicelle.ipynb b/module3/exo2/analyse_varicelle.ipynb index c012d1ee49d6d1e52a45eec90a7551954f745ead..4ccf83f1abc2ae41cb54b9c7320dd68a637ae5de 100644 --- a/module3/exo2/analyse_varicelle.ipynb +++ b/module3/exo2/analyse_varicelle.ipynb @@ -14,13 +14,6 @@ "import urllib.request" ] }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - }, { "cell_type": "markdown", "metadata": {}, @@ -32,7 +25,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 2, "metadata": {}, "outputs": [], "source": [ @@ -45,7 +38,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 3, "metadata": {}, "outputs": [ { @@ -1010,7 +1003,7 @@ "[1531 rows x 10 columns]" ] }, - "execution_count": 4, + "execution_count": 3, "metadata": {}, "output_type": "execute_result" } @@ -1029,7 +1022,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 4, "metadata": { "hideCode": false, "hideOutput": true @@ -1079,7 +1072,7 @@ "Index: []" ] }, - "execution_count": 11, + "execution_count": 4, "metadata": {}, "output_type": "execute_result" } @@ -1090,7 +1083,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 5, "metadata": {}, "outputs": [], "source": [ @@ -1117,7 +1110,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -1131,6 +1124,120 @@ "data['period'] = [convert_week(yw) for yw in data['week']]" ] }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + "periods = sorted_data.index\n", + "for p1, p2 in zip(periods[:-1], periods[1:]):\n", + " delta = p2.to_timestamp() - p1.end_time\n", + " if delta > pd.Timedelta('1s'):\n", + " print(p1, p2)" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmYHUW5P/55zzmzZF+HkA0SQghLhEBi2FSWgMQNVEADKly/eFHE3asXrveHaxS9LIoKCuIFvCogoCB7SMImISFsCUkIWUlC9n2dzMw59fuju7qrq2vrM2fOTDj1eZ55pk91dVV1d3W9+1vEGIOHh4eHh4eIXGcPwMPDw8Oj68ETBw8PDw+PFDxx8PDw8PBIwRMHDw8PD48UPHHw8PDw8EjBEwcPDw8PjxQ8cfDw8PDwSMETBw8PDw+PFDxx8PDw8PBIodDZAygXAwcOZCNGjOjsYXh4eHgcUHj55Zc3M8aabPUOWOIwYsQIzJ07t7OH4eHh4XFAgYjedqnn1UoeHh4eHil44uDh4eHhkYInDh4eHh4eKXji4OHh4eGRgicOHh4eHh4peOLg4eHh4ZGCJw4eHh4eHil44uBRM5i/ZgdeX729s4fh4XFA4IANgvPwyIqP/eZ5AMDKaz/SySPx8Oj68JKDh4eHh0cKnjh4eHh4eKTgiYOHh4eHRwpW4kBEjUQ0h4heJ6IFRPTDsPwHRPQOEb0W/n1YuOZqIlpKRIuJ6ByhfDwRzQ/P3UREFJY3ENE9YflsIhpR+Vv18PDw8HCFi+SwH8CZjLHjAIwDMJmITgrP3cgYGxf+PQoARHQ0gCkAjgEwGcDNRJQP698C4HIAo8O/yWH5ZQC2McYOB3AjgJ+3/9Y8PDw8PMqFlTiwALvDn3XhHzNcch6Auxlj+xljKwAsBTCRiAYD6M0Ym8UYYwDuAvBx4Zo7w+P7AEziUoWHh4eHR/XhZHMgojwRvQZgI4BpjLHZ4amvENE8IvojEfULy4YCWC1cviYsGxoey+WJaxhjbQB2ABigGMflRDSXiOZu2rTJ6QY9PDw8PLLDiTgwxoqMsXEAhiGQAsYiUBGNQqBqWgfg+rC6iuNnhnLTNfI4bmWMTWCMTWhqsm5k5OFRETDG8PC8tWgrljp7KB4eVUMmbyXG2HYATwOYzBjbEBKNEoDbAEwMq60BMFy4bBiAtWH5MEV54hoiKgDoA2Brpjvx8OggPPT6WnzlL6/itudWdPZQPDyqBhdvpSYi6hsedwNwFoA3QxsCxycAvBEePwRgSuiBNBKB4XkOY2wdgF1EdFJoT7gEwIPCNZeGxxcAmBHaJTw8Oh1bdrcAADbsbO7kkXh4VA8u6TMGA7gz9DjKAbiXMfYwEf2JiMYhUP+sBPBFAGCMLSCiewEsBNAG4ErGWDFs6woAdwDoBuCx8A8AbgfwJyJaikBimFKBe/Pw8PDwKBNW4sAYmwfgeEX55wzXTAUwVVE+F8BYRXkzgAttY/Hw6AxwvzkvzHrUEnyEtEdN4u+vrsEX7nzJqS73lvCkwaOW4LOyetQkvnnP6851fciNRy3CSw4eHh4eHil44uDh4QhvcvCoJXji4OFhQbW1Ss8t2YSfP/5mdTv18JDgiYOHhyNYlUzSn7t9Dm55ellV+vLw0METBw8PC7w52qMW4YmDh4eHh0cKnjh4eNgQGh28QdqjluCJg4eHBT4IzqMW4YmDh4eHh0cKnjgcYFi4dic27drf2cOoSXi1kkctwROHAwwfvuk5nH3jM509jJqCz57hUYvwxOEAxPa9rZ09hBqFFx08ageeOHh4WEA+0sGjBuGJg4eHhweAjbuacc9Lq/xe4SE8cfDwcIQ3SL+78aunluA/75+PV1dv7+yhdAl44uDhYUG8E1znjsOjY7Fs024AQKuXHAA4EAciaiSiOUT0OhEtIKIfhuX9iWgaES0J//cTrrmaiJYS0WIiOkcoH09E88NzN1G4iwoRNRDRPWH5bCIaUflb9fAoD3EQ3LufOoy46hFc/cC8zh5G5+Ld/5qd4CI57AdwJmPsOADjAEwmopMAXAVgOmNsNIDp4W8Q0dEApgA4BsBkADcTUT5s6xYAlwMYHf5NDssvA7CNMXY4gBsB/LwC91ZxPLFgPVZs3tPZw/CoMmrNlfWvc1Z39hA6BdzxwNOGAFbiwALsDn/WhX8MwHkA7gzL7wTw8fD4PAB3M8b2M8ZWAFgKYCIRDQbQmzE2iwU7td8lXcPbug/AJOqCezN+8U8v44zrnu7sYXh4eHQAvPowCSebAxHlieg1ABsBTGOMzQYwiDG2DgDC/weF1YcCEFmPNWHZ0PBYLk9cwxhrA7ADwIBybsjDo6PgF413N3I8waKXHQA4EgfGWJExNg7AMARSwFhDdRXHzwzlpmuSDRNdTkRziWjupk2bbMP28KgIfJxDbcBLDklk8lZijG0H8DQCW8GGUFWE8P/GsNoaAMOFy4YBWBuWD1OUJ64hogKAPgC2Kvq/lTE2gTE2oampKcvQPQ4gbNvTgl3NPgq80mgrlvDga++A+dXPCP90Arh4KzURUd/wuBuAswC8CeAhAJeG1S4F8GB4/BCAKaEH0kgEhuc5oeppFxGdFNoTLpGu4W1dAGAG68Iz+AcPLcDGnc2dPYx3LY7/8TSc/LMZzvVb2kqYcussvLJqWweO6sBfNG59bjm+fvdrePC1tfbKNYwuvPRUFS6Sw2AAM4loHoCXENgcHgZwLYCziWgJgLPD32CMLQBwL4CFAB4HcCVjrBi2dQWAPyAwUi8D8FhYfjuAAUS0FMC3EHo+dVXc8cJK/Nff51esvbe37MHm3T7Tqojd+9uc6y7fvBsvLt+Kq++v3DtJ4ADXKi1evwsbdjZjy+4WAPBZfTUgv6lTAgVbBcbYPADHK8q3AJikuWYqgKmK8rkAUvYKxlgzgAsdxttl0Fqs3Aw67X+eRj5HWPbTD1esTQ899rcV7ZUUOFAXjXN++SyIgMs/cBgAoK10gN5IB4PzAKUD9UVXGD5CukxUevoU/Qfbbrh6mTz31uZM7b4bguAYAwq54E587iAzPG0I4IlDmfB6yQMX3erz9koCumDITVko5ILP3UsOaoS000sOITxxKBN+/hy4yEocIlT5nVeaAYkkh5JaclD1d/2TizF7+ZaKjqOrgjMBnnYG8MShTBzIKoZ3G7LGITQUsk37zpIbKs2AFPJ6yWFfSxFrtu1Llf96xlJ8+tYXKzuQLopIfeg5PwAOBmkPNTpj/vhJWxmUG9R2oD99LjkUFc4UX7jrJfxraW1ICDpQpFbq3HF0FXjJoUzkc9XnJ9fu8LEVlUBWqa+zTA6VXqPykVop3XKtEwYRXisQwBOHMlHoBOJw/ROLq97ngQRXwepAEcAqLSnW5YM56/cr0MHbHER44lAm6vLVf3Tei6JzcaCr9XI5vvi53UepxlbJXJRbqbbuWwdPHMpEZxAHP2XNcFX/ZCWy7xa1Ui5jBHCxxhZJ8q6sCXjiUCYK+eqvGH7OVgblBhzarpq3ZjtGXPUIlmzYFZW1Fkv45+try+JGK/2+s0YA12pgpv/OAnjiUCa85HDg4hM3v5CpfrRDmOUFPDxvHQBg+psbo7Jbnl6Gr/71VTz+xvpsg4S7YfQPzy3HSs0OhUdf83h0nMvox19rxIG8zSEBTxzKRF2nSA5+1prQUY/HVa2kqrYu9DDburelcgMSsLO5FT95ZBEuuk0di7C3Jc4j9eKKwCPJWXKosfnm1UpJeOJQJjrDlbWW5+xLK7fioderk2qaMYaZizdGnPOMUBJoz+Mv5925XMONxnscstguXLszbNit/1ozSJM3SCfgiYMj5AnTKUFwNaxYuvB3s/C1v75alb6eXrwJn//fl3DL00uxt6WtXfsftNeYfd0Ti/GqYZ8KPg9d8j/F6SHs86i5tejVSjUOTxwcIX9PnTF/3m0Mza3PLsOidTvb3Y7LApyFC966J1ABLdu0J9oDIQsq9Z4YA34zc6nRRsK7chFkY4O0ve6R/9/jEXHoDCm5U+C3CU3AEwdHyPOlc9JnVL/PjsRPH30TH77puar01Zoh2VxdmHuppa2EHftajXUT6IQ1lEsBLpJDmJTVWad+xnVPAwDy75KstDbwuywyht8/swzrazwjgScOZaP6K/W7Ua1ULYLXlmFzpvrQE62lWMq2MVCF78VlEY/UShnadR3mntCYXSO0ISKwyzftxs8eexNf/L+XO3lEnQtPHBzRFWwOXhdqhunx6IiD6j3WF4JFoqWtBDHTRLUf/+urt1vrcIbBRXLg95rV4JqrEeoQqd3CD22XIDXWIqzEgYiGE9FMIlpERAuI6Oth+Q+I6B0iei38+7BwzdVEtJSIFhPROUL5eCKaH567icIZTUQNRHRPWD6biEZU/lbbB69WOrChUyupUJ8P9ntoLZaSex/Ynn+F19CL/zDbWofbBVzW75g4ZBtHJW0Oq7fujWw6XRWc0NaaK68MF8mhDcC3GWNHATgJwJVEdHR47kbG2Ljw71EACM9NAXAMgMkAbiYivrvKLQAuBzA6/Jscll8GYBtj7HAANwL4eftvrbJIG6Q7xSTdCX12DDrCXdC0hGXxvOFrYWux1OU9dq574i0ASNhGdOB30pnpQ97/i5k45drplWuwgpDvs9bjHazEgTG2jjH2Sni8C8AiAEMNl5wH4G7G2H7G2AoASwFMJKLBAHozxmaxYGW4C8DHhWvuDI/vAzCJXOTkKkImBl5yaB864l5MTeoykaqu4fQgUCsxoa7boFX1OurV/XNe4Gbb0maXjDhBzkrvKu2t1NyqH+u+liLOufFZvPz21or26QL5LjMIm+9KZLI5hOqe4wFwefcrRDSPiP5IRP3CsqEAVguXrQnLhobHcnniGsZYG4AdAAZkGVtHo0u4snZCnx2Fat+LqwQwe/mWKK6gpcgySQ6qTYQ6msMpL2dTtmuq6a20YO0OLN6wC1MfWVS1Pjlk24qXHBxBRD0B3A/gG4yxnQhURKMAjAOwDsD1vKricmYoN10jj+FyIppLRHM3bdrkOvQOgd8JLsCcFVvx9OKN9ooSqn0vrmv8p299EddPC1Q1LW1F5cY4Byr4I39q0UbMWua+uU81hfiWUMLrjNxlcZxD8KC6ukqxo+H0BoioDgFh+DNj7AEAYIxtYIwVGWMlALcBmBhWXwNguHD5MABrw/JhivLENURUANAHQEquZIzdyhibwBib0NTU5HaHHYTOsDl0xbn6qd/Pwr/970uZr6vkrbgsXTpiZCJSyzbtwd6WNqFu1pElOmrHxXpkmRMiJ6zLxaRCNWPguFfZ+p3N+Mer71SvYwFFTxwAuHkrEYDbASxijN0glA8Wqn0CwBvh8UMApoQeSCMRGJ7nMMbWAdhFRCeFbV4C4EHhmkvD4wsAzGBdjE1OjaZT0md0HtZu34eLbn0R89bY3Su7Isp9dtc/+Vb2vqr4oqqh+qhmhDT3Dnt7y158457XqtJnS1sJZ1z3NJ59K9BGcPNUrauVCg51TgXwOQDziYi/rf8CcBERjUPw3a0E8EUAYIwtIKJ7ASxE4Ol0JWOMRxJdAeAOAN0APBb+AQHx+RMRLUUgMUxp321lR1uxhGkLN2Dy2IOVYnTKIF2tgYl9duJknb1iC2Yt34I/v7gKx17Qt93tVftWyu0vGSFtrmvUvnSQaibLfZX7yKsZ59DSVv05vmFnM1YIKc+LIYGqdcnBShwYY89DLbk/arhmKoCpivK5AMYqypsBXGgbS0fijhdW4iePLMKvpozDeePSzlgpg3SNcRXc2LqvNUPEsAEdoZYzvxONWsnaaLmj6Xood87mqqj+b8voIrSzuRXb9rTg0AE9yu4z7cKa/F+r8BHSIXaGHOLyTepNU7qAVqlTXVmzMI8L1+7EG+/sMNY5UCSHd88AAhtKOaimt5LscmwjaOff/AJO+5+n29WnrCngEdJereQBAOjdrQ6APpioGukz2oolFAxeGl0ht5LLCHgyvZXXfqRjB5MBunHb3iNLHJsrq5bQSqyrnR3xU021kiw4lBhg2ldrycbd7e5Tbp57qNW6WslLDiEawkycolj71MINGHHVI1i2KT0BO2LatGgCtThGDixfdG4vOHdVKXVaJYlrltQR7UFnMZKdHQ2a68SU3dXg3uUeuLdSjQsOnjhwcCZB5JIemR/sCfz66u2K3EqVnzm2KNc+oXRz3LA+Fe/bFZW46217WvDHf62oQEsBXF5FuVKXSHg6a63o7GQB1aQNnZHCQt7rg//uCpJ6Z8IThxB8EqpEaMY6JkJ65eY9iV2+9luIA5/D1Z6yb7yzI96FrQKdf/f+efifJxa3v6EM0K0x5aiK7H1VNn1Gp0sODsSpVGK456VVTmk8TJC7cqUN7dnSVCZAXJ1U65KDtzmEUM0tPk/3t5Wwcae08UcFJs7p4WYqHLYPq9RJ4u7fhWAk02K6cVczejXUWdvbvrf6WTnLfWY2VZ8I1RqqSqlxoMElzuGf89biP++fj7Xbm/HNs4+oWN+u762txFBfpogj2xb4d+YN0h4AYm5P9YH/19/np+t3AP9ulRy6uLg7cep0vGdo56m8TCj3mYlJ4spZKyrxrkyM+/knDMP9r6zBheOH6Su1Ey6SA3fkaG86bpmYmhZosS+T8fiuWSsxe/lW/PYzJyjPy33s3BdExXfNr6x68GolCa6eGe1lKq5+YF6qzC45hP+rnC1SvFfbfc+3uLCq22/fw3S5+tHQfpTuu11dV6VNk/RRHzpS9GzsOD4vC0NeLjFcu30fgGw2hxN+PC06Nu3Xcc2DCyL7oQoyXZmzMsjcU+OCgycOHLHNwa1+eyfOX+esTpXZVBiRWql9XbcLlfhg5MVuwdqd7W8U+udy16yV+O3MZZnaGn9oP3slCeqsrJXwZVUXM8bw1zmrwuP2d6NDR6fP+Mer7+CUa2dg9vJ0MkD5tkZc9Qi+/+AbqXrFDNvApq6tcZdVHTxxCKHyVjKhI1Q7Nsmh3G0e24uOVmNl0eurYHsc1zy4IHOb6sCvzllEdDOyUtHqNphibzZItrhyiOHLbwdOGW+u35U6xxRT485Zb6fKXHb6G/29R7FkQ7oPTxzU8MQhRCS+VkFy0E1GmwGsswxkCbVSJRbIA8BGq0oZ0WlxDprnJU6jjmQY6jXE4dVV23DiT6fjb3PTUrAKunlv4sdc57zLAt9aZEqJvdYNzzp44hAipg3CTDVM2vZMJ52E4EocOnMud8XvqCMkG5UqxbUXVb2OeG6iKqUjX0tBE6K8ZEMQHDpnRZxd3/Qu7pq1MnPfrvfV5qhWKiokDC85qOGJQwiW0eYwbeGGsvvSEgHLHOXal870VuqIntu7cLbnet215aSMULqyVsTkoG5EjOYvlhi2tdNTSIdChTLvbdq131pHfh+uXL1zPISiXlbasH1vC6597E20tVMd2tXhiUMIPrmqkUdGNxdtk5R1AcmhVlDNPQxs0E1J0ZXzz7NX4fgfT9PmBmsP6kzJjQD87eU12LM/sH+0xwDPGEsRA1fi4Kx+UtTLqlaa+sgi/O6ZZXh8wfpM19mwtAJ5oiqJmicOb23YhcP/61Gs2roXQFJy6KgAJt1kdFUrVVtHKuqzO6Zr90aV0ccdMKaCSq1UBhdbGe8uNc6+8dlU2c4OIA4mgzTHDdOCiHdRqp27cmvCAKx7FOL9pVPju43R9TGrIqmzqpV4PFIl1VHTFm7AWTc8g3++vtZeuUqoeeJw95zVaCsxPP5GwAVUI49NmVqlTkufkYR777qUBu15wpUmBDoVXTUzkXYW2pNyAkDiRbYqdP4X/G5WgoC5ZMCVq7iri8o3XJf7HB6bv75ii/lbIRFdtK4ybt2VQM0TB74G8MlVlTVBMxedvZWqSB1Wb92bcB3UEjbFiVueyRZb4AKlsbcDHogyx5bjteKlWeeT7BoatJFu5MHX1PsrZ+mvI6aRSdrWvSfx/uR5VGmbg1qt5HYtBx/u4wvW46s859i7EJ44hP/5BKlGLhzdhLepLVxpQ6nE8NuZSzHiqkfwo38uNNZdv6PZaMj8zB9mJ37PUgQqAeoPbLbgxWJCFmnAuFhUcLXLW/TsJrRHuvnvf6QDvOSRtBZL+Prd7d9f2V1NVnkyojPmyj25EwdH24RKreSNeEp44hB+ecUqSg66qWibo3G2SHPFmYs3RllPbamxT/rZdBwvpCGQISfJ29sSB16JInoWO0h7nrGqm47wVipn9zPTFS6L18tvb1MTaqnhJwyG0CxqUVeO2cVWENd1I94/ejjNtDCFXqnSNgduWxSRVa3UkUtEVyJTVuJARMOJaCYRLSKiBUT09bC8PxFNI6Il4f9+wjVXE9FSIlpMROcI5eOJaH547iYKZzIRNRDRPWH5bCIaUflb1d4fgPI8gcrlqHTX2eaoa/oMlf63kuDpu3//bKw2UhEH148oy2ir5carjHOwdK3K35NlITn/lhcw9+1tqXK5jd6N+sy3WfqrumODcDxz8UYAQcI+MeGk/H4rbXN4ZdV27NnflijzcQ5quEgObQC+zRg7CsBJAK4koqMBXAVgOmNsNIDp4W+E56YAOAbAZAA3E1E+bOsWAJcDGB3+TQ7LLwOwjTF2OIAbAfy8AveWCRFXLkxOGxNW7pzSXWeb4Py0rV4l3TBVnOhDoRFuhbAnsSp7ge75uaruKuWZdNJh/TNfkzV5xtY9LWXv0Wwdi/QgebI9FTIRWtfKGRp1fbfcpnPcD5+M80MpxqSb6/KCnmVeyMQhK5HsSKeVrqThshIHxtg6xtgr4fEuAIsADAVwHoA7w2p3Avh4eHwegLsZY/sZYysALAUwkYgGA+jNGJvFgq/+Luka3tZ9ACZRR74BAbLNoWL6bwN03K+tOdcIaZUbZkdAHEaWZ+HK/VdKhTRmUC/DWNzhqn5Q3d/Lq7Zn6MmMlxXSRdR3R7wHB0OyS5vi2HQz1GRzEK+fvigZhJplXsip8btC+oyu6ByXyeYQqnuOBzAbwCDG2DogICAADgqrDQUgJjBZE5YNDY/l8sQ1jLE2ADsADFD0fzkRzSWiuZs2bcoydMNNJX9mWizKVivpyt28lWzddsaev+1RK+ng6pnUId+2YvCthohY27trj8ujvHCYdtHLxtyUOaAyIY5N5yqclhzU5/704ttSPTcJA0gneswa6NyxNofOJ1QczsSBiHoCuB/ANxhjJmdcnURuktSdpHjG2K2MsQmMsQlNTU22ITtBFoNduJu4bnl96q6zfayuk7gcY6oOpqbE+1DupKe5OP3M1e2rPnhVP7YPynRWt6ir1CNG4mAcQfvQUXaESruJ2vsTfujSkKcjHRRHwHNLNjv1qSQOkuTQFbyVuuKOgU7EgYjqEBCGPzPGHgiLN4SqIoT/N4blawAMFy4fBmBtWD5MUZ64hogKAPoAcPODbCfK3bMWMH9cjDGMuOoRXP9kmsvTqpVsi1wkOZjrVSgVjhXJyOl2GKQzRIxXK125iq61Gah3Uv1R6bFU3gMJsI/zsKYeOGJQT+151ahMi5z4jNR7tTOL5OD2/G3lMnHIPKe63jreIXDxViIAtwNYxBi7QTj1EIBLw+NLATwolE8JPZBGIjA8zwlVT7uI6KSwzUuka3hbFwCYwTpwFfj072fhV08tUZ7LItaZPkQuuv5m5lLn6yrlrVSJJ5f18as4NK1Buh1p0ZWSg2WolZpJg3o3GjqpTB/tRSbbj0vIMjJ68BkehDg23RSYtyZpm0kQXVO/GcaYVivFF5982AAM7dvNvbEKY9ueFqxWuNt2Blx4zFMBfA7AmUT0Wvj3YQDXAjibiJYAODv8DcbYAgD3AlgI4HEAVzLGuHP8FQD+gMBIvQzAY2H57QAGENFSAN9C6PnUEWCMYfaKrbjxqbc0593bMn2InDtR6svKDIJri+IczOMycbgca7fvw8w3NyrP7dnfhpFXP4qbn04TNhFiL397eY2iRvtYLOV9Gm5t+Wa1t5DRSKopV428W11eUWodVoSZizdixFWPZE6wluUpZiHq1kSPCCSBSunBbZIDANw7NzmPxNsx3Vp7JAeRONQZPME4OkIFxB/HvXPX4P2/mFnx9suBdeNZxtjz0M/PSZprpgKYqiifC2CsorwZwIW2sVQCuyU3NvnGbJ/B5R84DLc+uzyoq1FBtxZLmLcm2EtZLT6rr3MNgrOm2XAgDqdcO0N7bn2YwuFvc1ULfgxxIXpREzntAt1olWollUHa1r5Q4YrTR+GWp+1pPYiCJIwJtYahvu6diCqhO/61EkCQkO7wg/TqGtVYXJFFreQiZVTSi0bszrXdpOSQncirblFmnsQ+6jopG29X1FTVXIS07I+fZUNzADhmSG9l3WKJRYv31EcWRWknVB9BFgOsiEhyMNZykxxM2N0cENAeDXnjpBW7WbYpzQ27qpW0xNLSpwojrnoEl981V9vOmUcehHOPG2JuBAF3KO9jYOLKmQMReeatwMPOZNjWjcYVWdRKG3fa91cA9O8nK+FI2KNdt+N1lByy2K1k5kn8WciT9Rm63veW3ftx4k+fwsIy9kj/w3PLM19TadQccZC5j5SIaPm2ckT40XnHAEhOvGN/8AROvy4QB19dFfuhZ/EFt3sruamV2hvxyYOEejaYBUvx/ldv3Zc6r9teUkZ7DdJy2ZOGjZgIyY/btPDJwYRy1dZiCSf9dDoef2NdpkW5JWMEe6Zkehma/vBNzxnP8/uvlDllx944nbha3Zoua69HlapY/j7E34V8rmI2qmeXbMKGnftx67PZE1D+n+Sq2xmoPeIgvfgUF2s4BwTEgS/4ogvcnpZitECyRP10G7q122pziDhOi1qpnbO7uS0wETUadOyAfSGyXR+1oymfsWgj/jZ3dSJJm+ud6fagyLLQysGE8vvZsrsF63c24/sPLci0oGTdQayzUmIUcpTJU8qGNdtjBsLVo87d5qC5XvGoZddVWa1ke4b7WovG8xxcpbxUIVXb0MPCmFUDtUccbOctEyOfixd8F/WQynilN0ibx+ZqkP7F42+aKzgi4LKzp2DmaKxzm15vrhc2hBFu7hv3vIbv3DcPd7ywMipT2yHSSC4WCZLvtNgSpYMJXdfdm6arPeE42koMN065Ley9AAAgAElEQVR7Cz97dJHQtlvjNntSJX386kLJr1Jqpf2WRVU1n8T3bXMdV0F1jUlyqMvnrBK87E2k3bckfEBvvLMTa7aZPZDkZ9mj3hOHqkOeRCmDtINaiXMELlya62KmqtvcWsSp187Ac0sCXXXRweaws7kVK7dUxhXOdneqfEoi6hzVSj9WZOgUsU3MDOu4+JkkhwTB07ZHCslB6gP2hUvnLfOr6Uvw+2djvbJpKonDdd7zowLI5zghLY+Z4fjQr57DXbNWSkyA27XuDgHqche1Uvpai81Bbk9zM2IwarOjtBH10QUs1LVHHOQC6S2IE0XFYwbEIV030YdFFNZKDtLv5Zv24J3t+zD1kYDDjCUH/eQ9/+YXtOcqjXJdHP+1VO/ZpLo10eNL1aPyGWvOE9zUNGqbQ9zQ8k27cd0Tb4VtknaxU9l/2hRU1fQkxXloJdgVlRzIuEi59rVo3U5c8+CCRJnqef300bTEawu0jM5p7XhmaQRIcv5tpXQgngxZmtYRZHH6ZN1ZsCvke6o94iDbHOTziXD99AvK5wTJQfN12NIE6N67Lf8O11WbPsolFdik3NWP26Y6L2eCq64QPyx1m/YFgEP+sKe/qTZeE4Dzxw9LlInP/czrn8H9r8SuvnppMF3WpjBImxY+lQG9LtyM6GuTRkv9qRfT75wzJnNCRu6tVa53nQnOhmbNcaqeA6PGkVIrSc/MNjb5MeokaHGu2TwI5W+uK2QRrz3iIOeLl88LBaoXGvi/U6puFuhe/AvLkhy1zGy4SA6VQCJtualeB6g4lGk4HLyL0u0Ix/L9CO3d/ZKYIzKJL31gVOI3X1TkDZCIgqBCFVQMhByhG4xRD/Ed8HvhYzEZzeVn9fvPjXceIxC4dRKZYgjKn4eqZ6CCODazK6uu3G5zEDewKjF7yg5ZotSplcRqKobABC85dAakZ140KM5VLzSfo8jTIs6SKhEcZl5cs6pjeHMuNgcd/vj8imgTc9ePmix92Vqp1Py2SQ5q1Z36mCjJpemIHxFA0tfBibOc8pmx5HaqYopw1cKhinMwPatd+9siqZHX4wuYvFCZOM5JRw3CGWPSCSs37krHO/RsKODa84/VSpFLN+7GEwv0LsM27Gtx08En7sdwb9qtdxVlMnHYuS92sS0xZiSwgLtaSXw3NvfyjbuSe4d7yaETID9zWToQJ4ZKN5wXDNLFiDjo+8sSBGe7NtrhrYyJ86OHF+KjNz2PZZt249v3vp69AQXsxlG3dv7tlBHRseqSrcLWmc6SgyabZyojrOZ6UiyLfIG2fegDetZHx6pnpFQrGV7qruY2fPtvwTuTm0vZRRISUww+l1QbQamI1c/PPzbKMaRiJs664RnM0KRfccFeR+LwwrI4+2o5EdJKbyWmJw6MJZ+hS8ZWneSVS6iV9Ezoj/65ELc9t8KpzWqi5oiDyRgFJCeZ6iMmiiNn+XlxAn33vtcTH5uK88oqMu5rLWL3/rZIysk6bfjH3VIs4eoH5uOBV9/JfK36nPVqpz5E1YjNIO3ai1FycFC9i+pDIFhUo/ctzRk5JYvNIUGpVrI8qgdfC+xR8gIpq5V0NgeOrHELJrVSe+DqvfPrGXF+L9Mz0koOimL5m+dxPUDwfMWzSkYgpX2wPyGVirpUYtjZ3Krc592rlToB8jOXv9O7Zr0tnFOrlbj/Pp/gYr17565JbBfZLskhJCyrtu7F2O8/UbbNQazep5t+/2FVfd2837iz2ei7XZcnq6urqg+1E4BY1+3eTfWcvJWQfHenjBqA1vCG5MU9RRyEe3DZT0CF75wzRllukxy0hvjwrlU2aaPkqz8V4ebPnOBQK8BBvRoAABdNPMT5Gg7VMG/hfWvuwSXOQazSWkwapF3ye+m+EbGfFYotZH87cymO/cGTmnGr26wmao84CMdnXPe0kmpzqKh9joCGQhD5y3XPWT+ucpmCOPFetuvECd5LE3m5fkes8xTb14m3E3863bhnciGXc7atzH073rpDKTlYJAurK6ucMsWReRalvkKOoudvy41kU0s0t9olhytOG5WqA2Sz86jUSspEkJZWbfN10lEHmSsI6NlQQCFH6F7vFj2fHEd6IENC1VcWyUFWCDAWP5+mXg2pnGm2Nl1iXL57/7xUjqXpBrWcVyt1AsQJtkKT4pnj2GF9UmU5leRgdEN0/xgP6d9dujZ5nqs1Mhu0Heqc9LPpcf3wfoio7F2yggRmbnXnrdkRPUvVhyYGE7mORkyb8MArsRrNNSkhUdLHP5/LRTYfG+cv9nCkYv9qlUrFRsDi2JqgXkOYWnrnvqTUYtt0KEuWYD4Q2xOrU+TC0Em3RcbKVlWprrF5DrqolRgYxgzqhb/8+4m4+kNH4rCBPYW6iuulRv+uUdPK9eRcVn2766X4rrBdaA0SB/e6clZOIFioeM4glVpJRhbJIeVmK9XjRq3m1pLzhiB5KVeMy+3zOoyxspP4BWkI3K/lUdCq7m4WUmw7J+PT3GlrseQcxyGuo3V5igzS1qyqQtf5fLovJXGQhiszFVESvLDeuOF9AQBD+iY3INK9Lt5a1shbW/WcIs2IaRyMmYMGTVB7DvH+7Nw7h8ggrNuxD08s2IBVW/filFED0auxDnddNhGnHRF4damYI7nkl5q9YWzTpF/3eu25rrBtaM0RBxdE23EqFpgcCZJDpFbKpldyFYFNus3XpR2zACg3kcnnkh+ii70iaXMojzjkc9nYw6KjPcXZW0lTr7Wt5GyQTqiV8rlojFl81lXjULqyWtqJOP6w4gePORgPXnkqLkgF6omMgHr+Zu3b9E50KVJs30Q5nLHqGumxKK5JQ3xGJ/8s2NdE9J4a2LMBZ4WqMhe1UqtmPtikbpPk0BVQc8TBZXHhddT6b6QipLNy1661TQZd1a5kZ93wTKrs2KF9yljgzTpXF7hktxTB79XenZ2TA4AnFqxXXt1aYm7EAUm1UiFHkUHaupjadPiKMntULs9yFKr8ABw3vC+ICH/70slx20I7TwpxCLHNQTEeo1rUOCwtcdDnOmIBya2U5MCzP2WQHFxApG9XLqlTSIeAOYYKcHMO6UzUHnFwmJVM+i9CTLwXRatm9Ixx1Y/eMG2xtl3TlpUiejQUnJOXcfD6MxdvKttrIp+366pF/HJ6IJrbo67TZapLVHt3A4HkkLxY3xd/d4c19UA+Ryhym4/lxmznXQP5RMhqJXGRF/c85uc3796Pr/711ag89lZKz0hTyhUbHS1oFkaTdFy2zcGoVrJfM6xfaLx2mNT8eSu/balMRyC/eY85nihrOpNqw0ociOiPRLSRiN4Qyn5ARO9Ie0rzc1cT0VIiWkxE5wjl44lofnjuJgpJMxE1ENE9YflsIhpR2VtMwvQR8skTqZVYmsusL+TiCVmyt6k0SDtyM08t0nszNDimw5YjPrm/vAmVcLEu5Oypj0Vwo7F1X+N2qpUO7tMIccnTMQtEgSrptksm4O7LT4KoJbOqvhLjSNY9fUwTlqu8vCz3JRukxXmVjCAP/uskPnk+rty8B1/808vGvk2367qhU2IMoLJSb6jeFV9gU1HPza14Z/u+qJ/rLzwOT33rNABu3n7cCUKpVpJ+H3Vw71Qd0fuvHBwoWVnvADBZUX4jY2xc+PcoABDR0QCmADgmvOZmIuIs7i0ALgcwOvzjbV4GYBtj7HAANwL4eZn34gTTvIjUSQBeWLoZv56xNPVhjBzQQ5AcAmxSpB/g2LGvNfUhLN6wS1nX5YMZ1LvBWkdGVu7fJF2dfNgAjD+0n7UNIncimOjbKjmoPla3fqZ/+zSMHdon8eHpDH+89OyjB+GgXo3IUawmy3JX8nDr8zlldLDtHuRd2RL3IBzzMcrxDzq10oad5kWMiIxj021KY5UcKqRWKoTESY5A/tAvn8Op186IRt6tPh95eLmomriRXaUZki+fs3JrKubn63e/ChVc1bSVYNDaCytxYIw9C2CrrV6I8wDczRjbzxhbAWApgIlENBhAb8bYLBZ8/XcB+LhwzZ3h8X0AJpGK3a4QXPSrjAH3CRk3RYieGXySnfub5419ymL7K2+njcmA26Jz1OCAS3HdTEyWHFxgqh4YuO3tkaUdHSonOaQrjmoKXBTFyWWSHJK/Y9dcu1oprrBaWjR00d6qNu/9YmxLkF02xVaSiQmZsZ8suZjkflTQxSuY2iXo5zr3wlJBdQ3X98tG4XfCRIiRpAWzHUEGF4hUaiUVcfnJw4sSvzfvVjOM8wRHkq5AAExoj83hK0Q0L1Q7cVZyKAAxzeWasGxoeCyXJ65hjLUB2AFgQDvGZYTL+wiMZvFnsehHScEpIhBhY7YPTJ5MbaUSRgzojt99Npkl02WycDHaOVK4BCzLuE2hqWVXfbHIaWfq26qyUemAnYoyQeZPcoIkZOPyxflw79wkk6HbHlPV4sSR/SOXylwkOXDqYFYr6Z6jTDOyeq/JkHcsY5b+gzHoXVk/cEQTThml/vxVbXJ9v869mHP+/H3myE2SzpnUSorrZQlqZ3NbuhKS362RgB4gaiUVbgEwCsA4AOsAXB+WK+2vhnLTNSkQ0eVENJeI5m7atCnbiHnDDhODi74c3STuiJ9yXfzS+W+CyScH2bmoR2x7ScgoMYav/fU1p7rROCwftuvHpbbluS+s6rHZ+7ahnA8vl0Fy4NjV3Joq0wnFtsVc3ppWnFJiizrV1z/nrVP2b3O3tD0rne1LbyBmoeSgrpAjfZ+qoXLioHMvjry7oudoVpNx5DVM2P62Ihau25mq36sxSRxU714eZ1fIn2RCWcSBMbaBMVZkjJUA3AZgYnhqDYDhQtVhANaG5cMU5YlriKgAoA80aizG2K2MsQmMsQlNTenUw46j154hzQcnQ7Y52CETB5ZK7AY4Sg55PmndemYIJnQWWFzUnQZKpJ781sW/jFQOqivM92CnDnKNHOkXXt14VLYFrVpJ0xZ/XvLGO+I9qCWHZDsrQunx6MFJ46mLDtz0LDnjw4cQq2bNk0h3mkDaZ6TsP1IrqSWH3z8TbMWaiyQHN+ZGZ5C+cZp6f3CZgVSlSJHba19Sy45HWcQhtCFwfAIA92R6CMCU0ANpJALD8xzG2DoAu4jopNCecAmAB4VrLg2PLwAwg3XgbjYuiwZjep0hYI/KlJFO9MWQo/QWjK7qmix9M8acc+fH49C37RrbRpq0C6bUxUD7iUfWelqobA6uiQ8jb7f0KZ33ot69OTiRl9RK4txJbKPK1GPkAZufOfEQfO/DR1n75SCkOe2756xK9Z2F0TEzXibpKl3GvaV+8sii9EkAD4W7K+YE4uXy7ZCGOOgM+HlHgpbYdc5QrytIFWpXAwFE9FcApwMYSERrAHwfwOlENA7B/a0E8EUAYIwtIKJ7ASwE0AbgSsYYX5muQOD51A3AY+EfANwO4E9EtBSBxDClEjemg8sjv/+VNZi5OKm2evwb749EQtm3ultdHvsMKYjlBbFY0hCHDDYH15xHJQbj2BJ1Swy5nDm1ATnaEgKDtEJysBjSV1ryXbnGOZjvwTwGIC1diGoy293zMaqek15yCOp+46zROH1MnMgudl0N60WSQ2KwcTuaMfJ8UESEsUNjdabtXcqeRYwxXPXA/MR5IFh8i0L/ps13TP4mqlQc8bXpNl1jBUTJxuXTkeNK5HYA4LnvnoH3/2ImAD3RlyGm7jA9+84nDQ7EgTF2kaL4dkP9qQCmKsrnAhirKG8GcKFtHJWCy6Lx4vJ4u07OmRwp+DLnpC/1iIN74fXVag8kIK0PLbHgI8i66TggGCadiQNzVkEVGUMOZrE75/hx5XLqejaidvldc43nXe+b11LuaeBwvSrxXcmROpjUT9pXHlYe2LMh4bHDOVdKVkvMHXFhincnTDZ/xelxlldxDIsU+nN5vCJBl72C+EIf/BcXPnV7TgZrTbnqEtH7an9bMcqYnG5TUCs5fBDcC6qlmGSsRKZhuJAo09XBslSyPyPAfZ53JHyEtAWyLhEQDdJhm5YXqdptLkeqj8A+tjjox1o17MutHiAuLKaL3HS2BLWEYZMMWi2Nq84adwhTnCrHUzqXEwzSlvfEEzOqFiFdTAyvKQ9NbiJqU6NW0o3xQiEHk9jFdU+qk8aJbYttyfarSF0jjd/0jEweb0EGAvU58Zph/brhqW99IPEuF69Xxw8BcRbbnEbdKaN3mNpCznqrG5sro1d0lRw6nzbUIHEwcsVpQ7MqTUXscx7UtInmbdJKzr2VTHraHXvV3g75nHsgjzhGFT523JDkuHjEt6E9U3BbU684QE9nm7jwd7MMrdvTYb/y9rZUfVUabpt6Q64nQ/7URV217dGb+nxuyWZlucrQHJzg/SfLxV+qIDj5kai2B3VBcN/xb3n/bNnmEEsG6vZ4biXd+ZzCUYNDzFU0ZlAvHH5QMh26yWW7MWTyXG0OfUPisGNf8jvU0QDXx5uwYZgkB035g6+9o10bKg1PHASo3m+jwlVPtjnY9OiyUasUSg4p4iAcy8FTHDw4x12tpD83cWT/5Dgd1Caqj+ATxw9NlelcXl3tHzrIhscj/vsxfP5/X1LU1Kt2XJBWK5F14QOCCHaXejJUhmYg9sZJ2Rw0cQ46g3Q5Kkx+ndiWTBx0kde6BfiTJwwzRl0T6e0OogeQSvrL64JIADSG6iadi7WMnqFr6q79suSgHps8ZrGaSJhNNod+3euinfJUz2/F5j34+t2v4ev3qKOvK43aIw4Zl4tGheRAkoRh40RkVUmJhW0IE+iiiYckPkKdax53aXRWKxnu91xJcojSZptUAgp10SdPCIiDWGySMKoBk+SQcFnWrJlpg7SbK+uRB/d2lihV45WHc92FxwGII4fFrKwcjXV53Pq58WGfyfY4xAXKRa126IBYn56QHCTiTpLkwKG69bFDe+M7HxxjlBz2t5a0XLi4PavqFkxZULtllByi9y/V1UsOyRNNPWMpWvRksgXB8dQ0qpxV3OuwvXmbXFF7xCHjevXri45Plcm+3G8adJ1AWq3EbQ7iKlMnZTHV7VimC87RwSTVyBOdj9NokFYYmuPJH59w5dB0OPLg9A5qWcC7VtocxMR7mjGKTglAMuLbRPQSCfpcByvUld/JkL7dMLRvt2jM/N3ITDInHrp5kSUD6PP/eUbk0WSTHKJFUfAG0o2jR30BuRxFNgfVc3zg1XdSKjBeb78gOSidHXjad8UE5uph18j92HYil2skB6l4QEgcBvSoT7wrUxAcA3DDp8ZhQI/6VFBdcL66zFbNEQcjFO+d5+MRIepXZy3bkjovQ96APlArUSK6VOamdJJDZuJgXMiSN9zstCc2YcuelmQ7ioWHMowRAI4YlHzOuo/QNVV5rF5RnHP4yOZKto1EbiXDdaKrbzbJgUsE6vvmY9bWk1Sdcte5hORgHsuwfrHUINuOZOLAbUAuNoe4Xz3j0FYqpd49r9tiEZe55NCsCPrk6mFXb7uI3ZGfo+bZqQjaxJH9MeM/Tk9IDqJbu2oc3erzOHpIb0tqjfJUhFlRc8TB6XsV6ugWPiD4+LfvbUmdlzFt4YbEbx7n0FiXx3FhCg2SODRd5OrwfuYN1WVksbHw7St1i2c+dLHaKRnpVPv4qiQH0/aaY4ckU4nopv9l7xupbUNEzL2n78VlJzeevp0jTl/BLE4NdqOsClFdxY3ncohuSCdhyO7V8vwQF6gsS4sc1yKrlXiCO/kzUd07H2PwT/0cS6V0QBnvX+5bBp9equjkWK3k6G0nOZ1wuNpuGAtsCH261SWuaRHmnk4C1XlUVVtLa41zeLfBtKi6fjSiV5PLXHl43jr85mJxDLFa4IEvn4oSY5j6yKKkWkmYRDO+fRqG9euOeWu246BejQAWVMTmIIOL7bpHlA89rAp5Sqi9IiO5WFmh2/3F428q2x3SpzFVV2db1G2sIiO2OaTPuQQQfv9jxyTHQ6JaxzCHRAkjw9fMuWLVdBLtPHovn6TaSa6WkFIzUAeZ035VE88jc7MmWw+XklW3wp01RPB6SckhfTWXHFROD1yHLyZQNCGtKE2Wy0g7l7DYk0u4IXGzqTtnvZ245uTDBjiNsTpyQy1KDoZzruKaqFctR8Tj6TOAgBuvy4cbCAmDExffw5p6or6Qw4QR/Z1Td/CPIcteDlwcNyWBI0WbsmtvUJZ+1rOWq1VwDXX5VJs69UrWtALKbR4d2pBpkLjZjk0aY4xh+abduP35Fc7j5ITzlVXpxVeMC3h+yaZwPGmDOR8fkLzHC8cPQ/f6cvnAJKd97WNqAu/irfSvpcH7p5DgqN4Dl6rltl5+exv27FdnOhWvBaBMFxMF62nib9L1g/9yVd33ropPUXlymaTnGz89LqyvHmPsqaYfdyVRe8Qhi8JRd1qYOGJVWW8+YkB3nHXUoNQeszzOIdll8rfO8yIO6zffx5fPGIWPHjs402LK1UpPCmqwi088JDou5IIAJdngp4oPkQ2ZgDoRHaD2INFJPM4R0ly9oniMLsnmUnEFEWcej+ysowalruPqtM/8YTbufml16rwOq7YG6hlVTi/RHvWDfy4Mx6Mer0pyOHKwvFOZ++riymnbJJdk73pX1v496lPDW7B2J86/5QX8WNgzQTUkzlA1G9RPcsru+oJ6CYwdAJId6RZmeTxMYBxFF1sTceCekUTqeatzd+4o1B5xqEAb4qbm4iL/5DdPS/U1qHdDaiNxnpXVNDadt1KcZ166VpqdbUUGU958VWpkrlYSA7VEQhB4mlBKLaPKQ0OUljB0CQBV9gnduF0locMPCgj17c8vL6sN+fWIdhU+NpVxkhM62XCrgsg08GPZs433nVZvqCUHlSG+PS7FnMvXYeKI/lG9ZP8m1ZterfSrKbF3IN9IiEeVi44dYsAlh4sjgPxN6FJpkGZldM2qKzKOohTa4mDv0iWtjCSHKimWao84KJ76hNC3+D1D+6RPKiDm1jdR8Za2UmqR3NdSxKurtqcWSjkuQMfdciZE/gDktAw79rVGXJ8Y7Pbts4+IjmX9vYrjEvvJ50jpo65WK7lLDjly3z7RJglx//xxw4N3+vqaHZnbCMakV9twDk61UPCF3MWrShwGdzVVMgUKyUpnkFapldpjyLTtf3DzZ08Ih5hkEFx2glONS1z0h4deU3K9H557DK752NGpa/mzs7tixxV0cyGyOaTmurpdea4zxNKcaGB3k8I06lA+Ni85dBTSD/1jxw3By/99FiaM6KerkoAowpu8F5pbi5A3U//z7MAIlXKVlLrVedToUnb/19/nJ353b8iD564XuSPRfzpFHBQugBeMj7fnyBOl1C2HH9RTm5pB/khVXPHDX32fUsd6zJD0pu1ym3IsAhA8x96Nhag91dDE5yH2+lchFbV+8Y0XjPPHD0MK4UKuysn1kWMHJ36L98yDG1VqB5IHCr1kwy8Xq8uLu2lx+ekn3pNqV7fYDuxZj4GhP3/a5pGuz+cJ596VWVYFaUqWRjg+c+IhShuKLqW6OD/luaa7tzjQVa1C5Xj4q+9TXi8a1sVvxoUJ0tscQrWStYXKoOaIg27SDujZ4CyuxbuEMuOb2tdaTBlmdb7asvium0ScC5HFYdl975tnHRGpOBKeRRJBEDk1lQvgxJH98ZcvnAgAUQCTDGnX1LAsKRrPWrZFee9jh/ZRLkCfPEGx8CL+QDbsbMaUW19MnSeiMO24nrvXSfZXi6mopRcrOgLwy0cO7A4ZudB6rJIcTh01MHkvwnEhUiulB6dKNyETab4AcltVUq0ktZfqIcBvLz4hYWPilZ0krSj+BuF/8zUMDKu3plPEiMF6sh0FCOZaQeOxppMcGgrJeKI47Y1B/cTHmdbnJTBiYI+wz7jiuh378PaWvVE/InFy8ZRTqWQBYb74OIeOgerVxByNWxuuksOQvt0gbhID6PWF8gLAJ/pvLz4hUU82/C3ftBs79rXi+aXJhG6NdflQakkSGjlS9qlvnRZxPzo/cv4BBK6s0riB1J7awf3EC/kNTy7GRbe9mEr3HN1TKOr/dubSRLsq8I9Qa9wGsH1vK+6c9TZ2NbcqvUvED1nXj3yfkeqsJHKm6aspbF8VH3P0kN74vqAOET12Du7TCAAYGT5reSwqg6eIumhHtJB7Fl5GyhMsw+LCiR0AzHxzo7aezOUbOd9QTP7544tTdQLVZfg98uuEe9Gl5AaE9C9S36LRWeTKzbYJ3ncSaUeSsJ5Q8ZM3vwAAeHUVDxCMz4n1dKpHOeaJg68j1ZIcai7OQTUf+MtLSAQWUKh8172o2y+dgLFD++CWp5c5GcFlXT7nAE88LJkcjySbw5nXP4PDFAsKEOsuTcShT7c61DUF1zdrjKj8evHDFcGlmZTkEBbcNGNp6hpIdUuM4X+eiBcL3eJvk8rFb/exN9anXFIBvVop2ZA8Rt6/QFg0UhRLXw4gSHGxtyXpjsntVscf0g8PvPIOvvXBI1LXqdwv5XxdRIR8jqIIXNHbpT2753FiBwCfv0OV4DCA7LFmYpD5+1blCCoInj38UJy/Ou8isZ7ctdim6BZstoskbShxufRbIiLz1+zAuvC++LUJyUHotG/3OuzbkZ7ndTlSStnR3h5Vog41SBzSMyKK3MxAk7l4qpMcJoVujrInjvbFSuonLjnIi3mkVhIaXa7ZI4Gra8QJqbIPxCoJ9dcS6e9zKhfK+HzPhkLkUZLYHMcCVdSqTv3margGQqnGIjnor01eF6tN4jgH1avkEcU6w6XcbryYBkeqhGsqjyGVq3MhR5Fa6pVVsU1rYI+kZ0+WtcW8LWxSlw+YbQ5imyWmfg+yfUBuq8FAHKJFX5oj4iMVHSXcJAfzXJGJCI8YB4Bdza1h//E9Je0d6rZ7d6vDruZ0TEeWmKVKwKuVRESisb0d7sXxy6eSm6X07Z50W5W9kPjiLqe5Jok6iNy63G9w3j7G2OYQV1YRB5W3kYhoLAqDNBCnJvjECb0BQyMAACAASURBVOI9uaUpCPpPfyhnHzUIN110PO794smJcpu3x7nHxWNoKzE89sb6VB2XcaXUZ6JBOvI3V6iVwoXcNc1CpOLQvG8OecgqFV0hF0Sub9vTgv/+R7Ct+7nHDcEFkuE8C+dp2hZWfBe8TT5XjDm9coE7tO07U6mVVJLD/1xwbKJPXZAmP+afg046FSGPUXYYkImI6J7MGbxrPhpH27vsBNersYBdzW2puR5tGWsddWVgJQ5E9Eci2khEbwhl/YloGhEtCf/3E85dTURLiWgxEZ0jlI8novnhuZso/LKIqIGI7gnLZxPRiMreYhKqCWniBHXgRiPZ62jmt0/HtG9+IPqtM0j/7JNJr5BA3I1rxpJD8hXpXFk5+navi2wInCMXJ6Eq/YRsx+DghuhDB/TAuccNwW8uPiEtOYAwuE83PPfdM/CdD47BJScfit99drxz8BTvP+UymCOce9yQVI4j28Iubof5qsA9H9YUq95ctomUF35VbiWl5BDaebR5/6Vi2fVTRXBUz+eUUQNS9Qr5HNqKpUSMxUePHWzcl9kGldTCIS6UqcR7ivpMqMuYnifnwaRD+nZLtAmoJYcLJwxHY11O6DvZsvguOMO0t6UNJ/x4mmYEegKq24yK9y0ayzmhjLwgkZy/uu+jV2MdiiWWSgMSq5WqQx5cJIc7AEyWyq4CMJ0xNhrA9PA3iOhoAFMAHBNeczMRcavLLQAuBzA6/ONtXgZgG2PscAA3Avh5uTfjAtWU5B4E8X649nZ0AWb9etRj9KBeiXpy4jKi9CRP2xzMkoNugTt2WF8h3TIAuEgOwf8SY3hs/rqo/JTDB0bX3HTR8aFnkbJbDO/fHbkc4UfnjcXksQdbg6fk/nXqoh6Sy6JNJSSOT7xX8bIjB9vTgesM0qIrqzxP7r78pIgo6uaQLgeRyfVWlD4H92nEpyYMS2RO5ajLE1pLSZWWKs+QzqirerSc2KkgesHxPlujtO/695QnSrxvmdD9+/sPw31fOhnvH90EIDk3dGMX1UVpRkO4n5Bh2q1Q24gQA11FyEFsMhNQUMw5sY7oraTjUbhqUSZE0TplHHnlYCUOjLFnAWyVis8DcGd4fCeAjwvldzPG9jPGVgBYCmAiEQ0G0JsxNosFT/su6Rre1n0AJlFHkkbFC5E3cXdBsJi76K6Tk3Xjrv1oLOSt1J/rjmWbg47Lj8/Hx6u27sXm3S2JSabK6y+qTKwbzju7+5qDp+T+dYuJnNfexvTLKgSOUYLk8M2zjrAGqcmvJxkEF9YBJbZaHdKnW2Tn0b3eY4f1wQmH9E3FPKgWEnEskT6dpXMPcWze3YK/zF6F9TtjQ6/quY7JsFeGSQIU3W75vIg9huJ6PRsKiTK+H7eo+//bl07Gc989Izo/YUT/RLAph84gLbpDmzLScoYpS4ocEfKCHbk4K5g5VbqLpFpJ/Vz5Pcp2N/68ZW1FR6Fcm8Mgxtg6AAj/HxSWDwUgJpRZE5YNDY/l8sQ1jLE2ADsApGXmCkH1OiJdXgbq4LpRuRgKP2fFVtz90molNycuAEBgcCRCSiWg8pqRx8XxWpjEbcPOOHpap9Pm3lc2P2yVQVpXz9XmkNdIYUD6/vmCsk2RKv3Z75yRqC/e639/JHYhLeRzOOHQvpZRJftN5FYS5suvLzoeQ0P1B1GsthCJ6HUXHoelUz8EIHAxfuDLp0bR+GnJIf1A97UU8faWPWE9u1ph4dqYwOtsUx84oilKT8GhegWmILhWQSIdFLri8ud//i2BO+fnTjoU//v594ZtIfpfKsXquRID3juiP4b3T0pDKqOwdv5CbwxP2RwYLMZH0d4RY+ueFuwPA0Uf+sqpynqqADtxxInzwoP9479NiI655CDblWz7q1calTZIq94cM5Sbrkk3TnQ5Ec0lormbNm0qa4CqRagcFzFZFfKL0CgmQ1QJiB4kqXpSJHVbiWm5fNFDSDUuju995KjUedmGEV+Xtk/o6rlA56vNcfqYJvwuTL2QyyXvhy+2HIPDhQeI612mcKs8ZEBycRHH2lAnq/HM96FTK7UVGf5v9qrEuSjgLsw9xZBUZeRz6cAtmSuObQ7psSzbtAfLNu2J+rKZEBL7FOtyB0H9LaTqGeaaWPyrMKPoCOkd/McHx2BAj/qwrWDg+RxpF8lk30k7BgC8rOGauRqzVGKpRVS0wXDnB5kJmnTkQYnfct9rt+/DCT+ehofnrcOYQb1w7LC+ynrimqAKxBSJtTiEM4+MkzjWFYL6Ynrvx+avU+YJ60iUSxw2hKoihP95dMwaAMOFesMArA3LhynKE9cQUQFAH6TVWAAAxtitjLEJjLEJTU1NZQ1caXOI1Erc5mBfAPO5pN5UVFuI4B/hH55bjn+8+o62vbTkwLRcUt6ghjn+kNj4pdr/mi9aqsWqxJiTsdYFtm1C7/j8REweOziqK97PL6eM017Hh7dtb6u2TlxX1Imrn6U+PblaYvvby6vxerinQaRSYHEdoiDPldifKsOmnPXTJDkcJWRVNblPnxkucKKhWJs7iNLfgupZ6JLAAUmnin496tGroZCSVHo1FqLxiukkiizuXyetRp55DlSMR8X/+11z8YW75ibOibEl3OYgqsSmfmIsbv+39yauiSWCoN4qIZq7Xoq4FusVE2oj3mfcrktep7pIcogf5hV/fkWZzr0jUS5xeAjApeHxpQAeFMqnhB5IIxEYnueEqqddRHRSaE+4RLqGt3UBgBmsA3emV7V8dJjHh79EeScqFQr5nGREVV/D1U8/eWSRca9pmZMLJAc9l69TF1xxWuyto/JSyRPhm2cdgX98+dRk/5HkkE2tpK2H8uMcTPYAUa1jg8hBuwxb3LdaJznsa0l76IgLO2cAEpltlYtu8J+fibnM9LhOPmwAeoV6+2JJLzl85czDAcjEQV2XoCZaqnHqnvdFE5OpNnK5NNOSy1EqjihPSa8v0xhN5+W6JQZMV0Rxiyli+P3YGAfZ5pA0iqc3TlLVi+sIkoODQZoTB9u2qB0NF1fWvwKYBWAMEa0hossAXAvgbCJaAuDs8DcYYwsA3AtgIYDHAVzJGOMK9isA/AGBkXoZgMfC8tsBDCCipQC+hdDzqaMgvo+eDQXM+d4knCLlvNFx7HIdXVptESaxPFUxxJ9nv43nl2zWjkPOLJk8Jxvf5G4IXz9rdEQQxbpBNHXwW7chfSp9gMFl05XCy0ZPEwFiLMnFmVAUdbZa24hOhy3fZ/BftaOaaLNS+c7rdPkAwMLnzZ+7aqHK58S8QfoNpvi1oq5ax3WrUq+r21QzVN86Ox3JXRAitJN98QNE41uyYXe0+OnnMoznk+PUO0Dc+f8mJuuxJOOgYgbjxHvAC0s347E3Yi++ZHLAuB6gT7X/+8+NBwA8+NraqExrkNZ4K1Ub1ghpxthFmlOTNPWnApiqKJ8LYKyivBnAhbZxVAryIhRsu8l/c/HXgTgQSVGq+o/QiTYI4/ve34OQEq6rlcHVMPJHw/cwEMfoCt4mn7CFvG7RdIPJA0ndt9iHvhddm985Z0yqLOlqqW5TXMzFptPeSkFBo+BKycdZF21B6a66kgmLySkin8tFCzljeuaFF4uLik5NGEifDsxNmLpj6cbdifKvTRqdbjOnlmg5w8LHN/+dZBp1m3OFLidXYpwGw/ngPrENK7I5lJLrgBaM4eI/zE4Uvbg8qfXmzhyA/nmfc8zBqqaVqNMYpKuN2kufIRzrkmi5bFMsSw66F+2+mKbb0QUvcbWS3OeVZ4xK/M6aYE1UK9XpVFqOAVXcQKiC7CUjR0jrhl3I6RcAZfZVgXjLZ3l1m+uuXJ8bC8Wy//vCifjn62u1xFz1HNK5iFjo7aSWHMTIY90riBfT+L6zpBtRtxmMcelGvUo0GieRcnGMBQf1wHXaE/4sdLsiJus6upaHTItoc9ARddkOqG0TsYrIRZvAobc5pN9jZ6DmiIP4ttMpEtybKeTJKYGbu44+uVgAJtUOElw+R15a0FUcpm483EYQBd85Sg6629MFT7149SQM6i3l+nHkYgt5vTSiIugmyUHVjOldqfJv8aNRTT3xjbPSahYO1a2lJQf9IpXPBfYtris3LWaAm0H6qUUblOXpNoN5rrN/JccZq6oKOcLHwxQxskFahsmjCnBbcHXqL1U9l3xjvH9Vm+8fnVRDixmVn17s5kV5/8trtPeli3OwjbfSqDniYPKZjj1I7O24Sg7urp+8Hfuk5YY/uUtZjaS0ORj6Zyy+d+1i4DgvdcFTBwtuqcm6+raYMCbGmJIDVT1ncZHMQvhVUMaXOLapzHYquUCaJII42aI5zoE/g4RaqZ2+HXxe6JiFRP+5eKEv5CmSpOL4BnUbNrWSiwedboMcGQHTUkpsOaplmDR2jLsEG0bQZvweF2+wS6JtxRK+/bfXtecjtZLG5qBKztgRqL3EewnJQb2YukzGguTKqnWJdByXqp5OcuCurC+tSOo+5TmTJacOdwXkIry7QVrdHlGGIDjJy8XYpkbSVgaPtYpqJbdx60CKhco1WlwZrxL+Z4LkoFv0uf2HSw82tdL2fbGbb3s1E3kibcyNri4QuIqKO78FB+rrzj56kLI8Tsnhsug7JlTMBWvARbfFG0VpJTHoJMz0XOLVxg23BVemnQQe+VpyNzkuOejuu72qQlfUHHEoGYiDLOqbkM/lkpKDpl7WRUhsZ+UWtVcOhTYH2VCWJnbuaiVuc+D3pFsLXG/HlZMLxkRYIhg8TVKYS/AfR7NgbNZtGK+D3E0UBOdqyAzxobEH4xPHp3e1U9kcdA4EfJGds2IrWot2tdLm3XFEfBbJQVW1oS6HlmLJzUkjyrbKEgQlZ6ANH37PwfiPD6adCYD4fpwM5wrHj8a6HH5z8fGJMtUcMmUNcLM5xH2biPGXThuF+kIuNc5jhiT3rlfFOYhodfFBrgBqjjiIHH56TrirldKSg7qeq1FYZZDWQaeySSfpU/akbVM01OmG4S45uKUXUbWpr6d/N/L2pwCwt1VQHUjn5Ky4Ln0DSLhq2kY9sGcDbvnseGU+IJW3ku4xcKLx2dtnJ8aiG+POffF9f1DDlaugUqE0FvIolhiahZQvYwapczPlc4F9IlJNhu8k8gJUDHxAjwathBsZpF1cWRXu3Q9/9f346LFDEmUqryY9I+TmaRhk2w/fo2Hx4I4FInE6dIA6gSKgtjn0aiyAseoYq2uPOAjHp6dC5oP/bpKDe5yDC1QbmsuuqWLfqjHKH1kWV9YoclSxB3Gynmt7wce6wyGSWf445cju44YHnFX3+kIQPKVYxBpUxEGQHGQCNKRPN7l6Arq0z0nJwfwwzAbusB/By0UbES/bxrSSQ1DOF47FP5mcyBBcDnhch6ijP0hyKODg7rG8fz7uKBeV8hp93/HeJW62hFQAnqJtFWNl9FbKGAsiErIfnntMoh7PRivez42fTmcDqDdIDgN7Bs/eZS+K9qL2DNLhe7ntkgk4fUwyBQefIi7EIZAc4penC8Bx5Yrl8QHA+w4fqKyji5CWiUE2V9bgQ2hT7EEsQs4PpSNA/IPZuCvOEHpI/zSXFNRNtiHvo3zjp8fhrQ27cc2Db6DImFJ6UHGXCbWSNEzbo5Gbi1QciUykaQzsWY/Nu1u05+P2uJQa6+hVe20AaeJgswdFi3N7rfCICXXSgKsnYqLHGx9nnF5EpebUjzG2Odi5ZJVUqewPbvWCsWVRK8WSw7B+3fD8f56ZHmMuTex45LuIOikITkz/0adbXVTGjzsKNSc58Nc9tG+31McYGx3treRyhNY2YeroOG3HUanmpzY9sSZCWl5EVJyobptFrovlnLGOPv5r6ZboePyh/fCbi09Q1uOcnJj0jGfnTNUVhjnh0H6p893rCxg3vK8xxYeqbK+gCkkbpLMtnCrJQdWnKZAu0b9Uv7VY0jsgyMRBQ0RkA25Wl0elzSGcL3sE4qD1qgql6WijqnCccRZbSo3LxDxltTnI0dmqplWR1HrHvHTySJ6NVe5HlBxMOdGApKFZ9S65AwJn1I6+5onoXO+QIIjvo6NQc8SBv0TVxOFFLjrOQo6wvyhKDmo4q2GizUWSfaiQI3XqA533FQB8d/IY/PDcY3CMlDYj7j/gqPhiwPfANuGajx6dSrMs9r1tb2uC4+zdqOZ0dHswyIhz8iiIg2IB2WuQHEQ8Mm9dqkxujV//6uo4+ZlKrSiWmLyZ5J3TWisoObQVSwiSALZfcuCb6+zZHz/L4YqNhnj/xRJDWzHp8davRz2aejVEapb7vhRv/WokoFnUSgS0tEmLvkZSkZk/s1op/j24T2OUjVVEXT4XMUHFkt6xgEsOIhGrU7gIyzm7RPQO9zephlqp9ohD+F/1/mKKbRcd8jnC/lZ16gURrmol/iGJi/7bGm+lvMKoxseU6Fv4PbBnAy49ZYRRX11iDA2FHBrrcvjRecco6316Qpx018SZzg7dbG+aviQqk9NmR+MUm7HooAO1RfqcuID0CCOwRX9/02u48i+vYJ/lY+M7h70uEAfVYu7ikgvE+1E8uTDY47qtVFIuFEBaPWTckwNAq2PQmowjD04zDnwOiRKgzruoIKuVwvupy+fw0vfOijY4EsfmYnNoc3BlzSkkB3W99KKr2hsECF1Zhd+6hJDD+3eL8n2VGNMa2Pl7s90Pr6f6xrmaL0skdrmoPeLAJQfFKpQlp0khR4nFx3XXMxumC5GrTy9OZ5gEeAyBSq2U/C0SJpufei4HgAX3Pqxfdy0X+7mTDxX607fJM9Bu3BW7VTZqtnhMqhkMYwy5vp89tih1bmRTbLy/67ITldeacNQ1j+Od7fui37J00qr4GPsr0mWY4mhE8Odz+/MrAAQLhk5dJDejIyKizaEM2pBKxgjEhIlvcgMYiHzowdcm2RxSbTqqlWIPMReDdNo2oVqkVankRako0aZUt0FDHHo31mFvSxu27N6PR+ev1/I3/FmKXkh9u6fnkMkxpjF89tWIdag9g7Ri6z4OPpldsiHmc7kEN/XeEf2V9VxFe65++frdr0VlqgUJ4IukXa0kcpw2/XPAeTFs3dNiJCQi0XDRafcQ8ijpFjVR3WRSxXCiKKcouO9LJ2OC8Pz7dEtPa5e3sEvYV1h+us2K3ftUcE0nLkdcZ7E5yGlS4jZjzrQcyUHdd/A/uWGOXsIplhjue3mNcZyqrKYqxFuPukgElGLqlBkCKC0p6PT3geQQtzm8n9rDLZ8jNLcyXPPQAgBIxOwkxiNJDj/75HuiLVTlewHUqlL+/Xni0AGIJYc0OOfmkke9kKOIm/rxeccoN9YBzJywDboJoHNlTbs86s/JyBHhodeDdMKyt5AI0UjuojJLbIyiqS+m1DCtabpMojLHq+onq9eYjP0ScXj4q+9T1kvYHBxUJvx2WoslrbQmj12XMZe/4pZiCY0ax4Os4H3vbxWJg7puPvwmbpj2FgA9M1BwlhSD/+Ka/4njhyrrEgEbhb2zATWj0dJWSkizADB5bDpjathAYoEe0ldNHPi8tDGV/HFwCcclXYqM+ioSh9pTK4X/VR8un8wuXGI+T9EHo+OQAHdvJdWrvmjicEVpHCEtIyU5CLPPplYSL91t8IQQP3gXyYEv1J/UfNRAMq7BaMTNqUVtXXbd5DgsA5UgZ1jl/uW63xyu3kqxVxiLfmvTpLuOnathiiWtiorjtksmGM9z8HcsqpV0hDaXo8RCrpsfSZuDSVJM2wDPHTdEWffN9btSGQVU/e9XLOCHDtDv4tgi3JApDqXI7ClG+PUxcTC/7yJjWCuoOoF4nfDEoQMQi/3pF8Mn7eqte9GroYDXrjlb204hR2gOPxjTpFDpPQ/qpVhYFIveF953mLLNvKMrq6sXkHze5B8vSgKm+/7LFwK9P+e4/0Ox3wKH2IqN21YRB5njzupPr4K8YHDukhNH3a0n9gsxELoBPQPi0ytUqbUWS/o06dLYdQuDKI3Y3recGVcHPn9djPt5idPWqp8SaiV93/xcMsGl+6Joy0o8sGcDbvjUcYb+k3ZFk9qvWNK7GHPE3krBPeiJA0UBeKdcOyNxjn9PLh6V7UXNEQcOk+Swp6WIkU09lMYijnwu1nHqOD5AzcXefflJqTLVq9arD9SurLIXhDhXTWMM2hSv09dtyMdcvimx3+GDAgPxvnAy6+Irgr7diJhuQxeTOq1SICL071EfvXMdsfnTF2JjuImR/GaY4punf97fVtIaeuX703nTZXFAcE0aGBukReKgXyTFjXw2Seob1djMQXBp3XuWNVHF5IjP6ANHDMQnT0jnvYrHljRy674hvo+FTo0mjyeSHAyrr26f+OZQW1GswkZANUccXGwOgF1l4mrslSd/fSGHw5rUaTFS49HMHiLC1j2xUY1Pyu4N+VQ9jiySg2nSipvdmCQM3h6fzLqAvmCgwqFFB+26ZWRHwEVXfsIh/fDxUPVhWvga6/Lo170umo/NraUopkCG3Jdu0XUl8MHYjKdT7Tzzln2fgqcWJb3rdC6i7t5pwX+X1PjK6xVTTrxvnZ0wqoukxKRTH3O1kk7yi8fDiYNZcuDnVDwAl2CrITnUnkEaes5P/PhtnJfIRZg8Q2Ruom+GkHd9pCWwbme8QHzptFG4cPxwHCIl8cqT2xjT7RvUSgIBtRmPgdh+YyIOrpKDvJ1otZEkDnZVom39FaN197cVtZKDPFcvGK+3RXG4OCC4oD2EdtKR6kBK1+cY2RwEb6UsWWZVc16c2yZplve/T7A/aoMPc4HkYNvzgvfdZrE5BH2nGaFnv3NGlEyyy9sciGglEc0noteIaG5Y1p+IphHRkvB/P6H+1US0lIgWE9E5Qvn4sJ2lRHQTVSK0UwOT5JDFTVM8b6orTtDxh/ZTqpTEcSXHo1crNUsBeDJh4PVcxijXNamLEvft4KPOPy7TBiWkOVa1KS8OqsCkLPtYZIH48ZsWguhZWoYhqsn2t5b0cSDCc77k5EOV7zroV7jGJa7FAeXuOvbZkw7Be4b1UZ5zlRwim4OgQsmyJKruUezbRXIQ3VxNzFpbiVk34clnkBxUHomHDOgezYUuTxxCnMEYG8cY4+4PVwGYzhgbDWB6+BtEdDSAKQCOATAZwM1ExN/OLQAuBzA6/JtcgXEpYUqfIU4m20chLvomKUOUML42abRWpaQKojPtBCdyNLoAvEyurMK969QbQZtu0ghfoJtbi6jP58y6ZaGZfQZPMdnmMKh3Axb+6JxUPVE6u/jEQ3D/FSen6pQDV6M9P2dbVjl3+NLKrXhn+z6tv734fGwLCkfFbA5lrhBmadru3gzE9zpb2NRKJzj87rPpHF+qdyR+j3bJIbA/cmiZtTD4zx5ompSEbAGfolqJp17n7/hAIQ4yzgNwZ3h8J4CPC+V3M8b2M8ZWAFgKYCIRDQbQmzE2iwVy1F3CNRVH7KuUfjNJTtuNCwDMxt6ENGL4EFSTXregL16/K9LlA4Y9DhylGyB5740a9UbqGge1UmuRoVejWXspvosXlm3R1pNtDpe9b6RycenRUIhcUZt6NmD8oeoAxawQ05mYnmekVrLaeYL3fuuzywEgYcwVIbZjarJ7ffycrZ4zjgJBuWolVyJmszHJqC9oFmhFQ6p3JH7XNskBoERGVK3NITQe25552iBtlrxFyYEHyx1IxIEBeJKIXiaiy8OyQYyxdQAQ/uebJgwFsFq4dk1YNjQ8lstTIKLLiWguEc3dtMltI+/UgKMMkelzWbw9XPWmdY46etWr1ompO/Yl90jQbeiSRa20TIjq1Kk3ZBgXSOGUKagOcDeOyq6svTSJ/MQ2y1GL2DhKW7v8lK1rnrmW1zNlzFUdqzA+zGpr85xx1dyWq1YyMUx5x/uRGbjvTh6D0484SFlX7baqkByEejZ3XqJkgjuzKyvD755ZZmwvi1qJb9vLwY3yEXE4AFxZT2WMnQDgQwCuJKIPGOqqngQzlKcLGbuVMTaBMTahqalJVcUK0yMV373tI3T9aMSPxCQ5fHjs4MTvldd+xMqJAMDfv3wKjj8kneYaSBIjG7ETxeexQ9WZW1PtO6gEADuHlsU4WmJBdkwA+NQEtWEWiCUx03s6YlBaxTe0b7fUnr4qGNVKkUHaLjmUWDxWbVbWDES+e5iuxO6dZjwt1CuPOLhIVrZxyF1fcdooLbftamcSx9XUs9FQM1iYRLdhXReyfUCnxpQT79nUSmL6nKJMHLq65MAYWxv+3wjg7wAmAtgQqooQ/uf+bWsAiF/zMABrw/JhivKOgcHmQGVKDiaI7m2mD+Y9w/pEEsBH3jNYWw9IbppjUtlkkRxEfHfykU71zC688bEtxsJ1/SEKfN4ZAz41YZixf/7pmN7T1ycdkSr77uQxOPwg++5pZpUAVyuZ2+CJ3biaoU6rMhGvMbfJ16jXhOyxur5FTHmvmtCKz3hgz3r866r0JjbK6yrgDZXekta9rg4Jj0Sbd1EuGQRnUt+Ki7VOjclpv5taidAq9M3tFAcEcSCiHkTUix8D+CCANwA8BODSsNqlAB4Mjx8CMIWIGohoJALD85xQ9bSLiE4KvZQuEa6pOEyurOK7mvv2ttR5ETabBIdrkjEA2BLGLuxsNm+t+WnhQx7QQy8aZ+E4Reg4WFP7MpIqOpv+O65rInacQ2sr2dNDcE7OdN+qJmz2ERe4LlTcID0zTCSoe+5Z4lVcIT6Wh75yKq49/1hlPfH59W6sw1BNfiEgqQZ1nW+ujgo2OBOjDHa4+kIu4SChq54jck4rDsSpSExxETlKBuBxYtC3Wz3+8u8n4iyH/Vbai/Z8CYMA/D18uQUAf2GMPU5ELwG4l4guA7AKwIUAwBhbQET3AlgIoA3AlYwx/uSvAHAHgG4AHgv/OgQmV1bxw9u8Wx1oxOHqxZElF9HLIUF6bslmYz2R+zFtFViu5OAKF+4ZcNB/C8d3/r+JxjZLLNDZ2twGSyU7cVAttK72FhP4k/yabQAAEo1JREFU0GxENohziKG7J1fXT8BdvSLeu0n9JT4P2/pblydw7aQzcTCec5+zroTE1VYIBO8jwaBr6udzyY2/dODPhDuTmL6LHBHe2hDbAbnNob6Qwymj1NsHVxplEwfG2HIAqcQkjLEtACZprpkKYKqifC6AseWOJQv4u1ZNjCxc2WurY88SE8/g6q2UBXySdavLmxdo0TW340JH1H0n1Er2YCOOIX30nCkRop3GbGo9F7WSchtVR08tE/g8snuIJT1SdGnfxWZ0G0BlhTgdTFPjIMFoa/s+gvdsD3oUYda7OzUBwH1+Z/Hgk+9BV72+kHNL8y9JDqbvIp8jLNm4K/rtIplUGjWXPqNk8FbKsn5u16QGkJHIQFmhp80ndb/u5mjrjpYcTBAX/Dqbz314emjfbon03TLyFHhwtJbsboOcOpjUfyrCaorxcAVv154JN6mr/ubZaRsIr8exs9m8d7D7fhJxmybvLNGZwDaHRE7YNWV4naGeeN8PfPkUYzuuEpP4TmwERX4uOjXQO9v2Kctl8DH+dmbg1WQMDKXk2uGyy12lUXPEwVWtZIPrYps1xbUL+ARvrDcvZGJ/Lu6ZHQabKiR87qMV3kNyvRILPEhsqqrY5qCvo1ocTDEervEfvF0XQ3wyd4+6fkIKs0lM4fx+z1B1dLLYN4eNy7/qQ0eG47CrYTh0u6aZrpEh3qo1b5Gz5ODOrInP5ayjDtKmC7epoDnknepM8yNHyUDXahigU2Ooeo+djOgRK95LlrVbfLHjFJuOx/XcctdnAedAbO2JpyvBEZeLWYbANsAeScyRywUcVInZjdz8PZskB9Vi3LNBL425PkPXGAuC28ZSWSKfOb502ihzm8LksBGHnOP9iHPdlZCa+hYlB5u6z9m1PJNaKX7f155/rNYl29WJQY6AN9mk5DarsWe0jJojDpy1skVIW9PvCotOH4N6J4sBzBW8TZtYnM/wcWXB9z92NK443bz4JMZh0707Di1HFKWOtnHlXHIw7rWheH4mtdavLzreZZjRvtD/WmomijkiJ121mC7bdt/cG89WT1zIbcZ9/pxsC34hoVYyE1IegGbqW+S0D7dkMi7L5mC5RtziNmt2AxV2SSpB03cxqHdyHnrJoQqI0mdYbA4jNLtDcbhycCJn5Mrd2JP+BW3avgdXvXJWfP7UkfhPx1gIwL74cEJt+8jExdSuVgr+m11Zk+dMG78AwAeOcAu8dN1v2pU4HDU4Dkq0SkxRQJ2NOIiMg1uQYpNqkyoB5wt7I9gCH7lLrEly2Be6Ph07rI/VpuDKdyU8vyxtim67pkSLI5vMawXH6WOS0d2mBb+3LDl4g3THw2RzEMXY31ycTuRVDrIYwFzr8bU2i1qpAxPdWmFbfFxiEoBAvcHVMNYU5E7EIfnbtPELx9fOPBwnHWbO1XTzZ4K5o9o8XgSRm1pJbMcqOYT3bXs+og7fLjkE/3sb0pUAwJcFabKH5d6jvg3Ege+Wd84xmj2eBbiqdt43OnYDtX1nfQWNgIkZ/PbZ+l0ORTT1asDEkfHc6WawGcoqzM6QHGpvP4fIW8k8Mfr30O8CBwCnjBqAGW9uxMeOUxupOBI2B0dSbOPQuORga68cgqDcwrSd6G4xnPN57+LD78oZ80U3y0ZMLvjWB+0LwUGhSsDmOSRuQ/neEeoUKDJsEitfRFyifzlsz9I1kaB43vbOOUyEaVi/7pjzvUlo0uzXLWKwwQVaxAmH9EO/7nXYtrfVwfvKzV7o6rYLxASpR33eGFDI1cDcTbYzvJVqjziE/23Lgm3iXPa+kfjYcUNSukEZ4odn4/R7Nxaws7kNt14y3liPLxCuNowejh8qALzgmB4hC2y6ar6g2RYfca11yTsFmLnD/7+9Mw2Wo6oC8HfekoSsj0ACZCNhEQghJhCWQMKSEsIDJAguLJIAKiCLuFUBAi6lFGApJYtVCJoUKApliQoIUkAJiCBLFJAYw2oJyFpgSAKELMcffTuv82am752X7tfd752vamp67ty5c8+cnj59zr333Lxc9fiC52td6BpP+Ow+2wa1ffacHVPfX+kGPYelDKzDxudlqJHs6YSNesTJI9O24gUYPSz9/xXT/QK9YGbj37Pr5ircOISGkX3E15WZnoVscRh4bMdmvPjWqoYbPOVJ/zMOKbmVkvhcThHxGgbY2L0PHXNIWwgGbEhHHeJqPvmtQ4Lv4iCnldTei74LK3nqvbb8gw3HoX/WtF8oJN7fE+ILo3cMpaXLOISmLBk/sv5GPzGxcRjqS5PehNcUstq8O9ttmT6AfMnRU7nuzy+ww+iwLXObYcl35qae87HoPnHaezBe+InpdRNK13z3dp5xijistPngdu6+uLPX1ylBfxxzcM++pfm+Lf9CacZzCK2327gRbD9qSNDd5ojN2oMvPJDP2ITvxI7TD/vCZMnVxKEypRnQeKVq1oT2rUVkQ/irGR2lcX7nznQMbmdMR9gddwjrNoT9ws8NX6hlr0kjuW7+jFwueoPaWz1J+sLaiW9W5uw8Ovh/cfzeE1Lfj6ezbu25sYw9BxGhzbNZVl70Q88hZUQ6QVbpJpJK9Q1QDh3Y5l0BC9GmLvd+7cBN7Vqv4RvoXR948UkaB1/YIiYt7u8b2+kp8WCvL5wmdHkvofL46NxtGzo9WX2bJf4Ns5qKnTf+9SVhnl1MyJhHzK5j0tPdx9NZOzzZDeIxhyIGomP6nXGI8U4DzcGnSpudAHDjF/bhriWvpa6bqCJHegbtQy8+axJjBN7ZSo60/9bUcR1sMWTAhmy4WRGLUW9/643rdQ1I+2YMFUnXXgIFdyQjYtsRetkNsYmjhw1kzbr1G+3GV4847JeWMBO6wkpFGoc+ou5wAh2HTBPVXXj4LkFz5CdtOcS7srUvEjqVNTlG4JthE+PbMSu5hiAr4pDKnpPSp7yKJJKwFRBTBv8dLHQZ2CKnQ2fJEe5mxTf9tZkZQg+dN4dHL/iYt158g+NbbR+Hlcxz6EXS9nNIkmUs9POzt+Pzs7fLrL2+Rvwf9F17kmMEobOV1nv+XBp8/xjOVsMHccsZ+zLZY3jiXFHgl+fCw3fh3qVvpNZplnu+eoB3yjZ0GeK8BvB7m3MP3ZkzDtzeu24jzpIbsv4l9Hxc505235jMADMOvU+o59BX7pKqwLrAsFJywZgv02vMeo/nsP+Oo7xpLnrC7g22bk2SlMDnCeVxgxE6UygOlbz/YT4D+L1Na4t4p9BC5FW+eMlhmV4L1iX2ZUgj/splr69IrZcn/S+s5J7t2r8xi07ek28eMbmQ7+4ac0ivt3pNl3Hw/bni6b6zdkifT37q/sV5dMndBleu9k9EKIohA6MQyHsBaUGGD2oLXtBXBbK+SdxgHDyexlsrsh0H6wn9znPonLI1H9lqaKFZSsvIQTuN5qCwLADB/O7M/YLqhc5WSiag860xWXzRwUHfLSKMH7kZL70dlpM/L8Zvnr5+oUjm7ro1R08fy9fn+k+Qp749txd6VJ+T9p2YaYLJPFgb6Dn40vH3Bv3OOGy7xRC29STV891tGmFMG984lXmSeMGSL/6djHmPSUk90Cy3nzWbt1aF5eTPgwFtLd7FbUUyqL2Vyz8zrehuePn2kbsW3QUvg9pbWb12vTcR5nF7jue91Wu96ybypDTGQUQOBa4AWoGfquqlRfTjme8VsxqxP/PxqWNY/v4aPj0jPUXAeZ078+WbnwCynTAwYnB7odOHv9gPZ6j1V6aN7+D+Z970hrXbWls4reDzohQ+mIi0Aj8GOoHJwHEiUkgAfEBbixmHXqalRZg/c6I3zfNRntQEVeWsOTsU3QWjl7jq+Olcffx0xpU4jBhTFs9hL+A5VX0BQERuAuYB/yy0V0bpuPOc2d4FRFUjq9QZRvkZPqidI6amLwotC2UxDmOBlxKvXwb2LqgvRonJY9FaUdx+9iweefHtorthGHUpi3GoF8epmaAuIqcCpwJMmFDcQI1hZMGUsSOYMnZE0d0wjLqUxZ99GUiORo4D/tu9kqpeq6ozVHXGqFFhWzYahmEYzVMW4/AYsKOITBKRAcCxwK0F98kwDKPfUoqwkqquFZGzgLuIprIuVNUlBXfLMAyj31IK4wCgqncAdxTdD8MwDKM8YSXDMAyjRJhxMAzDMGow42AYhmHUYMbBMAzDqEHSNmAvMyKyAlhW560JwH8CmhgBLM+wXh5thsqSx3fn0abJk129PNosUp485M76WtBM3TLLs5OqDvO2oqqVfACPNyh/M/Dz12ZZL482Q2UxeUyeqsuTk9yZXgv6ijyNrp3dH30xrPS/wHq3ZVwvjzZDZcnju/No0+TJrl4ebRYpTx5yZ30taKZuVeRpSJXDSo+r6ozQ8irSl2QBk6fsmDzlJit5QtupsudwbZPlVaQvyQImT9kxecpNVvIEtVNZz8EwDMPIjyp7DoZhGEZOlN44iMhCEXlDRJ5OlH1URB4WkX+IyG0iMtyVDxCRRa78SRE5MPGZPVz5cyJypYhvF9d8yFCe+0RkmYg84R6jC5BlvIj8SUSWisgSETnHlY8UkbtF5Fn3vHniM+c7HSwTkbmJ8sL1k7E8ldOPiGzh6q8Ukau7tVU5/XjkqaJ+DhaRxU4Pi0VkTqKt7PUTOt2qqAewP7A78HSi7DHgAHd8CvBdd3wmsMgdjwYWAy3u9aPATKKNhe4EOisuz33AjIJ1sw2wuzseBjxDtAf494HzXPl5wGXueDLwJDAQmAQ8D7SWRT8Zy1NF/QwBZgGnA1d3a6uK+kmTp4r6mQ6MccdTgFfy1E/pPQdVfQDovpfiTsAD7vhu4Bh3PBm4133uDaKpXzNEZBtguKo+rNEveQNwVN59r0cW8vRCN4NQ1VdV9W/ueAWwlGjL13nA9a7a9XT91vOAm1R1taq+CDwH7FUW/WQlT+/2ujHNyqOqq1T1QeCDZDtV1U8jecpCD+T5u6rGm6AtAQaJyMC89FN649CAp4Ej3fGn6NpF7klgnoi0icgkYA/33lii3eZiXnZlZaFZeWIWOZf4oiLc/CQiMpHozuYRYCtVfRWiPwCR1wP19wofSwn1s4nyxFRNP42oqn58VFk/xwB/V9XV5KSfqhqHU4AzRWQxkTv2oStfSPTDPA78CHgIWEvgHtUF0qw8ACeo6m7AbPc4sVd7nEBEhgK/Ab6squ+mVa1TpinlhZCBPFBN/TRsok5ZFfSTRmX1IyK7ApcBp8VFdaptsn4qaRxU9V+qeoiq7gH8iijWi6quVdWvqOo0VZ0HdADPEl1gxyWaqLtHdVH0QB5U9RX3vAL4JQWFM0SknejEvlFVb3HFrztXNw5JvOHKG+0VXhr9ZCRPVfXTiKrqpyFV1Y+IjAN+C8xX1eddcS76qaRxiGcWiEgLcCFwjXs9WESGuOODgbWq+k/nmq0QkX2c+zgf+H0xva+lWXlcmGlLV94OHEEUmurtfgvwM2Cpql6eeOtWYIE7XkDXb30rcKyLk04CdgQeLYt+spKnwvqpS4X106idSupHRDqAPwDnq+pf4sq56WdTR7TzfhDdSb8KrCGykJ8DziEa2X8GuJSuxXwTiTK1LgXuAbZNtDOD6AR4Hrg6/kwV5SGahbEYeIpoYOoK3CyZXpZlFpH7+hTwhHscBmxBNJD+rHsemfjMBU4Hy0jMqCiDfrKSp+L6+TfRhImV7vycXHH91MhTVf0Q3TiuStR9Ahidl35shbRhGIZRQyXDSoZhGEa+mHEwDMMwajDjYBiGYdRgxsEwDMOowYyDYRiGUYMZB8PIARE5XUTmN1F/oiQy9RpG0bQV3QHD6GuISJuqXlN0PwxjUzDjYBh1cInQ/kiUCG060QLF+cAuwOXAUOAt4CRVfVVE7iPKfbUfcKuIDANWquoPRGQa0ar3wUSLlE5R1XdEZA+i/FnvAQ/2nnSG4cfCSobRmJ2Aa1V1KvAu0f4aVwGf1CgP1kLg4kT9DlU9QFV/2K2dG4BzXTv/AL7lyhcBX1LVmXkKYRg9wTwHw2jMS9qVw+YXwDeINlm522V4biVKhRJzc/cGRGQEkdG43xVdD/y6TvnPgc7sRTCMnmHGwTAa0z23zApgScqd/qom2pY67RtGabCwkmE0ZoKIxIbgOOCvwKi4TETaXW79hqjqcuAdEZntik4E7lfV/wHLRWSWKz8h++4bRs8xz8EwGrMUWCAiPyHKkHkVcBdwpQsLtRFtwrTE084C4BoRGQy8AJzsyk8GForIe65dwygNlpXVMOrgZivdrqpTCu6KYRSChZUMwzCMGsxzMAzDMGowz8EwDMOowYyDYRiGUYMZB8MwDKMGMw6GYRhGDWYcDMMwjBrMOBiGYRg1/B8LPzmx1+2eBgAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXeYZFd95/09N1UOnaenuyfPKGehDAgJCRljgpHWMmsjr7GFeWWDDe8+wHq9Zl9edg2vQfviJZgFjIQxoBVgYUsIsEAIBSRGI2mSZjQ9uXs6p8p109k/7j23buXq7kpddT7P00/X3Kpbfaum6vzO9xcJpRQcDofD4bgRWn0BHA6Hw2k/uHHgcDgcThHcOHA4HA6nCG4cOBwOh1MENw4cDofDKYIbBw6Hw+EUwY0Dh8PhcIrgxoHD4XA4RXDjwOFwOJwipFZfwFrp7++n27Zta/VlcDgczobixRdfnKeUDlR73IY1Dtu2bcPevXtbfRkcDoezoSCEnK7lcdytxOFwOJwiuHHgcDgcThHcOHA4HA6nCG4cOBwOh1MENw4cDofDKYIbBw6Hw+EUwY0Dh8PhcIrgxgHAs8fnMT6baPVlcDgcTtvAjQOAj33vAD7/xLFWXwaHw+G0Ddw4AEipBubi2VZfBofD4bQN3DgAyOoGFpLcOHA4HA6DGwcAqm5iIaG2+jI4HA6nbeh640AphWqYWEypMEza6svhcDictqCqcSCEjBFCfk4IeZUQcogQ8iH7+CcIIZOEkJftn7e6zvk4IWScEHKUEPIW1/GrCCEH7Ps+Twgh9nEPIeS79vHnCSHb6v9SS6MZFJQClAJLKa4eOBwOB6hNOegAPkIpvQDAdQDuI4RcaN93P6X0cvvnMQCw77sbwEUA7gDwRUKIaD/+SwDuBbDb/rnDPv4+AEuU0l0A7gfw6fW/tNpQDdO5zV1LHA6HY1HVOFBKpyil++zbcQCvAhipcMo7AHyHUpqllJ4EMA7gGkLIMIAwpfQ5SikF8CCAd7rOecC+/TCAW5mqaDRZzXBuLyR4UJrD4XCAVcYcbHfPFQCetw/9KSFkPyHk64SQHvvYCICzrtMm7GMj9u3C43nnUEp1ACsA+lZzbWvFrRzmk1w5cDgcDrAK40AICQL4HoA/p5TGYLmIdgK4HMAUgM+yh5Y4nVY4Xumcwmu4lxCylxCyd25urtZLr0hWc7uVuHLgcDgcoEbjQAiRYRmGb1FKvw8AlNIZSqlBKTUB/C8A19gPnwAw5jp9FMA5+/hoieN55xBCJAARAIuF10Ep/Qql9GpK6dUDA1VHoNYEjzlwOBxOMbVkKxEAXwPwKqX0c67jw66HvQvAQfv2DwHcbWcgbYcVeH6BUjoFIE4Iuc5+zvcCeMR1zj327TsB/MyOSzQcVXcZB14Ix+FwOAAAqYbH3Ajg9wEcIIS8bB/7TwB+lxByOSz3zykA7wcASukhQshDAA7DynS6j1LKor4fAPANAD4AP7J/AMv4fJMQMg5LMdy9vpdVO1k9F5Ce58qBw+FwANRgHCilT6N0TOCxCud8CsCnShzfC+DiEsczAO6qdi2NIGsrB0UUeMyBw+FwbLq+QpoZh00RLxZ4thKHw+EA4MbBiTkMR7w8IM3hcDg23DjYxmEk6kMiqyPjKorjcDicbqXrjQNzKw1HvQDA5zpwOBwOuHFwlMNQ2DIOsYzWysvhcDictqDrjQNLZe0NKACAZJa7lTgcDqfrjQNTDr1+yzgkslw5cDgcDjcOzDgEmXHgyoHD4XC63jhkdRMCASI+GQCQzOotviIOh8NpPV1vHFTDhEcSEfBYxeLcOHA4HA43DshqBhRJQECxjEM8w40Dh8PhdL1xUA0TiiRAFAj8isiVA4fD4YAbB2Q1Ex7JehsCHglJlRsHDofD4cbBVg4AEPRIPFuJw+FwwI0DVN0KSANAwCMiwSukORwOhxuHrJ5TDgFF4hXSHA6HA24coOoGPKL1NoS8EhI8IM3hcDjcOGR1Ex6ZB6Q5HA7HTdcbB1U3oYgu48CVA4fD4XSfcUhmdUwspUApBWAHpOVctlKnF8Fphom//8VxPtSIw+FUpOuMw4PPncZNn/45MprVcC/rVg6KhKxuQjfMVl5iQ3npzDL++4+O4Olj862+FA6H08Z0nXHw2Sohbe+cVVe2UtDL+it17q56JW2l6i6necoup3VQSvH3vziO2Xim1ZfCKUP3GQfFqmlgbpWsbjh1DkGP9TvRwUHpGDMOKbXFV8LpZqZjGfz3Hx3B4wenW30pnDJ0nXHwypYBKKUcWGfWRAfHHdgY1JVVKIf9E8tYSXGlwakfadX6/nWySt/odK9xsD+cqlHCOHRwxlIsbb225RoXe80wceeXn8PXnj7RyMvidBlsc5bqYJW+0ek64+CTc24l06TQDOo03gt1wUwHphxqjTnMxDJQdROTy9w3zKkfLCEkpXLl0K50n3FwYg4mVDsrqVA5dLRxWGXMYWrFMgo8cMipJxmuHNqerjMOXikXc8jau5dcQNoe+NPJxqGGmINumHjL/U/hJ4emcW45DQCYi2ebcn2c7oAZBx5zaF+6zjj4lFwqa9awPpjdpRyqxxwSWR1HZ+J4enwe07Zy4MaBU094zKH9kVp9Ac2GBaQzqgFVt5WD0z7Duq+jjUOmulspa78vJ+aSEAgBACwkVWiGCVnsuv0EpwHwmEP7U/WbTggZI4T8nBDyKiHkECHkQ/bxXkLITwkhx+zfPa5zPk4IGSeEHCWEvMV1/CpCyAH7vs8TYq08hBAPIeS79vHnCSHb6v9SLZyAtG44iyBrn+GRRCii0NEDf5hxiGX0spXgzN12Yi7huJUAYCHBayM49YEphyQ3Dm1LLdtAHcBHKKUXALgOwH2EkAsBfAzAE5TS3QCesP8N+767AVwE4A4AXySEiPZzfQnAvQB22z932MffB2CJUroLwP0APl2H11YSFpBOu5SD4toNBzwiEtnOzemPpXUIxL5dpp4jq1tf2HMrGZyYT0K0T+BBaU69yDK3Uger9I1OVeNAKZ2ilO6zb8cBvApgBMA7ADxgP+wBAO+0b78DwHcopVlK6UkA4wCuIYQMAwhTSp+jVte7BwvOYc/1MIBbmaqoN3kB6QLlALDOrJ25mzFNinhGw+aoD0B51xJ7XwBgfDaBPUMhADzuwKkfrM6Iu5Xal1U5kG13zxUAngcwRCmdAiwDAmDQftgIgLOu0ybsYyP27cLjeedQSnUAKwD6Svz9ewkhewkhe+fm5lZz6Q6CQKBIAtKaWzmIzv3WHOnO3M0kVR0mBbb0+gGUr3VgyoFx+VgEADDLjQOnTmR0HpBud2o2DoSQIIDvAfhzSmms0kNLHKMVjlc6J/8ApV+hlF5NKb16YGCg2iWXxSeLyGqmswiybCXAUg6d+oFlbiRmHMq1xGDBQsalo1EAXDlw6kdatT5jPObQvtRkHAghMizD8C1K6fftwzO2qwj271n7+ASAMdfpowDO2cdHSxzPO4cQIgGIAFhc7YupFZ8s5sUcPC7j4FfEjnUrsQK4MUc5lHMr5b/+Lb1+9PhlHnPg1A2mHNQOb5G/kaklW4kA+BqAVymln3Pd9UMA99i37wHwiOv43XYG0nZYgecXbNdTnBBynf2c7y04hz3XnQB+Rtk0ngbglQvcSm7loHSucogXKIdytQ4sW6k3oAAAhiNeDIa8mI1x5cCpDxmXYkjxwVNtSS3K4UYAvw/gFkLIy/bPWwH8DYDbCCHHANxm/xuU0kMAHgJwGMDjAO6jlLL//Q8A+CqsIPVxAD+yj38NQB8hZBzAh2FnPjUKryzmB6TdysHTRcqhnHGw35cLh8MAgOGIDwMhD+YS3Dhw6kPGpU5THfp92+hULYKjlD6N0jEBALi1zDmfAvCpEsf3Ari4xPEMgLuqXUu98CkiMl2oHFiNQ49fRtgrlW2hwdxKf3jTNrxxzwB8iojBkAfPn0w27Vo5nU3arRxq/L596/nT2NobwE27+xt1WRwXXVchDVgxh4xmOB9KVhgH2MqhQ4NkTDmEvTKifqVqKutlo1Hccv4QADjKgVKKBmUZc7oId9JDrems9//0GK7b0cuNQ5Poyl4IzK20ZLtVIj7ZuS+gSFB1E1oHBslYtlLIKyHql8unsrKGhC6j2R/0QNXNjk3z5TSXtGZAsosra2lXY5oUSyk1rwaH01i60jhYysHESlpDyCtBEvOzlYDOLM6JpTUEFBGSKCDikx3jWAhzK7ljMSF7vna8g6fkVePsYqrVl9AxZDTDSXio5bsWy2gwTOp0cwWsz+l3f30Ghtmw3JWupiuNg9dOZV1OqYj65bz7WGfWTow7xDIawrZK6vErWKngVhIInJ0dAAS9nT8lrxLHZuJ4/Wd+jr2nGpZh3VVkNAN9QQ+A2ozDQtL6rGZd7qiHfn0WH/3eAfya/580hK40Dj5FQEYzsJzW0ONX8u5jyqETM5Ziad1RAL0BxfnCFZLRDHgkMS+2EPJaRiWe6dy+U5WYtBsQnl7g6qEeZDQTfbZySNawEVtkxsGV5fT9lyYB8OLMRtGdxsGOOSyntLx4A2DFHIDOUw4zsQwOnltB1DaGPX4F8YxeMraS1c28flOAaxBSl7qV2OteqnGCHqcyac1AD3Mr1aBGWUfgjKtj8EtnlgEA8zzFuiF0pXHw2tlKllupQDl4Ok85LCZV/PYXn8ViUsWHbt0NAOgNWq+71GKX1cy8eAMAhLs85sBed6HaMk2Kbz53ii9QqySjGY5yqKUIrlA5/PNLkyAEEAhvJd8outY4mNRqJBftAuXw8tklTC6n8fm7r8CNu6w0wF7bKC4li91EWd1whiIxuj3mwNxpiwUL0aMHpvBXjxzCv7xyrtRpnBKYJkVWNxHxyRBIbUVwi0nL+DLl8IvX5vC6rb3oC3q4YW4QXWkcWF1DSjVKBKRt5dBB2UqqbmVzDEe9zrGegPW6F5LFX6ysXqwcuj3mwJTDoktpGSbF5584BoD7vVcDS0f1KSICilRTzGGhQDnEMzoGwx70BRTMc+XQELrTOCi5XXFhzMHPlEMH7ZBZXME91IilEZZWDiY8Ur5y8MsiCAESXetWspWDy6302IEpHJtNAOB+79XApsB5JQF+j5hXLV0O9r4z5ZDWDPhkEQMhrhwaRVdWSHtdwdbCmANLZe0k5eAYB6nYOCyWijnoRpFyEASCoEcqOz2u02Gv220c/u3VGQyFPejx893ramC1CjnlsArjoBuglFrGQRHRF1BwaoG3dWkEXWkc3O0yevyFysF2OXWgcpBdyoGl8Bb60AE7IC0Xi8pQBw9CqkYp5bCc0rAp7EXEr/Dd6ypwlIMswqeIq8pWohTQDIqUaikHWRQwH+eGuRF0pVvJHWwtjDnIogBFEjpKObAGg27jIIsCQl4pL1vp2fF5ZOxutYVuJcCKO3SrW4kph5W05hjb5bSGiF9Bf1DBPI851EzGZRyqxRy+v28CBydX8owym8XiU0T0Bz1Iu/qkcepHVxoHt3KI+JSi+wOK2FEfNtWwAtLumAMA9AUU50s3E8vgPV99Ho+8PFnSrQRYGUvxbHcHpIFc+m8sbdXJDAQ9mE+oaOAIko7CbRwqxRxMk+Jj3z+Az/7kKBaTKgK2qmdDqnyyiD47JZurh/rTlcahknIArKB0J9U5lIo5AEBPQHEWOpZtM59QkdHMolRWwOqv1K3KIZ7RnFgVC+KvpDVEfBIGQh6ohtm18ZjVMBvLOCNCfbJoTV4sYxymYhmouolnxhegGiY2R30A4PQE8ykiBuwWHPMlsu4466MrjUOlbCXASmftJOWgOW6l/FbbvX7F8eW6d8NllYNH6uoiuK29AQBW+i+lFCtpDVGfgn62QPG4A352ZAYvnCzd62h6JYPr/+ZnePTAFAArMcSvSGVjDqfmrUCzam9uhh3jUEo58Pe+3nSncbB3xUGPlOeHZ/hrzKDYKGiGCUIAUcg3Dm7lwNxLK2mtZJ0DYMUc4l0YkKaUIp7RsKXPmqC3mFSRyOowTIqIT84ZB75A4W9//Bq+8PPxkvedW0nDMCmePT4PwPoeVtpwFGYhjdh1OmwOCYs5AMWV65z105XGgblMSrmUAFs5dNAimDVMyKJQNKSHxRwopc7I0FhGs7OVSruVurEILqkaMCmwzTYOS0nVeb8ifhn9IXv32mHprGuJoWR0o+xnhA2bYs0LvbKIsV4/4lkdCyVU16n5JDySgGu39wKwxtUCObeeTxadlOxS53PWR1caB+ZWKmcc/EpnpWxqOi0KRgOWcsjqJlKqUaAcSruVQh4JGa0zByFVgi12W+zZ2wtJ1RmxmqccOmiB+uovT+C2+59a9awEVTfLKoHCmIxXFrFzwHLVHZ8rrlU4tZDC1j4/3nT+IIDc+8+GVPlkEV5ZRMgrdZxhbge60jh47YUvWiJTCWDZSp3lVioMRgO5/kqLSdWR6gsJFSZF2WwloPuqpNliF/UriPhkLCZVZxcc8cno8SsQSC6o//CLE3jXF5/Z0NlLx+cSGJ9N4KUzS6s6r9K0wFjB5EGvLGDnQND5e4znji9gIZHF6YUktvYF8Ac3bMM3/sPrsK3fMiRutxJgTSnsJMPcLnSlcZBEAbJIECmnHDxSZwWkDbMoGA24WmikVCzabhK2wJWrcwCs5ntmF03fYsoh5JXQa7vi2O416pchCgS9gdwC9eBzp/DSmeWyIy2/+avT+OWxuaZc+1phbSp+fGh6VedlKyqHQuMgYiTqg1cWcNxuQ2KaFPf8wwv4L48cwumFFLb1+eGVRdx83qCzYXFnKwGWe5R3Zq0/XWkcACsYzVoGFxJQxI5KZVV1s2TgnfXTX0yqWLLdSiywV6pCms10+OEr53D5//MTPDM+36hLbitys7dlxzi43UoArEK4RBZnF1PYP7Fin1fa9/7lJ4/ju78+24QrXzus9uDHh2ZWpYCYcijljoqlc0ZDEghkUYAgEOzoDzrKIZ7VoeomHj0whaxuOmoByMUKl13ZSoDV8qaTNnPtQtcahy+850rc+4YdJe/zKxLSmtExs2lVwywZcxgMWb7y2Vi2aK6Dt6RysIzDwy9OIJbR8ccP7sUrZ5cbcMXtBdsJR3wS+oMKZuPZXEDaNg4DIQ/mEioes9M03ecVktUNZ2fermTs7qdnFlM4Mh2v+TyWdlqq6nklraEvoMCviHmFqDsHg07MYaVgrvm2vpxxyCmHfLeSVxba/v3ciHStcbhhVz9Ge/wl72Ntu9M1DCHZCJSLOQxHvJAEgtOLSUc5MEr2VrKNw8n5JK7cEoUsCnjg2VMNueZ2gvnKQ14ZW/sCOLOYwnJKhSIKziI3EPTg9EIS3917FixjuNDHzsjqplMl3K6kVQM7BgIQCPC1p0/WdI5umM6GqpRhjGU0RPwy9gyF8rLhdg4EcHYpZY/utT6HIVulllIO7mwl9rtTvqvtRNcah0p0WttuzaAl3UqSKGCkx4fTCyksFYxMLVcEx7j1giFsjvq6ou4h7riVJGzvD0DVTRyeiiHsk5304Jt290PVTZyYS+KOizflnQfkssAAy/XS7sYhoxkY6/HjAzfvxMMvTtQUe1BdWWyl0lljaQ1hr4ybdvVjz1DQOb5zIAhKrU0HU2Qf/Y3z8f437sDmSG4GCatQZ24lZizYZEdOfenKrqzV6LSBP+UC0oCVHnh0Oo60ZmDPUBCv2P7ySgFpALhySw9+dmS2K3y98YwGUSDwySK22zvZV84uYzCcW7h++8pRvPPyEcwnslhIqnjswHSecXjXF5/Bb126GX/+5t1QDbPtd7oZzYRPFvGhW/fgF6/N4T//80HcfuFQUa2MG9UVgC+V0RbL6Ij4ZHzk9j15z+POWGLhjWu39+L3rtuadz77TCZVA7JInA2PlyuHhsCVQwmYckh2yK44WyYgDQBb+/wYt4OBW0v4d90wt5IoEFw+FoW/w1J+C8loBr71/GnMJ7IIeSUQQrDDNg5soXMjCASDYa9TP+PePZ9bTmMmloFuUlCKtt/ppjUDXtnqUPwbFw9jLp4tm33FcBuHUm6luN2osNDAMIN7aj7pZIGVyiQUBeJsctwxC68sIstjDnWHK4cSsLnSix1Skq8ZZp5LyM3W3oCzW2MVwEDpmINHslKALxgOO4NaZmOdm1/+7PF5/OUPDkIUCDbbrRsGQh4rm001SvblAnIKi2UrmSZFRjOR1U1nAW33AGrGHqYD5NyJyaxesiEjw208SmVqxTIawt7iz6FPERH2SnmjVsu9tx5JhGboef3RvLIA1Y53FLaI4awdrhxKMNJjlelPLqdbfCX1QSuTrQRYyiF3260cihcBQgi29QVw83lWxarVUXPjq6uVtIYfvnKu6DirujVMipDHWqwIIdhuV/VGyyxgAUWEQHK7Z7ZoZnXDub0RlAP7DDjTEaukd7uNQ2EhHKUUsbSOcJn3bCDkcbLAfLJY8vMH5OIObuXAbrf7e7rRqGocCCFfJ4TMEkIOuo59ghAySQh52f55q+u+jxNCxgkhRwkhb3Edv4oQcsC+7/PE1paEEA8h5Lv28ecJIdvq+xJXz6awF6JAMLnUIcZBLx2QBvINwrZ+l3Io4VYCgEc/+Hp86NbdAFC2F39GM/JcDO3O//fjI/jgt1/CxFIq7zjL4Dp/UyjPiLL0ynILHSEkr6Ec84dntZxyYMd0w4Tehu1Isprp7M7ZHIVqLWUquZUymgnVMBH2ln7PBkNezMWzVqfbMsWpQG7T4i1wK1l/gxuHelKLcvgGgDtKHL+fUnq5/fMYABBCLgRwN4CL7HO+SAhh/4tfAnAvgN32D3vO9wFYopTuAnA/gE+v8bXUDUkUsCns7SjlIJdZ7Fm/GgAY6/WDuYPLuQ8USXCku7/MFK97v/ki/uK7L6/zqpvDbCyDh/ZOALAGHrlZtNNVf/inN+HvfvcK5ziLO5RzfQCW4WCprMwQqEa+caCU4o8f3Iu//MHBss/TCgyTQjVMp9YlN1e9inEwygekmZsp7Cvt3nSUQ1qr+L4yd6dfKVYOPChdX6oaB0rpUwBKN2gv5h0AvkMpzVJKTwIYB3ANIWQYQJhS+hy1yi0fBPBO1zkP2LcfBnArKYxYtYCRHl/RTnKjYgWkS7+lPkV0iuF6/IqzsyunHNz4FREZzcwrFpxYSuGp1+YwtbIxDOtXnz7pLNiF8ZOlpIqegAxFEiC5lJfjVqqwww15ZaeymqmrrGZCNazblFqL6amFFPatsn9Ro2E7cJ9iveaAp7YEjXzlkB9zYIayvHLwYC6exXJKrfi+MoPljjkwg9HucZyNxnpiDn9KCNlvu5167GMjANx9ASbsYyP27cLjeedQSnUAKwD61nFddWE06usct1KFmANgxR1CXmu2BdvZlfP5ugnYWV3uHdsjL5+z/2b7V5enVB3/9PwZ3LjL+rjNFsxjWEpp6PEXt1jZ3m+lXlY2Drn25myxdcccACCjWq0mTi+k2qoaP+0a4wm4A9LVYg65+wvdSjnlUD7mkNYMnFvO1KQceMyh8azVOHwJwE4AlwOYAvBZ+3ip7SmtcLzSOUUQQu4lhOwlhOydm2ts47KRHh+mY5mOaE9drkKaceWWHuwZCgHIuUpKZSsVwnZvrFiQUopHXp50/ma789iBaSSyOj54y26IAsFsPN+ttJRUneaEbi4dieCT77gIt124qexzh93KwTEOZr5x0A2ksjpUw8S5NnJhZgqMg1P3sxrlUPBY1lep3MI/GLbU6+Ryumy3ZMCtHHLuKR5zaAxrMg6U0hlKqUEpNQH8LwDX2HdNABhzPXQUwDn7+GiJ43nnEEIkABGUcWNRSr9CKb2aUnr1wMDAWi69ZkaiPpjUGm240Xh2fB7js7l+OOUqpBkfveN8PPT+6wHkZH8lpcFgiwardXhtJoHXZhKQBLIhAtIP7T2L7f0BXLO91+qZVOBWWkypTnNCN4JA8PvXbyubHgwAYZdycNxKrlRWwFpsWaHlyfnieQatosg42AtxrQHpUkOhHOVQIpUVAAaCuYLCim4lRznkPp9sk8JjDvVlTcbBjiEw3gWARdR+COBuOwNpO6zA8wuU0ikAcULIdXY84b0AHnGdc499+04AP6Nt0Ah/I6ez/seH9+PzT+RGNapG+SI4wFrsWJA54pOh2N0yq+GT8wOVJ+xiut1DobzgZDtycj6JF04u4q6rR0EIsbJlEsUxh94SbqVasBbIfOWgFhgHd7PDdjIOadW6RnfXU6C6cmCqqD/oKXYrpSu7lZhyAEoXwDGYuzOvCE5iyqG9P3MbjapFcISQbwO4GUA/IWQCwF8DuJkQcjks988pAO8HAErpIULIQwAOA9AB3EcpZeb8A7Ayn3wAfmT/AMDXAHyTEDIOSzHcXY8Xtl5YU76JDRh3iGU0Z7YApRSqbkIpE5AuJOKTawpGA8XKgVW3DoY8WEy2d3Hcz47MAgDeebkV+hoMeTDlUomGSbGc1koqh1oIeWXEMxoopXkxB7dxcE8vayfjwDqysl26IglQRKFqOxn22voCSpGhXXGaF5ZTDi7jUCHmwK7J685WsgPnXDnUl6rGgVL6uyUOf63C4z8F4FMlju8FcHGJ4xkAd1W7jmYzbDf82mhBaUopklndqe7W7UBnpZiDm3ddMZLXCbMSToNCZhzspmkDIQ/2T7T3Lo7tglmm1mDY4/SVAqzFjFKgt8IuthJhnwSTWn2A8rOVcu/LQpsaB3a97t15wCNWVw72a+sLKjhR8HpiGR1eWSib6BD1y5BFAs2gFWMO7Hy/LBUd4zGH+sLbZ5TBK4sYCHkwubyx0lnTmgGT5lp/sMBwJbeSm2t39OHaHbUli/kLAtLLKRWKJCDslds+W0nVTYgCcVJUB0JeLCSz0A0Tkig47996lANgpXS6A9Ju5bBg7659sohTC+1jHApjDoDlWqo1IN0X9DiqiWWls46s5SCEYCDowbmVTG0xB6U45pDlxqGu8PYZFRjt8ZUcfN7OsOKjpZQKSik03VqkazUOq4EFKpMu5dDjt+oC2j3mkNWNvKD7QMgDSnOT8Fg8oFS2Ui0w90k8o+cVwbldH+xvXbg5jLOLqbYJ4hemsgJWOmu1gDRLZe0PKNAMikPnYjg2YyVGrKSn2K6kAAAgAElEQVS1svEGxoDd5bayW6lEzIEXwTUEbhwqcPOeQew7s7ShgtIshVAzKBJZHVm76KpchfR68LOhSHZAejmtIupToIhWtlIb5BWUJaubeem67ql4QE55lapzqAWn+V5aQ8blq3fvvllc6OLNYZjUmrrWDrAOp+5Cs1r6aDHjxgzq737lV/jo9/YDAM6tZBxXbTlY3KFy+wymHFyprBIvgmsE3DhU4LevHAGlwA/2TVR/cJvgXnwWk6rj3qk1IL0amFuJKYellDXpi8U39DYq7CpE1c28wLtjHOxaB9ZXaa3KIVxCOQC5edRALuZw/nAYANqmIt9RDq73x3IrVQ9ISwJxFEI8qzuxhzMLSYz1lp68yGAZS5WL4IqVgyRa3YK5cqgv3DhUYKzXj2u39+J7+ybbehfsxt3TZjGpQrN3c7UGpFeDVxJBSC4gvWK7lZgLq50L4bJ6fmEgG9zD2kYvpuqkHFwxByC/rcSCndHFht0U9nZqFbn2GflupVpiDh5JcF57yCthOaXh7KI1aXBrFeNw0eYwNke8FetHcsoh//Pslfg0uHrDjUMV3n3VKE7OJ3FwMtbqS6mJRJFyWF1AejUI9nQ0JyDN3Er2F7hdfOilyOpGXuZMf9AyAqyFxlJShU8W8xbI1cCUQyyjO3UDQH5bCaYc2LCb6ZX2SP/NKYfaAtKmrRCZwX3dth6859ot+Mu3XgAAeHp8HkB+k8dSvOeaLXjmY7dUnDbnHg2ad1zpHOPw0Yf34yMPvdLqy+DZStU4z24rMZfIwCrebm8KjcMm28/bCOMAsM6suYB01KUc2jkoXehW8kgion4ZTx+bx83nDWAxqa3ZpQTkir3iGS1v0YqlNUgCAUVOnUT9MvoCCqbbRDmkNQOKlF8IWSkg/e4vP4vrdvRZ9TSSgKhfwX971yU4Om0Fo395zGp1U82tVEu/TfZ/5lfyly6vLHRMzOHV6fbYiHLjUAXJ9tWr+sZwK7l3d0sp1dm919IOYy0EPCLSqo60ajWVi/hl52+1t3Io7jf1zstH8I+/Oo23/89nEPJKeTMcVovHLhxbSee7lRJZHYokgMCK1SiSAFkUsCnixXSbdLLNamZevAFgAWkjLz0VAM4upvDSmWX0Bz0IeqS895QphaeP2cphHe8ng7n5Cgct+eTOUQ6xtJbXBbhVtP4K2hy20Olm+y50bli2kkCsVEkWkG6UcvDJ1qKxnM756NkC0c61DlnNLKoE/8TbL8KLf3Ub7rpqFPGMvuZ4A2DtgnsCMpaTWt5ApHjGMg6FIzg3hb2YbpORq2nVKHKnBTwSDJMWzZH+xWuWKoilNbsSP7/+YFPYi1hGR9QvV6xzqJU3nT+I733g+qJCTa8sdkxA2nJFtv61cOVQhY0QXHWTzOoQBYL+oIIlV8yhEQFpwFo0UqruVEdHfTJYklI7v2dZw0RUKV6sIj4Zn7nzUlwyGqnqBqlGj1/BQlJFWjOc6t94RoNHEiAJ+UNrhiJevHR2eV1/r15kdKPIp19ujjQzDvGMjpBXKqqA3trnx3QsUzXeUCuiQHDV1t6i494OUQ7WOFWtLRJguHKoAnMraRvErZTI6Ah6JPQGPFhMqo7fv9ywn/XiV0SkVMMpGovYbRCANncr2X71UhBC8N7rt+FN9qzstdIbULCUUpHRDETslhCVlMNiUs2bidAq0qqRlyoKlJ4jreomnrWDzfGsVtJVx0aq1ss4lMNSDu37eauVjGZCN6mTAdhKuHGoApPJ2gZyK1nGQc5LZW1cQFpEKmtgxVYObrfSRgpIN4KegKXe0prhFHbFMzoUUXDaQARcxgEonkbXCjK66dQTMIKe4jnSL55eQlI1sDniRSytOwFpN1vtueSNNg4+WeiI9hmstXlWN1s+AIobhyo4bqU23gW7SWZLK4eGuZUUCSlNdzqyRl0B6XZ+z7K6WdO0u/XQF1CwmFKRVg0ngKoaJhRJdHbmzK3EssqmWjg/ZCGRxan5JDKqkTcvASg9R/rXpxZBCHDbhUNWNX4Jg9tc5dABxiGdq4Np9evhMYcqOG6lNg6uuklkdQQ8Inr9cl6dQ6OylXy2cmBupahPgSxZWTftrBxKuUDqTY9fwXJKQ0AR81pCKJJQNIKTGYdWprN+9qev4RdH59AXVIrSeFnqqDsbbv/EMnb0B7A56rNanKdUp16EcflYFKM9Ply9rQeNpFOylWKuIsmUqlcsCGw0XDlUQd5gbqVE1kDQK6MnoCCW0R3fZSN6KwF2cZSqYyVlBVp9iphTDm1tHIyGu5XYAptUczEHwEpz9TrKwfryD9lupZkWKofFhIrJ5TQml9JFMYfCOdKUUrwysYJLR6NORfRCQi0yuJujPjz90VuwazDU0Gu3AtLt+3mrFXd7lVZnLHHjUIWcW2mDKIeMhpBHQp+9MDEfdiMD0hnNxEJSdXbHThFcG79nakHjvUbgbvft7hfkkQRn8WW+/LBXgk8WW6ocUvbOeyGpFmUrFc6Rno5lMBfP4tLRCMI+u49UVm+YQq2GRxZa7oapB263UquD0tw4VEEUCASyceocklnDcisFrCZmbAZ2o760zGc+tZIbDK9IdrZSmyoHSq18fU+DFzL3iNE8t5IrIO23d+SEEKsQrpXGweUyKpfKygLSr5y1BiO5lQPQuNhWNXyyCFU3nVYeGxW3cuDGYQMgi+0/n4CRyOoIemQM2F1Gz9lVt41snwEApxdSzgKoiNbC0q4BafZ/WZiRU296ArlF060clDzlkPMpD4Y8mGthtpJ7MfIWqKrCmMOByWVIAsFFm8N5oz8bHeQvBzNmmTZIBV4PeQFpbhzaH1kUNoRbyTQpkqqOoEd0WlCzGdiN2tEx5TCxlMYOu7uo3ObKgdVfNDrm0BfIzUUOeCTHtadIgjMD2e+qRA77ZKfCvRW43TKFMQdFEqBIAhJ2ttL+iRXsGQrBK4t5lc+tVA7Axp/pEM9TDq37LAA8W6kmZJFsCLdSSjNAKRD0Sk5vfDaoSBIaE3O4cVc/fvvKEdx+4RDefMEQALR9QDrbwDbmbtyuJJ8swiOJ0AzdCkjbO+yASzmEPBISWa3oeZqFezEqdCsBVlyELV6vTsVwy/mDznFGq4wDUzobPe7gzlZq9WvhxqEGZFFo24XODZP8AY8EvyI5nTQVUaip4+VaGAp78bl/d3neMZYZ1a4V0tkmKQevLCJgN6zzKdbOG1nkVUgHXN1Fg14pbx5Hs0mpBvqDHswnskXKAbBmVMQzOiilWEppjusyL+bQooC041ba6MYhrcFn12zwmMMGQBaFts68YbBdHfNjM9dSozKVyqG0Ycvul88uOzvjnFup8f5xlrHkkyXHGCmi6HQ9ZVlAQK4tdqv66qRVA1dtjQLIjYB1E/JKiNvDiwyTOu4kryw4yrTRGWDl6BjjkNGdmhduHDYAG8WtxDJJWICQ7eyaLfXbLf33pTNLeOcXnsHnfvIaADj9i5rxvrBaB58i5oxDid5KgKUcNKO482kzUHWrp88lIxHc/zuX4TcvGS56TNBjuZViafY5s4wDIbnRoK1SDr4OMQ7xjIYh2yXc6tfCjUMNbDi3ku2qYKMvG5WpVA5RIBAF0pL37MXTi3k9aSil+H8ffRUA8L19E8hoBrJac9xKQG7+AIs5APkV0u6hNaGCdNFmwlSVT5HwritGES3RrpwpBzbqlNU3sPuA5rynpWDvp3vq3kYkltbQF/BAFEjLA9LcONSAtNHcSt5Ct1Lz/5tlkTTdrXRyPol3f+k5/Ov+c86xxw9O48XTS3j7ZZuxlNLw40PTuVTWJriV+gIu42C7XDySgJvPG8RHbtuD8zblKofZ/1sr4g7MheGvMBaVxRxY0NQda2DGoVUB6ZEeHwDg8NRKS/5+vYhldIR9EvyyyN1KGwFlg7iVmHJgropWuZUAy73Q7IA0K/g7OJlbIJ46NoeIT8Zn/91l2NLrxz89f8ZRDs14X1jMwasIzq7aIwmI+GT82a27IeaN4rQW29Yoh1qMQ75byZ2lxOIPLTMOUR8uGYng0f1TLfn79SKW1hD2yvApIq9z2AhsFLfSbNwqoGILUqsC0oC1SDT7PWPN/47Ys4sBIJbW0RdUIIsC3nbpMJ4/uejI9Wa4QPrsRnR+JTdCs9wCyox6vAXKIe0Yh/IJjCGvjERWd6b+hX3FyqFVRXAA8NZLhvHKxArOLqZadg3rIatbo3bDPtmZk9JKuHGoAUkkbRNcrcSp+ST6AoqzixsMWTGHVuzm5BYoh8WktWgddRmHlbTmVCf3By1jOZ+wHteMzJrfuXoMX3jPlQh6clPSygVt2QLbyphDJeXAlMK5ZUuhuSujmYupVQFpAE4Q/UcHN6Z6YJuCsFeCT5G4cdgIyKKwIbqynlxI5s3WZYVwrYg5tEQ52MZhNp51bscymmMsmZGYjTe235SbvqAHv3mptWh5alYOzS+EY033CmdHuwk5xsEqrHRXRrfarQQAW/r8uGQkgp8cmmnZNawH1joj5LWUQ1rjAem2R9kgbqVT80lnuArQ6oC00PQZGIu2WwnIuZbcyoH9nrPdb43urVRILuZQ+u8GW6gc0jUGpAHLOLgzrqz7WhuQZly0OYxTC8mWXsNaYU33wj5pY7iVCCFfJ4TMEkIOuo71EkJ+Sgg5Zv/ucd33cULIOCHkKCHkLa7jVxFCDtj3fZ7YJbuEEA8h5Lv28ecJIdvq+xLXz0ZwKyWzOmbjWWzvz03civhkqydOK5SDKDQ9X38pqTq77yPTMQB2gM9OuYz4C4xDkxcydyprKVoZc3AC0nKlmIN13+RyOi8Y7b6vVamsjOGID/OJ9pjFvVqcOew+2aqSbnfjAOAbAO4oOPYxAE9QSncDeML+NwghFwK4G8BF9jlfJISw7cWXANwLYLf9w57zfQCWKKW7ANwP4NNrfTGNYiO4ldhuaXt/0DlGCMFA0NOSgLTcArfSYkrDzoEAegMKjk7HQSlFLKMXKQcWuG/2LpfFOMr9XY8kQBZJi+scyisHZrwml9J5LiUgF5xutXIYjrTPLO7Vwppkjvb4HeXw2kwcL55easn1VP2fpJQ+BWCx4PA7ADxg334AwDtdx79DKc1SSk8CGAdwDSFkGECYUvoctXoDPFhwDnuuhwHcShrVCGiNbAS30ql5K0NjW3/+rN53XzmCW+yGeM1EEUnTA9JLSRW9AQXnDYVwZDqOpJrf5oH9bpVyYAqunJIjhFgtNNq4zgGwJtuFfPnG4bLRKM7fFMJoT2NnRVeDtZ5gcZGNxJmFJDySgIGgxwlIf/JfD+Oj39vfkutZa+O9IUrpFABQSqcIIYP28REAv3I9bsI+ptm3C4+zc87az6UTQlYA9AGYX+O11Z1WuZVWUhoSqo6RqK/qY5lycMccAODDt5/XkGurhiIJTj1Bs1hMqtg9FIRXFvHjg9NOgK8o5pCwlUOT3W1OEVyFLKmgV2ppnUOphnuM/LqG/KXjvE0hPP7nb2jMxa2CzdHWz+JeK6cXUtjS64cgECsgreo4Mh1vWRuNen87Su34aYXjlc4pfnJC7iWE7CWE7J2bm1vjJa4eWRRaUgT3X//lEO75+gs1PfbkfBKDIU9eC+hW0orakKWUil6/gsGQBwtJFQuJ/Hx8ryw4xXmK1LhOteWolsoKWIVwralz0OGTRQgVWruHSmQntRubItZGaqqFs7gL+cpTx3HoXPXK7TOLlnEALAWXVA3MxbOIZ/SWGIi1GocZ21UE+/esfXwCwJjrcaMAztnHR0sczzuHECIBiKDYjQUAoJR+hVJ6NaX06oGBgTVe+uppRc4+AOw7s4Szi6maunSems9PY201cpMD0hm7xXFPQHHqO07MJwDkFIO7QVwrAqfuCulytGqmQ0o1KrqUgPzuq+6+Su1E0CMh5JEw1SZupYxm4L89dgQ/2DdZ8XGUUss49FnGoTD2w1yhzWSt35AfArjHvn0PgEdcx++2M5C2wwo8v2C7oOKEkOvseMJ7C85hz3UngJ/RVvUsLoMskqanZcYyGk4tpJDVTSSrZC1QSnFiPontfe1jHJpd57CcshbU3oDipPCOz1rGwb3LjfhaV8lbrc4BaJ1bKa0aFYPRgGVcWVZSqE2VAwAMR71toxxYYLyaGlxIqkiphqMcCt17rDanmVQ1/4SQbwO4GUA/IWQCwF8D+BsADxFC3gfgDIC7AIBSeogQ8hCAwwB0APdRStnK9gFYmU8+AD+yfwDgawC+SQgZh6UY7q7LK6sjrXArHZqMObfn49m81s6FHJtNYDGp4kq7F387oDS5zoFVR/f4Faen1LGZfOUAoKXKoSegQCCVF9agR8KJudbEHKopB8C69qWUVhRzaCc2RXxtE3Ng1xGrUth4esFKKHG7ldy0QjlU/R+mlP5umbtuLfP4TwH4VInjewFcXOJ4BrZxaVdYQReldF1+ajZhq5a4gNtHuZDMVnQZ/fKYFbu/aXfzXG3VkJucrcRyxHsDilMZPj5nKweXCyTSQuPwGxcPY+dA0JnxUAqmHL71/GksJlT82a27m3JtKc2Ar0JfJcZGUA6bI14cPher/sAmMGMbh3LKgVKKk/NJpx/UVsetZL3PmyNenFvJbCi3UlfB6gTWuxP+va8+j888fqSmxx5wdRZlvYDcZHUDv/V3T+Obz53CM+Pz2N4fqCmrqVk0263ElENvQEZ/0ANCrDgMkL+QRVqYj69IAi4eiVR8TMgeqPO1X57E9/ZNVHxsPUlldQRqUg7WotWuMQfASmedT2TbYkxtzjiUVg6PHZjGLZ/9Bb7x7CkAcFKB/bZb6dodfRBIrjanmXDjUAOs/cR6XEuUUpyYS9Ysdw9OruASeyGZT+Q+GI8dmMJsPINHXjqHA5Mr+JsfHcGvTizgxl19a762RiCLQlPnOTDl0OO3OrD2+hXoJkXII+W1xW6lcqiFoEdCVjdxYj6JpVTzAtOrcSsB7ZutBOQK4WbawLXE2sjHyiiHJ49auTwvn13GprDXNQTK+n3+phD6gh6uHNoVqQ5jL5dTGlTDrKlfSjKr48R8Em/cY7mJWErmXDyL/+tb+3DP13+Nr/zyBLb2+Z3nvGlX/5qvrREoUnMzvBaTKgjJLf4s7hAuKNbKGYfWtZauRNDly49ltLypdo0kXatbydP+bqXhNkpnna6gHCilePb4As7fFIIoECfeAFjqhxDg8rEoBlpkHNpXG7YRiu1WWs9OeMbONqjFOByeioFS4IotUYS9kqMcjs1azeRenbL8qff/zmU4NpPAPzxzCtfvaDPj0MQ6h3PLaczGs4j4ZMeQD4Q8ODIdL2scWt3moRzuxANKrd5QPRViFPUipeqOK6MSG8GtxArhJpdTAHpbei0sW6mUcjizmMLkchqffMdFCPvkvFjUjoEgnv/4rRgMezEY9rTErdS+/8NtRD3cSjP2h6QW48AmmV0yEkF/yOMoh+N2auZ9b9qJg5MxvO3SzRAJwR/etN1pKtcuyKIAkwK6YToLdiMYn43jzZ97CgCww92u3K51iBQsYswd0q5upVBBFtBynYzDJ//1MHoDCu5706684794bQ4LiSxSNaSyWtfX/m6lsV4/BAKcnGt9d1amHFTdREYz8jrZPjO+AAC4fmc/dg0Gi85lM+AHgh4cmYoX3d9ouHGogXq4lWZjTDlUT1M8MLmCgZAHg2Ev+gMep93D+GwCAUXE/337eXlZU2yITTvBduaaQdFID85Zu1nZrecP4ubzB53jjlupTIO4Zgz6WQtsVOhlY1G8cnYZSykV27H++pUnj87CpCgyDl9+8jhePrsM1TBrijkMR73wygKibbYZceORRGzp9eN4i40DpRTTsYzVYVUzEM/o+cbh+DyGwh7sHKj8/zsQ8mA+kYVp0ooV7PWmPb8hbYZcB7cSk4W1KIdDkzFcvDkMAOgPKVhgxmEugV2Dwaa3fVgLTG01OijN+id9/K0X4Pev2+ocZ4VwkXJupRZOLKvEtn4/AoqIO6+0Wo+t1CkonVINnJxPYiWd/3yTy2mkNatBYS3G4a6rxvDTv3hjxXGi7cDOgSCO26nMrWIlrUHVTUcVuOMOGc3AL1+bw027Bqp+nwdDHugmdZIumkV7fkPaDKUubiVbOVSpfk2rBo7Nxp1Mpf6g1ScIsJTDzhLysx1x4jQNDkoz41DoA2e1DhstID3a48fB//oWp2alXgtC0v7cHZjIpUibJsXUSq7NRC0LviIJGOttbefVWtg5GMSJ+WTTAvqlYC6l3fZ31h13+PGhacQyOn77ypGS57oZsF2kc4nmxh24caiB+riVbOWgGRV7Jb06HYNJgYts49AX8GA5pWExqWImli3pm2xHcm6lxhqHlYLOq4xczKHAOPjb260EWG0qeuzrXK6DcqCUOi1YXplYdo7PxrPQDOq8R7Uoh43CzoEAVN3E5FLreiyxNNZdQ8XK4X/vncBojw/X76iegs5cpM2eUdG+35A2oh5uJZatRCmQqdDK+pAdjGbFUv0hKxj5wkmrF+GugY1hHJhbqdHGIZbR4ZWFIiUw6MQcCgPS9jjLNnUrMUJeGYQAy3VQDlnddHbQ+13GwcrmAf7ghm0AULFye6Ox0/6etNK1NOMohxCAXJX02cUUnh6fx11XjdUUQ2DZVyfnmxtDae9vSJvguJXWE3NwWf1KQekDkyvoDSjYbBfy9AWsRe75k1Zmw0ZRDk7MocFupZWUVjJzZmufHx//jfPx1kuH844HPRKiftnZjbUrokAQ8clYTq9fOSRdrsz9LrcSmzz2tkuH8a9/dhNucQX0NzrtYRyy9rVYAWemHJ46Zo0beMflm2t6npGoDyNRH5493twRN+0dVWoTHLfSGttnUEoxG89gIGQVs6RUA+XE5MHJGC7aHHaCVAO2cvjxwWkoopBXKNPOMLdSowPSK2mtyHUEWK6Z979xZ8njj3/oDW2dbcOI+uS6VEmzJIjzN1kT8mZjGQyGvZi021qP9PjaPsC8WnoCCnoDSkuNw1w8i6hfRr+9EYmlLSM9E8uCEGC0p7Z2N4QQ3LirD48fnIZh0ryK/0bClUMN5HorrW2hW0pp0AzqtNROllEOGc2aGXuJq/8OUw4z8Sz+41vOa2jNQD1R1mlQayWWKW0cKrEp4s1LKWxXon6lLm4l1gL89guHQAhw/7+9BsCaBd3jlzvOMDB2DgRwfLZ16axz8SwGgh4EFQmE5JTDXDyDvoBnVd/lG3f1I5bRaxoaVC82xkrTYtbrP2e+RzbfuVw662szcegmzWvOtrXPj4/ecT7+959cjz9+w441/f1W0DS3UlorykjqFKJ+uS4BaebGvGpbL/7kjTvx7RfO4tH9U5hYSmNzGzVrrDdbegM4u5Rq2POnVB3//qu/cuaGFDKXyGIg5IEgWLPBWbbSbCzrxMRq5fqdlq+BFc41A24cakBe5y6YGYft/ZYfNJUtbRwO2jMcLt6cMw6EEHzg5p24ckvPmv52q2Bupaze2PGG5dxKnUCPX8Fyuh7Kwfo/CHpEfPi2PbhsNIJP/uthnFlMtVUn33oT8kp58ZZ6c3ohhWfGF/DcidIL9lw86xSohr2yM9NhNp51Uq1rZTDkxZ6hIJ4Zb17cgRuHGlivWylnHJhyKP2BPTC5grBXwljvxv/CjvX6oEgCvvDz8Yaqh1gHG4eIT8Zysg7KwV4g/YoEWRTw52/eg+lYBifnkxip0e+9EQl4rDnMjRosyTwAsyW6v1JKLbeSrRBCXsnJVpqNZ1atHACrCd/Rmea10eDGoQbW61Y6u5iGKBDssDMoyrmVDp1bwcUjkQ1RAV2N4YgPn373JfjViUV8usYZFqvFNCniWb1j3Uo9fgXxrL7udGAWc2BN/d64Z8DJoOlk5RDwSDBM2rBZ5hmNGYfi+oOkaiCtGXltXOJ2l935hOrU4ayG3oAHKymtYcauEG4camC9bqXTtnxnKZeljIOqmzgyFa86DGYj8a4rRnHL+YP4ud2zvt7EMzooLa5l6BRYRlVhy4vVwj5vrMhNEAj+6PVW/GojVDuvlYAdaG/UTO40Uw4l5juzFtsDwZxyiKV1LCZVGCZdtVsJsD4Ptbb9rwfcONTAet1KZxZT2NLrh99jfTkL3UrfeOYkbrv/F1ANs6OMAwAMhb1Oi4v1cHwugd/5++fyFspy1dGdQrROVdJscXSPp73zqlH87V2X4U3ndU5tQyHs9ZaL8a2XlK0cZkooB9ZmP8+tlNUcQ7IWt1LU/pzXo/alFrhxqAF5na0gziwksaXP7/TLL7T8D++bgG5QfPDW3XjzBZ31ZY36Zayk1y+FXzi5iOdPLuKVs7kKXxbg61Tj0OO3alzWm86aUnWIAslrUy6LAu68arRt51rUg6C9GWuUcsg4yqHYODjKwTV0Kp7RnccOrMGtFLU/D0vJ5jTg69xPRh2RhbW7lWIZDUspDVt7/ZBEAYokFNU5zMdV3LCzDx++bU/H5ZxHfDI0g65bCi+6mg8yVpyme51pHOqlHJJZawRoJ8SyVgNTDuXqitZL2lYOC8ks0qqB+/5pH556zap+LjQOLCA9s7IO5VAnN2OtcONQA+txK51ZsPKst/bZg8MV0fFVAlZQdT6RdaooOw22q1/vB9oZeDRXbBw6VTnU671LZvW8CXPdAttoNSqdlW14KAV+eWwOj+6fwh89uBdPHp3FXDwLUSCO+ov4ZBgmxZFpK9toLe1bckqSG4e2gZWrr6W30mnbOGzptbJDAoqEpMsHupLWoJu0LQf21IN6LXCLydzAI0asw40Dm7pWav7wakipRkd1XK0VZhCT64g5aIZZdlPIlANgTdQDgOGIF3/27ZcwsZRCX0Bx1o4bdlpjfL+/bwIRn7ymCn2mHJo114EbhxoghEARBahrcCudXrTK97fYysGniEhruZ0MC1z1BzunI6abaL2UQ5Iph1w7hE53K7GRofES84dXQ6JLlUPAjjmsRzn88YN78b2izw0AABk6SURBVJc/OFDyvozLODx5dA4BRcQn3n4R4hkdPz08k6cOLtocxo7+AGIZfU0uJaB+G61a4cahRmSRrMmtdHbR2kGwL2dAEfN2MmyAx0CHKge2cK9XCrOYw3wi60xHW0lrEAWCQIfuimVRgFcWEF+nWySl6h0Xy6oFlsq6npjD0ek4Xjlbup+RO+twcjmN8zaFcNOufkR8MpKqkecNIITg7XYX1rWksQKAVxbhk8W69NuqBW4cakQShTW7lZhqACw/qDvmMG/70js95rDedNbFpIo+e97AuB13YE33OjnQGrKLp9ZDImvkpbF2C05Aeo3GlVIrHnhmMVUy2y6tmhgKe8A+fhcMhyGLAm6/cAhAcVzh7ZfZxmENmUqMqL8+nXprgRuHGpHX6lZaSGFrr9s4iHk7mfk4cyt1qHGoQ4YFpRQLSRVXb7P6Sx234w4rab1j4w2MkFfKGy+5FlKq7rhYuglFEqCIgtNbarWspK1uymnNcLKP3GQ0A0GP5GxaLhi25r6/9RJrhkihcdgxEMT737DDMRJrwerUy41DW6GIZNXKYT6RteVm2Dnm9xQqByurIdqhi1zII0EUyLqMQ1I1oOomLhuLQhEFJ2NpJa11bHU0w1IO6zMOyazelcoBAPweseJwrUq4DcLpxRS+8PNx/N0Tx5xjzF3HahaYcbhxVz9u2tWPG3YWT235+FsvwJvWMVQp6pO5W6ndkERh1TGHvaes0Z7XbM91VPXLBcohkUVfQKlpXOBGhBCCsFdal3FYtF1vA0EPdgwEcHgqBkopjs3EMdrB7R8AqzXIet1KyazRsXGZagQUac1FcCweCFgegAeePYWHXjzrHEtrBnyyiMGQ5Vo6f5M1DlSRBPzjH12L1+8eWN/FlyDqr890wFrozu3EGrAC0qtzK71wcgkeScAlI1HnmN8j5pXzzyfUjnUpMdY77nLBTmPtCyq4dnsvHto7gRPzSUytZHDt9t56XWZbEvJKmFop7t1TK4ZpuUW6VTkEPWtv2+1WDs8dX8Bs3JrgltEMeGURac1ExCfjstEIEk1SZxvGrUQIOUUIOUAIeZkQstc+1ksI+Skh5Jj9u8f1+I8TQsYJIUcJIW9xHb/Kfp5xQsjnSRtGGOU1KIcXTi3gii3RvBYFfkVESsu1Ee7kAjhGxCevTznYmUq9AQ9u3NWPtGbgy08eBwBcu73cwNXOIORZX0CauVQCXZitBDC30tpiDixZJOyV8PjBKQBWwRtza2ZUA35ZxIdvPw8P/8n19bngKlgDoNSmdGath1vpTZTSyymlV9v//hiAJyiluwE8Yf8bhJALAdwN4CIAdwD4IiGEad0vAbgXwG775446XFddWa1xiGc0HD4XwzXb8ne2fiW/jfB8PNuxNQ6MiF9Zl3FgNQ59AQXX7eyDQIDvvzSJHr+M3YPBel1mW+KeA7AWWNp0NyuHNbuV4lkoooCLRyJIugwMK8RMaTp8truuWfvZHr8M3aR519MoGhFzeAeAB+zbDwB4p+v4dyilWUrpSQDjAK4hhAwDCFNKn6OWOXzQdU7bsFq30r4zyzAp8LrthcbB+jCl7SEk8wm1Y2scGBGfvK5UVtY6ozegIOyVcdlYFIZJ8bptvR0bq2EEvRJSqrGmNGogl+PfjdlKAOtIsDbjMJ+wNm5b7dnvl41FIZBcIWZaNZs+izzqa17zvfUaBwrgJ4SQFwkh99rHhiilUwBg/2ah+REAZ13nTtjHRuzbhcfbitUqh8PnrJGfl49F8467C3NiaR2qYXZBzEFadYbFc8cXMLmcBmC1zvBIgmNYb9pltSK4dkdnu5SAXAuNte5+WXyrG4vgAMuttNb2GXNxy+XL+qJdt6MXY71+J5U6ozW/LUkzm++t1zjcSCm9EsBvALiPEPKGCo8ttcWjFY4XPwEh9xJC9hJC9s7Nza3+atfBao3DYjILnyw6X24Ga4mwnNKcbIj+UIe7lXwyYhm9Zj8ppRT3PrgXX3pyHIDlVuoLKI50v+PiTejxy7hlHSmBG4X1ttDIzXLoTuUQ9EhrrpCei2cxEPRgm20crhjrwa6BIMZnE6CUIqXq8DVbOTSx+d66jAOl9Jz9exbADwBcA2DGdhXB/s3GgE0AGHOdPgrgnH18tMTxUn/vK5TSqymlVw8M1D9NrBKrdSstJjX0+ItrF3YPWT7yo9NxV1+lzlYOUZ8Cw6Q1734TWR3xrO5kiywmVfS64jIXbY7gpf9yO7b3Bxpyve0Eq+OIrTEo3e0B6cA6spUst5IHN583iP/8mxfglvMHsWswiJPzSaQ1AyaFE3NoFs1svrdm40AICRBCQuw2gNsBHATwQwD32A+7B8Aj9u0fAribEOIhhGyHFXh+wXY9xQkh19lZSu91ndM2RP0KFhLFVZLlWEqp6AkUK4JtfQF4JAGvTsVw1G7fu62vsxe51TYMm3UZBfa7N9DZBrQcuc6sa1vg2A6zU5sTViPokaAZFOoq50ibplWV3x9S4JVF/NHrd0CRBOwcDEI1TLw2Y7mWmq8cmjcNbj3biSEAP7ClvgTgnyiljxNCfg3gIULI+wCcAXAXAFBKDxFCHgJwGIAO4D5KKXMGfgDANwD4APzI/mkrtvcH8IOXJmtuYraUUp3+624kUcB5m0J4dTqGuUQWQ2EPRns6d8g7kFuYVtIaRnuqPBi5ge0sS2khoWLnQGdnJZVjvW6lqRUrbjMcWXs/n40MiwkkszoUqXb37VLKmvVcmCzCPoeHzlnN+JqtHPoCHvzLn96ELU0o/lyzcaCUngBwWYnjCwBuLXPOpwB8qsTxvQAuXuu1NIMdA9bu/uR8Ehdtrj7neSmpYrSn9H/gBZvC+OmrM/ArKVy1taejG8cBLuVQo5+UzdldSFj53HPx7Jo7WW501jvTYXI5g76A0vSsmnaBpfAmsnpJJV8KzTCdwsPCGqTNUcvInrQzlpqtHESB4JLR5syZ5+0zamRHv7VjOOGaJ1CJpZSG3hIxBwC4YDiExaSKiaU0rtxSw1Z6g5Pzk9ZoHGzlsJK2gvaqYWJoHZ0sNzLrVQ7nltPYHO1sZVqJ4BpGhd755edw55efBVDcSn8g6IFAgFMLtnHo4LYk3RmlWgMs+FmLcdANEytpzcksKIQ16AKAK7d2vnEY6/WDkPwRn5VgygEAXpu2zule5ZDb+a6Fc8tpR/V2Izm3Um3prClVx/6JZWyO+KCbJnYVFFlKooD+oAcn5lujHJoJVw414lNEjER9ODFffYFjgdfeMjL2fNs4KJKAizaHSz6mkwh6JGzrC+DwuRhMk+KLT447vvBSzLp62rw6ZdWLDIW7Uzl4JBGKJKwqWymjGfjViQVQSrlyWOVMhyPTcVAK/PVvXYjn/9Ob0Vcik3A44nVmw3eycuDGYRXsGAjUpBxYmlm0jFsp4pMx2uPDpSMReKTO/XC5uWA4hMNTMbw8sYzPPH4U/+Onx8o+djaWBSt8ZsZhraMVO4GwV8LkUhpv/twvcM/XX8C/HZ6p+PiHX5zA3V/5FV6ZWEFSNTDSxcZhtQN/2OfNre4LGQp7oZtWWjtXDhwAwI7+AE7MJaoWcy0mKysHAPjbuy7DJ95+UV2vr525cDiMM4sp/OSQtbD988uTZaumZ+IZ7LCzQg47xqE7lQNgBaV/cngG47MJ7J9Yxl889HLFx782Y6VI/8srVrlQNysHVt9Rq1vu1akYQl6pYgbhJlfmF1cOHADWJKekauS5PUrBlEOpVFbGdTv6cPFIc7IO2oELbffZPz1/Gv1BD7K6ie/++mzJx87Fss7O7fhcAmGv1NFfwmqEvBJU3cTWPj/+4IbtiGf0itX6LLbz2AGrk2g3G4fBsAchr4RH7feiGq9OxXHBpnDFDMI848CVAwfIpbOyrozlYE2xak2d6wbYYh/L6HjXFZtx7fZefKeEcUirBuJZHecNBUEIoBkUg10ab2CwoPTbLh12XJWVGhky1ydLx2Tpl92IVxbxwVt248mjc3jy6GzFx5omxZGpGC4YDlV83CbX57HZvZWaCTcOq+CSkQiCHgn//78dg2GWdy2xlM3eCsqh29gU9jrtRNgYxZPzSWS0/CwSlqk0FPY6ymuoSzOVGCGP9b791mWbq1bIJrI6plYyTiBWEQX0d2l1OeOeG7ZhW58fn378aMXHnV1KIakaFeMNQL5x6OT6EW4cVkHUr+ATb78IL5xaxFd/eaLs45ZSKjyS0NWukEIIIbhwcxiySHDNdqu7JQBMLKXyHjdj1zgMhb1OzKZbaxwYl4xGcP2OPpw3FKraioQVZ/2mPeR+OOrt+Lbm1VAkAXdeNYpXp2JlA9MvnVnCZ35sGY9qxmHIdisRAnikzl1CO/eVNYh3XzmCN+wZwD88c8o5puom7vgfT+HR/ZZf0+oFxFVDIf/hhu34yO3nwa9IjnE4u5if0sqUw2DYgz77PRzocuVw35t24dv3XgdCSNVqc5Zq/e6rrF6WmyPdG29ww9penJwvnW34Rw/sxROvzuBtlw478bFyMOXgl8WO7m7AjcMqIYTg6q09mI5lHJfI2aUUjkzH8Xc/OwZKKZZTatkCuG7mzRcO4U/euBMAMNZrLVpnFvOVw7TtJx8MedEX5MqhEKdlc7p0ptfx2QQEAlw2FsHlY9GmtVpod3baxWylCjGzuoGFpIo/fdMu/M/3XAlZrLwsBjwSQl2QJMErpNcAW9gml9PYORDEabuU/sh0HPvOLNnKoTu7YNbKQNADryzgbIFx2HdmCZsjVnyCqa9urY4uRbSKcjg+l8SWXj88koiH/+R6iF3uUmJs7fPnTXFzw2ZFD6yilmZT2Iu01vhRna2EK4c1wBrqsYXt1Lz12ysL+OZzp7Gc0iqmsXIsBTbW489TDqZJ8dzxBVy/sx+EEKdNd7dWR5eCdbgtFZCmlOLYbNxxoUii0NFuj9XgkUSM9vhLKgc2N2RVxiHi7fgZGZ396hrEGDMOS5a//PRCEiGvhHdfOYpvPHsKhFgZOZzKjPX6nfcQsJTXUkrDDTut8Z/sy7qJGwcHUSAIeaWiSWCmSfHJRw/jtZkE3n3laJmzu5udZTocOMYhWPvn7MO37VlzM8SNAjcOa2Aw5IEiCZhgymEhhW19AXz0jvMR8kr4ylMnsGeoO+cPrIaxHh9+fXIRzx6fxw/2TTo96m/YZRmHt1+2GUGP6ASvORZRv1xU5/Dgc6fwD8+cwh/euB1//PodrbmwNmfHQBDPnViAadK8DK61KIcruqCbMjcOa0AQCEajPpy10zBPLyRx8UgEPkXER24/Dx+8dTdELuerMtbrRzyr4+PfP4DTCykQYrUoGbYzbCI+Ge+6gu+CC4n45Dy3EqUU//j8GVyxJYq/etsF3JVUhh0DAWQ0E+dW0nmzVphx6AtyV7AbHnNYI6O9fpxdTEMzTEwspbG1L/dhk0Wh63PLa4EpgtMLKbxhzwAozakGTnmiPiWvL9W+M8sYn03g7teNccNQARaLKXQtzSUy6A0oVbOUug2uHNbIWI8P+yeWcW45Dd2k2Nrhc6AbAYvdRP0y/v73rsK+M0u4uIYpe91OxC/ntTz/7q/PwK+I+M1LN7fwqtofd/ubN+wZcI7PxbNFQ3043DismbFeP5ZTGg5OWl1Dt3HjsGq29PnhkQT8+2u3wKeIPIhfIxGf7FRIGybFo/un8NZLhp2WGZzSDAQ96PHLODodzzs+F8+uKt7QLfBP0xphu96nx+cAANv6eNB0tQQ9En7yf9q73xi5qjKO498fu22h/0u7SwoUWqSxrUQKbRQUIQFJrNFggolUQyu8UPwTRN9YjIm+8QVECYGaYCMYUINoEC2KEEQE8Q/a2oZSakuhjbQUu8VSulvZsu3ji3umDJ3dsjN7Z++dzu+TTHbmzMzJc5/cnWfOnXvP+erFbb3eQCOmnjSG1w68eWQxn76Dh1jUBisKjpQk5s+czKZXXn9be09vP4vO8P/v0XyQrUGVC+EeWLeT8WM7/M2jQWdOn0Cnj/XWZer4MQwcDg4cPHRkOojKMrZ2bPNnTmbzK/sZSFOeR4RHDkPwyKFB7+qayNzuiZw69SSuvuBM/xBoo2ZK1YVwleJwlovDsMyfOZn+gcNsf7WPs7sn0ds/wBtvHnZxGISLQ4MmjOvk0a9dUnQY1oamnJSdcrnvQFYcJnjkOmyVtRqe27Wf3v5DR1Z1dP5quTiYtZi31nQ4yLY9fczpmuCR6zCd3T2RzhPEr9bt5PHNu9+a+beOq6PbhQ/2mrWY6mm7t+3pY84MX40/XOM6Ozi7eyJ/+NduOqSGJt1rFy4OZi2mMnLo6e1nx94D/jG6TpXFfFYsmXdkHq9uF4caPqxk1mKmjR9L5wli9fqXORz+MbpeHz93Jr39Ayy7cDZXLDyNp7e96vXeB+HiYNZiThzTwbUXzWHVk9lStR451OfSeadw6bxTgOxw0sd8ZfmgfFjJrAXd8OG5nD4tu9ZmtouDNYFHDmYtaPzYTlZ++nye3NJz5AdqszyVZuQg6SOSNkvaKmlF0fGYld3CWVO5/rK5RYdhx6lSFAdJHcD3gSXAAmCppAXFRmVm1r5KURyA9wFbI+LFiDgI/Ay4ouCYzMzaVlmKw2nAS1WPd6S2t5H0OUlrJK3p6ekZteDMzNpNWYrDYNf+R01DxKqIWBwRi7u6ugZ5i5mZ5aEsxWEHMKvq8enAywXFYmbW9spSHP4BzJU0R9JY4CpgdcExmZm1rVJc5xARA5K+DDwCdAB3RcTGgsMyM2tbpSgOABHxEPBQ0XGYmRmosthFq5H0PyDP0cUUYF+J+5sB7MmxvzzjK3vuKvLKYStsbzvtf63QXxnyV4nhzIh45zN6IqIlb0BPzv2tKnl/a8oaX9lzl3cOW2F722n/a5H+Cs9fvTGU5QfpRryWc38Plry/vOUZn3NXrv6a1Weeyp7DdstfjVY+rLQmIhYXHcdoabftbQbnsHHO3ciUIX/1xtDKI4dVRQcwytpte5vBOWycczcyZchfXTG07MjBzMyap5VHDmZm1iQuDgWRNEvS45I2Sdoo6Sup/WRJj0p6Pv2dltqnp9f3SlpZ1c8kSeurbnsk3VrUdo2mvHKYnlsqaYOkZyQ9LGlGEds0WnLO3adS3jZKurmI7RltDeTvcklr0z62VtKlVX0tSu1bJd0mabC55kZfnqdX+VbXaWUzgfPT/UnAFrK1LG4GVqT2FcBN6f4E4CLgOmDlMfpdC1xc9Pa1Ug7JLgbdDcxIj28Gvl309rVI7qYD/wa60uO7gcuK3r4S5u884NR0/xxgZ1VffwcuJJuA9HfAkqK3L6K1T2VtaRGxKyL+me7vBzaRTVN+Bdk/GOnvJ9Jr+iLiKeCNofqUNBfoBv7UxNBLI8ccKt0mpG9tkznOJ37MMXdnAVsiojKH/u+BK5scfuEayN+6iKjsUxuBEyWNkzQTmBwRf42sUtxTeU/RXBxKQNJssm8WTwOnRMQuyHZAsg/74VoK3Jd2srYykhxGxJvAF4ANZEVhAXBnE8MtlRHuf1uBeZJmS+ok+2Cb9Q7vOa40kL8rgXUR0U9WUHZUPTfoWjZFcHEomKSJwP3ADRHx+gi7uwq4d+RRtZaR5lDSGLLicB5wKvAMcGOuQZbUSHMXEXvJcncf2Yh1OzCQZ4xlVm/+JL0HuAn4fKVpkJeV4sudi0OB0ofS/cBPI+KXqfk/aahJ+rt7mH2dC3RGxNqmBFtSOeVwIUBEvJBGXT8HPtCkkEsjr/0vIh6MiPdHxIXAZuD5ZsVcJvXmT9LpwAPAsoh4ITXvIFu/pqI0a9m4OBQkHdu+E9gUEbdUPbUaWJ7uLwd+Pcwul9Jmo4Ycc7gTWCCpMhnZ5WTHkI9bee5/krrT32nAF4Ef5htt+dSbP0lTgd8CN0bEnysvToee9ku6IPW5jOH/zzdX0b+It+uN7MyPIDuEsT7dPkp29sdjZN++HgNOrnrPduC/QC/ZN44FVc+9CMwrertaNYdkZ+FsSn09CEwvevtaKHf3As+l21VFb1sZ8wd8E+ireu16oDs9txh4FngBWEm6OLnom6+QNjOzGj6sZGZmNVwczMyshouDmZnVcHEwM7MaLg5mZlbDxcGsCSRdJ2lZHa+fLenZZsZkVo/OogMwO95I6oyIO4qOw2wkXBzMBpEmU3uYbDK188imZF4GzAduASYCe4DPRsQuSX8E/gJ8EFgtaRLQGxHflbQQuAMYT3ah07URsVfSIuAu4ADw1Ohtndk782Els6G9G1gVEe8FXge+BNwOfDIiKh/s36l6/dSIuCQivndUP/cAX0/9bAC+ldp/BFwf2ZxEZqXikYPZ0F6Kt+bB+QnwDbKFWh5Ni3V1ALuqXn/f0R1ImkJWNJ5ITXcDvxik/cfAkvw3wawxLg5mQzt6bpn9wMZjfNPvq6NvDdK/WWn4sJLZ0M6QVCkES4G/AV2VNklj0vz8Q4qIfcBeSR9KTVcDT0TEa8A+SRel9s/kH75Z4zxyMBvaJmC5pB+QzbJ5O/AIcFs6LNQJ3Eq27OOxLAfukDSebPbca1L7NcBdkg6kfs1Kw7Oymg0ina30m4g4p+BQzArhw0pmZlbDIwczM6vhkYOZmdVwcTAzsxouDmZmVsPFwczMarg4mJlZDRcHMzOr8X8djlb+iR/ouQAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [], + "source": [ + "first_september_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", + " for y in range(1991,\n", + " sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "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": 14, @@ -1138,1036 +1245,85 @@ "outputs": [ { "data": { - "text/html": [ - "
\n", - "\n", - "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_nameperiod
020201474624260266467410FRFrance2020-03-30/2020-04-05
1202013773715268947411814FRFrance2020-03-23/2020-03-29
22020127812357901045612816FRFrance2020-03-16/2020-03-22
3202011710198756812828151119FRFrance2020-03-09/2020-03-15
420201079011669111331141018FRFrance2020-03-02/2020-03-08
52020097136311054416718211626FRFrance2020-02-24/2020-03-01
6202008710424770813140161220FRFrance2020-02-17/2020-02-23
720200778959657411344141018FRFrance2020-02-10/2020-02-16
820200679264692511603141018FRFrance2020-02-03/2020-02-09
920200578505631410696131016FRFrance2020-01-27/2020-02-02
102020047799158311015112915FRFrance2020-01-20/2020-01-26
1120200375968410078369612FRFrance2020-01-13/2020-01-19
12202002765344530853810713FRFrance2020-01-06/2020-01-12
1320200179835701912651151119FRFrance2019-12-30/2020-01-05
142019527794152461063612816FRFrance2019-12-23/2019-12-29
1520195175823367579719612FRFrance2019-12-16/2019-12-22
16201950764244276857210713FRFrance2019-12-09/2019-12-15
17201949766214540870210713FRFrance2019-12-02/2019-12-08
1820194875542338377018511FRFrance2019-11-25/2019-12-01
192019477753650581001411715FRFrance2019-11-18/2019-11-24
202019467263813163960426FRFrance2019-11-11/2019-11-17
2120194574492261563697410FRFrance2019-11-04/2019-11-10
2220194475728362778299612FRFrance2019-10-28/2019-11-03
2320194374834275169177410FRFrance2019-10-21/2019-10-27
24201942762793989856910713FRFrance2019-10-14/2019-10-20
252019417413020306230639FRFrance2019-10-07/2019-10-13
262019407421122186204639FRFrance2019-09-30/2019-10-06
272019397313713104964528FRFrance2019-09-23/2019-09-29
282019387307814164740528FRFrance2019-09-16/2019-09-22
2920193779701621778102FRFrance2019-09-09/2019-09-15
....................................
15011991267176081130423912312042FRFrance1991-06-24/1991-06-30
15021991257161691070021638281838FRFrance1991-06-17/1991-06-23
15031991247161711007122271281739FRFrance1991-06-10/1991-06-16
1504199123711947767116223211329FRFrance1991-06-03/1991-06-09
1505199122715452995320951271737FRFrance1991-05-27/1991-06-02
1506199121714903897520831261636FRFrance1991-05-20/1991-05-26
15071991207190531274225364342345FRFrance1991-05-13/1991-05-19
15081991197167391124622232291939FRFrance1991-05-06/1991-05-12
15091991187213851388228888382551FRFrance1991-04-29/1991-05-05
1510199117713462887718047241632FRFrance1991-04-22/1991-04-28
15111991167148571006819646261834FRFrance1991-04-15/1991-04-21
1512199115713975978118169251832FRFrance1991-04-08/1991-04-14
1513199114712265768416846221430FRFrance1991-04-01/1991-04-07
151419911379567604113093171123FRFrance1991-03-25/1991-03-31
1515199112710864733114397191325FRFrance1991-03-18/1991-03-24
15161991117155741118419964271935FRFrance1991-03-11/1991-03-17
15171991107166431137221914292038FRFrance1991-03-04/1991-03-10
1518199109713741878018702241533FRFrance1991-02-25/1991-03-03
1519199108713289881317765231531FRFrance1991-02-18/1991-02-24
1520199107712337807716597221529FRFrance1991-02-11/1991-02-17
1521199106710877701314741191226FRFrance1991-02-04/1991-02-10
1522199105710442654414340181125FRFrance1991-01-28/1991-02-03
15231991047791345631126314820FRFrance1991-01-21/1991-01-27
15241991037153871048420290271836FRFrance1991-01-14/1991-01-20
15251991027162771104621508292038FRFrance1991-01-07/1991-01-13
15261991017155651027120859271836FRFrance1990-12-31/1991-01-06
15271990527193751329525455342345FRFrance1990-12-24/1990-12-30
15281990517190801380724353342543FRFrance1990-12-17/1990-12-23
1529199050711079666015498201228FRFrance1990-12-10/1990-12-16
15301990497114302610205FRFrance1990-12-03/1990-12-09
\n", - "

1531 rows × 11 columns

\n", - "
" - ], "text/plain": [ - " week indicator inc inc_low inc_up inc100 inc100_low \\\n", - "0 202014 7 4624 2602 6646 7 4 \n", - "1 202013 7 7371 5268 9474 11 8 \n", - "2 202012 7 8123 5790 10456 12 8 \n", - "3 202011 7 10198 7568 12828 15 11 \n", - "4 202010 7 9011 6691 11331 14 10 \n", - "5 202009 7 13631 10544 16718 21 16 \n", - "6 202008 7 10424 7708 13140 16 12 \n", - "7 202007 7 8959 6574 11344 14 10 \n", - "8 202006 7 9264 6925 11603 14 10 \n", - "9 202005 7 8505 6314 10696 13 10 \n", - "10 202004 7 7991 5831 10151 12 9 \n", - "11 202003 7 5968 4100 7836 9 6 \n", - "12 202002 7 6534 4530 8538 10 7 \n", - "13 202001 7 9835 7019 12651 15 11 \n", - "14 201952 7 7941 5246 10636 12 8 \n", - "15 201951 7 5823 3675 7971 9 6 \n", - "16 201950 7 6424 4276 8572 10 7 \n", - "17 201949 7 6621 4540 8702 10 7 \n", - "18 201948 7 5542 3383 7701 8 5 \n", - "19 201947 7 7536 5058 10014 11 7 \n", - "20 201946 7 2638 1316 3960 4 2 \n", - "21 201945 7 4492 2615 6369 7 4 \n", - "22 201944 7 5728 3627 7829 9 6 \n", - "23 201943 7 4834 2751 6917 7 4 \n", - "24 201942 7 6279 3989 8569 10 7 \n", - "25 201941 7 4130 2030 6230 6 3 \n", - "26 201940 7 4211 2218 6204 6 3 \n", - "27 201939 7 3137 1310 4964 5 2 \n", - "28 201938 7 3078 1416 4740 5 2 \n", - "29 201937 7 970 162 1778 1 0 \n", - "... ... ... ... ... ... ... ... \n", - "1501 199126 7 17608 11304 23912 31 20 \n", - "1502 199125 7 16169 10700 21638 28 18 \n", - "1503 199124 7 16171 10071 22271 28 17 \n", - "1504 199123 7 11947 7671 16223 21 13 \n", - "1505 199122 7 15452 9953 20951 27 17 \n", - "1506 199121 7 14903 8975 20831 26 16 \n", - "1507 199120 7 19053 12742 25364 34 23 \n", - "1508 199119 7 16739 11246 22232 29 19 \n", - "1509 199118 7 21385 13882 28888 38 25 \n", - "1510 199117 7 13462 8877 18047 24 16 \n", - "1511 199116 7 14857 10068 19646 26 18 \n", - "1512 199115 7 13975 9781 18169 25 18 \n", - "1513 199114 7 12265 7684 16846 22 14 \n", - "1514 199113 7 9567 6041 13093 17 11 \n", - "1515 199112 7 10864 7331 14397 19 13 \n", - "1516 199111 7 15574 11184 19964 27 19 \n", - "1517 199110 7 16643 11372 21914 29 20 \n", - "1518 199109 7 13741 8780 18702 24 15 \n", - "1519 199108 7 13289 8813 17765 23 15 \n", - "1520 199107 7 12337 8077 16597 22 15 \n", - "1521 199106 7 10877 7013 14741 19 12 \n", - "1522 199105 7 10442 6544 14340 18 11 \n", - "1523 199104 7 7913 4563 11263 14 8 \n", - "1524 199103 7 15387 10484 20290 27 18 \n", - "1525 199102 7 16277 11046 21508 29 20 \n", - "1526 199101 7 15565 10271 20859 27 18 \n", - "1527 199052 7 19375 13295 25455 34 23 \n", - "1528 199051 7 19080 13807 24353 34 25 \n", - "1529 199050 7 11079 6660 15498 20 12 \n", - "1530 199049 7 1143 0 2610 2 0 \n", - "\n", - " inc100_up geo_insee geo_name period \n", - "0 10 FR France 2020-03-30/2020-04-05 \n", - "1 14 FR France 2020-03-23/2020-03-29 \n", - "2 16 FR France 2020-03-16/2020-03-22 \n", - "3 19 FR France 2020-03-09/2020-03-15 \n", - "4 18 FR France 2020-03-02/2020-03-08 \n", - "5 26 FR France 2020-02-24/2020-03-01 \n", - "6 20 FR France 2020-02-17/2020-02-23 \n", - "7 18 FR France 2020-02-10/2020-02-16 \n", - "8 18 FR France 2020-02-03/2020-02-09 \n", - "9 16 FR France 2020-01-27/2020-02-02 \n", - "10 15 FR France 2020-01-20/2020-01-26 \n", - "11 12 FR France 2020-01-13/2020-01-19 \n", - "12 13 FR France 2020-01-06/2020-01-12 \n", - "13 19 FR France 2019-12-30/2020-01-05 \n", - "14 16 FR France 2019-12-23/2019-12-29 \n", - "15 12 FR France 2019-12-16/2019-12-22 \n", - "16 13 FR France 2019-12-09/2019-12-15 \n", - "17 13 FR France 2019-12-02/2019-12-08 \n", - "18 11 FR France 2019-11-25/2019-12-01 \n", - "19 15 FR France 2019-11-18/2019-11-24 \n", - "20 6 FR France 2019-11-11/2019-11-17 \n", - "21 10 FR France 2019-11-04/2019-11-10 \n", - "22 12 FR France 2019-10-28/2019-11-03 \n", - "23 10 FR France 2019-10-21/2019-10-27 \n", - "24 13 FR France 2019-10-14/2019-10-20 \n", - "25 9 FR France 2019-10-07/2019-10-13 \n", - "26 9 FR France 2019-09-30/2019-10-06 \n", - "27 8 FR France 2019-09-23/2019-09-29 \n", - "28 8 FR France 2019-09-16/2019-09-22 \n", - "29 2 FR France 2019-09-09/2019-09-15 \n", - "... ... ... ... ... \n", - "1501 42 FR France 1991-06-24/1991-06-30 \n", - "1502 38 FR France 1991-06-17/1991-06-23 \n", - "1503 39 FR France 1991-06-10/1991-06-16 \n", - "1504 29 FR France 1991-06-03/1991-06-09 \n", - "1505 37 FR France 1991-05-27/1991-06-02 \n", - "1506 36 FR France 1991-05-20/1991-05-26 \n", - "1507 45 FR France 1991-05-13/1991-05-19 \n", - "1508 39 FR France 1991-05-06/1991-05-12 \n", - "1509 51 FR France 1991-04-29/1991-05-05 \n", - "1510 32 FR France 1991-04-22/1991-04-28 \n", - "1511 34 FR France 1991-04-15/1991-04-21 \n", - "1512 32 FR France 1991-04-08/1991-04-14 \n", - "1513 30 FR France 1991-04-01/1991-04-07 \n", - "1514 23 FR France 1991-03-25/1991-03-31 \n", - "1515 25 FR France 1991-03-18/1991-03-24 \n", - "1516 35 FR France 1991-03-11/1991-03-17 \n", - "1517 38 FR France 1991-03-04/1991-03-10 \n", - "1518 33 FR France 1991-02-25/1991-03-03 \n", - "1519 31 FR France 1991-02-18/1991-02-24 \n", - "1520 29 FR France 1991-02-11/1991-02-17 \n", - "1521 26 FR France 1991-02-04/1991-02-10 \n", - "1522 25 FR France 1991-01-28/1991-02-03 \n", - "1523 20 FR France 1991-01-21/1991-01-27 \n", - "1524 36 FR France 1991-01-14/1991-01-20 \n", - "1525 38 FR France 1991-01-07/1991-01-13 \n", - "1526 36 FR France 1990-12-31/1991-01-06 \n", - "1527 45 FR France 1990-12-24/1990-12-30 \n", - "1528 43 FR France 1990-12-17/1990-12-23 \n", - "1529 28 FR France 1990-12-10/1990-12-16 \n", - "1530 5 FR France 1990-12-03/1990-12-09 \n", - "\n", - "[1531 rows x 11 columns]" + "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAD8CAYAAACyyUlaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHgNJREFUeJzt3X+QVeWd5/H3B5uAOmLAgBEQYSLjBswOBqox4242LhPAZEuwRjM9OkplqMIxmElSW5Vo6RaW8sc4lawbytKVxETUiLJsKNmNRFvc1Li7TANGE0HC0BkJIkh3qolgtujY8N0/znPldNt03+57u++P/ryqTt3T33uew3k8Vn/7+XGeo4jAzMysFKMqfQFmZlb7nEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWckaKn0B5faxj30spk+fXunLMDOrKa+88spvI2LiYMvXXTKZPn06O3furPRlmJnVFEm/KaV8Ud1ckr4habekXZLWSxor6R5Jb0t6LW1fyB1/p6RWSXslLcrF50p6PX23RpJSfIykZ1K8RdL0XJllkvalbVkplTUzs6HRbzKRNAX4O2BeRFwOnAU0pa8fiIg5aXsuHT8rfT8bWAw8JOmsdPzDwApgZtoWp/hy4GhEXAo8ANyfzjUBWAXMBxqBVZLGl1ZlMzMrt2IH4BuAsyU1AOcAh/o4dgnwdER0RsSbQCvQKOkiYFxEbItsqeLHgaW5MuvS/kZgQWq1LAKaI6IjIo4CzZxOQGZmViX6TSYR8TbwbeAAcBh4NyJeSF/fLumXkn6QazFMAd7KneJgik1J+z3j3cpERBfwLnBBH+cyM7MqUkw313iylsMMYDJwrqS/Juuy+gQwhyzJfKdQpJfTRB/xwZbJX+MKSTsl7Wxvb++jNmZmNhSK6eb6c+DNiGiPiPeBHwN/FhFHIuJkRJwCvkc2pgFZ6+HiXPmpZN1iB9N+z3i3Mqkr7Xygo49zdRMRayNiXkTMmzhx0DPbzEastmMn+NIj22g7fqLSl2I1qphkcgC4UtI5aRxjAbAnjYEUXAfsSvubgaY0Q2sG2UD79og4DByXdGU6zy3As7kyhZla1wMvpXGV54GFksanFtLCFDOzMlqzdR879new5sV9lb4Uq1H9PmcSES2SNgI/B7qAV4G1wPclzSHrdtoP3JqO3y1pA/BGOn5lRJxMp7sNeAw4G9iSNoBHgScktZK1SJrSuTok3QfsSMfdGxEdpVTYzE677O4tdHad+uDnJ1sO8GTLAcY0jGLv6msqeGVWa1Rv74CfN29e+KFFs+K0HTvB6uf28MLudzjx/inGjh7Fotkf564vfpJJ542t9OXZMJL0SkTMG2x5r81lNoJNGjeW88Y00Nl1ijENo+jsOsV5YxqcSGzA6m45FTMbmN++18lN8y/hxsZpPLX9AO0ehLdBcDdXTtuxE9y+/lUevPEK/2VmZiOKu7nKyDNazMwGx91ceEaL2UC5FW89uWUCvPzNq7l2zmTGjs7+c4wdPYolcybz8reurvCVmVUnt+KtJ7dM8IwWs2K5FW9n4pZJUpjRsukrV3HT/Etof6+z0pdkVnXcirczccskeeTm05MYVi+9vIJXYla93Iq3M3EyMbMB8XMp1hs/Z2JmZn7OxMzMKs/JxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnEzMxKVlQykfQNSbsl7ZK0XtJYSRMkNUvalz7H546/U1KrpL2SFuXicyW9nr5bI0kpPkbSMyneIml6rsyy9G/sk7SsfFU3M7Ny6TeZSJoC/B0wLyIuB84CmoA7gK0RMRPYmn5G0qz0/WxgMfCQpLPS6R4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictMzOrDsV2czUAZ0tqAM4BDgFLgHXp+3XA0rS/BHg6Ijoj4k2gFWiUdBEwLiK2Rba65OM9yhTOtRFYkFoti4DmiOiIiKNAM6cTkJmZVYl+k0lEvA18GzgAHAbejYgXgAsj4nA65jAwKRWZAryVO8XBFJuS9nvGu5WJiC7gXeCCPs5lZmZVpJhurvFkLYcZwGTgXEl/3VeRXmLRR3ywZfLXuELSTkk729vb+7g0MzMbCsV0c/058GZEtEfE+8CPgT8DjqSuK9JnWzr+IHBxrvxUsm6xg2m/Z7xbmdSVdj7Q0ce5uomItRExLyLmTZw4sYgqmZlZORWTTA4AV0o6J41jLAD2AJuBwuyqZcCzaX8z0JRmaM0gG2jfnrrCjku6Mp3nlh5lCue6Hngpjas8DyyUND61kBammJmZVZF+X9sbES2SNgI/B7qAV4G1wB8BGyQtJ0s4N6Tjd0vaALyRjl8ZESfT6W4DHgPOBrakDeBR4AlJrWQtkqZ0rg5J9wE70nH3RkRHSTU2M7Oy82t7zczMr+01M7PKczIxM7OSOZmYmVnJnEzMzErQduwEX3pkG23HT1T6UirKycTMrARrtu5jx/4O1ry4r9KXUlH9Tg02M7MPu+zuLXR2nfrg5ydbDvBkywHGNIxi7+prKnhlleGWiVkR3JVhPb38zau5ds5kxo7Ofo2OHT2KJXMm8/K3rq7wlVWGk4lZEdyVYT1NGjeW88Y00Nl1ijENo+jsOsV5YxqYdN7YSl9aRbiby6wP7sqwvvz2vU5umn8JNzZO46ntB2gfwS1XPwFv1oe2YydY/dweXtj9DifeP8XY0aNYNPvj3PXFT47Yv0CtPvkJeLMh5K4Ms+K4m8usH+7KMOufu7nMzMzdXGZmVnlOJmZmVjInEzMzK5mTiZmZlczJxEYsL5FiVj5OJlYWtfiL2UukmJWPnzOxssj/Yl593acqfTl98hIpZuXn50ysJD1/MRdU8y9mL5Fi9mF+zsQqqhaX4fYSKWbl128ykXSZpNdy2zFJX5d0j6S3c/Ev5MrcKalV0l5Ji3LxuZJeT9+tkaQUHyPpmRRvkTQ9V2aZpH1pW1be6lupavUXc2GJlE1fuYqb5l9C+3udlb4ks5rW75hJROwF5gBIOgt4G9gEfBl4ICK+nT9e0iygCZgNTAZelPQnEXESeBhYAfwT8BywGNgCLAeORsSlkpqA+4G/lDQBWAXMAwJ4RdLmiDhacs2tbGpx7apHbj7dml+99PIKXsnQaTt2gtvXv8qDN15R9cndat9AB+AXAL+OiN+kRkVvlgBPR0Qn8KakVqBR0n5gXERsA5D0OLCULJksAe5J5TcCD6ZWyyKgOSI6UplmsgS0foDXbUNoJPxirkW1NCnCat9Ak0kT3X+R3y7pFmAn8B9Ti2EKWcuj4GCKvZ/2e8ZJn28BRESXpHeBC/LxXsp8QNIKshYP06ZNG2CVzOqLZ6tZJRQ9AC/pI8C1wH9LoYeBT5B1gR0GvlM4tJfi0Ud8sGVOByLWRsS8iJg3ceLEM9bBbCSoxUkRVvsGMpvrGuDnEXEEICKORMTJiDgFfA9oTMcdBC7OlZsKHErxqb3Eu5WR1ACcD3T0cS4zO4NanRRhtW0gyeSvyHVxSboo9911wK60vxloSjO0ZgAzge0RcRg4LunKNB5yC/Bsrkxhptb1wEuRPQDzPLBQ0nhJ44GFKWZmffBsNRtuRY2ZSDoH+Dxway78D5LmkHU77S98FxG7JW0A3gC6gJVpJhfAbcBjwNlkA+9bUvxR4Ik0WN9BNjZDRHRIug/YkY67tzAYb2Zn5kkRNtz8BLyZmfkJeDMzqzwnEzMzK5mTiZmZlczJxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnErIa0HTvBlx7ZRlsNvBrZRhYnE7Makn8Vr1k1Gehre82sAvwqXqt2bpmY1QC/iteqnZOJWQ3wq3it2rmby6xGFF7Fe2PjNJ7afoB2D8JbFfGbFs3MzG9aNDOzynMyMTOzkjmZmJlZyfpNJpIuk/Rabjsm6euSJkhqlrQvfY7PlblTUqukvZIW5eJzJb2evlsjSSk+RtIzKd4iaXquzLL0b+yTtKy81Tczs3LoN5lExN6ImBMRc4C5wP8DNgF3AFsjYiawNf2MpFlAEzAbWAw8JOmsdLqHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqsNAu7kWAL+OiN8AS4B1Kb4OWJr2lwBPR0RnRLwJtAKNki4CxkXEtsimkD3eo0zhXBuBBanVsghojoiOiDgKNHM6AZnZCOZ1yqrLQJNJE7A+7V8YEYcB0uekFJ8CvJUrczDFpqT9nvFuZSKiC3gXuKCPc5nZCOd1yqpL0Q8tSvoIcC1wZ3+H9hKLPuKDLZO/thVk3WdMmzatn8szs1rmdcqq00BaJtcAP4+II+nnI6nrivTZluIHgYtz5aYCh1J8ai/xbmUkNQDnAx19nKubiFgbEfMiYt7EiRMHUCUzqzVep6w6DSSZ/BWnu7gANgOF2VXLgGdz8aY0Q2sG2UD79tQVdlzSlWk85JYeZQrnuh54KY2rPA8slDQ+DbwvTDEzG6G8Tll1KqqbS9I5wOeBW3Phvwc2SFoOHABuAIiI3ZI2AG8AXcDKiDiZytwGPAacDWxJG8CjwBOSWslaJE3pXB2S7gN2pOPujYiOQdTTzOqI1ymrPl6by8zMvDaXmVkxPJV4aDmZmNmI4KnEQ8vvMxlB2o6d4Pb1r/LgjVd4sNJGDE8lHh5umYwg/svMRiJPJR4ebpmMAP7LzEYyTyUeHm6ZjACD+cvMg5VWTwpTiTd95Spumn8J7e91VvqS6o5bJiPAYP4yy3eJrb7uU8N4tWbl98jNp2e8rl56eQWvpH45mYwQxT7k5S4xMxsMP7Ro3bQdO8Hq5/bwwu53OPH+KcaOHsWi2R/nri9+0n3MZnXMDy1aWXmw0swGw91c9iFe98jMBsrdXGZm5m4uMzOrPCcTMzMrmZOJmZmVzMlkGPhpcjOrd04mw8ALLJpZvfPU4CHkp8nNbKRwy2QIeelrMxspnEyGkJ8mN7ORwslkiA106WsP1ptZLSoqmUj6qKSNkn4laY+kz0i6R9Lbkl5L2xdyx98pqVXSXkmLcvG5kl5P362RpBQfI+mZFG+RND1XZpmkfWlbVr6qD49Hbp7H6qWXM2vyOFYvvbzbUti98WC9mdWiopZTkbQOeDkivi/pI8A5wNeB9yLi2z2OnQWsBxqBycCLwJ9ExElJ24GvAf8EPAesiYgtkr4C/OuI+FtJTcB1EfGXkiYAO4F5QACvAHMj4uiZrrVWl1PpOVhf4MF6MxsOQ76ciqRxwGeBRwEi4g8R8bs+iiwBno6Izoh4E2gFGiVdBIyLiG2RZbDHgaW5MuvS/kZgQWq1LAKaI6IjJZBmYPGAa1kDPFhvZrWsmG6uPwbagR9KelXS9yWdm767XdIvJf1A0vgUmwK8lSt/MMWmpP2e8W5lIqILeBe4oI9z1R0P1ptZLSsmmTQAnwYejogrgN8DdwAPA58A5gCHge+k49XLOaKP+GDLfEDSCkk7Je1sb2/voyrVze+pNrNaVcxDiweBgxHRkn7eCNwREUcKB0j6HvA/c8dfnCs/FTiU4lN7iefLHJTUAJwPdKT453qU+VnPC4yItcBayMZMiqhTydqOneD29a/y4I1XlK314PdUm1mt6rdlEhHvAG9JuiyFFgBvpDGQguuAXWl/M9CUZmjNAGYC2yPiMHBc0pVpPOQW4NlcmcJMreuBl9K4yvPAQknjUzfawhSrOM+6MiuOp7uPDMUup/JV4EdpJte/AF8G1kiaQ9bttB+4FSAidkvaALwBdAErI+JkOs9twGPA2cCWtEE2uP+EpFayFklTOleHpPuAHem4eyOiY3BVLQ8vkWI2MPk/vFZf96lKX44NEb9pcYDajp1g9XN7eGH3O5x4/xRjR49i0eyPc9cXP+nBcrMcT3evLX7T4jDzrCuz4ni6+8jiVYMHoTDr6sbGaTy1/QDt7gs2+xD/4TWyOJkMgmddWbkMxazAauI/vEYOj5mYVdDdm17nR9sPcFPjNA9OW0WVOmbilolZBXhWoNUbD8CblVkxz1V4cNrqjZOJWZkV80CrB6et3riby6xMBtp15cHpD6v3CQn1zAPwZmXiB1pL5wkJleMBeLMq4a6rwfOEhNrnMROzMvJrBAbHExJqn1smVtVqrQ/dD7QOjlt1tc8tExt2A1mS3Ev9jxxu1dU2D8DbsCtmkNUrzpoNr1IH4J1MbNgMJEF4ZpTZ8PIS9FYzBjLI6j50s9riAXgbNgNNEH6oz6x2OJnYsBpIgvDMKLPa4TETMzPzmImZmVWek4mZWQ8DeRbKMk4mZmY9+GHZgSsqmUj6qKSNkn4laY+kz0iaIKlZ0r70OT53/J2SWiXtlbQoF58r6fX03RpJSvExkp5J8RZJ03NllqV/Y5+kZeWruplZd5fdvYXpd/yEJ1sOEJEtODn9jp9w2d1bKn1pVa/Ylsl3gZ9GxL8C/hTYA9wBbI2ImcDW9DOSZgFNwGxgMfCQpLPSeR4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictM7Ny8oKTg9dvMpE0Dvgs8ChARPwhIn4HLAHWpcPWAUvT/hLg6YjojIg3gVagUdJFwLiI2BbZFLLHe5QpnGsjsCC1WhYBzRHRERFHgWZOJyAzs7Lyw7KDV0zL5I+BduCHkl6V9H1J5wIXRsRhgPQ5KR0/BXgrV/5gik1J+z3j3cpERBfwLnBBH+fqRtIKSTsl7Wxvby+iSmZmvfOCk4NTzEOLDcCnga9GRIuk75K6tM5AvcSij/hgy5wORKwF1kL2nEkf12Zm1ic/LDs4xbRMDgIHI6Il/byRLLkcSV1XpM+23PEX58pPBQ6l+NRe4t3KSGoAzgc6+jiXmZlVkX6TSUS8A7wl6bIUWgC8AWwGCrOrlgHPpv3NQFOaoTWDbKB9e+oKOy7pyjQeckuPMoVzXQ+8lMZVngcWShqfBt4XppiZmVWRYtfm+irwI0kfAf4F+DJZItogaTlwALgBICJ2S9pAlnC6gJURcTKd5zbgMeBsYEvaIBvcf0JSK1mLpCmdq0PSfcCOdNy9EdExyLqamdkQ8dpcZmbmtbnMzKzynEzMzOpApdcTczIxM6sDlV5PzC/HMjOrYZfdvYXOrlMf/PxkywGebDnAmIZR7F19zbBdh1smZmY1rFrWE3MyMTOrYdWynpi7uczMalxhPbEbG6fx1PYDtFdgEN7PmZiZmZ8zMTOzynMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMieTGlfp9XjMzMDJpOZVej0eMzPwQ4s1q1rW4zEzA7dMala1rMdjZgZOJjWrWtbjMTMDd3PVtGpYj8fMDLw2l5mZ4bW5zMysChSVTCTtl/S6pNck7UyxeyS9nWKvSfpC7vg7JbVK2itpUS4+N52nVdIaSUrxMZKeSfEWSdNzZZZJ2pe2ZeWquJnZcKvn58IG0jK5OiLm9GgGPZBicyLiOQBJs4AmYDawGHhI0lnp+IeBFcDMtC1O8eXA0Yi4FHgAuD+dawKwCpgPNAKrJI0fRD3NzCqunp8LG4oB+CXA0xHRCbwpqRVolLQfGBcR2wAkPQ4sBbakMvek8huBB1OrZRHQHBEdqUwzWQJaPwTXbWY2JEbCc2HFtkwCeEHSK5JW5OK3S/qlpB/kWgxTgLdyxxxMsSlpv2e8W5mI6ALeBS7o41xmZjVjJDwXVmwyuSoiPg1cA6yU9FmyLqtPAHOAw8B30rHqpXz0ER9smQ9IWiFpp6Sd7e3tfVbEzGy4jYTnwopKJhFxKH22AZuAxog4EhEnI+IU8D2yMQ3IWg8X54pPBQ6l+NRe4t3KSGoAzgc6+jhXz+tbGxHzImLexIkTi6mSmdmwKjwXtukrV3HT/Etof6+z3zK1NGDfbzKRdK6k8wr7wEJgl6SLcoddB+xK+5uBpjRDawbZQPv2iDgMHJd0ZRoPuQV4NlemMFPreuClyB6AeR5YKGl86kZbmGJmZjXlkZvnsXrp5cyaPI7VSy/nkZv7f6SjlgbsixmAvxDYlGbxNgBPRcRPJT0haQ5Zt9N+4FaAiNgtaQPwBtAFrIyIk+lctwGPAWeTDbxvSfFHgSfSYH0H2WwwIqJD0n3AjnTcvYXBeDOzelWLA/Z+At7MrMq0HTvB6uf28MLudzjx/inGjh7Fotkf564vfnLIxln8BLyZWZ2pxQF7L/RoZlaFam0hV3dzmZmZu7nMzKzynEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWcmcTKxu1NK7H8zqjZOJ1Y1aeveDWb3xQo9W82rx3Q9m9cYtE6t5L3/zaq6dM5mxo7P/nceOHsWSOZN5+VtXV/jKzEYOJxOrebX47gezeuNuLqsLtfbuB7N64/eZmJmZ32diZmaV52RiZmYlKyqZSNov6XVJr0namWITJDVL2pc+x+eOv1NSq6S9khbl4nPTeVolrZGkFB8j6ZkUb5E0PVdmWfo39klaVq6Km5lZ+QykZXJ1RMzJ9andAWyNiJnA1vQzkmYBTcBsYDHwkKSzUpmHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqkMp3VxLgHVpfx2wNBd/OiI6I+JNoBVolHQRMC4itkU26v94jzKFc20EFqRWyyKgOSI6IuIo0MzpBGRmZlWi2GQSwAuSXpG0IsUujIjDAOlzUopPAd7KlT2YYlPSfs94tzIR0QW8C1zQx7nMzKyKFPucyVURcUjSJKBZ0q/6OFa9xKKP+GDLnP4HswRXSHLvSdrbx/XVgo8Bv630RQyxeq9jvdcP6r+OI61+l5RysqKSSUQcSp9tkjaRjV8ckXRRRBxOXVht6fCDwMW54lOBQyk+tZd4vsxBSQ3A+UBHin+uR5mf9XJ9a4G1xdSlFkjaWcp871pQ73Ws9/pB/dfR9RuYfru5JJ0r6bzCPrAQ2AVsBgqzq5YBz6b9zUBTmqE1g2ygfXvqCjsu6co0HnJLjzKFc10PvJTGVZ4HFkoanwbeF6aYmZlVkWJaJhcCm9Is3gbgqYj4qaQdwAZJy4EDwA0AEbFb0gbgDaALWBkRJ9O5bgMeA84GtqQN4FHgCUmtZC2SpnSuDkn3ATvScfdGREcJ9TUzsyFQd8up1ANJK1LXXd2q9zrWe/2g/uvo+g3wfE4mZmZWKi+nYmZmJXMyGSaSfiCpTdKuXOxPJW1LS8z8D0njUvwjkn6Y4r+Q9LlcmZ+lZWpeS9ukXv65YSfpYkn/S9IeSbslfS3Fy7bsTiWVuX51cQ8lXZCOf0/Sgz3OVfP3sJ/6Vd09HET9Pq/s2cHX0+e/z51r4PcvIrwNwwZ8Fvg0sCsX2wH8u7T/N8B9aX8l8MO0Pwl4BRiVfv4ZMK/S9emlfhcBn0775wH/DMwC/gG4I8XvAO5P+7OAXwBjgBnAr4Gz0nfbgc+QPWe0BbimzupXL/fwXODfAH8LPNjjXPVwD/uqX9Xdw0HU7wpgctq/HHi7lPvnlskwiYh/JJuplncZ8I9pvxn4i7Q/i2y9MyKiDfgdUNXz3SPicET8PO0fB/aQrVZQzmV3KqZc9Rveqx6YgdYxIn4fEf8b6PYmsnq5h2eqX7UaRP1ejfQMIbAbGKvskY5B3T8nk8raBVyb9m/g9MOevwCWSGpQ9qzOXLo/CPrD1LT+T9XQfdCTslWfrwBaKO+yO1WhxPoV1MM9PJN6uYf9qdp7OIj6/QXwakR0Msj752RSWX8DrJT0Clmz9A8p/gOyG7gT+C/A/yV7Zgfgpoj4FPBv03bzsF5xPyT9EfDfga9HxLG+Du0lVvQSOpVShvpB/dzDM56il1gt3sO+VO09HGj9JM0mW6n91kKol8P6vX9OJhUUEb+KiIURMRdYT9avTkR0RcQ3IlvyfwnwUWBf+u7t9HkceIoq6jqRNJrsf+IfRcSPU/hIajYXuj9KWXanospUv3q6h2dSL/fwjKr1Hg60fpKmApuAWyLi1yk8qPvnZFJBhRkgkkYBdwP/Nf18jrKla5D0eaArIt5I3V4fS/HRwH8g6yqruNTMfxTYExH/OfdVOZfdqZhy1a/O7mGv6ugenuk8VXkPB1o/SR8FfgLcGRH/p3DwoO9fpWYejLSNrOVxGHifLPMvB75GNuPin4G/5/RDpNOBvWQDaC8Cl6T4uWQzu35JNmD2XdIMoUpvZLNeIl3ba2n7AtmrBLaStay2AhNyZe4ia43tJTdbhGyywa703YOF/y71UL86vIf7ySaWvJf+v55VZ/fwQ/Wr1ns40PqR/QH7+9yxrwGTBnv//AS8mZmVzN1cZmZWMicTMzMrmZOJmZmVzMnEzMxK5mRiZmYlczIxM7OSOZmYmVnJnEzMzKxk/x+qftJYVEr3+gAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "1992 False\n", + "1993 False\n", + "1994 False\n", + "1995 False\n", + "1996 False\n", + "1997 False\n", + "1998 False\n", + "1999 False\n", + "2000 False\n", + "2001 False\n", + "2002 True\n", + "2003 False\n", + "2004 False\n", + "2005 False\n", + "2006 False\n", + "2007 False\n", + "2008 False\n", + "2009 False\n", + "2010 False\n", + "2011 False\n", + "2012 False\n", + "2013 False\n", + "2014 False\n", + "2015 False\n", + "2016 False\n", + "2017 False\n", + "2018 False\n", + "2019 False\n", + "dtype: bool" + ] + }, + "execution_count": 19, + "metadata": {}, + "output_type": "execute_result" } ], "source": [ - "data" + "yearly_incidence == yearly_incidence.min()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": {