diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb
index 0bbbe371b01e359e381e43239412d77bf53fb1fb..c9e6511e17c54607680a3e078b05f3a83b59d8b9 100644
--- a/module3/exo2/exercice.ipynb
+++ b/module3/exo2/exercice.ipynb
@@ -1,5 +1,1348 @@
{
- "cells": [],
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Analyse de l'incidence de la varicelle"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "%matplotlib inline\n",
+ "import matplotlib.pyplot as plt\n",
+ "import pandas as pd\n",
+ "import isoweek"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Les données sont récupérées sur le site web du réseau sentinelle. Les données recouvrent la période de 1991 à nos jours."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "data_url = \"https://www.sentiweb.fr/datasets/incidence-PAY-7.csv\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "La première ligne est un commentaire, elle est skip"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "
\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " | \n",
+ " week | \n",
+ " indicator | \n",
+ " inc | \n",
+ " inc_low | \n",
+ " inc_up | \n",
+ " inc100 | \n",
+ " inc100_low | \n",
+ " inc100_up | \n",
+ " geo_insee | \n",
+ " geo_name | \n",
+ "
\n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 | \n",
+ " 202041 | \n",
+ " 7 | \n",
+ " 4094 | \n",
+ " 1977 | \n",
+ " 6211 | \n",
+ " 6 | \n",
+ " 3 | \n",
+ " 9 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1 | \n",
+ " 202040 | \n",
+ " 7 | \n",
+ " 2090 | \n",
+ " 662 | \n",
+ " 3518 | \n",
+ " 3 | \n",
+ " 1 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 2 | \n",
+ " 202039 | \n",
+ " 7 | \n",
+ " 1058 | \n",
+ " 240 | \n",
+ " 1876 | \n",
+ " 2 | \n",
+ " 1 | \n",
+ " 3 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 3 | \n",
+ " 202038 | \n",
+ " 7 | \n",
+ " 2253 | \n",
+ " 782 | \n",
+ " 3724 | \n",
+ " 3 | \n",
+ " 1 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 4 | \n",
+ " 202037 | \n",
+ " 7 | \n",
+ " 1584 | \n",
+ " 405 | \n",
+ " 2763 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 4 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 5 | \n",
+ " 202036 | \n",
+ " 7 | \n",
+ " 919 | \n",
+ " 100 | \n",
+ " 1738 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 6 | \n",
+ " 202035 | \n",
+ " 7 | \n",
+ " 828 | \n",
+ " 0 | \n",
+ " 1694 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 7 | \n",
+ " 202034 | \n",
+ " 7 | \n",
+ " 2272 | \n",
+ " 371 | \n",
+ " 4173 | \n",
+ " 3 | \n",
+ " 0 | \n",
+ " 6 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 8 | \n",
+ " 202033 | \n",
+ " 7 | \n",
+ " 1284 | \n",
+ " 177 | \n",
+ " 2391 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 4 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 9 | \n",
+ " 202032 | \n",
+ " 7 | \n",
+ " 2650 | \n",
+ " 689 | \n",
+ " 4611 | \n",
+ " 4 | \n",
+ " 1 | \n",
+ " 7 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 10 | \n",
+ " 202031 | \n",
+ " 7 | \n",
+ " 1303 | \n",
+ " 100 | \n",
+ " 2506 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 4 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 11 | \n",
+ " 202030 | \n",
+ " 7 | \n",
+ " 1385 | \n",
+ " 75 | \n",
+ " 2695 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 4 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 12 | \n",
+ " 202029 | \n",
+ " 7 | \n",
+ " 841 | \n",
+ " 10 | \n",
+ " 1672 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 13 | \n",
+ " 202028 | \n",
+ " 7 | \n",
+ " 728 | \n",
+ " 0 | \n",
+ " 1515 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 14 | \n",
+ " 202027 | \n",
+ " 7 | \n",
+ " 986 | \n",
+ " 149 | \n",
+ " 1823 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 15 | \n",
+ " 202026 | \n",
+ " 7 | \n",
+ " 694 | \n",
+ " 0 | \n",
+ " 1454 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 16 | \n",
+ " 202025 | \n",
+ " 7 | \n",
+ " 228 | \n",
+ " 0 | \n",
+ " 597 | \n",
+ " 0 | \n",
+ " 0 | \n",
+ " 1 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 17 | \n",
+ " 202024 | \n",
+ " 7 | \n",
+ " 388 | \n",
+ " 0 | \n",
+ " 959 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 18 | \n",
+ " 202023 | \n",
+ " 7 | \n",
+ " 558 | \n",
+ " 1 | \n",
+ " 1115 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 19 | \n",
+ " 202022 | \n",
+ " 7 | \n",
+ " 277 | \n",
+ " 0 | \n",
+ " 633 | \n",
+ " 0 | \n",
+ " 0 | \n",
+ " 1 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 20 | \n",
+ " 202021 | \n",
+ " 7 | \n",
+ " 602 | \n",
+ " 36 | \n",
+ " 1168 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 21 | \n",
+ " 202020 | \n",
+ " 7 | \n",
+ " 824 | \n",
+ " 20 | \n",
+ " 1628 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 22 | \n",
+ " 202019 | \n",
+ " 7 | \n",
+ " 310 | \n",
+ " 0 | \n",
+ " 753 | \n",
+ " 0 | \n",
+ " 0 | \n",
+ " 1 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 23 | \n",
+ " 202018 | \n",
+ " 7 | \n",
+ " 849 | \n",
+ " 98 | \n",
+ " 1600 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 24 | \n",
+ " 202017 | \n",
+ " 7 | \n",
+ " 272 | \n",
+ " 0 | \n",
+ " 658 | \n",
+ " 0 | \n",
+ " 0 | \n",
+ " 1 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 25 | \n",
+ " 202016 | \n",
+ " 7 | \n",
+ " 758 | \n",
+ " 78 | \n",
+ " 1438 | \n",
+ " 1 | \n",
+ " 0 | \n",
+ " 2 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 26 | \n",
+ " 202015 | \n",
+ " 7 | \n",
+ " 1918 | \n",
+ " 675 | \n",
+ " 3161 | \n",
+ " 3 | \n",
+ " 1 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 27 | \n",
+ " 202014 | \n",
+ " 7 | \n",
+ " 3879 | \n",
+ " 2227 | \n",
+ " 5531 | \n",
+ " 6 | \n",
+ " 3 | \n",
+ " 9 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 28 | \n",
+ " 202013 | \n",
+ " 7 | \n",
+ " 7326 | \n",
+ " 5236 | \n",
+ " 9416 | \n",
+ " 11 | \n",
+ " 8 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 29 | \n",
+ " 202012 | \n",
+ " 7 | \n",
+ " 8123 | \n",
+ " 5790 | \n",
+ " 10456 | \n",
+ " 12 | \n",
+ " 8 | \n",
+ " 16 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ "
\n",
+ " \n",
+ " 1528 | \n",
+ " 199126 | \n",
+ " 7 | \n",
+ " 17608 | \n",
+ " 11304 | \n",
+ " 23912 | \n",
+ " 31 | \n",
+ " 20 | \n",
+ " 42 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1529 | \n",
+ " 199125 | \n",
+ " 7 | \n",
+ " 16169 | \n",
+ " 10700 | \n",
+ " 21638 | \n",
+ " 28 | \n",
+ " 18 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1530 | \n",
+ " 199124 | \n",
+ " 7 | \n",
+ " 16171 | \n",
+ " 10071 | \n",
+ " 22271 | \n",
+ " 28 | \n",
+ " 17 | \n",
+ " 39 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1531 | \n",
+ " 199123 | \n",
+ " 7 | \n",
+ " 11947 | \n",
+ " 7671 | \n",
+ " 16223 | \n",
+ " 21 | \n",
+ " 13 | \n",
+ " 29 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1532 | \n",
+ " 199122 | \n",
+ " 7 | \n",
+ " 15452 | \n",
+ " 9953 | \n",
+ " 20951 | \n",
+ " 27 | \n",
+ " 17 | \n",
+ " 37 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1533 | \n",
+ " 199121 | \n",
+ " 7 | \n",
+ " 14903 | \n",
+ " 8975 | \n",
+ " 20831 | \n",
+ " 26 | \n",
+ " 16 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1534 | \n",
+ " 199120 | \n",
+ " 7 | \n",
+ " 19053 | \n",
+ " 12742 | \n",
+ " 25364 | \n",
+ " 34 | \n",
+ " 23 | \n",
+ " 45 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1535 | \n",
+ " 199119 | \n",
+ " 7 | \n",
+ " 16739 | \n",
+ " 11246 | \n",
+ " 22232 | \n",
+ " 29 | \n",
+ " 19 | \n",
+ " 39 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1536 | \n",
+ " 199118 | \n",
+ " 7 | \n",
+ " 21385 | \n",
+ " 13882 | \n",
+ " 28888 | \n",
+ " 38 | \n",
+ " 25 | \n",
+ " 51 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1537 | \n",
+ " 199117 | \n",
+ " 7 | \n",
+ " 13462 | \n",
+ " 8877 | \n",
+ " 18047 | \n",
+ " 24 | \n",
+ " 16 | \n",
+ " 32 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1538 | \n",
+ " 199116 | \n",
+ " 7 | \n",
+ " 14857 | \n",
+ " 10068 | \n",
+ " 19646 | \n",
+ " 26 | \n",
+ " 18 | \n",
+ " 34 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1539 | \n",
+ " 199115 | \n",
+ " 7 | \n",
+ " 13975 | \n",
+ " 9781 | \n",
+ " 18169 | \n",
+ " 25 | \n",
+ " 18 | \n",
+ " 32 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1540 | \n",
+ " 199114 | \n",
+ " 7 | \n",
+ " 12265 | \n",
+ " 7684 | \n",
+ " 16846 | \n",
+ " 22 | \n",
+ " 14 | \n",
+ " 30 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1541 | \n",
+ " 199113 | \n",
+ " 7 | \n",
+ " 9567 | \n",
+ " 6041 | \n",
+ " 13093 | \n",
+ " 17 | \n",
+ " 11 | \n",
+ " 23 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1542 | \n",
+ " 199112 | \n",
+ " 7 | \n",
+ " 10864 | \n",
+ " 7331 | \n",
+ " 14397 | \n",
+ " 19 | \n",
+ " 13 | \n",
+ " 25 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1543 | \n",
+ " 199111 | \n",
+ " 7 | \n",
+ " 15574 | \n",
+ " 11184 | \n",
+ " 19964 | \n",
+ " 27 | \n",
+ " 19 | \n",
+ " 35 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1544 | \n",
+ " 199110 | \n",
+ " 7 | \n",
+ " 16643 | \n",
+ " 11372 | \n",
+ " 21914 | \n",
+ " 29 | \n",
+ " 20 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1545 | \n",
+ " 199109 | \n",
+ " 7 | \n",
+ " 13741 | \n",
+ " 8780 | \n",
+ " 18702 | \n",
+ " 24 | \n",
+ " 15 | \n",
+ " 33 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1546 | \n",
+ " 199108 | \n",
+ " 7 | \n",
+ " 13289 | \n",
+ " 8813 | \n",
+ " 17765 | \n",
+ " 23 | \n",
+ " 15 | \n",
+ " 31 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1547 | \n",
+ " 199107 | \n",
+ " 7 | \n",
+ " 12337 | \n",
+ " 8077 | \n",
+ " 16597 | \n",
+ " 22 | \n",
+ " 15 | \n",
+ " 29 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1548 | \n",
+ " 199106 | \n",
+ " 7 | \n",
+ " 10877 | \n",
+ " 7013 | \n",
+ " 14741 | \n",
+ " 19 | \n",
+ " 12 | \n",
+ " 26 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1549 | \n",
+ " 199105 | \n",
+ " 7 | \n",
+ " 10442 | \n",
+ " 6544 | \n",
+ " 14340 | \n",
+ " 18 | \n",
+ " 11 | \n",
+ " 25 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1550 | \n",
+ " 199104 | \n",
+ " 7 | \n",
+ " 7913 | \n",
+ " 4563 | \n",
+ " 11263 | \n",
+ " 14 | \n",
+ " 8 | \n",
+ " 20 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1551 | \n",
+ " 199103 | \n",
+ " 7 | \n",
+ " 15387 | \n",
+ " 10484 | \n",
+ " 20290 | \n",
+ " 27 | \n",
+ " 18 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1552 | \n",
+ " 199102 | \n",
+ " 7 | \n",
+ " 16277 | \n",
+ " 11046 | \n",
+ " 21508 | \n",
+ " 29 | \n",
+ " 20 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1553 | \n",
+ " 199101 | \n",
+ " 7 | \n",
+ " 15565 | \n",
+ " 10271 | \n",
+ " 20859 | \n",
+ " 27 | \n",
+ " 18 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1554 | \n",
+ " 199052 | \n",
+ " 7 | \n",
+ " 19375 | \n",
+ " 13295 | \n",
+ " 25455 | \n",
+ " 34 | \n",
+ " 23 | \n",
+ " 45 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1555 | \n",
+ " 199051 | \n",
+ " 7 | \n",
+ " 19080 | \n",
+ " 13807 | \n",
+ " 24353 | \n",
+ " 34 | \n",
+ " 25 | \n",
+ " 43 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1556 | \n",
+ " 199050 | \n",
+ " 7 | \n",
+ " 11079 | \n",
+ " 6660 | \n",
+ " 15498 | \n",
+ " 20 | \n",
+ " 12 | \n",
+ " 28 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1557 | \n",
+ " 199049 | \n",
+ " 7 | \n",
+ " 1143 | \n",
+ " 0 | \n",
+ " 2610 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "
1558 rows × 10 columns
\n",
+ "
"
+ ],
+ "text/plain": [
+ " week indicator inc inc_low inc_up inc100 inc100_low \\\n",
+ "0 202041 7 4094 1977 6211 6 3 \n",
+ "1 202040 7 2090 662 3518 3 1 \n",
+ "2 202039 7 1058 240 1876 2 1 \n",
+ "3 202038 7 2253 782 3724 3 1 \n",
+ "4 202037 7 1584 405 2763 2 0 \n",
+ "5 202036 7 919 100 1738 1 0 \n",
+ "6 202035 7 828 0 1694 1 0 \n",
+ "7 202034 7 2272 371 4173 3 0 \n",
+ "8 202033 7 1284 177 2391 2 0 \n",
+ "9 202032 7 2650 689 4611 4 1 \n",
+ "10 202031 7 1303 100 2506 2 0 \n",
+ "11 202030 7 1385 75 2695 2 0 \n",
+ "12 202029 7 841 10 1672 1 0 \n",
+ "13 202028 7 728 0 1515 1 0 \n",
+ "14 202027 7 986 149 1823 1 0 \n",
+ "15 202026 7 694 0 1454 1 0 \n",
+ "16 202025 7 228 0 597 0 0 \n",
+ "17 202024 7 388 0 959 1 0 \n",
+ "18 202023 7 558 1 1115 1 0 \n",
+ "19 202022 7 277 0 633 0 0 \n",
+ "20 202021 7 602 36 1168 1 0 \n",
+ "21 202020 7 824 20 1628 1 0 \n",
+ "22 202019 7 310 0 753 0 0 \n",
+ "23 202018 7 849 98 1600 1 0 \n",
+ "24 202017 7 272 0 658 0 0 \n",
+ "25 202016 7 758 78 1438 1 0 \n",
+ "26 202015 7 1918 675 3161 3 1 \n",
+ "27 202014 7 3879 2227 5531 6 3 \n",
+ "28 202013 7 7326 5236 9416 11 8 \n",
+ "29 202012 7 8123 5790 10456 12 8 \n",
+ "... ... ... ... ... ... ... ... \n",
+ "1528 199126 7 17608 11304 23912 31 20 \n",
+ "1529 199125 7 16169 10700 21638 28 18 \n",
+ "1530 199124 7 16171 10071 22271 28 17 \n",
+ "1531 199123 7 11947 7671 16223 21 13 \n",
+ "1532 199122 7 15452 9953 20951 27 17 \n",
+ "1533 199121 7 14903 8975 20831 26 16 \n",
+ "1534 199120 7 19053 12742 25364 34 23 \n",
+ "1535 199119 7 16739 11246 22232 29 19 \n",
+ "1536 199118 7 21385 13882 28888 38 25 \n",
+ "1537 199117 7 13462 8877 18047 24 16 \n",
+ "1538 199116 7 14857 10068 19646 26 18 \n",
+ "1539 199115 7 13975 9781 18169 25 18 \n",
+ "1540 199114 7 12265 7684 16846 22 14 \n",
+ "1541 199113 7 9567 6041 13093 17 11 \n",
+ "1542 199112 7 10864 7331 14397 19 13 \n",
+ "1543 199111 7 15574 11184 19964 27 19 \n",
+ "1544 199110 7 16643 11372 21914 29 20 \n",
+ "1545 199109 7 13741 8780 18702 24 15 \n",
+ "1546 199108 7 13289 8813 17765 23 15 \n",
+ "1547 199107 7 12337 8077 16597 22 15 \n",
+ "1548 199106 7 10877 7013 14741 19 12 \n",
+ "1549 199105 7 10442 6544 14340 18 11 \n",
+ "1550 199104 7 7913 4563 11263 14 8 \n",
+ "1551 199103 7 15387 10484 20290 27 18 \n",
+ "1552 199102 7 16277 11046 21508 29 20 \n",
+ "1553 199101 7 15565 10271 20859 27 18 \n",
+ "1554 199052 7 19375 13295 25455 34 23 \n",
+ "1555 199051 7 19080 13807 24353 34 25 \n",
+ "1556 199050 7 11079 6660 15498 20 12 \n",
+ "1557 199049 7 1143 0 2610 2 0 \n",
+ "\n",
+ " inc100_up geo_insee geo_name \n",
+ "0 9 FR France \n",
+ "1 5 FR France \n",
+ "2 3 FR France \n",
+ "3 5 FR France \n",
+ "4 4 FR France \n",
+ "5 2 FR France \n",
+ "6 2 FR France \n",
+ "7 6 FR France \n",
+ "8 4 FR France \n",
+ "9 7 FR France \n",
+ "10 4 FR France \n",
+ "11 4 FR France \n",
+ "12 2 FR France \n",
+ "13 2 FR France \n",
+ "14 2 FR France \n",
+ "15 2 FR France \n",
+ "16 1 FR France \n",
+ "17 2 FR France \n",
+ "18 2 FR France \n",
+ "19 1 FR France \n",
+ "20 2 FR France \n",
+ "21 2 FR France \n",
+ "22 1 FR France \n",
+ "23 2 FR France \n",
+ "24 1 FR France \n",
+ "25 2 FR France \n",
+ "26 5 FR France \n",
+ "27 9 FR France \n",
+ "28 14 FR France \n",
+ "29 16 FR France \n",
+ "... ... ... ... \n",
+ "1528 42 FR France \n",
+ "1529 38 FR France \n",
+ "1530 39 FR France \n",
+ "1531 29 FR France \n",
+ "1532 37 FR France \n",
+ "1533 36 FR France \n",
+ "1534 45 FR France \n",
+ "1535 39 FR France \n",
+ "1536 51 FR France \n",
+ "1537 32 FR France \n",
+ "1538 34 FR France \n",
+ "1539 32 FR France \n",
+ "1540 30 FR France \n",
+ "1541 23 FR France \n",
+ "1542 25 FR France \n",
+ "1543 35 FR France \n",
+ "1544 38 FR France \n",
+ "1545 33 FR France \n",
+ "1546 31 FR France \n",
+ "1547 29 FR France \n",
+ "1548 26 FR France \n",
+ "1549 25 FR France \n",
+ "1550 20 FR France \n",
+ "1551 36 FR France \n",
+ "1552 38 FR France \n",
+ "1553 36 FR France \n",
+ "1554 45 FR France \n",
+ "1555 43 FR France \n",
+ "1556 28 FR France \n",
+ "1557 5 FR France \n",
+ "\n",
+ "[1558 rows x 10 columns]"
+ ]
+ },
+ "execution_count": 5,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "raw_data = pd.read_csv(data_url, skiprows=1)\n",
+ "raw_data"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "On vérifie qu'il n'y a pas de données manquantes"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " | \n",
+ " week | \n",
+ " indicator | \n",
+ " inc | \n",
+ " inc_low | \n",
+ " inc_up | \n",
+ " inc100 | \n",
+ " inc100_low | \n",
+ " inc100_up | \n",
+ " geo_insee | \n",
+ " geo_name | \n",
+ "
\n",
+ " \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
"
+ ],
+ "text/plain": [
+ "Empty DataFrame\n",
+ "Columns: [week, indicator, inc, inc_low, inc_up, inc100, inc100_low, inc100_up, geo_insee, geo_name]\n",
+ "Index: []"
+ ]
+ },
+ "execution_count": 6,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "raw_data[raw_data.isnull().any(axis=1)]"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Aucune donnée n'est manquante, donc pas besoin de retirer des lignes."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Comme pour le cas de la grippe, il y a besoin de reformater les données qui ne sont pas dans un format standard interpretable par pandas. "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# ce code est repris tel quel de l'exemple du cours.\n",
+ "def convert_week(year_and_week_int):\n",
+ " year_and_week_str = str(year_and_week_int)\n",
+ " year = int(year_and_week_str[:4])\n",
+ " week = int(year_and_week_str[4:])\n",
+ " w = isoweek.Week(year, week)\n",
+ " return pd.Period(w.day(0), 'W')\n",
+ "\n",
+ "raw_data['period'] = [convert_week(yw) for yw in raw_data['week']]"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "De même, nous réorganisons le jeu de données pour que l'index soit les dates, et qu'il soit trié par ordre chronologique."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "sorted_data = raw_data.set_index('period').sort_index()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "On visualise les données, avec un zoom sur les dernières années:"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 13,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXmcHUW59vOec2Ym6yQhOwmQkIQt7ER2EAhLwAVc0LjivXhRLnrdrvcDlSsqKIqIooCioIAXAVEERQSEgIKBJCwBAoQEEkhISEL2bSYzc+r7o7u6q6vfWvpMz8yZTD/55Td9qqurq6ur6613JyEEChQoUKBAARWlnu5AgQIFChSoPxTEoUCBAgUKpFAQhwIFChQokEJBHAoUKFCgQAoFcShQoECBAikUxKFAgQIFCqRQEIcCBQoUKJBCQRwKFChQoEAKBXEoUKBAgQIpVHq6A7VixIgRYsKECT3djQIFChToVXjqqafeFkKMdNXrtcRhwoQJmDdvXk93o0CBAgV6FYjodZ96hVipQIECBQqkUBCHAgUKFCiQQkEcChQoUKBACgVxKFCgQIECKRTEoUCBAgUKpFAQhwIFChQokEJBHAoUKFCgQAoFcSiw0+PRV9Zg2bptPd2NAgV6FXqtE1yBAr4458Y5aCgTFl12Rk93pUCBXoOCcyjQJ9DWIXq6CwUK9CoUxKFAgQIFCqRQEIcCBQoUKJCCkzgQUT8imkNE84loARF9Kyy/hIjeJKJnw/9nKNdcRESLiWghEZ2mlB9GRM+H564mIgrLm4jo9rD8SSKakP+jFihQoEABX/hwDq0AThJCHATgYAAziOjI8NxVQoiDw/9/BQAi2g/ATABTAcwAcC0RlcP61wE4D8CU8P+MsPxcAOuFEJMBXAXg+51/tAIFChQoUCucxEEE2BL+bAj/27R7ZwK4TQjRKoRYAmAxgMOJaCyAZiHEbCGEAHAzgLOUa24Kj+8EMF1yFQUKFChQoPvhpXMgojIRPQtgNYAHhRBPhqc+R0TPEdGNRDQsLBsHYJly+fKwbFx4rJcnrhFCtAPYCGA404/ziGgeEc1bs2aN1wMWKFCgQIHs8CIOQogOIcTBAMYj4AL2RyAimoRA1LQSwJVhdW7HLyzltmv0flwvhJgmhJg2cqQzkVGBAp3Ckre3Yv6yDT3djQIFegSZrJWEEBsAPAJghhBiVUg0qgB+CeDwsNpyALspl40HsCIsH8+UJ64hogqAIQDWZXqSAgVyxok/fARnXvN4T3ejQIEegY+10kgiGhoe9wdwMoCXQx2CxPsAvBAe3wNgZmiBNBGB4nmOEGIlgM1EdGSoT/gkgLuVa84Jjz8I4OFQL1GgQIECBXoAPuEzxgK4KbQ4KgG4QwjxFyK6hYgORiD+WQrgMwAghFhARHcAeBFAO4ALhBAdYVvnA/gNgP4A7gv/A8ANAG4hosUIOIaZOTxbgQIFChSoEU7iIIR4DsAhTPknLNdcBuAypnwegP2Z8hYAZ7v6UqBAgQIFugeFh3SBPoWL//QCrpm1uKe7UaBA3aOIylqgT+GWJ14HAFxw4uQe7kmBAvWNgnMoUKBAgQIpFMShQIECBQqkUBCHAjs1eotF9FOvr8c3736h1/S3wM6PgjgU2KnRW9bas3/+L9w0+3VUe0l/C+z8KIhDgZ0axVpboEBtKIhDgZ0avU1M09v6W2DnRUEcCuzUKJbaAgVqQ0EcCuzU6G0b8V7W3QI7MQriUGCnhuhly21vI2YFdl4UxKEX4Y2127D07a093Y1ehd622PY2YlZg50URPqMX4fgrZgEAll7+rh7uSYECBXZ2FJxDgZ0avY5z6GX9LbDzoiAOBXZqFGKaAgVqQ0EcCuzU6C078V7SzZ0K9z63EotWbe7pbtQtCuJQYKdGb1t0ewsx6+1obe/ABbc+jbN/Mbunu1K3KIhDgZ0a1V6y2lL4txCDdQ+q1eDvhm1tPduROoaTOBBRPyKaQ0TziWgBEX0rLN+FiB4kokXh32HKNRcR0WIiWkhEpynlhxHR8+G5q4mIwvImIro9LH+SiCbk/6gF+iJ6CW2I0N39/f28ZZhw4b1Yt3VH9964h1EQYTd8OIdWACcJIQ4CcDCAGUR0JIALATwkhJgC4KHwN4hoPwAzAUwFMAPAtURUDtu6DsB5AKaE/2eE5ecCWC+EmAzgKgDfz+HZckW1KnDz7KVoaevo6a4UyIJesgb0VDd/G2bGe31t3/Kf6W2bhp6AkziIAFvCnw3hfwHgTAA3heU3ATgrPD4TwG1CiFYhxBIAiwEcTkRjATQLIWaLILrYzdo1sq07AUyXXEW94G8L3sL/3r0AP7x/YU93pUAG9LYdYk/1ts4+ty5H75oVPQMvnQMRlYnoWQCrATwohHgSwGghxEoACP+OCquPA7BMuXx5WDYuPNbLE9cIIdoBbAQwvJYH6ipsaW0HAGzYXsgoexN6yw4x0jl0c4d7yfAU6AF4EQchRIcQ4mAA4xFwAftbqnNbEGEpt12TbJjoPCKaR0Tz1qxZ4+p2gQK9ZvHrqX5KWtS3+IYiNLoPMlkrCSE2AHgEga5gVSgqQvh3dVhtOYDdlMvGA1gRlo9nyhPXEFEFwBAA65j7Xy+EmCaEmDZy5MgsXS/QC7C1tR1rt7Tm2mZvWwTy7u2yddvw5GtrnfX6mFSp12waehI+1kojiWhoeNwfwMkAXgZwD4BzwmrnALg7PL4HwMzQAmkiAsXznFD0tJmIjgz1CZ/UrpFtfRDAw6JOv+rX127F9+57qdctOr0Bp171Dxx26d+967e0deDsn/8LC1ZsNNbpbW8p72l13A9m4cPXP2G+X68boXxQfL5u+HAOYwHMIqLnAMxFoHP4C4DLAZxCRIsAnBL+hhBiAYA7ALwI4G8ALhBCSBOf8wH8CoGS+lUA94XlNwAYTkSLAXwZoeVTPWLu0vX4xaOvYUlO0VG3trbj5bc25dJWb8ebG7Znqv/0G+sxd+l6fOcvLxrr9JZFINq459TfJ15bi7aOqrNeLFYqWIcCSTijsgohngNwCFO+FsB0wzWXAbiMKZ8HIKWvEEK0ADjbo791g7zm1nm3zMPji9di8WWno1IufBLzRm/ZGefZy+eXb8TM65/Ap4+d6H1N3xMr9Y550ZMoVqMakdeOdO6S9QCAjt6yxa1DWIeulw1rHovW2q2B3mahR9ygvjrt+upzZ0FBHGpERzWn2RXu2IrJ2jXobcOaxzwohWyAT1u9bXwKdB8K4lAjWtvz8ZSO7dtzaa6Ahr44rlJEpHIhJgMKWS6veX75RnznLy/u9AYXO/fT5YOCONSIHe1uZZ8P5C6vtwSIqyf4KFF7m2w5j95ynENbR7rlNzdsx/YwHIwcyw9c9y/c8NgStOY0v+sVOzvxywNFmtAakdfHUwrXt4I4dA1627DmsWhJkqnOqfZqFY3aXvCYyx+Or9HmYWkn11D3smnRIyg4hxqRl84h5hzs9fISY/U19LZFIJf+MnosjnPg0FcMI/rIY3YKBXGoEbnNrehDtrc4b+n6vO6408FqrNQHVwFerGTndEkjKL1NHJcVO/vz5YGCONSIvBadWASQS3MFNPQ22pCntZIqVqo6Jpiuv+lt41YgfxTEoUbk9e2USoVCut7RndxHHjvaUmStFKOviIu8UQyHEwVxqBU5Ta5ol1ewDpnhozNV18Rax9i1rk679O+4+E8vJMr+8coavJ1zEEFf6Mrl4Nh+TV8Ts/Stp60NBXGoEXl9TLG1Ui7N9Sn4OXnFlf779/Nru4/j/NtbWnFLmFENCIjQJ2+cg4/+0hzwruabITCG+NGDr2DDtnRqz7VbWvGB62YDANYrqT9dhFEfy52d0djZny8PFMShRuQ3uQqxUldCtSr74zNv1tRGVrGSrL149RZrPdu1Nsx6eTWufmgRvnnPgtS515SAkEvXbouOXfOrr02/vsYp1YKCONSIvD6mws/BjhdXbMKtT77BnvMRK7VnYMlWb27Bc8s3AEjuumt9M7Vc5zMNdoSWR61taQsk05C4TK/1xXJnXzyLz82NgjjUiLzmlm5CWCCJM67+J7521/M1X9/uad8PAGf85DG892ePAwA+9Zu5UXnWd9NZBfaP//4K5i5N5bpK9YcjjqZc0C4a2dLWkRA97ezzcSd/vFxQEIcakZcFi1RI5xbIrxdg1sLVuPvZ2kQ8WdFe9fdkVxXI85dtiI6z7qI78yYFBH7890U4++ezrXUAnjiUiP/t4kw/cN1sfOvPaTFVgb6LInxGjcjNlLUPxlb6t18Hu/IzDx6XT4OWofP1DFahE+rsnEPmW2a6NuYc0tRBLysRoSqE1/y6aXasVN/ZZ6O6uXvitbXYsK0NM/Yf04M9qj8UxKFG5L2W9yHGoVtRC0fW0tb9oUqyzCe50HMCJL1M0oqs47Cze5arjzczTKO69PJ39VBv6hOFWKlm5CRWCt/Azv4x9hTaPVJl6tiuEYfMnEMn5sYzb2xwVwrBcQ56wDzKkNuhQAEVTuJARLsR0SwieomIFhDRF8LyS4joTSJ6Nvx/hnLNRUS0mIgWEtFpSvlhRPR8eO5qCmcuETUR0e1h+ZNENCH/R80XeX1sFJmy5tNeX4JP3NC2GgZ2+w6NOGTVOdTwLuWafsGtTzvrynDxun5BbUeiVCvnkKm2Gas2teCtjS05tZYfCmLphg/n0A7gK0KIfQEcCeACItovPHeVEOLg8P9fASA8NxPAVAAzAFxLROWw/nUAzgMwJfw/Iyw/F8B6IcRkAFcB+H7nH61rkbe1Ul/SObiQJxeVhXNoKAcvQxcr1dqdrnqlX73zOQB+xLFWnVZefT/iuw/hyO89lE9jOWJnN9XNA07iIIRYKYR4OjzeDOAlADZN4pkAbhNCtAohlgBYDOBwIhoLoFkIMVsEX//NAM5SrrkpPL4TwHTieOY6Qn5+Dn3PWsmFPIcii59DUyXYw+hipe4g3LVxGx5ipfBvvW4+Nre04eQfPYrnl2/s1vvW6XDUFTLpHEJxzyEAngyLPkdEzxHRjUQ0LCwbB2CZctnysGxceKyXJ64RQrQD2AhgeJa+dTfy2nn4+jn0pcmc3SPZXN+H6AohcOuTb0TOZWmxUjZ017vidk8lwxetD4NzjLvpGeYuXYfFq7fgygcXds8NC3jDmzgQ0SAAfwDwRSHEJgQiokkADgawEsCVsipzubCU267R+3AeEc0jonlr1qzx7XqXIK8FoBxSB5c9fr2xwdWqwC8efRWbWtrybzvHR/XZMc9Zsg5fu+v5SJbfEwrpWvhk1pTVIGz6xl0vYLPyrty0oXvm24724D4VE1XrItTX11Sf8HojRNSAgDD8nxDijwAghFglhOgQQlQB/BLA4WH15QB2Uy4fD2BFWD6eKU9cQ0QVAEMApFxEhRDXCyGmCSGmjRw50u8Juwh5Ta6GcvAKXPb49SZ1eujl1fjefS/ju/e+lHvb3IJeqx7C57JtGjG4+9kVyQrdMPa1PB63gOtKaklAFq7ajOseebWWrnUp5KbojXVbcd/zK7vtvoV1oBs+1koE4AYALwkhfqSUj1WqvQ+AjFl8D4CZoQXSRASK5zlCiJUANhPRkWGbnwRwt3LNOeHxBwE8LPrI22uoBK9ghyMndU/KjO96ZjnOu3leQkQjd6Hd5RPAPb6PWqqWUbtLC9DXHdZKtYCLtKoPifoZqTkdXF3srmeQ4U1eWbUF5/+f21IrD9z42BJ84oY53XKv3gwfJ7hjAHwCwPNE9GxY9jUAHyGigxHMs6UAPgMAQogFRHQHgBcRWDpdIISQK8j5AH4DoD+A+8L/QEB8biGixQg4hpmde6zs2NrajideW4vp+472qp8X7WoMLWR2dNgX2Z6klV+6PQh1vaWlHUMGNACIiVljJX9xAMs5MPV8xiSPccsuVuoedIabdI1Ldz3Djhr8UDqLb//lxW6/Z2+EkzgIIR4DrxP4q+WaywBcxpTPA7A/U94C4GxXX7oS37xnAe58ajn+9sXjsM+Y5m67b6Mv59D931AK6s5T5iSWYrE8wS16VSFQ9jLezA5Xq/XKwmblJk36iJ5ElsCIAPDWxhZUyoQRg5q6qEcFJAoP6RCrNgWOOr4OO3lt5OXi2uogDvWQ5lFdjGR/pfknEPgU3PXM8k5nteN2tTWLlXIYtszWU930rnjdTPL3VsXyShWPucVK3fMMWQIjAsCR33sI0y79exf1poCKgjiEaO4XiEs2tbR71c/LmqOx7Mc51IMKhiMOqljphseW4Eu3z685qU58n3RZreOdx3vKbMra6Tv6gTPT9b13HUwnAPXl31NPfakHFMQhhFzkVI/a7/zlRUy48F62fm6cgxQrOWSv9TBv1U2eFCtJnQkArNkchLxWE+XUBI449KCXcndGZe3K+2QRK3XXdKsnQVe9Ogr2FAriECKKdKnM1hseW2Ksn9c8avLkHOph4qp9sIWNtuG3Sq5ln/vo9+NQ6zlf1JuPiQTLYXl21fVM3ebIV0eBEArOIYmCOISQ80IPP2BCZ6dRtSrw5/krUA4N093EoZM3rBFqApjOfjwrN27HN/70grMeb61Uq1gpB9SpXCnrOKlTO4/F/9FX1uANJU91b0c9bMDqCUU+hxAi4hzSxGHjtjZs2t7G1q8Vv5v7Br5+V7xQ1qvO4dePL1X6wNdp76hi1eZW/mSint8zcLVqpUs+4+bavWa9d3dxGl25mPk8wzk3Br4CncmDUEeMQ8E5aCg4hxAi4hzS5w769gO4VPME7uw0Wr0puZi6dQ49P3FNTlSX3vsSjrn8YawNdQ2dXRx9PaR9FpY8Rq3exEq7DukHADh28ojUudymSXeJlTLUXbFhu7OOEAIf+sXsmryt68FcvJ5QEIcQcXYtz+naiY/nydfW4icPLUqU9QY/B27RJgL+8UoQ52rDtmyKaKPimrVWMsN0rq2jinv0UBg1oN4U0iMHBzb+ZWYn461z6GF6t3Kje6HXcfTlDzvrVEUQK6sWb+t6MBevJxTEIUSsYPWs3wnq8OU75qfKeoOfg1FEIyPLRj/9BvFLdzzLlrOK1hqI47uvfgyPLX7bWU/t7dRd0w6Q9aRyaGnrwPwwvDX3OnznpVMhnbln/vjbC2/hqO89jEcWrk59cOoc27BtByZceC/umLtMb8KIzoiG6oE7rycUxCGE/Fi6Y+elR/0E3GKlevBz4LooRLy4Zu2iiXPISyG9cNXmzNfwu/H6cYKbuzSOR5n1LlkWP1PVjdvaUiHNs2L+8iAV6oIVm1LbCHVtX7Yu4C5ufmKpd9sqcZh5/exM/eqs8+bOhoI4hJDzwn/nVTu2tqYd7VwLSj3MW3VxUT/qKE9xTvfh2mE9pHO6nwrOWq0O6DILX09yiQ7FIKDWZzro2w/g1B8/WtvFIWzvrbOEVeWwn3gtFdjZ+9oCBXGIICdld3AOHOvr0inUA8trYtn1j91EYH1FdtwOzjcYX2fBcQ5Z0fNvikdLuxpKww7bJknu6F2Yv2yDs44+Jzq7CeqMWKmwVkqiIA4hYs7BD1+76/ma71XiRBeOO9fDvFXXZ7U7cTa7jOIXj/u46uaNcg6cQ3fR8db2aqZkS61t/oqbPJ5hs0coGl0/lXC0VGNBeXaoUzoHVmwq8NOHFmHZup3Hn8MXBXEIEXMOXf9lVxji4JrTar96Sv/Asd1EsSjG1S1fb1iOUHabWIn5IurJlHWdoqe54v6FOPCSBxLnbe9ANXro7Bx61oMrcHth+/WBQN7EqlOcA3OTNzdsx5UPvoJ/+83cmtvVUa0KvLpmS27tdRX6PHGYu3Qd9v7GfVi3LdiBdce6y+1OXWIjVdTSUxImVx/zWkT5sBBd99Dq6+iMeWhUvwuJyRdu4y28fO7dmkmsZMdZ1zyeKnttzRY88dpax5XJ8U6Llfg7+4pV87ZWkkWdVcKr+OU/X8P0Kx/FC29uzK3NrkCfJw7XzlqM1vZqJB/tjl0iJ1Zy3VZ1Lu6pfaxpgZYcgWTL12/rXF5pVtHaqRb9wSqkszZSP4xGAm0ZFNK1EOOTrnwUM69/IsM90mWmtd13ze+MUtlmrfTmhu34qeabVCuefmM9AGD5+voWVfV54qCLOrqFc2DFSv4seHcpp2e9vDrx22RtK59G9suUq1hfcEyPkTWgXJ5cBSfyy0uXYoKv3uAWj6CFdWC34ISqZ9D7a5rbvpu2zpijmsSmElc++EqfUloXxEH73R0fF7c7dc05LiIqW68qcM2sxZhw4b347l9fMlcEsPTtrVZ2WZezyoRIOuTjvPyW3a9A77fpg+cW464kiOpixRPuLrs1AODnBmKqYtuOdlzsEbTQ1tXEuNawq88KtQlTrnH9vZru69uf9qrenv+D+Cz8WZMT9WYUxEFPyN4N92R3p45rqonv2lz70UVrcMX9CwEA1//jNWubJ/zwEfzHzfMcd47xz0VrouPX126NjuUYbtzeOXGSBPd0PSlWynr3LAvriys2YSWTfVDvxm/+tTRTHziIxHHXj6i6MH/q13P4OpZrkuV+99QX+E3b/ZJ3ASZrpeTvrGlNbah3Ls9JHIhoNyKaRUQvEdECIvpCWL4LET1IRIvCv8OUay4iosVEtJCITlPKDyOi58NzV1Mo0yGiJiK6PSx/kogm5P+oxidM/OoOSyBu/XEqpD05h6yT1ye8hMQd85ZHSrS7w5hFFP7zge/QcmPRXd6rlXI2hfQrjBd2loX3jKv/ibuYzHl6LwY0lFN1ONSDJz0H6ZC2uaUt4iIEkBpcs87BU6yk1bvorue8+8iJlXRio3MmtaAec3lz8OEc2gF8RQixL4AjAVxARPsBuBDAQ0KIKQAeCn8jPDcTwFQAMwBcS0RyZl8H4DwAU8L/M8LycwGsF0JMBnAVgO/n8GxeSHEOPabtdZz27FfZkxesdRH556I0Mck7HlVPcu5qTmwJW69PveofXdIPXRc2oNEvur5drMQfu+rmiQMueQC/UpJo+XIOazxCwgPpzdFfn3/Lu2+sc6rWn3ZHmJudCc6lRAixUgjxdHi8GcBLAMYBOBPATWG1mwCcFR6fCeA2IUSrEGIJgMUADieisQCahRCzRTADbtaukW3dCWA66V9HN6FbrJW60JTVO1lRjY/JNe/7olI6B5N82dPPwfUIuwxs9OuYgqZK+pOoByc4NaZSHnB1sTtybwvBKaST5yUu/KMfB9DWicXbR9eVp0K6Pnm8GJmS/YTinkMAPAlgtBBiJRAQECIaFVYbB+AJ5bLlYVlbeKyXy2uWhW21E9FGAMMBJLapRHQeAs4Du+++e5aum59J+12vPgS+OocK58XFINfHzDl7Hm/iaL7adKaW3UUDw3plXSjzGFu9779/ajlbL3XvelpxPPqStmBLX/R8Bn+Axatrdy7jFn6d1uQiVooiCnS6qS6Ft0KaiAYB+AOALwohNtmqMmXCUm67JlkgxPVCiGlCiGkjR450ddkLPaGQ5uCaKKo81DY/PWlDzWIl7kX5cw5+9/TlEvLyyFarcTvP7JxDfgtIdvgRUVcf87FWcjeSEitZzvnAFfreBi+dQw4K6XrKfmeD11JCRA0ICMP/CSH+GBavCkVFCP9Ko/jlAHZTLh8PYEVYPp4pT1xDRBUAQwDky0cboCuH8qLmv3j0VUy48N6EV2p0T1YhbW/PN3yGt3LYq5YbRFl0DtrvDMrHWhbcWmLotTEff0/s8LpCaZmYQ4Y633jXvrnf1wZfPwdfdOZ6TteV0jkUpqwxQtn/DQBeEkL8SDl1D4BzwuNzANytlM8MLZAmIlA8zwlFUJuJ6MiwzU9q18i2PgjgYdGFZhdn//xfkbdjmnPI57bXPRrYrm9t9SMO7pDd7g/bp524XufbkNjGPGPWe6pgrZU4bsLRoK/+RQWXVyOzWKkHxQW+9zbVyzv8urUPEFFuB4nOSm2yXn7FBw+MjrvLWkli5cbtWG3wHaoH+HAOxwD4BICTiOjZ8P8ZAC4HcAoRLQJwSvgbQogFAO4A8CKAvwG4QAghV4/zAfwKgZL6VQD3heU3ABhORIsBfBmh5VNXQAiBuUvX48oHX2HP5/Xu5TzzXqwd59U1y9akT/iAxas340/Pps0nAeCvz6/ExIv+iiVvb2XPc9Zdvkl1OjcWzI4+/PvMG3wQuFrY96P2HJ6+TyfmxJbWdky48F7cqFjoeCFD3304AsBvV12K5OE5KKQ9mrhbS+PaWZPlrP1uVAwQuHvr31MuYqXw5V5670s4/LsPdbq9roJTIS2EeAzmqTrdcM1lAC5jyucB2J8pbwFwtqsveWAzk2hH64z19G/PPQIfv+FJa53l67dFDmHcTqMWa6WECZ2lqs/O5uQfmc0v/xAqPhfVkEXNBX+FdD6iHXWc/2v6FFztERtnz5EDs9/IgsdCx8EbH1+Cfz92ovd1WeiaEH6EMOkgzQ9oLdyWCS7Lnq7gsLK22agYIPj41xRipZ0YHRrlz6qQnjBiAFuuKjKP/f6s6HgHoyDjPj/Xmq4u+jYxhzqZa/nOt4TEc1A/ft+gy8Kv+jvPgXFIh8/wqwe44y1NuPBefOvPCxLn1ef/7Dv3xK5D+rH3cw1TZ0xZP/vbINF91h1nlnfn6yD5r1fXxoTXUC/iHPxvb+lXLdd0r85B5Rx4a6VkmemZWto6cPwPZuHxDE6lEmrkgXpCnyMO+rvNqpBuqpRxwt4jcdD4IVHZ62u3YsrX78Ndz6TNDX3trl3scIenn4Nar5Y94NYdIXFoMhCHTmwsvZ3gOELgce2vH1+a+D1QcRwjkJf1Ek+4k/e+ZtZifOgX5vzEXF+z2t+r89IlaklY+Djm0R+f5sWJEtIJkLUYy7jwbtzO5wiP2mPKOq1zyMo5VOycgy5WMnFDr67ZgjfWbcOl99rjmQFITbIHX1zlvqYH0PeIg25XrU1Rt5IzYL3VOSLtsB9YkH7JnJiHW6RckzrJOZiRyPNcw0reEmYLy1O8IOH74f7hqeX4y3MrEtFKefNWe4Mq9+P7OOy70X5fcf9CzFliNqbj+toZ56wWxuJNRZbd8psbghSfpisaGSdAiawL75vr/dKJJu/ha1DB18vMOZRVziF9fouWzc50X/m91CKOHWjYiPU0+h5x0H7rOwHX1CqXCCUc+pGIAAAgAElEQVRKTkKZfrEfE/+GEyvx/XJxDnE7tg/gygdiMU9Pxtnx0Ru8tJJ3l7l93jJ87tZn8JU75kdleURl9SEQtViS+aCtQ+De51bi3wwB6Gz9cM2hpC7BD6ZHinONMPL3jOPQWgNB9OUcTPU6wzlwxhyvr03mXDBxDjKab3tVRPkafDGw0S9uVnejPklWF0J///r8dU1OokA8odaTOzuOOLCcA9Mup+f6zC3zsM+YZnzplL0SMmvTB7CltR2LFA/RzrDoWT6yhjKl/AOqAtBj2GU1CV25cTuIwIZZCBv0hjfn0InbCCFARGz99moVF9z6tGdLyX7kFc5dheldRKTBU+/D4YJbn8akEQPxi0ftUYHZfmXiHDguL3l9vwb7/jertZJpDFSfmrVb7OI0vdc9FCnIiT7IOegvO5tYSXIOaj2ZE4GbiKw4wTMq6/0LVuEnoYWNj0L647+yW1FlQZaFnNtNceP4uzlvZOpDZ0Vbeh9qbU5txpTISK3HPTvnXGeDumDo7e09enDid5ZNgGsIbGPkyznc+9xKXP3wYmc9znrM91F8M8ZxGzYVTQ6FtL+TXjxwQ/o3WO+ZTjDWcxy+DX2OOOizL6tTS6xziK+TH34jE5eHm3Dc9/fkknXWiI+Jdgxd9kn67guXk5QKbgi5st8+kY04EMVqWdZr2nG9el4PLb5FNWlWHskmVlr69lZ8/28vG++XZ0IitRt6q58+LmkSm3SQjI//dMExxvadHvlcWTesYUkuyHxDc6Ko5G+XOW1jOSYe/BzzU0h3BvWaXK7PEQf9PeispOsDKBGlFNJyAvELp59CGgBmLTSbtKn21Z2ZS6aMXDp6MoyBRJkUj90amlOv0Yf8rY0mZSlhgCYDlhuIE374iP1+2t9OQemv/ux6cECh7inUZwbwqaMnsM1f+cBCw23N490d6Wl93QhMXdH7qH/f+vxX5wUftsXePlee1Reiu9L+ZkXfIw6OnYVLnFIuEYh4OS8bM8lPqgQA6G9hgX1NWXUIEaQNXbYuUKz5xsV3M8/5Q0/RSUTRc9eyk01yDvr74Z+EKL34mnaLKzUCI+dEHt/6ltb2iJPU56TeP9vicsl7p7Llf56/IlX2tTP2MYqVOqoiu5d3DfBdKH3HWNcZbE5ZHyl12UxwdjF0fK1CHBwiRPktRm3WKevQ94iDziamdA726yPOwVfRnGGlaO5vtg9QZdZZ9AGrNrXiivsX4pxfz8FTr6/DxXe78xADNrGS331r2Q3p6VNLjl1dFqgiKhf0BcFkhqovAvFl9r76yJiFAL54+7Nscw2apl8Yjm3viuvCecdPUtpJVrhj3jJjyJk88YwiGrWNkukb0L9LfXOmmkfr7XDWSinOwcAUqMTBJnr64m3PpMS/dUob+iBxSL1sfWdgv75EQSrJxGIdNnrrnDdSjnC8WMnUNn9iw7YdNXMOcsJv39GBj1z/JB6xiK6S6NyMrWUt14lDUinL3MPVR9XnA2Zxnko2uBqmnaBu5uhtKeRZ7y/PrQSQnpMNFTPnkBClWcihaa6ZrJV0e/+uwsV/ijcvtnEymrJqv/UFXxcrqae5DZ+rPa4dTo/ZURXY3NKGPz2b5tgKsVKdwPWyr51lt7IgIjRVyom48bKJDdva8KXb5yfqs8TB8NGa5sjB334Q23bEH2eWySQJF8HOmfj2hQA8v9ydfKWWCa+LlcqJBcy9q0udV46zOMHpRMQkQ9ZjbElilaVfKv7zhElsuU4EdcOHWsbaNB6mcrUPnz9pcub7HTt5ROZruJH671P3Cs54OsGlxMYWHYKXtZKBKqnrCJdX/PL7XsIBlzzgvLae0PeIg/Ii3nHZ31MRPZ2B+RCYrLYqOxDbq2V1Dh4foI7tbWli5IK+2JpCYnBOVibTy6Vrt+E9P3vMeM+ssXnU95HWOcTHtbDeSYV0kiQbF0fmnK9Fm6zG1VYdnUyL+aeOmcCW69Vtec/VOWQjiC4zYdscmzB8IL562t7W61WMG9ofo5qbvOtLcMMuTVN9TVkBu0xfPeVjrWRqSiUsV9y/EM9pochtuawLnUOdQH3/LuXsgUr8JBX9GnjOgUOWXYGtaqunlZGKMpnFNCpO+3EcpVUu1qbd8rqtdgefSrirFZ4GG3cq6S/LWhq7klOsZMamlrZ0ekmLBVBUhRki36B5Nl2COg9M1UwLtl5dV6rW4gRnToQUWitZRrdSJowY5J+fmwg1SSk5ohiNkak9ZgDaGZHsd87aH3/+3LEJ3yQfa6Uv3fGs4bbJiu/92eOJ30MH2H0f6hF9jjhkwVhDFM+mSgk7OqpeNs/cglHRXYc9sF0hDrMWrrbUjFEqqZZUZFy8uNwNpgXRJZ6RegNfpbma81dfsB5TIlxmdRRasSF7XB8TfOMiRZwD09WWNnf4k7TnrGwvqL/X6EEAgGEDkgtzLRvPkoE6GLla5R4N5VKmTHVEtWmw9HdeJor6ZzYp5criwmtCsfGuQ/rhgPFDMH7YAPxk5sEADNZK2m9TKBPXWmAjDoWHdJ0gmxkoXy5ZW5kC1BpCmznlu0NUsXJjnDHqjnnLUueXMgt8RduJZ/lAuYxoPpDEoZYFy+6dmy7LSjBMYqXEMSiS6cvQCt6OT5FYqTYxgb5IyHkiH/PTx+2Juy84BodP3CVRTxVLuPQssqwzYqVKiTLZNBPMGxNrH7TfpRLFY2K8Jn1GEod5S9fhbwsC8Y76+DP2H5Ool2jPs98uCcHQAWZOq/CQrhNk+XBNa4IU18hFw7Z2+Ci54nK/vsmwyio4B639xzV7tceh1oxXkVipFiWpZcXxfW+2oGdeIbspjpL5i48fBgBo89Y5+CmkTbtefTMfZ2WTvwkH7TYUAHDJe/aL6snzbR1V/OkZe0huCV2/I+Gz5gecgz9q5xySv2vlHGRVVYyszgX9e7b1wQTXBqK5XyFWqnu4Xnbyo7GLVnzmDZ8Pmb9SLbVZBLmCiUmo+QyAbFyTKUicqw1pg++72f7FP+LgbGY5ONjB5m7xuyfNITp8F7Txw/oDAEYODpSoHb45OTzbN1uCGTgH6YGvnJuixFeS53/+yKu4SyEOJmK7bUe7UXfkQ0ArZb/cGHE/ajNt5jYEtnAq0O5zRMhhRUSbaQeIv3mWOHi+1R/ez3ucS+hm2ip6rViJiG4kotVE9IJSdgkRvanllJbnLiKixUS0kIhOU8oPI6Lnw3NXUzgiRNRERLeH5U8S0YR8HzEJ26seNbgJwwfG7J9ZaRl+tFVHRfgpubjyf1iyQw1o9DNJrQoRtfnmhu14Q/PMZPvgOK+HltAhRVl5s8oZJTudwk8/cgiu/dihEZHwbdPEOUzfZxRbTwdpX6NctGK9kVKXseRa6zAWkPj338x11tG7qP5szMw58JFq3Z1I/uyoishoQbeXWLZuG7a0tkMIgaZKCc9dcipOnSrFRXyf1GM96kHUBY+Ob25pw3yHeXedrv9W+GxBfwNgBlN+lRDi4PD/XwGAiPYDMBPA1PCaa4lIribXATgPwJTwv2zzXADrhRCTAVwF4Ps1PosXbItWpRRP4useeRUPvcwrfmNzzaD2ghV8XgIAeFYzlW3rqOJFQx4D32XohL1HetXLc3k+PZTLNlmSwQCxsr2We9sIALuDs4gQ1F3z4xeeFJQldAs8iAL58BkHjI0WEG/CZNA56ClX1bMqsdX7VNbk66qeQD2Wi5rNFFjFE6+ZExVFTnCWNziwqZJpsQs4h87rHNqr1Wh+6UYCx/1gFs7++WwIBGPT3K9BcehLE229+3owTVMfAOBfryZTgf7MEIFW7aOVKe6tOgchxD8AmGdTEmcCuE0I0SqEWAJgMYDDiWgsgGYhxGwRjMTNAM5SrrkpPL4TwHSirqOztg+9VKJoAsnAZNxOOWZtg78mIgIAN81+PfF74VvmTFEJe3Vmwkg/BV9LkarIrhw1texSBErECunsE952jW9z3POOGxpwAD7jptaJZf5+N5f1Nm1PmprqJsWqma/ui3HXfx4d318bSxNxk/fVreBq+Yps/h8S/RvL2dq26BxserG0w1ostuR8T15auQnVqlCU7vF1OnSFfJnIEFspXfbRXyadH00cm29u6Dp1c+iUzuFzRPRcKHYaFpaNA6Ca0iwPy8aFx3p54hohRDuAjQCGd6JfDpjfRLDzSn6Ix04egfn/e2qiXili90Vmqi93E5e9b//ULtzVUiQbzbBYvfCmmavJApciUCIWK2W/h+0j8XFQCgszw2TFVKJshE7W+sgvn0iU6zt6td/qcYmAQ3YfhlP3Gx39BvixVE1R5bjZ5NpZYRMrDWysZDNl1RtQcNI+oxP52JP3TF8kgw6azIurIn5v6neqt6cTt1IpwxzTsGl7G1uudtG2361TxqFm4nAdgEkADgawEsCVYTk3AsJSbrsmBSI6j4jmEdG8NWt8YwRpDVteRJkoJZYY0r8hFXZC5RyyJnGRH/L4YQNS9uq+ynL/naxZsZwVsdLOXq9SNi+oUjRlhoVzcFzpg6TM3mCtwxIHv/ZNRETf0SfCvSdEHdpuVkvZSQmxUvq+uulylgU8usbjkv4N2TiHQOdg0LPAYpbKnJDPaCIOAjHnYBMLciI8TiH90EtpqYBO7FWnRPVch2fo7p0qtpIQYpUQokMIUQXwSwCHh6eWA9hNqToewIqwfDxTnriGiCoAhsAgxhJCXC+EmCaEmDZypJ/cPdWG5VyplFacNVZKqYUkUkhDRFngvO8fToQS2Xcl3HyxWVVwqGXSma6QC6Vrwlc0JaoKV7+tpy36ha6Czi2ZPOZd/dEXk2RynvT9ol2T5uegzcK4vfCV6JyDXEBr8c61DW1WJ06btVJJ2ZD59CESKxk2ZXOXrotGRhcL6iK8RD9KaeKwqaUt4aQpoYfW39wacw5jmmPH2STnwHY30b96Q03EIdQhSLwPgLRkugfAzNACaSICxfMcIcRKAJuJ6MhQn/BJAHcr15wTHn8QwMOiC0fLzTkItHdUIycwNXb+qNC0MfZcBba1ZYtWKecf91E4U5TmtJOtBfKZXXGGYj+H7P2xnfd9ls4SxKTOIbk41+q/kHJGVKtx4xSOMWlVSg7Ooawt2lvDOGGP/PcJ1n6riJP9xB174c2NuFuJJhrkNMkgViKb5Z95k8R9D/J7/I+b57HXvPDmpkicJJ+Ft1ZK/i6X0gppEwHSpXetive7+qq9xb9etbofTptIIvodgBMAjCCi5QC+CeAEIjoYwXMtBfAZABBCLCCiOwC8CKAdwAVCCLm1Ph+B5VN/APeF/wHgBgC3ENFiBBzDzDwezATbbr0ccg4X370gKpNesr899whMCcMXqItGW3tWsVKs09CvdKwZmcVKeSq6IhGLizhYwmfYCMv2HR3YsI2X3QL84sKWhX9rNWlIipWCv/KZXfJn2+44WS+pc/j4kbvjvQeNU4LKJRviFNLsrlojQlvDSL6qd65rUed8eN7902SgxXIpm8CKYBYrlSyEgyuWXMtqxaEtlStcaTtoJ92Sj7WSSgT+deFJ+I+b52HBik2p0CMqZ6gaH6hctu2TrVexkpM4CCE+whTfYKl/GYDLmPJ5APZnylsAnO3qR16wvQcZi+ivz6+MyiQbe+yUOORwrHMQmV+srM+FL1Cb4kQw8sPIGs4hD8hJ7+Ic5M6Oq2br921z7fmluSvZMkv31Lg4SWsfsOW6zsH1qk3ndTFMYhMgAu9ZNSSGTuA4sVIrE6tJF19N3dUuBuMQi7IsmyjFU9mvTRvnYGlLu2ZMc79UFjwgPdfkM9h0RvoC31Cm1EZP7fOuQ/vHim5DmBP92DfKwM5ordQrYfvAZQwY9d1zk1E163QRh2l7DEv8lpuJElEq0mrCioVptyfFSnJT6rRWMiikhRD456K3uUvC8/b7Z30Wrvr2jJFtdZ1DXn3UdQ764tgRiZU0nYNSr6U9fpZYVJlsZ3RzOnAkl2tAhc1CKqpToszKblNzRLZYY/FV5xy1B+75/DEJvYrk6PRNh7QCjN4fk2pWv2Nzv4ZUljgTB6ePc4JzSFiRCfZYR50yDn2QOFi201LUo757jjhAmXSuhVo/He/ygE16di2l8tWMY42ceHnI30/SvHa5PqiQH7BpN7TH8AEAzAppPWuajvteWGk9/4gWiXZHe5VVjovob7qfJuJgMnGUu9rYFNIO03ndt0UdGyFEaqHViUHcv7ieurEwcQ4ctrbaCaRNTp+ol9VayTKvTKIuda4duscwjBrcL/GMG0ITUn2eS/GcrtBP9in5e0j/hpRY05RCmFNmSyStleI6duKQPrdtRzt+P29Zjyqr+x5xyPiF61m3AD3PgL1BXQwTy48ZsZKjaxFx8GQdbLUmDB/o1YaEr82/KXzGqUrOCA5zl5oD5gHA7+YkI9Hu9Y378Nnfps10bbv8Fq8ETWmRQSxWsj+7MSwGmetxnIMUZeoKcXXtV82ghU5NLHDRD1/HP1/acNyUEVGGCNP9TMZPXGRg9dvrMHAOkjjoOgfbMw3qV4l0NBIpgxGlzyrU/lOif6roL3nNEYoYkZs33/nLi/jqnc9h9mtrjX3uavQ54mADQX6s8QvWk7nLeoDUOdjbbNcmePyhx+2eNnV04pwJ2e3uzef2HjPIr5EQg8MQECqxO2T3odGx7HvZED7DFAc/d1j0Ay2JbHrq9j0+1NdXgpsgfmjaeOM9AU4hnTzWZ9jlHzgQACI9BLdjPWLP4fivMF2nSazEwUQ/IvPMSJRWWzsqBjVV8OtPvcOqc2hpqxrFSqpeRYIT2+jcrAxMmeV7IaT7yIlG1XYlBikRV9W9JJdkSMWuYb4Yrn9vhSH6s5rK54k+RxzcSkWR+Fg/OG23VB0pfxfCnU1On7ic3FIG0nNZw9TiIZ0XxoVB6FTnI25Slz05jK6CHMNa789bsYRtG5qsxTkxedPkXUc398O4of1TYiW9b0fuGQQSkJykK0dD0Ea6zvWfOAxzv3Fyog23Z7D7Xv0ayqiUS9aQ3XfMW2bsN8c5qMRBLr4/m7UoUad/JFYKfvvMBc5Rz2RAofe3QekTx9kA6bkhAPz1C8eF52z9sna7S9H3iINV5xAsBPKFDB3QwOZdjuWyIpVkXofOqqp+DhK+u+pAeZfFlLWGBdKxs3x8cczm6iaZgN0Jrp5g6p6ul1Dt8E1zJyIOhjb1Hb0u5uCGvFSKGzRKjbRdvh7DyRfDBzUpcbuS9zTBh0uRdWzJfjqqAiXDKpTI0858N5Io6oEv0zqH9IZB7w7HOZh+6/0VACaOGIinLz6FDYioH0sMHdCI3Xbpz45NPXw+fY84OKi0UFTSph0NZwtuwvL127FRib2i2qzf9O+BY7kUz4hoMeBbPmDcEGP0SA4uyyz+Ir6YU3ZyH5stfAYHUypWHe8/ZJy7ktIP1/0NUqVUxi7VWdHIOTi5peTY6ZwIN80I8XuWrer1UvkefBZs5otX323Jk7j7Jk6Sf03NBWG4O8c56Hq9/o2azkHOCcsejBN96ZyDyQy9KgQG96tgl4GNbP/UPqTuC0M0WBGf7yn4JQbYieBasoRwf2Q+tuAqvv3nF3Hlhw5KXFMiwjv3GomXvzMDC1ZsTPRNnVRPXDQd/RvLeG3NFuy3azN+P2+5M76RRBb76cCE1/zgrH05048yE3ivtZ2Xm542dTSeXbaBPadDD3ttgmshV2qypTKCq0SJeFPIRB3HgqoPK7eT5a4RWn1XGlEfcIuNynGo+jSgczJvyRFzu3IJNVCeDlXnIAlg0slMJP5KDIgU0kmdg10cmw6dk9I5hH85HZJ8N2r/VMOR3z+1PHHNIWFGv5KBcEZlhVip++Ba0GO+wQwfW3AAeNeBQZQRTk4vJ1jAAieJjTrZxwzphyH9G3DI7sPQVCmHjnq+nIP/qtHuWABdnIMEJ1baZjCfHNTU0GU+G1nq26rq1kqTRqatvKRFW0dV4EcPvpI6r4+cbE+GpLhFC+sur5H9mhWGhE/rQ2R7UkwVn/vD+Uexz/PgS6tSZUlT3qSI7OK7X0jV5/rCYb3F412iymxKhBCYs2Qdb62khqcwEAedc5DjY7PyU02W1b7p/ZJ1U8/A9M/mMPrV0/YO+0jsN2ATOXYX+h5xsJwjrYJpceVswU/ff0xk0SNx/jsnAQAOUxzhOIW0LqayTaoS8dEjVfzHcRMxY+qYTDvKlrYOCCGw5O2tUZnqtctZbelWNwDvi2HyL6iUyNss198pXFjvGdXzbI8ICfHOPmPT+QckV/XPRW/j6ocWpc6nN8ZBe29vCYwZuHwAJYp3stLnxRQAMiJeyuTdbZcB7PNw+i11vulK3DcM/imZ1BtkzgQ3jAkKeNczb+JDv5iN2+fG5sv6/AIsxEHTOchnsXEO3OPoc87EOQT3Sp+zfacyBhnIxUH2HHnoe8TBsihIiwXX++Bitpx1yDj86pPTEvU4KxY5X4hh5WVzHRa3+7Jhp6GirUMY49MD4aKvPWNLWzWl2FMX7jIjrOYUbhVGObvNIJqolLPoT/zqjRoc6DC+cNuz9vYSxxZiXFKUqYJfRCRxMClp0zLq4G+Fc7CUYBaNLJyDHuzPBtWirqTtVIyhtjOsWQEXxLfz048cGh3LubM0JEhq/CSpMOfESronuOQcdOW6m3NIlumLe0SgGLGSHDeVeHGcj45Ap5XuV6xz6Dn0OeLA8Q6TRwU2/+OH9Q/tzu3WJ7FbvtKqSMdribOixWV3hwngk5xDUrHYbtGcEbNo6G7/m7a3RYqufZWd7seP3D06PmBcMu5Oa3uHUQGnPouKz0+fEh3La2WgQvXabTv4yLWVkpvQRX1xfGeSyxk7NFgoXCbGvkiIlcDv5BoqqoiQb0OFHBqOG5PQuVjAopDWFNdANssldb7pBEfHBSdOCvtn3j37lgPAqOam6HjyqEEpp9PmfhX84AMH4hSZAEnlHMI+ThqVFPX113QO8lFsO3kuOKA+BvKbGTG4MVUvEislHGSNt4tgCzzY0+hzxIF7Eee/cxKevviUgDggg0Jam0z6+llidtEypWiJ4Rx8JnGJCS189d+TooxB/SqRQlPdLQ1WnHX21cQjnDPSobvH4jCOOBw+Ie3lKYmD2kUuIdLPP34o+yzDBzam6qrtA8CTjNeoXGh9PzS13qdv4sM/A6FCOtqZC3YnJxe0RkN+7ZGDmxK/fcJdcHb3poBvcoOqPpMevtuECcMH4MDxsTOjLjLVx/ML0/eSFVMwZaKzKaQrSoRXYizxBvdrwIfesRur8I08kLW2ZWpfPR6YOg317nCcg74h+faZUwEAuwxMvk/ViEUdAp8Amebc1SLqV0+h7xEHpqxSJuwysDEOvCfrGt4tfz6tWNN3dtw5ID0BbDqHMjOZ9PoXnr5PZIKpylltu8mWto5EPw7abSj+3+n7xOw88+ET8yHExCH4vfTtrXh97dbUtTP2H8vmtPjRhw9m+ycfcd3WHfjw9U+kzvuGFJdQF95Wi5+J9H0JrgmeeZ8xgxN12PhbCvRFM+IcLKIfbkepv4E4+VPQfxenx+GRr56Y8OWJpEoGOb28J9e6ObteQOhWbWpJnSuVKLGw6m9PT6Wb4BzC1+aKreStkNbK9HaHDmjEniMHJsq3tLZj9mtro5hVpsB7NvAK6bBfPShY6nvEgXkRUj4rJ4jLVJUzIVTljhI2IsOtCzrnoMZfidqkOBH6snXb8PaWVvzmX0sTdQY0VqIdr7p70cVec742HT8OF+PW9mpisS9TsOjtGVrneBMHmewHwPqtO3DCDx/Bl++Yn35YxLvyx5Rorf0Mu2/5LloMimbZvysffMWYQjLZnrNK1Ec1kxj3qdp8O6btMQxfOHlKokzWG9AULGJcpjbW/l27ueSWJGem1vbxluYQK6SBDdt24Jk3kqbGkWMb175lMyUE8PNHX02dUzcs3C5a58aSfgSSKCbb1HUOPqasnKMe9z51cd8l9wS5X15cuSl6huh6T86B61ZEHHqQc+h7fg7Mm1B3Q4kF39CGbtEh68p5O2xAA777vgMwIhQncEq9pFgpKX6SnMDMw9OhO8qKKetxP5hl6GHMoqvPo6/vo5r7YUzohNba1pFIf6inJOUUnOquRn4Hkc6hKnDhH58z9i/oT9DHK+5/OdFvDq5dmDqer65Jp3asFeqCJULuUN/NRrolZjH41plTE+I8FQPDsCnXfPTQ1DlOzKHvIiuKCW3Qweycg454HAW+8ae0GWucKyF9rTmhTzCG6xirLDXXBffM+lgn/QjC+2rXNEYGAvGzAA4xjwfnEPQxFve1tnfgTs1/IekhHZcPG9DAmvY2lIlVXJtCpnQn+h7nwJRV4u2Q+sdIHdR8DlG7Il5Am/s34PQDxlr9ITixUuznEEwWzkLI10NafmjqB8GJleSzd4hkGAO5CMjrSyVGp6L8lru1BoVz2KyHJE/1MRDZVC0ETMIZGl0ZEx9W3J9zSMZWIqSnhXxPbIImi5WX/KubQEd9dPRNvru2jiqEEInQJjqX6AuVc7C9P26MTe8oMKLgzzVVykq99Aejnlf7B8QLqP49NGjEwSeqLqcX4RhQtR73vhPcNGOlqKO5fwO2tKSJhmIg12Poe8SBGW2pvIsVY45GWM5BYMKIwLZ879GDw2rmycktCDrnwO3+SopYyYZI52ARKwFJW3mVYPXXUlZyeYNLRJEPhww+1hhlghPOBViKbNRxHNK/Ab/85DTM0vIeuwjiKfuNiY7Xb0vvUHW4A8sFSPg5CLAmpjGB5a9P3Tusx8XZkuDEDfpiJHfd7VWBW+e8gb8teAsAcPVHDrE9khWqyDQvx8NyyRxbSYUuBgLSYqWmSikSt5oW3EoU8jxsrypQrQqr7wtnUqpHVAaS3A1HHP5LseBLmrDznR3cr5LO66KgJ1OIOokDEd1IRKuJ6AWlbBciepCIFoV/hynnLiKixUS0kIhOU8oPI6Lnw3NXU7jSEFETEd0elj9JRBPyfcQkuLkpdXoAACAASURBVEWhrHMMUV0enC5BiEDWf+9/HRuFyohNxoOKkn388il7YSCrBAz+yknHyfldHtL3fO6YqI9VIRLKarY9hdCp8vwrzg7CRv/4w4fg3QeOxeSRg1J7RSLglnMPx9yvn4xrPnYoTt9/TGRKCuFegEsM5wAAp+w3GhOGJ524XN/IRxQR3I2PLYmOLz0rlZk2E/RFmpAmzlIBmjUgohpnSwdRur39dk1amEmupL2jijfWxc5qx00egVoRL9D2Z7ERPR3ScdN0fr/Qcm5UKIZV760TByLC/zt9n0S9NLGOHcyC88B/3zkf18yKdR4pkR0jVmrliINi8soRB9UKMBmVNVUVADC4qQFbGOIgq2cJgZM3fDiH3wCYoZVdCOAhIcQUAA+Fv0FE+wGYCWBqeM21RCT5wusAnAdgSvhftnkugPVCiMkArgLw/VofxgvMYEv2z9cyQO6uOG/TqbsOiWTM+qLfsiOoL03tJOL7JiediXMwfbgjBjVFZokx55C2YU/cW1G+f/aWp6Jy6Uy2367N+NlHD0WlXEor3IkwoLGCkYObcMjuw3Ddxw+z5pBOP4vkMHjZrgr3YhXXN1mMqCalWRTSCVNWZtEe0j9436yYgWlT1jLlBwCSi1Vzvwo+dfQE7KKZ+cr50V4VyfhIzE3HMClDOajhM6yeu9pvG2Esaf4su2ve218+ZS/84fyjcVAYb0htSdc5AGlHvbRYiRL1BAT++PSbxv4BvFiJ+75dnIMKn0xwjZUS7yxneLbuhJM4CCH+AWCdVnwmgJvC45sAnKWU3yaEaBVCLAGwGMDhRDQWQLMQYrYIZtHN2jWyrTsBTCd9ZcgR3FBLRaK+mBvDZ4T1VOczrmacAzg4K6Oz6s5S+n3bbZwDESu+COrHx0ve3oq3NrXw3q+J9uJ7r9iYNjXUHojtN1fFR6wk00f6GAFk+UhU8dkohSDc/8XjvduQSJmyAvjme6ZG/hj9GkoJOb0OSZDef+g4vDuMtRUH3gv7yxEHxXpGCL6OFJ98688vJu7NbXJ+/W/vsD5ndK2i/7I6Hmq3sK2TksCqj/CH84/CP756IoBAsX7YHsNYQw/Od8TmGS7bC+oRe54D51fCE4d4PriIg49YqaFcQhtzH7lhnbtEX3q7D7XqHEYLIVYCQPhXJiQeB0DN57g8LBsXHuvliWuEEO0ANgIYXmO/nODeUZzQXatraENOOjViJbfD0Rf946+YBQAJq6BEPa0/nDJT3cma+gUAc18PJtXm1phlNREbwH+nb7qfhLrzdDUZOXEpz6PvKiVk/9Tw5xL//J8TE7/VXfRpU2NdhLrz9qU1KqcgRPD7qEnD8dBX3hk9gymF6sJLZ0QhwH/0oYPx4Xfslrg3F2dLYntbB14PRUVVIdg6qn/FZlWpydTVQ0yYoC6oNrGgTFAl8dU7A3PlQ3Yfir9/+Z2Jc9I3J1ayChy2xy7YXRMdcos550Oi6/L0uSs5Kpe3d7LN5H23trZHCvkvnbxXop6c2TZ/JMCskL7tvCOj44YKb60k2772kbT5b3chb4W0iYu2cde2c8nGic4jonlENG/NmjU1dZCb8CbZr2lOyXrqzuKUfUen6nFWTQDDOWhDIHf7psXcZD+tLtbffu/U1Hmbr4LfB6SJlbg6antMk4P7VfCJI/cI+xuUqYo/dWxO3z9e2OUCff5vY9GXhB5kjgtNUiukeO6ZN9bjrU0t0RjIvyWKnbj0naRuaaMTYi7OlsTi1Vvw2pqtUT3OmEAVO6p3tnF0Lqicn23tO2h8EEpCiquk2OZjR+wRhaOJ+xPoaSI5uoEjie6t3PjhMKJAsr3grxACbR3VhOgUiAmKjaObOEKLrkvJ733qN+/H1+56HgBwdpgGVrbpK1ZST6tchMzgBwQGHLpfzjWzFmO+Zyj7rkStxGFVKCpC+Fe+weUAVOP88QBWhOXjmfLENURUATAEaTEWAEAIcb0QYpoQYtrIkSNr6ji3BuofqWknKCF39KrCirUECv+u2LAdF9z6dFQ+YhAfIkKfdBUmBEKZCTkhcczkeNI1MXF+JHFQFxWbF7cOfdExWdkEDfKE+NxjJ+I7oZK4pMjMAeD4vczvVL6jpYYooYk+eJhx+lorSU7tfdf+C0B6DIhiYqRyNS6RG6DqHNJ191MUm7pIRkIl9qpowhY11IXE7l0bIpU7JiK868CxkSOfxMQRac5P+uYI7bmN91bKON8IddGf8vX7cMe8pK+BXGx1D2l57dLL35UKaWLTNzYmnjvN4ZugEjlT1YZyCVWRbOuK+xda2+0u1Eoc7gFwTnh8DoC7lfKZoQXSRASK5zmh6GkzER0Z6hM+qV0j2/oggIdFnsmPNXANTxie3EXoDmA65PlWR1hoOYlvm7sM9z63MirXWXzdqqk98nPgxTamiXbxu/eL+8isBkdM3AWfPnZiZE0F8GIlPSifXlfvd6Is/GvSOSR9GpLjLBW7HEyKaw4+QedceiW1j8nFRS5gIjovn/oGxUqKa5bUhReqWIl5V3vuEpk7m3QORITvvf8AAMmcISxH58k7qJyfTkDv/OzRid9lhovVuSVAsVYKf5s8lX25WN25TYe8XFdcq+e4e5vmQoIoKrogVy53LmqxDsnl+Hj1dzd8TFl/B2A2gL2JaDkRnQvgcgCnENEiAKeEvyGEWADgDgAvAvgbgAuEEHIFPR/ArxAoqV8FcF9YfgOA4US0GMCXEVo+dRXUCTC6uQlzvj49MhGUc6ms7Wh1RMTBkfvZJNJImedFfQOeen1dtMjw1kpmt3xVPsuF56iUSvjGu/fD2CFxtjPug/z0cRPZ9vXecI+XlFmnofY9EitJRzvL+lUVwKbtdqc6CVN0VA4+SvOEslcTIZXI3u9kW/Ke4eISTh+Xn4NJ5yDvDySDwbHTLmsfkZ7fU0YnxUXlEqUWSFPYjkWrt0TcjWkd1ImnCS492f7jmsN6COt56hwQcCo/1XJyqPoVH87hxk8FofvveyHwOxEW4wxpWeUT3ru74QyfIYT4iOHUdEP9ywBcxpTPA5AyOhdCtAA429WPvKC+ozJRZLIJxLsr186z4kkcgKSMMrqv9qWrH+QHrpttrCfL5GQf3dyEVZuC0NQDGsuJRZH7SG1+DgmLIdO3pItULGGbTTt99UPVg+XZRt20q/vaGfukymzRTiVcme/iPibvre9Z+zWUvfUautjE5udQLsVJnQLiwN9D3luGDQ/66Pfube0JIfDc8o1R+ZkH75oiuoH+K3k9J66Ti+QDL66K2uYQz0X7W7FxGDOmjmHEw9bmojaFAC78w3NRPyXUsSOFaLcbzAZP2iepf7Q9jtwochZLPY0+5yGtzt10Zq3kXxNizsGdX5fdxKUKebk/G8+IYlNWddJ/5dS9E/V85c6cXNbHGir4zbUX/GVE1gACnYNeV+4+bcHifPsE+LHoF4dxg3xEGBznMHxgI7548hT89tNHZOYcdJ2D2WRZEgczFyqfv609KVfX4RtryRTyhbu6UkpnJRw5qImpmYTrXboWc1tYGpvBiQ3SuY2zhkvWMzvfmWCrF4uV7PqYnkCfIw7q5NHXXjnpXLusiHNocy9C7A7eQ3Zv6kdZ2cmq80j/+P05hzQr77LSin/b2kuy0vO/eSpe/PZpGK4sHpHOoUPuoC3EwTDUtRIHGUXTTRySdX43J7DUJiJ88eS9MGnkoAzOk8HfWFwky7l3FXBUNqU1EC/aCZ2DRdzn7qP/brtUoogDO27KCIwb2h+jPExmTeIYf8u59IZGQi3K5ueA0GnUj2sBgAUrNporKrh93jLjOV3noN9/hmKK3d3og1FZ4+PUIi3FSg7i4KtzAMwig+R9030D3B7S6i5Dt9DhUgxw4rIsctksnENV0zlwymZd52BbvzgFqakPLR5EWyKrzsEVbM3RGoD0zpN7BvmebQQEiDc4SYV0J8RK4V+f+VAuJZ9FRvh1wbz5CDcLjgXazjko9UqyngfnEOoS9Oe+/hOHJesp9120yh39d2trO75+Vzq6rYSuc9D1nC6ld1ei73EO6s7CIPt37bIiU1YvsRIn3tHvK2XRyYlgFDdUBRav3pJITJ/mRphrPZ3gPFUODs4BzpVX9qfDQyEtBL+b5cbIlJbU1K61j+ReXHzFONGiFv62+TlIcZZrbOR4b9hmN6P1FX3p7yRuk+eAZb32DuFNgN590Fj+3uHlnGdysl5QkTMYUbsg6y1fv92jVzI5VrJUNd4AkPCkHjs0ec4HD3wp6aXfqHEO+jNxwf+6C32POCjHJiLgks/KHUnenIOrnrx3VQCf1ZzBdE6B5xL8ZPqusCE2ZNl5ygVH7ppsprLBLprbtac7taXVnzj46RzcdXT89QvHpcr0JFI2kZF897PDlKi2LGsA8PaWOGc2V9NXaS7NNvUFWkbfVVEqxaasHVVhzYktcezkEfjOmXwwRLmRavMU7eiitFGDmyIfmqC9AG9ucBOHoE2RcqhrakhbFkZiQUs/P3/SZBCl59deo/ksglJnpBMDlxd2V6LvEYeEtUzyXGTl4NQ58B8QB5Y4GD5Ul1UTIE0c0wtlSuTDvFlWfFFKLliAhXPwWGCSgdvsddPjb65rIg4cId+2w8zRffCw8Ynfrk+Ps8hx4Y7PHJVaBIK2wntqiwu7Kw8rn3PjnKgffP+Cv6oitdYscEBsBqxzxR87Yvd0HxWleXtVsPlHdAwb2BjFPtIhn8WURzyuFxIRZSGdNHIQ5nz95KT1YYZxkIu+vlFv1PqqWh/aRD7SFFmdO4fuPjRVr0ESY4Zz6NeQ9p7uTvQ94qAcHz8l6ZErp5LLlFV+uKaUlSpYG3Zt1ElbNCS4eDim8Mdp81hGnMXsKVn5rUkmzBdr95XtCazebA/kp4+NbukyZVSwwA5oLAc6DKZf+s4OsHMOx2rhrH1CU2eNjGmaPnpMIFuARb0Nl1ipraOKfcYMxtLL31Vzoh8gILYlSnPFJgIWiZWqVS+LKFsVLjKBTa8lw4uY6mUZhkjnoO1ouJDhsn+q6O284/dM1JPvtE2hDpeedUDqvnGq15A4KHKt4QObEvHbuht9ViF9/ScOw0n7jEqck5PO21rJh3NgylLik0Q4rwAjDCaBMvyxi3Pw8RJWr0vqHPjFcPXmVraca08I4O0t9qQ7+jB/7qRkruXPnzQZR08ajusefRXrt+5wWqdIbLUQB/3ZhOMVSmJsg4/ZJ5B2gpPEgRPH6O/PNCdlcVuHwIDG2olC3EdCU6Xst/FRfG58dQ7WGsqzqP3h+ggk4y7ZQrn4JRoKOPJ2jU0cNqBRq6fojML3t+iy01MBAiNHWuVZBjalnTN1ayU1UdXgfhVs7UHi0Oc4B/lqxw8bkGJv5fRy7bxk5jg/nYP7w9cXDbVMhzStTBEHD1NWjuvnzAd9NspfOWUvtlzVOchFdZxBcaeOTaVEqV1apVzCEXsOj5Sz3BrNLdzbWs0flP5szoREJTfnoJ836wdk/eCvlC9z/iz6QstFJw1bDdqqVtlYXLWgX0PJa25XFFNWm87h8Am7RMc2kVccTsWeg0SW7aV5bOvIIl2LrZXisg9P2y1Kf5toUxMr2fR7qliIe4cVhYis3dKK03/yz+hcc/+GgnPoTsjvmJs4usLQhNjPwcNaySLnT/VNvc7Qnh7+WC1P3iM+vuj0fTB0QAOb6J7bXelRNXXsP64Zn58+hT0n23t1TWzmN7rZwAUpfbZ9yBFBZAgB96q2WKyV1CYWr96cCDvBYf6yjZkU3ICNsMsTQSd2dJg5B53AuDiH9g6zF3VWNFZK3j48QgRzp0PTObxjwjAcOzkQ2/7owwfh2O/PAmDXA8hnUcXsrLVfWLbDmaskJDZenEPaz+FjR6b1LA3lEja1BfNBzkfue5Z0QCUOHPFWRWlrtSCDzf0asGjHZmffuwp9jnOQr95mRWRyi5ewZYLT4fO5ykmjLn4mEQ6FCtK0zoHvIxBwSR9+R3qiq/WqIgi/fMzk4Zim7PQ42BSP8rZqSkbTrlf9pmyLhvQ14Hbw3IdvWwvUNj7yyyed5sg+hCGrWOkHYdTN9o5ATs9zl8nfJnm+ukP1NSV1QVU0A+YscmqAyjZN5/D7zx6NL5w8Jey7EvPLQ+egcg42s9z2RFTkdD2K6vlZzgmF2wX4GF3jh/XHsjDPRocwi9I4c1sbh1gVaQLXv7GcELF1N/oecZCcA/MJN0QJ2+2LfrbYSu4PVu76//eeBc66RrGSzjkov22LhioGaq9WsYcWoVaFtFixSS+4x+Wyeen9so2S9DX46cOLU+f2UHI5/GTmwZZWQijDtmZzK4763sPuaxT84IMHMk36iZXkO5GK1PaqMIqC9HdmtPCJdqj+fgYSn3nnnmx5sAGJn0kPuKf3sUMIdFjur5bbuJtoMa/aOQI5vuqizxtbSA7Dz+JHIMk59GMizDb3a8DW1na0dVRxzaxXnZGbVc5BTw8MJMW6+vg1VUpOfVdXos+JlWxeqb4mqrqH9L8dM8FY1+d7le35cCIy8J5OHPSJpa4ltkVDfkCt7VVsamlPme6pkIu8T3uJ6wxtypDULkhfgzufSsbtv+MzR+HwiTGXM9whIgLcOgYX9h3TnCrz5Rz0D72to4oGAxemiypMnEO0UFar3kYIEnpGQgk9Z4iL2D31+nqs2Nhi1Dmofbd1kfOQtlkrqYu+rd4aD0MKqUtQCdPYoYy1YDg2f3luReqcipg4BO1d8cEDWU5EDT6pP0NTpdSjHtJ9jjjYxEqVyJXdZd5IKJco0jl8/iRe/i7rusCxmyZI6xl9Q5H29o5/20wMKbz15fe9DICXf0vIRd62++NYeJNYaUxzrKi2jbge/E5ChlqP61kaCdGZjdi0PYbhgPHpXBd6k6bh0TnN9g4z55CyPnOKlcxtmWAiJmVF0Wy9PnytH/3lk9Y+lssqcXDrHNR7f5/h1OQzJ4LkWXQO/1z0tvGe0eWgVPgMbt5Kr3CX70skVjI4eEbtWcRKTZVyj3IOfVCsFJGH1Dk5GeSin0olqKBcouhj9xHb2JCBNqBUipWAif5YTFlt1lf6hDSJLwA/zoHbzX3rzHTKUiCwiolg+QZMvgapxc1jsDuzEdvfkAQpZa1k6IhuHhpYGBlEbtqzmYg2Fy7ChH3GJB3zTPOiREku1uVjIWEa22TmQXP/5LipC6KaUlNvY8nbsZ8DtwHihkM3X1fr6joHDtKU3Df+muQcTN+4bKYqRMKIAwh8eAri0ANgOYfwTb2xbhv2HduMv3z+WOP1lVKcGNy2266Vc9A9eSWkzkG3JrHpIGycg37KbDJpPycxaWRAUHcJvVw/d+JkY3J7dWxs4h418Y0KfafsY61jM0v9wQfSu9TzT5gUHZuaT4mVDPWaNN1LW4dAg2mB9tQfqITItWAdMTFpaGDTEahcjqlVfV69bkjhmlRIu+eiuiBytdnAghbOQeKEvUeyOiN5HwF7SAwgfmanubsi7uP6Et03LF+3dQc+9eu5iXOSU3f1qavQ54hDrJBOQy5+7VWBvUcPwsAms9RNnYxWzkE7ddn70nFluOtNa5iM879Zs6JRA68ByZ2Kj225hGmxAmLOwSZykLmFm/tVUBVmZXRw7/jYtqOXBFGHTqx8llPbZzZicDpsw4ThscLbJIY5cPwQDFSUjabhPmT3YVF9IBAzcbm+gbT1mWnM1TF0EQd9o2KzgPIJKqnfb9Fq3uzSm3PQFlTZl1Q9LjQMGy4mee3Xz9jX6FwacA5us1f5zC7aXYo4B4dYKSzf3JLOIxEbyBTEoVsgd6jcjl7dibrixKhyVJvOQJ9EHztij3RbGaxMTLF+RmkJ0xOcQwZZdINlMW/w3MmUiLA9FKHYiIO6A3TrHNwfiM9u2+bDwp1S54HpPTWUS/jZRw+NftvyOxy15/CIg2ht60hxExL6YmJSqqrz2KWQ1k+bFqztbR144rV17D0S12vjYUqUUypRdG+7yXLwt93hv8C1YMtyKGFLHyujrbrDhVPYB5dYKfgbiZUcSn32XKSPKIhDt8DGOaiLvCtOjO9uyEfZbLO04OqqsutjJg/HP756Io7WYgapky4L52B7brnQu3ZXQR8DCmazflJvbVu0yUAQufu6YPeBSJclLW3cH3JQz3yPUinuQ4uFc9Dfy+n728NcA8kNi0+bpvEyiYd06MToywaveSAeR+vYUHox5LgEW6j4xD21uWcixEDsBOcyT5dNuszdY4W02ToSiMdDt7ZdeOmM2Hu6N3IORLSUiJ4nomeJaF5YtgsRPUhEi8K/w5T6FxHRYiJaSESnKeWHhe0sJqKrKUs4xYyweUirugOnTNFz0VDbvPQsQ6hi5nqbWGm7QhyEAHZXRB9x/+LjvHQOjYyznqmPLT6cg0ocrO2liQffrnva2HZhHIFKcpP2Pka9sIlOEHNBNs5BnV/nHjvRmEhHnacuzkF/11kTANnuPWXUIJx3/CRDTVUc456L7Q6dAzdm3LPoGxMTIZY38lmCI9Nvhwd55FAbEhHTGhGZsXckxXhNlbLiFNgLiUOIE4UQBwshpoW/LwTwkBBiCoCHwt8gov0AzAQwFcAMANcSkXxb1wE4D8CU8P+MHPrFIrZVYnYfygRzcw5+Qyd3L4ftMQwfPzItUsqKUokSnINprUuIGzJYK3GOOhL+nEOs0LTrHBSxklXnkDZlXfid9BRRo7pecOIk3MfkVLB9Z9ypiufiq56ziRwCy6sgx8CTS9YlIouqSHJ+lj4rnXaJD9Nh3T2Jg6Gar94NiL8Xu84h+PvMGxuiMu6dNJRLqRDi3PeqG4rYOQc/6hD5ODkc6+TY7HBYNKqmyBI/+tBBAOJn2pkU0mcCuCk8vgnAWUr5bUKIViHEEgCLARxORGMBNAshZotg63azck3ukLtDbsL7ehX7nJeQLzirg5Kp+totrYkFzrQTrlVhbnNMixT2zvAi8XGzpT3fIdGD3332nZPYndi4YbHfxB7DB2LfsZzDWu06B9tTk+diXqLAnv7hl1YBSCbpSdZTji0NHrRbbF5rcqiLO5n86T8nDQubyp06CJOcg64wKb7Q57TNk1rCShzIz0HS12G1FHEOdrGSHEO1vQGNlcS9eqVYCcH38gARPUVE54Vlo4UQKwEg/CsNi8cBUDNtLw/LxoXHenkKRHQeEc0jonlr1qypucMmZLH88CUOckHN4ssAAGccwCcWf+DFVYnfnziK50Z8TVn1D8hmoRWZ1nlwDhITLL4ivoHi9DzOXD7qoL342Ce3gA7unaq+GA9qY2+61ipWCkVkroXQN+REU6UccXvdzTkkzGgdz5NF5yDx3fcdgGYmWCRX12cuWQkT/Hxg5H1ufGyJtV5kyuqwVuJCfMiqvV0hfYwQ4lAApwO4gIiOt9TlRkdYytOFQlwvhJgmhJg2cuRIroobFp2D74Jqup6D/GCzWCTN/99TcdI+o531fjLzYLz7wF2d/cuS/EXNpKWjwcOUNbhffMzFp5Hw5hzCBXVKGC3208dNNLTn5pY+ymQ0A4D3HzKOdZBSQ0zYPtLk7exiPOHhRFVKEAdr1aiPLj+UlM7B8wUcwmQwA5Lvz3cz5aNzkDC9K64d165/6ACeyEiQn1QpWhdk6tGL370fW093gjMbmKQ5EVkm79UrdQ5CiBXh39UA7gJwOIBVoagI4V+ZkWM5gN2Uy8cDWBGWj2fKuwQ2U9YsC6ovMW+IZK329q768EHR8RDLRJ75jngIbSIg9WPNsovee0w6vaVEk6cpayJ0h81B0MszIfbt6BAC7zpwrJcznmmxGtBYYdM1fu8DB7DXqPH8bebNvtZKhIDIOO3kPa3NAEShnme/uta7TQCYNiGdF1rH187YB581KJqTIVrs7yRyHrOOTRaxknfV8L4O4kXupE5Ael14z4EGKzLNz8HocEjJekFfk33udcSBiAYS0WB5DOBUAC8AuAfAOWG1cwDcHR7fA2AmETUR0UQEiuc5oehpMxEdGVopfVK5JnfYTFnVCWTLJgb4vzC5ODonp6fc6YS9Y47JFmguy+IisR8jo1fRkMGUNbrG8hX70qzAt0Og3eJRrCOLKA0IxDMcpPwXsC9ISYW0GRRyDv/vD89bamnteb6/havssf9lKzKdqC0CL6K6zdYwG9GxY/pKw4wsnIMNnUmFysFmOKHCN0OfrCcNM0zEkxidgxyjcg9zDp0JvDcawF3hxK0AuFUI8TcimgvgDiI6F8AbAM4GACHEAiK6A8CLANoBXCCEkGY35wP4DYD+AO4L/3cJImslh1hpqyWbGGB2+NEhPwrf1KMuqH0cPsiciD2LslDCJduUOge3yaQf16IuenpoB709IcL8B55bRmvOCa8WAqhiJSvnoDyLiyD6OfQlr8kDcryzKH5tcyeLAUekc/Donw90InP/ArM+yHVfIANx0J7TFUpdWhaaY2MF5fe98FZUJqUbnEK6pa0D9y94CweOH2qN/ZYHaiYOQojXABzElK8FMN1wzWUALmPK5wHgnQByhi2fg/reXTJMX+LQELHTfjJZF9SP1Zdz8G3btWbFyvUMxMGTc5ixP6+Al/U6hEC7JRWlDhsNyWI4poqVrP4iniHSA4V0/Hv3XdI+KkCSAK/c0OLRUzdqyRTny/nZnB2BOHmVj4e0D7Ja/7ng6n9UT09l68k5mMaRew5pDRhHbI0nzKbtbfjCbc/i0rP273Li0Pc8pGFWEHWF7120oDqaroVz0PPbqkiYsno+l2tH6ysiU09bOQeFQJ9z1ARjPRlDv70qvP1L7JyD/3sekNA5+O2i7U6HlNh4/OITh/H1lDaypik13zv46yKwv/7UO6JjO6GLz9lCUwDxM/j6vfzzf060tpdVrKSn4NTRaDGcUKGbw7oSHN0QWjWZODD5yOqckWIkTiEtuQjfTVJn0PeIg1XnkP/9fK2VfCe77+Kofmi+LLOLOMjTmcRKNoW0KjpxLKhVESjtfEVkVmKbBZkjQAAAEQ9JREFU4T2rOz5bm+r7dXEOqgi52WiWm53zu+Q9vOWMem/AbdV03JQ4FIvNd0Ltlc2HQIVth66GpMhqeXXQbrxFlS/Ub2Rwvwru+s+j2Xp6v0xzQlcTmJ6bSyfaFvlGMMQh4iq6funue8RBHjh0Dk6HIk/4imJ84R3yQOm+SdmqwyVWksTDxYgkFNI5yP6lKWt7h/CyVAI6n2OD74cn52CNJ0UJRzyTgj1hbeZJEE805CqI+ujpkOlL3EsZOAcJ2453+45Yz+ciNvoznHkQb9LtC5U4fPHkvaIIujp0YwzTPNMNWkxzoqFMqWeVObQ5hXRb1Z0mIC/0OeIgV0BOtKDOtwFNfpP9fw12zhLyxeb1Kn2Jg/rx+E6k197mQzlITBo5CMdNGYHvM3kPVCQspTxFMa72OqpBjmtf8VsWr3BVjGKDTQ/lL0pLEuFBBnPkDmUX7fvMLmV9ZFbt+OrVdzaw0ayWVNt5YcVGa5vH7xVY2dnEN2rMsGEDzcYWeh+BznP9Az11S7oZt0kUrYsCbWle9Xwnus5BJUhyPcmaL7wW9Dni4GutZEpQY2rPhHvmBy4bf3/Jbk3hC99JoT5LXrqUxkoJt5x7hJuF97ydN3EoSbGSyGCt5P/Mx04ZYT1/66ePAABrZFh10XDpJlTxnSmP816jY38T32d2mfn66oxUjGo2Gz1MHhn3cZvDum9M2I6ti5JzeP+hbICEBFLhMzq5WO46NA69YmvLRbQkjp6UzGBn48p1fyWdACQ4hw67aWye6HvEwaJzUL+Zfz9momd7fjbIrlgsvvAmDj34Zn0DhWXxkJbIxc9BeftNlZJTVLX/+CEY3FTB56dPNtaZoi7mVuKQ3FCYCPdgJWxEXpxDFOcrw0JqExep0YB99Vo2TAo94I9zEGsgHebFNI5cfC0Oqge1Tex26O7DvN7H8EFNOHLP2Dy7qcE8PrpYqd1CHDoKhXTXIQ68lx7cWpS4Ljlv1vZcqMWqyRefOX7PzNdwaHME5pPwJw5+8nzTNTrUj81nnJr7NeD5b52Goye5Fy3AbeGTNVaO77t0zTE5dl1hleerkLbd+tDdh2HO16bjfYfwKXJVjFd2+oCZI7nXkupXhbpB2LbDbh2mOqLaIOfB4KaKVRKh6wQjzkFTSM9ftiHKm90dYqXOOMH1SsQhu9PIsqC++O3T0NYhjEHgdNRqR60jq85hrCEPAIeLztjXu64NvtyUr0mp+sy+OyabIlVNgbnDEXq5Frhs+aV46gvTp3i1d54n0R5oMW0G4o1FYxfsOq25EhCnsR3UZP9eRnmKc/XQ8ocaFMi+4iaVOOgpd3X4bn6kRdGRmohJh/zmxzT3w1ubWqJn0TmHM695nO1vV6HvEQdr4D3/dgZYFHUcfDkHXVapY7inzLNUIsz/5qnOBaMr4Ovt7zvevvJ8FTYOS8321d2hCQgUyY19Cd2u2i7Z2LZjcyMXlK5YWC48fR/r+S+dshda2qs4cZ8aA2ZqUNPZLr7sdG+O0gR187bnSLtz2Qtvbkxdw0G+jT0dzmqS6zpp31H49nunpqIqcOFqaok6nBV9T6wU/uWtlfIf8MveJx2/PU1QHedHDm7C4RN3wedONMu/JYb0b+j0R1MLfMUmvuOt7u7zECvlpf+pBSWKnyevd/O5EyfjHR5B9KS1VV4iThX7jrHL9vcd24yb//3wzJsqE1QT6TzGUSXUZx1sV4hLh7rp+9pFytKc1ZTFT0JyXZUSJZ4lIg4Mp5IlL3yt6IOcg1kj3RXEWMYMcslk5U7AZS9ORLjjM0fl07kuQt6cg7qY+yukzeN96B7DnGa7XYUSUfQ8ee3g//u0vb3q5X1fFd0hA1eR9+KYyJ/h+Sz/M8POLW1uCYiDK1y4XBtMCYx4zqGwVuoyuExZ80JLmGvWFl4bAI6ZNAKfP2myMc90b4IkwN97/wHWer6cgyoG8jZltSwePTnGROoi3b0LqpSlmxLodAbdIeZQkffYZZEavCd0uHPFNpK+Dq7xlpycTrRtmRd9nQ47gz7IOQR/O6uQ9sWeIwdit13645vvmWqtVyoRvnKq3w6w3hF5UufUnkocfBcF21rVHR+WCUQUmSp2x+5PxYgwx7YpeU9nkHcIbRf+f3v3HmNFdQdw/PvbZQGBhVWKlJdAraCIgrC1JcUuNcVKa0qjEjUqVJoUG41a07TaR9rUtNWGGlGbIFGI9mXTWFOstAbaoqW1VagiIiAopoJE8MnDCC7++secYWfvvXNn7u7ceez+PskNs2fPPcxvzu6ee86cOafW3k+SMS+5ZCqL51V/EBQ67mdFrVAQ1nPw048cPVo2yaN/lamxSel9jUO1zX7qcL0H9O3DP751bvIF55g/rJRUYxucXRT3D2oWOytGLRYHnXusUQ3dzXNOZd2ON7p7WsdcOG0UY4cOoHVc+PLoRVHLfYZNPzwv1n2WCcMHccao6EakoUHo1xD9AcMfDor6v/3Go7T35T8bcfiDD4+NQPis51AHafcceqO4azDFdaTTsFK8QqNuio9qOe7YVo9JGROy/HZQ8Pc/6r7HoraTWdRWeRe2rhCRHtEwQMe9pzg/Y80xh9Eeu6HaLse183sOUY3D6/u9JdlfefO9Tun+bKjD7R+ydtveTt+L+1xJd/S6ew7Vls/wk8KWNOipfjT3dO6d35pYef7wRdxnQKIEh5XizrQZ1hy+7APAmhvbunVOXfXi6wePHR94P96eIHl22dnh+zzXlftlbUnoZwy8xjPJGYvtblZa1JTXVZv2APDIxs67I/uzmA63H2Xnm50/SMRdTLM7el3PYc7kjzJh+KCKF9dvOFIePs3c/Cp7KXTFjbMncPKwgcyeNDyR8oI9h6gG55VbvxirTH8vjKiZJEnb/XZHb2XWhHhP12dh8bwpx3Yxq+anF54ROfGgHo4f0JfTRw6OPVMrC3F7DlPGtPDUzrf4bMmT137v4Certpa9x+451MHYoQND9871W/jMPg31EP2bGrnkE/GvYVQXOdg4TB45pMvnVWrtN2eF7qdQL8ENmj6XUONZDxdPj17CIktNjQ08et05WZ9GVf2aGjl05Ghkz+GeK6bz8DO7uXLG2E7pYTPAHvr6jLo8k1UqN42DiJwPLAEagXtV9da0z6Fvnwa23nJ+7KUuTPctvWJa5OJoi9o+xpMvvwkk+wDXuDpvs1jNiRHDXqb4po5p4W9b90beFzl+YF8Wzixf6LO0AThz9BBWXhtvragk5KJxEJFG4BfAbGAX8LSIrFTVF9I+lyynOfZG508eEZln1sT8Dr/Uyv91D9se1PQcSy6dytpt+2JNVIjjgjOjf1eSlIvGATgb2KGqLwOIyIPAXCD1xsHk05ob23pEj25ky3Fs33sw9TWdTPqa+zcde2AuCXHX2EpKXhqHUcCrga93AZ/M6FxMDn3crfVfdIvnTeG+dTtDt6E0JsifYPHCa/s5bURzRO5k5eWjWKVRubKPViLyNRFZLyLr9+3bl8JpGZOsYc39uGnOqamvRWSKbdLIwanchA7KS+OwCxgT+Ho08FppJlVdpqqtqto6bFgyS/8aY4wpl5fG4WngFBEZLyJ9gUuBlRmfkzHG9Fq5uOegqu0ici3wGN5U1uWqujnj0zLGmF4rF40DgKquAlZlfR7GGGPyM6xkjDEmR6xxMMYYU8YaB2OMMWWscTDGGFNGSrefKwoROQBsq/Ctk4D/xShiCPBuBvlqyWuxpJevlrwWS3r5asmbVSz1KLOesUxU1ejHrVW1kC9gfUj6vpjvX5ZFvhrLtFgsFosl57HU6frULZawv52lr544rPROzHyPZJSvlrwWS3r5aslrsaSXr5a8WcVSjzKzjAUo9rDSelUt29syLL2ILJZ8sljyyWJJtuwi9xyW1ZheRBZLPlks+WSxJFh2YXsOxhhj6qfIPQdjjDF1kvvGQUSWi8heEXk+kDZFRJ4UkU0i8oiIDHbpfUVkhUvfKCKzAu+Z7tJ3iMidkvbi6MnGslZEtonIs+6V+j6aIjJGRP4uIltEZLOIXO/STxCR1SKy3f17fOA9N7vrv01EPh9Iz7RuEo4l07qpNRYRGeryHxSRu0vKKlS9RMRStHqZLSIb3PXfICLnBspKp17iTtXK6gV8BpgGPB9Iexpoc8cLgVvc8TXACnd8IrABaHBfPwXMwNtY6M/AnALHshZozbheRgDT3HEz8CIwCfgZcJNLvwm4zR1PAjYC/YDxwEtAYx7qJuFYMq2bLsQyEJgJXA3cXVJW0eqlWixFq5ezgJHueDKwO+16yX3PQVWfAN4qSZ4IPOGOVwMXueNJwF/d+/biTQdrFZERwGBVfVK9q/sA8OV6n3upJGJJ4TRjUdU9qvpfd3wA2IK33etc4H6X7X46rvNc4EFVPayqO4EdwNl5qJukYknznMPUGouqHlLVdcD7wXKKWC9hseRBF2J5RlX9Dc82A/1FpF+a9ZL7xiHE88CX3PE8OnaR2wjMFZE+IjIemO6+NwpvtznfLpeWB7XG4lvhusffT7u7X0pExuF90vkPMFxV94D3C4HX64HK+4SPImd1081YfLmom5ixhClivUQpar1cBDyjqodJsV6K2jgsBK4RkQ14XbQjLn053sVaD9wB/AtoJ+Ye1RmpNRaAy1X1DOAc97oy1TMOEJFBwEPADaq6v1rWCmlaJT11CcQCOambGmIJLaJCWt7rpZpC1ouInA7cBizykypkq0u9FLJxUNWtqnqeqk4Hfos35ouqtqvqN1R1qqrOBVqA7Xh/ZEcHiqi4R3UWuhALqrrb/XsA+A0ZDWmISBPeD/qvVfUPLvl11/X1hyb2uvSwfcJzUTcJxZKLuqkxljBFrJdQRawXERkNPAzMV9WXXHJq9VLIxsGfaSAiDcD3gKXu6wEiMtAdzwbaVfUF1107ICKfct3J+cAfszn7zmqNxQ0zfcSlNwEX4A1NpX3eAtwHbFHV2wPfWgkscMcL6LjOK4FL3bjpeOAU4Kk81E1SseShbroQS0UFrZewcgpXLyLSAjwK3Kyq//Qzp1ov9bjLneQL79P0HuADvFbzq8D1eHf7XwRupeNhvnF4K7VuAdYAYwPltOL9QLwE3O2/p2ix4M3I2AA8h3ejaglupkzKsczE684+BzzrXl8AhuLdSN/u/j0h8J7vuuu/jcAMi6zrJqlY8lA3XYzlFbyJEgfdz+WkAtdLWSxFrBe8D4qHAnmfBU5Ms17sCWljjDFlCjmsZIwxpr6scTDGGFPGGgdjjDFlrHEwxhhTxhoHY4wxZaxxMKYORORqEZlfQ/5xElit15is9cn6BIzpaUSkj6ouzfo8jOkOaxyMqcAtjvYXvMXRzsJ7SHE+cBpwOzAIeAP4iqruEZG1eOtffRpYKSLNwEFVXSwiU/GefB+A9+DSQlV9W0Sm462h9R6wLr3ojIlmw0rGhJsILFPVM4H9eHts3AVcrN5aWMuBHwfyt6hqm6r+vKScB4Bvu3I2AT9w6SuA61R1Rj2DMKYrrOdgTLhXtWNdm18B38HbeGW1W/G5EW85FN/vSgsQkSF4jcbjLul+4PcV0n8JzEk+BGO6xhoHY8KVri1zANhc5ZP+oRrKlgrlG5MbNqxkTLiTRMRvCC4D/g0M89NEpMmttx9KVd8F3haRc1zSlcDjqvoO8K6IzHTplyd/+sZ0nfUcjAm3BVggIvfgrZp5F/AYcKcbFuqDtxHT5ohyFgBLRWQA8DJwlUu/ClguIu+5co3JDVuV1ZgK3GylP6nq5IxPxZhM2LCSMcaYMtZzMMYYU8Z6DsYYY8pY42CMMaaMNQ7GGGPKWONgjDGmjDUOxhhjyljjYIwxpsz/ASpQurFCxL3QAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "sorted_data['inc'].plot()"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 15,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 15,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEVCAYAAAALsCk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXl4Y3d57z+vLEvyInm3Z988njUbM5NJJmSDkAVoCaRwG0pJStMb4NLSlvY+QJfLbW+BUlpKaRuWsCSBFhIKNGFNQhLIvsxM1pnJzNizehYv40XeJGv53T/OObJsy7Zky5ZkvZ/nmcf2T+cc/eyR9D3vLsYYFEVRFCUZV643oCiKouQfKg6KoijKJFQcFEVRlEmoOCiKoiiTUHFQFEVRJqHioCiKokxCxUFRFEWZhIqDoiiKMgkVB0VRFGUS7lxvYLbU19ebNWvW5HobiqIoBcWePXu6jTENMx1XsOKwZs0adu/enettKIqiFBQicjyd49StpCiKokxCxUFRFEWZhIqDoiiKMgkVB0VRFGUSKg6KoijKJFQcFEVRlEmoOBQh8bhBJwAqijIdKg5FyLvueIp/eeRwrrehKEoeU7BFcMrsOdw5yLLqslxvQ1GUPEYthyIjEoszPBqjd3g011tRFCWPUXEoMgZCUQB6hyI53omiKPmMikORERyxREEtB0VRpkPFochIWA7Do5qxpCjKlKg4FBnBkGU5RGKGodFYjnejKEq+ouJQZDhuJYDeIXUtKYqSGhWHIsOxHEDjDoqiTI2KQ5ERHIkmvu9Ry0FRlClQcSgy1HJQFCUdZhQHEVkpIo+JyAER2Scif2yv14rIwyJy2P5ak3TOJ0WkVUQOisj1SevbReRV+7EviYjY614Ruddef05E1mT/V1XAijmUuATQWgdFUaYmHcshCvyZMWYzcCnwERHZAnwCeMQY0wI8Yv+M/djNwFbgBuAOESmxr/Vl4Hagxf53g71+G9BrjFkP/DPwuSz8bkoKgqEoSwI+XKKWg6IoUzOjOBhjzhhj9trfDwAHgOXAjcDd9mF3A++0v78R+J4xJmyMOQq0AjtFZCkQMMY8Y6wE+3smnONc67+AaxyrQskuwZEIVWWlVJd7VBwURZmSjGIOtrvnDcBzQJMx5gxYAgI02octB04mndZury23v5+4Pu4cY0wU6AfqUjz/7SKyW0R2d3V1ZbJ1xWYgFCVQ5qamvFTdSoqiTEna4iAilcAPgD8xxgSnOzTFmplmfbpzxi8Y8zVjzA5jzI6GhoaZtqykIBiKEPCVUqOWg6Io05CWOIhIKZYw/Icx5of2coftKsL+2mmvtwMrk05fAZy211ekWB93joi4gSqgJ9NfRpmZ4EiEQFkpNRUeTWVVFGVK0slWEuAbwAFjzBeSHnoAuNX+/lbg/qT1m+0MpLVYgefnbdfTgIhcal/zlgnnONd6N/Co0cY/80IwFLUth1K1HBRFmZJ0hv28EXg/8KqIvGSv/QXw98B9InIbcAJ4D4AxZp+I3Afsx8p0+ogxxmni82HgLqAM+Ln9Dyzx+baItGJZDDfP8fdSUhCNxRkMWzGHUrfQOxzBGIPG/hVFmciM4mCMeZLUMQGAa6Y459PAp1Os7wbOS7EewhYXZf4YDFvV0X5fKWWlJYxGrcE/FV4dCKgoyni0QrqIcFpnBHxuaso9wPzUOnQOhPjYvS8xPBqd+WBFUfISFYciwmmd4QSkYX6qpH91sIsfvniK/aenS2pTFCWfUX9CEeG06w74SiktsVtozIPlcOLcMAADYbUcFKVQUXEoIoL2FLhAmRtfqdXRZD7E4di5IQAGQyoOilKoqFupiEi4lewiOJifgT8nemzLYRpx6Bka5aY7nuLZI+ey/vyKoswdFYciIuFWKiulqqwUEegZzn7M4bjtVhoMT33tz/38dfae6OPxQ9oGRVHyEXUrFRHBUBQR8HvduFxCVVlp1i2H/uEI/bYITeVW2nO8h3t3W+23jnYPZfX5FUXJDioORURwJEKlxxIGgNp56K90vGfswz6YQhyisTh/9d/7WBLwsba+QsVBUfIUdSsVEcGQ1VfJoaZiHsTBdinBWNFdMvc8c5wDZ4J86je3sHVZgGPnhojHtVOKouQbKg5FRHAkit83ZizOR9tuJxi9srZsklspGIrwhYcPcdWGBm44bwlrGyoIReKcDYayugdFUeaOikMRMclymAe30rHuIRr8Xhr9PgYmBKRPnBtmMBzlvTtXISKsrasANO6gKPmIikMRERyxZjk4zItbqWeY1bXlVHrdkywHJ1BdU27tYW2DioOi5CsqDkWEMwXOoabcQygSZ2Q0Ns1ZmXHi3DCr6yrw+9yTKqT77LTZKlscmvw+fKUuFQdFyUNUHIoIZwqcg3MHny3rIRSJcTYYYnVduSUOU1gOVbZry+US1tRVcEzFQVHyDhWHIiEeN/Ysh/FuJSBrE+FO2sHo1XXTu5Wqkvag6ayKkp+oOBQJA+Eoxljtuh2y3bbbSWNdVVuO31fKSCRGJBZPPN4/EqG0RCiz+zqBJQ4neoaJJh2nKEruUXEoEpI7sjrUVjhupdmls/7itTO8+R9/lbAInIZ7q+sqqLQHCA0lxR36RyJ2246x2VFr6yuIxg3tvSOz2oOiKPODikORMDbLIYXlMEu30s9ePcuR7iG+8+xxwKpx8Pvc1JSXUmlbKMlxh6AtDsmsrdeMJUXJR1QcFjGhSAxjrOrjsSlwYx/Ozgf1bNxKxhheONYDwLeeOkYoEuP4uWFW15UjIgn3VbI49Ks4KErBoOKQx4yMxnhuli2tQ5EYl372Eb77vNXgLnkKnIO7xDXr5nvtvSOc6Q9xw9YldA+G+cHedk70DLO61vqwr/Raz5PcQqNvZHSSONRWePD73CoOipJnqDjkMT968RQ33/ksXQPhjM890x+ibzjCT145DaSOOYCVznpuFuKw+7hlNXz0mhYuWFHF1x4/QnvvMKvqygESbTqS23anshxEhHWasaQoeYeKQx7TOzyKMdAxi95DZ/utc54/2sNAKDJuClwy6xv9vHqqP+PrP3+0F7/PzcYlfj50VTPHzw0TiRlW11rikCrm0D88WRwA1qg4KEreoeKQIXuO9xKOZq+ieDqGR60P1s6BzMXBEZRo3PDk4W4GbLeSk0XkcFlzHcfPDSdqFNLlhWM9bF9dQ4lLuH7rkkTsIGE5eMeLQzxuGAhHU4rD2voKTvePEIoszN9VUZSZUXHIgN6hUd79laf57M9eX5DnGwpbH5azcSs5nU7LPSU8drCT4EiUSq8bd8n4//LLW+oBeLqtO+1r9wyN0to5yMVragEocQl/+Kb1eN0uNjT5AfD7xsccBkJ2ncUU4mDMWEdXRVFyj4pDBpwbstw8//ncCdp75/+DLGE5BGchDv0hKr1u3rSpkccOdtE3MjquXbdDS2MlDX4vT7VOHfjuGgjzWpLrabedpbRzbW1i7be2r+DF/3Mt9ZVeAHylLkpckrBYnFqIajt9NhnH6jjSpa4lRckXVBwywMn4GY3F+eIvD8/78w3ZDfG6BjMXh45giKaAlzdvbKRrIMxzR3omBaPBCghf1lzH023nEmmvE/nCw4e46Y6n2X86CFguJU+Ji/OXV407rtwzJj4iMq6FRt+IFfSeKuYAY0V0iqLkHhWHDHAyfnasruGHe9s53DEwr883HJ6D5RAMsaTKx1UbGxCBU30jk4LRDm9srqd7MMyhjsGUj5/oGWI0FuePv/cioUiM54/1cuHKKnxJbTBSkdyZNVVfJYeArxSv25W1Hk+KoswdFYcMcDJ+Pv7WTZR73PzTQ4fm9fkSMYfZWA79IZoCPuorvVywohqYnMbqcNn6OgCeak0ddzjdF2JVbTmHOwf51P372HeqPxFvmI5K71hn1unEAaxYhCO+iqLkHhWHDHA+vFbXlvM/r1jHL/adTfjf54MhO+aQaUA6Hjd0DoRZEvAB8OaNjUDqYDDAippyVteVpwxKG2M43TfCDect4QNvXMO9u08SjZu0xCHgK024lWYUB5874bZTFCX3qDhkQHKV8W1XrGVZlY8PfOsFHj/UNS/PN2zHHDoHQlPGA1LRPRQmGjcsqbLFYZMtDikC0g6XNdfz7JGeSd1Re4ZGCUfjLK3y8fEbNrGxyY9LYNvqmhn3UelzJ0aFpmM5TJz/oChK7lBxyICBUBRPiQuv20Wl180P/tdlrKgt5wN3vcB/Pnci68/ndDQNReLj2lDMREe/ZWk02ZbD1mUBrmip5+K1U9/tX76+nsFwlJfbxxfEne6zUmKXVZfhKy3hmx+4mK/87vYpP+STSQ5I949E8Lhd+EpTv+QCPnUrKUo+oeKQAcGRCIEyd6Ll9NKqMr7/oV1c2VLPX/zoVb6/+2RWn294NEadPZCnMwPXklPj4LiVXC7h27ddwm9csGzKc3Y1W3GHpyfEHU71Wa20l1eXJb5et3VJWvvw+9wJUXOqo5PbdU88NqiWg6LkDSoOGRAMRScFdSu9bu68ZQcNfi/PHc1e/MEYw9BolNV2xXEmcYeEONhupXSorfCwZWmAZ4+Or3c4bYvD0gyu5VCZ9IGfqq9SMhqQVpT8QsUhA4IjEfwpPuDcJS4a/d6spmKGInGMGasByEQcOvpDlLgkUZCWLluWBSals57pH8HrdlFbMbl4bSb8Xjej0TjhaGxmcfCV2lXU6cdWFEWZP1QcMiAYikwZ1K2r9HJuFimnU+FkKq2ps8QhU7dSQ6WXEldqF85UrG+spGsgnAgegxVzWF5dNqU7aDqcFhpD4TTEoczNaCxOOKrjQhUlH1BxyAAr5pD6A66uwjOr1tdTMWzXOCyrLqO0RDKzHIIhmmbhBlrfUAlAa+eY9XCqb4RldrwhUyoTzfciaVkOgLqWFCVPUHHIgFQxB4e6Cg/nBrMnDo7lUOktoaHSm1Fn1rP9IZYEMnMpgWU5ALQlicOZ/pFZxRtgbKbDQCg6ozg4x2qtg6LkByoOGeBkK6WittLDSCTGyGh22k47aazlHjcNAV/GloOTqZQJK2vL8bhdtHZZ4jAajdM5EJ695eB84I9EGAilbtft4FhkmrGkKPnBjOIgIt8UkU4ReS1p7f+KyCkRecn+97akxz4pIq0iclBErk9a3y4ir9qPfUlsJ7aIeEXkXnv9ORFZk91fMTuEIjHC0fi0lgPAuaHsxB2cpnsVtuWQrjiMjMYIhqI0zkIcSlzWVDbHrdQRDGHMWBprpvjtUaFOOqy6lRSlcEjHcrgLuCHF+j8bYy6y//0MQES2ADcDW+1z7hARpzvbl4HbgRb7n3PN24BeY8x64J+Bz83yd5lXBhKT1KYSB8uNky3XktN0r8LrpjGQvjhMrHHIlObGyoQ4JNJYq+fmVmrvnVkcqsoct5JaDoqSD8woDsaYx4F0E/hvBL5njAkbY44CrcBOEVkKBIwxzxgrV/Ee4J1J59xtf/9fwDUym9SYeSbROmOKbKXaSstyyFY6a8Jy8LhpqPTSMzxKJDZzJo8zHjSTGodkmhsqOdk7TCgS43S/9aE+V7eSWg6KUnjMJebwhyLyiu12chrtLAeSy4Tb7bXl9vcT18edY4yJAv1AXaonFJHbRWS3iOzu6pqffkZT4XxoTeVWqrcth+4spbM6g37KPSU0+L0Yk55V4owHbZql5bC+sRJjrME7idYZVXPLVjrlWA7l0wWknZiDioOi5AOzFYcvA83ARcAZ4J/s9VR3/Gaa9enOmbxozNeMMTuMMTsaGhoy2/EcCSbcSgtkOYSdmIObRr8lPOm4lmZTHZ2Mk87a1jXIqb4Rais8lHmmn9swFb7SEjwlLtr7rKl501kOvlIXpSWizfcUJU+YlTgYYzqMMTFjTBy4E9hpP9QOrEw6dAVw2l5fkWJ93Dki4gaqSN+NtWDMZDlUeErwZHFgzfBoFJeA1+2iwRaHdNJZnfGgzl17pqxrqEDEqnU40zf7NFaHSp+bM7YFMp04iIg231OUPGJW4mDHEBzeBTiZTA8AN9sZSGuxAs/PG2POAAMicqkdT7gFuD/pnFvt798NPGrysIdCcrvuVIgI9RUeurMUkB4Kx6jwWE3+nMyjdCwHZzzobPGVlrCyppzWrkFO94VmHW9w8PvcROPWf+dMnVwDZaUakFaUPGHG20sR+S5wNVAvIu3Ap4CrReQiLPfPMeCDAMaYfSJyH7AfiAIfMcY4if8fxsp8KgN+bv8D+AbwbRFpxbIYbs7GL5ZtgiO2W2kKywEs11JPllJZh0ejlHstd059ZfqdWZ3xoHNhfWMlbZ2DnO4bSXRrnS2OBeN1u2YcKxrwudVyUJQ8YUZxMMa8N8XyN6Y5/tPAp1Os7wbOS7EeAt4z0z5yTTAUobREppxHAFY6a7ZaaAyGo1R4nA/WEqrLS9OzHPpDXDrHD/T1jZU8fqiLaNzM3a1ki0M68x/8vlIGNCCtKHmBVkinyUAoQsA39TwCyG4LjeHRWMJyANIqhJs4HnS2rG+oTLiC5u5WskQhHXEIlOlMB0XJF1Qc0iQ4Ep0y3uBQW+HJXoV0kuUA0BgY668Uj5txzfEcugbHjwedLc12jyXIhjikbzloQFpR8gcVhzSZrl23Q12ll1AknqhRmAvDozEqkjKOGiq9dA2GiccNn/zhq7zlC7/muSPjB/O8dsoa8bmxyT+n53bSWQGWzbI62iETt5IVkFZxUJR8QMUhTaZr1+2Q6K+UBdfS0GiU8qT6gga/l85gmL++/zXutceRPt02Xhz2nujF7RIuWFE9p+euKi+l3p4H0eifmzgkLIdpCuAcAj43oUicUZ3poCg5R8UhTaZr1+1QV+k035u7OAzbqawOjX4f4Wic/3juBB+6qpnNSwPsPj6+HGTP8V42Lw3MumgtmfWNFSwJ+DIeGDSRygzcSk58QoPSipJ7ZlcpVYQERyKJu+CpcEZpZiOddSgplRXGfP9/cPlaPn7DRkZGo3x/TzuRWJzSEhfRWJyXT/bz2xevnOqSGfEnb9lAbxZEzp+RW2ms+V5dhiNOlfklFInREQyx2p5MqCx+1HJIk2AoHbeS018psw/VeNwQTWqqZ4yxYg5JlsO1W5q474O7+Mu3b0ZE2LGmluHRGAfOBAF4/ewAI5EYb1g1N5eSw6Xr6njr+UtnPnAGMspW0uZ7ect3nj3ODV98gnA0O/NKlPxHxSENwtEYoUg8jYD07Porfe7B17n+i48nPV+cWNyMsxw8bhc719YmUml3rLF6Hb5wrBeAF09YX7etqiGfyDQgDdp8Lx85dm6IkUiMvmH9vykWVBzSYKZZDg7lnhK8GfZXOjcY5u6nj9HWNZS4K3OmwCVbDhNZWlXGipoydh+z4g57T/TR4PeyomZuqafZxglE15R7Zjw2kIg5aK1DvtERtFylKg7Fg4pDGszUdM9BRKiv9GbUtvvup48RilgupU77DThsz3IonyGwfPGaWl441osxhj3He9m2qnraIr1csG1VDf/vxq28cX39jMf6k8aKKvmF0wq+bzh7c9KV/EbFIQ1matedTG2FJ23LYTAc5e5njif1TrLegEN2ncRMnVV3rKmhezDM3hO9nOgZzjuXElijR9+/aw0e98wvNXUr5S8JcVDhLhpUHNIgXcsB7CrpNAPS33v+BP0jET527UZgzHJwZjmUzyAOO9fUAvDVXx8BYPvq/BOHTKjwlOCSsSaHSn4Qi5tE6xa1HIoHFYc0mKlddzJ1lelZDuFojDufOMKudXVct7UJGLs7cyqsK2ZwKzU3VFJdXsrDBzooLRHOW1414/PmMyJCoEyb7+Ub3YNh7FZbGnMoIlQc0iCddt0OdXZ/pZlGUtz/4mk6gmE+fHUzteUe3C5JtOROWA7TBKQBXC5hx+oajIEty6pmbIldCAR8OtMh33BuWkDdSsWEikMajFkOM8ccxvorTZ8P/p/Pn2Bjk58rWupxuYRGvzeREZKwHLwzf9jvsF1L27JU35Br/DrTIe9wXpegbqViQsUhBaf7RrjyHx7jcMcAYMUc3C6hLI0787Eq6anfRCd7hnnpZB/v2rY8kV3UEPAlBaTTsxwAdq2zZjdcsrZ2xmMLActyUHHIJ5y55AGfW91KRYSKQwoOnh3gRM8wP3rxFGBZDn6fO600Uaf53nTprA+8bI3PfntSBXKT3VgPYDicvuVw4cpqfvJHl3P91iUzHlsIBMrcGpDOMzqDIVxiDYFScSgeVBxS0Gubzr880AFYRVnpBKOBRE+g6SyHH798mm2rqllZW55Yawx46UiyHETA504vhnDe8qq8q2+YLQGdBpd3dARDNPi91FZ4NeZQRKg4pKDXvjs61DHI8XNDVrvuNILRMHPb7sMdA7x+doB3XLhs3HqT30ffcIRwNMZQOEp5aQmuOXZELUSsmQ5qOeQTZ4PWdMHq8lKNORQRKg4p6E96Azy8v8Nq151GMBrGYg5Tte3+8cuncQm87YLxTe2a7NGencEww6PRcYN+igm/z81gODquEaGSWzqDIRoDPmrKS9WtVESoOKSgdzhCTXkpG5v8/PJAR0aWQ7mnBF+pK2XbbmMMP37lDJeuq5s0RKchYLmjOgdCDIVjRSsOzt95MKzWQ77QEQzRFPBSXe5hJBIjFNHOrMWAikMKeodHqS738JYtjbxwrJez/aG0xUFEqKvwpmzbve90kKPdQ5NcSmC5lWDMcpipr9JixYntaPO9/CAUidE7HGFJwJforNuvcYeiQMUhBf0jEarLS7l2yxJiccNAOH23EkC935toN5DMj18+jdsl3HDe5MyiJtty6AjalkMaaayLEactun4A5QfO69hyK1kuU3UtFQcqDinoHR6luqyUC5ZX0ei3PrTTtRwAGv3eRM1CMi+e6OOildVUp2hfXZNUJT08YQpcMaHN9/ILp8ahyQ5IgxbCFQsqDinoG45QU+7B5RKu2Wz1PUo3lRUccZhsOXQMhBLjPieSXCU9NFq8lsNY2251K+UDTuuMZLeSprMWByoOKegbjiTu7q/d0gik1zrDoTEpLdXBGMPZ/lDCfZTyPLtKejhcxDEHn1oO+YTTOsMKSKvlUEyoOExgNBpnMBxNvBGubGngr39jS8KCSIfGwORZ0sGRKOFoPJGymvI8u0p6aLSIs5U0IJ1XdARDeNwuqspKNeZQZKg4TMAJhNbY4uAucXHb5WszjjmAlR/u4FQ/TycOTQEfHQMhqwiuSC0Hv9eNiE6Dyxc6giGWBHyICOWeEkpLRN1KRUJx3p5Og2MyV6Ux83gqnBqG5LhDR3BmcWj0exN3ZcVqObhcQsBXytHuoVxvRWGsxgGsNO2qMo+6lYoEtRwm4LTOcCyH2dCYKGhLFocx3+1UJAvHTIN+FjO/tW0FP3nlNAfPDuR6K0VPRzBMY9LrUqukiwcVhwk4d0U1c7Ac6io8iEBXslvJ/n5iZXQyjUnCMdOI0MXMR69ZT6XXzWd+diDXWylqjDEJt5JDtYpD0aDiMAHnhV+VQerqRNwlLuoqvJPcSgGfm7JpLIJk4SjWVFaA6nIPH72mhV8f6uLxQ1253k7RMhiOMjwaG2ftVpV5NOZQJKg4TKBvxLYcKmZvOcDkWoeOYIglVVNbDTDe5VSsRXAO79+1mlW15XzmZweIxacfuarMD6niZNqZtXhQcZhA77A19W2uPv/GwPgq6Y5geNpgNFiurNISq013MVsOAF53CR+/YROvnx3gB3vac72domQsTqYxh2JExWECfXbTvbkOz2lMmuwGdtvjaeINYGXqNNjDgoo1lTWZt52/hPWNlfz4ldO53kpRktpy0M6sxYKKwwT67Hbdc6XB76V7MEwsbojHDZ0D4WkzlRyczJBiTWVNRkS4aGU1B85o1lIuOJtIokiOOWhn1mJBxWECVrvuuYtDo99H3FjjQs8NjRKNmxndStZ51huxmFNZk9m8NED3YDhlI0NlfjnePUx1eem4GxWtki4eZhQHEfmmiHSKyGtJa7Ui8rCIHLa/1iQ99kkRaRWRgyJyfdL6dhF51X7sS2L7bUTEKyL32uvPicia7P6KmZHcV2kuJKqkB0JpFcA5OMcUcyprMpuX+gHUesgBe0/08oaV1ePWJvZXMsaw73T/gu9NmX/SsRzuAm6YsPYJ4BFjTAvwiP0zIrIFuBnYap9zh4g4t8BfBm4HWux/zjVvA3qNMeuBfwY+N9tfJhv0DUeonkMaq0NyIVxnonXGzG6lLcsCNPq9lJeq5QCwZWkAgANngjneSXHRPxLhcOcg21bVjFuf2Jn1gZdP8/YvPakFi4uQGcXBGPM40DNh+Ubgbvv7u4F3Jq1/zxgTNsYcBVqBnSKyFAgYY54xxhjgngnnONf6L+AamWs0eA70Do/OOY0VxmoWuoLhlFkfU3HzxSt56hNvxuXK2Z8gr6gu97CsyqfisMC8eKIXgO2rx4vDRMvhkQOdgIr3YmS2MYcmY8wZAPtro72+HDiZdFy7vbbc/n7i+rhzjDFRoB+oS/WkInK7iOwWkd1dXdkvjgpFYoSj8azEHBpSuJUa/DNbDiJCaYmGgpLZsizA/tP64bOQ7D3Rh0vgwglupeSYQyxueOKw9T5s6xpc8D0q80u2P4VS3e6aadanO2fyojFfM8bsMMbsaGhomOUWp6bXvhuqLpu75eArLSHgc9M5YFkO9ZUe/dCfJZuXBjjSPaTpk1MQisR4z1ee5pm2c1m75t7jvWxaEpiUNZfcmfW1U/2JXmStnSoOi43Zflp12K4i7K+d9no7sDLpuBXAaXt9RYr1ceeIiBuoYrIba0HoHZp7071kGgM+OoNhu7PlzC4lJTWblwaIxQ2HOtSvnYpDHQO8cKyXZ45kRxxiccOLJ3onuZRgfGdWp7XJtlXVKg6LkNmKwwPArfb3twL3J63fbGcgrcUKPD9vu54GRORSO55wy4RznGu9G3jUjkssOE7rjKpsiYM9S1rFYW5oUHp6nGBwV5bSfQ91DDA0GkspDjBWJf344S7OWx5g59o6jp0bIhqLZ+X5lfwgnVTW7wLPABtFpF1EbgP+HrhWRA4D19o/Y4zZB9wH7Ad+AXzEGOP4Aj4MfB0rSN0G/Nxe/wZQJyKtwMewM59yQV+iXffc3Uow1l/Jap0xc7xBSc3ADKk0AAAgAElEQVSq2nIqPCWazjoFjjh0BCfPLZ8Ne45bweiJmUoO1eWlnOwdZu+JPq7a0EBzQwWRmOFEz3BWnl/JD2ZMpjfGvHeKh66Z4vhPA59Osb4bOC/Fegh4z0z7WAh6s9CuOxnHrRSJx2dsnaFMjcslbFqqQempOGi727JVKLj3eC/1lV5W1palfLyqzMMLxzoAa4yux23dY7Z2DrKuoTIre1Byj0ZIk3Ash2xkK4FlOYzG4hiTXhqrMjWbl/o5cCZIjjyOec1cLYenWrv54d6xZMK9J3rZvrp6yv5izvuj0utm2+oamhstQWjr0ul9iwkVhyT6hkfxlbrwZakALTl1dUmVupXmwualAQbCUdp7R3K9lbyid2iUzoEwlV433YPhjP3+o9E4f3rvS3zsvpf57M8O0DUQ5ti54SldSjCWsLGruY7SEhcBXymNfq8GpRcZKg5J9A5HspLG6pDsSlK30txwgtL7NSg9DseltKu5DmPg3FBmsxZ+9uoZOgfC7FpXx1cfP8IH7noemFz8lozTXubKDWPp5OsbK2nVWodFhYpDElZfpey4lGC85aBupbmxcYkfEc1YmojjUnI+qDuC6ccdjDF866mjrGuo4D/+4BI++ub1vHYqSGmJcN7yqinPW1FThtslXD1BHI50DqrbbxGh3d2S6BsezVowGsb6K5W4hLostOQoZso9btbWV/DaKRWHZA52DFBVVsr59od5ZwZxh70n+ni5vZ//d+NWXC7hY9dtZEVNOX0jo9O6Vt9+/lK2raphZW15Yq25oZKBcNRuTa83QouBoheHWNxQYvcx6huJ0NKYvWwLv9eNr9RFTblHeyVlge2ranhof8e4/7Ni5+DZATY2+ROp0h0ZZCx986mjBHxubto2Vp/6Py5eOc0ZFu4S1zhhAMtyACtjScVhcVDUbqVHX+/gor99KNEXxpkCly1EhEa/T98sWeKN6+vpH4moa8nGGMOhswNsXOKnvtKLSPqWw+m+EX7x2llu3rkqK4OlksVBWRwUtTg8fqibgVCUz/z0AMaYrMccAK7b0sQ1mxpnPlCZkcuarX6MT7V253gn+cHp/hAD4SgblvgpLXFRV+FNu9bhnmeOY4zhll2rs7KXRr+XSq9bG/AtIorarfTqqX5KXMIjr3fy89fOEo2brPVVcvir39iS1esVM40BHy2NlTzVdo4PXtWc6+3knINnLQtq0xJrIFKj35tWrcPxc0Pc9fRR3nb+UlbUlM94fDqICM2NlWo5LCKK1nKIxuLsPx3k5otXsqq2nP9z/z6ArLqVlOxzWXMdLxztYTSqfXwOnrU+iDc0WeLQFJjZcjDG8Mkfvkqpy8VfvT27Ny7rG1QcFhNFKw5tXUOMRGLsWFPDX7xtM92D1h1XNqbAKfPHZevrGYnEEsNoipmDZ4MsrfIlprM1+n0zWg7f393O023n+MTbNrGkKruxsObGCjoHwgRDOl96MVC04vBKex8A5y+v5vqtTexaZ/mzszEFTpk/Ll1Xh0vgqSzOLihUDnYMstF2KYFlOZybpkq6cyDE3/10PzvX1vLei1dlfT/r7b5Kr2uDxEVB0YrDq6f6qfCUsK6+AhHh/71zK2/Z3Jjw3yr5iZPT/0xbcQelI7E4bZ2DbGwae702BHzEp6mS/vRPDxCKxvnsTefPS2r1pc11lJWW8IM97TMfrOQ9RSsOr7T3c97yqsSbZH2jn6/fejF+n7qV8p3L1tfz4ok+hsLRXG8lZxw/N8RoLJ6INwA0OaNpU7iWjDE8+nonN71hOc3z1Dk14CvlnW9Yzv0vn6J/WF1LhU5RikMkFufAmSAXrJi6RYCSv7yxuZ5o3PD8sZwMDMwLWjutDqgtTWMf9E49TaoWGp0DYQZCUTbbParmi9+9dBWhSJz/2qvWQ6FTlOJwuGOQcDQ+bf8YJX/ZsaYGj9vF00Vc73Ck28oKSp6f0DhNlfThDuv4bHYASMXWZVVsX13Dd549TjyufZYKmaIUh1dPWcHoC1ZU53gnymzwlZZw0cpqXjhW2BlLZ/pHeGjf2Vmd29Y5RFPAKjxzmK5K+nCnFSRe3zT/w3jef+lqjnYP8bQmDRQ0RSkOr7T34/e5WV2bnQIgZeFZUVNG10B2xmLmim8/c5wPfmfPrGInbV2Dk2IHVpW0J2WtQ2vnIFVlpTRUzv9ckbeev4TaCg/ffvbYvD+XMn8UpTi8eqqf85OC0UrhUVfh4dxQYYtD7/AoxpBxywljTEpxAKvWIbXlMEhLY+WU092yidddwm9fvJKH93dwuk+HMxUqRScOo9E4r58Z4HwNRhc0tRVeQpE4w6OFm7HkjKV14gHp0jVoBZebGyomPdYY8KaMObR2Diaa4y0EN160jLiB548Wb9JAoVN04nCoY4DRWDzR/14pTJz5GOcGM5t8lk/0j1jicKgzs6KxNjtTqTnFh31TCsvh3GCYnqHRBRWH5dVlAAXv+itmik4cXmnvB+CC5RqMLmRqbXHoyXAsZj7hiENrhpaD44ZK5VZqCngnzZJ2+h21NC1cgWelPcsk3S6xSv5RdOKwtNrHOy9axsraslxvRZkDdZWFLw4Jt1KGzeraugYp95SwJMWckFRV0s715zuNNRlnlolaDoVL0bXsftPGRt60UecrFDp1FVbWjdMwsRAJjkQQgZO9w4yMxijzTD2aM5m2riHWNVSkTKhIrpJ2iuJaOwep8JSwNMuN9maiwe+lU8WhYCk6y0FZHNQWuOUQixtrUE+jP+OMpbbO1JlKYM28gPFV0oc7B1i/QJlK4/ai4lDQqDgoBUmFpwSP21Ww4hC04w071tQAVqKEw3NHznHZZx/h608cITKhw+rIaIxTfSNTioMzSzr5Q/lwxyDrGxe+oWSD36tupQJGxUEpSETErnUoTHHos8XhwhXVlJbIuLjDfbvbORMM8Xc/PcDbv/QEzyRVGo+1zZicxgpWlXRpiSRa0vePROgcCI/rwbRQNPq99I9ECEViC/7cytxRcVAKlrpKT8FaDk6mUl2lh7X1FYlah2gszqOvd/DOi5Zz5y07GB6N8d47n030kTrSZaexTmE5lJa4+J2dq7hv90lePxscy1RawGC0Q4O/8ONCxYyKg1Kw1FZYw21yyed+8To/feVMxuc54lBVVkpLoz/R+2jviT56hyNcu6WJa7c08dCfXsny6jI+/bMDxONWZbQIrK1PbTkA/Om1G6gqK+VvHtjPYdtdtZA1Dg6Nfiv+oXGHwkTFQSlYcu1WisUN33jiKHc/fSzjcx1xqC4vpaWpkhM9w4QiMR7efxZPiYsrNzQAUO5x8+fXb2Df6SAPvHyatq4hVtSU4SudOrOputzDx67byDNHzvH1J4/idbtYUbPwfcQcy0HjDoWJioNSsNRW5Nat1N47zGgszsvtfYxGU4/mnIr+YWvfAdtyMMZKOX14fweXNteN67Z644XLOW95gM8/eJADZ4JpDev5nZ2r2LTET6ud2VSSgz5ijf7JwXGlcFBxUAqW2goPw6OxnAU8nfTTcDTOvtP9GZ07zq1kB4sf3HeWY+eGuXbz+Docl0v4i7du5lTfSOLDfiZKXML/fcdWIDcuJbD+f0TUcihUVByUgqXernXIlWvJCQ4D7Dme2WyJ/pEIvlIXXncJa+oqcLsk4Z56y5amScdftr6eqzdarqZ0x3xeuq6Oz950PrdfuS6jvWULd4mLugovXdpCoyBRcVAKllq7SjpXQem2riFqyktZUVM2SRx+8spp7n/p1JTn9o9EqC6zxM3jdrGmvoJgKMp5ywMsrUrd2uUv37aZ5oYKLllXm/Ye37tzVU4nHmqtQ+FSdO0zlMWD03wvd5bDIOsaKllZU8ZTbecwxiAijEbj/PV/v0ZpiYt3XLgsZWVy33CEqrLSxM8tjZW0dg7yls2TrYbEMU1+Hvmzq+fjV5k3tEq6cFHLQSlYnLbdPTlq293WNcS6+gq2r66hayBMe6812Oaxg530DlvFZ1M11esfmSAOdsfU6cShEFHLoXBRy0EpWHLZXykYitA9GKa5sZLtqy03z57jvaysLeeHe9vxe90MhKM8ebibDSlaZfePRMall/7OzlU0+L1sXRZYsN9hIWi0xSEeNzp5scBQy0EpWPxeN6UlkhO3khOMXldfwcYlfiq9bnYf76F3aJRHX+/kty9eyZq6cp6yK5snEhyJUF0+ZjksqfLx/ktXL3hzvPmmwe8lGjeJdiFK4TAncRCRYyLyqoi8JCK77bVaEXlYRA7bX2uSjv+kiLSKyEERuT5pfbt9nVYR+ZIstneIMi9Y/ZVyUyXd1un0OLJqCN6wqprdx3r58SunicQMN21bweUt9Tx75Nyk5nlg9VZKdistVsaqpDVjqdDIhuXwJmPMRcaYHfbPnwAeMca0AI/YPyMiW4Cbga3ADcAdIuKUeX4ZuB1osf/dkIV9KUVArgrhjnQP4nYJq+ss19C2VTUc7BjgO88eZ/PSAFuWBbh8fT1DozFeOtk37txILM7waKwoxEGrpAuX+XAr3QjcbX9/N/DOpPXvGWPCxpijQCuwU0SWAgFjzDPGGAPck3SOokxLXWVuWmgc6RpiVW05pSXWW2jHmhqMgUMdg/zWtuUA7FpXj0vgycPjXUvJBXCLncak4UNKYTFXcTDAQyKyR0Rut9eajDFnAOyvTrnncuBk0rnt9tpy+/uJ65MQkdtFZLeI7O7q6prj1pXFQK4sh7auwXFtsy9aWY1LrMrkGy+yXr5V5aWcv6KaJ1tTi0NyzGGx0qAtNAqWuYrDG40x24C3Ah8RkSunOTZVHMFMsz550ZivGWN2GGN2NDQ0ZL5bZdGRC3GIxQ3Hzg2zLqlS2e8rZduqGq7b0pT4QAS4fH0dL53sYyA0FpB1ZkcHisByqPC6qfCUqFupAJmTOBhjTttfO4EfATuBDttVhP210z68HViZdPoK4LS9viLFuqLMSH2ll8FwdEH7K53qHWE0Gqd5wsCdb992Cf/82xeNW7t8fQOxuOHZIz2JtWARuZXAmSWtAelCY9biICIVIuJ3vgeuA14DHgButQ+7Fbjf/v4B4GYR8YrIWqzA8/O262lARC61s5RuSTpHUabFqZJeSOuhrXssUymZMk/JpFba21ZX4yt1jUtpLaaYA1gZS2o5FB5zKYJrAn5kZ526gf80xvxCRF4A7hOR24ATwHsAjDH7ROQ+YD8QBT5ijHFu9z4M3AWUAT+3/ynKjCSLw7Lq1D2Jsk0ijXWagTsOXncJO9fWjRv12We3664uEnFo8Hs5cCaY620oGTJrcTDGHAEuTLF+DrhminM+DXw6xfpu4LzZ7kUpXupy0F/pSPcQ1eWlCWGaiQuWV/Hl1m7C0Rhedwn9I1GgOGIOYInD44fUcig0tEJaKWjGLIf5+/AxxvDd50/wgz3ttPcOWw336ivSrmbeuMRPLG4S85z7RyJUeEoSabCLnQa/l4FwlJHR3MzdUGaH9lZSCpq6Sqdt9/xZDm1dg3zyh6+OW3v39hVTHD2ZzUut3koHzw6wdVnVpKZ7i53GpEK4VXULP650MRGLG778q1au2dzE5qXz24dLxUEpaAK++e+v9ND+DgC+fdtO2joHeelkH/9jx8oZzhpjTV0FHreL188OANA/MkpVeXouqcVAY2CshYaKw9xo7x3mHx86RKPfp+KgKNMhItSUe+a1bffD+zs4f3kVV7Q0cEVL5vU17hIXLY2VSeIQoaqseN56DZXaQiNbHO6wXJPrm+Z/9GtxOD2VRU1txfy10OgcCPHSyT6uTTG6MxM2LQnwup2xU2xupaVVluVwqm8kxzspfA51WjcYCzEXXMVBKXjqKj3zFpB+9EAnxpAFcfDTORCmZ2i06MShpsJDXYUncderzJ7WjkGWBHwEfPP/+lFxUAqeRr+P1s5BTvYMZ/3avzzQwfLqMjYtmTywJxM22UHp188G6RuOUF1EMQeAlqZKDnYM5HobBc/hzkFaFsClBCoOyiLgg1etA+B3v/EcncHstWkYHo3yxOFurt3SNOchPBttcXn5ZD/haLyoLAeADU1+WjsHsRovK7MhbqdDL4RLCVQclEXApiUB7vr9nXQNhHnf15/LWiuNJw53E47GuW6OLiWwgrJ1FR6eP2pVShdLAZxDS5OfwXCU0/3aY2m2nOobYSQSSzl2dj5QcVAWBdtW1fD1W3ZwvGeY2+/ZnZVr/nJ/BwGfm4vX1s75WiLCxiV+dh/rBYqnr5LDBvtu95C6lmaNU0TZopaDomTGZevr+bNrN7D7eO+s4w/RWJzDHQP8+OXT/PJAB2/a1Ji1SuZNSwIMhK3WGcXSV8nBuds9rOIwaw4vYKYSaJ2Dssi4emMjn/356zxz5BwrazMruOoIhrj2C78mGLI+wD0lLm7aln4l9EwkB7WLzXKoqfDQ4PdySDOWZs2hjkEa/N4FS2ZQcVAWFRuaKqmr8PBM27mMqpgBXjrZRzAU5eM3bOKqDQ00N1bgdZfMfGKaOBlLUHziANb/jVoOs+dw5+CCuZRA3UrKIkNE2NVcx9Nt3RlnxrR1WXe1v3vpKrYsC2RVGABaGv04SU/FKA4tjX4Odw4Sj2vGUqYYY2jtGFBxUJS5cFlzPR3BMEe6hzI6r61ziKaAF/88FRiVeUpYW2fNgCi2bCWw4g7DozGtlJ4FZ/pDDI3GWL9AmUqg4qAsQi5rrgPg6aQBO+nQ1jVIc8P83pltWurH73NT4ppb3UQhssEu3nICq4uFT/zgFf76v1+b1+c4bGcqbVDLQVFmz+q6cpZV+XimrXvmg22MMQsiDh+6qplP/ebWeX2OfKXFvutdTEHpSCzOAy+f5oGXT8+ru8yJ1bSo5aAos8eKO9TzTNu5tN+wXYNhBkJRmhtmHv05Fy5YUZ3RLIjFRFVZKU0B76KqdXjtVD/DozH6RyKJpnjzweGOQeoqPGlPH8wGKg7KouSy5jp6hyOJNtkz0dZpxSeaF9BsL0Y2NPkXlTg8d7Qn8f3zSd8DHO0e4ukMrNd43PDaqf6Ujx3uHFiw+gYHFQdlUbIrEXdI783pZCrNt1up2GlptHosLZaMpeeOnKO5oYJlVb5xQgHw1//9Grffs4dYmr/rPc8c4zf+9UleOtk3bt0Ys6AN9xxUHJRFybLqMtbWV/BMmkHptq5Byj0lLLGnlinzw4amSkKROCd7s99Bd6GJxQ27j/Vyybo6dq6t5fmjPYn06XODYZ5u62YwHOVgGtZrKBLjy79uA+Cx1zvHPdYRtFyeLY0LF28AFQdlEbOruY7njvYQicVnPLata4h1DRW4ijCLaCFZTEHpA2eCDISjXLK2lp1r6+gaCHPsnCV6D+7rwDEY9hzvmeYqFt/ffZKOYJjq8lIeP9w17rEn7J+3rarJ7i8wAyoOyqLlypYGBsNRXjg685uzrXP+M5WUsXTW/aeDOd7J3Hn2iGWVXrK2jp1rrQ9up+vuT189zdr6Chr9XvYc7532OuFojDt+1caO1TXcculqXj7ZR/9wJPH4g/vOsry6jPOWz+/M6ImoOCiLlis31ONxu3j4QMe0x43YhVkqDvOP31fKecsDibvhQua5oz2sritnSZWP5oZKais8PH+0l+7BMM+0nePt5y9l++oads8gDv+1p50z/SE+ek0LV25oIG7gKTtWNhiO8vjhbq7fumTOM0UyRcVBWbSUe9xcvr6eh/d3TNtK40i3BqMXkqs3NLL3RO+4u+NCIx43vHCsh51rrHbuIsLONbU8f+wcD+47S9zA2y+wxKG9d4SOKYZQjUbj3PFYG29YVc0VLfVctLIav9fN44cs8fzVwU5Go3FuOG/Jgv1uDioOyqLm2i1NtPeOTJvS2tblpLHOb42DYnH1Ruvu+InWwrUeDnUO0Dcc4ZJ1dYm1nWtrOdkzwl1PHWNdfQWblvjZvtpyN+2dwnr44d52TvWN8NFrWhAR3CUuLltfx+OHujDG8OC+DuorPYnrLCQqDsqi5ppNjYA1uGcqWjsHcQmsqVNxWAguWllNwOfmVwcLSxw6B0KJKYPPHbHiWJckDYLaaX9/uHOQt1+wFBFh67IqvG5XStdSJBbn3x5r5YIVVVy9oSGxfuWGBk73h9h/JsijBzq4dktTTtqtaMtuZVHTGPBx0cpqHj7QwR9d0wJYeeOdA2Ga7LTVtq5BVtaW4yvNbhdWJTXuEhdXbGjg1/bd8UL70mdD92CY6/75cfqGI2xs8hOOxlheXTZuZsjmpQH8XjcD4ShvO38pAB63iwtXVKcMSv9o7ynae0f4m3dsHfc3uLLFEorP/ux1hkZjXL914V1KoJaDUgRcu6WJV9r7OdsfwhjD3/x4P5d85hH+/bFWq6eSZiotOFdvaKBrIMz+M4WRtfTZn73OUDjKH715PY0BLx3BMNdtHT9bvMQlvHF9PRub/OMGO21fU8O+0/2EIrHEmmM1nL+8ijfb1q3Dytpy1tVX8GRrN36vm8ua6+f3l5sCtRyURc+1W5r4/IMH+eWBDroHw9z1tOUT/vyDB+kbHuVo9xBXtOTmDVisXGW7UX51sIuty6pyvJvpeeFYDz/Y286Hr27mz67bCFgB6VQ1MZ9/zwVEY+Otoe2ravhyzPBKe3/C9fTfL57iRM8wd96yI6XldEVLPUe6h7hmcyMed27u4dVyUBY9LY2VrK4r54u/PMQXf3mYd29fwcMfu4pbdq3mzieOEo7G1XJYYBoDPrYsDfDrPIg7jEbjfPf5EykziiKxOH/1o9dYXl3GH715fWJ9qmJJv6+UmgnN8Zxg8m67GC5qWw1blwV4y+bGSdcAuNq2Jt5qu6dygVoOyqJHRHjL5ia+8eRRrt/axN/fdD4lLuFv3rGV6rJSvvRoK+ctz++718XI1Rsb+OrjRwiGIgQmDFgaCEUoLXHNexwoGIrwwXv28MyRc9RXevm333kDlyZlIN399DEOdgzw1fdvp9wzu4/LmgoPzQ0V/PpgF35fKQ/tO8vxc8N89f3bp4y3XL2hgR98eNeCV0Uno+KgFAUfvGodDX4vv3fZGtwllsEsInzsuo3cdvk6qsqLbzJbrrl6YyN3/KqNpw53j7tDDkVi/Oa/PklVuYcffviycZk6n7r/NV441stnbjqfi1ZWz+n5z/SP8IFvvUBr5yAfv2ET3999kvd9/Tn+9/UbqfCU8OOXz/DC8R7etLGB67Y0zXzBadixupZ7d5/kuaM9LK8u48NXN097TRFh++raKR9fCCTTObv5wo4dO8zu3btzvQ1FUWZJNBZn52ceYU1dOd//0JgI/PtjrXz+wYMA/O2NW7ll1xoAHj/UxS3ffB6v20UkFufDVzfz0WtaZjXr+0z/CDfd8TQDoShf+d3tXN5Sz0Aowp9//2Ue3GelPbc0VvKbFy7j1l1r5nzz0N47zJOHu9m5tpa19RU5zdASkT3GmB0zHqfioChKrrj/pVP88fde4hNv3cSHrmqmMxji6n/8FW9cX8/IaIyXT/bxyJ9fRYXHzfVffByP28W9t+/i8w++zn2729m6LMA9v7+Tukpv4pone4Z59PVO3vmG5VSlmNUdixve9/VnebW9n/s+tGtcQNwYw68PdbG0qoyNSxa2C+pCka44aEBaUZSc8Y4Ll/HW85bwhYcOcfDsAJ9/8CCRWJy/fNtm/vbGrYSjcT7z0wN84eFDtPeO8Pc3XUCD38s/vPtC7rxlB21dg7z3zmfpGggDsPtYDzf++1N86oF9XPG5R/n3x1oZCkfHPeedTxzh2SM9fOodWydlSokIV29sXLTCkAlqOSiKklPO2QVmlT43J3qGuf2KdXzybZsB+MJDB/nSo62IwPsuWcXfvfP8cec+3dbNbXftZlm1j1svW8Pf/eQAy2vK+Iu3bebeF07wywOd1FZ4uO3ytdyyazXHzw3zrjue4i2bm7jjfdsKogAv2xScW0lEbgD+BSgBvm6M+fvpjldxUJTFwy9eO8uHvrOHugoPj/3vqxPZS6FIjOu/+DjhSJyHPnblpKwmsKaxfeCuFxgejbFzbS1f/d3tiXTSF0/08i+PHOZXB7sI+NxUeN0YA7/4kyuoLl+4ecz5REGJg4iUAIeAa4F24AXgvcaY/VOdo+KgKIuLu58+xoYmf2LEq0PP0CixuKHB753iTHjpZB9PtXbzP69Yl7Jo7NX2fv7tscM8drCLb/3exbxxffEWPRaaOOwC/q8x5nr7508CGGM+O9U5Kg6KomRKLG5y0sQunyi0gPRy4GTSz+322jhE5HYR2S0iu7u6cl9ZqShKYVHswpAJ+SIOqf7HJpk0xpivGWN2GGN2NDQ0pDhFURRFyQb5Ig7twMqkn1cAp3O0F0VRlKInX8ThBaBFRNaKiAe4GXggx3tSFEUpWvKit5IxJioifwg8iJXK+k1jzL4cb0tRFKVoyQtxADDG/Az4Wa73oSiKouSPW0lRFEXJI1QcFEVRlEnkRRHcbBCRAeDgDIdVAf0ZXDbT4yeyCjgxj8+n+5vb8+n+5vZ8+b6/uZ67mPdXD3Tb3682xsxcC2CMKch/wO40jvlahtfM6PgU53fN5/Pp/nR/ur/5O3cx7y+dz8uJ/xa7W+nH83z8RPrm+fl0f3N7Pt3f3J4v3/c313N1f0kUsltpt0mjP8hCko97Skb3Nzd0f3ND9zc35rK/2ZxbyJbD13K9gRTk456S0f3NDd3f3ND9zY257C/jcwvWclAURVHmj0K2HBRFUZR5QsVBURRFmYSKwywQkcFc72EqRORdImJEZFOu9zIdM/0NReRXIrLgwUERWSEi94vIYRFpE5F/sZtBTnX8n4hI+QLuL29fe6Cvv7mST68/FYfFx3uBJ7E626aNPaq1qBFr2vwPgf82xrQAG4BK4NPTnPYnwIKJQwGgr79Zkm+vPxWHWSIilSLyiIjsFZFXReRGe32NiBwQkTtFZJ+IPCQiZQu1J+CNwG3Yb04RuVpEHheRH4nIfhH5ioi47McGRTRA6OQAAAa/SURBVORvReQ5YNdC7HHCfq8WkZ8k/fxvIvJ7C72PJN4MhIwx3wIwxsSAPwV+X0QqROQf7f/rV0Tkj0Tko8Ay4DEReWyhNpmPrz1nX+jrby7k1etPxWH2hIB3GWO2AW8C/slWfoAW4N+NMVuxCld+a4H29E7gF8aYQ0CPiGyz13cCfwacDzQDN9nrFcBrxphLjDFPLtAe85mtwJ7kBWNMEKtlwR8Aa4E3GGMuAP7DGPMlrKFUbzLGvGkB95mPrz3Q199cyavXn4rD7BHgMyLyCvBLrJnXTfZjR40xL9nf7wHWLNCe3gt8z/7+e/bPAM8bY47YdyLfBS6312PADxZob4WAkGI8rb1+JfAVY0wUwBjTs5AbS7GffHvtgb7+5kpevf7yZp5DAfI+oAHYboyJiMgxwGc/Fk46LgbMu2kvInVYZul5ImKwhiYZrBkZE19wzs8h+w2bK6KMv0HxTXXgArGPCXfaIhLAGmF7hNRv3FyQV6890Ndflsir159aDrOnCui035xvAlbneD/vBu4xxqw2xqwxxqwEjmLdpe0UawSrC/htrIBhPnAc2CIiXhGpAq7J8X4eAcpF5BZIBEn/CbgLeAj4kIi47cdq7XMGAP8C7zPfXnugr79skFevPxWHDLH/c8LAfwA7RGQ31p3c6zndmGXC/2jC2g+A3wGeAf4eeA3rDTvxuAXF+RsaY04C9wGvYP09X8zlvozVLuBdwHtE5DBwCMu//xfA17F8v6+IyMtYf1ew2hL8fCEC0nn82gN9/c2ZfHv9afuMDBGRC4E7jTE7c72XdBCRq4E/N8b8Rq734lBof8N8oRD/bvr6K1zUcsgAEfkQVkDtr3K9l0JF/4azQ/9u2UH/jumjloOiKIoyCbUcFEVRlEmoOEyDiKwUkcfsqtN9IvLH9nqtiDwsVv+Th0Wkxl6vs48fFJF/m3Ct37YrG/eJyD/k4vdRCotZvP6uFZE9dhXtHhF5c9K1ttvrrSLypaSiOUVJiYrD9ESBPzPGbAYuBT4iIluATwCP2P1PHrF/Biuz4K+BP0++iJ0D/nngGrtytUlEcp02p+Q/mb7+uoHfNMacD9wKfDvpWl8GbseqoG4BbliYX0EpVFQcpsEYc8YYs9f+fgA4gFWNeiNwt33Y3VhtAzDGDNltAEITLrUOOGSM6bJ//iUL29ZAKUBm8fp70Rhz2l7fB/jsHP6lQMAY84ydLnmPc46iTIWKQ5qIyBrgDcBzQJMx5gxYb2CgcYbTW4FNYjVGc2O9MVfO326VxcYsXn+/BbxojAljCUp70mPt9pqiTIm2z0gDu9vkD4A/McYEM3XXGmN6ReTDwL1AHHgay5pQlBnJ9PUnIluBzwHXOUspDtM0RWVa1HKYAREpxXpj/ocx5of2codtqmN/7ZzpOsaYH9vdJ3cBB4HD87VnZfGQ6etPRFZgVSDfYoxps5fbgRVJl12B1c1TUaZExWEa7IyObwAHjDFfSHroAayAH/bX+9O4VqP9tQb4X1jl8IoyJZm+/kSkGvgp8EljzFPOwbbraUBELrWveQtpvGaV4kaL4KZBRC4HngBexXIHgdXn5DmsniyrsPqdvMdpoWt3yAwAHqx++tcZY/aLyHeBC+1r/K0xxmltrCgpyfT1JyJ/BXyS8VbpdcaYTrFGXt6F1aX158AfGX3zK9Og4qAoiqJMQt1KiqIoyiRUHBRFUZRJqDgoiqIok1BxUBRFUSah4qAoiqJMQsVBUeYBEfmQMws4zePXiMhr87knRckEbZ+hKFlGRNzGmK/keh+KMhdUHBQlBXaju19gFZy9AWvY+y3AZuALQCVWi+zfM8acEZFfYfXMeiPwgIj4gUFjzD+KyEXAV4ByoA34fbvf1nbgm8Aw8OTC/XaKMjPqVlKUqdkIfM0YcwEQBD4C/CvwbmOM88H+6aTjq40xVxlj/mnCde4BPm5f51XgU/b6t4CP2v22FCWvUMtBUabmZFKPou9gta44D3jY7oxaApxJOv7eiRcQkSos0fi1vXQ38P0U698G3pr9X0FRZoeKg6JMzcTeMgPAvmnu9IcyuLakuL6i5A3qVlKUqVklIo4QvBd4Fmhw1kSk1J6dMCXGmH6gV0SusJfeD/zaGNMH9NvN9QDel/3tK8rsUctBUabmAHCriHwVq9PpvwIPAl+y3UJu4ItYIzmn41bgKyJSDhwBPmCvfwD4pogM29dVlLxBu7IqSgrsbKWfGGPOy/FWFCUnqFtJURRFmYRaDoqiKMok1HJQFEVRJqHioCiKokxCxUFRFEWZhIqDoiiKMgkVB0VRFGUS/x+G+ecoXNZq/QAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "sorted_data['inc'][-100:].plot()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "On constate que la periode creuse se situe en fin d'été"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Analyse statistique"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Nous utilisons des periodes annuelles allant du 1er Septembre au 1er Septembre."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 19,
+ "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": [
+ "Puis nous calculons les incidences annuelles en utilisant la même technique que pour les états gripaux"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 20,
+ "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",
+ " 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": [
+ "On a alors les incidences annuelles"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 21,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 21,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAD8CAYAAACyyUlaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHgNJREFUeJzt3X+QVeWd5/H3B5uAOmLAgBEQYSLjBswOBqox4242LhPAZEuwRjM9OkplqMIxmElSW5Vo6RaW8sc4lawbytKVxETUiLJsKNmNRFvc1Li7TANGE0HC0BkJIkh3qolgtujY8N0/znPldNt03+57u++P/ryqTt3T33uew3k8Vn/7+XGeo4jAzMysFKMqfQFmZlb7nEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWckaKn0B5faxj30spk+fXunLMDOrKa+88spvI2LiYMvXXTKZPn06O3furPRlmJnVFEm/KaV8Ud1ckr4habekXZLWSxor6R5Jb0t6LW1fyB1/p6RWSXslLcrF50p6PX23RpJSfIykZ1K8RdL0XJllkvalbVkplTUzs6HRbzKRNAX4O2BeRFwOnAU0pa8fiIg5aXsuHT8rfT8bWAw8JOmsdPzDwApgZtoWp/hy4GhEXAo8ANyfzjUBWAXMBxqBVZLGl1ZlMzMrt2IH4BuAsyU1AOcAh/o4dgnwdER0RsSbQCvQKOkiYFxEbItsqeLHgaW5MuvS/kZgQWq1LAKaI6IjIo4CzZxOQGZmViX6TSYR8TbwbeAAcBh4NyJeSF/fLumXkn6QazFMAd7KneJgik1J+z3j3cpERBfwLnBBH+cyM7MqUkw313iylsMMYDJwrqS/Juuy+gQwhyzJfKdQpJfTRB/xwZbJX+MKSTsl7Wxvb++jNmZmNhSK6eb6c+DNiGiPiPeBHwN/FhFHIuJkRJwCvkc2pgFZ6+HiXPmpZN1iB9N+z3i3Mqkr7Xygo49zdRMRayNiXkTMmzhx0DPbzEastmMn+NIj22g7fqLSl2I1qphkcgC4UtI5aRxjAbAnjYEUXAfsSvubgaY0Q2sG2UD79og4DByXdGU6zy3As7kyhZla1wMvpXGV54GFksanFtLCFDOzMlqzdR879new5sV9lb4Uq1H9PmcSES2SNgI/B7qAV4G1wPclzSHrdtoP3JqO3y1pA/BGOn5lRJxMp7sNeAw4G9iSNoBHgScktZK1SJrSuTok3QfsSMfdGxEdpVTYzE677O4tdHad+uDnJ1sO8GTLAcY0jGLv6msqeGVWa1Rv74CfN29e+KFFs+K0HTvB6uf28MLudzjx/inGjh7Fotkf564vfpJJ542t9OXZMJL0SkTMG2x5r81lNoJNGjeW88Y00Nl1ijENo+jsOsV5YxqcSGzA6m45FTMbmN++18lN8y/hxsZpPLX9AO0ehLdBcDdXTtuxE9y+/lUevPEK/2VmZiOKu7nKyDNazMwGx91ceEaL2UC5FW89uWUCvPzNq7l2zmTGjs7+c4wdPYolcybz8reurvCVmVUnt+KtJ7dM8IwWs2K5FW9n4pZJUpjRsukrV3HT/Etof6+z0pdkVnXcirczccskeeTm05MYVi+9vIJXYla93Iq3M3EyMbMB8XMp1hs/Z2JmZn7OxMzMKs/JxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnEzMxKVlQykfQNSbsl7ZK0XtJYSRMkNUvalz7H546/U1KrpL2SFuXicyW9nr5bI0kpPkbSMyneIml6rsyy9G/sk7SsfFU3M7Ny6TeZSJoC/B0wLyIuB84CmoA7gK0RMRPYmn5G0qz0/WxgMfCQpLPS6R4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictMzOrDsV2czUAZ0tqAM4BDgFLgHXp+3XA0rS/BHg6Ijoj4k2gFWiUdBEwLiK2Rba65OM9yhTOtRFYkFoti4DmiOiIiKNAM6cTkJmZVYl+k0lEvA18GzgAHAbejYgXgAsj4nA65jAwKRWZAryVO8XBFJuS9nvGu5WJiC7gXeCCPs5lZmZVpJhurvFkLYcZwGTgXEl/3VeRXmLRR3ywZfLXuELSTkk729vb+7g0MzMbCsV0c/058GZEtEfE+8CPgT8DjqSuK9JnWzr+IHBxrvxUsm6xg2m/Z7xbmdSVdj7Q0ce5uomItRExLyLmTZw4sYgqmZlZORWTTA4AV0o6J41jLAD2AJuBwuyqZcCzaX8z0JRmaM0gG2jfnrrCjku6Mp3nlh5lCue6Hngpjas8DyyUND61kBammJmZVZF+X9sbES2SNgI/B7qAV4G1wB8BGyQtJ0s4N6Tjd0vaALyRjl8ZESfT6W4DHgPOBrakDeBR4AlJrWQtkqZ0rg5J9wE70nH3RkRHSTU2M7Oy82t7zczMr+01M7PKczIxM7OSOZmYmVnJnEzMzErQduwEX3pkG23HT1T6UirKycTMrARrtu5jx/4O1ry4r9KXUlH9Tg02M7MPu+zuLXR2nfrg5ydbDvBkywHGNIxi7+prKnhlleGWiVkR3JVhPb38zau5ds5kxo7Ofo2OHT2KJXMm8/K3rq7wlVWGk4lZEdyVYT1NGjeW88Y00Nl1ijENo+jsOsV5YxqYdN7YSl9aRbiby6wP7sqwvvz2vU5umn8JNzZO46ntB2gfwS1XPwFv1oe2YydY/dweXtj9DifeP8XY0aNYNPvj3PXFT47Yv0CtPvkJeLMh5K4Ms+K4m8usH+7KMOufu7nMzMzdXGZmVnlOJmZmVjInEzMzK5mTiZmZlczJxEYsL5FiVj5OJlYWtfiL2UukmJWPnzOxssj/Yl593acqfTl98hIpZuXn50ysJD1/MRdU8y9mL5Fi9mF+zsQqqhaX4fYSKWbl128ykXSZpNdy2zFJX5d0j6S3c/Ev5MrcKalV0l5Ji3LxuZJeT9+tkaQUHyPpmRRvkTQ9V2aZpH1pW1be6lupavUXc2GJlE1fuYqb5l9C+3udlb4ks5rW75hJROwF5gBIOgt4G9gEfBl4ICK+nT9e0iygCZgNTAZelPQnEXESeBhYAfwT8BywGNgCLAeORsSlkpqA+4G/lDQBWAXMAwJ4RdLmiDhacs2tbGpx7apHbj7dml+99PIKXsnQaTt2gtvXv8qDN15R9cndat9AB+AXAL+OiN+kRkVvlgBPR0Qn8KakVqBR0n5gXERsA5D0OLCULJksAe5J5TcCD6ZWyyKgOSI6UplmsgS0foDXbUNoJPxirkW1NCnCat9Ak0kT3X+R3y7pFmAn8B9Ti2EKWcuj4GCKvZ/2e8ZJn28BRESXpHeBC/LxXsp8QNIKshYP06ZNG2CVzOqLZ6tZJRQ9AC/pI8C1wH9LoYeBT5B1gR0GvlM4tJfi0Ud8sGVOByLWRsS8iJg3ceLEM9bBbCSoxUkRVvsGMpvrGuDnEXEEICKORMTJiDgFfA9oTMcdBC7OlZsKHErxqb3Eu5WR1ACcD3T0cS4zO4NanRRhtW0gyeSvyHVxSboo9911wK60vxloSjO0ZgAzge0RcRg4LunKNB5yC/Bsrkxhptb1wEuRPQDzPLBQ0nhJ44GFKWZmffBsNRtuRY2ZSDoH+Dxway78D5LmkHU77S98FxG7JW0A3gC6gJVpJhfAbcBjwNlkA+9bUvxR4Ik0WN9BNjZDRHRIug/YkY67tzAYb2Zn5kkRNtz8BLyZmfkJeDMzqzwnEzMzK5mTiZmZlczJxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnErIa0HTvBlx7ZRlsNvBrZRhYnE7Makn8Vr1k1Gehre82sAvwqXqt2bpmY1QC/iteqnZOJWQ3wq3it2rmby6xGFF7Fe2PjNJ7afoB2D8JbFfGbFs3MzG9aNDOzynMyMTOzkjmZmJlZyfpNJpIuk/Rabjsm6euSJkhqlrQvfY7PlblTUqukvZIW5eJzJb2evlsjSSk+RtIzKd4iaXquzLL0b+yTtKy81Tczs3LoN5lExN6ImBMRc4C5wP8DNgF3AFsjYiawNf2MpFlAEzAbWAw8JOmsdLqHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqsNAu7kWAL+OiN8AS4B1Kb4OWJr2lwBPR0RnRLwJtAKNki4CxkXEtsimkD3eo0zhXBuBBanVsghojoiOiDgKNHM6AZnZCOZ1yqrLQJNJE7A+7V8YEYcB0uekFJ8CvJUrczDFpqT9nvFuZSKiC3gXuKCPc5nZCOd1yqpL0Q8tSvoIcC1wZ3+H9hKLPuKDLZO/thVk3WdMmzatn8szs1rmdcqq00BaJtcAP4+II+nnI6nrivTZluIHgYtz5aYCh1J8ai/xbmUkNQDnAx19nKubiFgbEfMiYt7EiRMHUCUzqzVep6w6DSSZ/BWnu7gANgOF2VXLgGdz8aY0Q2sG2UD79tQVdlzSlWk85JYeZQrnuh54KY2rPA8slDQ+DbwvTDEzG6G8Tll1KqqbS9I5wOeBW3Phvwc2SFoOHABuAIiI3ZI2AG8AXcDKiDiZytwGPAacDWxJG8CjwBOSWslaJE3pXB2S7gN2pOPujYiOQdTTzOqI1ymrPl6by8zMvDaXmVkxPJV4aDmZmNmI4KnEQ8vvMxlB2o6d4Pb1r/LgjVd4sNJGDE8lHh5umYwg/svMRiJPJR4ebpmMAP7LzEYyTyUeHm6ZjACD+cvMg5VWTwpTiTd95Spumn8J7e91VvqS6o5bJiPAYP4yy3eJrb7uU8N4tWbl98jNp2e8rl56eQWvpH45mYwQxT7k5S4xMxsMP7Ro3bQdO8Hq5/bwwu53OPH+KcaOHsWi2R/nri9+0n3MZnXMDy1aWXmw0swGw91c9iFe98jMBsrdXGZm5m4uMzOrPCcTMzMrmZOJmZmVzMlkGPhpcjOrd04mw8ALLJpZvfPU4CHkp8nNbKRwy2QIeelrMxspnEyGkJ8mN7ORwslkiA106WsP1ptZLSoqmUj6qKSNkn4laY+kz0i6R9Lbkl5L2xdyx98pqVXSXkmLcvG5kl5P362RpBQfI+mZFG+RND1XZpmkfWlbVr6qD49Hbp7H6qWXM2vyOFYvvbzbUti98WC9mdWiopZTkbQOeDkivi/pI8A5wNeB9yLi2z2OnQWsBxqBycCLwJ9ExElJ24GvAf8EPAesiYgtkr4C/OuI+FtJTcB1EfGXkiYAO4F5QACvAHMj4uiZrrVWl1PpOVhf4MF6MxsOQ76ciqRxwGeBRwEi4g8R8bs+iiwBno6Izoh4E2gFGiVdBIyLiG2RZbDHgaW5MuvS/kZgQWq1LAKaI6IjJZBmYPGAa1kDPFhvZrWsmG6uPwbagR9KelXS9yWdm767XdIvJf1A0vgUmwK8lSt/MMWmpP2e8W5lIqILeBe4oI9z1R0P1ptZLSsmmTQAnwYejogrgN8DdwAPA58A5gCHge+k49XLOaKP+GDLfEDSCkk7Je1sb2/voyrVze+pNrNaVcxDiweBgxHRkn7eCNwREUcKB0j6HvA/c8dfnCs/FTiU4lN7iefLHJTUAJwPdKT453qU+VnPC4yItcBayMZMiqhTydqOneD29a/y4I1XlK314PdUm1mt6rdlEhHvAG9JuiyFFgBvpDGQguuAXWl/M9CUZmjNAGYC2yPiMHBc0pVpPOQW4NlcmcJMreuBl9K4yvPAQknjUzfawhSrOM+6MiuOp7uPDMUup/JV4EdpJte/AF8G1kiaQ9bttB+4FSAidkvaALwBdAErI+JkOs9twGPA2cCWtEE2uP+EpFayFklTOleHpPuAHem4eyOiY3BVLQ8vkWI2MPk/vFZf96lKX44NEb9pcYDajp1g9XN7eGH3O5x4/xRjR49i0eyPc9cXP+nBcrMcT3evLX7T4jDzrCuz4ni6+8jiVYMHoTDr6sbGaTy1/QDt7gs2+xD/4TWyOJkMgmddWbkMxazAauI/vEYOj5mYVdDdm17nR9sPcFPjNA9OW0WVOmbilolZBXhWoNUbD8CblVkxz1V4cNrqjZOJWZkV80CrB6et3riby6xMBtp15cHpD6v3CQn1zAPwZmXiB1pL5wkJleMBeLMq4a6rwfOEhNrnMROzMvJrBAbHExJqn1smVtVqrQ/dD7QOjlt1tc8tExt2A1mS3Ev9jxxu1dU2D8DbsCtmkNUrzpoNr1IH4J1MbNgMJEF4ZpTZ8PIS9FYzBjLI6j50s9riAXgbNgNNEH6oz6x2OJnYsBpIgvDMKLPa4TETMzPzmImZmVWek4mZWQ8DeRbKMk4mZmY9+GHZgSsqmUj6qKSNkn4laY+kz0iaIKlZ0r70OT53/J2SWiXtlbQoF58r6fX03RpJSvExkp5J8RZJ03NllqV/Y5+kZeWruplZd5fdvYXpd/yEJ1sOEJEtODn9jp9w2d1bKn1pVa/Ylsl3gZ9GxL8C/hTYA9wBbI2ImcDW9DOSZgFNwGxgMfCQpLPSeR4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictM7Ny8oKTg9dvMpE0Dvgs8ChARPwhIn4HLAHWpcPWAUvT/hLg6YjojIg3gVagUdJFwLiI2BbZFLLHe5QpnGsjsCC1WhYBzRHRERFHgWZOJyAzs7Lyw7KDV0zL5I+BduCHkl6V9H1J5wIXRsRhgPQ5KR0/BXgrV/5gik1J+z3j3cpERBfwLnBBH+fqRtIKSTsl7Wxvby+iSmZmvfOCk4NTzEOLDcCnga9GRIuk75K6tM5AvcSij/hgy5wORKwF1kL2nEkf12Zm1ic/LDs4xbRMDgIHI6Il/byRLLkcSV1XpM+23PEX58pPBQ6l+NRe4t3KSGoAzgc6+jiXmZlVkX6TSUS8A7wl6bIUWgC8AWwGCrOrlgHPpv3NQFOaoTWDbKB9e+oKOy7pyjQeckuPMoVzXQ+8lMZVngcWShqfBt4XppiZmVWRYtfm+irwI0kfAf4F+DJZItogaTlwALgBICJ2S9pAlnC6gJURcTKd5zbgMeBsYEvaIBvcf0JSK1mLpCmdq0PSfcCOdNy9EdExyLqamdkQ8dpcZmbmtbnMzKzynEzMzOpApdcTczIxM6sDlV5PzC/HMjOrYZfdvYXOrlMf/PxkywGebDnAmIZR7F19zbBdh1smZmY1rFrWE3MyMTOrYdWynpi7uczMalxhPbEbG6fx1PYDtFdgEN7PmZiZmZ8zMTOzynMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMieTGlfp9XjMzMDJpOZVej0eMzPwQ4s1q1rW4zEzA7dMala1rMdjZgZOJjWrWtbjMTMDd3PVtGpYj8fMDLw2l5mZ4bW5zMysChSVTCTtl/S6pNck7UyxeyS9nWKvSfpC7vg7JbVK2itpUS4+N52nVdIaSUrxMZKeSfEWSdNzZZZJ2pe2ZeWquJnZcKvn58IG0jK5OiLm9GgGPZBicyLiOQBJs4AmYDawGHhI0lnp+IeBFcDMtC1O8eXA0Yi4FHgAuD+dawKwCpgPNAKrJI0fRD3NzCqunp8LG4oB+CXA0xHRCbwpqRVolLQfGBcR2wAkPQ4sBbakMvek8huBB1OrZRHQHBEdqUwzWQJaPwTXbWY2JEbCc2HFtkwCeEHSK5JW5OK3S/qlpB/kWgxTgLdyxxxMsSlpv2e8W5mI6ALeBS7o41xmZjVjJDwXVmwyuSoiPg1cA6yU9FmyLqtPAHOAw8B30rHqpXz0ER9smQ9IWiFpp6Sd7e3tfVbEzGy4jYTnwopKJhFxKH22AZuAxog4EhEnI+IU8D2yMQ3IWg8X54pPBQ6l+NRe4t3KSGoAzgc6+jhXz+tbGxHzImLexIkTi6mSmdmwKjwXtukrV3HT/Etof6+z3zK1NGDfbzKRdK6k8wr7wEJgl6SLcoddB+xK+5uBpjRDawbZQPv2iDgMHJd0ZRoPuQV4NlemMFPreuClyB6AeR5YKGl86kZbmGJmZjXlkZvnsXrp5cyaPI7VSy/nkZv7f6SjlgbsixmAvxDYlGbxNgBPRcRPJT0haQ5Zt9N+4FaAiNgtaQPwBtAFrIyIk+lctwGPAWeTDbxvSfFHgSfSYH0H2WwwIqJD0n3AjnTcvYXBeDOzelWLA/Z+At7MrMq0HTvB6uf28MLudzjx/inGjh7Fotkf564vfnLIxln8BLyZWZ2pxQF7L/RoZlaFam0hV3dzmZmZu7nMzKzynEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWcmcTKxu1NK7H8zqjZOJ1Y1aeveDWb3xQo9W82rx3Q9m9cYtE6t5L3/zaq6dM5mxo7P/nceOHsWSOZN5+VtXV/jKzEYOJxOrebX47gezeuNuLqsLtfbuB7N64/eZmJmZ32diZmaV52RiZmYlKyqZSNov6XVJr0namWITJDVL2pc+x+eOv1NSq6S9khbl4nPTeVolrZGkFB8j6ZkUb5E0PVdmWfo39klaVq6Km5lZ+QykZXJ1RMzJ9andAWyNiJnA1vQzkmYBTcBsYDHwkKSzUpmHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqkMp3VxLgHVpfx2wNBd/OiI6I+JNoBVolHQRMC4itkU26v94jzKFc20EFqRWyyKgOSI6IuIo0MzpBGRmZlWi2GQSwAuSXpG0IsUujIjDAOlzUopPAd7KlT2YYlPSfs94tzIR0QW8C1zQx7nMzKyKFPucyVURcUjSJKBZ0q/6OFa9xKKP+GDLnP4HswRXSHLvSdrbx/XVgo8Bv630RQyxeq9jvdcP6r+OI61+l5RysqKSSUQcSp9tkjaRjV8ckXRRRBxOXVht6fCDwMW54lOBQyk+tZd4vsxBSQ3A+UBHin+uR5mf9XJ9a4G1xdSlFkjaWcp871pQ73Ws9/pB/dfR9RuYfru5JJ0r6bzCPrAQ2AVsBgqzq5YBz6b9zUBTmqE1g2ygfXvqCjsu6co0HnJLjzKFc10PvJTGVZ4HFkoanwbeF6aYmZlVkWJaJhcCm9Is3gbgqYj4qaQdwAZJy4EDwA0AEbFb0gbgDaALWBkRJ9O5bgMeA84GtqQN4FHgCUmtZC2SpnSuDkn3ATvScfdGREcJ9TUzsyFQd8up1ANJK1LXXd2q9zrWe/2g/uvo+g3wfE4mZmZWKi+nYmZmJXMyGSaSfiCpTdKuXOxPJW1LS8z8D0njUvwjkn6Y4r+Q9LlcmZ+lZWpeS9ukXv65YSfpYkn/S9IeSbslfS3Fy7bsTiWVuX51cQ8lXZCOf0/Sgz3OVfP3sJ/6Vd09HET9Pq/s2cHX0+e/z51r4PcvIrwNwwZ8Fvg0sCsX2wH8u7T/N8B9aX8l8MO0Pwl4BRiVfv4ZMK/S9emlfhcBn0775wH/DMwC/gG4I8XvAO5P+7OAXwBjgBnAr4Gz0nfbgc+QPWe0BbimzupXL/fwXODfAH8LPNjjXPVwD/uqX9Xdw0HU7wpgctq/HHi7lPvnlskwiYh/JJuplncZ8I9pvxn4i7Q/i2y9MyKiDfgdUNXz3SPicET8PO0fB/aQrVZQzmV3KqZc9Rveqx6YgdYxIn4fEf8b6PYmsnq5h2eqX7UaRP1ejfQMIbAbGKvskY5B3T8nk8raBVyb9m/g9MOevwCWSGpQ9qzOXLo/CPrD1LT+T9XQfdCTslWfrwBaKO+yO1WhxPoV1MM9PJN6uYf9qdp7OIj6/QXwakR0Msj752RSWX8DrJT0Clmz9A8p/gOyG7gT+C/A/yV7Zgfgpoj4FPBv03bzsF5xPyT9EfDfga9HxLG+Du0lVvQSOpVShvpB/dzDM56il1gt3sO+VO09HGj9JM0mW6n91kKol8P6vX9OJhUUEb+KiIURMRdYT9avTkR0RcQ3IlvyfwnwUWBf+u7t9HkceIoq6jqRNJrsf+IfRcSPU/hIajYXuj9KWXanospUv3q6h2dSL/fwjKr1Hg60fpKmApuAWyLi1yk8qPvnZFJBhRkgkkYBdwP/Nf18jrKla5D0eaArIt5I3V4fS/HRwH8g6yqruNTMfxTYExH/OfdVOZfdqZhy1a/O7mGv6ugenuk8VXkPB1o/SR8FfgLcGRH/p3DwoO9fpWYejLSNrOVxGHifLPMvB75GNuPin4G/5/RDpNOBvWQDaC8Cl6T4uWQzu35JNmD2XdIMoUpvZLNeIl3ba2n7AtmrBLaStay2AhNyZe4ia43tJTdbhGyywa703YOF/y71UL86vIf7ySaWvJf+v55VZ/fwQ/Wr1ns40PqR/QH7+9yxrwGTBnv//AS8mZmVzN1cZmZWMicTMzMrmZOJmZmVzMnEzMxK5mRiZmYlczIxM7OSOZmYmVnJnEzMzKxk/x+qftJYVEr3+gAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "yearly_incidence.plot(style='*')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Et un histogramme:"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 22,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 22,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFGtJREFUeJzt3XuQJWV9h/Hnxy4oMLggDKgLukbIRGFVZL2gVTqDxhKXaGm8IZpotNYqFTFlyqCoaLytF7zESyobJZB4nShUkDUaDI6oUZQV4oI4amBVlighKjCI4Oovf3QvGadmdmZO95lzOu/zqdraPrfu73mn53v69Ok+E5mJJOn/v70GHUCStDIsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhVq/kwg488MA88sgjV3KRrbn11lvZf//9Bx1j2bqaG8w+KF3N3tXcsHj2bdu23ZiZo02Xs6KFf9hhh3HZZZet5CJbMzU1xfj4+KBjLFtXc4PZB6Wr2buaGxbPHhE/bGM57tKRpEJY+JJUCAtfkgph4UtSISx8SSpEo8KPiLGIuGLWv5sj4uVthZMktafRYZmZOQ08GCAiVgE7gfNbyCVJalmbu3QeC/xnZrZyvKgkqV3R1t+0jYizgW9l5vvnXL8J2AQwOjp63OTkZCvLW2kzMzOMjIwMOsay9Zp7+86b+pBmcevXrrlzuqtjDmYfhK7mhsWzT0xMbMvMDU2X00rhR8Q+wPXA0Zn504XuNzY2ltPT042XNwhdPYuv19zrTt/afpgl2LF5453TXR1zMPsgdDU3LOlM21YKv61dOidSbd0vWPaSpMFqq/BPBj7e0rwkSX3QuPAjYj/gD4HzmseRJPVL42/LzMxfAge3kEWS1EeeaStJhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUiMaFHxEHRsSnIuK7EXF1RBzfRjBJUrtWtzCP9wKfy8ynRcQ+wH4tzFOS1LJGhR8RdwMeDTwPIDPvAO5oHkuS1LbIzN4fHPFgYAvwHeBBwDbgtMy8ddZ9NgGbAEZHR4+bnJxsFHhQZmZmGBkZGXSMZes19/adN/UhzeLWr11z53RXxxzMPghdzQ2LZ5+YmNiWmRuaLqdp4W8Avg48KjMvjYj3Ajdn5mvnu//Y2FhOT0/3vLxBmpqaYnx8fNAxlq3X3OtO39p+mCXYsXnjndNdHXMw+yB0NTcsnj0iWin8ph/aXgdcl5mX1pc/BTyk4TwlSX3QqPAz8yfAjyNirL7qsVS7dyRJQ6aNo3ROBT5aH6FzDfD8FuYpSWpZ48LPzCuAxvuWJEn95Zm2klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqROM/Yh4RO4BbgN8AuzLTP2guSUOoceHXJjLzxpbmJUnqA3fpSFIhIjObzSDiWuDnQAJ/m5lb5ty+CdgEMDo6etzk5GSj5Q3KzMwMIyMjg46xbL3m3r7zpj6kWdz6tWvunO7qmIPZB6GruWHx7BMTE9va2F3eRuHfKzOvj4hDgYuAUzPzkvnuOzY2ltPT042WNyhTU1OMj48POsay9Zp73elb2w+zBDs2b7xzuqtjDmYfhK7mhsWzR0Qrhd94l05mXl//fwNwPvCwpvOUJLWvUeFHxP4RccDuaeDxwJVtBJMktavpUTqHAedHxO55fSwzP9c4lSSpdY0KPzOvAR7UUhZJUh95WKYkFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSIRoXfkSsiojLI+LCNgJJkvqjjS3804CrW5iPJKmPGhV+RBwObAQ+1E4cSVK/RGb2/uCITwFvBQ4A/iIzT5rnPpuATQCjo6PHTU5O9ry8QZqZmWFkZGTQMZat19zbd97UhzSLW792zZ3TXR1zMPsgdDU3LJ59YmJiW2ZuaLqc1b0+MCJOAm7IzG0RMb7Q/TJzC7AFYGxsLMfHF7zrUJuamqKL2XvN/bzTt7YfZgl2nDJ+53RXxxzMPghdzQ0rl73JLp1HAU+KiB3AJ4ATIuIjraSSJLWu58LPzFdl5uGZuQ54FnBxZj6ntWSSpFZ5HL4kFaLnffizZeYUMNXGvCRJ/eEWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCtGo8CPirhHxjYj4j4i4KiLe0FYwSVK7Vjd8/O3ACZk5ExF7A1+JiH/JzK+3kE2S1KJGhZ+ZCczUF/eu/2XTUJKk9jXehx8RqyLiCuAG4KLMvLR5LElS26LaSG9hRhEHAucDp2bmlbOu3wRsAhgdHT1ucnKyleWttJmZGUZGRgYdY9l6zb195019SLM8h+0LP71t0Cl6s5zs69eu6W+YPZjv57wS496P57yUdX2Q6/WenvNi2ScmJrZl5oamGVorfICIOBO4NTPfOd/tY2NjOT093dryVtLU1BTj4+ODjrFsveZed/rW9sMs0yvW7+Ks7U0/ZhqM5WTfsXljn9MsbL6f80qMez+e81LW9UGu13t6zotlj4hWCr/pUTqj9ZY9EbEv8Djgu01DSZLa1/Rl/J7AuRGxiurFYzIzL2weS5LUtqZH6XwbOLalLJKkPvJMW0kqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCNCr8iDgiIr4YEVdHxFURcVpbwSRJ7Vrd8PG7gFdk5rci4gBgW0RclJnfaSGbJKlFjbbwM/O/MvNb9fQtwNXA2jaCSZLaFZnZzowi1gGXAMdk5s2zrt8EbAIYHR09bnJysqf5b995U/OQDRy2L/z0toFG6ElXc0M52devXdPfMHsw3+9VV8d92HPv6ec8MzPDyMjIgrdPTExsy8wNTTO0UvgRMQJ8CXhzZp630P3GxsZyenq6p2WsO31rj+na8Yr1uzhre9M9YCuvq7mhnOw7Nm/sc5qFzfd71dVxH/bce/o5T01NMT4+vuDtEdFK4Tc+Sici9gY+DXx0T2UvSRqspkfpBPBh4OrMfFc7kSRJ/dB0C/9RwHOBEyLiivrfE1vIJUlqWaMdXpn5FSBayiJJ6iPPtJWkQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqRKPCj4izI+KGiLiyrUCSpP5ouoV/DvCEFnJIkvqsUeFn5iXAz1rKIknqo8jMZjOIWAdcmJnHLHD7JmATwOjo6HGTk5M9LWf7zpt6TNiOw/aFn9420Ag96WpuMPugdDX7sOdev3bNgrfNzMwwMjKy4O0TExPbMnND0wx9L/zZxsbGcnp6uqflrDt9a0+Pa8sr1u/irO2rB5qhF13NDWYflK5mH/bcOzZvXPC2qakpxsfHF7w9IlopfI/SkaRCWPiSVIimh2V+HPgaMBYR10XEC9qJJUlqW6MdXpl5cltBJEn95S4dSSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVonHhR8QTImI6In4QEae3EUqS1L5GhR8Rq4APACcCDwBOjogHtBFMktSuplv4DwN+kJnXZOYdwCeAJzePJUlqW2Rm7w+OeBrwhMx8YX35ucDDM/Ols+6zCdhUXzwGuLL3uAN1CHDjoEP0oKu5weyD0tXsXc0Ni2e/T2aONl3I6oaPj3mu+51XkMzcAmwBiIjLMnNDw2UORFezdzU3mH1Qupq9q7lh5bI33aVzHXDErMuHA9c3nKckqQ+aFv43gaMi4r4RsQ/wLOCC5rEkSW1rtEsnM3dFxEuBzwOrgLMz86o9PGRLk+UNWFezdzU3mH1Qupq9q7lhhbI3+tBWktQdnmkrSYWw8CWpEBa+JBVi6As/Ik6IiPsOOsdydTU3mH1Qupq9q7mhvOxD+6Ft/Z08nwB+AfwWODMzvzTYVIvram4w+6B0NXtXc0O52YdmCz8iDo+Iu8266pnApzPz0VRP7uSIOH4w6RbW1dxg9kHpavau5gaz7zbwwo+I+0fEZ4GvAH8VEbu/fO1XwH719CTV90w8IiLm+zqHFdfV3GD2Qelq9q7mBrPPNZDCj4j9Z118MHBdZq4DLgbeWV//M+D2iDggM38GfA84DFi3glF/R1dzg9kHpavZu5obzL4nK1b4EXFQRJwTEd8ENkfEaP2K9EDgqxERmXkB8IuI2Ej1JA4A1tez+D7VN8rdsVKZu5zb7GYvJbfZl559JbfwHw3sAp5I9S2brwbuVme4R/7fp8fnAs8GvgHcQvXHVcjMrwEnADevYGbobm4wu9mXp6u5wexLy56Zrf6j+k6dFwFfovoe/EPq6yeBl9XT9wU217c/lGof1ar6thHgv+v5rAWuBl4K/D3wQWC/tjN3ObfZzV5KbrM3z96PLfyTgCcBbwCOB95eX38R8Mh6+sfAl4ETM/ObVK9qEwCZOQNcCjw0M3cCz6XaP/UT4DWZ+cs2w876oOOPupR7jk6NOTjug8jumLuu9/RtmfU+pYyIh1K9xfgysDUzbwd+H7gmMy+OiGuBd0TE44FtwFMi4pDMvDEivg/cGhH3Bt4HPCciDqX6Tv3/oXrbQmZeBlzWS8495N8AvJDqbdE7gBuA3xv23HX2To55nd1xd11fSuZO90v9HIZy3Je9hR8Rd6t/GOPA2VSHCD0OeGt9l98C34uIfTPzWqq3HQ+sn/j1VMeQAvyG6q3JXsCnqb4e9BTgOGBLZv52udkWyT0SEXeNiHPr5V0LvDczb4iIvaheSYcud5394Pr/RwLn0JExrzMfFNVxxFuA8+jWuB8aEQfXv7zn0pFxj4hDIuIhEfF+4Hy6Neb71P3yGKpdFZ0Y8zr73hGxb0Scw7B2zBL3Pe0H/Cnwb1QH/AP8OfCSevog4NvAsXXozcC6+raT6sCH1NPbgTVUHzJ8Fthn1nL2arKPbA+5LwY+WV/3NuBFs+6zuv7/pcCbhiT37jOgn061D28KGANe3oExn539q8AXqFbqoR/3ep77A8+jept9E7CxI+O+O/cXqLYCn9qhMd8beDHwGeBvgCOB04Z9zOdkv5Bqa/zoOt9QjvuiW/gRsTdwFfA04B2Z+cf1Tcfuvk9m/hz4Z+BlVL8ohwL3r2++BHgMcEdmXgh8GPgU8AGqraZfz5pPa6+6c3K/PTN3v3puBx4QEW+tX4n/LCLuDnwOuMegc9fzy4hYAzwDeE9mjmfmNNUr/O77DN2Yz5P93Zn5uMz8Nh0Y94g4kmr3wWOB1wA7gR9RfZC2ql7m0I37nNxnUO3T/WGdZ/0wj3ntJVQl9x6qv4n91Pryb+plDt2Yz5P93VTryFOo1pmxiNg8dOO+xFex84BT5lz3TODSWZfvBVxfT7+E6pTfg4B9qV657z3rvoe0/Uq7jNz3rrN9EjiZaiX73DDlrpf1YuCN9fTureanDvuYz81eX94XuGedb3JYx53qF/Yusy6fTbXB8ORhHvd5cn+Iamv/0I6s658B/qSefgFwat0v3xjWMV8g+/OBV9Ud88lhHPelPqmTqA72P4tq98LrgPtQnfF16Kz7XQQ8vJ5+E/CvVG8v/3KlfgB7yP1F4PVUH3ysnXWfvevndEJ9+c2Dzl3neEad6xTgW1Sv+E+m2s1wyKz7DdWYz5N9G/ARqiMUZq/YQznus/IdTLUP+XFU+1J/RnVM9NCO+5zcj68vd2FdfzHV7uJ/ojpK5UNUu0Z+DowO85jPyf4jqo2E+wNHDOO4L+eJfZ7qhIAjqF65TgO+RlX+Ady9/kHde9aTPAa466B+GPPk/ijV6clHzbr9HlQfhD5wyHIfVY/vX1NtDTwbeBfVcbivrEtoWMd8bvanU23Znzjs4z7neVwOPLWe/gjwtnp6KMd9Tu4n1dN7zbp+aMec6oiWs6n2Z78eOBOYrn93h3Zdnyf7GcDfAY8YxnFfzlE6T8nMt2Tmj4G3UB0e9T6qtyUXUL3KrcrMHwFk5q8z88rM/NUyltEPs3O/lWrwD4yIB0XEGVS7fX6Z1X7mYcr9I6pTpVdntQ/zwvq6C6k+jB7mMZ+bfSvwHeDYiPiDIR936iMqoPrA/Mh6+s3VTfEZhnTcF8i9V33EzqsZ4jEHHgBMZeaNVO9QVgEfZ/j7BX43+7nAbcD9IuLoYVvXl3wcfv7uQf2/oCr812bmxyLiOcBVmXl52wGbmpP7FqoTFa6meiXem+pIgGHMfXtEbKZ66wrVVs7xwObMvHzIx3xu9n2o3ua+keqoqX0Y0nGH6gOyiNj9bYQ/qK+7GnhlRJwCfGcYs8/JvaO+bldEPItq/RnKMY+IVcB1VF8x8A/ArVRnmb44M68d5nV9nuy/onqH+3aq3T13YYjGfcl/ACUi7gI8gersrqOpDp/6YGbu6l+85ubJvSUz3z3YVEsXEW+mOmHjWKr9fmfWW81Db072LwCvzP6ehdmqiJgGXpeZn9x9MtCgMy3F7NyDzrJUEXE/qkMU76BaZ84H3pTV2aVDbZ7sFwBnZOaKfxHbYpb1F68i4kVUJz7845C8lVqSrubeLSLGgB+afWXMOtPzwVSH9u7qQtl3Nfdu9ZmlRwH/npm3DTrPcnQl+9D+iUNJUrsG/hevJEkrw8KXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhfhfHCg3qEW8+c0AAAAASUVORK5CYII=\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "yearly_incidence.hist(xrot=20)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": []
+ }
+ ],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
@@ -16,10 +1359,9 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.6.3"
+ "version": "3.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
-