diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index de710cddb58269c98814d00280d718ebebce5a85..0f0f501901eeb7cd1b7fcd62df76f5134415f1be 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1082,6 +1082,295 @@ "Deuxièmement, nous trions les points par période, dans le sens chronologique." ] }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Nous vérifions la cohérence des données. Entre la fin d'une période et le début de la période qui suit, la différence temporelle doit être zéro, ou au moins très faible. Nous laissons une \"marge d'erreur\" d'une seconde.\n", + "\n", + "Ceci s'avère tout à fait juste sauf pour deux périodes consécutives entre lesquelles il manque une semaine.\n", + "\n", + "Nous reconnaissons ces dates: c'est la semaine sans observations que nous avions supprimées !" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "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": 12, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ/vudc7ck5GbfE0ggYUkCBIjsKrJGVEABJ+oAjvjDYcBlGHFAR2VUFEaRGVxg2GQZlU0UEAHZ15BwCdlDyCULuUnIvueu59Tvj+46XV1dW5/lLrn1Ps99bp/q6qrq6qr66luLGGPw8PDw8PAQkenqBnh4eHh4dD944uDh4eHhkYAnDh4eHh4eCXji4OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZGAJw4eHh4eHgl44uDh4eHhkYAnDh4eHh4eCVR1dQOKxdChQ9n48eO7uhkeHh4ePQpvv/32ZsbYMFu+Hkscxo8fj4aGhq5uhoeHh0ePAhGtdsnnxUoeHh4eHgl44uDh4eHhkYAnDh4eHh4eCXji4OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZGAJw4eHh4eHgl44uDRK8AYw6Nzm9Dcluvqpnh49Ah44uDRKzBrxRZc9dB8/PjJJV3dFA+PHgFPHDx6BXa3dAAANu5s6eKWeHj0DHji4NErQERd3QQPjx4FTxw8PDw8PBKwEgciqiOiOUQ0n4gWE9F/hunXEdFaIpoX/p0tPHMtETUS0TIiOktIP4aIFob3bqFwO0dEtUT0YJg+m4jGl/9VPTw8PDxc4cI5tAI4lTF2JIBpAGYQ0fHhvZsZY9PCv78BABFNBjATwBQAMwD8loiyYf5bAVwGYFL4NyNMvxTANsbYRAA3A7ix9Ffz8PDw8CgWVuLAAuwOf1aHf8zwyLkAHmCMtTLGVgJoBHAsEY0CUM8Ym8UYYwDuA3Ce8My94fUjAE4jLyT28PDw6DI46RyIKEtE8wBsBPAsY2x2eOtKIlpARHcT0aAwbQyANcLjTWHamPBaTo89wxjrALADwBBFOy4jogYiati0aZPTC3p4eHh4pIcTcWCM5Rhj0wCMRcAFTEUgIjoIgahpPYCbwuyqHT8zpJuekdtxO2NsOmNs+rBh1oOMPDwSYCaeV4MVm3ZjQdP28jfGw6MbI5W1EmNsO4CXAMxgjG0IiUYewB0Ajg2zNQEYJzw2FsC6MH2sIj32DBFVARgAYGuqN/HwMKAUGeWpN72Mc379etna4uHRE+BirTSMiAaG130AnA7g3VCHwPFZAIvC68cBzAwtkCYgUDzPYYytB7CLiI4P9QkXA3hMeOaS8PoCAC+EegkPDw8Pjy6AyxnSowDcG1ocZQA8xBj7KxHdT0TTEIh/VgH4GgAwxhYT0UMAlgDoAHAFY4wHtLkcwD0A+gB4KvwDgLsA3E9EjQg4hplleDcPDw8PjyJhJQ6MsQUAjlKkX2R45noA1yvSGwBMVaS3ALjQ1hYPDw8Pj86B95D26JW46K7ZeHLB+q5uhodHt4UnDh69Eq8u34wr/jC3q5vh4dFt4YmDR6+Ct3Lw8HCDJw4evQLe397DIx08cfDw6EZoz+Vx1YPzsGLTbntmD48KwhMHD49uhPlrtuPRd9bi6kcWdHVTPHo5PHHw8PDw8EjAEwcPDw8PjwQ8cfDoVejuUVm6d+s8ehM8cfDoFfDWSh4e6eCJg4dHN4KnYR7dBZ449DA0rNqKvW0dXd0MDw+PfRyeOPQgbNzVggtum4WrHpzf1U3xqBC8zsGju8AThx6E5rYg8vmS9Tu7uCU9G91dKe3h0R3giYOHRzeC1zl4dBd44uDRq8BQ3DnSnYVu3LR9Gs1tOdz/5mrsbvX6PA5PHDx6BaiH7cl7Vmt7Ph6btxbf/8si3PvGqq5uSreBJw4evQ49YXfeE9q4L2FHc3vsv4cDcSCiOiKaQ0TziWgxEf1nmD6YiJ4louXh/0HCM9cSUSMRLSOis4T0Y4hoYXjvFqLANYmIaonowTB9NhGNL/+renh0f3Qlx3D1w/Mx/ponu7AFXQ9vrBDBhXNoBXAqY+xIANMAzCCi4wFcA+B5xtgkAM+Hv0FEkwHMBDAFwAwAvyWibFjWrQAuAzAp/JsRpl8KYBtjbCKAmwHcWIZ3KzueXLAea7bu7epmeJSI7rwAdGXLHn67qQtr71pwD/puPDQ6HVbiwALw4PLV4R8DcC6Ae8P0ewGcF16fC+ABxlgrY2wlgEYAxxLRKAD1jLFZLJid90nP8LIeAXAa5yq6E674w1yc8+vXuqx+P3B7D7rd4N/HwXVSfopFcNI5EFGWiOYB2AjgWcbYbAAjGGPrASD8PzzMPgbAGuHxpjBtTHgtp8eeYYx1ANgBYEgxL1RpbNvrZZI9GYz1jAWgJ7RxX4LnHJJwIg6MsRxjbBqAsQi4gKmG7KpNDzOkm56JF0x0GRE1EFHDpk2bbM328CiA+eXWwwF+nERIZa3EGNsO4CUEuoINoagI4f+NYbYmAOOEx8YCWBemj1Wkx54hoioAAwBsVdR/O2NsOmNs+rBhw9I0fZ9A9xO0VQZrtu4tu15ALK4n7A4r8alnvb/F68w04FLsnjA2Ogsu1krDiGhgeN0HwOkA3gXwOIBLwmyXAHgsvH4cwMzQAmkCAsXznFD0tIuIjg/1CRdLz/CyLgDwAuvGWsPv/2URtu9t6+pm7JNYvG4HPvpfL+Lu11c5P/Pa8s24+O45yOf1Q6b7jqbOwxfueBMf/a8Xu7oZ3RK9ZN+VCi6cwygALxLRAgBvIdA5/BXADQDOIKLlAM4If4MxthjAQwCWAHgawBWMsVxY1uUA7kSgpH4fwFNh+l0AhhBRI4CrEFo+dVfc/+Zq/PyZZWUr770Nu7CzxesyAOCDLcHOds7KLc7PXP5/b+OV9zZhtyFabU+jDeVq786WdrzzwbYylbbvItI59LSRUjlU2TIwxhYAOEqRvgXAaZpnrgdwvSK9AUBCX8EYawFwoUN7uwzyoMkZdqlpcebNr+DwMQPwxNdPtrShbFV2W1SKveff7+X3NvUqufJX72nAnFVb8d5PPtnVTenW8JxDEt5D2hHyYlXuxWvh2h3Oefdl3UNhB1fEs6Zv0tPIQbk+8bym7QCAfG/YWZQBvpcieOLgCHnQdMVk600Dt/ycQ+XK7s7IeBNNJ3iFdBKeODhCFit1xRjibdiXB3DlmKKe0Wnl/rbcuauYzcwL727Ab19qLG+DuikKRLSHjJPOgFXn4BFAHjJdsUD3hmEbOcaneFsHirIvE1QTMo5iusaNuxNpX7mnAQDwL6dMLHOruiE855CA5xyKRFdYNfSGgcvX+TLq+wH0HMJabn0SJ7YmzmFB03ac/suXy1txD0MRW5J9Hp44OEKeW5lM52uF54YmifuycrFSynbeZRnqHUSWo6Dgz+vzfOAd4wroTWPDBk8cHCHLIquznd9133lkAQBgV8u+e1pVSfbmRmul4Ga/mu4tSS334pRx4Bw8xE2J7ycOTxwcIc+t6mzX2ZPu06asFYqOyb8fQ+9SOvKxYiIOPe2UvErA90ESnjgUiaqM77qKoBTTS8P87inkoOw6h/C/SYfTm4ilDZ7BiuBXOEd0J86hNwzg8nMOQYntuXzF+u+CW9/AF+94M5a2sGkHGjfuci6jUmIlHxbCDE4gfTdF6N4C2G6MbBcopDn2ZflxwWqkQu/Y2pHHY/PW2TMWgYbVyRhGnwkPh1p1w6dSleXKQby9ehtWbt6DC44Zm7j3mxcbsWVPECByb1sucb9QlxepFOC5qAieODhCHjTlGkLFLIL7MG1ASQcAGvpFjIX13T8vLK54xkprX6q63PKdf+sbAKAlDhzn/uZ1fV1+QYx0Ur4rCvBiJUdUKrZSMeX0Ds4hnn73ayvx/qako5YrOnLufbZ2ezOWrNsJANiyu7WQ3hndXk7akxEK29Hso/6awKT/Hp44OCPhId2Fw2gfpg3KxTGfZ/jRX5fgPMPuFzB/k/a8wdBfwkk3vICzb3kVQHAGQlR+5SF+24ZVW/HLZ98ruqxi6czu1n3XVFqLXhCaJi08cXBEQvxTLs6hiGf2Bc5hx952/OSvS9DWEV+0I1PW6B35lW7R0nEbIooNsf7ehohb6cx+JwIuuG0Wbnl+eUllFIOpP3wGTy5YX3S9PRER59Dz51a54ImDI5KcQ9dhH6ANuOHpd3HnayvxxPy4clh10Hs5lNPtKcRKOnRmv5ejrlL0I1f8YW7pDehBEPv7gTkfYOn6nV3XmG4Cr5AuEuWypimmnNw+QB3acwHHIL+LiQuwvbbpdkfOXaykL79n9fu+7CxZbhTmIQOueTQwWEhrYbavwXMOjqiYQrqIZ8p5Cl1XQdt/itDJrm9rIrQdRfSZfCZ1Z9JklVlsWnjakB49f2aVD1biQETjiOhFIlpKRIuJ6Jth+nVEtJaI5oV/ZwvPXEtEjUS0jIjOEtKPIaKF4b1bKOR7iaiWiB4M02cT0fjyv2qJqIzKoddDt4CV+3Ce9iI4h9aO0rmNrkSmm7EOO/a2Y+XmPV3dDCX4ENsX9Hnlggvn0AHg3xhjhwE4HsAVRDQ5vHczY2xa+Pc3AAjvzQQwBcAMAL8lomyY/1YAlwGYFP7NCNMvBbCNMTYRwM0Abiz91cqLhJ9DF5qy7ssoePQKaa7iHFOuYritlva441h3+1a2d+osnwxXnH3Lq/jEL17q6mYo4f0ckrASB8bYesbY3PB6F4ClAMYYHjkXwAOMsVbG2EoAjQCOJaJRAOoZY7NYwP/fB+A84Zl7w+tHAJxG3WxkJ42V/CgqBbr+KyU4pmliF6OQbumQiEMnfHOVaEwnLuuwmOd2rxkU+I/Y8PU/voM7X13RCa2Jw8/mJFLpHEJxz1EAZodJVxLRAiK6m4gGhWljAKwRHmsK08aE13J67BnGWAeAHQCGpGlbpZGwViqbzqF3D0vdHiCmcyhDF+Uc/BxyeYZ731hV+N3SHn+mM3aVqip09dpeqQsjvBSNJ+avw0+eXNrp9RaO4O30mrsvnIkDEe0H4E8AvsUY24lARHQQgGkA1gO4iWdVPM4M6aZn5DZcRkQNRNSwadMm16aXBZWK9dPT2diHGtZgRQmeyzJKYRhNhNZFqvTgW2vww8cXF34nxEpFt8wdqvGgk4Pb5OO6mEliWA2POHyAwghOxIGIqhEQht8zxh4FAMbYBsZYjjGWB3AHgGPD7E0AxgmPjwWwLkwfq0iPPUNEVQAGANgqt4MxdjtjbDpjbPqwYcPc3rBC8IMo6IPvPLLA6rlcXNnueQsExfCMS3m7WuIhJp5aGHcE66pvriNsNpNmHefw82eWxX53x6F8xysrymJ+7ArxvA+PAC7WSgTgLgBLGWO/FNJHCdk+C2BReP04gJmhBdIEBIrnOYyx9QB2EdHxYZkXA3hMeOaS8PoCAC+wbrb6dicnuK7Ai8s24kt3vok9gpcyl+PvLOPJdCr2vhwjoRjx3S0vxHfYncI5KGrRtV02tZVh4sK62fRK4Pq/LcVfKhQ9V8YX73gTv3t9ZafU1ZPg4gR3EoCLACwkonlh2ncBfIGIpiGYM6sAfA0AGGOLieghAEsQWDpdwRjj/PnlAO4B0AfAU+EfEBCf+4moEQHHMLO010qP9lwezy3ZgBlTRyonVaX8HHoKrnt8MVZv2Yumbc04ZGR/AEBrqLCtKeORqYUwBkIHl8NaqSwEpjM2simc/+zWSvp7uTxDVXgmSXdRXMsEq7PMSt94f4u2Db0ZVuLAGHsNap3A3wzPXA/gekV6A4CpivQWABfa2lJJ3P7KCvz8mWW47R+PxoypoxL3kyG7y+UhXZZiKg7eTnFB4n4ANVVFEAfNe6vYe9c+KrUvbYtkp1grKdJ0i6RNrGR6n448Q1VoYN5dxqBsUda/1rw85fIMyz7chcmj68vWhhTxGfd5eA/pENvCQ1FWb9mrzlAxD+luMjMtUJ1FrCMOzy3ZgM1CqGtjuZr0ci9Y5dgRdtUiqmMQ7NZKeurQHZ29WiXT4b4W4nDzs+/h7Ftexbsfli8OUnfsl66CJw4hBvSpBqCPe9/bdQ4c4txpDa15agXi0NyWw1fva8DFd80prnx+XGMsLd2z6nulo6uslXSEzco5GO4VE06k0pDFZDZT3PlN2wEAH+5oKVsbumG3dBk8cQjBd7/ipHls3lqMv+ZJNG1LchNd4SE9bnAfAMDhYwaUp/IU4PM0p+AcqgWdA7+/eos5TIL2tSOlQ5TkKD4xZSuLzqETdpUqAqfnHGzWSvrVNVeGKLXFQtfutHo9flRvOT+L1zlE8MQhBB+v4nz6yztrAQDvbdilGICdP4imjg6IwqTh+3V63VxJL7Ld/Dqts1XDqq14+b1NYbnxe0onMEt5LvO5HOK7bsc5lKKQ7sJFUCe6SdsiTvzKGYjSi5Ui+JDdIaKFTm2pVInYSkvW7URzu7sZKJ8EnT2Af/NiYyFgmrjrM+7UDeVdcNssa53FvGHFrZU6oduVxLFIPweTKWu5I/u+8t4mjBvcFxOG9rPmzTGmXHhkImhrId+UlJPQdSFD1e3giUMIPjDF6cQnV3N7Dht3tkr5S6+TH0XJkc8zZAzbcD6fO3v8ik5TpjVl9ZY9GNyvpqS6VAHQbH0diZVK0znoPIqjMrpm5dDutLuRzuHiuwMdk8sZCLpmJ/R6jsSvFFFQgiB5zqEAL1YKwceEarN15R/ewbmSF3AlFoo2i0dowUGsC8dvTsM5vLliCz7+85fwUEMQPqtY0/lIIS1ShyILE8vtIRpp1eKk9ZAuxVqpCzWvemIn/TaUsbetAxt3toTl6fO9unwTPnXLq9qQ7XKdvfL8bA08cQjBx4hrDPxSF5tL73krkWY7P4BPqq7c2+ic05ZvDOIrvfNBukNqXmvcLJUf/5+uber09Tua8ery0mNxdYrOQZlWfp1DV1orye3euLMF7bm8QnSrb+OXf/cW5jftAGAWs37vz4uweN1OrNNEhJVFUu98sN3Y9t4ETxxC8AHmGvitlKnFGMPz725MpLdZiUP0fFdBnNcRtxUJZNLqQx6du1aZHhMrlbgsn/CzF9C0zR4uWvz0R+0/MHG/U3Q9SoW0Ouujc5vUNxzgEqW2UhDHUHNbDsf+9Hl899GFVl+iz/zqNZzy8xcBAHNWRqHXTESSm1nrNl77wqmKlYInDiH4GHG1vCllndjTllOm28RK3YFzUCn/GGORk1z4CsW2kUn/g/KLLKwEZB1CqHQWdETpztdWFl1mOTiH7XvbEpFrXSBubvjzf1+yITFm5CYuXLsDq0InVfHzmL5LbXWwxJ158yt4csH6xH1vnaSHJw4hIoW0K+dQ/KDaqXG0s8mBCwO5C8dzPiZWilAwKyxxskV6FXU95mdLqjoGlWFA54iVVMS3/PWIXGqxxU/70bP49K9eS/2cuFs3LfKmhXs/wXvalK+OxwgBcN0TixP3PeOghycOISKxkuMDJQwq3W7LtgjwXXlX7nbipqzCJA//84lfvEI6PXhd5TQSUHMOle93VRWV+N7NGu41LRo3pj/LQ7cgpzEXzwrE2yQa4sEFdfm8WEkPTxxC8IEobhhNC9yj76hl5S5wnRzJ57reWimmcxDSVbGXygXbolyJ7siqOIdOsVZyS5OxbU9bqvMPWgTOoTOCsorf0NlaybHDzaHJo2tVvWmttlo7cvjp35b2CqsmTxxC8CHSOUdXqwekbZxGUqXuYWkizjUujit5JxY+rhNfGR+tsFipq6Ba1OR+PurHz+L7jyXFJjqI3Gup3TZ/jd3CJ29ZpBlLjmrTUCqGwVdyDikHzUMNTbj9lRX4n+feS/WcDY0bd3U7H4teTxyWrNuJCdc+WbBmcTVlLQXaWDmWwdEdOIf4AI5RBwClK82LkblHYiUz0jjoZRXDoDt5SP+3YnH664L44TimjU4ximQdZB8gAFi7vRkvLYss8sSFOa/YVABJTsE0H8R3M85YkSgpJl5aTpdzZ3J48VIwe8UWnP7LV/D72R+UrcxyoNcTh9/PXg3GgjDTgCRWqhCd0HqI2nQO3cBayWrKGkor9hYp0y7Fz8GGNJ9TKVbqjPMclE5wybTXJf+QtGitsFhpxn+/gi//LvLlicXkErlPIQpvGie4eCQDfT7xm6m4hGItel9r3Iz731xd3MMSeGiahaHfRndBrycOsqy8M6RKukXGxlbmK7hw6vCGtAiJEyyuc0haK/EzMopBnD9xe+FysuW6GFuVxt2K4yrlDe/OlnbMVThrpRm6xViDpcEu6ejYOOegFk3KMH1P13ka0zkoCEFasRKvtnHjbnz/L4uMeXs6PHFAPOxvmpg+xUK3W7GJ66NJZc64fW8brn10AcZf8ySeXpS07Rbx/qbdRhHDF++cHfu9anMUijuucwggOifJC4QLou+g0Xybni3xvogqhVypM2jymyu2KlLjNV/14DxFniRM62dniyZFgiCemaKK8qv7HQcJ+fS5YjqHMiikexM8cZAiO3bGWNFyDo7WSrY2/vDxxfjjnDUAgL+8oz+kfU9rB0676WX828PzzQUK+PuSDwvX4mlv5eK4VOSvK6avmnOoXEu27WmLEVYR8vd+bmnSux5I6hhMre3sNVHcEJ3z60hHwYSLhJ+DQeQT94/Qv4x4L5dnqfQa6norJ1robqdCWokDEY0joheJaCkRLSaib4bpg4noWSJaHv4fJDxzLRE1EtEyIjpLSD+GiBaG926hsKeJqJaIHgzTZxPR+PK/qub9wv+RPD/6QJ2tc7DJPwvex5YBvac14gRM78A5hlnCAes2LFq7E39bGHAj//L7uYV01WJaTP+VFmHTfD/NQlCl0DmYFtQ1WzXHyzriJ08uxef/Vx3K3LXdafq7s31lVLv2lvZcjIuQc7jqHEyvovKyjrWrG3AOnWIgWQRcOIcOAP/GGDsMwPEAriCiyQCuAfA8Y2wSgOfD3wjvzQQwBcAMAL8lIu6meCuAywBMCv9mhOmXAtjGGJsI4GYAN5bh3VKhkopQXV0ynK2VLOULB7NVZOA93LAmkVb2elKI96LdnKX/rIHqopeoFTxrlY2ScP6tbxjLtuG9Dbu091zHpPwJzGKlTiYOir4/59ev4bSbXg7ao2iTaT5s3BVxrWl23Jt2xUPvp+ccUmXv0bASB8bYesbY3PB6F4ClAMYAOBfAvWG2ewGcF16fC+ABxlgrY2wlgEYAxxLRKAD1jLFZLBgF90nP8LIeAXAaUed8hk6qJoZi2UdXAlaViT6raziQNJD7TFdDMbtT/kTcz6E8C1maUnhMntjzhgK2WJTvtsXY5FYh9oUu9DSQbiyn6YtyEBLVWHhvQ9y7OsE5ONZrytYhseMykeoGjEO3RSqdQyjuOQrAbAAjGGPrgYCAABgeZhsDQNxaNoVpY8JrOT32DGOsA8AOAEMU9V9GRA1E1LBpU+khmFWID8jKEI6S/Rws5Wdc3bwrjGJY9lLWIduzacpW6hwM+fvXmc/NWm4LM+Ho5Sv7MhSLNIrYcjAZLmMhoXPQPHLT35c55QOAlZviZ5nL7Ug7RisxnWSjmO4CZ+JARPsB+BOAbzHGdpqyKtKYId30TDyBsdsZY9MZY9OHDRtma7IT5DnZKQppzSiw1c3ltrYdVbbCtMG1zOLk2kkC6FqM1VrJUpD4Xpt2tybumx639Ukp98V6d7fqLcsSYiUjN2JpUJnhNhbieXSP/OqFRukpfdmys5oc+bhUnUNZxHP8JMPSSyornIgDEVUjIAy/Z4w9GiZvCEVFCP9zE4omAOOEx8cCWBemj1Wkx54hoioAAwCoTTfKDFnsUi7qvWNvO8Zf8yT+T+Eoo6vCGkPIsW1ZUazUhULSUsJCF2HJakWa5qi+hWkRsvWz7TM4i5UsZ364Ik581e/1r6cfnMibuh5uYWdpNmMsVVTWeD739sjnOpjqUHJX0ofsykOTKg0XayUCcBeApYyxXwq3HgdwSXh9CYDHhPSZoQXSBASK5zmh6GkXER0flnmx9Awv6wIAL7AKasw+f9usQvgBedKWy1ppbXjylJI4FMk5RKasFs7BkR80lWISO6j6RRXlsxSxkvgdyiF7lsu0QdX2kjgHK/HQ3xe/t+nMjzTj1aVPXc7mttcT/HdxNtu4Sz6n3ZVldP8w8oFa4mcePaAOMz8S7Wtd2txRxjAa2/e2xfyIuhouy8hJAC4CcCoRzQv/zgZwA4AziGg5gDPC32CMLQbwEIAlAJ4GcAVjjK8clwO4E4GS+n0AT4XpdwEYQkSNAK5CaPlUCTDGMGfVVvz3c8s199XXurK09wwLke4x22TgAzHNXNWtF4vX7cBDodWRvBC++O5GHPjdv2HpepP0MI5rHl2YSCsmNAFvSZp33B062/1O4V0ca0+KMs+aMjKRZlyDLAszATjs+09rzVVNj4vtNp8W6E4dRIKjI0zlsHbiuWwbBQbgorviDpeu3ytNgD65/8R2VWUzse+oarNcnqzwLgXPLd2IU37xUtnKKxVmLRoAxthr0I+T0zTPXA/gekV6A4CpivQWABfa2lIOyKF25RezjcerzzoEP38mUIjlmTpAW3NbDss+DEwTVYpNfchuM/hgtc1LcXOpW7Q+dUt0SMsO6fChJ0Kl56K1pcV6KWbiqN7N9r58N/3AW2tww/lHGAoHzj58JP628EN9nhBTxwxQPF78LnHWii1obs9pHd3MAR8FziGFWMkYcyi2CbIv3DJcF+6gbHISEcllunOMxXOEYruqs5Q6QkI5OIfuah3b6zyk5fUqIVayjIiDR/SPypJMDPnAu+IPc3HVQ/OV5ZvqsFmQcPmmbZFSHcKTBnwnrrfAiZe6QsMKl+JoVeyj4695Em+vVi/AOcbws88aiEdJbTL39Pf+bInD46g8NpuymqsQMc8hzHYkVlK1KZ0+oJixUBbOQeoTWVQkEovqbCbWTpVYKWnAoq58b1sHTrrhBbzxfvoAibc8r5ZqdDZ6HXGQF9ZEyAEm3ks+n80E3AMQHxiTvvcULrgtcIR6c0XkcZzGJNJqrRRSNrtsPUIxCuk9bQFx6FdrZSyNSHH2TAGFKJ2in0PKdeXC29Sim1yeOVNLpfmc1I7/e3M1TrvppVRtS1Mfh7ga6jplAAAgAElEQVRpkBWqrmXIeGpRxD1pxUqOehAT+Pfctkd9NK4J7nXoIRucmDiHKolzUImV5CSdXqJx426s3d6Mn/3tXUPr1FAFX+wKlDb7eyDkb5kUK5kHJBEVFny5rHfCSJliusoKRXsSlqXugs7BmEs6o9eSV4WW9mABqqtWeQm7oxSx0rodLXh79VbsbcvhgMH9UpURDyteHAui5Piknv8PISpnqUZhJrGSWGu5FNKx8q1ipXTKeVU+bqBhy1cM0kRvlblzcYhWZTLxQ4kUxGG3FExSVzX/nssMnu869KvpHstyL+QcLPctGbJEhQVft8jnbeyHViFtrrtgNmfJ9/h8wVGqhEXLbbnSw5U2iKISseTzb52Fi+6aU5KsP7EZMPRH/J6C43NUfMo6HBnrtjfjtJtewvod0YJp9kmIKu6KKKKliJU4ijlcyJlzSKGQToiVEjoHs1hp9RazU12h3rDito485n6wTd9AJDm0frWlbcrKhd5HHCysg3hXNeiymYhzcJmnant5NeTJ0LRtL4647hms2BR41/KBaJo0G3a22BvliDTiL9Upa66T+++LN5jbkXI9FLk1+VFXWpmCpidgipMEAH+c8wHe37QHDzdEAQNclcem40vLHSqlHE503//LIlz3+GLc+FR68Yq7zkGfUebITGKlddtbYt9YRYhlvwYdcRAPi9qy2xxaRe7mzjiN0gW9jzhIv+UJJQ4WlViEKJo0PK98sHtsgCkGrusB60/MX4+dLR14MDQ55e0xzZnjfvp8vL3KHXBpOzJV+1WKUtXju1qSu2pxIpXDvUXcicnlleIUKJZ156srpHKj6z5FiOOMYiWmvgaA86aNVrYBqExcLQ5Xwv/w2024541V2NVqPtujWcFZlEPnIMfIShAH4ffa7c1WzkFOsomVAGBAn2pDC5PoDpFigd5IHKxypehS9ZGyos4hXA9lJWHcWshdPPHye5p4USwYxNFJcO6DR7XmmOTWsWpTOOupDvZRTe4PFKGt6xRB7mLtMN5NwuRt7ApVEfy9m7btxU+eXCrlj57oU2MmDmnpX7wfGfoJ5Y8f6qaPefRfTtTe0x3paor5w8pn3q+Faz+ZiIh8y2StRBQf20pHSGk0ujjKpd3wdHY4dR16H3GQPq78IcRfKhvmTCapc5AHiO3b6j7+Xa+tjP0WF3aRnU0zdFSL3FbH4zu1/hhltj8Xd9pqP4d0kyV2aH3inuk54VpBVflicd5vXk/c+7BEcZ6rZRBjcU4r72h8MHlUPYbulxT9AcC1CgfGoE368jpjASuHZ7x8SxYViRxLv5oqSb9jr0u3y4+ZxKbkBLoJbeh9xEEeLTIBEAekSqyUISrIfQvhLKSPL/4qRXYt5hcHWJrBQxSIfG546t2ColSOaa+v1z7wjc8rsqk4qZoq8zBctSVdSIH4KWH6e8YyFGlcfLjZIkO2LWoqEZyJ29kpcGUyceiI7XzjhYg/66qzeOLrJxvbJeJXXzgqqlMaB20dedz52gr5kbIjnaOdrgyZU4jfl40HbCJhuS5d3eJ8tcVf2rKntDMmKoVeRxzkbs/lZZFQdK36qNkMFSahy9Giyjnv+O3FZ9vzaoseexmEvy5Yh9tefh83Pv0unl+6Af/19DL7g4aKXK2QVP2iNIu03L/i9++4VRhClPeW8+jFds2HTi4Y5nJU8nUTzfrGH6P3Z2CxwIpplJfZhE+PvqGfOXJ0oU1ytttfeR+/efF953qLxX2zVjvlM3IO0j154d0pnWdt1TlIv3ViJXGOmDiHXzyzDD+VfCFcRFWdge5hUNuJSOwkDGIl1UfNEFAd7tw41yHm+8Fji6ymrGl3Brta2rFHVOileL66itDeEeRv68jj0nsbnJ/V1eLOOaTPV455YeQcHJW0qjVX3khwrJQ8xG07XpVZp6uiPM/igRVF6yVr5NeMTBzSmPZG2KnQL1UCmxVh01Uw9XeSc4j/5j49PK+YXWWt5CpWEtcVpbEGY9jZ3IFfv9iYuFfGcE0lofdxDomPG/8tyv3lWPBAsFPjzmF8kosD8L5Zqx0C9rm1lU/OP85Zg9PD4xSBdJxDdTZT2D277jH5Ym2yVlru4NxTTPgDpSgqpYLZtJt2FyslM6rGAwCcKnwbwE48xQWJQ1y3p40biCPGJmM7AaFYSXgJmRsQId+Rz8V2HUdyPpGYf/vMgx1LiVAOgwEAuOPi6UF7UgS5TG4GBfFPjlnDZyT7Ql2vSDRUps33v7kaR/7o75o2dw/OofcRB+H65BtfwB/nfKDNq9opZohQG8rIuZWS2c46mabLb7La2SNYlKThPKod43e7WGYU0hlwxs2vOJRq33kBgVMYRznkrcWuPeLuXbXmugZZK0asJLZ68uh6XP7xg9Rlg8Xamea88CTnYG5oZK2kz3fAkH74/qcnmysW8PGDh+FwRVDDYnBkSEBNm5DEudSyflD4OWlEf7u1koUTUeX7xd/fwzuSI9zrjfqYS91FrNT7iIPQ8U3bzC79J09MnjaXzSQ5B5NMURmjR5P3qHGDpGfVsz3N2KlWhY1V4Et3vulcT7kPYfnmA/OsdaaBuEi+8O5GfcaUcA0HYuccksQh5rjHDGImSRRkcoqTkdA5WPK7cFnVWcIQhQOkqcxyLX1EFLTRqJCOrkcPqEvoERmAmmwGj11xEu79p49gzMA+0bMqayXp969eUAfJk9eEz/72jdjvQX31fdZNaEPv0zmk6XjVIfPZTJTOxQNpZYSlfvzF63aiI5dHlQNXkKUomJhpsr+5IhnJVK9zcGgk1IukzaxP9YzODl8HccF8uGFN7F4pzqc6sVJaqIhDreQ4p1vzGeJiM3HBt+lTsgqdgwtM2aoyGRClCCOeol4bspngjY2cg9D6TIZinMOOve249aVAsX7kuIEAgGs+eSh2trTj0blr1Tt4Kem5perNh233P6Cv3jGumzhI9z7OIQ3U8m+Bc+hI6hxU+ZPlaljRRMRYfdtUC5UqhouYK63HrLadJZiy2sz6yrFuiG8pl1eMQnr8kL4A9AppGTbOQfXt+sQ2IkyrN8kzBiLgjWtOxVvfOx2jhZ2uDQni4Njb8uuIv6stZsgyiKhsFmQZCgilqTxxuGUofq7ERXfPTuSvq87inCMDr3N1VFZHrtkyVEycQ3dBryMOaXYtqkGXJSrs1vguxLRLSCNWKrVtn5NYVyC0wHAvNiw7gE7GXgrnYJtc5dA5bNvbrlWqOyukiTAo3N3d/A/TALhzDnaDhKSBgGieCgAZzcxkLHhu9MA+GNa/Fp8+YlT0jFBg48ZdmN8UP6wpacpqbqeLBVW1YNrtgnJyDpkMJbyaE4gRh/jcW9CkPswqip3mpjNTwSaCLDUcfmeg9xGHNEslQ0H5zJFRhOxOe96yboDJyXKIBhHuDkLmtpjw1fvUZq+lHKquEyutCcNqlCusTKkLEAHYf0gQmoLv8lw9Xe0EUJUaJQYEQKNvgqw4V1+f/sukwYBKP7Fxl4Nnt+F1qrKZVPwoURnFSkQgyylz4r2Ac3AoNxPf/IlwXT+s5zh0F8WCAVbiQER3E9FGIlokpF1HRGulM6X5vWuJqJGIlhHRWUL6MUS0MLx3C4UjmYhqiejBMH02EY0v7yvGYfomYyQWPc+S7H1NVXTObMFDOuV31g7mFOW4L1TR9QNvrdFnTIFJw/dzrNtd5/DOGn4WRvpJoyJ6vO6kWMm9nLsvmY7bLzqmEDgt7bnJMsYNDsbX4nXJHatYtMn/gDFWtnB6jAHHXv+89n7BCc4wMMX54AYqm0I6Q4SspEcAAr+Cxo2B+WieMXzt4wfi3R/PCLkMe+18zrsE3pPNg4FAl2E7x8Hsm2FtYqfAhXO4B8AMRfrNjLFp4d/fAICIJgOYCWBK+MxviYhr2m4FcBmASeEfL/NSANsYYxMB3AzgxiLfxQmmfhdlss8v3YA7Xl2ZMDscUV8bHZ8Ypq3bobd6emtVUg+wSLE4BOXZRwWPQ5TGwaxcm5QTDxoCAM5yblW9OuLAu76Ytp43bUwiTQ5SOH5IX8z9/hnOIhAiYMh+tThzyshUIdrFOmVwIqP0HBfS8oyhSpIrMYHY6V4hLdGwH2xlL6NfTTaVLivgHMqkc8gEp7fJ4r7rHl+M03/5CjbsDEJw12YzqKvOBvoJh7ojziF5T366I8+wZN3OWNoNT6u5htaOaC0xtaPH+Dkwxl4BoD6UN4lzATzAGGtljK0E0AjgWCIaBaCeMTaLBW9+H4DzhGfuDa8fAXAapRFipoTLqVGMMfxxjnqXTSErK5b1T797y1inHNJ78dqdynwuY+JjBw8FkE7EUc4QEhnH3VdQdzJNp5/hfVquuDJyOYP61WBwvxr38xyEnJRRl6kDzyZ/d/Npb1HZHXmGEw4agjMnj0iUyRhLJeN3aWcp+frWVqXiHExZbTG2ZGSIUJPNJOT7/Jjenc3tQdvDBhK5WRZyI0AXzgEALrs/Ln7dukft2f20cDTrvsI56HAlES0IxU7cQH8MAHFVbQrTxoTXcnrsGcZYB4AdAIaU0C4jXPpdZuuX/ijOOPFNnbtyKsn2fmT8IPzo3Cmp21ZTFXAOaXQOL5bJ1j9DgfLR1RlMKVbSPBuJ6tK3S3mgkqyITllmzJdA0jHZwN9BPEZULEcFseyW9hyyGcKPz5sa3RfyaUtJ+ZK214nESnoEnIM7TDqH4f1r8a3TJzmXlSUKOIeOeIFy8Zwrla2VdChwisrBmEzbT1Iuq8LXA/H5YA4z3j2oQ7HE4VYABwGYBmA9gJvCdJ1xjslox9mgh4guI6IGImrYtElz9oEFLv0uN1iOzx/tct3qlHf5jAX24ZNH1Uvp9gK5U5trdMg809tipwVRMNFsJp0F8VwqziF8pMycQxprJdJc246F1dX9hHhcK6D1XQDiXcX9IER5tmh9VS6+2uoh7VBRXXU2JeegNz3NCFy5U1kURABol8cji/0rlBmY0drBx29y3jLlXKqvi/ssiMShf11EOJw97J1yVR5FEQfG2AbGWI4xlgdwB4Bjw1tNAMYJWccCWBemj1Wkx54hoioAA6ARYzHGbmeMTWeMTR82LOm97Nh67R1xgaLCbiOZL1r73D6jPCjyjCGTcQ+2JqIm5HlNh5bIdZUTBLJaKxVktipTVu2Zu+l25yJUj+jeu5g+l3UOx00YbG5PmK+lI4VYSWhuc4E4RNMzOj6caRfQ9DoHx3zyRkS4DsZjWp2D/n6auEtEhOpsJhl2P/z/w8cWx8rMOOo7dArpvy5Yr8wvO8vuFE47HFlfV7gW52iaA4q6CkURh1CHwPFZAJx/fhzAzNACaQICxfMcxth6ALuI6PhQn3AxgMeEZy4Jry8A8AKrIF9lKrmwQAHYGJ55oDryMa0IRJaJcisoeSK4iZUyyro7LR6Lg86B73hVuXSERbYAKxXiYhqUn3LpFLLLbbO1UOXHAJhDXYgbjfGhCW2VEPqE3zdxDi7v+MrVnxDaac4rG14AQMOqrTHZOfc1cIUpfAYDSxUOBAgi4j4+f13sjBLe/7NC3QMvs1RT1u3NySNuxfwcrUJgRXFDIM7RNPGgugpWTwwi+iOAUwAMJaImAD8EcAoRTUPwnVcB+BoAMMYWE9FDAJYA6ABwBWOMq+gvR2D51AfAU+EfANwF4H4iakTAMcwsx4vp4NLt981ajXc+CEwruUjpr18/uRDELtrlun1EeVefYzwuTHxQuRTHOQd54OokPcUstrpHOEGzKcMLToIqzkFTOJ9Ef5m3TnnfBFWRshNcap2D8ERG/t6WLi2IM6RKTRFUwYCxg/rgv84/AkcfEKjwYsQhLDTP9ERAZVYpY//Q2zvWUA1U3NwFt81K5rPWKuY1WwwVKzJbuHY7Tj00UODrhmdaU1a5HLF7X7/mVJx0wwtBudLzIrEQiV1P4xysxIEx9gVF8l2G/NcDuF6R3gBgqiK9BcCFtnaUC0bOIfz/hhAxsTZUAE8VIknKCsoJQ/slYvqLkA+JYYwhS0kWOh3nIBOcyouViBzFStnkosKhIyy8K5auV1typYVcTWrGQaGQlrkRfd3BfflddV7PQZlBnSdOHFpIq1Y+oPdzqFNwuSaoo8NGcNUDpeLKDJwDgcwE1IA9rYKZqCYMDTlyDlyv15aL94+4YRB9omRxofjZxPBnovjL1KXdhDZ4D2kVxG/dV3FYPL8dOcGZy5QtdLhYKSE7dlJIa8RKWp2DtchUSMM5qBaVxzScgWnhtMGkc0hFG2NB7CLICmlbmeIuX4RZ55DUJYi7TrFuvVhJnc79K2T8+K9LtO0BhE2QMVdazkFfIANLBJPU6aj+dPmJsd8vLosUxfL34Rs8V51DfdhfO5vjVkc6xkwWhenO24ifBW7iHLoHeeh9xMHU75TMI1sqAUnOwRpoTdY55LlYSWqbcL1ordpRjnMO8gKtm0Smpp03bbT+pgKEoN06zoHb5Wc1OgfGGF5+T21lRiA0O0RflU/Vam7LKSeT/E3SBh2MPZvgHMywic5U4JyD9j6z59PRbJ24ab3BeRNCPbbxnU7nYLYYkkPMP7FAvZk45oB4ePu4ziGel+sNAyc4exs5MZXPl9a9p9y9bmIlff26e4/NW4tte8znl5cTnjgIUH37uiqTQjrkHCyONYmFnDHlDlxs2yvL1Ysoj/WUOHRE82KmXcjAlJEhAwfAJCHiLDafFNUavchDDfrwHUTAq5p3FrF6y97C9RuNm3HYD57G7JUK4zZHEZC+PdIuXth12nZ2urtGU1Zm3oGzQr5kSBduZ69rV1olb+E5xWZJBdWiKcckA4J2Esz9J3uG63wGZAzvX6e916cmKNPVgbM6m0FddQa7W+N168Rn8vcQiYPIOXRYdA48LI3q3pqte/HNB+bhyj/Otba/XOh9xCHlYlGnEitJVhy2xULe7TIWDCjxnIJTDx0e95LV2ETzXaBMDHScg2kyDNYc0qJ7ol/oDSsTNa5A5ZMkMmWNP7/sw93atrhC5OS4NYq4a+SIwmeECcUbKwGIO1Dp+ufz08eGdRbLOSTvjxpQFyszryAiL377FIwf0lf7rdPK8c8IOUBXr3UVV6Y6gXDWtacarZXWbG2OKeEBwfLNSpCj+2u3xzkivsEjRyc4QK041/Wi/F1FMV7MWkmYEKpmHL1/MI9E3wgOHnpj/Q6HQIllQu8jDg5jQ8zys88dnrgvWisxxrDO8sF0fg5726KdCUltk0MvcGQLu/J4ul4hbWxaKowd1EcpVuLDn+9Qdaasrk5gAPDZo5LxkuQyTN9SXshTWytJD4hmkLp6o7Aq6vtGU1amVjRfevIEAJJIS2rcsP61mDi8v5aDlU0tC3Uq0l79zifw6y8eBbEa65xRFK+qsqYqYw3ZXSMRFa6DaNPMBxfwDR7BfT6o/DF0xF1OrhWkDSIjFCcOyYb86LwpGD2gDiPqk1xQsVZ3paDXEQcTVDs3OVIrEGe3H5271lruHok9DQ5sIRwysr9Qd/wZ2cKJQ2cmqj/L1tq8BHQDsLYqgwwlj97kbeebvoxmt2daHGXOR5ezlLMkSoFoBqkr2eb/YtrBB4u+qsw4m6ojIiYzzTTK/hH1dYXFzdUxUdUe1SmFBLIe9qPjHNo6LMTB0Ma6qkjn4Co4IEWRun6UiUaeMRw6sj/m/eCMGGGOh89IllNblcWUMQOMY7yCYecS6HXEwUXn4Hrwep4Bm3erg2yJeFmSpedDsdLYQX3Rv7YKowuigyiPbrHn1lPyQqDbNRptym0Nl1BlOdglwTlIVW81KNMSr6tTujpSB7lulzklZpGDwMWUmZo+jXbabos0k36omqjyxle9i4krSyNWEpXXrpEAVGNCxa1kKMkhJ+uPdxIvp9VCHEwt5HMmk3HfNJBCea3jHBIm6SwQDQ3sWxN7RuS4dX2q08mUd6vjhu5/HFGZYT7S060McdK4PPO/L6/AtZ88LNYGvkF65wdBGOmv3f92bACIxOG5qz6O0QPrsHLzHmzc2RqWEa+jnH4OWoVqRm3zUxArJXQO8ZIeebsJOqjOzlC2LSXnEHkrpyOFfWviUyNDEWGykVsd/ZLfT1ws5q3ZnrCOiUqM61BU/WQKKnfMAeZwH7FyhJVO5oS271UTdyXnoCAOFDjKGL+hbK1U8DmwEQdDoQXikErn4K6flDniYF2IzwUgPqd/8+L7sWdOmjik0EZVE71YqRNg+tzOC4gwaUwLmu5QHHEhrMpmgoPSST/AJw7fD31rqjBl9IDChE14XUu/jxgbOO2lIQ02jolHZZUhp5nCZ8jQKRz1R5Tad19BPofKU0DWORwqiAQ5xHe55/WVyjJEiDv6tdubE9YxQFy/tW1PG15r3Kw0JTaFhvjp5xK+pwBQiAKgg+wZftn9byvzqaaAnnMwzzFZHJUhwnsbdlnNbk2fmzsHyk5wxhDhKXQOSbFScsMEmAPv3XHx9CC/hruJwsDom1xu9D7i4LBzsMn1Ch9cKkteMMYP7Ycjxg7A/oP7xtK5n0OsTmurAmQ1i6k8oG48/wgcOLRfqkWytSMPxhjmr4kWjc8JiuFs6Jshz/todxtUltFwDiroLJtkCy+OtJFRuSWTk1jJkCeuc2DKA49E/47rnkg6mMn91q/W7tEsqhz+NDfgvMTvo2qfjFqFOXYa8G+zXHO6mWrBV/Wl65nUIqqyhDNvfgXn35oM2yGi4HOkGPB9BIW0OG8OHNpPW57rwUBBXrkx0RohEgfT2OWcqs6iKuIcvM6hYjByDukYhwTn8PS3PhaviwWmiLKXNQv9HOS6VWNHZs91cV9kzoHIvGB8+cTxifdt7cgbw4BkMhTaiks3wnL4zoibMbrMrSoNIdGJEdKexrYq9IvQfVvnBSATLRY6nwQuAtErhiXxg0PVolJYZbJbKNvRwSsNZM6hpV39TVR9q7XoMnDIx44fnOAcsgotsIorMUUriMRK8flv1D9Scq1wdW4MDE6Ca3H+ygE4lfVq2lUgDp5zqBzU9sUDASBxvoIO4qQxfay2XD60l47SNu9uxbodLQkdgS7OvTwRdGKlBU3x3eTetlyoVIvnO+fIyCv6+APjZyq1tueMu5usQixwxuQRhTS+29c56qkwZlCwA5ft93Wmi7Yy//OcKYXyXOC6oMpiJbUoxUwU5TAWTqFcCu00RyxNc0KfK+TwGfKY432tapXe3Fe/QTt87ICEtZL8vf/fRyfgpW+fknjWdJ67aK0k9pEpkrFqkZb7l68XchQFBoFzSLnCZkjtQV7O0xyd29LpNXY5kp382aPH4p3vn4EjxwVEwh5oLPhv0zm0tOfCHUhU3o9CccNLy+IWTDLnwK91cmq5jf/+p4Wx35w7kV9FdLCZPj6uqGxpzydY/8PHRgEHuW5ExD9MH1dI48ShRhP/ScZxEwbjFxceqXwfG+ewq6Udb65IekZzb2HX8Bk8Vz7P8L0/L1LmAeKLbyj9TeTRcUEcl3/8oNhvFwImipVMUVczRM5Hx+pw6qHDlXUXAglK73XJiePDjMmydHNIxyEDwfiS/RzkvGdMHolxkpgWiLhW1SLKiSoRxaz6zMYpyc2avPG/558+AiD5XUTOQRx3LrRbR+QjzsGLlSoGnehmUL8ad7GSsKMyPdPSnkuw+9zTMeF9KZXDB0g/6QhCPtBt5zccvf8gpXWGkZh15GLzfPyQvvjyieMLnFUgVtI/364RK+1qacfbq7cl8p931JjCblpuZ4smYijPd+Uf3lGWqTsc3vZtbU5WFOMcgsn/D9PHxfKYDjkCkl7DTmIlwbFO58zG6y6VONz95Y/EfsvhM/QRdVViHnUdnENWEX+ipJ+D3Jc6JbJLUER54TWKlRT35Tk3vD4QGYvJ2/e24Z0PtkfjV3gdF85Op3Pgfd+JUqVeSBwUaXzSuXa8aM9uouTD+9cC0oDkE0m1yIpt46dLPfS142N5ZJ3DorU7NBNNrZuQF5jZ3z2tsHtvbc/H2pUJ/RomhlZXKnv50DoRgMA5SGHFD7/u7zj/1jeS7ytwXvk88NX7ooPaN+9Wm03yMt/9UB3am7/fL599T3k/2QZHnYMkKyckJ7uOMHEMksKVuJ0tELYTzMg5VGUzBeJcPkTErmnbXn0uRbP0J/EF333uBwrCTpTwc5CLUcVsAiKzYJupunjbRExVOgfVWJHFZN94YB4A4K1VwfvFFNLCuOBhUWxtLLQ1TPQ6hwrCqPR17HkxKqvuiVu/dDRuunCaVsaatFaK6wdyOYZ+NVkcOCxuDiuGj16zdS8+/avX8J9PLNa2U7bekInSiPo6jKivBRByDsJtTgz4pjqTIaUMVRYr8d0fg3vkT3lX9uFOdUgS23rK3++5pRvMGXl5KcqVw2brCK+Kq7vwmLGJNHmxOf2w4Yk8MRGmgThUZymm8DxCEAcWC5Fz+Mo9bxnymTc6IkxnKhw4bL+EoYYr58AXejPnEBcVmYlzUsGvYi5FZ7lNu1rxihR1WHwfcVzsV6t2MavOZJSbPT6PPXGoIFTUPzXnEP43OW598vBRGNC3OiFW4tlrJPYZ0k6lI880tuJ8d8oKTlNzNfbqmUxSXKKIaqA9UD3iUoL0LCXPoCACtu8N2sHtyaszmcJO23aoeuKUtRAjFfFlVG2UITtRcejWAZ5u28XHFNIIzl6QnylYKynaKJ8zLLeprjqT2AgAoljJwjlkorOU1+9oxvINpQc5FC2lVA56Ub5kmj74YHBPdfv8o5PxtOSu1HEOLt9R9gWxWSvJJE5F9EWLQKWfinAdD5+hrry+T5UyEi1vtzdlrSCMewXHfhetOH7wmF6JycsUJ0p9XSBjLyjzeD6pcbl88uATIO4XwHfosk8A35VkiBKLs2r3KRIBccz2De3w+YKsUkgDKCjyDxkR+HlUV1FBduq605cXgYf/+QQ8eNnx+MKx+8fSbYv4foLC/XXhRL8Lpyd37mnKjfk5sPhvDm6tpFxELHJ5rXGDsHs36bhjhbEAACAASURBVHuqBM7hhJ+9gObQGOLV73zC+F4miCKttOoMbXypDCHHmLK/VSJamcjInMODlx1faKOp3qD86Ju1tOeMeiaVzqFdJb4VrlXE+8pTJxauxU2Dbrj1r6tGc3suMadL1ScVAytxIKK7iWgjES0S0gYT0bNEtDz8P0i4dy0RNRLRMiI6S0g/hogWhvduoXAkEFEtET0Yps8movHlfcU4VB+FDxhXqixaccjeqnO+dxpevvqUKC+SuoT6uip8TbJckQ9B0XEOoikrH4yNG6Nd4sUnHIBn/vVjhbrlQaYLvQAEfSOGEf/tl44GAFzzyUMxY8pInDllROJ5AuGfP34Q3rjmVHzj9En4+MHD8PVTJxUspWwmeKKYTMTogX1w3IFDUN8nzn7b5kj/2shc9Et3zi5cnztNHeWVt8829UQOkCH4DvKEHd6/1qmNct2A2veF18vB+4hHahVRlUluBBiD0rLHFaJIy6a8lWEKW57PJ+cNByfuh4VmonIu2aHvuAOH4CPjB0XyfJtYKbx/6PefNvqNqGT/KsdMlbm4iInDI8fY+GZA/Qy3JtzdkgzWydvVWXDhHO4BMENKuwbA84yxSQCeD3+DiCYDmAlgSvjMb4mIf81bAVwGYFL4x8u8FMA2xthEADcDuLHYl3GBarHi39zdWiksS/F9h/evwwFDIs9L2WKopT2XsEACkt6buXxeuRMRPaRVC/0JBw4pRJIlokR0V5VSWVygxcNERg0Iyhk9sA9uu+gY9K2pUi4E2Qxh9MA+qK+rxr1fORYj6utCsYsiswAxBo02r5RuUyDzg11cUVjwLf5JMVNWxhL+K0B0OplKrKRW2savVd9TFGHy/N8+85BEvqpsBh155qxgd4Ecml6fL5lmii+VZ+o+AoL588SVJ+O/zj8iLCeeTyVWEs1ObQrpVOc5SINPZ2HFc1mPC46JldR5OGckczXd0lqJMfYKANmg/FwA94bX9wI4T0h/gDHWyhhbCaARwLFENApAPWNsFgtG2X3SM7ysRwCcxrmKikDFOaTseFEO7JA5toh8uLNVeRC8bB1h1TkwnadodD1vzfZEqAWVWEncIa7YpPeQFvPaELyPi1gp+K/ry3rJccxEcEbU1yr7VsZr/54UtaTTOQAg4OunTcSEofGNAKAWAZxyyDAAwPWfnYqLTzgg6EdJBq3iHMSNiGn3yJ/9w5wPjO8BBN7xLogIk42zUo1D9RPZMHaQ2EcvX30K/vwv0ZnQh48dUHAsk/tS9oMAOOE21xvkc/ciV3EOKjGU+Bl13BCHSBB17eTvJxMiLjKc36Q+PrgSKFbnMIIxth4Awv/czGIMAPEsyKYwbUx4LafHnmGMdQDYASDuultGqD5JWjMx2f7bBHGH+VDDGrzy3ias3pJcgGUZpyg2UtWd08jzbWIctTmqO7FL0G1Nn/HJZSsxMv+Mcp59+MjC9cUnHBDLzyfVhp1JkcCzV31cuXjI/Th2UFLUYmtnbNfJgtc+dGQ9XhS8dUVxo4j3fvJJnHpocLral447AD86d2poChsWF35L1Z6IE5yOfF7rGAkAH2wNTE1NjnwcR4V+KzZE9ah1BBxyeJhfPLMMu1o6cNCwfnjyGycnysxJYqUDhvTDUfvHz4TmVbvozERLP9N3TMc5xMvavretIHL9w1ePi9ojcC02vYDOx+I3Xzy6cM05B1mEZT3PogIot0JaKX40pJueSRZOdBkRNRBRw6ZN9vOGVVCNjYhlc6MOhR2iUNi9XzlWkzdadF8OzdxUY0heGGycA9Mo9cSyv3DsuMR9k1zbZd6Yzg2QyzSJI6aOqcfZU0cpxUr/evrBhev+dTLnoG9kfV21cvE4dFQygipHGmslxhgWNG3His17lAu5qA8SoTK/FBcV06K/cVdg0vvSsk0FAqr6BmtC4jAs1HuY4MqYixyl6RwNrh8YEvpx/PrFRgDA544eiymj4ya1QYwqt/4G7M6eQZnhRoQxrfMkL1NVnOpAL9FElTGGaT96Fne9thL9a6tw4sShQr5o7NqIg7jei+8vboa4o6Tos/L0ovX431dWGMuuBIolDhtCURHC//xosCYA4oo0FsC6MH2sIj32DBFVARiApBgLAMAYu50xNp0xNn3YsGFFNVy1sy4sYOFksJ0rkAmtdsQJM2agwakFwI1Pv4vnLbb3YttyOZZwCBLblsureQQxAqzKlpovVCLh0SmFVUiYsmryFWTLmiL/+vWPYlC/GmXdk0bYF/NkfcF/FWdkMqd1kVUH5QehF8759esA1O/N+2ajQdFZyKsQhagW/RMPChai+j7Vhfyq8XldGOfoSAf/BlfxqbhpkHtHrue8aaMTurSDFKa5GQoWfJsIJqMhtLp25hnDXa+txMk3vmgsU/7OXzlpAl6/5lRlfj42xIOGkoQ+IiK2toobpZhjrPA9I+IQ1fnP/zfXGl69EiiWODwO4JLw+hIAjwnpM0MLpAkIFM9zQtHTLiI6PtQnXCw9w8u6AMALrJxaNQmqkg8Ndz78E5mcjTiqMvJZyupnuD38rS+9X4hqOahvtSJfvG06zkE0ZVV107RxkchA9R5fOG5/fPnE8fjm6ZOifIrduy6csbPOAVyJ6rhDdDTxyTOmDB0dBTqLGsjb6nT+sKV6WSSh5sCC/3e9ttJanbgzLSz6ikIHCmPFpHOYMnoADh3Z37qx0T2vzCfUK3/GRKgNx/Ad2XAht53oxwmtjYhw5Bnw53fMR/aqfFNUfj9A2Edh1lYhGq2sEBf9Iaycg4NCunC4UQlnZpcLLqasfwQwC8AhRNRERJcCuAHAGUS0HMAZ4W8wxhYDeAjAEgBPA7iCMcb5vMsB3IlASf0+gKfC9LsADCGiRgBXIbR8qhTEb1JfV4W3vnc6PhIGoCuc3uQwe+RYNjp6wlleEXKMnaDyKN//PLcci9ftSMSZCdoY/E/jQyCib3UW150zpeBvIZfJcdnHDlSWmfDs1vQV1zm4Eoc0i8CXf5f01r3p80EIEJGg9gtj5K8yhCGXF2hTO8U8LocemcuLiLtp0RcdFHlMJ109pgN/Yu10NtmOOIdmSVwzsG88FEhVRhUTKNkYThR5jDF93cH/nMMiqYtkCgBPffOjUT7FXNQRU87xL1q7A/e/uaqQ3iZxoS4KaS5yFoNt6va/OoV0V8B6TChj7AuaW6dp8l8P4HpFegOAxJFUjLEWABfa2lEuiB8lk6GYjJYPE1OAM46qTCY2GPQLQ3LgKv0X+E4pl8fNzwVxgUYq4q+IHtLyWJTt39Vmq2Y9Bodusrkuf9z+27ab0snpddARmxPC8OPiO48f2heL1u40Lpj8lpWIZeL9o+oHV31M8Hw8HEfwvJ5TzIXf28QZ8DZOGV2PxevUsaeCMh3bGFYlf8dvnjYpMYazCU5ajTdXbAGQjCIsg3NRLmWavK4PE8Lw89AdeYd5y5Xcn/7Va7F0+cx40apJN4Y/fnBSBK7lHDQK6a5Ar/aQlicaF/G7iJUyBCfOgSi5S1Das4f5RNbT5OcQmBfGy/2cFH5AKfqwmMdyaGX7jiOGH6yiCkInWmekFStpnasyXKzk1j6OhxsC4zo3U1Z1nh98ejIeu+KkVKENxHDqJp1DlkTioDZ3ldto60vVATq68oL2ueWVRUUThiZ1Dk3bzLG2ovKC/y7jQhfJVFUmk/QdJrGSy4hMY60kwha51xOHroDwTVTevkBwqI0NVdlMbDDoFgbVmQpqziFoGjO0DxB22owlFl5ZgZ3eKkXkHNSD11UhzSesSOy+cdokzP3+GfjUEaMKaWkWAUBNbIBoERU5B5f5xVl929oS6FCi348K8u2vnDwBR44bCEoxm8QAdBFxSPYm/6acczB9U16mbbftyjmI4TNs4GExAODTR4xC/9oqHKI6Z9tlciGdziHyxrflSxJPrVgJrqbq0bh4e7XSjiYGxhgem7e2EI9MRiRWqpja1Rm9jjiIA12eJHycuCyqMhute4RcxUrhABcHrmoiiYf9yJNWfh8X8RggWitFabqJkcYXRH6f6gxhsBS22hTJVAXbUY3iO+ccjmXksnQrcXDYnaaQKhWcBIGI4Cn9HMJvynUOZs4hDOlioYouSmveRsCdc+DfOs+YUiQKmA8sipcX/HfnHOxl8oVcNHfV68zUeowF152ZzBdmXL7RHuwwz4BvhmG9VaipCtqj4xx0UWkrAavOYV+DOL91imcXY6mqDMUWHy1xoOSCpsrL3fXFvCr2X9Q5/HnuWuW96LfxFYT2pNA5JBTS2lJDU1YzAeXluesc1Om8q8T26c48FtEcOjbZxUr2trkuukB8ZxpFvU3mqxKC+ZmiAPP6GQN2KqJ6inDfNPBvE+9HneKci5Xac2pLOyDJ3erPNXDnHIK+dCMijDFc/ciCqN1GziFZZl+FBz4n8hM0Fn4i5HH+xJVxJ0GbWMlG+MuJXsc5mCxOrHF+BGQzFHNUMVmQJMRKBnGROA/l2PC8PADIMeBOyWRSnpAuuhOxzFg7NZPNdfnj8t2Y6M2wsN3u6ORjCugmw2YRAwicgyWfSedgasP/XXqcIid3Bot22jxNBv+mO5rb8fDbTcb+JwSLyra96oOS5DJt4EHuZMuZqaOTvhRVglgpiCisrkOse3j/2piHuYjIUMHRWkn6NH1rsrjiEwcp84lnfWi7QqNzkCMliyavJvPc78wI4mHJY+hwyV/ERhxsDonlRC/kHAR5oyxWCv+7KLeqMuRsbqY7MSzZNjermSBvMl+Sc0grVrJzDnLT9boWSoiVyoFA7p5MV/WpC+cgikJM4E5wJsjt+sl5U3HypKHqvBCd4KI0Gfy9bn3pfXPlYRt3t3YUFkpdDCXXcVEXnkEh9+Ppk0ckyxRMuzvyagdOIH7exn51VdpYWGlMnDNhvCZx03XD+UfgnCNHS2Umv7NuA0WAk0ZaVFybRKOiYYHYHhmcOMgmswDQv7YKu1o70JbLoy5jjyFWKnod5yB2+amHqA9UdxF/ZzPktDNV7Wq01kqID7DPHZUMM22y7kmYFzoThyTHpNc5uMur84zhwx0tsbRSoVvEVdE6TaEUPjl1ZKw8m1hCJR5U5TH9FqGyclGLEY1VJurjO87rPjO54DUtw1G9g9pw4XYZ51lB59CRU0cUDtoYpZtFZMH/nMPRp9ws+N0PI+dInaWePG10bXDhFHk+VuCYgrRshvAtwcmUpwHx+S07EgJmP4fB+wX6OtO4Lid6IecQ/L/z4umFSJkcfBfsstutymRiOyrd4kKULE/NOSRjzsgsp/isqolyua6LsdJaSfM+b6+On/2rq4M7JomhJA4fYw/tcP7RyUN5Xr76FDRta8aX7pytnbAqotVq4Oxkeb9tHRAVrjokLbn0H0BUovIFXXWKHRElHC5NbeShQlxFRybUVak5B2Xd4cLLTUV1YiWblzlHKp2DQgSkC3+e4BwMusK01kp5xlBXncG7P/5kIl9hUycQOzluGBAclAVEY2JvW6Q/GhBGKN7TlsPA4o/pcEav4xz4MBo9sI9afgg3sVJG4hy0O20ojt/UWCsBcfNUlSe1qY3yhEireHRRSIv4p5PG4/gD9QF084zFdkAfnWSPhzV1TH0i7YAh/QrB0fKMlZUDiYiEOT8XXRjzpOAciKL+joiDejo6f0fhWFiTL4OLaSoQcQ6yd7QKnEvNMx5RWF3/zI9EJ/uZOIdI5+CqaJbSdPmkNL0pa/I8hzsunm4sM5dnWm6df0NRl6COIBykccXz5B88U7jHoxrsVRxHWgn0Ws5BbTEUwFXnYNqZcmQy8YPfAbVVCrdeEetW7SRFayUZCYW08JI3XXik8pChWJlCkSobdRH9a6vww8+oxRZAIIqYvWKr9ixoHUxsPqD3c0gL/q5N25qxbU+b9ZvPX7NDeUawiISBgyGv6P/Cd8daJW4K8WDEORgyuoqVQs7BZTHi9eXygSltVhhrnzt6TEG38MXj9sf1Ty7BnjYzwUmlc1CI/NThz5McscnKUB4Spx06PJGvKkNoDYlnLs+0OgyeLh6+paKfKrNyDn4qoq3vyoVexznwPlcNiipFuFwdspJC2sQ5yB9atRPk1h7iTknFzmcVC7l8j0Os5rBR9ZgxdSRUELmRMQP74MSDhhSigepg66HNu9uwdnsztu4xW87o2mJqY3tZHISiMn785BKrcYGNMADpOIcNO1vxwFuBdzbfTep2266cA1G0ETFzDklMHpXk2Ph4EjdBZyiU0WJ9/OhccWPzy89Pw08/e7iQN7jnpHMo0lpJHRgxqXN4cuF6bbli1myGlAv/uMF9C2dp5JnehJf3pWiKulexyEdi4+RX4gS2s8xZex9xKFiGKBbobJL100HmHLSsumaQytjdmsP2ve1YIJz0tHZ7MtRAgd1WiZUy8u+oHpNnqqiIb8/lY2G/ZXz37EOD+h2tkFQTwAStSXCGi76AD3e2KPNwfEwRy0aG2H2Pzl2Lc3/zunsjAfz7jEMTaa5BCWXw3b6KUwzKcWtThiJFpqsXNMevvnhUsrxMkjgcoBkbvL6OfLDBMRE0vgkzRfFIwzmodAk6zkHOJ8dKEp8Xs6pEQEAQvHNvWw572zpw36zV2vDw/HFx4zl+SNIvwmROzzm5clsA6tD7iEO4iKsmXHXGPSJiNkMxqwHV6WKAOQSGiPc3Bd6V1/wpctDRHWiuYo95m3R125yngMAKYkdzu1b2DURycdcBmtayQtdK/mryYlFTlcHs78ZjQE4anozpI8P1RDAdxNDoOriqRsqmc6DIa9/IOSheXbf4yVZ5JgMEAFjQtB3vfrhLywXxMsVnTChW56C2+Ujm0zEm4Vak8HvsoOSBQEBkwsv9dHQcJn9X/q2/d/Zh2oOZMqQWG4sOkZ2B3qtzUNxLxTlkI87hf2ZO005g1wWiVnGwuG5i6MzsEgpp4bcpbAF/7idPLg3yGrgM7r7vOkBdlJmqtujSResNABjYpxojJL2Gy1pa7OZrQJ9qjKivxQkH2U+ydeUc+G5S1+9p/FX42DXpKfarS057k75DPM9A9058/H/xjtlhPn07q8O8pv5Jq3PYtCsZLVWVT8a1Zyc5QP68OMTlMcbhej4F7x/OWZj6J5tRz++0G7NS0Qs5hwBKziHLd9AunEOmoIgy7ZLUiu9kIm+XOCZOP0wt3w3MKlVtksUa+nvx8uK/dbtIIP0A5QvL3V9OWnoU6hN8FGw6B9kIQNm/DotpsdOrI5dPHH2pg8uSHph+mjkHd6vUyGvf9L1VXI8ufyZjNgnWPd9sECdms5xz0Jen8nM4eaJaD7Zw7Y7E4Tim0DMidHNMNo/VDSkeY80W6bZgrRR+a7OlVqAbWSeJlTnB8MShQojEMSqlcNAdq7fuwdD9ajD/h2cm8kR5I87BvPAm78kHsgcNC/6JOwbdZJDPFuCQd4tiu0xtlBdTE+egcjZTgVt2tHTkcNioepx6qHoSAvFFQtdM3o+yyE9FmF022ianN/GAGI4rPzExqD+XT+0/YkJrR76wm3RxHDNBfNwmipooid60ynCK69b0Yr/4HdEhTQavy9RClZ+DThSoqkt9lonYBsJ/fOowrYc2P88helbPMeXzrOCfoEOBEwq/tUsAxRNveCGWzr+DJw4VhtpaKeIcJgztV3A6UUGMymoS2aju3HjBEYk0rgvpiJm66cUMKosqOb84oM1ipfhvF52DDdyqpaU9ZyUoYjtNMaqA5PGJtdUK4uCwZzfpHFTfnSvp23Pm4HeiuaMp31lTRhTawWX6PJaRjLSe7oA9+ql8V885SJEANMXK9amMKeS6XKyV4o6Z2uzJ51WHLQr1TR5dj69+VH3aYZBXPvtFnY+HKq925Rxy9g1llpJnYwCImcx2BnodcTDpHMSFz7bzEids1mgJFL93woFDMHQ/tSLKFVkivLo8CsrXXyFDBuLEwhSET56kJgJgEjmp6m5uz1nDDMeIgzZP8F/mHFQLhosYxrTQqG6J3JRpsf7BZyYXrk1r+jEHDALAQ0gH71SnIHSAu1hJXJ9sQRflb64jJrtaOvDq8s2F39pYWik8sqsciAOfN6L+z9V5D7CfglinIcSF+mE/W4WnmwINys/z97HpW1Tr/2fCWFGlGlO4ovcppAvWSiqxkrjzsuwEsmJe/YeWzRP7qERKSLcrIkLMf+DOi6crwwXLbLSpPBE6k0ogOsbQBj4ZmtvyVs5BrE1/pq9arNSqUHi7iZX092wOhqby3ZXHFLZD4Bx0Ig6hzH89/WBtmWI+K+cg3S413Ia8GH/91In6vAWFtLnMDMV3yWk2zKrFWhzXKo4zBsnow6SIz+eZtb9lhbT5uNekQvqlb59S4Jpdz1svFSVxDkS0iogWEtE8ImoI0wYT0bNEtDz8P0jIfy0RNRLRMiI6S0g/JiynkYhuoTQntaeEK+dg2wm5EhIxRMfHDh6GG89PipSAdArSTCZ+0FBtdRbDFdYU4gBMoxcxvY8z5xAW2dqesz8jVK83lQz+y8ThohPGK/LGC5EPGALS777EPjH2ZYyI2HfGMc5BQ0TFz3GlYdGNPWOZQq6cQ/I5dbrYJ2MG9sG/nXmItoyCn4OljUTSudSaT6ZSKqvKFjlYnQivUDeAPYJZqk2sZFdIB/87CgppfV6VP8b4of1i54l3BsohVvoEY2waY4ybo1wD4HnG2CQAz4e/QUSTAcwEMAXADAC/JSL+hW4FcBmASeHfjDK0SwlT+IyY05jjTgAwK3DFcq6ZcajWtjkNskQxE1G9KMaNOMh9sZ8mzAYQnVRlQ8F3osMuVjK1RS5P1Dm8+p1P4PJTDkrkFfPc9o/H4OlvJRXMRrGS4l5MrGSRF3OYeorfY4zhh48vBqDnHOL+Kvoyn1n0YeHaxP0BxXMO2u/jaPwARHPChXOYszI6elMnVrrunMmJNFUbxM2fjXMgijtwanUyFIRyt5l2R+M35BwsmzWROHw0DPuuCvtdSVRC53AugHvD63sBnCekP8AYa2WMrQTQCOBYIhoFoJ4xNosF5gH3Cc+UHZGtkkomGV27DnBbXlc9hstJVhxb9rTFdtB6eWh0nYZzUNnBc9Rk3eLI8yLbc6wQMMwFuh2dylppnMZb99jxgwvXR+0/EMP7J7kqk/xadU/83qZPJW4gzeKnZFk6ziFGcAyFflYI8S4HlUzWHy+nVGY967gREe+7cA4idAYiqvqUnIPQJy46hz1tIuegEysFhCFn8Y3inGd7h92UVQ4t3q+mKqyrZxEHBuDvRPQ2EV0Wpo1gjK0HgPA/N98YA2CN8GxTmDYmvJbTEyCiy4iogYgaNm1KnpLm1GCmd0KJ7bQtA1dkI806B5E4uLYyHSYMUx9PGDNltSjAROgC9AGwmuypyhxvOz5RGOt6c87gv1OwQ4eFyijyUWwcRJPHx+evc6rbReEq7hB1Czovxra5FyPk2jiHYlUMOoW0q9k04M45iLfPmjICP7/gSHXdKv2hSucgEN8R9WYOnohivhpaP4dQIf0/zy83lpdJK1YSCEDk9d6ziMNJjLGjAXwSwBVE9DFDXlV3MEN6MpGx2xlj0xlj04cNs8fPUZZhuCcOAJMFEuDOOYiD1LRYXPaxyKxuSL8arLrhU8b6ORb/51laMZC7ziH+2zRx3HUOUaE6K5w0kMVK3O9ABbGbdcTml59XLzRXn3UI9h+S5EhEQwITl+cqVopMNQ2ZCnmDzLZFV/SfsZkcF8spHKYI0AfIOrjycA7i/S8edwAGKXRHgPpdVASjRpiLNvEuIR4HSfeZ+CZx2952AMBt/3i0Ol/BWslBIS2JlXL5uPlrZ4XPKGnWMsbWhf83AvgzgGMBbAhFRQj/bwyzNwEYJzw+FsC6MH2sIr0yMOkcUnEO0X2zX4DbAj1j6qjCxLtgevLAGxUyBPTRyKmB9MpRjkNHqhcAwN3PIWYplYJl0g17UcENQHv8JhDvZ51sVxcL6woN0XFdeON9rs1W6HMX80w+Lm0Lep3QxrTWSi74n5nT8KkjRinvxfrcUnhkymquT7yvskpT1V14VpEm6r6sYreMHHXZHFCP4xOKsN5iG7kpq03nIEqpZM6hs86QLpo4EFE/IurPrwGcCWARgMcBXBJmuwTAY+H14wBmElEtEU1AoHieE4qedhHR8aGV0sXCM2WHyZRVHNSL1u5I3BdR5TgZRCsX26RZun5n0EbLt+fHhw7uV2MZZOZyVLCd1ubqIS22qzpFQ4Zodofy5DKJTdI4g7lCJMLuXJiJIAf/f/1Co7XujONiKrbRRsR5H335xPHOXKpp0xDj1ixcN99tp+EcTOJE1UZOlRYT8Vrqrq3KSEYf6vzy/NNZ+skRbk1zIpOJ+3cUjpEN29xZpqyl+DmMAPDncJGtAvAHxtjTRPQWgIeI6FIAHwC4EAAYY4uJ6CEASwB0ALiCMcZ7/3IA9wDoA+Cp8K8iMJmyit+raZvewxOQd0r6fOIi5uoo1Lhxt1PdpTrTqWAz8XTnHIQFOgXnMF1QJoso+DlYzj4I8kbX5TguE4iLlUx9EOco7QTsvlmrrXXzIm0xv8T67MQh+J+GgzAt+q66lqBtnBMy1yeORFOkZFK8qs1ayTYkZfGpSedgq1fMx6MUG7lPIqzcvKfwmxMDvn50ls6haOLAGFsBICG4ZYxtAXBa8gmAMXY9gOsV6Q0AphbbljTg3aoOpU1CPvMHeHbJhkSZKqTZrRTyOcpsh+yn3mVzFDOGbAPP1SxVfAWbctS1mRmKFgnXhcq1z23oWxNNFee6y8TVuTrWxYwkrKasbnJ/Ea6OlLYxwhdIm5hMDH/NPcpVUHIOiramEX3J76DL7zofshLnYPo+WSK8tyGKF9VR2BB1LnHodeEz8kZrJfdy5BDBOoi7ZoufTPSMI3EY3M/MORTjZm97pBjOwXX3/rt/+ojxfnC2gDmCabF12+AuVhJFWgbuRuJdP62R5QdlurQwnQiNZ7UtblefFTmzXmf6EwAAEXZJREFUmTjAuAGC2Uz0xWWBpeFGxzn0jdMmGS3eVN9DdY5ImrEgi0914811oeb9c+tL7wMwG3YQxfszJ1sr9QSFdE+EWazkPnhs1kwcomzRdRdrE8PwRcDkrAYUp7iyi5Xc3iFuMWR+H67sG26zIKFISegaSLBczvYxcZHDQTaA2+l7HLrzAgB3caRY30BD0EgRNuuzc8J4PoBFTi4SB8fdtKsuao/17O5kmsoHRhwzaTjkyaPqCycgynA5UhhIEiYboRX1HT1OId1TUehWJefgLrMVB9qYgepTogCJc3AVD1iy8cWixpKxmDFkIw7u4aPdZO+xsi3RVEWxkqvFUBq4PuYqLjKLYdzb+M4H253yxcxJLeOXL442ziGuaHbT8+g8vWW4zgddYEkOm/K5kE90ZLTUKfbLTz93OIZo9Huu3LnL6ZIc8vt2lUK61xEHzjqoFiJxjFkHpDDQTGx0VREKade4TjYCVoxY6f1Ne6x5LjhmLG6/6BhjnjQK6YL82zIaMyLn4Cj3T4NnvmVy04ngejZGGp2DqxWYuV3uZfD1xVav+A7Kc0hCiH3+RuNmbT4AuOFzhzu0MMIlivhZIlzFRSIHa5sa/QQdk4nj150ZLUM+PrTD4FEtc5GiQprIcw4VQyF8hoVzOEBx+LcIm6iEoyZmIeE6iG06h6BMW4TUNCE50uAXFx6JM6eMNOZx3UEDgte6hXMQD54x9X+xaoZJI/ob7z9+5UkA3ENHm9sYb6QpENyRDudVB/W5vzjfOKQ5a8O0CdpfEOPssoiB0lrZ6ZzfOFy5MHGPYts4jRakAaZ5O3qgXhwoYuqYuBmwafcvh5sRRWBZIs85VAqupqy/0Ljqczgv9I7nAMTbYRevAHbOweEo7IpB5H5sCxAf63YLm8iU1dXPwRX8cCITJg7fD0P3q8F3zlLLn2Wk0TmYvMi/fOIBJdcnI1+EWMmEkQOiRbLWMi4LmxqHsqeM1vtWpIUYudi2vA7sGy3Qpn694Bg3h9X+ddWx6LGqSMEcclBAHnIDCOaVV0hXCFFspeQHFxeVAX3NCj3uKfoNSwjlmLOcY2/bRFp8aNh0DsUMois+kYxyWgzSWK/wXZyNeGYcZerFSJVcuqpvTRUa/uMMfOxgt9AtLjoHHqrkS8frCYC7Kav7i/OxkeYgJlfYIp666qAWXHcm/nT5ianr16G+rhpfOm7/4Iflg9c46gqJKMY1mSAO2akGZ1O+meIWcgJtQFVGfUpcJdALD/sJoPrcaZynvn3mIfjKSROsMVpix346TrSrztDHwgeixdSmmzh4RHBO8CUnuO08AeBqx12xDeJuy0Yc+Dy19fn2MH4NYF54izFfjR2FWSYYdRPh/741VThsVK3V8swFrqJOAGjvYIX6TSiKOFginu5sDr6j7Z3TRPMVcehIvXiwcMiSpYw0p0K2O7LovJzjJqgdPTl4/w2vr8XqLXsx8yNR1KHOFCv1PuJgiK0UUyY6xFZyOZuBT0KXMjl0p8Vx8NDFu1vMst1DR9Zj3g/OwMC+ZpltJSDuvOxipaCP0ohFXP0cXOES7TUtXNrY1pG3cop8XOriGkVlBv9dCA0Xz5mUzEBxToQ/tSicTzhoKD46aSiu+8yU1GXb8PA/n4DpBoc5/jq23beoz7P15/odLQCAI8eaQ8/wb36gJooyh8g5rPjp2YmAoF4hXSGYznMQUS7nqcmCzLRYE0sZZ08dhXGD+ziJN7qCMADprLS4iCNNn5fL+5gjjamhK8wWVcH/1o68c6gLm4FBv9oqjKyvw88crIH4+9qIgyo0hQ1HWOJzDehTjfsvPU4Z/bZU1GQzlnDsAeycQ5CztirjfECXeJ6GCrzOkfV603fg/7d35rFWVGcA/31vYXks78mmIFDQWsIiiyCVCNqauEDU10RrtFZQ2lgTrdbGtoptTdqYqGlNRZtSYqHazabpIlat1aZqXShK3UCKKDUVQa0BWbRC0a9/zBkY311m7ntn7sy8+/2SmzvvvHPPPd98M/eb7zvnfOdgqO9DVTdD6aMPreY5pIRWG5GO4Ms4+GonyrhhA/jb10/y3q5Poj94h8SM34QqqSUsUp1k5/yIYQPY7HLYpOE5VAuLhPf7vv0fxGdQdfIkWb2+eknZzDUlhJ5D/1b/YaWkC0TTIGnakLhzGXq+I2L2fYhyRmTBYDnCxXwdMfdD6DmUW6hXbn/ptGg4zyEk7ppP4Tc9lrsvPZ6bzi6/x3TRCFcRt/Vpjo1rh/gypEmn8P7BTU2F4EfaN9WeyqOJBOMHhYN3n78JO934zbCY/FzdUYmvTLjdIW7Fd3jfJx1ziIswwMGdBystlAsJjcPg/tXvh3ARYTnj0NokiVdl95QG9ByC9ziV+0q7AMHGNG/uej+23rQxHYnntOed8AnOV2K3Wkjqdkef7H16DqMP6c+WHf+teg2F52Xv/g8T773g84lx+pgO1ry6PXbNQbc8hwyNQ1yIbs4RQ1n5+KtMjRkfqEWGuy4+LtHMwPBHPW7A/oDnUKbNvq3NqYRAy9F4xqHKfg5pcdWp1Wcf9UYOPHnVcJp9/ah0J2ulzxtu1WVzeXtP9aRyB2LfGp+fZ84RwxjV3o8vn3SUpx7C7RfO4s2d78eOByXd+S9KtdxTaRMXVjpl8mGJJmmM6uiPCFx5cvw5b2oSmhJ4GAdSlsSc09B4fFDmuujb0pTKzLpyNJ5xSOg5GD0j6VaQUertOUTZ53HF4JABfaoucoKPrnl5LCbdRHtbK09ck2wsISmD+7UmmiranUkUviZedIdBCWRKMkljyIA+bqaQP1mS5rMK752tO0ujDX1bmmL39fBFw405VEuf0ch8p3MyP41JmV0L4QU+roYZKb48h/Eu9UkSY7Pm2uBHd2+dbriQJLHsohHdBz0rBsdMC64F39GFpMYhXAdSjr4tzeY5pMWCKSOZcNig2Lhfo7EwJrlZrbT3b+VH5x/DsTELfqL4uhnb21oTb3051O2JMbuGfvogKmpcqvKs+fbpkxjZHp9DaMmCiSxZMLEOPSrlprOm8ucX36hruLhWwjQYccah2v+bmmD15u2oauqyNpxxGDu0LXZ+dZI8O0Y884+uvmgrZHC/FnbFLOiL4jNq0dwkPPTVEz6SaK0eRMNtj3zt03X97lpZPHd81l2I5Zxjx3BOZCVxHgkzBcRdvmdOG8Xbe/aWlWf15u0A3L/uDRYkvL+6S26Mg4icBtwCNAO3q+oNWfRj0/XzvW0taSTj3svnsX7rrsT1Pxfmx/HEx0dUz8aaBtFLLG5FvNE7mDq6nfVbd8VOo21qEr44r3qIbsd7+/x1rFI/Uv+GBIhIM/BDYD4wCThPRCZl0ZfW5qZMB9QakTFD2jhtSvUU4ADLPh/sIXHk8IFpdyl1wtlcE0f6yzpq5JvrzpjM0vNmMMPDdPV6PMDmwjgAs4GXVXWzqu4D7gI6M+6TkTNOnXwody6eHbv5SxEIF5/Fpa8weg/9Wps5c9qoHo0V+Paaq5GXsNLhwGuRv7cAn8yoL0ZOEZHE6bLzzpHDB3LVKZ+gc3r1fDyGEWXJgols37OvLmNkeTEO5UxpSWhORC4GLgYYO7Z+FtQwfCMiXOZxUZvRGAzs28KymC16fZGXsNIWIDo0PxrY2rWSqi5X1VmqOmv48N7xBGkYhpFH8mIcngKOEpHxItIHOBdYlXGfDMMwGpZchJVUdb+IXAY8QDCVdYWqrs+4W4ZhGA1LLowDgKreB9yXdT8MwzCM/ISVDMMwjBxhxsEwDMMowYyDYRiGUYIZB8MwDKMESbrfbt4Qkd3AxjL/Ggv8O0ET7cBOj/XSaDOpLGl8dxpyN6o8WV5DtXy33Tv1rVdLXZ+6maCq8dkmVbWQL+DpCuX/Sfj55T7rpdFmUllS+u405G5IeTK+hmr5brt38iu3N91U+u3s+uqNYaV3Eta7x3O9NNpMKksa352G3I0qT5bXUC3fbfdOfevVUjeNa60qRQ4rPa2qs5KWF5HeJAuYPHmnN8nTm2QBv/IkbavInsPyGsuLSG+SBUyevNOb5OlNsoBfeRK1VVjPwTAMw0iPInsOhmEYRkrk3jiIyAoReUtE1kXKponIkyLygojcIyKDXXkfEVnpyp8TkU9FPjPTlb8sIkulJ9sx9QCP8jwsIhtF5Fn3GpGBLGNE5K8iskFE1ovIFa58iIg8KCKb3Pshkc9c43SwUUROjZRnrh/P8hROPyIy1NXfIyK3dWkrU/14lqWIujlZRNY6HawVkZMibaWjm6RTrrJ6AScAxwDrImVPASe648XAd93xpcBKdzwCWAs0ub/XAHMINha6H5hfcHkeBmZlrJuRwDHueBDwEsEe4DcBV7vyq4Eb3fEk4DmgLzAeeAVozot+PMtTRP0MAOYClwC3dWkrU/14lqWIupkBjHLHU4DX09ZN7j0HVX0U2N6leALwqDt+EDjLHU8C/uI+9xbB9K9ZIjISGKyqT2pwNu8EPpN238vhQ546dDMRqrpNVf/hjncDGwi2fO0E7nDV7uDgue4E7lLVvar6L+BlYHZe9ONLnvr2ujK1yqOq76rqY8D70XbyoB9fsuSFbsjzjKqGG6CtB/qJSN80dZN741CBdcCZ7vizHNxF7jmgU0RaRGQ8MNP973CC3eZCtriyvFCrPCErnVv8rSzCMFFEZBzB083fgUNVdRsENwGB1wPl9wo/nBzqp4fyhBRNP5XIlX56KEtIkXVzFvCMqu4lRd0U1TgsBi4VkbUELtk+V76C4OQ8DfwAeALYT8I9qjOkVnkAzlfVo4F57nVBXXscQUQGAr8FvqKqu6pVLVOmVcozwYM8UEz9VGyiTFkm+vEgCxRYNyIyGbgR+FJYVKaaF90U0jio6j9V9RRVnQn8iiDWi6ruV9UrVXW6qnYCHcAmgh/Y0ZEmyu5RnRXdkAdVfd297wZ+SUbhDBFpJbi4f6Gqv3PFbzp3NwxJvOXKK+0Vnhv9eJKnqPqpRC7040mWwupGREYDvwcWquorrjg13RTSOISzC0SkCfgmsMz93SYiA9zxycB+VX3RuWe7ReQ450IuBO7Opvel1CqPCzMNc+WtwOkEoal691uAnwAbVPXmyL9WAYvc8SIOnutVwLkuVjoeOApYkxf9+JKnwPopSx7040uWoupGRDqAe4FrVPXxsHKquvExqp3mi+BJehvwPwIr+QXgCoLR/ZeAGzi4mG8cQabWDcBDwMci7cwiuAheAW4LP1NEeQhmYqwFnicYnLoFN0umzrLMJXBhnweeda8FwFCCgfRN7n1I5DPXOh1sJDKrIg/68SVPwfXzKsGEiT3u+pyUB/34kqWouiF4aHw3UvdZYESaurEV0oZhGEYJhQwrGYZhGOlixsEwDMMowYyDYRiGUYIZB8MwDKMEMw6GYRhGCWYcDCMFROQSEVlYQ/1xEsnUaxhZ05J1BwyjtyEiLaq6LOt+GEZPMONgGGVwydD+RJAMbQbBAsWFwETgZmAg8DZwoapuE5GHCXJfHQ+sEpFBwB5V/Z6ITCdY9d5GsFBpsaruEJGZBPmz3gMeq590hhGPhZUMozITgOWqOhXYRbC/xq3A2RrkwVoBXB+p36GqJ6rq97u0cyfwDdfOC8B1rnwlcLmqzklTCMPoDuY5GEZlXtODeWx+Diwh2GjlQZfluZkgFUrIr7s2ICLtBEbjEVd0B/CbMuU/A+b7F8EwuocZB8OoTNfcMruB9VWe9N+toW0p075h5AYLKxlGZcaKSGgIzgNWA8PDMhFpdfn1K6KqO4EdIjLPFV0APKKq7wA7RWSuKz/ff/cNo/uY52AYldkALBKRHxNkybwVeABY6sJCLQSbMK2PaWcRsExE2oDNwEWu/CJghYi859o1jNxgWVkNowxuttIfVXVKxl0xjEywsJJhGIZRgnkOhmEYRgnmORiGYRglmHEwDMMwSjDjYBiGYZRgxsEwDMMowYyDYRiGUYIZB8MwDKOE/wOypFk5Txw8PQAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Un zoom sur les dernières années montre mieux la situation des pics en hiver. Le creux des incidences se trouve en été." + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEVCAYAAAALsCk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXeYJGd57n0/nXNPzpvzrrRKK2mFEkIGRDAiSLYAG4GxdcwFfD4OB6Nj+7O/jwMHH4MxJvlgg0WysGw4SGAJSSgQrLS70mpXm/Ps7KSe1DlWv+ePqre6OnfPdE9XV7+/69prZ6qruqu7eup+n0yMMQgEAoFAoMXU6hMQCAQCgf4Q4iAQCASCIoQ4CAQCgaAIIQ4CgUAgKEKIg0AgEAiKEOIgEAgEgiKEOAgEAoGgCCEOAoFAIChCiINAIBAIirC0+gSWS19fH1u/fn2rT0MgEAjaigMHDswxxvqr7de24rB+/Xrs37+/1achEAgEbQURXahlP+FWEggEAkERQhwEAoFAUIQQB4FAIBAUIcRBIBAIBEUIcRAIBAJBEUIcBAKBQFCEEAeBoANgjCGbFVMfBbUjxEEg6AD+/cAEbvjsU5CEQAhqRIiDQNABnJuLYiaURCSZafWpCNoEIQ4CQQcQS0kAgKgQB0GNCHEQCDqARFqIg6A+hDgIBB0AtxyEW0lQK0IcBIIOIOdWklp8JoJ2QYiDQNABcLeSsBwEtSLEQSDoAGIpWRREzEFQK0IcBIIOIJ7OAgCiKSEOgtoQ4iAQdABxRRSEW0lQK0IcBIIOIC5SWQV1UlUciGgNET1DRMeI6AgR/YGyvYeIniSiU8r/3Zpj7iei00R0gojerNl+DREdVh77eyIiZbudiP5V2f4iEa1v/FsVCDoXka0kqJdaLIcMgD9mjO0AsBfAR4loJ4BPAniKMbYFwFPK71AeuwfALgB3APgqEZmV5/oagPsAbFH+3aFs/zCARcbYZgBfAPDXDXhvAoFAIS7qHAR1UlUcGGNTjLGXlZ/DAI4BGAVwJ4BvKbt9C8A7lZ/vBPB9xliSMXYOwGkA1xHRMAAfY+x5xhgD8O2CY/hz/TuA27lVIRAIVkZayiKjNNwTbiVBrdQVc1DcPVcBeBHAIGNsCpAFBMCAstsogIuawyaUbaPKz4Xb845hjGUABAH0lnj9+4hoPxHtDwQC9Zy6QNCxcJcSICwHQe3ULA5E5AHwAwD/lTEWqrRriW2swvZKx+RvYOzrjLE9jLE9/f391U5ZIBAgVwAHCMtBUDs1iQMRWSELw/cYYz9UNs8oriIo/88q2ycArNEcPgZgUtk+VmJ73jFEZAHgB7BQ75sRCATFaC0HEZAW1Eot2UoE4BsAjjHG/lbz0CMA7lV+vhfAw5rt9ygZSBsgB55fUlxPYSLaqzznBwqO4c91F4CnlbiEQCBYIbw62mO3CLeSoGYsNexzI4DfBnCYiA4q2/47gM8CeIiIPgxgHMDdAMAYO0JEDwE4CjnT6aOMMb5c+QiABwA4ATym/ANk8fkOEZ2GbDHcs8L3JRAIFLhbqd9rx2Is1eKzEbQLVcWBMfYrlI4JAMDtZY75NIBPl9i+H8BlJbYnoIiLQCBoLNyt1OexYWIx1uKzEbQLokJaIDA4vMah32tHWmJIZuqPOyTSEv7Ld/bj/Fy00acn0ClCHAQCg8NbZ/R57ACWF5Q+ORPG40dmsO+8yBPpFIQ4CAQGh7uV+lVxqD8oPRdJAshPixUYGyEOAoHB0bqVgOUVwgXCsjjEhTh0DEIcBDXx0L6LePbEbPUdBbqj2K20HMtBznJKKHMh1OdOSbjj736BRw9PrfAsBXpDiIOgJr70zCk88Nz5Vp+GYBnEUxJMBHS7bQAaazn85NAkjk+Hsf/84spPVKArhDgIaiIYS6t+Z0F7EUtJcNks8NjlzPXlBKQDZWIO33txHAAwG06s8CwFekOIg6Aq2SxDOJnBXFgUULUj8bQEh9UMt13unL8st1K4WBxeuxTEwYtLAIDZsFg4GA0hDoKqhBMZMAbMR5Mo19UknEirrgeBvoinMnDZzKrlsCy3kmo55GIO//LSOOwWE27e0ieuvQER4iCoSjCeBgCkJab+rIUxht95YB9+54F9q31qghqIpSQ4rWa4VbfS8i0HnvmUymTx8CuX8PbdI9gy4MVMSLiVjEYtvZUEHU4okROEuUgSXS5b3uOPvTaNfecX4XWIr5MeiaclOG1mWM0m2CwmRFL1iUMyIyGUyKjPBQBL8RSiKQlXrvEjmpIQS0mIJDOqdSJof4TlIKiK1loIFMQdkhkJn33sOADZ/RSr88YjaD7xlASXTY43eOyWui0HnsYK5GIOMSWo7bJZMKDUT8wK68FQCHEQVEUrDoUZS999YRzjCzHceeUIAGA2JHzPeiOelt1KAOC2m+vOVuIuJSKNOKS4OJgx4HUAEEFpoyHEQVCVSuLw9PEZ7Bj24a5r5DlOwvesP+Ip2a0EAG5b/TMdeLB52OdQA9LxtPwcLrsFAz7FchDiYCiEOAiqUkkcFqJpjHY5MOgTq0e9wgPSwHLdSvI1HetxqTGHaFJrOQi3khER4iCoSjCehsVE6Pfai2odlmIpdLlsGFRcC8Jy0B/xdC7m4K5DHHhrb1Ucup0l3Up+pxU2i0mksxoMIQ6CqoTiafidVvR57EWWw2IshW6XFT6nBXaLSVgOOkR2K8lZRLWOCj09G8Zlf/k4DlxYRCCchM9hgd9pVS0HnnjgsllAROj32MW1NxhCHARVCariYMsTh3hKQiKdRbfbBiLCgM8uLAedkZGySEnZugPSr4wvIS0xPHLwEuYiKfR57XBYzSUtBwAY8NlFCw2DIcRBUJVgPA2f04p+jz0vrZHPI+5W6h4GvQ6RraQT4ikJs6GEutKv1610OhABADx5dAaz4QT6PXY4rWakJYaMlNVYDoo4eO3i2hsMIQ6CqqhuJa8dgUiuhUZOHKwAgEGfAzNi9dhU9p1fwKGJpar7fenpU3j7l36lVjQ7NHUOkVSmbBsUzplZeRzoZDCBVyeCiuUg3y4SmazGcpDdVQNeh3ArGQwhDoKqaN1KqUwWYWXluRiVs5i45TDgE6vHZvM//uMY/ubxE1X3u7QUx2w4iRMzYQCAS5OtxBgQTVV2LZ0JRHDdhh4Qya0yuOUAyLUO8ZQEu8UEs4kAyJZDMJ4Wk+IMhBAHQVVkt5JFHRbDi6JUy0GZEzDgdSCSzCyrd4+gNuKpDJZixf2tCuH77Dsnz3zm7p81PS4AwPm5aNljkxkJ4wsxXL+hB9es7QYgT5GzK+IQT0mIKs38OLzWQWQsGQchDoKKMMYQSmTUbCUg105hqTDmIIqhSsIYw78fmFBdPCshmcnm9boqx5JSm7JPGcLD3UrbhrwAgOPT4bLHXpiPQcoybB7w4E27BgEAfR6bajkkM5I6I4IjqqSNhxAHQUWiKQlSlhWIg3wDWFDcSl2amAOQX+uQkbII1rDSNTKnZiP4k397Fd96/vyKnyuZzpbsjFtIUBHuVy7K4sDdSut6XLBZTDgxHSp77JlZORi9qd+Dt+0ewVi3E5ePdqniEE9lEUtKeZYDn08dEDEnwyDEQVARfiOSA9KyhcDFYTGWgtdugdUsf4245aAVh88/eRJv+rufr+Yp646w0tH04YOTK36uZEZCKJ6uGlDmlgNvd8HbZ1jMJmwZ8ODETCRv/1gqgwdfGkdayuK0Ig4b+90Y7XLiV3/6Buwc8cHBxSEtIZaW4NJ0YB1Qr72wHIyCEAdBRfiq3++0osdlA1Eu5rAUS6nxBgDo564F5QaRlrJ4aN9FzISSVQOVwVgaf/3T40hL2Yr7tSPcnXRsKoSTM+XdObWQzGSRrRJQzmbluRvcogOQt8rfNuQtshyeOR7A/T88jB++PIEzgQhGu5x5biMAcNqUbKW0hFgyo1ojANDrtsNEIuZgJIQ4CCrCLQef0wqL2YRetw0BJeawEEuraawA4HNY4LCa1GKoX5wMYD4q7xuq4gr55ekAvvbsGZyo4AtvV7RtzB9ZofWQzMjiWenz5JP7XrepV93m1Nzotw16MRNKqjEjQJ7yBwD/8POzODkTwcZ+d9Hz2i0ayyGV71YymwhumwVR0bLdMAhxEFRE61YCgD6PXV0d8r5KHCKSax0Uy+GHL18qep5yJBX3R6wBQVu9wQvRxrqdePjVS1VdQuXISFlIWfnYSp/nUly+6d+wsRckZ5qq8QKgdFB6QRHxc3NRHJ0KYfOAp+h5uWsqkZbkfk0Fg30cNnNDgu4CfSDEQVARvkL1OWRxWNvjwvl5OQ1yIZpCjzt/Ktyg14GZUALBWBpPHpvBJmUFWlUcMlwcjLfy5DfM39yzBhcX4nhJSS+tF/4ZAZUtB57GOtLlxMY++fPXrvK3D/kAIM/FtRRLw2u3qPtv6i8WB4emziFa4FbirxEXdQ6GQYiDoCI8bdKvuI+2DHpwfi6KVCaLpVi+XxuQA5OvXFzC3f/7OaQyWXzwxg0AahEH+aZixCIqfsN819WjGPDa8ac/OFRTxlEheeKQKC+iPBjd5bJi14gfRIDdkvtTH/TZ4XNYiiyHXo8Nv//6TQCAHcO+oud1auoc4ikJLru56HEjWn6dihAHQUWC8TRMBHgUn/WWAS8yWYaTM2FEkhn0FMyTvu+Wjfj13SNw2y146+VDuGlzn/o8WqQsw8RiTP09ZzkY7+bC31Ofx46vvP9qTCzG8Sf/9iqy2frcS1xAgSpuJSWW4Hfa8N7r1uJ3b9oA4v4lyO6/7UO+vPjOopJccPc1Y3j4ozfi6rVdRc+rikM6W1QEB8huJyOKe6cipoELKsKb7pmUNgncF73vvOwa6SpwK+0e68LnfyN3Y+G+7MKb2U8OTeK//dsh7PvzX4PfaTV2zCElwaSs3q9d34P737oDn/rJUTz22jTetnu45ufhaalAZbdSUGM5bB7w4AZNYJqzbciLH70ixz+ICAvRFAZ9DhARrlhTLAxAzvoIxtPIMhRnMwnLwVAIy0FQktlwAsmMpPZV4mwe8IAoJw7dBW6lQnwO+QZSKA4zoQRSUla9yRndrcTnHgDAe69bAwAYX4hVOqwIreVQqUp6KZafRFCKrYMehJMZNXlgMZpSK93LYTIR7BYTFpTMpkLLwSUC0oZCWA6CkrzrK89hrNsJq9mUd5NxWM1Y2+PCS+fkyttCt1IhFrMJXrulSBz4TAEuBkZ3Kzk0wVun1QyziRBJ1hd3SOZZDhViDrE0PJrixFKMdDkBAFPBOIb8DizG0uhxVxZ6QL7+vDLeXWA5OKwiIG0khOUgKCItZXFpKY4Xzy3gV6fnilagWwY8apV0VxVxAOQaiUJxiBeIAl8VG1Ec4gX+eSKSW2dXCCqXQhuQrpbKWslqAPJbncRTcmpqt7v6tXRazarl4CxhORgx26xTEeLQgcRSGXz8wVdwdLJ0f51FJU4wqqwueRorZ/OAV/25u4bVpt9pLfKR85uIajkoq2KjupWcBWmfHrtFbatRK7W6lYIlssgKGfLL4jAdTKjddatZgQDgsJrUwsaigLRVuJWMhBAHDZ0y5vCbvzqHH786iWdOzJZ8nHdd/eRbtuOua8bwazsH8h7foimQquanBmRxKFzpcguBB1lL1TkYpWFfLCUVrbK9Dos6F6NWuIA6rebKdQ7x6uLQ47LBZjZhKpRQkwZqsQIdVjMWIlwcCttrWIRbyUBUFQci+iYRzRLRa5ptf0VEl4jooPLvrZrH7iei00R0gojerNl+DREdVh77e1Kic0RkJ6J/Vba/SETrG/sWa+ORVydx/WeewuRSvBUvv2osRlP43z8/CwBl5z3zVgpDfgc+d/cVeNdVY3mPbxmUxcFpNef50stRShz4CpOvhgvdSqdmwrjqU0/g1YvVp57pnURaKlplr8St1O+1V65ziKXQ5aweXB7w2TGjtRxqcCs5rGZV1EoFpNMSM2R/rE6kFsvhAQB3lNj+BcbYlcq/RwGAiHYCuAfALuWYrxIR/wZ9DcB9ALYo//hzfhjAImNsM4AvAPjrZb6XFfHt586DMWBeMyPZiHzt52cQTWXQ7bJiKlhGHJTPoLfMzYJXz1bLVOLUYzlwt9LEYhxZBhw0gDjEUsVuJa/Dgki9loMioANee9VUVn8N12bI58C0xnKoJSCtfR/uEkVwAIT1YBCqigNj7BcAaq33vxPA9xljScbYOQCnAVxHRMMAfIyx55ncWObbAN6pOeZbys//DuB2blWsFqdmwth/Qc6+SUnG/WLPhhJ44LnzePfVY9g91oXpMuLAg829yvyGQtx2C0a7nDUFMAG5urqc5VAYc+CiwVenp2bbvxFfvIRbyeOwLkMc5M9owFdeHBhjcuV6lYA0IFuGcgO+/HGvldC+D2eRW0lpryHiDoZgJTGHjxHRIcXt1K1sGwVwUbPPhLJtVPm5cHveMYyxDIAggOKqnSby4Eu5U9ZmhBiN770o9+v/2G2bMeyXV42lmIukYDWTWqNQil+/YgS3bu2v6XX9TisS6WxeQDWWVgLSZdxK3OVyqmDuQDtSPiBdbyqr/Nn0e+wIJzNqEz4t0ZSETJZVjTkAsuUwFYxjPpoCUeW6CI7DmrtluEsEpAFjZpx1IssVh68B2ATgSgBTAD6vbC+14mcVtlc6pggiuo+I9hPR/kAgUN8ZlyGRlvDDVybUzJyUQcUhlcniey+O4/Vb+7G+z41BnwNzkWTJ9zsfSaLXbUclA+6Tb9mOT9yxvabX9ik3Ha31EEvmWwyFbiVeA8AHz7Qzhe2tAbk4sP5sJW45yJlGpWIWvHVGtZgDIFsOiXQWF+aj8Cst2auRV69RIuYACLeSUViWODDGZhhjEmMsC+AfAVynPDQBYI1m1zEAk8r2sRLb844hIgsAP8q4sRhjX2eM7WGM7envr23VWo2fnwxgKZbGb+1dB8C44vDo4SnMRZK493XrAQDDfgcYK52hNR9NqVPfGgFfkWpdIWrMQbUc8t1K/MY3H01hPtLeA2Rkt1K+FeaxW5DMZOv6vvHPqM8jX5tStQ5qdXQNlgOvdTg2FarJpQTkxMFiItgKxITPqRaWgzFYljgoMQTOuwDwTKZHANyjZCBtgBx4fokxNgUgTER7lXjCBwA8rDnmXuXnuwA8zZbb8H4ZTCzK2Ul71suesZRBMy0eeO48NvS5ccsWWVQH/cXznjnccmgU/lKWg1rnULoITpvm2c7WQ0bKIiVli91KissuWkfcIZmRYDGReiMvVeug9lWqwUU0rHwHzgSiNScX8PfhtJmLLEuXpqW3oP2pJZX1QQDPA9hGRBNE9GEA/0tJSz0E4DYAfwgAjLEjAB4CcBTATwF8lDHGvykfAfBPkIPUZwA8pmz/BoBeIjoN4I8AfLJRb64WFqMpmE2EPiX4akTL4cJ8FAcvLuH9169VG+jxG0OpjKW5iNy+uVGUEge1QrpMEVwkkVFXpqfaWBz4+yyVygqgLtdSMp2F3WIqaYlxuOVQS80CtxykLKspjRXIxRwK3w+QczMJy8EYVO2txBh7b4nN36iw/6cBfLrE9v0ALiuxPQHg7mrn0SzmlYZjvOOkEcVhUblhaEc/DvvkGEthxhJjDHORpCqWjaBQHNJSFmlJNg5z7TNyRXCMMUSSGaztdWFqKd7WlgMXB0dREZz8mYTr6K+UzGRht5rVGE4py4FPgaslIM3FAagtUwnIWQ6FfZUAEXMwGh1fIb0YTaHHbYVNEQcjZivx1TmfAQwAPqc877lQHKIpCclMtmyNw3JQxUERKe3KMtd4T/4/y+RrEElm4HVYsHnQ29bprDxlt3BqmldxK9VTCJfMSLBbTCUD/JxaOrJybBaTGr+o3XLIuZXKPRYX/ZUMQceLw4JiOdgMbDkklPekTUMkIgz7nZgqiDnw4G8jLYdc2275phEvEAfGGJKZLLyKqyWRlhBOZOCxW7BlwNPW6axcCMu5leqpdUhmZLcS/zxLdWYNxtNwWE01Va4DOeuhFjcUkBOA0paDvE30VzIGQhxisn9ddSsZMCBdynIA5Dz3mQLLgfdVamTMwWI2wWO3qC6PqGZlmczILibGgC6lQjeWknKWw4AHs+Fk2/ZZKu9WWm7MwQyP3QITyW6lAxcW8NzpOXWf+Uj11hlahhRxqKU6GsgPSJd7LCbcSoag48WBDznhwU8jupVKWQ6AnOdeGJCea4LlAOS30Ci0HLhLid/UYikJEY3lALRvpXQ5txLPVqqn+V4iI8FuNYGI4HNacXo2gg/98z782Y/UtmcYX4hiba+r5ufk3VnrTWUtbJ0hPyZ/v0SFtDHoaHHIZhkWYyn0um0gkvO2DelWKmc5+B2YCSXyZhnPN8FyAPLbdufHHLKqIPMgaiItWw4euxUb+uQg+oX5+qam6QX+Xou6strl91pXzEHJVgLkNuqPvTaNUCKD8YWYKrDn5qLY2Oeu9DR55CyHGgPSNvn1ndZitxIRiVGhBqKjxYHPwuU9gmwWY4oDv/naCyyHYb8DmSxT+/MDuZhDrTeLWtFaDrzGwWUzI5GR1PPjq9dIMiOLg8OiTixr12655VJZHVZT3dPg5IC0/Dw84HzFmB9SlmF8PoZgPI25SEoV1FpY0yNbGQNeR5U9lfO2lLccAGVUqHArGYKOFof5aH6rYpvFZMjGe+ViDjwYqc1Ymo+m4HVYivZdKaXcSt0uG5LprHp+vBCLu7a8dgscVjP6PDZMBttUHBQhLKyQJiJ5pkNd2UpZ1XWzoc+Nq9d24a/esQsAcCYQwbm5qPpYrbxt9zAe/L29NbuieOykVMwBUEaFCstB5dJSHI8fmW71aSyLjp4hXdjH3qhupWSZmEOuEC6Oy8f8AOQbc3+D4w1AoeUg3zx63DZEUxn1/PyK5RAIy+LA/fIjXU61kr3d4DfKwgppoLaZDufmovA5LOj12JVsJfl5vvCbVyLLmPp9PROIqit2bT1LNaxmE27YVHufS245uEq4lQA+KlSIA+fbz5/HN355Dqc+/ZaKvcr0SEdbDryPPXdnGNWtlEhLIEJRLxwejPzHX57FQ/suIprMYC6SbHi8Achv282zWbrdiuWgupVky0EVByXdc8TvbFu3UqyMWwlQOrNWCEhnpCzu/ofn8DePnwCQq3MAALOJYDWb4LZbMOx34MxsBOcCUZgo5ypqBtxiEG6l2liMppDJMrVNTDshxAG54KvsVmq/i1gNnh9fuHLp99jxuzdtwMWFOD7xg0N4x5d/hXNz0Yb2VeJ0ueS23Ym0pLpaelxWJDOSxq0kX4fZAsthtNuJyaUEVrHlVsOIp2Rh5jd1LV5HZcth/4VFzEVS6ueRTGeL4kaAPHzpTCCCs3NRjHW7Gu4S1NLltMJiIvR7S39HhFspH74gqnd2hx4Q4gCN5WBQt1IiLZUsiiIi/Pnbd+L5+9+A73z4OgTjacyEmmM59CmCMxdJIqq06+5y2fKylfwFlgMvihvpciKeltTq33YirkyBK+VS8DqsFdtnPHFkBkCuDbfWraRlU78bZwJRnA1E64o3LIdutw2P/+EteNvlwyUfF5ZDPkux/CSMdqLjxcFly81BtllMhqxz0KZAloKIcPOWfvz44zfh7buH8eZdQw0/By4485EU4mnZPeKymZU6h/xspSLLoUt2f11qQ9dSrMT8aE6lmANjDE8clQOZfPWpdStp2TTgQSSZwYmZcF3xhuWyqd9TdvaD02Zuyxths+DXji+I2onODkhHU3kpm0YVh0SmtOVQyLDfiS+/7+qmnAMvqpuLJBFLZVRRzmSZejNx28ywmU0IKDMmPBrLAZDF4bJRf1POr1kkUuU/e0+FOdLHpsKYWIzDbTMjGM+oLUZKioMy01vKsrpqHJqB02ppS/96syhM324nOttyiOWLg93AAelKlsNqoLUc5MloFvWceHGc3WKGw2pSU4x5odhoG9c6lJoCx6mUyvrE0WkQAW+5fBiheBopKQvGAHsJoeHiAAAb+jxFj68mTpupLW+EzYK7laJtGIfpbHGIdoY4yPnxzQtS1gK3HAKRpDIZLefO46sru9UEl80CHnfmGTE9brn31aU2TGeNlZgfzfFWmAb35NEZXLO2Gxv73UhJ2dxnVELkB312dZ7zhlVwK1XCZbOImINCMiOpn0U9Q530ghAHV75byYjZSnqwHBxWuWFcznIwq3UX2hsfX2U7rWbVr01EGO1ytlUh3MUFud1HQhHCUpTrzBqMpXFkMoRbt/arldCzITkOU+o6EhE2DXhgt5gw7Kut0rlZOKxmJNLZvJYsnYq2pboQhzZjIZpSW2cAxs1W0oPlAMiupTnFcnDZzGrmTVDjVuI3Uh6M5ox2O3FpqXhqnR45PBHEzf/rGRy4sIhYOqO2si7E4yjdX+nl8UUAwDXru1Vx4Blc5dJUb9zch9dt6lUn/bUKLu58NngnU2pmejvRsQHpRFpCLCUVBaSNKA6JdBa97taLQ5/HjvloErF0BgNeh2o58LkENotJdcHwNFbOiN+J49Ozq3vCy2RamZGx7/wCYqkKbiW1M2t+Ouv+CwswmwhXrunCgQuyUMwqQfpSdQ4A8Kd3bG/Iua8UtW23ElfqZLSp19E2jMN0rOVQ2DoDMK5bKam0em41vW4b5sIpxJKyq8WuiTlYzQSzicpaDiNdTgTCybYYXs9dCIcnghXdSt4yc6T3n1/ErhEfXDZLTW4lPcHfqyiEy3crxdowlVXf37QmwltTa/vY28xmQ1oOyXRW7YnTSvq8iuWQkuCymtVzCsbTqruEuyU89mK3ElA881qP8FXiqxNLFQPSnhKjQtNSFq9OLOGadd0Act1Xee1HqWwlPcHfqwhK51sOokK6jeCWg7Ya2KhuJb1YDn1uGxaiKUSScp0DP6dQIq2uiPnNpVAcRpRCuHZIZ+WWw8RiHKF4ukIqqxJz0Nw4jkyGkEhnsWddD4CcOMworiq9Ww4uYTmocMvBYW3P9F59f9OaCJ+A1l0iW6kde/hUIqEjyyHL5Juh02YpsBwUcVD81IVuJT6UZjrUBpaDxoWQZeXbW3MB1Dbf239+AQCwZ71sOXABma0SkNYL2phDp7MUT4NI/u62Y51Dx0WMnj8zj3/4+Rn84lQAdosJg77PeJZXAAAgAElEQVRcAzHtHGm9/xHWQyKtD8tB29BPm8oaSWbUOgi+8iwMSPMVdD3zD1pFLJWBxUTIKOmc1QLS2qyWAxcWMdbtVGdtmE0En8OiyVZq/XWsBBfCdogNNZtgLAWv3QKvw4pYG7qVOk4cxheiODkTxsffsAW/sWdMXZkBuZbWqTINztqRjJRFJsv0YTloXHguTUCasdxNz1UmIM2vk/ZGqlciSQldLht8TgvOBqJl3UoOqxk2iwmhRO49vTy+iL0b8+cr+F1WNdZSOJNDb/AMJWE5yBZxl8sGl80sLId24F1XjeGua9bAXCIf3KbcoJKZLLyrfWJNotyI0FbQ69FaDhY4NKtgLg4ONeZgzTvWZjHBYc2/keqVWCoDt92M3aN+nA1EK9aYaGdrS1mGmVAS63vdRftcXJBjLXpftIiAdI6leBp+pxVuu0VNRW4nWn/HWGVsFlNJYeCPATBUUFqdAqcDd0S/p9CtlLvRFWcrFd8EfQ5rW7iVokm58O3ysS4AqJjv73NY1DoPLhLchcbR/t4ubqV4GwZgG41sOcji0I5dWfX9TVtl7G0uDom0hMePTEPStC7gvl89pED6nBZYFGF22sx5Nzpu2ZRzK8nHW9vCcogmJXjsZlyrBJUHfOWHJ2nHp/L/u1yVxKH117ESXByEW0lug+JzWuG2mUX7jHbHpglItyP/89Fj+C/fOYBHXr2kbis3P7oVEJGaOuyyyb2TuFhUcysBcgCXr7L1jNyS3ILdY114+o9vxR6lZqEUWsFbqsVy0MF1rIRwK+UIxtPoclrhslnaUiz1/U1bZbQBab1SLs324MUlfPuFCwCAf/rlOXU/bjnoISANoCgriYtBzq2kpLLaS1gODivCbWA5RJIZ9fw39nsqDpavxXLwOYuTJvSK2USwWUwdX+fAGMOS6lYyI5rKtF2KvL6/aauMNiCtR3748gT2/s+nkCxoapaWsrj/h4cx6HXgz966A0cmQ3jhrJwvr6eANJALSjut8s2TWzTccrhqbRfeevkQdo36io71OiwItUHModIMh0K04sDHgZazHGxmU8sb69VCv8eOM4Foq0+jpURTEqQsg1+xHBhD2w1B0scdQyfoPSB9dDKEmVASJ6bDedufPj6LY1Mh/MXbd+K3b1iHHrcN3/jVWQB6tBxybiUgZzFw8erz2PHV918Dn6PYreRzto/l4C5h+ZTC55CzlRhjmoB0/gxvLg56D0Zz3nLZEH5+clad0d6JcKHvctrU5Ip2a6HRHt+2VcKu85jDXEQuhHrtUihvO2+tcP3GHjisZvzW3nX42bFZXFyI6c5yUN1K9nxRqCXQKt9I9W2eM8YQS0nqoKJq+J1WtWo8WCXmoJdrWI33XDOGtMTwyMFL1Xc2KLyvkk+xHID2GxXaHt+2VcJmlv+g9Wo58PGZr00G87YvRhVftXITuXVrPwDgdCCSy1bSieWwrtcFp9WsWgbcoqllVex1WJCSsrp1+wGyG0/KsprbVfuc8n7BeBpLsTScSmGcli7FktDLNazGjmEfdg778IOXO1ccQpr4EV8otFs6qxAHDXp3K80pnWSPXMoXh6W4XKbPJ6d1KwHNpVgq51bSyarzN/aswc/++FY1EF0Yc6gED8zqOZ2VpyyWCqiXglsFoXhGzYsvt0+7WA6AbD0cvhTEyZlw9Z0NiDbzTFgOBiCXyqpPhedupWPTYaQ1rq9gLI0ud+6mwpsJLkbTObeSTladVrMJo11O9fdczKEWtxLvRaTfPzKeslhrQJoLXjCeVitqC8nFHPRxDWvhzitHYDERfvRKZ1oPwVKWQ5tlcAlx0KBnyyGbZViIprCmx4lUJotTMxH1scVYSnU9APINh0hevSTVIjh9Xuq6LAeH/i0HHnSsJyANyO8pWFUc9HkNS9HnsWPLoLcoeaJT4DEH3j4DaL850u3zbVsFeA65Hn3awXgaUpbh9VsHAOTHHZYK3BFyJ08rlmIpTRGcPleduTqHWtxKpSen6QnuOqhVHPwayyEYKy0OXocFRO0lDgAw2uXEpTaYv9EMgvE0bGZ57K3bJsSh7dGz5TAflV1Ke9Z3w20z58UdlmJy90ct3S4rFmNp3aWyFsJveLVmKwH67szKg47uWuscXLn3VC7mYDIRvHaLLlqg1MNYtxMTi3FdZ5c1i8mlOPq9dhCR6mJstyrpquJARN8kolkiek2zrYeIniSiU8r/3ZrH7iei00R0gojerNl+DREdVh77e1LKRonITkT/qmx/kYjWN/Yt1o5dx0VwgbAcjO732LFrxI/DeeKQUjOVOF0um2o5EAFWsz6Lp1TLoQa3F2/brWfLIVqnW8ljk62CUDyNpXiqpOUAyCLSbpbDWLczL0W3kzg3F8XGfrm7rupWSmVwcSGGLz99qi0Es5Zv2wMA7ijY9kkATzHGtgB4SvkdRLQTwD0AdinHfJWI+HLnawDuA7BF+cef88MAFhljmwF8AcBfL/fNrBQ9t8/glkOvx45doz4cnQohI2WRzbKSK85ulxVLiuXgsJgrtnBoJctxK+k55sCDju4aU1lNigswEEkikc4WWYCc3967Dr9+xUjDznM1GFPmfk8sdpZriTGGs4EINvV7AMjfbRMBsaSEB18ax+eeOKlOotQzVf8iGWO/ALBQsPlOAN9Sfv4WgHdqtn+fMZZkjJ0DcBrAdUQ0DMDHGHueyZL57YJj+HP9O4DbqUV3MpOJYDWTLovg5iO5mddbB71IpLOYCScRTmSQZSi6qXS5bFhULAe9BqOB+txKTqsZFhPp3K3ELYfaXUA+pwXjCzHl59KWw323bMI72k4cXAA6TxwC4SSiKQkb+mTLgYjktt2pDI5NyQWs7TDudrl3jUHG2BQAKP8PKNtHAVzU7DehbBtVfi7cnncMYywDIAggfxSWAhHdR0T7iWh/IBBY5qlXxmY26dNyiCRhIjlNdcivzFMOxrEU52X6hW6lfMtBr3A/emHhVymICF6HRd9upToD0oAclObiUHgd2xmestzOQenZUKJuS5X3leJuJUC2JKPJDI4q4jBjBMuhTkqt+FmF7ZWOKd7I2NcZY3sYY3v6+/uXeYqVsVn0KQ6BSAo9bhvMJsKQj4tDUk2ZK3Yr2RBJZhBJZnRtOdSTygrof6ZDLCnBVGdmkd9pxeRSQv3ZKHS55FkGE4uxVp/Ksvmtb7yIv37seF3HnJuTxYFbDoDcLubiQhwzIdk9bGTLYUZxFUH5f1bZPgFgjWa/MQCTyvaxEtvzjiEiCwA/it1Yq4ZexWE+kkSvW+5LNKxYDlPBOBZ5g68S2UoAMBNK6tpyyLXPqLFozGHVtVuJN92rxzPqc1jVAU2lspXaFSLCqJKx1I5IWYazgSgm67R8zgYisFtMGPHnij3dNgsOXlxSfzeyODwC4F7l53sBPKzZfo+SgbQBcuD5JcX1FCaivUo84QMFx/DnugvA06yFoXybxaTPmEM0pQ7K8TutcFhNmA4myk8PU8RiOpjQteWgNt6r8Rz17laKpTI1B6M5WmvBSJYDIMcdLrWpOMxFksgoCR/1cG4uig197rz26i6bWR2A5HNYjOFWIqIHATwPYBsRTRDRhwF8FsAbiegUgDcqv4MxdgTAQwCOAvgpgI8yxnhy70cA/BPkIPUZAI8p278BoJeITgP4IyiZT61CzzEH3tGUSHYtTYcSWIyWjjlwy2E2nNC15TDkc8BiIrXlRzV8Dn27laLJ2juycowtDs62dSvxWMlyxEEbbwByvbYGfXZsHfS2RbZS1SUOY+y9ZR66vcz+nwbw6RLb9wO4rMT2BIC7q53HamGzmHVZ5zAXyVkOADDkd2A6mCg7WpLfbNMS07Xl8IbtA/jFJ25Dv7f8nGUtfFRoNsvwxadO4d1Xj2Jdr7v6gatENFX7LAcOz1AiytVyGIXRLidCiQxCiXTJGR16ZlIVh9ot1bSUxfhCDG+5fChvu0v5Tuwc9sFltxQ1z9Qj+r1rtAg9upUSaQmRZEa1HAColsNSLA2vI9eRlaN1M+m5YRsRYUTTiK8afODPa5NBfPGpU/jJoakmnl39xJK1T4HjcHHw2i0wt8Gkt3rg6azt6FrKiUOqatHa3/3sJN779RdwJhBBJsuwoc+T9zivmN8x7MOw8rer90I4IQ4F2M0mtVmdXuBzHHrdWsvBiZlQAouxVEmXjHabXtp1NwKfw4poSsLPT8ipzIs6mzamnR9dK7zbbLkCuHZmtI0L4XgGWVpiarygFIwx/Ou+i3j+7Dz+5N9eBYAitxJv271zxIchvwOJdFbX3YUBIQ5F6NFymFdadWsth2G/A2mJ4UwgUjLDxWUzqy0z9Gw51ItXuZH+9Mg0AOhuFGUslal50A+HuwSNFm8AclXSl9ow7qDNUqoUdzg+HcZUMIF1vS51SuPGvnxx4HGoHcM+DPJUdJ1nLAlxKECPqax8joM25sC/YCdnIiVvKkSkrkQNZTko7/XIpPxHuBDTlzhEU1LdMQd+/YyUxsrpddvgsJra03II1iYOz5yQM/m/8zvXY8uAB/1ee5EVuHdjL35txyDW97pzRaw6F4f6vsUdgF1H4sCYbM7yCXCFlgMg94Eql+nT7bIiEE4aynLgLhhAbqehN8shmszU3JGVwwWvXOuMdoaIsKbbhfPz7Wc5TC3J1sCF+RiCsfLi8OzxAC4b9WFtrwvf+73rMRcu/k7euLkPN27uAwC1iFXv6azGWVI2CD25lb767Bns+svH8ZlHjwHItxy4OADlV5xGtBx4No/VTLh1a7/ac6qZvHRuAX/xo9fKBhATaQk/OTQJKcsQS0lqZkqtqJaDAcUBALYOenFqtr2G/iTSEuajKewY8gEobzkEY2kcGF/EbdvkDkIDXgd2jvgqPveAT17k6d1yMM5do0Hoqc7hTCACt82Cy0b8uH37QJ4vu9djVzNbyt1U+HZDWQ5KZ9ar13ZjrNupVog3k6eOzeA7L1zAxYXSrpH/ODSFj/3LK/iPw3LmlKfOOgefQ57cZ0S3EiCLw/hCrK1mKPN4w47hyuLwi1MBeQjXtoGSj5fCbjGj123Tfa2DcCsVoKeYQzCWxtoeF777u9cXPWY2EQa9dkwGE2WzXLoNaDnw93rT5j6YTIRYSpKbCzZxEE5Y6bS6/8IC1va6ih4/Py/30vnWc+cBoO6AtM1iwhfvuQrXrOuuvnMbsm3IC8aAUzMRXLGmq9WnUxM8U2nHsBdAeXF46tgMul1WXFnn+xr0OTAjLIf2Qk/isBhLodtdfjU5qLiWyrqV3O03e7gao11OfOE3r8C9N65XU3vnmxx3iCS4OCyWfJx3VD2gPF5vKisAvOOKEbWLqdHYPiTfYNtpnjQPRm8b8qrDmAqZjyTx6GvTeNvu4brrU3gRq54xzl2jQdgsJiR1EnOQZ0OXz30friIOOcvBOG4lAHjXVWPwOazoUcSh2bUOEcVyeLmMOFyYj2FEEwOqtwjO6KzpccFhNeHETBuJw1IcRMCw3wmfw1rScvj+votIZbK494b1dT+/sBzaELsSc9BD9eJSLF0xSDnkk1ea5QSEH2s0ceD0rLLlcGImXPImMb4Qw63b+nHFmB/A8iwHI2M2EbYOetvLcliKY8Brh81igt9ZLA5pKYvvvnABN23uw5ZBb93PP+RzYD6aQjKjr4JbLUIcCuBDZ9JSa8Uhm2VYKlP9zBnyy1kPZQPSyrFGcitp4eKwoIxQbRbhZAZehwWMAa+M51sP4UQaC9EU1va4cfceuVu936CB5ZWwddDbZpZDAsNKy+1S4vDEkRlMBRP44OvWL+v51/TIz81nP+gRY941VgAXh1ans4aTfPxn+RvNG3cO4Z5r12BNT3GQFMilvjoN6ubIiUNzu7RGkmncsLEXZhOpcQUOjzes7XHhnmvX4J8/dC12DldOZexEtg95EQgndVeXUo7JYFyNAZUShx+8PIGxbidu2157lpKWvRvlYZe/PDm3shNtIkIcCrApDexa3V9pqcwQHy0b+tz47Ht2w2oufRmvXtuNz7zrcrxuU19TzrHV+BxWmE3UdMshkshg0OfAjmEv9p/PF4eLijis63XBYjbhtm0DdQ366RS2Kq6X49OhFp9JdRhjmFpKqJXMpcThwnwUl4/6l90ocaTLic0DHvziVHPGHTcCIQ4F2JSagFZbDur4zxUURplNhPddv7am+cztiEmZA9FsyyGalOBxWLBnXQ8OXlxS25kAcjAaQFnrTSDDM5ZOtkHcIZzMIJ6W1Epmn9Oa17abMYapYE48lsstW/rx0rkFJHTW6JNjzLvGClDdSi1OZ+XFXZVSWQVAj9vaVMshmZGQkrLw2C2457o1YGD4yHcPqN+P8YUYulxWQzbNayT9Xju6Xda2iDvwthaDGsshFE+rSSrhZAaxlJTXpWA53Ly1D8lMFi+ea9lU5IoIcSjArhNxyI3/NF4b50bS47Y11Y/NM5U8dgu2D/nwN3ddgX3nF/GXjxwBIIvDWmE1VIWIsLbX3dIGfPGUhHd+5T/xqmaWcyl4W4tBZQCV32lFSsqqbbt5fcKQf2V1KXs39MJmMeGXJ/XpWhLiUAC3HFo9Da7c+E9BPr1ue3PFIZkTBwD49StG8Pu3bsKDL43jpXMLuDAvxKFWfA6L+nm2gslgHAcvLuH5s/MV95sJyZaoNuYA5BZsXBxWajk4bWZct75Ht3EHIQ4FDCirhe+9ON7SWody4z8F+XS7rU0VhzC3HDTdYP/g9i3o89jx+SdO4NJSXIhDjXjsFtUSawXxVP7Kvxy8OI23xS8nDjwmsRJu3tKHkzMRzIb1VxAnxKGAq9Z24yOvl1eGn3viRMvOo9z4T0E+PW47luJpSNnmCDlf6Xo1hW1Omxm/f+tGvHhuAVKWYV2JfkuCYjz21loOUeW1p4KVXVvTwQT8TqtaPKqKg5Ikwhvm8e6qK2GDMhRoNtTcjLvlIO48JfjEm7fhnmvX4CvPnMHhidYMAi83/lOQT4/LCsZyqb+NJlLCcgCA91+/Tp2vITKVasPjaK3lEFNiBtW6oc6EEnlWQZHlEIqjz2NrSLfjwufWE0IcSkBE+IDSL2WiReMNl2JpdItK26r0KDfoZrmWCmMOHKfNjI+/YTOsZsLmAU+pQwUFeO0WRFIZZJtk5VUjlqxdHLRWQeENvBFprBw+4KlUY79WI8ShDGr1bRPnBXzu8RP4/398tORjS7EU/MJyqEqzO7OWEwcA+MAN6/D8/bdjwNuYG4XR8SgtSGItyuvn8yTmIsmK2YgzoWRlyyGYUPuarRRVHBJCHNoG3rZiqcJ4wJXy7MlZ/PjQZMnHluLCcqgF7nprVmdWVRwcxeJARHmjWwWV8djl73OrXEsxJSDNGMoGgKUsQyCSzLMMvA5LXtvu6VBixZlKHOFWakMcVjOcVnNT20HPhVMIhJOYjxQHoxajKZHGWgO8f1TTLIdEBiaS51ULVoZXEdhIsjU3wqhmEh3POHr+zHyeUMxFkpCyDAMay8FkInjtFgTjacRTEpZi6Ya5ldw2M0wEhOL6m5InxKEC3S4rFptkOWSzTG3DUNjKWMoyhBIZUQBXA6thOXjsFtEvqQFw6yvUIsuBp7ICwGQwgUA4iff90wu488v/iVNK5TZPYy1MU+1y2RCMp9UCuUaksQKy9elzWoVbqd3octmalgUTjKeRUQJzxwvEgZuYwq1UHZvFBKfV3DSzPJzIwOsQ16ER8HTgVrmVoklJbZQ3HYzj8KUlMCZnBt71D8/jtUtB1aIYLEhT3dDnxr7zi+ps6Ua5lYDSjf30gBCHCvS4bU0LSAc0rqRCy6GWjqyCHP4mrrwiybQY3tMgPKpbqbHiwBjD+RrmIsTTGfS4bXDbzJgKJnBoIggi4OGP3gSLifDFp06VtRzec80YLi3F8YMDE/LjDRQHn8MqspXajS6XtWkB6bmwLA4OqwnHC5qRcVdWpVkOghzNXHlFkpmSwWhB/XiaZDk8c2IWt33+WZwNRCruF01KcNvM6vzmwxNBbO73YNuQF+++ehTPnpjFsekwzCZCb0GiwZt2DsLvtOJHBy8BaLA4OC0tc7VVQohDBbpdNrU7aqPhlsO163twaiacl/sdjAvLoR58TksTxUGCW1gODcGrZCuFG2w5nA1EwRhwbKpyx9dYSoLTZsFIlxOTwQQOXQricmW067uvHkNaYvjhyxPo99iL5jQ4rGa866pRZJm8GHHZGvedEG6lNqTbJV+0ZrRmCCiWw81b+hBLSbioKbZbjIqYQz3ILZWbs/KKJNJ5rTMEy8dtlzO+Gm05zCp/S9Ush1gqI1sOPgdOTIcQCCexe1QWhx3DPuwY9iGRzqqtugv5zWvlMbCNjDcAwq3UlnS5bGCsOdWLc5EUbGYTrl3fAyB/1bMoYg514XM02a0kxKEhWMxy8kCjU1l5nKDaPGbZcjBj2O9AIi0XwV0+1qU+/u6rRgHkWnUXsmPYh70be9TBRY3CJyyH9oMP2mmGaykQTqLPY8O2IS+I8oPSwXgaJoJYsdaIz9m8lVckIWIOjcTrsKidbhsFzzA6oxGHaAnXlWw5WDCszIY2mwi7RnLzvu+8cgQmkkd4luNbv3MdPnf3FY06dQCy5ZvMZHU3EU6IQwXUHPomiMNcJIk+rx0umwVre1w4MZObrTsfTaHLZYNpmfNpOw2/04pwMtNw95+UZYimJGE5NBCPw9LwmIPWrcQYw+GJIHb/f0/gtUv5TTNjKQkuJSANyHOtHZrixgGfA9/84LX4vVs2ln0tu8Xc8E7JPmXx0WjRXClCHCqQK7Bq/Ko0EE6iX8mI2DXiw8HxJXV+xIHzi9g57Kt0uEAD708TbnA6K6+o9QrLoWF4GzzTgTGGmVACDqsJ4UQGc5EUnj4+CynLcODCYt6+sZQEl92sxgx4vEHL67cNYLSC5dAMfDptoSHEoQJNtxwUcbhlSz8mgwmcmAljJiT/f/OWvoa/plHxq50tG7vy0o4IFTQGT4OnwUWUec48dnc2EMGL5+RJb8enQ3n7xlIZ1VLvcllx67b+hp3HStBr8z3xra9Al7s5zfeyWYb5aAr9SuDrtu0DAICnj8+qHT5v3qKPL2470KzmZfwmJlJZG4fHbsFcuHFt8Hkweu/GXvzy1BxOzITx8rhsMWg7D0hZhkQ6C5fNDJfNgpf//I26cdv6HPps270iy4GIzhPRYSI6SET7lW09RPQkEZ1S/u/W7H8/EZ0mohNE9GbN9muU5zlNRH9POmlk47VbYDFRwy2HxVgKUpahT2kaN+hz4PJRP54+Notfngqgz2NveEaEkeE+22aJgwhINw6P3dpQy4HPe75qbRfsFhN+9MolJNJZjHY5cWI6Vz8UV4K9LpscY9CLMACA39mc7+9KaYRb6TbG2JWMsT3K758E8BRjbAuAp5TfQUQ7AdwDYBeAOwB8lYh4NOhrAO4DsEX5d0cDzmvFEBG6mlAIxwvg+jVzAN6wfQAvjy/imeOzuGlzr66+vHrH72qOWc7dSiJrrHHI2UqNu07cchj2O7Ghz42Xx5cAAO/fuzavfiimCFIji9caRc6tZPyA9J0AvqX8/C0A79Rs/z5jLMkYOwfgNIDriGgYgI8x9jyTI7Lf1hzTcrpd1oYHpOfCsthwywGQxSHL5C+IcCnVR7PdSsJyaBx8jjRPvlgp3HIY8NqxsV+ex7xt0IvXbZJjdty1xGc5cMtBTxjSrQSAAXiCiA4Q0X3KtkHG2BQAKP8PKNtHAVzUHDuhbBtVfi7cXgQR3UdE+4lofyAQWOGp10YjWmjsP7+QV6ATiMirnX5Nsc3lo341QC2C0fXB/7gaLg4iIN1wvA4Lsix3s14pM6EEvHYL3HYLNvbJ41qv39iDrYMeEAHHleJSnnmmR8vBYTXDZjEZThxuZIxdDeAtAD5KRLdU2LeUn4RV2F68kbGvM8b2MMb29Pevzuq6Ec33/uD7B/G3T55Uf1ctB404mEyEu64Zw02b+/IGjQiq47KZYTFRw/+4eD4+7wkkWDmN7syqnffMLYfrN/TCZbNgnaZ+KK5jywFobmfh5bIiGWWMTSr/zxLR/wFwHYAZIhpmjE0pLqNZZfcJAGs0h48BmFS2j5XYrgu6XTYcvLi07OMzUhZTwThmQrnc6blIEjaLqciX/cm3bF/263QyRNSU5mXccuA9gQQrh1th4UQGgw0o5ZkJJdSittu3D+L3bt6ANyjZf9uGvBrLQRYHvV5Ln8Oiu2lwy7YciMhNRF7+M4A3AXgNwCMA7lV2uxfAw8rPjwC4h4jsRLQBcuD5JcX1FCaivUqW0gc0x7ScLrdsOSzXRxqIJJFluRbdQK4ATidJWYagGf1pwok0HFZTwytiOxlvwy2HJAaVxA6/y4o/e9tOOBXrYPuQD+fno4inJMQVt5LTqj+3EqDP/kor+aQGAfwf5QZnAfAvjLGfEtE+AA8R0YcBjAO4GwAYY0eI6CEARwFkAHyUMcYdjx8B8AAAJ4DHlH+6oMdlQ0rKLruNwuSSHF/QDvcJKK0zBI1DHrXY2JXXZDCOYf/qVssaHY/iomtElTRjDLPhRFk37I5hL7IMODkTRjSpb8vB77RioYnz6pfDssWBMXYWQFEHKsbYPIDbyxzzaQCfLrF9P4DLlnsuzUQ7o3g54sCbgoUTGSTSEhxWM2ZDSaztdTX0PDsdn6PxMx0uLsSxpkdcp0aiDvxpQGfWhWgKaYkVjfTkbOqXA9Tn56OIKXUOTp3GHHwOa03T7FYTYS9XgU9jW25QeioYV3/mMxwmg3GMNLgnfKfjb0Jn1ouLMazpFpZDI/E2qMmc3FNJ/nsaLGM58O6ql5biap2DW4fZSkBzB1YtF31+Ujqi2y1bDvPRZJU9S8MtB0AORHe7bQgnMhgS7oqG0mhxCCXSWIqlsVZYDg2lEeLw09em8QfffwXblC4C5cTBbbegy2XF5FIcvW7ZunBa9Wk5+BW3KGNMN7FIYWT0V8gAABR8SURBVDlUYUOfnB53aqbylKlyTIVy4hAIJzGtWBIjXcJyaCQ8oNeo4qqLC3JlrXArNRa3feUB6Z+fnAURMLkUh9VMFQV8tMuJyaUEYqkMHFaTbjsP+BxWtUW8XhCWQxX6PHaMdjnx6sTy0lmngwls7HPj7FwUc5GU6vMcErUMDcXvtCKTZYilGjPzmYuDsBwai9VsgsNqWpE4HJoI4tr1PfjnD16LhVgqr5i0kJEuJ8bnYxj2O3TrUgI0LTTiad0UXQrLoQauWOPHoQl5cMhcJImrP/UkHjs8VdOx08EEdil94wPhJKaU7KVK06YE9ePXtD2+uBDDrMZiq5WfHJrE3/1MLla8uCBbeGu6hTg0Go/dumy3UiIt4cR0GLvH/LCYTWoX43LIlkNcHRGqV5rVAmYlCHGogd1jXRhfiGExmsLTx2axEE3hb588qXZ8LIeUZZgOJbC2x4lulxWBSAKTilupnJ9UsDx4C43JpQTe+sVf4rrPPIU7v/wrPHdmrubn+O4LF/CVZ04jnpIwvhCD12FRm/oJGod3BTMdjk2FkMkyXD7aVX1nyOIQTmYwE0ro2nLg1s9seHmxzWYgxKEGdo/JK/9Dl4J45sQsTAScmo3gqeOzFY+biyQhZRmG/E70eeyYC6cwHUygz2OHzSI++kbCV17ff2kc4WQG996wDmcCUTy072KVI2UYYzgxHUZaYnjl4iIuLsaES6lJeOyWZScPcAv+ijXFU9xKwS30U7MRXVsOfDrd1FK8yp6rh7hD1cDlo34QAQcuLOKXp+bwnqvHsKbHia8+e7piAHRKyVQa9jnQ77UjEEliMpgQwegmwMXh4YOT2Njvxl+9Yxc29buxWGMK8lwkpe774tkFjC/EhEupSewY9uL5M/M4E6g/yePQRBB9HnvNMTv+txYIJ3VbAAfIngQiYDJYvzu0WQhxqAGvw4qNfW78y4sXEElm8Madg7jvlk14ZXwJL51bKHscz0wa8ivioGQriWB04/EpA1NSUhZ3XTNW9yyOkzNyDx6rmfDC2XlMLMZFoWKT+JM3b4PDasL9Pzxc1TVbyOFLS9g95q853XNUU6ei19YZgByo7/fYheXQjlwx1oW5SAo2swk3bu7DXVePwUTAc2fmyx7DLYeRLsWtFJED0iIY3Xi45WAi4D1Xy30ce9y1i8MJpe//Wy4bxr7zC0hlsqIArkkMeB3472/dgZfOLeDfDtTm9gOAaDKD07MR1c1bC31uO2xKbyw9Ww4AMNzlxPQyEimahRCHGuFfyOs39sBtt8BpM2Oky4nz8+VL3qeDCdgsJnS7rOj32hFLSQgnM2oXSUHj8DqsIAJu3dqvBvu7XbaaBzWdnAmjx23DWy8fBl/MihqH5vEbe9bgijE//vk/z9d8zJHJELIMdYmDyUSqa0mv7bo5I34HJnVkOejXztIZV66VR2Hftm1A3bahz12xH8pUMIFhvwNEpA7yAXLBJ0HjMJsI/+Odl+G69T3qtm6XPK84lclWTQA4MRPG1kEPrtuQO16IQ/MwmQg3bOrDN351FmkpC2uFzrfn56L4+6dO4Ren5MyzWjOVOPIiLqbLQT9ahvwO/PxkQDdV0sJyqJErxvz4yvuuxvuuX6tuW9/rxrm5aNmg9HQwocYXtIU6otNnc3j/9euwZdCr/s5bnyxVcS0xxnByOoztQz70uG3qFLFR4f5rKtuGPEhLrGrDuQdfGsePDl7CDZt68ZX3XV2x6K0U3I2rf8vBiVhK0s0saSEONUJEeNvuYTg0vVnW97kRSmTKZsRcWoqrLqR+YTmsOj2KOCxUEYeJxTiiKQlbFWF5084h7Bjy5V1rQePhn/fJKq1ppkMJjHY78aX3XoW37R6u+3VGVXHQt+UwrLi/tM06W4kQhxWwoU92O5ybK/5yX1yI4dJSHJcr1dF9XvlGRSQK4FYL3lG3Wp98nqm0bUhu8fxHb9yKn3z8puaenACb+j0wkezSq8R0MKEO9FkOo21iOeRqHfQRlBbisALW98pN+c7NxYoee/ZkAADweiVG0eu2w0QQBXCrSI/qVqoclOY3J+6SMplItw3ajITDasb6PjdOTlcWh5lQAoMrsLbbxa3E3c1TOql1EHepFbCmxwWziUr6TJ89Pos1PU5sUoaem02EHrdduJRWkR5lUFMlyyGSzODJozMY8TvUFhyC1WPboFe13ErB5zaspDZo25AXvW5bXjxKjwx45QWkcCsZAKvZhLFuJ84VpLMm0hKeOzOP27YN5GUdbB5wY5vOv6BGokszxa8UFxdieM9Xn8OhiSD+8I1bV/PUBApbB704Px9FIl26VXUokUE8La1IHPq9dhz4izfiyjX1ZTmtNryR4KRO3Er6jtC0ATydNS1l8dHvvYybt/RhXa8b8bSUl/YKAN/84LUw6SBFrVOwWUzw2C0lEwYYY/jYg69gKhjHAx+6Fjdv6W/BGQq2Dclznk/PRnDZaHH9woxSFLYSt1I7MdzlwHRIH5aDEIcVsr7XjZfOLeDhg5N44ugMnjg6g039btgtJuzd2Ju3r96zJYxIt9taskr65fElvHpxCZ9652VCGFrI1kE5CeDkTLiyONSZvtqujPidODYVavVpABBupRWzoc+NWErC3zx+HDuGfbhj1xDOBKK4YVOvrrtAdgrdZforPfDceXgdFrz7qtEWnJWAs67XDZvZVDZjiY/Z7ZSuAsN+ByaD8YZNNFwJYim7QtYrY0RnQkn8v2/fhTftGsSXnjqF128fqHKkYDWQW2jki8N0MIHHDk/hQzeub8jUOMHysZpN2NhfPmNJtRw6JP17yO9AIp3FUiytFnG2CvGXsUI2KOmsG/vduOOyIZhNhD9607YWn5WA0+O24WxBHcp3X7iALGP4wA3rW3NSgjy2D3nx/NnSDSynQwl0uawdU5DI024ng/GWi4NwK62Q0W4nXr+tH3/+th0wi9x43dHlsmKpoPnek0dncOPmPtE7SSdcPtaFmVCy5GjXmVByRQVw7QavnToTqNxSZDUQ4rBCzCbCAx+6Dm/YPtjqUxGUoMdlQ1hpvgcAGSmLs3MR7BqpvbOnoLlcwSctKlPetKy0AK7d2DzggdVMODJZ/FmsNkIcBIamsPnehYUY0hLDlgFPK09LoGHniA8mksfwFiI3r+yMTCVATr/eOujF0cnWZywJcRAYmm5eCKfUOpxSW2UIcdALLpsFWwa8ODSxlLc9I2UxF1lZdXQ7smvEh6OToZZnLAlxEBiabnd+871TSgfQzcJy0BW7x/w4PBHMuyHORVLIMmCgw8Rh57AP89EUZkLJlp6HEAeBoekpcCudmo1grNspChJ1xu4xP+ajKVxaiuPhg5dw/w8P4ZIyFa3jLAelGPDoVGvjDuIvRGBouFtpQSMOIt6gPy4fk/se/fS1aXz+iZOIpyUcV2ofOqUAjrNj2AcAOHIp1NJEF2E5CAwNn+mwGE1ByjKcCUR0352zE9kx7IXVTPjsY8cBALdt68cr43IMolMK4DgeuwXre1040uKgtBAHgaGxW8xq873xhRhSmaywHHSI3WLGtiEvMlmGP3zjFnzxvVdhrNsJq5nQ2+JisFawa8SPoy3usSTEQWB4ej02HJ8OqXMDhOWgT+7YNYTrN/TgQzdugM9hxQMfuhafu/uKjhy8tHPEh/GFGEIJOcsuLWXxxZ+dwmx49dp5i5iDwPC8//q1+Myjx9WJcCJTSZ987A1b8LE3bFF/3zzgxeaBzhTynSNy3OHg+BJu2dqPn742jS/87CRSkoT/9ubtq3IOwnIQGJ4P37QRezf24MhkCKNdTnhEsz2Bztm7oRfdLiu+88IFAMCDL40DAB49PL1q9Q9CHASGx2wifP43roTXYcH2oc5ciQraC6fNjN/euw4/OzaDp4/P4Lkz89jU78a5uSiOTVWeud0ohDgIOoLRLid+9NEb8al3XtbqUxEIauK3b1gPq9mEj//LKzCbCF9+39UwEfDo4alVeX0hDoKOYVO/R22JLBDonX6vHe+5ehTRlITbtw9gx7APezf24tHDU6viWtKNOBDRHUR0gohOE9EnW30+AoFA0Gp+7+aN6HXb8OGbNgAA3nr5MM7ORctOzmskuhAHIjID+AqAtwDYCeC9RLSztWclEAgErWVjvwcH/uKNuF6ZR//mXUMwEfCfp0sPR2okeknbuA7AacbYWQAgou8DuBPA0ZaelUAgEOiIfq8dv/jEbRjrbv6gKl1YDgBGAVzU/D6hbMuDiO4jov1EtD8QCKzayQkEAoFeWA1hAPQjDqVKIIsiLoyxrzPG9jDG9vT396/CaQkEAkFnohdxmACwRvP7GIDJFp2LQCAQdDx6EYd9ALYQ0QYisgG4B8AjLT4ngUAg6Fh0EZBmjGWI6GMAHgdgBvBNxtiRFp+WQCAQdCy6EAcAYIw9CuDRVp+HQCAQCPTjVhIIBAKBjhDiIBAIBIIiaLXavzYaIgoAuFDmYT+AStO5qz3eyOeq9bUa9Tx9AOaa/Fqr+Z5W+/Nbrc8GqH6t9PKeVvO1GvU6jXwtvVynRjxXHwA3Y6x6LQBjzHD/AHx9JY838rlqfa0GPs9+I72nFnx+q/LZ1HKt9PKe9Pj5NeKcG/U3paf3VG2/Wu4P/J9R3Uo/XuHjjXyuWl+rkefc7Ndazfe02p/fan02taCX97Sar9Wo12nka63W66zmdahK27qVBKUhov2MsT2tPg9BdcS1ag+MdJ3qeS9GtRw6ma+3+gQENSOuVXtgpOtU83sRloNAIBAIihCWg0AgEAiK6AhxIKJIlcefJSJD+BTbGXGd2gtxvYxNR4iDEan2hynQB+I6tQ/iWuXTMeJARK8nop9ofv8yEX2whackKIG4Tu2FuF7GpWPEwYgQkYeIniKil4noMBHdqWxfT0THiOgfiegIET1BRM5Wn2+nIq5T+yCuVQ4hDu1NAsC7GGNXA7gNwOeJiE/V2wLgK4yxXQCWALynRecoENepnRDXSkE3LbtXgQzyxdDRqhNpIATgM0R0C4As5Lnbg8pj5xhjB5WfDwBYv/qntyzEdWovjHa9jHyt6qKTLIcLAHYSkZ2I/ABub/UJNYD3A+gHcA1j7EoAM8j9cSY1+0lon4WAuE7thdGul5GvVV0Y+s0BABFZACQZYxeJ6CEAhwCcAvBKa8+sIfgBzDLG0kR0G4B1rT6h5SKuU3th4OtluGu1XAwvDgB2ATgDAIyxTwD4ROEOjLHXr/I5rQj+hwngewB+TET7ARwEcLylJ7YyxHVqLwx1vQx+rZaFodtnENHvA/h/APxXxtgTrT6fRkFEVwD4R8bYda0+l0YgrlN7YcTrZdRrtRIMLQ5GxIh/mEZEXKf2QVyr0ghxEAgEAkERnZStJBAIBIIaEeKgc4hoDRE9o1RnHiGiP1C29xDRk0R0Svm/W9neq+wfIaIvFzzXe/9ve/cWYlUVx3H8+yeNMjOjVAoTCaQ0K8WH7EZUKPUQBRUR0Yz20M0uDwVlFEHQQzcpfcgiLbUoDYusyCghyW4PJjRNkmhIWUN2mXRSCINfD2sdOLhntDPuc+bMOb8PDGfOOnsv1uI/s/9nr73X2nnW5zcRsT4iTh6KPrWqkmN1Q45Td0Q8ORT9aVWDiNOciNic/3c2R8RlVXXNyuXbI2Jx1YS5Yc/Jofn9C9wnaSowG1gQEdOAB4ENkqYAG/J7SDM8HwHur64k343xHHCppHNItx7e1ZgutI2yYnUS8BRweZ6NOyEihvv8gWZSa5x+B66SdDbQCayqqut54FbS7OkpwBWN6UL9OTk0OUk9kr7Ov/cBW0mzNq8GVuTNVgDX5G32SdpEOvBUi/xzXP52Mwb4pf49aB8lxup0YJuk3/L7j2nxpRoaaRBx2iKp8r/SDRyTJ/2dAoyR9IXSxduVlX1agZPDMBIRk4GZwFfABEk9kP7YgfGH2lfSAeAOoIuUFKYBy+rY3LZ2JLECtgNn5sXeRpAOOKfVr7XtaxBxuhbYIukfUkLZVfXZrlzWEpwchomIGA2sJd1ut3cQ+48kJYeZwKmkYaWFpTbSgCOPlaReUqxWA58CO0lDIVaiWuMUEWcBTwC3VYr62axlbv90chgG8oF9LfCapLdy8a/5tJb8uvsw1cwAkLQjnwKvAS6oU5PbVkmxQtK7ks6TdD7wPWlpCitJrXGKiInA20CHpB25eBcwsaraibTQUK2TQ5PL1weWAVslLar6aB3p4hj59Z3DVPUzaYG0cfn9HNJYq5WkxFgREePz64nAncBL5ba2fdUap4gYC7wPLJT0WWXjPPTUFxGzc50d/I/YDheeBNfkIuIi0tBCF2kJYYCHSGOka4BJwI/A9ZL+zPvsJF1wPpq07vxcSd/lmaD3AgdIq2nOk/RH43rT2kqO1evAubmOxyS90ah+tLpa4xQRD5OGYKvP3uZK2h3pGdmvAMcCHwB3q0UOqk4OZmZW4GElMzMrcHIwM7MCJwczMytwcjAzswInBzMzK3ByMKuDiLg9Ijpq2H5yRHxbzzaZ1aIdniFt1lARMULS0qFuh9mRcHIw60dekG09aWLUTGAbaQbsVGARMJq0lPM8ST0R8QnwOXAhsC4ijgf+lvR0RMwAlgKjgB3ALZJ6I2IWsBzYD2xqXO/MDs/DSmYDOwN4MT//Yi+wAFgCXCepcmB/vGr7sZIukfTMQfWsBB7I9XQBj+byl4F78vpJZk3FZw5mA/upai2dV0lLLEwHPsoP/DoK6KnafvXBFUTECaSksTEXrQDe7Kd8FXBl+V0wGxwnB7OBHby2TB/QfYhv+vtqqDv6qd+saXhYyWxgkyKikghuBL4ExlXKImJkXuN/QJL2AL0RcXEuuhnYKOkvYE9eBA7gpvKbbzZ4PnMwG9hWoDMiXiCtyLkE+BBYnIeFRgDPkh4deSidwNKIGAX8AMzP5fOB5RGxP9dr1jS8KqtZP/LdSu9Jmj7ETTEbEh5WMjOzAp85mJlZgc8czMyswMnBzMwKnBzMzKzAycHMzAqcHMzMrOA/aviBCcFgkToAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-150:].plot()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Etude de l'incidence annuelle" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Comme l'incidence de syndrome grippal est très faible en été, cette modification ne risque pas de fausser nos conclusions.\n", + "\n", + "Encore un petit détail: les données commencent an decembre 1990, ce qui rend la première année incomplète. Nous commençons donc l'analyse en 1991." + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [], + "source": [ + "first_september_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", + " for y in range(1991, sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "En partant de cette liste des semaines qui contiennent un 1er août, nous obtenons nos intervalles d'environ un an comme les périodes entre deux semaines adjacentes dans cette liste.Nous calculons les sommes des incidences hebdomadaires pour toutes ces périodes.\n", + "\n", + "Nous vérifions également que ces périodes contiennent entre 51 et 52 semaines, pour nous protéger contre des éventuelles erreurs dans notre code." + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "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": "markdown", + "metadata": {}, + "source": [ + " Voici les incidences annuelles." + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 20, + "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": "markdown", + "metadata": {}, + "source": [ + "Une liste triée permet de plus facilement répérer les valeurs les plus élevées (à la fin)." + ] + }, + { + "cell_type": "code", + "execution_count": 22, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2002 516689\n", + "2018 542312\n", + "2017 551041\n", + "1996 564901\n", + "2019 584066\n", + "2015 604382\n", + "2000 617597\n", + "2001 619041\n", + "2012 624573\n", + "2005 628464\n", + "2006 632833\n", + "2011 642368\n", + "1993 643387\n", + "1995 652478\n", + "1994 661409\n", + "1998 677775\n", + "1997 683434\n", + "2014 685769\n", + "2013 698332\n", + "2007 717352\n", + "2008 749478\n", + "1999 756456\n", + "2003 758363\n", + "2004 777388\n", + "2016 782114\n", + "2010 829911\n", + "1992 832939\n", + "2009 842373\n", + "dtype: int64" + ] + }, + "execution_count": 22, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Un histogramme pour montre la fréquence de l'épidemie" + ] + }, + { + "cell_type": "code", + "execution_count": 23, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 23, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFGtJREFUeJzt3XuQJWV9h/Hnxy4oMLggDKgLukbIRGFVZL2gVTqDxhKXaGm8IZpotNYqFTFlyqCoaLytF7zESyobJZB4nShUkDUaDI6oUZQV4oI4amBVlighKjCI4Oovf3QvGadmdmZO95lzOu/zqdraPrfu73mn53v69Ok+E5mJJOn/v70GHUCStDIsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhVq/kwg488MA88sgjV3KRrbn11lvZf//9Bx1j2bqaG8w+KF3N3tXcsHj2bdu23ZiZo02Xs6KFf9hhh3HZZZet5CJbMzU1xfj4+KBjLFtXc4PZB6Wr2buaGxbPHhE/bGM57tKRpEJY+JJUCAtfkgph4UtSISx8SSpEo8KPiLGIuGLWv5sj4uVthZMktafRYZmZOQ08GCAiVgE7gfNbyCVJalmbu3QeC/xnZrZyvKgkqV3R1t+0jYizgW9l5vvnXL8J2AQwOjp63OTkZCvLW2kzMzOMjIwMOsay9Zp7+86b+pBmcevXrrlzuqtjDmYfhK7mhsWzT0xMbMvMDU2X00rhR8Q+wPXA0Zn504XuNzY2ltPT042XNwhdPYuv19zrTt/afpgl2LF5453TXR1zMPsgdDU3LOlM21YKv61dOidSbd0vWPaSpMFqq/BPBj7e0rwkSX3QuPAjYj/gD4HzmseRJPVL42/LzMxfAge3kEWS1EeeaStJhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUiMaFHxEHRsSnIuK7EXF1RBzfRjBJUrtWtzCP9wKfy8ynRcQ+wH4tzFOS1LJGhR8RdwMeDTwPIDPvAO5oHkuS1LbIzN4fHPFgYAvwHeBBwDbgtMy8ddZ9NgGbAEZHR4+bnJxsFHhQZmZmGBkZGXSMZes19/adN/UhzeLWr11z53RXxxzMPghdzQ2LZ5+YmNiWmRuaLqdp4W8Avg48KjMvjYj3Ajdn5mvnu//Y2FhOT0/3vLxBmpqaYnx8fNAxlq3X3OtO39p+mCXYsXnjndNdHXMw+yB0NTcsnj0iWin8ph/aXgdcl5mX1pc/BTyk4TwlSX3QqPAz8yfAjyNirL7qsVS7dyRJQ6aNo3ROBT5aH6FzDfD8FuYpSWpZ48LPzCuAxvuWJEn95Zm2klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqROM/Yh4RO4BbgN8AuzLTP2guSUOoceHXJjLzxpbmJUnqA3fpSFIhIjObzSDiWuDnQAJ/m5lb5ty+CdgEMDo6etzk5GSj5Q3KzMwMIyMjg46xbL3m3r7zpj6kWdz6tWvunO7qmIPZB6GruWHx7BMTE9va2F3eRuHfKzOvj4hDgYuAUzPzkvnuOzY2ltPT042WNyhTU1OMj48POsay9Zp73elb2w+zBDs2b7xzuqtjDmYfhK7mhsWzR0Qrhd94l05mXl//fwNwPvCwpvOUJLWvUeFHxP4RccDuaeDxwJVtBJMktavpUTqHAedHxO55fSwzP9c4lSSpdY0KPzOvAR7UUhZJUh95WKYkFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSIRoXfkSsiojLI+LCNgJJkvqjjS3804CrW5iPJKmPGhV+RBwObAQ+1E4cSVK/RGb2/uCITwFvBQ4A/iIzT5rnPpuATQCjo6PHTU5O9ry8QZqZmWFkZGTQMZat19zbd97UhzSLW792zZ3TXR1zMPsgdDU3LJ59YmJiW2ZuaLqc1b0+MCJOAm7IzG0RMb7Q/TJzC7AFYGxsLMfHF7zrUJuamqKL2XvN/bzTt7YfZgl2nDJ+53RXxxzMPghdzQ0rl73JLp1HAU+KiB3AJ4ATIuIjraSSJLWu58LPzFdl5uGZuQ54FnBxZj6ntWSSpFZ5HL4kFaLnffizZeYUMNXGvCRJ/eEWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCtGo8CPirhHxjYj4j4i4KiLe0FYwSVK7Vjd8/O3ACZk5ExF7A1+JiH/JzK+3kE2S1KJGhZ+ZCczUF/eu/2XTUJKk9jXehx8RqyLiCuAG4KLMvLR5LElS26LaSG9hRhEHAucDp2bmlbOu3wRsAhgdHT1ucnKyleWttJmZGUZGRgYdY9l6zb195019SLM8h+0LP71t0Cl6s5zs69eu6W+YPZjv57wS496P57yUdX2Q6/WenvNi2ScmJrZl5oamGVorfICIOBO4NTPfOd/tY2NjOT093dryVtLU1BTj4+ODjrFsveZed/rW9sMs0yvW7+Ks7U0/ZhqM5WTfsXljn9MsbL6f80qMez+e81LW9UGu13t6zotlj4hWCr/pUTqj9ZY9EbEv8Djgu01DSZLa1/Rl/J7AuRGxiurFYzIzL2weS5LUtqZH6XwbOLalLJKkPvJMW0kqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCNCr8iDgiIr4YEVdHxFURcVpbwSRJ7Vrd8PG7gFdk5rci4gBgW0RclJnfaSGbJKlFjbbwM/O/MvNb9fQtwNXA2jaCSZLaFZnZzowi1gGXAMdk5s2zrt8EbAIYHR09bnJysqf5b995U/OQDRy2L/z0toFG6ElXc0M52devXdPfMHsw3+9VV8d92HPv6ec8MzPDyMjIgrdPTExsy8wNTTO0UvgRMQJ8CXhzZp630P3GxsZyenq6p2WsO31rj+na8Yr1uzhre9M9YCuvq7mhnOw7Nm/sc5qFzfd71dVxH/bce/o5T01NMT4+vuDtEdFK4Tc+Sici9gY+DXx0T2UvSRqspkfpBPBh4OrMfFc7kSRJ/dB0C/9RwHOBEyLiivrfE1vIJUlqWaMdXpn5FSBayiJJ6iPPtJWkQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqRKPCj4izI+KGiLiyrUCSpP5ouoV/DvCEFnJIkvqsUeFn5iXAz1rKIknqo8jMZjOIWAdcmJnHLHD7JmATwOjo6HGTk5M9LWf7zpt6TNiOw/aFn9420Ag96WpuMPugdDX7sOdev3bNgrfNzMwwMjKy4O0TExPbMnND0wx9L/zZxsbGcnp6uqflrDt9a0+Pa8sr1u/irO2rB5qhF13NDWYflK5mH/bcOzZvXPC2qakpxsfHF7w9IlopfI/SkaRCWPiSVIimh2V+HPgaMBYR10XEC9qJJUlqW6MdXpl5cltBJEn95S4dSSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVonHhR8QTImI6In4QEae3EUqS1L5GhR8Rq4APACcCDwBOjogHtBFMktSuplv4DwN+kJnXZOYdwCeAJzePJUlqW2Rm7w+OeBrwhMx8YX35ucDDM/Ols+6zCdhUXzwGuLL3uAN1CHDjoEP0oKu5weyD0tXsXc0Ni2e/T2aONl3I6oaPj3mu+51XkMzcAmwBiIjLMnNDw2UORFezdzU3mH1Qupq9q7lh5bI33aVzHXDErMuHA9c3nKckqQ+aFv43gaMi4r4RsQ/wLOCC5rEkSW1rtEsnM3dFxEuBzwOrgLMz86o9PGRLk+UNWFezdzU3mH1Qupq9q7lhhbI3+tBWktQdnmkrSYWw8CWpEBa+JBVi6As/Ik6IiPsOOsdydTU3mH1Qupq9q7mhvOxD+6Ft/Z08nwB+AfwWODMzvzTYVIvram4w+6B0NXtXc0O52YdmCz8iDo+Iu8266pnApzPz0VRP7uSIOH4w6RbW1dxg9kHpavau5gaz7zbwwo+I+0fEZ4GvAH8VEbu/fO1XwH719CTV90w8IiLm+zqHFdfV3GD2Qelq9q7mBrPPNZDCj4j9Z118MHBdZq4DLgbeWV//M+D2iDggM38GfA84DFi3glF/R1dzg9kHpavZu5obzL4nK1b4EXFQRJwTEd8ENkfEaP2K9EDgqxERmXkB8IuI2Ej1JA4A1tez+D7VN8rdsVKZu5zb7GYvJbfZl559JbfwHw3sAp5I9S2brwbuVme4R/7fp8fnAs8GvgHcQvXHVcjMrwEnADevYGbobm4wu9mXp6u5wexLy56Zrf6j+k6dFwFfovoe/EPq6yeBl9XT9wU217c/lGof1ar6thHgv+v5rAWuBl4K/D3wQWC/tjN3ObfZzV5KbrM3z96PLfyTgCcBbwCOB95eX38R8Mh6+sfAl4ETM/ObVK9qEwCZOQNcCjw0M3cCz6XaP/UT4DWZ+cs2w876oOOPupR7jk6NOTjug8jumLuu9/RtmfU+pYyIh1K9xfgysDUzbwd+H7gmMy+OiGuBd0TE44FtwFMi4pDMvDEivg/cGhH3Bt4HPCciDqX6Tv3/oXrbQmZeBlzWS8495N8AvJDqbdE7gBuA3xv23HX2To55nd1xd11fSuZO90v9HIZy3Je9hR8Rd6t/GOPA2VSHCD0OeGt9l98C34uIfTPzWqq3HQ+sn/j1VMeQAvyG6q3JXsCnqb4e9BTgOGBLZv52udkWyT0SEXeNiHPr5V0LvDczb4iIvaheSYcud5394Pr/RwLn0JExrzMfFNVxxFuA8+jWuB8aEQfXv7zn0pFxj4hDIuIhEfF+4Hy6Neb71P3yGKpdFZ0Y8zr73hGxb0Scw7B2zBL3Pe0H/Cnwb1QH/AP8OfCSevog4NvAsXXozcC6+raT6sCH1NPbgTVUHzJ8Fthn1nL2arKPbA+5LwY+WV/3NuBFs+6zuv7/pcCbhiT37jOgn061D28KGANe3oExn539q8AXqFbqoR/3ep77A8+jept9E7CxI+O+O/cXqLYCn9qhMd8beDHwGeBvgCOB04Z9zOdkv5Bqa/zoOt9QjvuiW/gRsTdwFfA04B2Z+cf1Tcfuvk9m/hz4Z+BlVL8ohwL3r2++BHgMcEdmXgh8GPgU8AGqraZfz5pPa6+6c3K/PTN3v3puBx4QEW+tX4n/LCLuDnwOuMegc9fzy4hYAzwDeE9mjmfmNNUr/O77DN2Yz5P93Zn5uMz8Nh0Y94g4kmr3wWOB1wA7gR9RfZC2ql7m0I37nNxnUO3T/WGdZ/0wj3ntJVQl9x6qv4n91Pryb+plDt2Yz5P93VTryFOo1pmxiNg8dOO+xFex84BT5lz3TODSWZfvBVxfT7+E6pTfg4B9qV657z3rvoe0/Uq7jNz3rrN9EjiZaiX73DDlrpf1YuCN9fTureanDvuYz81eX94XuGedb3JYx53qF/Yusy6fTbXB8ORhHvd5cn+Iamv/0I6s658B/qSefgFwat0v3xjWMV8g+/OBV9Ud88lhHPelPqmTqA72P4tq98LrgPtQnfF16Kz7XQQ8vJ5+E/CvVG8v/3KlfgB7yP1F4PVUH3ysnXWfvevndEJ9+c2Dzl3neEad6xTgW1Sv+E+m2s1wyKz7DdWYz5N9G/ARqiMUZq/YQznus/IdTLUP+XFU+1J/RnVM9NCO+5zcj68vd2FdfzHV7uJ/ojpK5UNUu0Z+DowO85jPyf4jqo2E+wNHDOO4L+eJfZ7qhIAjqF65TgO+RlX+Ady9/kHde9aTPAa466B+GPPk/ijV6clHzbr9HlQfhD5wyHIfVY/vX1NtDTwbeBfVcbivrEtoWMd8bvanU23Znzjs4z7neVwOPLWe/gjwtnp6KMd9Tu4n1dN7zbp+aMec6oiWs6n2Z78eOBOYrn93h3Zdnyf7GcDfAY8YxnFfzlE6T8nMt2Tmj4G3UB0e9T6qtyUXUL3KrcrMHwFk5q8z88rM/NUyltEPs3O/lWrwD4yIB0XEGVS7fX6Z1X7mYcr9I6pTpVdntQ/zwvq6C6k+jB7mMZ+bfSvwHeDYiPiDIR936iMqoPrA/Mh6+s3VTfEZhnTcF8i9V33EzqsZ4jEHHgBMZeaNVO9QVgEfZ/j7BX43+7nAbcD9IuLoYVvXl3wcfv7uQf2/oCr812bmxyLiOcBVmXl52wGbmpP7FqoTFa6meiXem+pIgGHMfXtEbKZ66wrVVs7xwObMvHzIx3xu9n2o3ua+keqoqX0Y0nGH6gOyiNj9bYQ/qK+7GnhlRJwCfGcYs8/JvaO+bldEPItq/RnKMY+IVcB1VF8x8A/ArVRnmb44M68d5nV9nuy/onqH+3aq3T13YYjGfcl/ACUi7gI8gersrqOpDp/6YGbu6l+85ubJvSUz3z3YVEsXEW+mOmHjWKr9fmfWW81Db072LwCvzP6ehdmqiJgGXpeZn9x9MtCgMy3F7NyDzrJUEXE/qkMU76BaZ84H3pTV2aVDbZ7sFwBnZOaKfxHbYpb1F68i4kVUJz7845C8lVqSrubeLSLGgB+afWXMOtPzwVSH9u7qQtl3Nfdu9ZmlRwH/npm3DTrPcnQl+9D+iUNJUrsG/hevJEkrw8KXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhfhfHCg3qEW8+c0AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.hist(xrot=20)" + ] + }, { "cell_type": "code", "execution_count": null,