Caractéristiques des épidémies

parent 08d98b0b
......@@ -2146,6 +2146,306 @@
"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\n",
"le début de la période qui suit, la différence temporelle doit être\n",
"zéro, ou au moins très faible. Nous laissons une \"marge d'erreur\"\n",
"d'une seconde.\n",
"\n",
"Ceci s'avère tout à fait juste, cela confirme qu'aucune période n'a été supprimée par manque d'information."
]
},
{
"cell_type": "code",
"execution_count": 10,
"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": "markdown",
"metadata": {},
"source": [
"Un premier regard sur les données !"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fb2ee2b06d8>"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXmcHkWZ//d533cymdz3DSSQAIZEAoQQkPuMoFwra7w4fqwoi6uIx4K6LrigsqIo64Iit8gl4nIj9x0Sw5GEJEJOkpCQg9zXTGbe+v3RXW9XV9fVPe8xk7e++Uze962urqrurq6nnpsYY/Dw8PDw8BCRq/UAPDw8PDw6Hjxx8PDw8PBIwBMHDw8PD48EPHHw8PDw8EjAEwcPDw8PjwQ8cfDw8PDwSMATBw8PDw+PBDxx8PDw8PBIwBMHDw8PD48ECrUeQFYMGDCAjRw5stbD8PDw8OhUePPNN9cxxgba6nVa4jBy5EjMnDmz1sPw8PDw6FQgog9c6nmxkoeHh4dHAp44eHh4eHgk4ImDh4eHh0cCnjh4eHh4eCTgiYOHh4eHRwKeOHh4eHh4JOCJg4eHh4dHAp44eNQNlq/fjpffX1vrYXh4dAp0Wic4D4+0OOYXL6DIgKU/P63WQ/Hw6PDwnINH3aDIaj0CD4/OA08cPDw8PDwS8MTBw8PDwyMBK3Egoq5ENIOIZhHRXCK6Kiy/kog+JKJ3wr9ThXOuIKKFRPQeEZ0ilB9CRHPCYzcQEYXljUR0f1g+nYhGlv9SPTw8PDxc4cI5NAM4njF2IIAJAKYQ0eTw2PWMsQnh3xMAQERjAUwFcACAKQBuJKJ8WP8mABcBGBP+TQnLLwSwgTE2GsD1AK5t/6V5eHh4eGSFlTiwAFvDnw3hn0m1dwaA+xhjzYyxJQAWAphEREMB9GKMTWOMMQB3AThTOOfO8PuDAE7gXIWHh4eHR/XhpHMgojwRvQNgDYBnGGPTw0PfIKLZRHQbEfUNy4YDWC6cviIsGx5+l8tj5zDGWgFsAtBfMY6LiGgmEc1cu9bbq3t4eHhUCk7EgTHWxhibAGAEAi5gHAIR0T4IRE2rAPwyrK7a8TNDuekceRw3M8YmMsYmDhxoTWTk4VE2vPT+WmzY1lLrYXh4VA2prJUYYxsBvAhgCmNsdUg0igD+AGBSWG0FgD2E00YAWBmWj1CUx84hogKA3gDWp7oSD48KYVtzK867bQYuuOPvtR6Kh0fV4GKtNJCI+oTfmwCcCOAfoQ6B4ywA74bfHwEwNbRAGoVA8TyDMbYKwBYimhzqE84F8LBwznnh988BeD7US3h41BytbcFUXLx2q6Wmh8fuA5fwGUMB3BlaHOUAPMAYe4yI/khEExCIf5YC+BoAMMbmEtEDAOYBaAVwCWOsLWzrYgB3AGgC8GT4BwC3AvgjES1EwDFMLcO1eXiUFX634lFPsBIHxthsAAcpyr9iOOcaANcoymcCGKco3wngHNtYPDxqAq4R89TBo47gPaQ96hJvL9uA0254BTta2qx1uVG1pw0e9QRPHDzqEj95bB7mrtyMeas2W+t6hxuPeoQnDh51CW/u4OFhhicOHh6O8AZ0HvUETxw86g6MMXTk4CwrNmzHv937NnbusutDPDwqBU8cPOoORZZNrFQtvuHKR+bh0VkrfUpTj5rCEwePukNa8ZAXJnnUIzxx8Kg7pF3svarBox7hiYNH3UFc7NPoHqpHJDw18qg9PHHwqDuwtIsvy3ieh0cnhicOHnWHtBxA9YlCBzal8qgbeOLQyfD+6i34aNPOWg+jLuF1Dx71BJeorB4dCCdf/zIAYOnPT6vxSDovUnMOVScKngp51B6ec/CoOxRZNkGRX7I96gmeOHjUHVKbslZkFCZ4nYNH7eGJg0fdIbUTnFc21A0efHMFlq/fXuthdAh44uBRdxCX+lR7dO/nsFtj5cYd+O6fZ+GSe96q9VA6BDxx8Kg7pDdl5Z9+0d6dwXOFL1rjc4UDDsSBiLoS0QwimkVEc4noqrC8HxE9Q0QLws++wjlXENFCInqPiE4Ryg8hojnhsRuIAv9UImokovvD8ulENLL8l+rhESLjGl896VLtdA5fuXU6Drzq6Zr1X0tw4u+3AAFcOIdmAMczxg4EMAHAFCKaDOByAM8xxsYAeC78DSIaC2AqgAMATAFwIxHlw7ZuAnARgDHh35Sw/EIAGxhjowFcD+DaMlxb2fHMvNVYtNbvKjo7GNKFZa0nU9ZXFqzDph27atZ/LcGfs1cxBbASBxaAr4gN4R8DcAaAO8PyOwGcGX4/A8B9jLFmxtgSAAsBTCKioQB6McamsUDDd5d0Dm/rQQAncK6iI+Grd83ECb98qdbD8GgnOr6HdIAO+Ars1vDiwzicdA5ElCeidwCsAfAMY2w6gMGMsVUAEH4OCqsPB7BcOH1FWDY8/C6Xx85hjLUC2ASgf5YL8vCwoZhxa1jtJcNbSVUX/H772x7AiTgwxtoYYxMAjEDABYwzVFdtd5ih3HROvGGii4hoJhHNXLvWJ0LxyIbU737VFwvPMdQSnjYESGWtxBjbCOBFBLqC1aGoCOHnmrDaCgB7CKeNALAyLB+hKI+dQ0QFAL0BrFf0fzNjbCJjbOLAgQPTDN2jE2HLzl3YsK2lYu1ntVaqHirX43PzV2N9Be9tZwZLfKlvuFgrDSSiPuH3JgAnAvgHgEcAnBdWOw/Aw+H3RwBMDS2QRiFQPM8IRU9biGhyqE84VzqHt/U5AM+zDsxT/+TReVi5cUeth7Hb4rCfPoeD/uuZVOd867638cislfaKiMuU08j1qz0ly61z2NrcigvvnIkLbp9R1nZ3F5QU0p46AHALvDcUwJ2hxVEOwAOMsceIaBqAB4joQgDLAJwDAIyxuUT0AIB5AFoBXMIY45nSLwZwB4AmAE+GfwBwK4A/EtFCBBzD1HJcXKVw22tL8P7qLbj7Xw4rS3vL129HYyGHQb26lqW9zo7tLW32ShIefmclHn5nJU4/cJi9svDuuyz4tdqmlIsYrdm8Exu278KwPsH8Wujt+DXwOgcRVuLAGJsN4CBF+ccATtCccw2AaxTlMwEk9BWMsZ0IiUtnQUtbsWxtHfXfLwDwkVarhdWbmzOdV701o7wcw+SfPYciA+ZeFbgctRb96qeCJwpxeA/prPATqdPi2qf+UfruIrphVd9RlrcjmRZ44qAGkz7rHZ44ZERWc0iP2qOxkG7a1+pRl1vnwC+jzRMHJSInOH9/AE8cMsNPn86Lrg15e6UOgHIvUlnae3z2Ktw1bWlZx9HR4d/tAD4TXEb43UXnRWNDLtUCsLvkc7Bdh0pRzSOUnnv4yPIPqIOh+uLDjg3POWSEnz+dFyP6NNV6CBZUZnaZFr1HZ63Eib+q79AwnijE4YlDRuRrEPfGOy+1D/26dwEA7DOoR6rzasUlVjO20rxVm6vWV0eFJw5xeOKQEYV89YnDHa8vrXqfuxPyueCZpfaQ3l0Wjd3lOioE7/wWhycOGdGQ97eus6JeF4F6vW5X7DabgDLBr3AZ4YlD5wPn9YrFZFk9wC9+HmngV7iMKOTqaVnZvZDWzt91Ud3e0oqRlz+OP89cHit/Zt5qbG1uTdVnJeBpgxmeeMbhiUNGNKR0pPLoOPj+X2ZXpF0eluO3LywslS1Ztw1fvWsmvvvArIr0CQBPvfsRXlmgDmH/zXvfLn1vLWPIF4/dH97PISMaPOfQ6SAyDGu3uMdXcpXVq2bEtpBjWLZ+u3N/afH1u98EoI7NJUaq/cFf363YGHYHeJ1MHH77mxEFr3OoKhau2Yo7XlvSrjbEkCerNu001v370vXYtD3IpfzWsg3t6hfoGCKdZ+ev1h7zIhV/D2T4FS4jauHnUM84/bev4spH57WrDVexyq62Is753TScf0eQ9+Db92cXCbV3mvzlzRV46K0V9oplxvaW2utIqg1PG+LwxCEjPAvafjwwczleX7TOqS7P8dAehzSVIlq1ePN6by/bmKmfcu5Av/PnWbisgvoKHcb++G+YuTSRjHG3hg+JE4fXOWSEn0ftx/cfDBTDafJYMJZ9N97m+NDEammUuLsbM3nebfWVMU6cHS+/vxZtRYbj9h9Us/HUGp44ZISnDbVBe+570XGdF4nIrjb3Hne3DcO2DBn5OjPE53duSBjrOQGXFytlxO62EHQWtCePRqsjdRD7cD2nEnjxvTU169vDw0ociGgPInqBiOYT0Vwi+lZYfiURfUhE74R/pwrnXEFEC4noPSI6RSg/hIjmhMduoDCyGBE1EtH9Yfl0IhpZ/kstL2qic/AUqV23wNX3rShUTOMwVy6xEr/GP01fVp4GOxA27diFJeu21XoYGvj3S4QL59AK4DuMsU8AmAzgEiIaGx67njE2Ifx7AgDCY1MBHABgCoAbiYhnV7kJwEUAxoR/U8LyCwFsYIyNBnA9gGvbf2nlRUJZVYN55BN4VScDn3ifO3pKzTkrNtV6CKlw1o2v4bjrXqz1MJTwe684rMSBMbaKMfZW+H0LgPkAhhtOOQPAfYyxZsbYEgALAUwioqEAejHGprFgpb0LwJnCOXeG3x8EcAJVM16xAzoAbfAWUlWCSICKGYiD6jmlsYRJM/Pnf9S5Qm0vXmvnGn7w1zn4zbMLqjCaOPzbFUcqnUMo7jkIwPSw6BtENJuIbiOivmHZcABicJkVYdnw8LtcHjuHMdYKYBOA/mnGVmnIE8ebvVUXfMGs9G1fsHoLnp8fyfrTcA6k8JFWlZUTHWoHVSbcM30Zrn/2/ar361/pOJyJAxH1APAXAJcyxjYjEBHtA2ACgFUAfsmrKk5nhnLTOfIYLiKimUQ0c+1adSyZSkEmBjXhHDrg5J23cjMem73SXrGdKEVUrfBNOOn6l2Oxl9IG6SsH0lxiLiOD/afpH2Q6b3eG3/DF4UQciKgBAWH4E2PsIQBgjK1mjLUxxooA/gBgUlh9BYA9hNNHAFgZlo9QlMfOIaICgN4AEh44jLGbGWMTGWMTBw4c6HaFFUIt5lFHFH+fesMr+MY9b9srlgnVvgXLKxgTqRzIKnz9YSeIs3TLK4uxc1f1zGk74OtVU7hYKxGAWwHMZ4z9SigfKlQ7CwCfbY8AmBpaII1CoHiewRhbBWALEU0O2zwXwMPCOeeF3z8H4HnWwch4QqxUkzHU7pZsbW7FV26dXjPzSq6Cqva0+OIt0+2VJLR3iB1L21Y7XP34fPz+pcXV6euxebi8QtF6OytcnOA+BeArAOYQ0Tth2Q8AfIGIJiBYJ5cC+BoAMMbmEtEDAOYhsHS6hDHGyf/FAO4A0ATgyfAPCIjPH4loIQKOYWr7Lis9WtuKeHreanx63BBl7t6EQroGtKuW5HLZx9vxyoJ12NbcimP3q77XaCRWKne75VuJa7Go7+6EZFeVwozf8mr7gjrujrASB8bYq1DrBJ4wnHMNgGsU5TMBjFOU7wRwjm0slcQf3/gAVz06D78850D80yEjEsflXXsl1um2IivlOe5o4Ne/IYxUWpU+GUsS6g7FT1YGaTYB5SJu1eZKlc9Wge6N5iWKMYa5Kzdj3PDe5RqaRwjvIR1i844gCqXOQSfxwlbgXWppNe+SailpS9P10nXbMH3xx2Xts2St1AmoQzUfU2flHHT3SJ7j3Rvz6oohHpu9Cp/5n1fx6KzKG0XUGzxxCNG7KdihbNrhtjOuxCJlJw5l79IZ3ErIhUAde92L+PzNb5StTyDaIXdEpbwJWRbvmoinqmwUqwuCKBfbuItFa7cCCEyQPcoLTxxCdCkEOxQxls6L763ByMsfx/uKiVeJhbq5zWyZ0bupofydOoJfbzXXZlVfKuLkQrA6mH1Dh0O1OTKdSXKi1PLceF6VzrZp6AzwxCEEn6ziTuWpdz8CALz5wQaFQrr8Y7BxDjxv9ZBeXcvfuQUR59D+tlpai7jpxUXOfQIoab1Ui4DLmKpBG8rVRyqdQyeVK7mKlWzIhTo613DsHu7wxCEEn5QqfTBjKoV0+yfjig3b8XchoUqzhThEBKzdXafCqk07cNaNrwMoz3Xf8foSXPvUP6z1VO971v5rvXSkWbvSLHQdiTS8+N4aLPvYzS9E51yYVrWXK3EOtX7Cux98PocQfK6Ksle+CLe0tuEjKedwOebikde+EPvtqnOo9oIwY0lEwEzXvX5bCwp5++i2Nrs5NikYB+Vq4fIoqiFWUhGuLIS8s4pIzr/978gRsPhn9hwIWrFSymvnm7lKPd5ikZW4k3qDJw4hmGFXrspdXIm5aCMOfLfVUUUJB//XM+hSKB8zGlNIG8RKLtCdVs5bqeR0Mow3DSFrz465Egufe1h0dXmCQze019JaxIcbdwAwhzn54ONtOP/2v+P+iyZjUEqRbBtjyHUo/qx68GKlEHxqucaqae9O5apH5ybKWiwOP7VincVubUOwETgVdI5OquvNGvH0lQXqWFyVuKXt5VLSPGdet6fFH0CFWoYjl69x7ZZmtLQWUz2Pqx6di7umfaBsT8Rd0z7AknXb8PA76c1daxFbq6PAE4cQJbGS8yahfZPm9teWJspsC+vuMk/lW/ynN9RB4MTL5eK+LOvuqk078P/umJnqnFEDuqNLPt3roRpaJrFSCvp6/TMLMvdTy4VPXswPveZZ/Nu9byXqyYT24rvfxH4/CgIrTBfEnabQ6l0bgufY3Jo+TlM96zI8cQgRKaSrwzmo0FGd4Cpt5qjLVcyE2xGJlVTchBlXPDQn9Zio9J87yvV80ixIy8LAgFnEQ+VIgbpxe0um4Hgqpfvf5q625k158t2PSoYbPbtG3JKJzjWGZurXPf0+bnguXZ4Izzl4RJZAjvUrMWVsi0IaR7RyIi5Wan/fiYgYmjaVhCCTDN9wTPckSW25ZuwnXXV9OxkayhK6O4sIUMaEnzyDL/4hvcOj1pQ1xV3s0SgSBzvnAAC/eiZdnog6pg2eOHCoxUr6F64SC7RtIvLjtZyv1ew7LlYy1KvQoLLmSlANJ82il8VmP0tMrh1lCof91rKNqc8Rd+Tiu5Tm0sVrNhGHQi77MpclE+DuAk8cQvC55bogvPBe9mRDOsJiIzh8otZSDFrNvuPWSnp7dtvCazqqux5CepPhclkrZZFzF3KEj7c2pzpn567qRDxVIZaKVeRMpXrut8KwkXMeVRI6Qn37a0swb2XnStGaFp44hOALTDXMRHUT3s45hMSh6hE0q9pdCeUTK6U/iYgyzIXyjNd1syrK+ldt2olDrn4W//f2h879lDORztyVm1LV14kq5WduuhXOIuB2mvuqcNWj83DqDa9kbleFNZt3Osd2qwbqnjgsWL0Fo3/wRMmzU+TOK0Un9A5ANp1D/LNaYLHv5e9ce9lCuSlNaCWIV3sefWzhy3C/XBezLyhk/dOXuEfDLSdxOO2GVxNlH27cgZfejzhs8bpiYiXhnKyhQ0zvanvmRzXDckz66XM45hcvVK0/G+qeONz/9+VoLTI8PmcVgOxy5jRwWAuVKGd8o6xI07dO4ekaAVQZR8m9eyfo2iMqj1gpC1zFSm9nkPWLKIdC2oQp17+M826bUfotEgTxGll89xGDiVC6G49kfzAqayVxTFc/Nq+sFk0bq5gvxYa6Jw6lPAFKhbQZWdlVvVjJIjsvHa4eddiwrQU3vrgwGkOKc7/751nt6lsVeK8jm53zoYlzKO14W9uKePdDuyxbFSk47N25r0pzoFuaW2O/2zTiI3HxTrOQi/e5UgYLKmtf8b7d8uoSzFyaSHe/W8ATB0nRmWa3mHXS6YiAs87BUo8xhlteWYyRlz+Ob9yTdCwSsW5rM9Zs2ak9/u9/mY3Fa6MESGu3uCs9n52/2qmeCyfFn0u1zHjbk9+gPUN8Y7HbQnPy9S9n7yREbCE2jLlcoeLFhVaUrZs88M330lGs5DY8JVRipXrxffDEIfzMEreo3N6TtoWPj9HW66K1W3H14/MBBJmyTJh49bOYdM1z2uNbpd1fOeB6i1XKQKWZqJVYmo6pD1IGcyVzP/bzF6/diqUfqzMRuiIN5+u6xpVL0ioutP9007TSd+ZGoxKIcw4ms3NpHCkWd1Vd+b3vqLHO2gsrcSCiPYjoBSKaT0RziehbYXk/InqGiBaEn32Fc64gooVE9B4RnSKUH0JEc8JjN1B4V4mokYjuD8unE9HI8l+q7gKDj0palcjQOgBZOQdez01xXSmcdeNrAIAX/rHGWM/dmsReh7+AmayVMu4d077y0xatS/adouspv34FP/q/d631yiGHB9RiO3V/KRo1QLXQNre2SVxECrGSYz35+T/01grnPlTj8ZxDhFYA32GMfQLAZACXENFYAJcDeI4xNgbAc+FvhMemAjgAwBQANxIRTwR7E4CLAIwJ/6aE5RcC2MAYGw3gegDXluHaUoHvalSiDB2ycg56sZJN58DCeub2y6lUVzXFFaGrNunFUcY2Heupo7KquAnzDTlqzEDXoQn9UeoIs+rove5zxBZ4kcP0/DPrzKqw3qkW2i/+YTom/yziXOUapvv39LxIbJnGWun1Re4WXUqxUoJzcG6uU8E6+xljqxhjb4XftwCYD2A4gDMA3BlWuxPAmeH3MwDcxxhrZowtAbAQwCQiGgqgF2NsGgtmyV3SObytBwGcwLmKSoOzo23Rttz53Kw7Kq2M3VnnYK6YxVs2C7I6n7netqxZ32T0MEQsNTXXkDLwXtRm0hJnwZqtmdpSwaQjSqMrcb2XWZ02Zah23G9+sME4JtcuTFfdHqsspVhJKttNaUM6nUMo7jkIwHQAgxljq4CAgAAYFFYbDmC5cNqKsGx4+F0uj53DGGsFsAlAf0X/FxHRTCKauXZtdg/leJvx32mme7U5B9fwGR0lX7Lz4qO5IvE6yFDXqnNwG0YMhOycUaXxrfveKUs7sTUug0I3rXTFxWdAfr66Pq7723ux36a95KK1ccKcRiykslaSz99dpUzOxIGIegD4C4BLGWMmWzvVU2KGctM58QLGbmaMTWSMTRw4ML2oQAW5Y5UoQ4esxCGzzsGROpRTJppG0Zc4rhmou1hJOMegc7BerWGg2vAZZdoOVmLdWGcIk5FOId0+sVLa+e8UBFbmHDQD++0LC5XlKshhQmQiZdpMuYiVWh3FgZ0NTsSBiBoQEIY/McYeCotXh6IihJ9cO7kCwB7C6SMArAzLRyjKY+cQUQFAbwBVMR5ORgh1P9e2Bo+8/HH812NJObTuRbSJaaLwGWZUy6uzXByKnlim5xLKhfYQB11oiHKh2RATqdym2Bd8aqR2wqUlDm6cg/TbsQuTKFV+BrJYyNSHWqwU/13LpEmVhIu1EgG4FcB8xtivhEOPADgv/H4egIeF8qmhBdIoBIrnGaHoaQsRTQ7bPFc6h7f1OQDPswrKRj7/+2n4zbNhkhTpdUrTqWmI/Nitry5JHMvu52DvFygP51CO29/eXblK6qEOn5FN91ENlKtv8Rq3t5THvNjlGRMos45M7sclwumG7S2x365RUU25y+UW5DkkE63ffvEgbV1V/XLkxRAxe0X7PN/LBRfO4VMAvgLgeCJ6J/w7FcDPAZxERAsAnBT+BmNsLoAHAMwD8BSASxhjPIjLxQBuQaCkXgTgybD8VgD9iWghgMsQWj5VAowxTF+yHtc/q47rXi7OwbSb0Mtw3RZ9K+fg8EJ9vLUZL7+v19uMuuIJ/PCv5iQ5Yi+yYlE+ngVZEvuoYL6tOtFXcrEZM6hHht4jrN/WgpGXP457ZyxLfa54DaZQ2+n8dOx1TM25cg68msu8/Na9cX1KOYhrghhIa7k8Lp4cSHUMSBKs1rbybj9O/+1reOpds39SNWBNPMsYexV6bvUEzTnXALhGUT4TwDhF+U4A59jGUg7ITl1JhbT5QZ9/xEjc8fpSAPqXo63I8M5yPfXXB94zdl3aobgSERPOvul1fBAGG9Sd/6fpy3Dk6AHWtoAggGECrjtLTXlRkQlOKWqytS+cIz4/E4iAS47bB//7wiLnftR9R9+5X8jdb3yAL0zaM107wvdyhdp2X9yzcbpiPzmQk1hp+Yb4nHQWXRmquUYe4BBNmJWcg0wcKiBWWqp5N6uJuvOQljlAmerZ5uK44b2jtoTKbUVWmjS/euY9nPO7aYlzbX04e0hbxuhCHHSEAQC27gwIaI7c7cdnKOLLaBXSrqlYFYIlFz3Pvj98MnYvxVPGDO6B7568r7VvApCXxukqamOaX98JY03tyqDAdF0k0+hK9PGZhPaQndPl4LVc7p9sdpqGAGn7t/glyH2IucOV1kqyDsPQ92X3v4P/fuof2uM6/PzJ9OeUG3VHHBILVsoFIJ8DfnrW+LBuVH7YT58tOfPMXmGOba9bNG0vQmsZxUombN4ZeKya/AOA+Evx0FvJPALtycAF6IigSucQ/93SVsTy9TuUbRLieRq0j5sokZdZVfXc22bgl0+/pziib39XBjGEM3FIoZL+n+ftFj/GzYEjjeNDd4k4Ku/CnQmyoZp8TL4k+X3pUohqqLgdWaxket8eevtD3PjiIu3xjoz6Iw6WiSIeV71oOaJSzgdxUqzb2lIKSmc18czIOZRkmzbOoZ3K5ObWQKbdtSFvqWlGY4Pj9NKM969vf4jHZq/E5p27EtFzbWhpi+TysWeaYmcte5qr+n75/bVOi6yILJxDLV1XdH27cw5BvRUb1AQ7S9/Jsdj750hwDtLJ4nNXKcTlaLMuYyxn7oxqof6IQzuPB8RBn7LSVF7qQ/uymfvmOgebXuTWV5aYG7JAnU87CdtL0egYfmL6ErXV8q2vLsE37nkb33lglmCtpBpIskiXetKVNhDiiZ8A912sGLlWdUZrG8NT767CuUKuAxucPYXL7K5LRAZON92uPssCmZYAKduw0OJkOAyBc1BMuKXr4sERXTj17z8421qno8GqkN7dYJM/2uZiPkepd7GJMTi+bG1FhpN+9RK+P2U/TBk31EnnsL2lFc9ZAuJZx8fin1nRs6tbqGcdceBYtSnacWYRM8Q8rSU9ikGq5CRWSjMOjl1tRXz9bnModRnuYqXywtSeq/TyqkfnorGQx6OzVtorZ+zDKFaSOQfpqhKpSUVdooNC2oVTf2OxezynjoL6Iw7Sb9NEUe24Xi/LAAAgAElEQVTCcgQr52CbK64TfuvOVixetw2XPzQHU8YNddI5XJxy0VHBNVicrdaoAd2V5b/4m15Gr0KOBIKsHIdCLqx5CBT+s4FgFitta27NpGgEMoqVhO/9u3fBx9sif4Ahvbrio82VC/WRVQzKce+M5fZK2r5dNwMGzkE6lFAoS4+DxY4pdFwp+ubo0608OTGqifoTK5n10VbkiEremLpF3jW6qowH31SHEubVuc7B1P60FBEndUjr3JT1uCtyFC3orjmktQua9GpzZ0gZRIRJo/rFysS+z71tBu6c9oFp2OE4kgPJYvoo9j28bxN6hsYCXfI5dG+MdEO6+TyoZ6PRSkurLCaTtZJxyGWBaxdGnUOCGJh3/jbOQW7Phdbn22mcUQt0vhG3E0kWUzouPHjVQiTuYssdW+kfH0mmhdLgWh3ESuVAXIGbXVBRrsUjNgRXwqWIjAoElkJie68uTOZg4Dh4z76x36I4QXb640r85DiSaBbMNZ13xsICxFgk8srnyPk+nzNxD+2x+2eqd/cmLqvcya6UfThenInblYcpL/g7JI9zsUt1Dmlze+oxVGcjVU7UHXGQ51DCdE74rpoY+VykkNaGMrYNwXUelGT/LByPwxZF8y7fN2OZ0TFPRCllKpknrd0qq3ycAxebKMN4K/sWj0c/drUVneTyqjomE9Q/vLxYWa66BeK8cl3YRa9oBoZCSBwKOYq1Z8rlMbhXV/z564e7dQjg9vMP5R0mx9PShttfW+LcVla43J8BPboY66l0eSI27dBbHym5Uvm3E3Fo3/FaoO6Ig/wM5Iki7zJl5ATioGUnbRPBWaYfrye66addeC9/aA7O/N/XsGbzTlx639upzs0K14Vv/yE9jcdFZzTne6d5wV3l/ao11hRDRzZvdIXrcxQT4hSLEeeQEwwkAGg3B7wX10RQI/t3w3H7Dwo2CIp7/sun38Mf2mkV54LbHAgQEVkU0nHI7/zmHXGRmuzcmmhP5hwcJrpp3i5csxV7/+AJaxvVRt0RB1vQLfEhqh56jqIgX6qF5jsPzEKzZQFy9/oMPne2FrFl564Yl5N1p3HD8wvwf++YrUacdQ6Whdp1iL0sCexj1kVK/UKyUHzO4lFZrKTtU7HKthk4h5Ub1Qphe6Rd+1iSbUYEs5CjhCe3CWkTQelqb9xhd2irFggpFdJSgWxiKzal1DlIz9SJOGiqbNm5K5MVVzVQd8RBfkiyTPPuN6KgaCrFYZ6o5BzG5czixPzLWyswyyK+cd0t8gWupbWI8Vc+HZuEWbnQPk1drHXEya+TK2/c3oIlkr23iIY8lU2sJK59aUM2yAjESvYFUl4ATh0/BLsMnIP2BbcMV3c9F3xqpH5sjJUW+Xwu7sltu7ZChiyBNoX/5Z/eP3Wb5cIPT/0EcjbOQVZIJzaE+vpKayVm/q2CKhPgzl1tGH/l0/jNc2qjCFd9S6VQf8RB+H78dS/i9xpZMaCW8RNRybmLx9VP+wyzxosRxRp6U00zenRVWy+LgfNsCjkAmPLrV4wmioVczpkDmWHxc8jFxEpJqHUOai5rV1sxk6NYjihT9E3VGf272wn0Zz45FIeO7Ks93iCEqHbhHPg9SJtfnAzWShyjB/bApSeOcW7zpLGDraJEV0zeuz+IbLGV4r9NomRAmv9KzkFqz3Gi3/JKfK2xWa2ZNiPVQP0RB+FBLjbsfAFgv8G9EmX5XMQ57Aw5h/SxjNT15UVDnnNxnUPKLi046fqXhbaDxgmk9S612dUX8pTKmmXxWn2O5ZgoJKVIDgCuFfwRXOMayTvwQo7QWmRgjLU7dlVLm53IQ4oBJaLIGPr3aAQQxMGKcQ6WtT+9WMlev5AnDAjH49Zm+UCh35HpiYj3+MARvSFPAc4l3vyVQ/DMt4+ObaDUnEO8zDVI3tWPz8cyIeCljbOuMW2oR+LgXleVQCSfA7qGMYN2ljiHdIuF3gY/jiTnEP2ev8qUqTV936pxMLDMcZoa8rlU90XM1ywvcK8siMxNXf0cdFQka7J5bqfeWmSpnNhUYxMzueluUeBsqW9zQI9gIzFxr34QgoiakroCQKyuK2yLWEM+l4ojCazg0o9DhcjnyI1z6NalEJMItBUZrno0yNa498DuGDO4Jybs0Qf/dea40nFTe2khmjzb9hiuxheVQt0RBxfw3YLqpSAidA2TgXBFVupcujrikJCNxo9vEpSAL76XTNSz7OPtMRt6ZR8OEy6mkMu4S87nzHJgE0zLjGubumEHYqX0e9deTcFusq3IUjmxqe63C+dApPfkZgju77OXHY3ffeUQ7Du4Z+yYCeUSK4nX1ZDPpRLVBdnlyrPw5YiQy8F44eI9LuTj3PCNLywUYmFFF3H6J4cBQILLAMpnom3lHGpLG+qPOLg8V75bVtXNCx7SfOFMu4C6vhgmhVRXRcTTo3/xQqJsSK+u7gMrwa6QtqEhJXFwdbyzLVQcukiprW0sk1hjWO8mAAFxKafDkq6lHAE6p1rGGAiE0YN6okdjAdecOV44FtV7e5noqBeavqYlDg51TCk6lW2WlXMIxZ+ODRZycRHUQkGcKXJq/N5nSS5lgniubcjVcDI0of6IQ4qds6qmKs5PWtmgq1jJZOPtGk57UK+4LNg03zgx4jRp9eZmfJghzDIA5PPpdoffEnwvzMHeVFu5ZJEuBWqWuEYAYhuCNK+s1flJMxwCaRdyhrjoralLXjgWdXjWja9LZ6mJgxjYUNmf5RoacrlUBNdFye0K/j66OsHlc7l4eAzhRHFTIm8ARbiu2e/JEQ8M41LBNV9GpWAlDkR0GxGtIaJ3hbIriehDKac0P3YFES0koveI6BSh/BAimhMeu4HCJ0FEjUR0f1g+nYhGlvcS4zA9jz36BbtD/oKpHl6XQq40iUqJ01OLldT15eLbX1uqbUPFObj0ZQp6p+KYxABvadCQy6Vii8V+TLvb9m6mxg3vnUoE8uDXD8fT3z66tKtkLN1LKw/31PFDpONRjZ5CcqUgeqyGODDzMRNUpx3+s+fTnSChSyGlzgHlM3Pmsc5U8ZH44lxkwJQDhmDuVaegQTKUiBEHqV3AzVpJh8seeCdRJja3O+gc7gAwRVF+PWNsQvj3BAAQ0VgAUwEcEJ5zIxHxbc1NAC4CMCb8421eCGADY2w0gOsBXJvxWpxgut2cA2AMmL74Y/xaEZRt9KAesYUCsC+g8ouwUGHzrKqnArdoGtjTzTokzTvYZtC1cJw5YRhG9G2ytlVoh5+DOTVp9pf1mW8fjQs+NdJtlxtWmjiyH/Yd3LNkEVRkLNVLK49Xzq4nLhADhGdKBHTRiGsCsZIaNtt4OQy5DS61u3XJpzNBKiPnkM8RGnI5tEoc4UNvfYhTfv0ynv/H6sAvJE/o3lhAjuKxqOKcQ7xdwM1aCQD+7+14JsTtLa2YuzJpNDJvVZQlstPrHBhjLwMwG6JHOAPAfYyxZsbYEgALAUwioqEAejHGprHgjtwF4EzhnDvD7w8COIFMQud2wnXB0nkR53NihNCg7J9/r88XDQAzpSBtM5duUNZzGRkPg+0qykozwfiOynRK3tF/oZCScxBhJA7ZmgQA7D2wR6DozTC9+DlFlo7grpc2DnJ0zrg/hriLJVxz1viSVVLsHOjvkd7YIew/46tlem+6W9LJyiCgbNSBCGgoJH1Q+MK8eO22gNMSOtdxDiL3w++TSgqpuhWX3h/nEuTnzvHt+2dF7ShrROjMOodvENHsUOzEvXWGAxA9o1aEZcPD73J57BzGWCuATQD6t2NcRrjc7oBtj37P/0mccSpxDmFruonAIcstW4tFDOrZiLsvPCz14FxM92JNMoY1jrH+XZIJ2YLxcbSHczBBqXJw7KY9O46IW0zH7F/2wKzYb9lDWZexrrVYxOBeXXG1oGwu1WP6a9HNCx77KS1tcEls1a1LPqUpq9kv4cefGevcVj5HKORyMQswIC6SYYgWftmbWmeqzTksl/AZKmzYZg8vkjVjZLWQlTjcBGAfABMArALwy7BcNUOYodx0TgJEdBERzSSimWvXqhWONrjZ+cfZdlHhFwwk+HDdGScWhGIwSccP7y31a0da4lBkzNlJJxKr6dvOOYoEZKsQETaHqTQRNtOAr19ZPaSD/ttnyig7oek8ubkPTYNCtFRkTK+s1oyN+3ekJg4OJLWxkFIhDfM9TJvnu0s+l+AcePNEgY6Bt5mTNjc6sRIQhkNXTEbR70aH9dujDeNp44cq61h1Dp2Rc2CMrWaMtTHGigD+AGBSeGgFADFo/AgAK8PyEYry2DlEVADQGxoxFmPsZsbYRMbYxIEDB2YZOkxLm7hLMk3QnMt2Sqyf2C0y5Agg6e67TAZboiEZjEWe3DaUFNKGOq5mgwWDE5xOll6CoXnVobSKuywcRK60IWifmjBBHGLfo188RLfKo5nptluwz4tuXbIlf0w0KxQQUaoF3WatlIYLyRGhkCetFdp/PTYP25pbY5yDTucg95unpKIbcCMOGwRpgi5UiE0/1OF1DiqEOgSOswBwS6ZHAEwNLZBGIVA8z2CMrQKwhYgmh/qEcwE8LJxzXvj9cwCeZxUkmTaRCRBM3HVb9KIiPoVcH5486YqMOzlJY3NoKwvnsKPFjTiUdA5WsZK9rSDwnvqYzYnM7O3aHs4hu2Ap0jmwdrH7SbGSmnPoEWZ30y2Uuh29bV70aCzgmW8f7TLUoB9BnMYxbdHH7cpTHnAO+uNpdOZEwKK1WzHzgw2Yu1Kt7F23tSXGNeruuXyribIHv9sqhHDvUsgmoKm1zsG6jSCiewEcC2AAEa0A8J8AjiWiCQjWs6UAvgYAjLG5RPQAgHkAWgFcwhjjK9PFCCyfmgA8Gf4BwK0A/khECxFwDFPLcWE6uNzuh9/5EE/N/ShW9rdLjy6x5nIO6UE9G7Gm5GWZhKxzKDKGXM6co1iHgsH+WgXG4oliTCh5hhu5K7O8mMOkkDZN+k07dhmJR3t0DiU4EAm5RpTgqX0mhkmxUvz76QcOw9kHD8chewUpStXhQvRhx11GNmawe9C7nLBh4vjCH95Q1Eunc0hzPLheDTEsMqzeHLx7j81ehQOGBaJaeQpxYhqY0Wr6lZ56kGVPf0df+f5xuPuND5TBO0XipCMOVp2D8WjlYSUOjLEvKIpvNdS/BsA1ivKZAMYpyncCOMc2jnLBuCsOJ4cqSuh+AmsoLhQAMLJ/dyNxkM3siowhT3onJxPkvm0oMlaSX9vAWWgT3XFWSOf0CmkTYfsfTfhijnLYfrvcdbkXUazUniGYOAcgCEVx7H6DtMf52PRxl9INzmZMQdJGSF8vVbfG5ygT0CIDdJLIghAsaruwW0+kAxZ0Drprke9pQz6XiMUl3t9BvRq1C7/YQ4MmoJVtf+dDdlcZaRcXlUIQ4kIB+4uT5BzintZpxsZDFThzDg7jk8dpVUgzYHgfs6+DLrbS5p27sEGXzB728MdKzsF4RhLtVkhb6irnTAj53DjnkOQIVCbLRcb0YiXNPuDMCcOU5V+9a6ZmpAFSqtecYBIrESUJqC5Y4t0XHoZ+QiTjh96KfA3k9nmYfWPuB+mW9moqYPNOfQpRExciLuyNWTmHzqhz6MxwkXWKdVRhKuQdhu0h71KIlQIPWP3YdC8Et5NPo3PQ1T1gWDwkeclaydCeq0I6n1OHz/jZE2bLqT/PXGE8/vjsVbHfrW3FBGdWDshLLwkbgva8tO9+uCn2O0YcFP2qujIZTGh3xRpW44OPzWHrORGyXXMqfY5Bb0VI7rRvemmRsu6RYwYAAJrCd3RLjHOIg9fJ5eL3SBy2THD7NHWJBbsE4puXQLGuvhCxVMtd2DgHRYV5Kzfj70td3c7aB08cBKgmuIo4yKy2bRMvJw1ioSliQucgfN+qyUmcT7mTKxb1u8njBPEFIEx8m1hJKjvv8L0S9WSrEI57ZyxLFgrQXTeHrAg9/OfP45hfvGg8R4aLeaaMSJxn9pC2KeyTzzyuHE0oRRVtMKR3gtM7v5nvhezTo4OqFdWiKDqRqvujRCC/j7fqRbZAxKnt1b9bqUx+Bvw9JmlemhTSPRoL2CpxDjHCArvTIQB00YiVVBz6K98/LjquOOfUG17BOb8zO92WC/VHHFIKIZoMnENpLbVxDpINNvdzEOfihUeOis0G3W6Ycw5p8izodpOXHDc69rvNQSGtYsu/cvjIZEWDbLecWGvQ9ejQfj8Hfb0vH7aX8bpV8nSOwL8mfvyY/QKT7V5CApqg+Xi9l753rDE1qz7Jj7o+13tEHJPm9BCqe9qg6PPNH51o1Fu1FlmCc7Dp5qJYZ1HZasnxk/sqBSKtqKKYXVHuR8UZiBstMoioxDmgV0gny0Rv81pbK9UfcUh5v2/68sGJMjl8xqwVmxJ1RMgOOlysJE7GLoVcbCLqlMiFDKasuqpySGj+othCcxQZw7aWaEelWnhyNmP2GsKFNkyXjBJc/RzyBuc/flyE7AQnr4MN+Rw+e+CwUuY3fo58y/fq3x3D+jRp54WOOKh0V9OuOB4/Ou0TAJKWeTqouAGVKKuxkA8WaENbss7GNXud+GSelzhMzjnkKC4WFTduKlGifNky5+Cy2dT5lqjuab/uXXDxsfsAaJ/ZdjmQzSNmN4VqCu6rMPuL/CGYFDNfjW3NMmuaVEjLSjqd41opAJyjQrrI9NncZFEDJ0imlvmiv1FQKqveXULtdz7lREmUWDS/tLYkRyqfFw6duEje7erq6UR5fFwqqIjD0N6RsYFqV66CajyyYpnXs4neCtKuxZXTM7XJoxjLHv4i55AU6SUJfULnYOEcnrr0KK1IT3cu1wV2Sie4zgyX9UqsonqwogWHbM2gwpPvxpWo3M+BiHBUqFSTN9o7Nb4J+wwMA++56hyYyXRPJg5Bn7bQBnJ7KrbfNcwGx6Ej+9orATh636ye8RGyiZWCT5tCmi+Iunso76aT4TPs95In+5FhYtb04TY0J/A2FeM01RMhBxkEQuJgyQQn6xzyZL6nLhZV3boIOgfh5RG5elnnqBJ/sYRYSaOQDov36NtNawxgey9rvbeqO87BuJtV7cZ0IhPwl9SOF6SUntzPAQDuuGASiozhN88uiE000Vrpxe8eiyG9u2LOh5swemAPXP34/FROcHqxUnz0PMWoqWWVEk61K+UxbUQ8Oksd6fbQkX2diV3/7skopWmhWrRkXHjkqNhv1xc2Cu2tts2Xb5XY3LqtzcowEPK9NHMO6gGOGdxDWW7j7mT9mk7H48w5gKycg6xz4PNL1t2lQVNDsNTJOjPR4VJ1DSbOAdBfBy/OEWk5hzteX6osj/nU1BD1xzkYjrluKHk9Lh5KCx4+Awjj0Yc5eOPsbvRr5IDu6NqQx6Ej+8XyCriAMeZMSFw4h5wisYrqHqgWAFV8eyCQB6fRobQXqkUrUUda2Tk9sXEOfGHb2tyK6xSJlZJipaAxbuL64JtJU15R5Lho7VZs3L5Lz61pBvfFSXsqy7dZQqvIlnln3/SarmaiRK2LssdWkokDEWHpum1Ytn67cawmlBTSEue7r0A0kzqHpIhQnn+66+D1iPQiPflZ33/R5FK/qr6qjfojDg433LZ0iLtI8R2VA2z1694Fpx+YdD5SKRTlPnULej4X9W3CGROG4dzD9zL6OcjgxOG1hR+Xyk45YHBinLK+Q7URV1k16URlJlm5jHLIYV0VnCJisZXCJeHbJ+6bqMeDCt704iL89oVkHmtdyBRZLxXvPKp39WPzAABrtiTDsOeItMYEWeNKyeEzlq9XpxRV6ko03A0MzmNAknj3airg2OtexIm/ekndd/hpere5WCnwkI7KRVP1hLWSok157ms5h7DcRBxk7NGvW6lfU9vVQv0RhzK0ITpEifPpqUvjAc0YYxjYs1GR/UsdclnH7org88xmytraxkoioF5dG4x1OZpDhfQf3/hAW4coyTko9TJI7nyadUp2w45XRtrd1ByFJZkL5yDL9FWmrEpFvEU+Lpu88+tpMARnE+fKW8s2AgC2NSfvpUqU1244WysloTolcv5Utzesd9cE52Dzxi/1ZzjWJFgrifqOokGspOJw5NdSn/I35BxAzsSB999RdA71RxwUN3zSyCDI2YEj+ji1UVoAYHaoamktJhbJltYiXlv4MTbvlEJISDNTdpzj0JkWbpHa27RjV0lhtme/yDnospOSu10Ol9DeKnGRzpQ1yTnor8ldh2KvN7R315Kif+Wm5E7X5WWVFwpVsh+TPqpR4R8DJOcLvxydo1RwTvS8+dhbFfND3hWXA6VL1LT7i899MhijcrOTPIlHI9Y9xoP36psQ6cl1px66B/4mbcRU9UQ0dVE7wYnf5WejGqe8MdKlzBU5B1fRcymsuCDCrCXqjjioZvmp44fgrf84CZNG9dPUSIJbMpjWmR272hKL6WOzA6Xs+6vjeaRl1lj2jeAoEQdpFZBFGL27NZQWaHFCd5MTFwlQiX3OPyKumFVdrmphUFk1qQjAPV89DLlcUqx00J5qQi2uibNXbFTWGdq7a+m72pxS/dDEgItyDVWyH1Ueb940N520gd8ieUGU2+T1SsRBMT9MCmlXnH3Q8Nhv8bpVOGdikL5FyTlo+jDpHHhOaBHyNR2978BYIMxos6a/9ohziBMtXSgN3m7SCS7++8uT90KvrgWMHSqFoilxl3rOQSda9jqHGkF1v/M5Qr/uXVKZOPKF1yTL5QprcYLpFn3ZHM+mc5APy+3+9KzxpQiU8WxX8fFO3CsyIVXt7A/fpz9e+t6xpd+qXZDOWkkc0VvLNihDYxyxzwClWOk/P3tAoi4Qf2FO/61aOSoukqrdve6lu/KRucpyIC5K5GerAqrx7hoLSSI8sn83NBSkXbHDVkQ0/TSFbM/nyJorw4ZffX6C1HeALFFZdecEAevUhhI5osQ9kqvpuCzTpXeLKaTVY0z6OZid4IBAZzFpVP/EU/z1c++X2tBZK40eJFmQUezD+zlUG6r7zReQNDF3OKtv4hwG9WwEZFZfU78U4Cz8zV/yq06PL5Kymduyj7dj/bYW3Prqkli93k0NpRchZq4n9Xvn/5uEl793HPI50uoE9urfPTpfMX6dLwhf8F9buA5n3/h6wms1uqZgMX/4nSiqZpYFgKM5FN0tXrtVyTnoHAhFIpcUK0W7ObM1NJ9LSUwZNxTfO2X/+FgYH5O+zVxO1HMELcvBHIHAs7jcQQhL8m8AS9bpg/SpiYO+LgMwf1XSei1HlHCCkzcOjRquzOQYqvOQjoXDkMVKChGqqgvZH2Lx2q1xsZJmlT1VSh8qZqsLUFvqUHd+DqoXmz+MaPdufyi5cGesYxxuv+BQ7De4J+6a9oHTM473Hcngj5GcviJP3eD40b94QSsq4lyLSenWvbGA7o0FdC3knPI+qDgltbVSdK+/dMt0Y5vcWunyv8wplenELC7PZk5oFnrN4/Nx4VGjEsdtSeV1Ywz6B/gD1YnTdPjmCaMToRTksO/X/tN4xZmR2C0SKyWfVSGXzKUMmEWJVvDNSJHhzP/VmbGqN1a6RxVYVamJ7L6DeyTCZ8j15I0Dr61LFQpEXJ5snmriHKDykFZQB/m043/5ktCmXqyUTEgUgM/9Zk1k5mqhDomDgh0vEYcUciXi1krqc3jEUzm5iK4Hea/Ad/uqiZWXZPTbNbbqXPkWCxGg6T+fc1MKq85X29y7y7+5fkLMWKdT4qUJOEik5mr0kUv1bak8pF31L1Eb+gWU36v+3XV6DLtYqZCnklHBorWRTishvkgBccxy+Or4IJNFJo9mxtRK9a8etTe2S7oveR7Jyn5O1I0cXWlXHo1Ntu5SBt6TTVkVndjmukg4e8YCKKrP4daFWxyiL1QSXqwkIA0zl6Og4o0vxGPNy2Z3svKNL/pyGIikziF4cVQ76By5LZJ8crss+jlD5jbVOEXwBfiksZFPRCB2szYX9K2wbBrSuyt+/5VD8Oq/Hxcrd2lTHOM1T8xPHNddpyzOiLcZ6XpYqSxZjy8+6mMq4sA5h+C3Un+DuPULoN4lFwQCf0K4e917YHfcccEkxRUloRqzs85BUaYNAkiBObRyoc2RIltevI7MOfzpXw4DAPTrYfeeFxXsu9qKsd25fA2qDb8qz4pK/CRC1E2JXLzOoa53UychDkR0GxGtIaJ3hbJ+RPQMES0IP/sKx64gooVE9B4RnSKUH0JEc8JjN1D4thFRIxHdH5ZPJ6KR5b3EOFQPUfRmdAVPevPs/NWx8ie+eVTMzC5H8UWXWwT9Rlb8SRYXJs7BtlN56F+PCOsFL4FJIR1vM152zVlRVtcvT94Td194mJpLyBGmXXE8/veLB+PSE8fg6jPHGePOJPtWxWsCTjlgCEb07RYrt7V5wxcOEl5GwmyFn4OOWIq0Qe5GzTmoF3IdVIsN78Y0B7kIEwCGhZuP7568X6JeIZ9LhJhoasjHsqWZoOQKw3tie5RKU1ZN3cA6jWmNM7ifw94DuoftmHUO+w7uiTMnDHPyCxCf47G/eBGvLFgnXEO8rirTW4syvIl5QymKK3VmtEB0jzl3IZunVxsunMMdAKZIZZcDeI4xNgbAc+FvENFYAFMBHBCecyMRcR7wJgAXARgT/vE2LwSwgTE2GsD1AK7NejEuUFmHtEkvpguNyGl2C727NcTN7BCfBFx00qSRA8vWSqrdLJfZqtC7qQEH79m3VE/mHHQEkC/Qc1dGi+mp4yKF2dVnjseRYwZo783Q3k3oUsjh0hP3xZcn72XdTcnXIy/YOrGSbQdbyJGQ2EVdZ4hg6irC5P+gSvYjt3/C/oME4wZ9GyL4c+RETxeKhF9376YG7D2wOz4tKTMBtUJ6oyEla7IfFbHjO+0MnIPWkTPw5NZGC84RHv/mkbjt/EPDvuPH1VZi5mB+pXEKhgUfbtyhPBb9Tq4Xu1ScA9w3QuI1JziH8DfP/6DLBlktWIkDY+xlAHJeuijG4RgAACAASURBVDMA3Bl+vxPAmUL5fYyxZsbYEgALAUwioqEAejHGprHgDtwlncPbehDACaTb3pYDKs4hnH2prJUUO21NxdhPHrhMl1eWg++qdPFpdH2L1ddsbkaR2bOrhQNFkcXDI6iSlJiUtvFxmPMaxHpWiJV0sOWaGNG3CV1DM1Jdm0fsMwBfnpyMNSQuzPKpkX+JIOIB8L1Toh38fkN6xhwkZfDmzzpoOP7p4BGxeqJdfOI84VoY0xPOLTtbsWDNVry+KNoN6/JZT967n7JcN2bb41ENSXdOPkzVKT7LV//9ODzwtcNLvw8Y1ru0g5YXXtW8JNKHDlGN02W+qTY4Ks4BFs5BhEgwdZZQnHMyKdh3WGJilQNZdQ6DGWOrACD85PkmhwNYLtRbEZYND7/L5bFzGGOtADYB6J9xXFaoHiLftaYSKyl2FSqICrC5Kzfh9teWhucndymqMalMMTlbru4vqv/+mi0AgEVrzXmCo3HG5cC6DFYuCDgmt1eGLxYiRLlyPsaWB/VUu6onvnkUPjmiT2ncopf5ryUx3sS9koujKOOVh85va0ysREE2Pe51Lmb3E4n3xL364v2rP1165td/fgI+d8iI2PXwRUMneuLKYJP59MylQW6RP71hTsUKAMfsOyhRZuJ2bDtj2Qrrur+9hy07WzF6UA889m9HJtpsYyymkB7Rt1vJCbU0Hskyj0Nl5qxSHquQJjSFnM9he0tr6TlMPXSPeJuum5uYGW38pMG9AmOEQi6IXNtiiEL70NvmXOvlQLkV0qr5xQzlpnOSjRNdREQziWjm2rVrVVWsUE2KLJGAZdHOT85QO22JWePmaaKSivX4+Kw6Bw3rIBKdfz12dPK4pn/O6ouTV5esxQU5S9Kbk8YOxn9+dmzUt6xzEPo+J1xIgej+XP/s+4k2eeDDKP5UdEyWu6sI+2aDNY6oyHx4FvfHiD+zHOk92GVCK4el4NVVTPOroVx89ead2rhcQCSyHNEvMorQiW6UymdDmY1LPjj0aB8Wiuy4x/5ZBw3HuOG9Y3W5AYJzuHC5XMNN83omcUyqcNgSwRn747/hG/e8DQA474iRYjXnjZDoYS/e0389dp/SsycKIjWLnMOGbS34zgOzhOuonHCl1EfG81aHoiKEn9y7aQWAPYR6IwCsDMtHKMpj5xBRAUBvJMVYAADG2M2MsYmMsYkDB2ZL+qJaFEpipZK5m/3GF6TQ1fsMVJsL8qbufH0p/jxTT+0j9j1os2St5GDKKoLHFAI0O3+tQtruTR3Uc5uUopxchZ+fPR4XfGpUqR/xer51whjtebzNDzckYybxRYOPsc2QH1gF8d7J84Srfjbv2IXfv7Q4HHd8TERRdj/RoVB1H0TrJ0DUOSTHdVYY0mJrcyvEcO8y/nDuRADxnbU2SquyTP+8VY/yZNE6jQinjR+a0KWNUZjR8jZbWi16DOkecajMk/kG49UF67Dvj57UtqmKTXb43v2x9OenJfuHniEQ360UjAOApPmysv18LqbfuP7Z9/GXt6L1Q+d1XU5kJQ6PADgv/H4egIeF8qmhBdIoBIrnGaHoaQsRTQ71CedK5/C2PgfgeVbB5KmqlnlaPn67XQKzyX4BuhHzpn7y2DzMWKqkebG+3TgH/Y7wSsGjWhct9T8+MxZ//vrh8fJwgc6ieFTXM3MO4iE5fIYcxVYEH98jmsRBQEQk4lm+pP4VYxvQI+ljEI0xuajIIpccEdZvawEAXPd0xNmoCHm0K47vIlW74r3DjUdbMXBo1E3Pk8YORkESOdqyp9kQicmS7fz+K4fEfps2LfF6wadJpg6I+oF4o+pwLcE9fGWBm0RBHKfOgpkMq76oMxR1Qi7g7654XfLzaMhT7P7I9z9L2Pm0cDFlvRfANAD7EdEKIroQwM8BnERECwCcFP4GY2wugAcAzAPwFIBLGGN8C3UxgFsQKKkXAeDk/VYA/YloIYDLEFo+VQriLe7eJY8ZPzwBR4yOUnUC7pyDuPjo9A+uunVR8ffgmyvw3Pw1yOdIu3vXvfSiAlIdUjrIcnboyLh8NwjRwKxKPfcIk2adjLzIljOfAyeK4sslj1pt0ixWiB/j91IViqRN2PWr8jKohizviosGzoEvpq1tzChWCs6P30utWEllhqucL3rOQZ6bQWyn+AQyXbsuvweHjmtRW3QFmxFVf5yjip0rEgfN/eScw8btLbjhuQWxY+IGRhWgT8azlx2DAaEfxpsfBLoh01xuyOdiyu8E91QF4mD1kGaMfUFz6ARN/WsAXKMonwlgnKJ8J4BzbOMoF8RFNZcjDOoZmTXyF8Yl3n8+78Y5uO/Qol3od/8cyBZ1CmExqUvProWSIrV7l3xMMehqWVRqkzGrc53z9ZB58rPYy+kus9URxbMPjqKJ8usWXy4XIi2n4hTBz1eZBXNiVMjnUFDshtWhq+N9iqIpGTytaZC4yfxcZXGebnOuasOsc4hfg+gDU2rT0WKIR9P93oOzjfV0+gHV5auCN3IcsU//WD25TS1xCNv84V/fxeNz4nngxffMhXMYPahHyeHuqXc/wuS9+0uRE+JjaMjnYmI3+drSvNtZUdce0ip3efHThEIuF9tF6uaGq3msymRQJ1eMiZWEE757yn5SvbTEwW7x4coJmbgbABgoiHDkpO+mLvShL5Ic065WfZuqZkzj5aercg5/7Zh9cP4RI3H+ESOV+alVzco72EiprSIOUd+2YI+yt7mO6OryO2vH6YB8LulgOGpA90S9VZuSWexU0IULV3uRU2hJpifEQLSoxomDun/OOSRyr0DWObiZYnNOlp9r2jx1KcQV0knuyd5fe1F3sZXi7GT8kEm+KiNHkDiH9sl2o3ai7zoORjRlNVkXqQKb6ogVgSukzeNz1zmoJ//ff3gi+nZriO180omV1BXzUntAnHNIpudMtsMMx1VKbp6NrXtjoaTrUT0ypUJaOmYWK3HTXGb0c+Dnm0I0cLh63+reCdU8ygtGGqeNH4qX31+LfQf3VNZLg0TfSlGr2wZNpeTWjYcv+rbwM6IZrSndK/deL+QIc1duwr0zIpPjtDqHaoTzrjviIMoG5UnharYHBJxDmwPn4PoekIJ1yGscmERT1hgnlLgeN9EBL2ewx2xyFytFYjIRqgQ5NrFSfCds7g+I7oPoKGTkRorM6DsSjDFoQOQcVM6Fao5E3x4/xueSknMQRFounIM4Rt3CprKhV92jKCuZvs/SOHPRvGwrslKoDxlyMh8dXH0sgIj7tFXllxhXBpt0DsnYZDz7XbxegAVr4km8RPB2GvI5nHbDq7FjQ3vH75Vsyirf/3KHZleh/sRKwk1OOqKFuwqHNyGRWEW3aLmKlRRlOs6BW4XMXbkpFpFVFkOlMXfjYiDbi5jGlBVws+LIWxZmES65mfltE3dx8qjFVngkU9Nj59fTZnGKUQ1PlRUuqXOIl4vgmxhOHEyiPaJ49FSXucyhi64b9B1fjJTjFAhTa7Go3ZHL5fqQ8wjbciEO0CqkxQ0hv56/vv2hUKZrFErOYY9+8XhfvG/AvMn5ZmiiPVxKLXr7+YfiC5P2iJUFCmk9B9jepE4uqDviEGMnZeJQqmO/8YU8xV4YvbVSuvGJ7azb2qJts40xTL35jVi5zDnoQjGoUIp348BCu0CX61rdJjnlkjC1J14r/761RSAOhoFz8ZDYdELGy81jLfdHngfD+zThf790cKKezChyoqdaUHlk3iXrtuGNxeuNnAMRlUK0ACZrJbcyntFOzi2gCgOeEziH1iLT5uQQy3t1LeDNH52kHqPCCEAHVbiWHo0FnDx2cClFaFAv+Hxt0cfWNrk5tnwPZUMRMVOfiRiffuAwAMnQOcftPygxP2U/B/kx6oIWlhN1Rxxi1kqyzkGjAFNB5hz01krl32lHu3xpTBoFu62Ml8tOcOp6btej82411XWBViEd0zmEfRs65y8qEOkmTFxTtIuNXlgVay+P78GLD0+IDIC4dZp4nkkp/IO/zon91o1TVKD+4NRPaOq53XTO9cixfGRTaCAKxQ0EC7pWZyb03bNrgzYIZcnT3Yk4xEObAMDXj9kbN587MS5yVIiq9KLK8LhUQV7cRc7BNNYoxa/9ehoKJCmko3OO2Kd/kGWywqhDnUOEY/ePx5dJxTnksifHMdUTW+yvCbWcp2SoCyC560yj+ONKYdvO2Pl6UnAO8kLVp1v8uvkutWdjQdueKEJSL3zx87oKu0n+8ptMWVU6B+W9clDcAkiEwuZER8k5OHCE4jEePuKPF07CUWPUkQRcOUDOObgEbxQdQ3e1FbX5MczZ18QxunMOPHHPba8tSZwfr5c8V2tMEo5Vfs5J4hBxLSKXMV4KGxLpjqKyE/ZPxrgCArHSViHWlziEe746WXlOuVF/xCG8yX84dyKO3U+dcMdlQXPlHNIqpMWJepxm4nCxkikeUdB3CuIQOsHZFF3LN2x3ao93vb3ZHj1SvkdnHzQ89vvCI0dh/IjeuGf6spKNvIw5H0Zxq0wZ11Tgi49JrKTSOahYe7lE9whkM03eliqKqs5wQgUS5P4mnZNSrKQo5JyDK3EoxjgHNXEQF3vbHM0REjkqVBB376a2VQTDxDmodA6yJ71orcQZywe+dngikGBJuS+01ytM7COjIZ8rifLaikxpTltp1J1Yib++w/s0lULjckRspL0VV2ulLB7SHPL4OPK5QKxkizmjOr2pi3o/wJ3gTJEgAeCV96Nw0PsP6ZmIuMnBd8w7hRhDv/vyIeq6wrgH9WxMErkcYfLe/UPRgXqn17tJdP5L9mG6Kr7bM21QXTmHBFHRtCdzqVyEoHrmMnEw7aSDxVTPhZigmqs8JafJRJND3DCZdA48XDkfr21MsjJcWS/UD4ixnJTWV4oNoG4zyHUJ4v0+bfzQBGdLEMRKJd1Rsr2SYYFgcbZqUzJGGBDm5gj7/c9H3sWL72ULNNoe1B1x4A9RKY9H9PBsyOfiiVXKFcNGbEYXi58nx0nGnJH7js7/t+NH44enfgKfUSSJ4f0XWSRHF71KRWwTlLzfPmnfRMTNaIzBp6gc7dNNvUtSKZOV9cKdqWpt7FIQlY7pOIdbXgkC6cXFShJXFrb52OwoppOKy0oYJmh1PCGnGP7eFV6UakGVF3mTMjJHZMwFkhZchLLNiQMMQ1iwIMubKA47+6DhpTDl/3TICPQN54Jt82SKIybXKzKGHkKOZrUndVAoOkhu1aTj5JyDOC/EaKwcfJcfhJ9hsX5EiCbJPDLDhm1qjkDMb/2XNz9U1qk06o84hJ+qOclfTBcb4kKOYhYc7eUcxBAJHHM1Ib5zOR4kTyo3mLL2aCzgq0fvrXW755OR7zp5Fi4ZXz9mn9J3U5gRbk75rfveKZXpOCGxGdPtMoV6Fln1A4YlCZZJ2XzvjOXY3tLqpJAWLchkk0Qg+Uy0Todh8dWPzQMQzTmVD4D8XOX4RfI4Wxw4B9WVniWJ84BoDolz/YJPjVS2KZrcBpxDdC2/+vwEXHfOgULd4Jjt9SAiJ8sclbWSKX6U6CA5M4x1pKrLEH+mouUTx4i+Tdja3IoN23eVuAxlwExBId3Dkgo0T5HXdTXiKKlQf8TBYBXCZaQuMs58jmJx47XWSo7j4s//USHa6JuaSct3SbJ4QZ5E4qJim2ClkBNtDL2bGmIKWxEnfiIK02yK77L04yDB0JJ1UaIhXfY72VvaNMYiY3j4nWREVp5PAAC+c/K+ieO2Jzr2x3/Dk+9+FNVP6ByS41Itpkm5t7o/fp1rQs6KL4AqzkFuwyRWChbTbGKl752SzEvNn40Ygny4xrlNFJu0FYvGzQPnil10Dq5+DklPalV7QaFbCs5ggRaJscpnhesNtjW34pZXF8f6ESFyDtzY5KtH763sWYzAUCPaUIfEQZP/F4gmrMvEyUucwzH7ts8qhM//qx+fb62rSo4DKBTSwtO1E4dAVPXxthatOAuQMrQZLo4HAxQJgo5z0CnlVGMsMoY/TlsaK7/unAPxrRMjgqBalFxMhNcIIrBk38kyFcGQxUom71sRu4ruOgfTYpnLRZsbI+egaMLUtzjXddfE685esQnvr95q7J8fs+oc4GgVSMn4Rrp83IA6XpK6bjxSsewAB8QX/TcWB2H5TZxDW5GhSyGHg/fsU8ppkqgrvOPVCLKnQt1aK6luN2eDlXliJQTEIdhN/cdnxhpstcv/YLkpq6pc17ct0myOCK87OAaJDkCml59bt3QRrC50UWaH9ooi45oiK/CcExuljG37D+kZGwtPuhOzPkqVjkXdt4irTldn/nNUOSRyGXAZuGqBlueQWSFNTgrpfj2SZtLKxFJcrCQ4Kepa5XXP+d00a/+FEnFw0Dk4+jkAsAZw5P2tEJJFjVCIB4FI0RznHJLvOef2xHVDde2in0Nb0RJ6PRcROxeLv0qgDjmHAEqdQ/jwbHHmeV2+6JkWXlfSoFq8ZFPbUpsaJZ3Jz8G2+3ClYfF8EfqTuoYKYm7tMmZQD+yl2HUB8RfOFG4kR4HuQF4snDzBU9KG/YbEA8bJt2+PfuoFxUW0ASDhEd5aLIJIn8gmXtcgVhKOm+blFyftmShTpt8MVwhRrKSDfL7pPeLXaVdIx3NEqGJzARGnOmvFplKZMv+6ojtd5kGuc7DtFXMKAmoWKwXzxPROcqK4c1eb02a1Eqg/4lB6eZMPhu/alq/fjm5d8njnx2q3fiBQqPHJYNohqSaJq9jjuycnZcC8P2VoYvm3UODCObhAVDLqTBUB4NdTJwAAhvcJuIKfnj3eoAwXx2EeY5ElPVa7FFSLQLwsLd8gJpAHgvj9ImHULWpJU1Z1vQZpzLvamDYgXVrOwRTEjyOfI21Mo1g9hYxe16w8Hd5frQ9Cx/V71llHcessndHAR5uTYcDVOSviZV87Zm+cqdAdBV1zk3Hz7OHvf0tbRAyVhKkUxDDMy2Gd6wyXCgYd1UbdEQcOk7XStpY2jBrQPWHPHKuboxJFN8noVbP/8W8elShTTT/d4isuACI2bI+LW0QxkyrPQKzNDDPBtPgM7tUVDfkoZpJOGQ3IxMG8m1KZsnZT+G7IzbjoHOLnJ8dx2CgxaYx6nF8/dp/Yb9Jc9v5DgtS03C6/ubVNe4/kzYcpvWZ8Q2B+qC7bAaXOwVKXQzRG0NW1zTtZv6aji9sVfhi6XNMcXQo5XPHpT2h1YZxzsFkvljiHVotYSdBNFB3ESkXmnva0Eqg74mDSOYg7N9tOWwynbVp4VRNAFllo+9ApMwmxaKwcAyQ5ctxaydxXFt2Ii5J7Ryha0OkbAGl3baKzFAR2k3dy3ZXEQeYc2qdzAOLEWnfpA3o04j8+MzYah6G9Y/cbWNJV7dxVLIngZMiPxqw4F0WJhs4dwe/jS++vTZTJsG1ARGS1VtJxDioPbqWYTCjqapiTQKRzsKk8SpyDQBxUfcvRdV02QnKdH52mjpVVCdSfQrpkraQQ98QWfPOkFYmHi8keRxfbKi3AZBWyWVDKnjlhGC47aT/s2T8u0xcnaJoX1xW2kOA5opLc2XTd4u2zxQ1S+Tl0a0wuqnIraTkHFcTnbBTZiLTO4XoAM+cg93XD1IOMbUbjtSx+ZTaWkB/xvx0/2lA3hc6hzc45rFUQTNV7KXIJOmLMEVhAMav1Yr6klxE4B6WlFJVMbovMrAfkHNMWgejtP6Qn/uUotelrJdCuFYOIlhLRHCJ6h4hmhmX9iOgZIloQfvYV6l9BRAuJ6D0iOkUoPyRsZyER3UDlnrUCzNZK7i9WPudGSMR2DhzRG49qwk2oVi/dXcgRYacwEYsMCcIQ1BPG6/ASpoVdJBApJY2cg9C3aRQ6sZJKLCATkN6O5rImuD5z8aV3uR4gUGaqbOiB+LMbO7QXpowbom1zV1HcvRo6t4zNBJvvBhAo7L+j0ZkB0XthNWWluLWQjnM49/CRiTLVMxLniu5+i2CwWy/yjZdIHHTWi3mKohtYdQ7SRP/xZ8dqalcG5dhOHscYm8AYmxj+vhzAc4yxMQCeC3+DiMYCmArgAABTANxIRPwO3gTgIgBjwr8pZRiXEqbwGbE8xJY748o5iATn4mNHa0VKaTa2OaJYCGXduVmc4NLARkADzsFsxir3rVIsluqFclhR37LkZ6cq63Jb/+P3H4Q7LjgUB+7RR1kvDURlvKsRgpnmBtezafsuPD5nFRatVcvo04iKFgttWDOuZaUOFj8HwH1zZRsCEeGd5VGwRd1cV3ltqzY8ovFCY8HGORg6FKBS2msTGOWCsOZtVrFSMnXuEfsMsA+mjKiEzuEMAHeG3+8EcKZQfh9jrJkxtgTAQgCTiGgogF6MsWks2BbcJZxTdkS2SmZllX1yxxOM69DgKKpKI/ZYtn5bSZYPuOVVLpe1Urx983GiYNdFBPRs1O/cxfu3RRPnhrdXFKxHDhvVT3vveUTMpi55HLufOrptWojB51wcvAA309y5Kzdp6wBxpbbtOYk2+yZrskpAlVNDh0IKnYNo9ttXYyRCRIk+Vc+oSz5atG2cQxB4zw6+FIjEQRs0k+vNijZrMrfo0JVEe4kDA/A0Eb1JRBeFZYMZY6sAIPzkb+ZwAMuFc1eEZcPD73J5AkR0ERHNJKKZa9dm0+JzttTkIAOk0zmY0BDbbdrHJUJn0y2bCPKAZjLEa7Rdz8fbophBrnTC+mKHfQ7t1VXLZqftj7HIL0J3f4CIKKdJlWqD+MxN7Yq32tQ9UbApcJG7R+eY635WSGKkW6BKbRmPpj8vzebK1QlOJK5jBvXAPV89TFvXJZ+JaEI8uGfXxPFY3+SWv5rPhWfmfWSpGVx3a9EtF3hnJw6fYowdDODTAC4hoqMNdVW3ghnKk4WM3cwYm8gYmzhwoNpBzAbT7RZ3aFYxTOy4vlXxJTG9CEdJ4TeW/OxU9Opql5PfccGhOE6zM86nIHbzhCB/k0epI7Im2nfkRnRxmuR6NnAZ/f6haO5nZ4+3tmki4hceOSpR1rUhp12ARAVmOcRKfAGw3cfYc7Tcqm7CGG3EISv2HtBdWZ5mM+IaeE9s5p8n7oERfdWOlEFdiTioxErCPTFtLoDQWsk8vKCfcJAvhGG1D9pTL8LkqVRt1kqBZZ4+jlU10K7ZwxhbGX6uAfBXAJMArA5FRQg/14TVVwAQPYtGAFgZlo9QlFcGBp1DZTgHtzb3GdijlOTm6H0HGneIk4QUjT276g3O0oTPELv77Rf11jC69tXHw74tK5rrDjZHgbx2V1sRh47si54G4snHZrrnsqMbAFzwqVFa2W6jowOgu1gpEFu4KGXFc0wQOTTbHBbnWF9NOHUZt50/EUeMVt+fLO+PlXMQjtsUwwniYFFI2+almNnOVk/EJcearbTaQi9/s2gy2AgdNiqZjrVayEwciKg7EfXk3wGcDOBdAI8AOC+sdh6Ah8PvjwCYSkSNRDQKgeJ5Rih62kJEk0MrpXOFc8oOkymr+KymLzbHGXKNeCkqMW0vwuZQ3r7OYMcOAJ+dEIkO+nfX735i0U5TEAeVU5kKtpcr2r3bFdcu4GKlIFeArc3g09WqiMNEbBscgw46K6RDTuhzYSwil/Zs98rGpalwxaf3x9s/PtmpLnfeUyFmHefoJ2S1eBOON1tNSuN9qp6vaC5sEzl2KeSc46zFfhvei8CB1ZwMiddTpSitJtrj5zAYwF/DRbYA4B7G2FNE9HcADxDRhQCWATgHABhjc4noAQDzALQCuIQxxrWqFwO4A0ATgCfDv4rAZMoqvnibDYpRwJ1zcLWNB4Bn568GAMxbpc7joGpzgGOicSvnAJGQODVZip+kA79eowc54gvov0oexqr2WtqK6NZonrounINqcTBdkyivdiU6pvueI3KzhhHasMmhbfdaBL/8NGoZ04JGGTgHE2clH7f5G8hdqvQFIoG3bZgaCzknQxH5vTZZieVzQfiX1jZmFPvxIJOPzKqcEMWGzMSBMbYYwIGK8o8BnKA55xoA1yjKZwIYl3UsacCftSmcrwve+2hL1KZhArlGMU0DcVHr7hAfB7C/CK6OaCLsuoTgs2BTjArdHWpgo3l7La1FNDiKTEyLs9KaxWBy2+Bsypoch66ei9JRbE+XmCYaV3phQBpLNdPCJ7Zic/bkO2Jb18vWRznLTxprtjqTn8nG7cmw3A0pOQcXyM6Ltg1JG2NoaSsaCXk5DSmyou7CZxQdrZVscE34HbNwcTD9dIEcntrpnBSyXdf7YIqXJLbpqrgGoFWuAxGBa24tOsiL+afBv0IxLtNYXYmD64tNiHv8/vHCSep6KealK0cLRAHtbIvgFyZFuhmbKITDZib6+OxVAKKkUDZ8etwQHLKXWf4uz1vVOyoSLTFstwqu0QxkfwnTPVq5aScefHMF1m5pxqzlehNm/hhdgiNWCnVHHFzFSja4cgENKXQOrtYlfPKlCcXRaLPpjnEObm3auBEuI262hEDnfdtEImIWLxs3EomV9HWUIQ4MbboaF7gmZwkU0hF10Nnwp0Ea3wYewM82j7502F6l76Y5Kt4Tm4NZqT1HTscljL6Mz3xyWKJMfIZvLzdzYWJe8j7dGvCXiw9X1pMJoev7YxIf8zmkiqFWLdQfceBflJyDezuiQnSSQRQSD9ZmYWMdF3veThrO0+oNKn4vE0u7bmugWBdj7KuQxpQVCIiOTawUEQcT55AsM+kyxIXRZoXkBEIsaVM5TE855zB5b7uVCxft2DgHZ6s34ZBLaArA/ZoHWXwSgKT6pl/3JLEV57YqYKMIkTP+5vFjtJxLMkZT+98fk2FEtVB3xIGzDurk49k4h/499EphkYjYuA1XZSJ/QdPoMGwvaxp/m99MnYB/nqh2vMsC18uIwke3OesxjKFNFNRhygH6uEX/PDESr5iivC4X5OQmyESkRxkWhDQ6By5itREH55hSwvXYdCO//vwElyGW8O2TknnBZbg4rIlQBWwU0b3RzSzYlRCmXCX1wgAADy9JREFUweBedmJYadSePFUZpkxwInSZvjhc6UiXCoiVovy7KYiDhXNQZZbT4YwJw3HGBHWClGxwlNELYiVXEVQa3cBBe/YxLpTD+jRhSK+u+GjzTiMxbVXlcFWOMa6Q7mGxwHKBqxUQEG0IbFylq4JddI5ThdAW0V+RptSEIb3Tcw42NFkMKoYJDmimeZRGvOsKFwfYSqPuOAeTzkHETV86xHjcxTkGSBcGPC1xSCVWsuxuXK+nEuC3xR5GIvhsaSta/Rz41Zh1A/HfJo9rjitPH4sBPRoxwMAtflGQ0ZtAIGfiYPK6FcHnW5q5YeMcXDnqPYQ0sLYFk891l6ZdHfRsGyAZNgIq6oBM84iIYqJl0zX94nOfLH3fZ6Da0xywv6/VQO1HUGVEsZXME8PmWj957yDExHGaPM8cceJgHpurMjGfQaxk2x3WMo4Lfxa2y+HXy5j9Xm0JLVVcA+QBbpY+U8YNxcwfnWhcUG1JZDhyubg4rxwWUJzGp+EqbQt5Fgts2+LmaiY658qT8drlxzvVHdTLzefnvosmA3DxGYnGaFOIT3b0ZBaf8c3nTtTWs1kCVgP1K1ay1LO9XF86bE+c+InBVna3wTF6KxC9pLedr580gBgLv3zWVVw5+cYVSheViqLkzZzC3NbGZb0YxrmZsWQ9LjlO16/cX3kU8e5EO3B06tm1oA2eyMHHetHRexvr8cixaURU3S2y9ywRe227eB5y3qZ4NYVIkZGW87Zth0TRpZyCt719A0B/hcKcw9Xaq5KoP+JgiK0kwiUujYscNJaNzfElMwUXA4ABPYNJ5bIGzbnyZKddGr8v1Q7zDESLj2usJsDdnt+UbzkRbqFMl+4qhgn6C2JFWReXsEmTHwgAbA09+7unIA42W/oszptXnX6A8fiBe/TBkaMH4MrTy5/A5pfnHIgzD9LrxPjV2BTYosOcLtggx87WiLMYP7y3tp44x00JqETO4W+XHl3KNV5N1B9xCD9t8sZyeTNXos19BvbAAcN64UwHpXCanRfgbndeCVgJtmhSmWGnZmov6L+6hJEndAlCKVg2I+GnLRf2UfsOQEOelMlvdGiymHRmuS3jRugXSCDgbO7+F3347fagsSFn1REA9tzQorjt9AOTPhMi3li8vvTdROhFS0PTfBOJQ0OenH1nyon6Iw6OGulKEAcbq89t7G1dN+RzePybR5VrWDFUgnO42BAvCYi4K9s9F99lm58DRxpntWq/f0SBIUBr0R5IkCLqYMTQ3k1YcI06O54O3coUUl1EGk/tcsP1XtrFSu5xmLjl01ePGmWsx9uxrQWiWKlSoddtqL3Wo0ZIs0stF2x5jG/80sG47KR9sc/A6rOQHJUgDt8/RZ9LGIjotG0REgOvuXIOtkXq68dEhMvF/LOcyBFFXsqOjmiVMCrrZZmX2bIE1o44iKlAVShdj0WslOYaTgnzel94pFknxNcVK3EQFPq1upd1yDkEn1aFdBnJ5tVnjsMjs1Zadx/D+zThmyeMKV/HGVAJsZJrpjN7OI5IruvqMGh7sUYNiPQ75d4PuFic8BAjrvk2bGKlNOjZtYAtO1trkl+8EuBSAdtOe2ioKzxp7OCy9f3lw/bEOYeMsAaj5PfatgkTRVq10AMC9UgcDPkcRGTZLenw5cl74cuT3Wzfa41ayDb5M7F13SzkEi6X2ELcPP7/9u4+2KqqjOP49wHElyvyJuAN5aVRSEB84cpkKUgjFWbShEyWIwjNpI1NJk2TZk425QyWOfnSDFFBWqlNY01YFoOWOpaJoKIiIWiOYoxSKiLmC+PTH3sd7/ael73PZe+z9+H+PjNn7r7r7rPOeu6656691tlrrSwbhzUXz2BIwlpJ0dr+UQGSekJLZk9gywsPcewR6eY7pPGXr56aOFkNaLjFaz1F9hySGqZRhxzAw5fPTuzJA7x/RMd7Nteqx8xS7aWR9jb09y77XswAT99rHFL2HIoa5+uLKkMlSUN58Y1X0g4rJU3fmBi2HIXuVUqzcNSoQYnnvGcp7mdeqrltacW0scNYe9lpWRTtXSMG7Z84nweS1yAqi/+Fi4c0t/EObXAbadxdS2ZmeqPCuz2HJhrPohraPvcfMO3yGX3NbV88iYtPS16/Jg/vvJOuNxffCSzpNsTKFqDnz2z8YfjxY4a+ewWZZtevLD3/SveS0S/tfqulr92MZnqT531oXH4FSbDz9eh3OLQju6Unsr6DLc2CkD01s4FTltrjkiBDc6YcxoRRB5dikkmZTBs7LHG9/GbduWQGL7zaeMtTSD/3JL70d9Kba+m8qSydN7XhORWfmNrJzQ882/K18/8baxDmJUyCK9qS2ROYMrr+FqEVV5w5mSsS5jjk5cpPH8PVqzdzWAkWraune1gp/XOKGsXoc43D2OEdjB3eeEJLmuWOJdmRIwdx5Mjk4ZXKnrudCZMK41f2ZxzbuXeFi7nik5P53PQxiZMPsxa/Ipxf8sah6Bsl0pg1cWTiJMGiVWaEv/F2ci/1psXTGdDP1DiY2ceBa4H+wE/dfWkR5Xjyu3MK/TCtLxo95EC+f9ZUZn2g8Rv7MyeO4Za1zwHZrlo5cEA/pjSY1ZqX+K25rZ6AJ8Wo9GoOTbEq7YwJjddty1spGgcz6w/8CJgNbAMeNLNV7v5Eq8uSdkEwydb82F4J9RyX4Z06ZfLzRScWXQRpkaEdA/nJgq7Uq+wWqRSNAzAd2OruTwOY2a3AXKDljYOU271fm9XyD47z0jn4QODlTO+SkvLLcn5FnsrSOIwGnot9vw3IZ+EVaWtjhrf2c4E8ffvMyXQOPoCZCcu+ixShLGMotQZcqy6nzOwLZrbOzNbt2LGjBcUSyc/QjoFcevrRmlMjpVSWv8ptQHzQ+XDg3z1Pcvfl7t7l7l0jRuhqS0QkL2VpHB4EjjKz8WY2EDgbWFVwmURE+qxSfObg7nvM7EvAaqJbWVe4+8aCiyUi0meVonEAcPc7gDuKLoeIiJRnWElEREpEjYOIiFRR4yAiIlXUOIiISBVLWhe/rMxsF7C5xo/GAM+myGIwsDPD8/LKU/Fk99pZ55k2ljxeO488FU925+WRZ1bxTHT35OWS3b0tH8C6Ouk7Uj5/eZbn5Zin4ilpPGljUTyKp0zx1Pvf2fOxLw4rvZLyvNszPi+vPBVPdq+ddZ5pY8njtfPIU/Fkd14eeeYRT13tPKy0zt270qa3K8VTXvtSLKB4yi6reNLm0849h+VNprcrxVNe+1IsoHjKLqt4UuXTtj0HERHJTzv3HEREJCelbxzMbIWZvWhmj8fSjjWz+83sMTO73cwOCekDzWxlSN9gZqfGnjMtpG81s+usoE17M4znbjPbbGaPhEchO6ub2RFm9lcz22RmG83sopA+zMzWmNmW8HVo7DmXhnrYbGYfi6UXWkcZx1J4/TQbj5kND+e/ZmY39Mir8PdPxvG0Y/3MNrP1oR7Wm9lHYnllXz9pb40q6gHMAE4AHo+lPQjMDMeLge+E4wuBleF4JLAe6Be+XwucRLSx0J+AOW0ez91AVwnqpxM4IRwPAp4EJgHfAy4J6ZcAV4XjScAGYH9gPPAU0L8MdZRxLIXXTy/i6QBOBi4AbuiRV+Hvn4zjacf6OR54XzieAjyfZ/2Uvufg7vcCL/VIngjcG47XAPPC8STgrvC8F4lu/eoys07gEHe/36Pf5E3Ap/Iuey1ZxNOCYqbm7tvd/aFwvAvYRLTt61zgxnDajXT/vucCt7r7m+7+L2ArML0MdZRVLK0scyPNxuPuu939PuCNeD5lqJtQvkziKYtexPOwu1c2QdsIHGBm++dVP6VvHOp4HDgzHM+nexe5DcBcMxtgZuOBaeFno4l2m6vYFtLKotl4KlaGLvHlRXTzezKzcURXNw8Ao9x9O0RvAqKeD9TeL3w0JaujvYylojT1kzKeekpVN7DX8VS0c/3MAx529zfJqX7atXFYDFxoZuuJumNvhfQVRL+YdcAPgb8De0i5R3WBmo0H4Bx3PwY4JTzObWmJezCzg4HbgK+4+6uNTq2R5g3SWy6DWKBE9dNEPHWzqJFW2Psng3igjevHzCYDVwHnV5JqnLbX9dOWjYO7/9PdP+ru04BbiMZ6cfc97n6xux/n7nOBIcAWon+wh8eyqLlHdVF6EQ/u/nz4ugu4mQKHM8xsP6I/7l+5+29D8guhu1sZlngxpNfbL7wUdZRRLKWpnybjqacUdQOZxdO29WNmhwO/Axa4+1MhOZf6acvGoXJngZn1A74JLAvfH2RmHeF4NrDH3Z8IXbNdZvbB0H1cAPy+mNJXazaeMMx0aEjfDziDaGiqiLIb8DNgk7tfE/vRKmBhOF5I9+97FXB2GCsdDxwFrC1DHWUVS1nqpxfx1FSGuoHs4mnX+jGzIcAfgUvd/W+Vk3Orn739RDvvB9GV9HbgbaIW8vPARUSf7D8JLKV7Mt84opVaNwF3AmNj+XQR/QE8BdxQeU47xkN0F8Z64FGiD6auJdwlU0A8JxN1YR8FHgmP04HhRB+mbwlfh8Wec1moh83E7qoouo6yiqUs9dPLeJ4humHitfD3OakMdZNlPO1aP0QXjrtj5z4CjMyrfjRDWkREqrTlsJKIiORLjYOIiFRR4yAiIlXUOIiISBU1DiIiUkWNg0gOzOwCM1vQxPnjLLZSr0jRBhRdAJF9jZkNcPdlRZdDZG+ocRCpISyE9meihdCOJ5qguAA4GrgGOBj4D3Ceu283s7uJ1r76MLDKzAYBr7n71WZ2HNGs94OIJiktdveXzWwa0fpZrwP3tS46kWQaVhKpbyKw3N2nAq8S7a9xPXCWR+tgrQCujJ0/xN1nuvsPeuRzE/D1kM9jwLdC+krgy+5+Up5BiPSGeg4i9T3n3WvY/BL4BtEmK2vCCs/9iZZCqfh1zwzMbDBRo3FPSLoR+E2N9F8Ac7IPQaR31DiI1NdzbZldwMYGV/q7m8jbauQvUhoaVhKpb4yZVRqCzwL/AEZU0sxsv7C2fl3uvhN42cxOCUnnAve4+yvATjM7OaSfk33xRXpPPQeR+jYBC83sx0QrZF4PrAauC8NCA4g2YdqYkM9CYJmZHQQ8DSwK6YuAFWb2eshXpDS0KqtIDeFupT+4+5SCiyJSCA0riYhIFfUcRESkinoOIiJSRY2DiIhUUeMgIiJV1DiIiEgVNQ4iIlJFjYOIiFT5Pw1TUawZ1zC7AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"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 démarrant en hier et chutant en été. Le creux des incidences se trouve en automne."
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fb2ec1d41d0>"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl4JFd97/09XUvvrX2bkWafsT3eF4zBmM14gQSMAScm9wZISEyIc1nfLLy5SXjzXt9Aci8kJIE3BBMMIbEdMNckNhgHs9h4m7HHnn2f0WgkjbaW1Hut5/2j6lRX7y2ppS51n8/z6JGmuktTXeo+v/P9rYRSCg6Hw+Fw3PiafQEcDofD8R7cOHA4HA6nBG4cOBwOh1MCNw4cDofDKYEbBw6Hw+GUwI0Dh8PhcErgxoHD4XA4JXDjwOFwOJwSuHHgcDgcTglisy9gufT29tItW7Y0+zI4HA5nXfHSSy/NUkr7aj1v3RqHLVu2YO/evc2+DA6Hw1lXEEJG63kedytxOBwOpwRuHDgcDodTAjcOHA6HwymBGwcOh8PhlMCNA4fD4XBK4MaBw+FwOCVw48DhcDicErhxKGLv2TiOTCaafRkcDofTVLhxKOLPvn8If/nDo82+DA6Hw2kq3DgUkVUNzKSUZl8Gh8PhNJV12z5jtVB0E1nNaPZlcDgcTlPhxqEIRTeRyGqglIIQ0uzL4XA4nKbA3UpFqLoB1TCRVPRmXwqHw+E0jZrGgRAyQgj5CSHkCCHkECHk4/bxzxJCxgkhr9hf73Cd8xlCyElCyDFCyG2u49cSQg7Yj32J2FtzQoifEPKQffwFQsiWxr/U+lB0EwAwl1KbdQkcDofTdOpRDjqAT1NKLwFwA4B7CSG77ce+SCm9yv56HADsx+4GcCmA2wF8mRAi2M//CoB7AOy0v263j38YwDyldAeALwL4/Mpf2tKhlEI1mHHgQWkOh9O+1DQOlNJJSunL9s9JAEcAbKxyyh0AHqSUKpTSMwBOArieEDIEIEYpfY5SSgF8E8C7Xec8YP/8HQA3kyY4/DWDglLr51muHDgcThuzpJiD7e65GsAL9qHfI4TsJ4R8nRDSZR/bCGDMddp5+9hG++fi4wXnUEp1AIsAepZybY2AqQYAmEtz5cDhcNqXuo0DISQC4LsAPkEpTcByEW0HcBWASQD/mz21zOm0yvFq5xRfwz2EkL2EkL0zMzP1XnrdKK4UVh5z4HA47UxdxoEQIsEyDN+mlD4CAJTSKUqpQSk1AfwjgOvtp58HMOI6fRjAhH18uMzxgnMIISKADgDx4uuglH6VUnodpfS6vr6aI1CXTIFy4DEHDofTxtSTrUQA3A/gCKX0C67jQ66n3QngoP3z9wHcbWcgbYUVeH6RUjoJIEkIucH+nR8A8KjrnA/aP78PwFN2XGJNUfW8cZhNc+XA4XDal3qK4G4E8OsADhBCXrGP/d8A3k8IuQqW++csgI8AAKX0ECHkYQCHYWU63UspZf6ajwL4BoAggB/YX4BlfL5FCDkJSzHcvbKXtTwUnSsHDofDAeowDpTSZ1A+JvB4lXPuA3BfmeN7AVxW5ngOwF21rmW1YcpBFnw85sDhcNoaXiHtQtEtgTPYEcAcdytxOJw2hhsHF8ytNNQRwHxGhe4KUHM4HE47wY2DC+ZW2tgZBKVAPMPVA4fDaU+4cXDhKIfOAABgNsmNA4fDaU+4cXDBlMNAzDIOiZzWzMvhcDicpsGNgwumHLrDMgAgzdt2czicNoUbBxdMOXSHLOOQ4saBw+G0Kdw4uFDtVNbuCDcOHA6nveHGwQV3K3E4HI4FNw4umFupM8iUg1Ht6RwOh9OycOPgQtFN+Aggiz6EZQGpHFcOHA6nPWk743BmNo2H9pwr+5hqmJBF65ZEAiJ3K3E4nLal7YzDjw5dwB9+9wCmE7mSxxTNgF+0xl2H/SJSKjcOHA6nPWk743DDNmv66AtnSmYJFSoHP1cOHA6nfWk743DphhgifhHPn54reUzRTfht4xCWRR5z4HA4bUvbGQdR8OG6LV0VjQNTDmG/yOscOBxO29J2xgGwXEunZtKYSRZOe1N1E7Jg3ZJoQESaxxw4HE6b0rbGAQBeOFOoHhTdhF9iAWkBaV7nwOFw2pS2NA6XbYgh6hfx4ItjME3qHFd1A36Bu5U4HA6nLY2DKPjwB2+/GM+cnMX9z5xxjqu6Cb9kZyvJIlTddKqmWxHTpPja06eR5K3JORxOEW1pHADgv752E267dAB/+cRRZ26D4oo5hP0igNbur3RsKon/8dgR/OeRqWZfCofD8RhtaxwIIbh19yA0g2IhbRkHVS+skAZauzPrYtZ63QsZrhw4zeUbvziDsXim2ZfBcdG2xgEAgrIVfM5qVuDZXecQYcqhhTOWEtw4cDxAStHx2X8/jO+/OtHsS+G4aG/jIBUaB7WozgFASxfCJezXxhREPRyeSGAupdR+IodTJ1nV+vy1sgt3PdLWxiHAjIP95ixsn2E91spupbxyUOs+57987Xn87VMnV+uSOG1Izt6cZVSeOu4l2tw4WC8/Z0+Aczfei/glAGjpWgcWiF+oUzkkchrmMxomFrKreVmcNiNvHFp3I7YeaWvjwGIOuTLKIWwrh1aWuoms9drqjTlMLlidbKeT3K3EaRzMrZvmysFTtLdxcMUcTJNCM2hJQDrZysbBVg61Yg53/P0v8MjL5zGxaCmG4rYjHM5KyGlWLVGmhT9r6xFuHGAZB9Ww3qDFAenWVg61Yw66YeLVsQU8fWLWUQ4zSQWU0orncDhLgSsHb9LWxoH1UcppJhS7EpoVwUmCD7Loa23j4FIO7jYibpjRPD2TwqStHFTDXFKGE4dTDRZzyHLj4ClqGgdCyAgh5CeEkCOEkEOEkI/bx7sJIU8SQk7Y37tc53yGEHKSEHKMEHKb6/i1hJAD9mNfIoQQ+7ifEPKQffwFQsiWxr/UUoKOcTCg2EFpZjAAINri/ZVYzMGkQLJCyq6iMeOQxsRCfnoedy1xGkXOUQ6t+1lbj9SjHHQAn6aUXgLgBgD3EkJ2A/gjAD+mlO4E8GP737AfuxvApQBuB/BlQghbcb8C4B4AO+2v2+3jHwYwTyndAeCLAD7fgNdWE0kgEHwEWdVweiixxntA6zffS+Q0+Ij180K2vGuJKaqkouPA+AIE+wQelOY0CidbqYUzA9cjNY0DpXSSUvqy/XMSwBEAGwHcAeAB+2kPAHi3/fMdAB6klCqU0jMATgK4nhAyBCBGKX2OWg7rbxadw37XdwDczFTFakIIQVASkNUMZxFkjfcAyzi0tFspq2FDZxBA5YwlpqgA4PhUCrsGogC4cuA0DuZO4qms3mJJMQfb3XM1gBcADFBKJwHLgADot5+2EcCY67Tz9rGN9s/FxwvOoZTqABYB9JT5/+8hhOwlhOydmZlZyqVXJCD5kNPyykF2KYeIX2hZ5WCaFElFx6buEIDKtQ5KUVfaq0Y6AADTyVy5p3M4SyZnv8cyqsETHTxE3caBEBIB8F0An6CUJqo9tcwxWuV4tXMKD1D6VUrpdZTS6/r6+mpdcl0EipQDy1YCLOXQqlWbKVUHpcgbhwoZS0zyM3b0RxGUBK4cOA2DKQfdpE4CBKf51GUcCCESLMPwbUrpI/bhKdtVBPv7tH38PIAR1+nDACbs48NljhecQwgRAXQAiC/1xSyHoCQUKAdWIQ0AYbl13UosjXXENg6Vso+KlcOGjgD6on4ec+A0jJzLdckzlrxDPdlKBMD9AI5QSr/geuj7AD5o//xBAI+6jt9tZyBthRV4ftF2PSUJITfYv/MDReew3/U+AE/RNdKXQVkoCEi7lUNIFlpWObDspBFHOVQwDna2UndYBgAMdQbRH/VjOsGNA6cx5FyfMV7r4B3qUQ43Avh1AG8lhLxif70DwOcA3EIIOQHgFvvfoJQeAvAwgMMAfgjgXkop+4t/FMDXYAWpTwH4gX38fgA9hJCTAD4FO/NpLQiIgl3nYKeyiu0RkGbKoTcsI+IXawakLxmyAtFMOczwzqycBsEqpAFeJe0lxFpPoJQ+g/IxAQC4ucI59wG4r8zxvQAuK3M8B+CuWteyGgRkAYmsVlU5UEqxBslTawpr1x0LSugISjVTWX/9hs143bYe9EX96I/68YuTs2t2rZzWJuuKa9Wr1L/70nl0BCW8bffAal1W21PTOLQ6QcmH6YThvCmDriK4sF90gmTuWEQrwJRDLCChMyRhsYZyuHgwhtsvGwIA9EX9SOR05DTDaXvO4SwXd9JDvYVwf//TkxjpCnHjsIq0dfsMAE6dw7ydrdMVkp3HQnbX1lYszmGtM2JBEZ0hqXIqq1Za/9Eb8QMA4un650BwOJXIagZEu7iy3s9aPK0W1OBwGk/bG4eAna20mNVACBAN5MVUWG7dUaGsdUbEL6IzKDvGsRilTBZXNGDNuqjUcqMdGItneE5+g1A000l4qOezphsmFjJaQaxCN0w8vGeMG4wGwo2DZGUrLWQ0dAQl+Hz52ELInunQihlLiZyGiF+EKPjqciu5A/UR24CmlPZsvje+kMWb/uon+OmxxhRitjtZzUCPrUbrSWWdt9+rbnfUf+yfxB98dz9+fpzHwhpF2xuHoGxlKy1ktQKXEuBSDi2YQZHIao5K6g5byqFcZ1a2O3MbB3Zeok2Vw4XFHEwKjM6lm30pLUFOM9DjKIfaxoG5M1VXDc4j+8YB8LYujYQbB0mAapiIpxV0BKWCx5yYQ4sph3haxb6xBXTaxrArJMOk5QvhFN2A4CMQXW1Fovasi1SbGoekHa+J1zlBj1OdrGagyzYO9aSyzqUtA8CUw3Qih2dOWCpulqdYN4y2Nw5sjvSFxRw6Q4XGoRUH/qQVHe/9yrM4F8/gU7fsAgD0RKwPZrxM3EHRzALVAPCYA3vd8XThQkQpxb++eI7P2F4iOc1ENCBCFnzIaPUrBxYP+/6rEzApIPoI5rhxaBhtbxxY6urkYg6dbaAcDk8mcGY2jc+/93LcYqcBMnfafJnsI0UvNQ7tHnPIG4fC+/XzE7P4zCMH8J2Xzpc7jVOBnGYgIAoIykJdyoHdd6YcfnpsBhcPRrGpJ4TZFM+gaxRtX+fA8vQzquG4WRiOcmihbCXN3m0NdQSdYyxTZK6scTBKajzCsgBC2lk52G4l1/2ilOKv//M4AO73Xio5zUBQ9iEsC3XFHOZShcohmdMwEAsgqxrcrdRAuHKQ8wtfxZhDC9U5sK6XkiuGwPy9lZRDQCp8mxBCEPGLbWwcSpXDz0/MYt+5BQDc770UNMOEblIERAEhv1hXthK777pJoRsmspqBoCSgNyrze99A2l45uCuii2MOoRasc9AMKyPJ7SrqDtWKOZRWQccCUtsah0QZ5fDk4QuI+kXsHIjwBWoJMNdQUBZs5VB/QBqwZkFkNQNBWUDEL2IuPbdq19putL1ycLd/KE5lFXzWpLhWijloZZRDUBYQlATEy/hrFd0oqI5mRPxi28cc5jOak/67kNHQF/VjqDPI/d5LgPVV8kss5lC/WwkAFM1AVrXauPREZCxkNOc9zlkZ3Di4jENHkXIAgLBfaKlsJZYbLgmFjQS7w3KBcnjxTBzJnFY2IA1YtQ6tOiWvFizmYJjUURGLWQ0dIQl9ET9mecyhblh7lqAkWPNTqiiHx/ZP4qXR+QLFltNNZFUDIVngbV0aTNsbhwK3UrDUOITk1poGVy7mANiFcPaHKqXoeP8/Po+H9ozZxqHUrRQJtG/MwV38x4L4i1mrwr43IiOp6CUT9DjlYcohIPlqxhz+5NGD+NwPjiCeVhG244E5zcjHHOyUbJ4Q0Bi4cZDdMQe55PGQ3FrKgUnuYjXQFZadoq65lALDpJhJKchpRgXlILVxEZzuBOnni4xDX9TavfK4Q22mkznHGAQlASGpcsxhMashnlbx8rkFzGdUbOi0su2SOR0mtT7HTDmUy7rjLJ22Nw7uTJxyyqHV5khregXlEJKcoi7WuyaR1S3lUCHm0K7tM5I5DZu7wwDyC9FCRkNnUHIWKB53AJ45MYunT5TvP7WY0XDT53+CB/eMAbDcuyF/5ZgDa1VimBQmtSYSAnAaRgYlwenPxN16jaHtjQNzKxFiDb4pJlRnBsV6gWUrSWWUw3zaMgpsN5zIamXrHAAgFmjfgHQiq2FTjzVeNZ62elIlcsytxBcoxt8+dQJffPJ42cdmUjkouolnT1mN8gKSlW2UVnUYZXp8nZ3LFPx7Y2cAALDAjIOcdyvNpfm9bwRtbxxYQDoWkCD4Sqe9hWWxRescCl9rT1hGStGh6PnZFomcVrZ9BmAph5xmtl1mCKUUKUXH5u68cUjmdFAKdIRk9LaoW2k57clzulkxLrVot4wftRf9gOTDSHcIJrXaoRdzdtZSDjft7AWQL+JkG5qgbVxk0VeQzcRZPm1vHPyiD4SU1jgwQn6hpbJynGwlX6lyAKwPW9zlR7eK4EqVA+vM2m5xh7RqwKTWNLyQLCCeVp2GhR1Byeku2krG4aE95/CGz/9kybMS1CrGgWV5MYKSgO19EQDAqZlUyfPPzqUx1BHA7ZcNAgA22caZDakKSAIIIeiL8PnmjaLtjQMhBAFRKBtvAGzl0FJuJROSQArmVgCuQri0igUn5sDcSmWUQ5s232NprNGAZKX/FhmHgCQgGhCdjJkfHryAt//N0+taYZ2aSWN8IYtnTy6twEzVjYobq0RRB+CAJGB7X9j+//LG4aXROKYSOYzOZbC5J4RfuW4EX//Qdbh8uANA3q3Euhn0RGSuHBpE2xsHwPJXdpTJVAIs5VBPv5f1gmUcSv/srL/SfEZ16h2YcqhU5wAASUUrOweiVWHGMBoQHeOwkLXuF1OffRG/E5D+1vNncWQyUXHS3r/tHcOPDl1YgytfPiwt94klXqeim0gp5WMIxckMQUlAZ0hGb0TGqen8nIzf/uZL+IPv7MfZ2TS29IQhCT689eIBR82y5AmWddgTlnnMoUFw4wDLf87cAcWEZRGq3jq+dVWvbhziadUJSM9nNFBqVa8Ww2Y6PHVkGlf++Y/w5OGpVbxq75BXDmJZ5QBYM7ZnUgrmUgqeOzVnn1d+B33/M2fw7RfOrcGVLx+Wbvrk4amyC30lmAuznHoopxwAYHtfxFEOhkkRT6v42fEZzKVVbOkNO89nG5YFV7YSYGcXtlCMsJlw4wDgf911JT52886yj7Va227VoGWNA8vPn0rkSna55d1KlnF4ZN84kjkd9/7Ly85C2MqwHW/MzkyaTuYcNxwzDn1RP2ZTCp44NAW2llYyDopuer5gLmcv8nNpFS+Nztd9Hkt+qGQc/KIP/fb7jr3HtvfnjUOxAdliZ4gBcCmHfLYSO+71+7le4MYBwPVbu7HVtStxw9p2t0rcQTPKu4k6ghKiARHn4hknA4RRqQgOAM7MprF7KIaukIT7nzmzOhftIdiCFQuI2NwdwlRCwVQiB8CtHGRMLuTwzefOgoV2ihc6hroOjENWNbCpOwRZ9OH+Z07XnbmkulpqF5PIaYgFJVw0GIVf9DkxsO19EcxnNMylFCfYzFRqOeXgzlZi37Mev5/rBW4casCUQ7pFpCoLSBdDCMHmnhBG5zKYz6gF7cvLts/w5xv63rSrF5t7wm1R95CPOUjYagdQXxlbgF/0ObvZ123vhegjOHoh6WTXuJVDMqc5BsFSDt52WSq6gd6IjE/dsgtPHJrCIy+P13keMw7llIOOWEDETTt7sWsg6hzPB6XTjsvok7fswm/euBU77GwmwCriFHykxK0UkHyev5/rBW4cahCWW085lHMrAcDm7jBG59KYz6gFEr5chTQLSAPAtZu6EJJbq3ttJdwBaaY2Xx1bKEiFvv2yQez/7K148Y9vxmfefol9Xt5w/tf7X8Rf/vAYACujx+s73axqtcT+7Zu24fqt3fjsvx9yVEElDJM68Yly6c5MOfz2Tdvw7//tDc5xdzori+VcOdKJP33n7oI55gAQEH1OsghzKzHlsJy6DE4h3DjUIORvLeVQKSANAJt6QhiNZ6AZFJt7SiW8G7/ocxTINZu77JTf1rhH5dANE99+YRRTiZzTyn2LfY8SOb1kUBQhBP3RgFM/4s7rn1jIOq4o1fC+WymnW2M8BR/BO68YQjKnOwt3JdzGo7imAbDcbB1BCYQUqtiNnUHIog9nZ9PO/1GpBsnv6m7A3qPsmFLDeHFq0/bDfmrRGcxn8bQCqkFLWmcwNneHwDZcBcqhjFuJEIJoQHJaRoTqnP+7Xtk3toA//t5BCD6CaEAEIQRhv4jBWAAXErkS48AIywJ8RSNVc6oBRbd2t6o9rMbLZFUDATmfDQQAaUV3khjK4TYOZd1KOb1gA8Lw+exCtqRSEugvJmC/j4N2ARyQD1TnNKNs8SanfrhyqMHGLqtMf3yhtKR/PaLpJvyV3EquD2uBcijjVgIsA/Lmi/oAsB5U3l7k6iGj6vjevvMlbok5u+rWMGmBS425ljqC5VOhy41UzWoGFN0aj2nS/EwDr5LTTATEQuNQq2uAu5q6UrZSLFh+b9oX9WO6DuPgdwWhGUHHOHj7nq4HahoHQsjXCSHThJCDrmOfJYSME0Jesb/e4XrsM4SQk4SQY4SQ21zHryWEHLAf+xKxTT0hxE8Iecg+/gIhZEtjX+LK6AhKiPpFjM9nm30pDUEzTEhiaUAaADa71MKmGsoBAB76yOvwx++wfOqVevErurHktgvN5G+fOolPPvQqjkwmC47H7ayYiwejTkdWIJ9BU2kBA6zgNctWYjOTFc10dteqYUK3v7xYT5PTDARla6lgMbhabeyVAuVQ6FailGIxqyEWKH/P+qOWcljMaoj4xYpuUOZKcisE1mXZ66669UA9yuEbAG4vc/yLlNKr7K/HAYAQshvA3QAutc/5MiGE/eW+AuAeADvtL/Y7PwxgnlK6A8AXAXx+ma9l1djYFcT4QgsZhwoftsFYALL9gesJy84OuVzMAbAyRliQMCQJUA2zJFD5qYdfxW89sLdRl7+qLGY0fOu5UQBwYgIMlk//3Y++Hl//0Guc49vqMA6xoOTURzAXklJ0r3K6iU89/Co+/uC+BrySxpLVDJdysGNwNRI0VJeRKw5IZzUDuknLdkEGmHLIYSGrVr2vzCiE5FLl4HVX3XqgpnGglP4cQLzO33cHgAcppQql9AyAkwCuJ4QMAYhRSp+jll7/JoB3u855wP75OwBuJsVRqiYz3BXE+RZRDkqVgLTPRzBiu9G6w7LzwQxUcCu5CdnuBrd6iKdVPHHwAiYXc5VO8xTffO6s4wKZThZeczytIiQLCNudPxnMrVQpaApYmU1s95yz74+iGQULaE4zcHYujZdHFxryWhoFpdRWDtaiG3FiDtUX32oxh4TdkbWycghgPqNhJqlUva9s0+Ie2OWOOXBWxkpiDr9HCNlvu5267GMbAYy5nnPePrbR/rn4eME5lFIdwCKAnhVcV8PZ2BlsKbeSXME4AFaswUesDy778FZyK7lhYxszWn4heGz/BHSTetJVUoxmmHjgubO4cYf11ptOFPbnmc+o6CrTf4vVOlRbxGIBsUQ5qHqhcsiqVpO6C4mcp9KmVcOESfOLrjsgXQ2lIFupyDjYhrJazAEATs+k61IOhW4lrhwaxXKNw1cAbAdwFYBJAP/bPl5ux0+rHK92TgmEkHsIIXsJIXtnZspPmFoNNnYFkVRqp++tBzSDFux8i7lmUyd29Efg8xHng1nJreQmWKZY8P+8MmH9n+sgrfCpo9OYTan48Bu2oiMolbR9nk+rTv8pN9t6w/iL91yOd16xoeLvjgUkRzlkXcVv7liMohtOT6Czs95JfmCB3WLjUCsg7TZ8xcWRiWz1QDNrqTG+kK1qdJmiDcmlMQevB/nXA8syDpTSKUqpQSk1AfwjgOvth84DGHE9dRjAhH18uMzxgnMIISKADlRwY1FKv0opvY5Sel1fX99yLn1ZbOy0grPrUT3sORvH4YmE8+9KFdKM333zDjz+sZsA5Hd29SmHQrfS+EIWL43OQ/SRAveJV/m3vWPoj/rxxp19ls+7SDnEM5pTs+CGEIL3X7+p7GMMy61kKwfmVtKNgt11VjWd3fiZ2XTpL2kSzD3DFt1wnR0DmHFwv3aGoxwquJXcKbKVssCA/PuyIFtJ5sqhUSzLONgxBMadAFgm0/cB3G1nIG2FFXh+kVI6CSBJCLnBjid8AMCjrnM+aP/8PgBPUY+VN+bTWdefcfjTRw/hr5446vy7WkAasOIOLMjsKId6Yg5yYaDyzIy1wO0ciNaspm0204kcfnJsBu+9dhiiYDWDK445zKdVdFfZxVYjaisHSmmBcihwK2lG/t7Nlg67aRbMOLAFWBR88Iu+mgFppor6Iv7KMYdKyiGWNw71KAe3cWCBcx5zWDk1i+AIIf8K4M0Aegkh5wH8GYA3E0KuguX+OQvgIwBAKT1ECHkYwGEAOoB7KaXsr/RRWJlPQQA/sL8A4H4A3yKEnISlGO5uxAtrJMO2cTg/7x25Xy+JrAbRNdinWkC6GGYcqsUoGKGiBoVsxkF/1O+MePQqPz0+A8OkePdVVhisP+rH3qLuo/Nptao6qEY0IMKk1hS5XAXjMJ9RnQ6uZzzkVso6yiG/AEf8Ys2YA3tt3WEZp4v+/ouu5oXl6Am7lUO1gLQdc5C5clgNahoHSun7yxy+v8rz7wNwX5njewFcVuZ4DsBdta6jmfSEZQQk37p0KxX7hit1ZS3HO6/cgJAslkyNK0e4qLU5K2Dqi/o971ZiCx3zdffHAphOKqCUghBijbtUdGda3lJhO+RkTkNWtWsbdLPAreSeXuYt5WBdo3t3Hq7HONh/856IjFfPLzj3EsjHHKIV3Eqy6HNmZVSa0Ai4Yg5llYO333PrAd4+ow4IIdjQuf5qHSilSCl6wQ5VqzDPoRxXDHfiiuHOup7LdmwsqMq6ZfZH/U4TNqEOI9MM2P1h7rP+qB+qbjo9k9hrWYlyACx3intH6zbcrAI7KAk4O+ch5aCWKoewX0SqRsyBGb6eiB+aQQtmkSdyGoKSUDUxoi/it4xDVbeSHXNwB6RlXgTXKHj7jDoZ7gqVHXzuZRT09BtPAAAgAElEQVTdhGFafu6sajiLdL3GYSk4lbPMrZTRnLoAAJ5OZ2ULGXOfsYDojB13YGNTy2Ur1UM04FIOrkXLPeNhzu7dtXtDzJ7j7Y1eXjm9MCANABG/UHcqa2/EupcHxhedxIjFKq0zGCzuUD0gXVohLQs+EMKNQyPgxqFO3ryrD8enUjg+laz9ZI/gDgTOZ1Rnga7UPmMlsO61jlspq6EzKDkLrpddS4puQHAF4plxYBlLrOliuTqHenDmbed0pwgOKFQOs7ZyuGxDDIB3MpZyRS2xASAki7UrpJlysA3qb/7THqf6e3Ixh8GOYNXz+2yjUo9ycKeyEkIQEPk0uEbAjUOd3HHVBog+gu++dL72kz2Ce3cXT6vOAl1PgHmpyPbwFScgnVHREZId14GXax1UvTAO0x8NAACmk9aCzaaNLVc5sJTNRLFyyLndSpYBunjIMg5eqch3lIO4vIB0T8S6Z0lFx9m5NAyTYnQug03doWqno89RDlUC0mUa7wGWIeMB6ZXDjUOd9ET8ePNF/Xhk3zh0D++C3aSKjANboKv5epcLIaRg4M9CRkNXSHJcWJrhqezkAhTdLLgnzKUxXeRW6govL5WVZeUkcoUxB3dDurm0ZYjYsJvi3k7NggXQ3coh7BfqrnNgmUfRgAjNoBiLZzC+kMXmGsbh0g0d6AnLVQ1yufYZgNXKmwekVw43DkvgfdduxExSwYtn6m011VxKjIO9QK9GzAGAPdPB5VYKSc6i6+VaB0UrVA5Rv4iA5MOMoxxW6lZyZyu5jUOpchjqCCAg+XDBI/2onCI4sb5sJdOkoJRC0Q2IPoIrRzrwa6/dhD+/41IAwHOn52CYtKZyeOcVQ9jzx2+rOpOhXPsMwEptbRXl8P/+x2Hc883mNK7k2UpLgM26LW6t4FXc3TAt42DHHFbJOIRdvuiFjIaOoOxUY3s55qAaZkEVOCEEfVE/XjgTx0ujccTTKqKByq2jaxGQrKl5iaxe4AtnAWm/6HPUScQ1QMgLOHUOsjsgbf2d3empjA99Yw82dQcREK1spJAs4n/eeTkm7Ey/p09YbW9GahgHQghqtd9kBj1UohwEKC1iHI5eSGBioTnvBa4clgBbHLy8C3bjVg7zGdXJIKnWPmMlhPwCsqo14WwhY6Uhyo5bybv3TNGNElfbu67cgKOTSbz3K8/hkZfPLzveAFgLXSwgYTGrlU1ljQUlZwJfyC9gsCPgGeWgaAYIKYxThWSrqK94dx5Pq3j6xAxOTqegGoWuOtYO/pkTswAK54UsF6bkihVdUBZaxq2UyOpNa8TIjcMSYMZBN73rP3fDFh8fsVIltVUMSANASLJ2lGnV6tfftU7dSgDw+7ddjJf+5G340Ou3IJHT0blMlxKjKyxjIaOWuJUEH3HaYEsCgV8UPKcc3GM4ASuVFSjtr/T0iRlQai1oql7Y/dfnI9jcHUIip0MSCAZjgRVf22u2dOHffud1uNTO8GIEJF/LuJUSOa1ps9m5W2kJsB23l3fBbphxGOoIYt5tHFYhIA1Yu153jn5nUHYFpL17z9QKVePRgITPvutS7N4QQ9S/so9Kd8iq+JVFy8WkGRTJnAa/6HO5R6z/Y6AjgOmEUtZts9bkNLPEp19pjvTPjlsuo6SiQdHNkp5cm3vCODGdwnBXqCEFkYQQvGZLd8nxoCQ4FfrrnURWKzthcS3gymEJiOvMrZRWdPiINY9ibg1iDixbyZn968pW8nLMQdHMqgbzV64bwdsvH6r4eD10hSXMZ1TkNMMp7ErmdMiir2SQzmAsANUwnfqKZsKUg5tybbtNk+Lnxy2XUTJXqhwAa+Y4gJrB6JXil1ojIE0pRTKnQzdpU9YcbhyWgLzO3ErJnI6wX0RPRMZ8WoWqr3a2koiMojvGoctV5+Blg6roRl1tyVeC1SvIijmwwq5kTocs+EpGcDKXixdcSznNKFEAkTIDfw5PJjCbUrChI4BkTrfjOIX3dLM9NW+1jUNQElpingMbpwqgKeqBG4cl4LiVPLzQuUkpOqJ+EV12EzOnCG4VKqQBq/leRjOcjqyFAWnvGlRFr78Z4XLpDsuYz6jIqIbTTI4FbZlyYG6lwQ7bODQxKL2Y0XByOmWNCK3kVnIFSvectdK7b710EIZJsZDRStTYWimHVok5sNbmQOGExbWCG4clwPykXvafu0krtnKwFyaW3icLq7NLDsoiMoqBeVs5dAYlp1WHl++Zqld3KzWCrpAMw6SYThTORZZFn7P4Om6ljuYrhy//7CTe8+VfIKMapTGHMgN/9p9fRH/U76R7z6aUEoN76YYObO4J4bXbSuMEjSQotUb7jISrSLIZQWluHJYAIQSy4IO2TtxKKUVHJCCiKyTDpPnmbqvRWwmwFg3VMDFrF491uJSDt91K5pq4lQDrb+JuJicLPsdtw/L1+yJ++Agw1UTlEE+pSOR0HJ9KVlYOLrfS/vMLuGK4w+kjNZdSS4xDd1jGz37/LXV3+l0uAds4eGxm2JJxV9Bzt9I6QBLIunIrReyYA5B3U6xazMFeNCYXswjJAvyisD4C0mUyaxqNu923u1+QXxJKlIMo+NAb8TdVOWTsnfdsSi3oyAqUBqSTOQ2nZ9O4YrgzP7tC0VctZboWAUmASb39nquHArcSNw7eRxR86yYgncpZxoHtWplxWLU6B3vnO76QdfzqTuM9D39QFd1Y9YXMPSgo4hec6Xx+wZfvLurP79AHOwK4kGheJX7GpQpquZUOjC+CUhQoB2D1UqZrwa53vRfCFbqVeMzB80iCb93sSJhyYLnoE4tWC4PVTGUFgNG5jFM0tm7cSqusHNwV1gFZcFwu7phD2FVL0R8NYLqZysG1Uy02DqLgQ0DKz5Hef34RgDUcKuaa7rbaQf5KMKWz3uMO7nkf3K20DlhvbqWwX3RaULM20KtWBGdn25yfz2Jbn5W2KHlcOVBK7ZbdaxNzAKyAKWs37c5WYgOTACAWFAsa86017myf4pgDYLnAmFvpwPlFjHQH0R2WC+ZCN0s5BB3lsM6NQ665biVeIb1EpHXiVqKUIq3oiAZEu3U2ccacrlZvpdds6cJ7rtmIN+3qwzvsojGvp7IyFbjau9yQbDWiU3XTMg5MOQj5Cmm3coi6Ft9mUKgcylePM+N1ZDKByzZ0OMcZzXYrrfd01gK3UhNeCzcOS0QSyLpwK2U1Aya1FhxCCPoifkysckC6MyTjC79yVcExZogUj6otdl2rbRwIIegOybiQyCEo5+cnFyoH10CdgGUcmtVCI6sa6I34MZtSyiqHaEB0smniGdVxXbIOtJpBVy1luhbBVok5ZHUE7WrvLI85eB9J8K0LtxJr180yYNw9cNYyi4QQYi8W3rlnB8cXnV2ZukbGAchnLLmVg190V0jn92oRvwTDpE1b4DKqjms3WymnoTJ9pSzjoDstHlisgRDiqIfVjuNUwt8qMYechgF78BTPVloHrBe3EnNJsOyRPjvuIPoIfA1oerYUZA8Z1BNTSdzx97/Aff9xBEBeOayFC6TbniQXtNN82f9b3FsJsJQDYDWxawYZ1cCWnjD+7teuxnuu2VjyeMRvKYeMasAwaUGWEvu5WamswRZxKyXtbsB+sTkV39w4LBGv7YIrwYwDC3Ky0Zer5VKqhiQ2J8Nr37n5kiyp+x4/AsOk+P6rE0jkNKdqfLUD0kB+7kBBzMGVreQeWsO6wKaaEJQ2TApFNxGUBfzyFRuchAY3LObA4g4xV+0GUxHNjjnkmtTNtFEkshqiAdGaic2Vg/cRBZ+n0zIZjlvJ3sX1R5lxWHv/tST41tygTidzeM9XnsVDe845x545MYufHpvBu67cgKxm4NFXJtYsIA0APbZbKSAJjutDFn24YVsPPvG2nbh6U5fz3EiZzqdrBcupL56w5oa5lZh7rpxyaFYq62AsAB+xmgGuZxI5DbGghJAkcLfSekBeZ24ltsiw3V8zdnOy4HM6wq4V1jwE4OB4foF4+sQMZMGHv7rrCuweiuFfXjjn+PTX4r44MQeXW8kvWG6lT7xtV8E1MKPeDOXAdqlBuXK+SjQgIeXqwOuub3DcSk0yDl1hGa/b3oPH9k+u6xYaiawVy+HKYZ2wXtxK03Z/I7Yg5ZVDE4xDE9xK8/bAoaNTSedYIqehIyTBLwp455UbcGQy4fSBWgu3Uk/E+huEZdHxx1daQJlRTzZFOVgLUbiKcmD1DJN2YWU5t1KzlAMAvOPyIZyeTePohWTtJ3sUSzmIVit8nq3kfdaLW+nsbBp+0YchezYAizk0YzfXjMJBNijnxFQSpq30FrOa09eo1+43xYzoWmTW3HHVBvz1r16FwY6A8/9VMkrRJioHZhxquZWAfGFloVupuTEHALjt0kH4CPD4gcmmXcNKyGkGVN10lAN3K60D1otb6excGpt7Qk5mUl+7KQfbOGRUA2PzGQBMpluLGDMS08nV7TflJhaQ8O6rrcwfd0C6HM2MOWTt2QG13EoAMGEXVrrdSrEgy1ZqTp0DAPRG/LhhWw9+dGiqadewEpxAf0BESG7OZDtuHJbIenErnZlNY0tP2Pl3b8QPQpqUrdSEgHTcNUOYuRbcyoF9n1lD5eDGncpaDifm0ES3Uj3KgRkHrykHALhsYwfOzqXXZdyBBfpjQckZv7vW1PzrEUK+TgiZJoQcdB3rJoQ8SQg5YX/vcj32GULISULIMULIba7j1xJCDtiPfYnYZZ+EED8h5CH7+AuEkC2NfYmNRfRQzn4lDJNiLJ7F1t68cZAEH7pDMuQmZCvJTXDFzadVZ3E7ZhsHlv0BWLMmAJdxWIOYgxt3+4zyjwuQBV9T+iuxhahcZTSDGYDxhSxk0VfQnK/Z2UqMoY4AFN10hk+tJxbsmFksKCEoiZ4NSH8DwO1Fx/4IwI8ppTsB/Nj+NwghuwHcDeBS+5wvE0LYu+YrAO4BsNP+Yr/zwwDmKaU7AHwRwOeX+2LWAmkdDPuZWMhCNcwC4wBYrqV2cSvFMyoGOwLY3BNyjEM55cBiDmu9y3WnslbCaqGx9gtbPamszO01Pp8tcCkBza9zYAzZE/VY0Hw9wWI5I11BWznoODObxvOn59bsGmr+9SilPwcQLzp8B4AH7J8fAPBu1/EHKaUKpfQMgJMArieEDAGIUUqfo5bG+2bROex3fQfAzaQZzWTqRF4HbqUzs2kAwJYi4/Duqzfilt0Da349zXArzadVdIdkXDQQxdELCZgmRSKrOQsX+55XDmu7kDF/fFXj4BebHJCuHHNgsZu0ajgxBsalG2LYNRDBtqL331oz2BEEAEwuNK/1+XIZnbPiZMNdIcet9FdPHMXHH9y3Ztew3MZ7A5TSSQCglE4SQvrt4xsBPO963nn7mGb/XHycnTNm/y6dELIIoAfA7DKvbVVpllspmdMwn9awqaf2cPazc5ZxKFYOv/Om7atybbWw2mesrdqKp1WMdIewtS+Mnx6bQUrVYdK8YgjJ1sCdmVRzjEM+W6mGcWhGQJoZB39tt1LxzwAw0h3Cjz75ptW5uCXgKIcmzsVYLqNzGQzGAghIAoKyAEU3cWQyiemkAt0wIa6BB6DR/0O5HT+tcrzaOaW/nJB7CCF7CSF7Z2ZmlnmJK6NZbqX/9cQx3PUPz9b13DOzaYRkwaltaDbNaJ8xn7GUQ380ANUwMRa3dmJsl0sIQUdQcmIha+5WqpGtBFhupWbGHEJVYg4ByedMs3PPcPASvRE/RB/BBQ+5lR549ixeGp2v+byxeAabuq2NIHPvnZlNg7pmwa82y/1ETNmuItjfp+3j5wGMuJ43DGDCPj5c5njBOYQQEUAHSt1YAABK6VcppddRSq/r6+tb5qWvDOZWWusMiJfPLWAqoUDRawemzs6msbkn3JRWz+WQBLKmAWlKKeYzGrrCsmMgT06nABTOb3b/vNZN4pwK6SrGoVkzHTKqNTa12u7U6r5qGQV3AZyXEHwEA7GAZ9xKlFLc9/gRPLxnrOZzR+Npx0tQnFI8vUbjY5f7ifg+gA/aP38QwKOu43fbGUhbYQWeX7RdUElCyA12POEDReew3/U+AE9RD+eeiYIPlFoZQWuFqptOUDVex67h1EwaW3tru5/WCr+4tjGHjGoVEHWHpRLjUJiPn6/kXWtDms9WqhL0DTTLraQ7nWKrEXXiN95UDoA1i3ty0RvGYSGjQdXNmp12c5qBqYTiKIfirDFWm7Pa1JPK+q8AngNwESHkPCHkwwA+B+AWQsgJALfY/wal9BCAhwEcBvBDAPdSStlW96MAvgYrSH0KwA/s4/cD6CGEnATwKdiZT16FZfusZSHc8amk45aZTVY3DhMLWZyLZ3CNq4lbs1nrudvMgHaFZKf478SUbRyCpcahGVk1XeHCrKlyNDMgXS1TieEoh4A3lQNgxR0ueCTmMGUv6ols9b8pc4EWu5Wc37NGyqGmyaeUvr/CQzdXeP59AO4rc3wvgMvKHM8BuKvWdXgF1tVUNcySwetLIZ5WIQqkrg/WoYlF5+fZdPU3xi9OWnH8N+zsXfa1NZq1HpDE+ip1h2X02+1DTs5UdiutdY0DALxpVz8evffGqgkGkYCIpKLje/vO49iFFP7o7RevybVlNKNO5SAWfPciQx0BPHl4qmkT9dxcsBVMMldeOVBKcWY2jXPMODhuJetvMdwVxPn57JopB+/+VT0KUw4rXew+8q29GOkK4Qu/elXN57o7i86lSpWDYVLc/dXn8JaL+3H8QhK9ESuF0yvIom9NZ0g7yiEsI+K32g+ctdN7YwXGoXnFWoKP4MqRzqrPifpFqLqJf/jZacwklbUzDoruzAGphuNW8mjMAbDSWRXdxIIdg2omUwlmHMorh6dPzOIDX38RV2+y3heOcrA3obuHYsiohlObs9pw47BEGuVWOj2Trrsg7cD4Ii7f2IED44uYTeXfGE8ensLuDTHsH1vAnrPzeHVsEUFZwJsv6mv6LskNcyut1e7NUQ6hfEfas3MZEJIfogO4lYM3u8iwQrOjF5IQfGTN7l9GXZpy8LJbaYNTCJfzgHGwPruJCsrhp8esDMx95xYQlgVn/gerN9k1EMW5eMbzAem2RWRupRUoB1U3MZdW6+qXohsmjkwmcP3WbgQkH+Zs45BRdXzkW3vx6197AX//05PY2BmEKBAsZjXcuMM7LiUgv/iulXqIp60PH1sMWNwh6hcLRqR2NDHmUA8R16JrmBSJNYo/ZLU6Yw5+77uVBj1UJc1iH5X+js+emsWO/ghk0YeR7pCzEeiP+SH4CK7e1In+WAAz3K3kTVjK40qyb2ZcC3wtTs2koegmLt/YgZ6wH7O2W+n0TBomBU7b7pI/v+NSpBUDX3zyON64szlpvpVgcRrNMFd9IZ5czGIqkYPgI04WDRt0VOz+WC/KgbGY0aoGsBtFRjUw3LWEbCUPu5U2dFpV0uMLzTcOU3bMQdVN5DSjIGY5m1Jw9EISv3/bRdjRH4HgUogDsQCe+6O3oi/qxw8OXsCJqbWZUcGNwxJphFuJ+R7rUQ4Hxq1g9GUbO9Ab9TtuJZaa+ZE3bcOB84t437XDCEoC3nftsLNT9grsnqm6ifAqXtrkYhY3fu4pmJR1oS1sV168sOaH0jSvtXQ1infk8xm1rgr5WnzhyeNQdbMkhvH86TmcnkkjqxoISvXEHLzvVuqP+hGSBZyeSTf7UpxsJcCKO7iNw3OnrJ5Jr9/eUzAulsESK/qjfswkFZgmLVDBqwE3DkukEW4l5jOsxzgcHF9ESBawtTeM3rCMCXv3cXI6BR8BPnXLroLFzWuGAci7bVa71mFiIQuTAm/c1Yc37OhxjrN7UryIOcphjdt11wtTDleOdOLVsQUsZBvThO/nx2dwfj6DP7z9ooIYxtefOYOnjk5DFEhdbqWhziBEH3EGJ3kRQgi290Vwys5WayYXFhUEJWs2QzKnFXxWnz01i6hfxOUbO6r+jv6oH7pJEc+o6I2s7mfdm58KD9MItxJLRavHrXRwfBG7h2IQfAS9Eb8Tczg5ncKWnrBnd71uHOWwysaB5Y9/8m07cc8b832k+ispBxZzaEKn2noY7goiGhBx17VWcwHWxnmlZFQdsynV2Wgwxhey0E2KnGbWZRx+6fIh/PjTb3LGn3qV7X3hpisHzTAxl1awoz8CoDBjyTApnjo6jddt76nZM4kpiLUISnvzU+FhGulWymlm1Uprw6Q4PJnAZfZuojcqI55WYZoUJ2dS2G6/0byO7HIrrSaL2fyAFDfsA1XcPdTryqEn4sf+P7sVt182CMCqsG0EacVSrPvHFgqOu/3y1TqyMgQfweae5nZerYftfRGML2SbMoeZMZNUQCmw0/7MujOWfn5iBlMJBXdevbHS6Q4D9rjftah18OanwsMwt9JK6hzcVr/a+L8zs2lkVMMxDj1hS1LOpVWcnU07uxCvI69RthIzDsUKoZJyYAN/vKy+CCHotK97vkHKIW0vkq+ezxdXphUdC66Adz3KYb3ANlHNVA8sU2nHQKly+Le9Y+gOy7j5ktrt9FlyxVrUOnDjsEQa4SKZcv1hM1V65xx0gtExAECvvci9NDoP3aTY0bc+jIPUAFdcPSSYcgiUNw7FxyOyCEK861ZiiIIPUb/YMOWQYcrhfF45MNXwwddvASFWdXmrsN3+nDQz7sAylXb2W8WprEp6LqXgycNTuPPqjXVl8vVF/fjQ67dge9/qKzYekF4ibCHRV7ALnnb1eqkWlD44vgi/6HOMQK/9gX3hjJXZsF6UA0tlVdbArRSUhJIPWU/Ejz975268rWhn5vMRbOwMejKIX0xnWGpIzEHVTWdjc+D8opP1Mm5PHnvTrj68/bJBbFuDxWet2NwTgo9YaeHNgrmSi2MOz56ag2ZQvPuq2i4lAAhIAj77rktX5yKL4MZhiYiunP3lMpXIoc9OSUtX8YMenFjEJUMxJ0jFlMMTBy8AwPqJOaxRtpJ7DGgxv3Hj1rLHH/no60uG1XiRzqDckGwl5ne/eDCKoxeSODOXxva+CM7bymG4K4gBO0bTKgQkASPdoaYqh5mUAtFHsKk7BELyKpcZjUakKDcab+tpD7JSt5KiG5jPaNhqB/IqKQfTpDg0nihIbWPl9JOJHD5+886SIimv0ogMr3pI5JZeJNYfC9TVKqLZdIYkzDfArcRagN98ST8EH8EXnjwOSinG57OQBII+j2ceLZftfRGcmm6icUgq6I1Ylc4RWXSqpKeTCvyiz5Ntz7lxWCIrdSuxmcVb7HkLlYzDuXgGSUV34g2A5Qf+k1/ejX/5rRvwyVt2Lev/bwbSGmYrFWcktQqdIbkhbiX2frtkKIZP37oLj+2fxEN7xjC+kMVQR3DVC6uaxabuEM7Pr16VtKqb+NA/vYgDriC/m5mkkq+3CUqOW2k6kUN/zO+pXmiM1vwkrSIrdSux5ltbey2XUKWA9EG7TfelG/LKgRCCD7+hvHvEyzC30urHHHRs7GwtlwijKyQ1JCDNlEPYL+J33rgdTx+fxV/84Cg2dAax0W410YpEAyLSqr5qzQunkzn89NgMrt/ajcuHSwvZZlKKo8qiAdFJZZ1OKk4GktfgymGJrDTzhvkYt9ZQDgfGFyELPuzyUOvt5TLUEUBIFvAPPzuFXJXU3ZWSyGqe7vOzEjqDEhI5bcUTCFmmUli2mhB++tZdWMxqODKZwMau1jUOYb8ISqunjq+ErP05rlSc5lYO0YDoZCtZxsGbrjxuHJZI3n++vA8pm/LEFv1KhTmHxhO4aDDq2Y6hS6EzJONv7r4a+8cX8aePHly1/ydRJSC93ukMyaA0H8hcLnnlYMVZrt3c5cyVaGXlELbjc6s1dpUZnXLFaaZJMZtSC9q4FLiVuHFoDVbqVhqNZ9AVkpyMkHLKgVKKA+OLBfGG9c4tuwfwris34Cd2z/pGY5gUSUX3dBO4ldAZakwhHNuMsGE+hBD89k2Wq3Kk23sZM40ibCcdsOrwRlNNOcxnVBgmLXEr5TQDiZzuVPB7DW4clshK3Upj8Qw2dYecoeHpIuPw8J4x3PyFn2ExqzmV0a3CYCzgVDGvhPGFLO76/551gvtAfkfdqsqhyx5ctNJ01rQr5sB4x2VD+OKvXol3XD64ot/tZdjrTa+ScsjYymGqjHJgbfb77NhC1FYO7P3r1TobbhyWSH42wfLcSqNzGWzqCcPnI1aHxiK30vf2jSOZ03HvW7bjl6/YsOLr9RIdIcnpZb8S9p2bx56z83jxTNw5xgJ8rWocmHJYacYS24wwtxJgFQPeefVwXf2U1iuRVTYOOZdyoLRwbSg2ArGgiGROd+KP3K3UIhBCIPrIspSDZpgYX8hisy3fw36hRDnMphRcu6kLv3/bxS230LHXs9KsGzYj+qQrb71S071WoZMphxXeu4yigxA4yrVdcJTDKjXfYzEHRTexmNXwyYdewROHrGLVmZRlBPIBaQmGSXF2zoo/8mylFkISfMsyDhMLWRgmdaohg7Lg+CoZsykFvdHW6WvjhhmHlbqW5myZ7q54rdR0r1Vo1L1LKQbCsujJvPrVZLVjDu7Y4Ytn4vjevnH87rdfxuMHJkuUA/tbst5p/TGuHFoGUSDLciuN2jsFRznIYoHM1QwT8xlt1Yd4NItGLXDllAOb5dCqxoFNXUuucI50RtVbquNqvTQi5qAZZsVNodtV+rPjVtLFcFcQn3zoFZyaTiMoCY6Bet02axDVIy+fh+gj6A55czPIjcMykJepHEbtNFbWAz8oCwV512zRa1Xj0Bm0PgSNMg6nZ1Mw7bz/vFupNf3mkuBDQPI5+fHLJaXo66btSiNpRCrrpx5+FR/7131lH3N7AH56bAaCj+Bz77kCim7iP/ZPoC+ar4Le0hvGlcMdSOR09Eb8nq1K58ZhGSzXrTQWz0AWfU4Aqlg5MPnZqsYhH3NYWVB1Lm3dp5xmOq2mW92tBOSzXFZCRjUQ8rehchxf5Z4AABgxSURBVLB37fWM5q3EsQsJvFo0IImRcW3yxhey2NYbxmu3dmOoI4C0apSMUn2X3YXVqy4lgBuHZSEKZFm9lUbn0tjUHXJ2CiFZKHizzqSYb9KbMnOlNNKtxJoQsrhDIqdBEkhLB1qtytqVGYeUojs1Du2EKPjgF30rcivNJBVMJnJQ9FIDk1UNRAOiY4QuGYrB5yN4+2VDAErTVX/5iiEQ4t1MJYAbh2UhC75ldWUdncs48Qag1DjMtrhyiAbEgnbFyyWeVnHdli4A+bgDa9fdyoHWaEAqGC+5HDKqXlDj0E5E/OKy3UosHkgpMBYvbeCX0wwEJcEpaLtkyCpgZbUjxcZhIBbAx966E3dePbys61kLuHFYBpLgW7JyyKg6Tk6ncNFgvldSyC8WGodUa8ccfD6CWEBakXIwTYr5jIad/VF0hSRngMtiVmvZ6mhGrAHKIa0YbWscQn5h2W4lliEHAOfiafzTL87gcz846hzLqAZCsuAYgUuGrM/5NZu6cPPF/XjDjr6S3/nJW3bhl64YWtb1rAXt+S5ZIVa20tKUw75zC9BNitds7XaOhSShoLfSbEqxshpa+MPbEVyZcVjMWs3nusMydvZHcWQyAQA4MZXEcAu3fwAs5TWxsLK202lFd1wf7UZYXr5ymE3lq/FH5zL45+dHMZdW8Ye3XwRCCLKagYAkOG6i3bZy8PkI7v/Qa1Z+8U2AK4dlIC3DrfTimTh8xGp0xmDKgWXctHKNA6MjKK2oBcScnanUE5Fx/dZuHBhfxLm5DI5PpfBal+FtRaL+lQek00p7u5WWG3Nwt2o5cH4Rp2bSWMhozvsxpxkIygKuGO7ARQNRz7bEWAorMg6EkLOEkAOEkFcIIXvtY92EkCcJISfs712u53+GEHKSEHKMEHKb6/i19u85SQj5EvG441hehlvpxTNxXDIUK3B9sHzznB3gmk0pLetSYqxUObA01u6wjBt39MIwKf72qRMA0PrGYYVuJdOkyGhG2yqHkF8s6UhQLyxZJBYQncpnIB/zytpupXveuB0//MRNLRH7aoRyeAul9CpK6XX2v/8IwI8ppTsB/Nj+NwghuwHcDeBSALcD+DIhhL1LvwLgHgA77a/bG3Bdq8ZS3UqqbmLf2Dxes6Vw8Squ2pxNqq1vHEIrNQ7Wh7Q7LOOazZ0ISD48sm8cftGHK4Y7G3WZniQakJDVDOjLbPqY1QxQijZWDsKKlcM1m7sKDAwzDhnVcDLlWsEwAKvjVroDwAP2zw8AeLfr+IOUUoVSegbASQDXE0KGAMQopc9Rq2PVN13neJKl1jkcnFhETjNxfdHONminFLICmnZRDivJVmJB+56wH35RwPVbe2CYFNds6mqJ2RfViARWVsjF+gqF2tQ4FNcVLYXZlIKIX3QSSi4ejCIkC04qdc6OObQSK/00UQA/IoS8RAi5xz42QCmdBAD7e799fCOAMde55+1jG+2fi497Fss41O9WOjxhBU2vGinc2TrKQdWhGybiGRV9kTaIOWS0ks6V1XhpNI6zs1ZWEnMrdYUt99wbdlitCF67rbVdSsDKW2jkp8C11iJWL+EVxhx6IzI2d1vdDa7Z3IVtfeG8W0kzWq4tyUqNw42U0msAvB3AvYSQN1Z5bjmtRascL/0FhNxDCNlLCNk7M7M6Q2PqQVqiW4ktaMVBqmggP8AlnlFBKdDbAoGsanQEJegmXVJK4cf+9RV88T+PA7DuZdQvwi9aH8Rbdg+iKyThbZcMrMr1eomYbRyWW+uQKjPLoZ1gXZCXsjFhsDGfW+ymmVePdGJHXwSnyriVWoUVGQdK6YT9fRrA9wBcD2DKdhXB/j5tP/08gBHX6cMAJuzjw2WOl/v/vkopvY5Sel1fX2ne8FqxVLcSW9DYoCDGzoEIAOD4hSRmk61d48DoXGKVtGFSXEjkHJ/vXFpFt0tdbe0NY9+f3tpyg5HKwTYTy1YObJZDG1ZIA5ZRNEwKRV96zGY2ZRmH67d247//0iX45Ss2YEd/BBOLOaQV3Upl5crBghASJoRE2c8AbgVwEMD3AXzQftoHATxq//x9AHcTQvyEkK2wAs8v2q6nJCHkBjtL6QOuczxJZ0jCXEqteweykFHRFS51F/VH/egOyzgymcSxKcv1tMVuyteqLLWFRjxtjVhk6iueVtBd5l62Ayt1K7GeVq3anLAWKxn4Y7mV/BAFH37rpm0IygJ29FubuxPTKai6iZDUWvd1Ja9mAMD37Mi8COBfKKU/JITsAfAwIeTDAM4BuAsAKKWHCCEPAzgMQAdwL6WU+RY+CuAbAIIAfmB/eZatvWEkFR0zKaWuQR3xjIauUGn1LiEElwxFceRCArLoKwh4tSpLNQ5sYDvLJ59LqRjuCq7OxXmcvHJYnltpctG6l0Md7Xn/2KS7tGKgJ1L/eYpuzXruK1L12/usX3JowprLEJRbKyFi2caBUnoawJVljs8BuLnCOfcBuK/M8b0ALlvutaw12+w3xemZdF3GYSGjVtztXjIYw7eeH4Wimbh6UycEj7bvbRSxJU6DYwPb42kVpkkxnVRwjauQsJ1YqXKYWMhCFn1O08J2I2J3o11KtpdumLhgG9XieOBQp2Vkz9gtXHjMgYNtvZbr57T9pqhFPK1WHOhxyVAMim7i2FQSV29q/UVvqbOQmXIwTIrZtIJ4WsWAR8cqrjZ547A85TC+kMWGjoBn5wesNssZFfob39iD2/765wBQohwifhERv4izc7ZxaLFYTmu9mjViY2cQftGH064xldWYT6vODOBiWPdGoLC1RqsyGAvAL/oKRnxWgykHADh+wTrHyz3wVxO/KEAWfStSDhs629OlBLjdSvXdP90w8eKZOHrCMlTDxMVDpS7fgZgfp2e5cuDY+HwEW3vDzpuiGopuIK0a6A6X7xi6oz8CSSAgpLQOohURBR8uGozi8GQClFJ87enTGJ2rfB+nbOUAwGmyN9CmxgGwO7Mu0S3y7MlZAMDEQq6tjUM+IF1fGvXZuTQU3cSnb70Ie//7LRjuKm3sONQRxDl7/G+rxRxa69WsIdv6wnUpB+Zbr6QcZNGHHf1R7OqPtvQUMze7h2I4PJHAqZkU/sdjR/D5Hx6t+NzphALmBWHGoZ44T6sSDUiYSSp4+988jf/ytefx2P7Jqs9/7MAkfu1rL+Dlc/OYSra3cQj7Wbua+ozr4ckkgEJ1X8xALADdbpwZbLFsJW4clsm23gjG5rNQa+RMz2fyjeIqcd+dl+Ev3nt5Q6/Py+zeEMN8RsN3Xx4HADxxaAqTi+VbUU8nFScB4DAzDm2sHKIBET87PoMjkwkcnkjgYw/ug2FWTqk+PmUtcI/tnwSlwMbO9jWsrL6j3oD0kckEJIE4KavlGOzIvxeDvM6BA1jKwTApzsWru5ZYfn5nmVRWxjWbunBNGwSjGWwn9u3nR9EdlmFSim8/f67sc2eSCi6203tPzaQg+Ah6wu1tHFTdRG/Ej3vfsgOGSav2qjo1bb0/Hz9gKYx2Vg6xoIS+qB+PH5isq0bpyGQC2/siVXt2DbrSgnnMgQMgn87KeqtUgrmV2rVwqxxssU/kdNy6ewA3XzyAB/ecK/nAUkoxncxhuCuEqF+EZlD0RuSWT/etRtRvbTLecfmg856qVjNyetZ6f7Iah3Y2DoKP4NO37MLe0Xk8dqC6Ow6wjMPuKi4lwEqwYPDeShwAwEUDUXSHZfzdT05WdS058wcqxBzakWhAwiZ7atvrd/TiDTt6MJtSnUI3xnxGg2ZQq5LcbpkxEGtftwiQT2d955Ub8mnBFYyDbpg4O5txArEAsKFNC+AYd103gkuGYvjcD45WVQ/xtIqphFI13gAUGgfelZUDwPIv/s87L8fB8QS+9OMTFZ/H8vkrBaTbFbYje/32HozYhuJcPFPwHFbjMBALOLvkdg5GA8Dlwx24crgD127qQkewunI4P5+Fapj4pcutOcXdYbnl/OJLRfAR3P2aEZyfz2LaNd3NzaGJRdz32BEA1YPRADDAYw6cctx+2SB+6Yoh/NMvzji7EEop7vj7X+ChPZYPPZ7WEPGLLT9rYKn8+us24+M370RvxO+oiLEi4zBl1zj0x/xOnKGdg9EA8IHXbcGjv/cG+HzEyW6rVFDIXEp3XLUBkkCwoY2D0W5Y24tKtTb3fvtl/Pv+Cdyye6Bm7VFv2A/RdnPymAOngOvsyVDMJTKTUvDq2AL+7icnYZoUCxm1ajC6XblxRy8+ecsuAHDyx4uNwwU7g6k/6ndaPrRrdXQ52PuqknJgwejdG2K4YVsPLm+DzrX1sL3f6nBwqkyHA0opJhdz+I3Xb8E/fuC6mmrA5yMYiAUgi76Wi4W1VmJuExhxLWy9ET9G7YKYsXgWPzsxg3iVvkoci6AsoDfix1i8MJ315dEFxAIihrtCTsyh3ZWDG6eJYYU+VadmUugJy+gMyfjGb1yPFlu7ls1gLICQLJStU0oqOhTdLJm9Uo2BmH/Z0/m8DFcOK2S42wrwjc1bCxubWBaQfPjn50Yxn9HQxeMNNdnUHSyJOTx7ehY3bOux01dZQJobB4Yk+BCWhYoB6RPTKceFIvhIy8w2XimEEGzrC5dVDmxuyFKMw1BHsCWn63HjsEKYcjg/by1so3MZCD6C37xxK358dBoHzi+UbdfNKWSkO4Sx+bxxGItnMBbP4vXbrTGg7MPa7tlKxXSG5JIOt5RSfOFHx/DS6HzJ3HKOxbbeSFnl4BiHJQzd+t23bMf/c8e6aSpdN9w4rJCwX0R3WHZcImfn0hjuCuITb9uF/+vWXQhIAnYOtPaMhkYw0hXC5GIOe8/G8amHXsEThy4AsGITAHDr7kH85fuuqJl33m50BKWSmMMjL4/jS0+dxK9eN4JPvG1nk67M22zvi2B8IYucVthnaTnK4dINHbhld+uNqeUxhwYw0hUsUA6be8KQRR9+76078Ttv2g4fl/M12dQdgmFSfPrfXsXoXAbfe2UcvRG/07ogKAv4letGavyW9sMyDoXZSv/8wih29kfwufdezl1JFdjWFwalwJnZdEG66nKMQ6vClUMDGO4OYSyeAaUUZ+fSzhBywOpC2q7985cCi92MzmVw085eUArcuKOHL2416AxJBW6l41NJ7Du3gF99zQi/d1XY7hrY5WYmpUASSNs0wawGVw4NYKQrhB8duoDZlIpkTsfmFp8DvRqw2E1QEvB3778Gx6aS2NrL72MtOkOFbqWH9oxBEgjuvHpjE6/K+2ztDYOQ0vY3M0kFfRE/N6zgxqEhjHQHoRkUL56JA0CBcuDUx1BHAFG/iDuv2YiOkMQDqXUSC0pYyGqglIIQgn9/dQJvvbgfPUsIqLYjQVnApu4Qjk0lCo7PJBXuUrLhxqEBsF3vMydnAIArh2UgCj48/vGbeB3DEukMylB1EznNhKIbmE4qbTFRsBFcMhjDEXtmA2MmqWCog2fEATzm0BBYb6Dv7RuH4CMY6W7v5mbLZaQ7BL/Yevniq4m7SvqMXWOztbfy/AFOnkuGYjg7l0bGNVN6JsWVA4MrhwawsTOIyzbG0BGU8GvXb+YLHGfNcPorZVWXceDKtR4uGYqCUuDohSSu2dQFw6SY48bB4f9v715jpDrrOI5/f+yyVij3i9KCIpVYsGqxRKm2JampEd9UUxPZGMH6ohetl3e2xsS+MbHENpVigqTFtGpaNNUIWIvYaLXewRLouikFREtLgE2RyxJF7N8X5xmcMrvA7Jzdcw7z+yQnM/PM2Wf+zz9n5z/nzJnnuDjkoKtzFBs/d23RYVgbmlg3hcbevn5GidMTGdrZ1U5h7d1/lK6OUYzuGMWr4dNYa1wczCpsQt01Hfb09TNr8hjPAHyeZk56PeMu6mRTzwHuXt9zei+smV9HX8i8FZlVWP3ke3/r6/chpSZIYt4bx/PrnYeIgL7j2Y8JveeQcXEwq7DaRaRq3zm4ODRn3oxsaps7rn8rH0xTYHj+rowPK5lV2NiuDsZ0dfDzngOcOPlf5rg4NGXJO2bw91dOcMt1czh56lWWPH/w9NmH7c7FwazCJHH74su4d/NOwKexNmvRnCksmpPN/DumCz66YGbBEZWHDyuZVdytiy9jbpqgcPZUf+q1fHjPwaziujpHsbJ7ARu3v8ylE/0DTMtHaYqDpA8B3wQ6gAcj4usFh2RWGfNmjH/N1NNmrSrFYSVJHcC3gCXAfKBb0vxiozIza1+lKA7Ae4BdEbEnIk4CjwE3FhyTmVnbKktxuBR4se7xvtT2GpJukbRF0pZDhw6NWHBmZu2mLMVhoCtrRENDxJqIWBgRC6dNmzYCYZmZtaeyFId9QP0FgmcCLxcUi5lZ2ytLcfgzMFfSWyR1AUuB9QXHZGbWtkpxKmtEnJJ0B7CJ7FTWtRHRU3BYZmZtqxTFASAingCeKDoOMzMDRTR871sJko4Bz6eHE4AjOXZf9v6mAn059pdnfGXPXU1eOazCeNtp+6tCf0Xlr/a6b46Ic5/RExGVXIAtdffX5Nx32fvbknN/ucVX9tzlncMqjLedtr+K9FdI/pp93bJ8Id2qDW3WX97yjM+5K1d/w9Vnnsqew3bLH1Dtw0pbImJh0XEUoZ3HnhfncOicu9YUlb9mX7fKew5rig6gQO089rw4h0Pn3LWmqPw19bqV3XMwM7PhU+U9BzMzGyYuDiUgaZakX0rqldQj6QupfbKkzZJeSLeTUvuUtP5xSavq+hknaVvd0ifp/qLGNZLyymF6rlvSDknbJT0paWoRYxopOefu4ylvPZJWFDGekTaE/N0gaWvaxrZKur6ur6tS+y5JKyUNNO/cyMjzlCovQz4VbQbw7nR/HLCT7LoWK4A7U/udwD3p/ljgGuA2YNVZ+t0KXFf0+KqUQ7Ifhh4EpqbHK4C7ix5fRXI3BfgHMC09fhj4QNHjK2H+FgCXpPtXAC/V9fUn4GqyyUh/BiwpalzecyiBiNgfEX9J948BvWRTlt9I9g9Guv1IWqc/Ip4B/jVYn5LmAtOB3wxj6KWRYw6VlrHpU9t4LvBJIHPM3RxgZ0TU5tP/BXDTMIdfuCHk79mIqG1TPcBFkl4naQYwPiJ+H1mleKT2N0VwcSgZSbPJPln8EXhDROyHbAMke7M/X93AurSRtZVWchgR/wFuB3aQFYX5wEPDGG6ptLj97QIulzRbUifZG9usc/zNBWUI+bsJeDYi/k1WUPbVPTfgdW1GiotDiUi6GHgc+GJEHG2xu6XAo61HVS2t5lDSaLLisAC4BNgO3JVrkCXVau4i4jBZ7taR7bHuBU7lGWOZNZs/SW8H7gFurTUNsFphH+5cHEoivSk9Dnw/In6Umg+kXU3S7cHz7OtdQGdEbB2WYEsqpxxeCRARu9Ne1w+A9w1TyKWR1/YXERsi4r0RcTXZ3GcvDFfMZdJs/iTNBH4MLIuI3al5H9m1bGoKva6Ni0MJpGPbDwG9EXFf3VPrgeXp/nLgJ+fZZTdttteQYw5fAuZLqk1MdgPZMeQLVp7bn6Tp6XYS8BngwXyjLZ9m8ydpIvBT4K6I+G1t5XTo6ZikRanPZZz//3z+iv6m30tAduZHkB3C2JaWD5Od/fEU2aevp4DJdX+zF3gFOE72iWN+3XN7gMuLHldVc0h2Fk5v6msDMKXo8VUod48Cf03L0qLHVsb8AV8B+uvW3QZMT88tBJ4DdgOrSD9ULmLxL6TNzKyBDyuZmVkDFwczM2vg4mBmZg1cHMzMrIGLg5mZNXBxMBsGkm6TtKyJ9WdLem44YzJrRmfRAZhdaCR1RsTqouMwa4WLg9kA0gRqT5JNoLaAbBrmZcA84D7gYqAP+FRE7Jf0K+B3wPuB9ZLGAccj4huSrgRWA2PIftz06Yg4LOkqYC1wAnhm5EZndm4+rGQ2uLcBayLincBR4LPAA8DHIqL2xv61uvUnRsTiiLj3jH4eAb6U+tkBfDW1fwf4fGTzEJmVivcczAb3Yvx/7pvvAV8muzjL5nSBrg5gf936687sQNIEsqLxdGp6GPjhAO3fBZbkPwSzoXFxMBvcmXPLHAN6zvJJv7+JvjVA/2al4cNKZoN7k6RaIegG/gBMq7VJGp3m5B9URBwBDku6NjV9Eng6Iv4JHJF0TWr/RP7hmw2d9xzMBtcLLJf0bbKZNR8ANgEr02GhTuB+sks9ns1yYLWkMWQz5t6c2m8G1ko6kfo1Kw3Pymo2gHS20saIuKLgUMwK4cNKZmbWwHsOZmbWwHsOZmbWwMXBzMwauDiYmVkDFwczM2vg4mBmZg1cHMzMrMH/AGaYKaljhBOXAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sorted_data['inc'][-200:].plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Etude de l'incidence annuelle"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Etant donné que le pic de l'épidémie se situe au printemps, nous définissons la période de référence\n",
"entre deux minima de l'incidence, du 1er septembre de l'année $N$ au\n",
"1er septembre de l'année $N+1$.\n",
"\n",
"Notre tâche est un peu compliquée par le fait que l'année ne comporte\n",
"pas un nombre entier de semaines. Nous modifions donc un peu nos périodes\n",
"de référence: à la place du 1er septembre de chaque année, nous utilisons le\n",
"premier jour de la semaine qui contient le 1er septembre.\n",
"\n",
"Comme l'incidence de la varicelle est très faible en automne, cette\n",
"modification ne risque pas de fausser nos conclusions.\n",
"\n",
"Encore un petit détail: les données commencent an décembre 1990, ce qui\n",
"rend la première année incomplète. Nous commençons donc l'analyse en 1991."
]
},
{
"cell_type": "code",
"execution_count": 13,
"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": "markdown",
"metadata": {},
"source": [
"En partant de cette liste des semaines qui contiennent un 1er septembre, 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": 14,
"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": 15,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fb3247501d0>"
]
},
"execution_count": 15,
"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": [
"<Figure size 432x288 with 1 Axes>"
]
},
"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": 16,
"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": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"yearly_incidence.sort_values()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Enfin, un histogramme montre bien que la force des épidémie est plutôt répartie. Il y a quand même plus d'épidémie de force moyenne."
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fb2ec0d50b8>"
]
},
"execution_count": 17,
"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": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"yearly_incidence.hist(xrot=20)"
]
},
{
"cell_type": "code",
"execution_count": null,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment