diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb
index bc059cd02dabe9586b4c321343c8972f7b3053c7..73b5965aba51a969eab94ee9d696b5d29d876751 100644
--- a/module3/exo2/exercice.ipynb
+++ b/module3/exo2/exercice.ipynb
@@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
@@ -14,7 +14,7 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
@@ -51,14 +51,1263 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 3,
"metadata": {},
- "outputs": [],
+ "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",
+ " 202531 | \n",
+ " 7 | \n",
+ " 5873 | \n",
+ " 0 | \n",
+ " 13316 | \n",
+ " 9 | \n",
+ " 0 | \n",
+ " 20 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1 | \n",
+ " 202530 | \n",
+ " 7 | \n",
+ " 7110 | \n",
+ " 3586 | \n",
+ " 10634 | \n",
+ " 11 | \n",
+ " 6 | \n",
+ " 16 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 2 | \n",
+ " 202529 | \n",
+ " 7 | \n",
+ " 6385 | \n",
+ " 3384 | \n",
+ " 9386 | \n",
+ " 10 | \n",
+ " 6 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 3 | \n",
+ " 202528 | \n",
+ " 7 | \n",
+ " 5584 | \n",
+ " 3123 | \n",
+ " 8045 | \n",
+ " 8 | \n",
+ " 4 | \n",
+ " 12 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 4 | \n",
+ " 202527 | \n",
+ " 7 | \n",
+ " 5667 | \n",
+ " 2850 | \n",
+ " 8484 | \n",
+ " 8 | \n",
+ " 4 | \n",
+ " 12 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 5 | \n",
+ " 202526 | \n",
+ " 7 | \n",
+ " 5872 | \n",
+ " 3285 | \n",
+ " 8459 | \n",
+ " 9 | \n",
+ " 5 | \n",
+ " 13 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 6 | \n",
+ " 202525 | \n",
+ " 7 | \n",
+ " 5953 | \n",
+ " 3698 | \n",
+ " 8208 | \n",
+ " 9 | \n",
+ " 6 | \n",
+ " 12 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 7 | \n",
+ " 202524 | \n",
+ " 7 | \n",
+ " 4580 | \n",
+ " 2558 | \n",
+ " 6602 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 8 | \n",
+ " 202523 | \n",
+ " 7 | \n",
+ " 4911 | \n",
+ " 2663 | \n",
+ " 7159 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 9 | \n",
+ " 202522 | \n",
+ " 7 | \n",
+ " 6837 | \n",
+ " 3940 | \n",
+ " 9734 | \n",
+ " 10 | \n",
+ " 6 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 10 | \n",
+ " 202521 | \n",
+ " 7 | \n",
+ " 4693 | \n",
+ " 2653 | \n",
+ " 6733 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 11 | \n",
+ " 202520 | \n",
+ " 7 | \n",
+ " 3083 | \n",
+ " 1535 | \n",
+ " 4631 | \n",
+ " 5 | \n",
+ " 3 | \n",
+ " 7 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 12 | \n",
+ " 202519 | \n",
+ " 7 | \n",
+ " 5084 | \n",
+ " 1997 | \n",
+ " 8171 | \n",
+ " 8 | \n",
+ " 3 | \n",
+ " 13 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 13 | \n",
+ " 202518 | \n",
+ " 7 | \n",
+ " 5003 | \n",
+ " 2718 | \n",
+ " 7288 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 14 | \n",
+ " 202517 | \n",
+ " 7 | \n",
+ " 6246 | \n",
+ " 3424 | \n",
+ " 9068 | \n",
+ " 9 | \n",
+ " 5 | \n",
+ " 13 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 15 | \n",
+ " 202516 | \n",
+ " 7 | \n",
+ " 6151 | \n",
+ " 3193 | \n",
+ " 9109 | \n",
+ " 9 | \n",
+ " 5 | \n",
+ " 13 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 16 | \n",
+ " 202515 | \n",
+ " 7 | \n",
+ " 5557 | \n",
+ " 3262 | \n",
+ " 7852 | \n",
+ " 8 | \n",
+ " 5 | \n",
+ " 11 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 17 | \n",
+ " 202514 | \n",
+ " 7 | \n",
+ " 4984 | \n",
+ " 2858 | \n",
+ " 7110 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 18 | \n",
+ " 202513 | \n",
+ " 7 | \n",
+ " 5964 | \n",
+ " 3608 | \n",
+ " 8320 | \n",
+ " 9 | \n",
+ " 5 | \n",
+ " 13 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 19 | \n",
+ " 202512 | \n",
+ " 7 | \n",
+ " 3855 | \n",
+ " 1995 | \n",
+ " 5715 | \n",
+ " 6 | \n",
+ " 3 | \n",
+ " 9 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 20 | \n",
+ " 202511 | \n",
+ " 7 | \n",
+ " 5878 | \n",
+ " 2747 | \n",
+ " 9009 | \n",
+ " 9 | \n",
+ " 4 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 21 | \n",
+ " 202510 | \n",
+ " 7 | \n",
+ " 2921 | \n",
+ " 1421 | \n",
+ " 4421 | \n",
+ " 4 | \n",
+ " 2 | \n",
+ " 6 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 22 | \n",
+ " 202509 | \n",
+ " 7 | \n",
+ " 3381 | \n",
+ " 1468 | \n",
+ " 5294 | \n",
+ " 5 | \n",
+ " 2 | \n",
+ " 8 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 23 | \n",
+ " 202508 | \n",
+ " 7 | \n",
+ " 2835 | \n",
+ " 1286 | \n",
+ " 4384 | \n",
+ " 4 | \n",
+ " 2 | \n",
+ " 6 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 24 | \n",
+ " 202507 | \n",
+ " 7 | \n",
+ " 4502 | \n",
+ " 2382 | \n",
+ " 6622 | \n",
+ " 7 | \n",
+ " 4 | \n",
+ " 10 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 25 | \n",
+ " 202506 | \n",
+ " 7 | \n",
+ " 3455 | \n",
+ " 1958 | \n",
+ " 4952 | \n",
+ " 5 | \n",
+ " 3 | \n",
+ " 7 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 26 | \n",
+ " 202505 | \n",
+ " 7 | \n",
+ " 2087 | \n",
+ " 1056 | \n",
+ " 3118 | \n",
+ " 3 | \n",
+ " 1 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 27 | \n",
+ " 202504 | \n",
+ " 7 | \n",
+ " 6895 | \n",
+ " 4466 | \n",
+ " 9324 | \n",
+ " 10 | \n",
+ " 6 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 28 | \n",
+ " 202503 | \n",
+ " 7 | \n",
+ " 2462 | \n",
+ " 1161 | \n",
+ " 3763 | \n",
+ " 4 | \n",
+ " 2 | \n",
+ " 6 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 29 | \n",
+ " 202502 | \n",
+ " 7 | \n",
+ " 5966 | \n",
+ " 2757 | \n",
+ " 9175 | \n",
+ " 9 | \n",
+ " 4 | \n",
+ " 14 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ "
\n",
+ " \n",
+ " | 1779 | \n",
+ " 199126 | \n",
+ " 7 | \n",
+ " 17608 | \n",
+ " 11304 | \n",
+ " 23912 | \n",
+ " 31 | \n",
+ " 20 | \n",
+ " 42 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1780 | \n",
+ " 199125 | \n",
+ " 7 | \n",
+ " 16169 | \n",
+ " 10700 | \n",
+ " 21638 | \n",
+ " 28 | \n",
+ " 18 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1781 | \n",
+ " 199124 | \n",
+ " 7 | \n",
+ " 16171 | \n",
+ " 10071 | \n",
+ " 22271 | \n",
+ " 28 | \n",
+ " 17 | \n",
+ " 39 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1782 | \n",
+ " 199123 | \n",
+ " 7 | \n",
+ " 11947 | \n",
+ " 7671 | \n",
+ " 16223 | \n",
+ " 21 | \n",
+ " 13 | \n",
+ " 29 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1783 | \n",
+ " 199122 | \n",
+ " 7 | \n",
+ " 15452 | \n",
+ " 9953 | \n",
+ " 20951 | \n",
+ " 27 | \n",
+ " 17 | \n",
+ " 37 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1784 | \n",
+ " 199121 | \n",
+ " 7 | \n",
+ " 14903 | \n",
+ " 8975 | \n",
+ " 20831 | \n",
+ " 26 | \n",
+ " 16 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1785 | \n",
+ " 199120 | \n",
+ " 7 | \n",
+ " 19053 | \n",
+ " 12742 | \n",
+ " 25364 | \n",
+ " 34 | \n",
+ " 23 | \n",
+ " 45 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1786 | \n",
+ " 199119 | \n",
+ " 7 | \n",
+ " 16739 | \n",
+ " 11246 | \n",
+ " 22232 | \n",
+ " 29 | \n",
+ " 19 | \n",
+ " 39 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1787 | \n",
+ " 199118 | \n",
+ " 7 | \n",
+ " 21385 | \n",
+ " 13882 | \n",
+ " 28888 | \n",
+ " 38 | \n",
+ " 25 | \n",
+ " 51 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1788 | \n",
+ " 199117 | \n",
+ " 7 | \n",
+ " 13462 | \n",
+ " 8877 | \n",
+ " 18047 | \n",
+ " 24 | \n",
+ " 16 | \n",
+ " 32 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1789 | \n",
+ " 199116 | \n",
+ " 7 | \n",
+ " 14857 | \n",
+ " 10068 | \n",
+ " 19646 | \n",
+ " 26 | \n",
+ " 18 | \n",
+ " 34 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1790 | \n",
+ " 199115 | \n",
+ " 7 | \n",
+ " 13975 | \n",
+ " 9781 | \n",
+ " 18169 | \n",
+ " 25 | \n",
+ " 18 | \n",
+ " 32 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1791 | \n",
+ " 199114 | \n",
+ " 7 | \n",
+ " 12265 | \n",
+ " 7684 | \n",
+ " 16846 | \n",
+ " 22 | \n",
+ " 14 | \n",
+ " 30 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1792 | \n",
+ " 199113 | \n",
+ " 7 | \n",
+ " 9567 | \n",
+ " 6041 | \n",
+ " 13093 | \n",
+ " 17 | \n",
+ " 11 | \n",
+ " 23 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1793 | \n",
+ " 199112 | \n",
+ " 7 | \n",
+ " 10864 | \n",
+ " 7331 | \n",
+ " 14397 | \n",
+ " 19 | \n",
+ " 13 | \n",
+ " 25 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1794 | \n",
+ " 199111 | \n",
+ " 7 | \n",
+ " 15574 | \n",
+ " 11184 | \n",
+ " 19964 | \n",
+ " 27 | \n",
+ " 19 | \n",
+ " 35 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1795 | \n",
+ " 199110 | \n",
+ " 7 | \n",
+ " 16643 | \n",
+ " 11372 | \n",
+ " 21914 | \n",
+ " 29 | \n",
+ " 20 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1796 | \n",
+ " 199109 | \n",
+ " 7 | \n",
+ " 13741 | \n",
+ " 8780 | \n",
+ " 18702 | \n",
+ " 24 | \n",
+ " 15 | \n",
+ " 33 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1797 | \n",
+ " 199108 | \n",
+ " 7 | \n",
+ " 13289 | \n",
+ " 8813 | \n",
+ " 17765 | \n",
+ " 23 | \n",
+ " 15 | \n",
+ " 31 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1798 | \n",
+ " 199107 | \n",
+ " 7 | \n",
+ " 12337 | \n",
+ " 8077 | \n",
+ " 16597 | \n",
+ " 22 | \n",
+ " 15 | \n",
+ " 29 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1799 | \n",
+ " 199106 | \n",
+ " 7 | \n",
+ " 10877 | \n",
+ " 7013 | \n",
+ " 14741 | \n",
+ " 19 | \n",
+ " 12 | \n",
+ " 26 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1800 | \n",
+ " 199105 | \n",
+ " 7 | \n",
+ " 10442 | \n",
+ " 6544 | \n",
+ " 14340 | \n",
+ " 18 | \n",
+ " 11 | \n",
+ " 25 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1801 | \n",
+ " 199104 | \n",
+ " 7 | \n",
+ " 7913 | \n",
+ " 4563 | \n",
+ " 11263 | \n",
+ " 14 | \n",
+ " 8 | \n",
+ " 20 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1802 | \n",
+ " 199103 | \n",
+ " 7 | \n",
+ " 15387 | \n",
+ " 10484 | \n",
+ " 20290 | \n",
+ " 27 | \n",
+ " 18 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1803 | \n",
+ " 199102 | \n",
+ " 7 | \n",
+ " 16277 | \n",
+ " 11046 | \n",
+ " 21508 | \n",
+ " 29 | \n",
+ " 20 | \n",
+ " 38 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1804 | \n",
+ " 199101 | \n",
+ " 7 | \n",
+ " 15565 | \n",
+ " 10271 | \n",
+ " 20859 | \n",
+ " 27 | \n",
+ " 18 | \n",
+ " 36 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1805 | \n",
+ " 199052 | \n",
+ " 7 | \n",
+ " 19375 | \n",
+ " 13295 | \n",
+ " 25455 | \n",
+ " 34 | \n",
+ " 23 | \n",
+ " 45 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1806 | \n",
+ " 199051 | \n",
+ " 7 | \n",
+ " 19080 | \n",
+ " 13807 | \n",
+ " 24353 | \n",
+ " 34 | \n",
+ " 25 | \n",
+ " 43 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1807 | \n",
+ " 199050 | \n",
+ " 7 | \n",
+ " 11079 | \n",
+ " 6660 | \n",
+ " 15498 | \n",
+ " 20 | \n",
+ " 12 | \n",
+ " 28 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " | 1808 | \n",
+ " 199049 | \n",
+ " 7 | \n",
+ " 1143 | \n",
+ " 0 | \n",
+ " 2610 | \n",
+ " 2 | \n",
+ " 0 | \n",
+ " 5 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "
1809 rows × 10 columns
\n",
+ "
"
+ ],
+ "text/plain": [
+ " week indicator inc inc_low inc_up inc100 inc100_low \\\n",
+ "0 202531 7 5873 0 13316 9 0 \n",
+ "1 202530 7 7110 3586 10634 11 6 \n",
+ "2 202529 7 6385 3384 9386 10 6 \n",
+ "3 202528 7 5584 3123 8045 8 4 \n",
+ "4 202527 7 5667 2850 8484 8 4 \n",
+ "5 202526 7 5872 3285 8459 9 5 \n",
+ "6 202525 7 5953 3698 8208 9 6 \n",
+ "7 202524 7 4580 2558 6602 7 4 \n",
+ "8 202523 7 4911 2663 7159 7 4 \n",
+ "9 202522 7 6837 3940 9734 10 6 \n",
+ "10 202521 7 4693 2653 6733 7 4 \n",
+ "11 202520 7 3083 1535 4631 5 3 \n",
+ "12 202519 7 5084 1997 8171 8 3 \n",
+ "13 202518 7 5003 2718 7288 7 4 \n",
+ "14 202517 7 6246 3424 9068 9 5 \n",
+ "15 202516 7 6151 3193 9109 9 5 \n",
+ "16 202515 7 5557 3262 7852 8 5 \n",
+ "17 202514 7 4984 2858 7110 7 4 \n",
+ "18 202513 7 5964 3608 8320 9 5 \n",
+ "19 202512 7 3855 1995 5715 6 3 \n",
+ "20 202511 7 5878 2747 9009 9 4 \n",
+ "21 202510 7 2921 1421 4421 4 2 \n",
+ "22 202509 7 3381 1468 5294 5 2 \n",
+ "23 202508 7 2835 1286 4384 4 2 \n",
+ "24 202507 7 4502 2382 6622 7 4 \n",
+ "25 202506 7 3455 1958 4952 5 3 \n",
+ "26 202505 7 2087 1056 3118 3 1 \n",
+ "27 202504 7 6895 4466 9324 10 6 \n",
+ "28 202503 7 2462 1161 3763 4 2 \n",
+ "29 202502 7 5966 2757 9175 9 4 \n",
+ "... ... ... ... ... ... ... ... \n",
+ "1779 199126 7 17608 11304 23912 31 20 \n",
+ "1780 199125 7 16169 10700 21638 28 18 \n",
+ "1781 199124 7 16171 10071 22271 28 17 \n",
+ "1782 199123 7 11947 7671 16223 21 13 \n",
+ "1783 199122 7 15452 9953 20951 27 17 \n",
+ "1784 199121 7 14903 8975 20831 26 16 \n",
+ "1785 199120 7 19053 12742 25364 34 23 \n",
+ "1786 199119 7 16739 11246 22232 29 19 \n",
+ "1787 199118 7 21385 13882 28888 38 25 \n",
+ "1788 199117 7 13462 8877 18047 24 16 \n",
+ "1789 199116 7 14857 10068 19646 26 18 \n",
+ "1790 199115 7 13975 9781 18169 25 18 \n",
+ "1791 199114 7 12265 7684 16846 22 14 \n",
+ "1792 199113 7 9567 6041 13093 17 11 \n",
+ "1793 199112 7 10864 7331 14397 19 13 \n",
+ "1794 199111 7 15574 11184 19964 27 19 \n",
+ "1795 199110 7 16643 11372 21914 29 20 \n",
+ "1796 199109 7 13741 8780 18702 24 15 \n",
+ "1797 199108 7 13289 8813 17765 23 15 \n",
+ "1798 199107 7 12337 8077 16597 22 15 \n",
+ "1799 199106 7 10877 7013 14741 19 12 \n",
+ "1800 199105 7 10442 6544 14340 18 11 \n",
+ "1801 199104 7 7913 4563 11263 14 8 \n",
+ "1802 199103 7 15387 10484 20290 27 18 \n",
+ "1803 199102 7 16277 11046 21508 29 20 \n",
+ "1804 199101 7 15565 10271 20859 27 18 \n",
+ "1805 199052 7 19375 13295 25455 34 23 \n",
+ "1806 199051 7 19080 13807 24353 34 25 \n",
+ "1807 199050 7 11079 6660 15498 20 12 \n",
+ "1808 199049 7 1143 0 2610 2 0 \n",
+ "\n",
+ " inc100_up geo_insee geo_name \n",
+ "0 20 FR France \n",
+ "1 16 FR France \n",
+ "2 14 FR France \n",
+ "3 12 FR France \n",
+ "4 12 FR France \n",
+ "5 13 FR France \n",
+ "6 12 FR France \n",
+ "7 10 FR France \n",
+ "8 10 FR France \n",
+ "9 14 FR France \n",
+ "10 10 FR France \n",
+ "11 7 FR France \n",
+ "12 13 FR France \n",
+ "13 10 FR France \n",
+ "14 13 FR France \n",
+ "15 13 FR France \n",
+ "16 11 FR France \n",
+ "17 10 FR France \n",
+ "18 13 FR France \n",
+ "19 9 FR France \n",
+ "20 14 FR France \n",
+ "21 6 FR France \n",
+ "22 8 FR France \n",
+ "23 6 FR France \n",
+ "24 10 FR France \n",
+ "25 7 FR France \n",
+ "26 5 FR France \n",
+ "27 14 FR France \n",
+ "28 6 FR France \n",
+ "29 14 FR France \n",
+ "... ... ... ... \n",
+ "1779 42 FR France \n",
+ "1780 38 FR France \n",
+ "1781 39 FR France \n",
+ "1782 29 FR France \n",
+ "1783 37 FR France \n",
+ "1784 36 FR France \n",
+ "1785 45 FR France \n",
+ "1786 39 FR France \n",
+ "1787 51 FR France \n",
+ "1788 32 FR France \n",
+ "1789 34 FR France \n",
+ "1790 32 FR France \n",
+ "1791 30 FR France \n",
+ "1792 23 FR France \n",
+ "1793 25 FR France \n",
+ "1794 35 FR France \n",
+ "1795 38 FR France \n",
+ "1796 33 FR France \n",
+ "1797 31 FR France \n",
+ "1798 29 FR France \n",
+ "1799 26 FR France \n",
+ "1800 25 FR France \n",
+ "1801 20 FR France \n",
+ "1802 36 FR France \n",
+ "1803 38 FR France \n",
+ "1804 36 FR France \n",
+ "1805 45 FR France \n",
+ "1806 43 FR France \n",
+ "1807 28 FR France \n",
+ "1808 5 FR France \n",
+ "\n",
+ "[1809 rows x 10 columns]"
+ ]
+ },
+ "execution_count": 3,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"raw_data = pd.read_csv(data_file, skiprows=1)\n",
"raw_data"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "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": 4,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "raw_data[raw_data.isnull().any(axis=1)]"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "data=raw_data"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "def convert_week(year_and_week_int):\n",
+ " year_and_week_str = str(year_and_week_int)\n",
+ " year = int(year_and_week_str[:4])\n",
+ " week = int(year_and_week_str[4:])\n",
+ " w = isoweek.Week(year, week)\n",
+ " return pd.Period(w.day(0), 'W')\n",
+ "\n",
+ "data['period'] = [convert_week(yw) for yw in data['week']]"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "sorted_data = data.set_index('period').sort_index()"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "periods = sorted_data.index\n",
+ "for p1, p2 in zip(periods[:-1], periods[1:]):\n",
+ " delta = p2.to_timestamp() - p1.end_time\n",
+ " if delta > pd.Timedelta('1s'):\n",
+ " print(p1, p2)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 9,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ/vudc+/Nzb4DIQkkQFgCmAAxBgEFAhIVBVQcdBSc4Scu4CjOOIIjjAs4oCKKCiqgsiiLqICyKPtOIJCQELICIftG9u0u59Tvj67qrq6urc9y7wm33ue5z+1TXV31dXd1ffWtRYwxBAQEBAQEyCh0NwEBAQEBAY2HwBwCAgICAjIIzCEgICAgIIPAHAICAgICMgjMISAgICAgg8AcAgICAgIyCMwhICAgICCDwBwCAgICAjIIzCEgICAgIIOm7iagUgwbNoyNGTOmu8kICAgI2K3w0ksvrWeMDXfV222Zw5gxYzBjxozuJiMgICBgtwIRveVTL6iVAgICAgIyCMwhICAgICCDwBwCAgICAjIIzCEgICAgIIPAHAICAgICMgjMISAgICAgg8AcAgICAgIyCMwhoGFRLjPcOWMZOkrl7iYlIKDHITCHgIbFXS8vx3/fNRvXP/VGd5MSENDjEJhDQMNi844OAMCGbe3dTElAQM9DYA4BAQEBARkE5hAQEBAQkIGTORBRKxG9QESvENFcIvouL/8OEa0goln870PSNRcT0WIiWkBEp0jlRxHRHH7uGiIiXt6LiO7g5dOJaEztbzUgICAgwBc+kkMbgBMZYxMATAQwjYim8HNXM8Ym8r/7AYCIxgM4C8ChAKYBuJaIirz+dQDOAzCO/03j5ecC2MgYOwDA1QCurP7WAt4pYN1NQEBAD4STObAI2/jPZv5n+15PA3A7Y6yNMfYmgMUAJhPRCAADGGPPMcYYgJsBnC5dcxM/vgvAVCFVBPRchBEQENB98LI5EFGRiGYBWAvgIcbYdH7qAiKaTUS/JaLBvGwkgGXS5ct52Uh+rJanrmGMdQLYDGCoho7ziGgGEc1Yt26d1w0GBAQEBOSHF3NgjJUYYxMBjEIkBRyGSEW0PyJV0yoAV/HquvUes5TbrlHp+A1jbBJjbNLw4c6NjAJ6IKa/8TZWbd7Z3WQEBOz2yOWtxBjbBOBxANMYY2s40ygDuB7AZF5tOYDR0mWjAKzk5aM05alriKgJwEAAG3LdSUAAgH/5zfM4+SdPdjcZAQG7PXy8lYYT0SB+3BvASQDmcxuCwBkAXuXH9wI4i3sgjUVkeH6BMbYKwFYimsLtCWcDuEe65hx+/AkAj3K7REBAbmxr6+xuEgICdnv47CE9AsBN3OOoAOBOxtjfiegWIpqISP2zBMAXAIAxNpeI7gTwGoBOAOczxkq8rS8B+D2A3gAe4H8AcCOAW4hoMSKJ4awa3FvAOwRhmRAQ0PVwMgfG2GwAR2jKP2u55nIAl2vKZwA4TFO+C8CZLloCAgICAroGIUI6YLfA32evxGdvnO6uGBAQUBP4qJUCArodF/xxZneTEBDQoxAkh4CAgICADAJzCAgICAjIIDCHgABPPPf627j8vte6m4yAgC5BYA4BDQ/WIKn3PnX987j+qTe7m4yAgC5BYA4BDYuQezEgoPsQmENAQEBAQAaBOQQEBAQEZBCYQ0BAQEBABoE5BAQEBARkEJhDg4ExhqcXrUe53BgeOgEBAT0TgTk0GO6fsxqfuXE6bp3+VneTEhAQ0IMRmEODYeWmaBezpW/v6GZKAgICejICcwgICAgIyCAwh4CAgC7BkwvX4ZVlm7qbjABPhJTdAQEBXYKzf/sCAGDJFR/uZkoCfBAkh4CAgICADJzMgYhaiegFInqFiOYS0Xd5+RAieoiIFvH/g6VrLiaixUS0gIhOkcqPIqI5/Nw1xJPnEFEvIrqDl08nojG1v9WA3RU9bQ/pr9w2E2Muuq+7yQjo4fCRHNoAnMgYmwBgIoBpRDQFwEUAHmGMjQPwCP8NIhoP4CwAhwKYBuBaIirytq4DcB6AcfxvGi8/F8BGxtgBAK4GcGUN7q0mYIzhlueWYEd7Z3eT0uPQU9Pu/e2Vld1NQkCAmzmwCNv4z2b+xwCcBuAmXn4TgNP58WkAbmeMtTHG3gSwGMBkIhoBYABj7DnGGANws3KNaOsuAFOpQVJyPrZgLS65Zy5+cP+87iYlICAgoMvgZXMgoiIRzQKwFsBDjLHpAPZkjK0CAP5/D159JIBl0uXLedlIfqyWp65hjHUC2AxgaCU3VGtsbysBADbu6OhmSgIaBayn6bkCeiS8mANjrMQYmwhgFCIp4DBLdd2Kn1nKbdekGyY6j4hmENGMdevWucgOCKgLAm8I6AnI5a3EGNsE4HFEtoI1XFUE/n8tr7YcwGjpslEAVvLyUZry1DVE1ARgIIANmv5/wxibxBibNHz48DykB3Qh1m9rw872UneT0dBYvHYrZgWf/4AGho+30nAiGsSPewM4CcB8APcCOIdXOwfAPfz4XgBncQ+ksYgMzy9w1dNWIprC7QlnK9eItj4B4FHWYLL7kwvX4Y/Tl3Y3GbsFJl32MD523bPOejvbS/j4dc9i3qotXUBV7VCLgXnST57E6b98pgYtBQTUBz6SwwgAjxHRbAAvIrI5/B3AFQBOJqJFAE7mv8EYmwvgTgCvAXgQwPmMMbGM/BKAGxAZqV8H8AAvvxHAUCJaDODr4J5PjYStuzrxrb/OqaqN7W2dmL9695oIK4XPhD/jrQ146a2NuPy+nmPsX7hmK9Zs2dXdZAQEOOGMkGaMzQZwhKb8bQBTDddcDuByTfkMABl7BWNsF4AzPejdrfHFW1/CU4vWY8Fl09Crqei+oIeAGdbiwl/tsQVrk7qMdfve0pFQWxkNH7j6ydoSExBQJ4QI6S7Ei0siM0op7NWQC281WIba8PYCegICc+hCEF9tNpY1JSAgICCLwBy6EAWuiQi8oT7oKh+GrmTuW3Z14MI7ZmHzzhBnE9C1CMyhCyF05eUgOnihs5R9To3w6Ew2klpgBd/sSeB3Ty/BX2euwA1PvVG3PgMCdAjMISeq2dtZ2FEbYYJrBJDDqFss5DP6vhOe6zFXPJr6LRYShcbIJhPQgxCYQ05Us+oXn7dN/SF75vR09G+NnOkmjB7kVb+reENXMiEWmENANyEwh5yoZl4oFNwG6Wdff7uKHt5ZEI+ppUiZsp4CIagG3hDQ1QjMISeqWTWK7zvYHDyR8zE1WFB9TZColbqZkIAeh8AccqIaY6QwSPu0EFaK9TX8VoNa8iAXQxOSQyFwh4AuRmAOOREkh66DeEwdkteSbTL1eaq/fGxxape1jlIZf5+9stukDle3gi6X8T4goNYIzKELEad9CLzBC+Ix1TJ76Y/+sSD1+5pHFuGCP87Ew/P8HQFcEs1Dr63B4xrHgjN/lU1G6FooiLNBcAjoajhzKwWkUZXkwD/wkD0jDdMzzSthVfJuVm2OkuBt3NGe/2IDPn/zDADAkis+nCp/ccnGTF3XWBCu08FbKaCrESSHLkQhZg49izssWL0Vtzz/Vu7rdDmo1JJSmeGxBWvBGMPSDdsrpFDTsK1qDV+faywEb6WA7kJgDjlRlUEaPTNC+pSfPolL7n4193W6CGkV1z/1Bv7tdy/i4XlrcdJPuibj6eotu3DRn2ejvbNcdVu2obCzvRSC4AK6DYE55EQ183ohR4T07sA/ymWGKx+cj7V12p/AJ3vtSp5uYsXGrsvc+u2/vorbX1ymtSvkRcnyog+59ME4p1KwOQR0NQJzyIlq5uw8uZV2A96A6W9uwHWPv45v3DW7Lu13ejCHlmI0hDs8pAwZiRdQfoj3V4t9JVxj4a8zV9Ssr0bB3TNXhC1SdwMEg3ROVOPymMcgvTvs+SAmto5S9eoVHUrlbLvq429uiphDe4U0VPKUBQ21mK6ZJ9nvJMnha3fMApA12Ac0FoLkkBPVSQ7Rfy/JYTfQK9WbxDySQ1tO/X81tC9et63yixX42p/eSZJDwO4BJ3MgotFE9BgRzSOiuUT0VV7+HSJaQUSz+N+HpGsuJqLFRLSAiE6Ryo8iojn83DXERzwR9SKiO3j5dCIaU/tbrQ2qC4Lz3+xnNxAcqoZrvvORnlqahFopJ3MQNOS6KsKG7ZHbay3ma1/mUK1BevbyTbvFgiOgceAjOXQC+E/G2CEApgA4n4jG83NXM8Ym8r/7AYCfOwvAoQCmAbiWiMSGydcBOA/AOP43jZefC2AjY+wAAFcDuLL6W2s8JCm73R9pT/No0kFnR1C9xYTkkNdzSH0HtUrVsXpzPuO87yKgGrXS/XNW4aO/eCa2XwQE+MDJHBhjqxhjL/PjrQDmARhpueQ0ALczxtoYY28CWAxgMhGNADCAMfYci77MmwGcLl1zEz++C8BUalQ5uipvJWGQrms3XQ7fN2ViiiY+qLM5qKhWcqgGuvtetXlnttBGRxdIDm+uj+I/Fq3Vq8O+8adXcO3jiytuP+CdiVw2B67uOQLAdF50ARHNJqLfEtFgXjYSwDLpsuW8bCQ/VstT1zDGOgFsBjA0D21dheriHCK8Y2wOOZ9F3lty2RzmrdqC6W9GKc6rNYpXkrtId03eNY2v5FCLpZLp+f/ppeX44YML9CcD6o5tbZ0xA28keDMHIuoH4M8AvsYY24JIRbQ/gIkAVgG4SlTVXM4s5bZrVBrOI6IZRDRj3bp1vqTXFFXFOXDdgI8u3WPRXHPMX70Ff3tlZe7rfCfWvI9O96zlsg/+7CncP2c1AKC9szpGVRHT19x23jn8Z48slGgy01CN5JCk9Oq+BcfusNjpLnz6+udxwo8f724yMvBiDkTUjIgx/IEx9hcAYIytYYyVGGNlANcDmMyrLwcwWrp8FICVvHyUpjx1DRE1ARgIYINKB2PsN4yxSYyxScOHD/e7wxqjmiHexJmDzyq3O2wO0376FL5y28y6tV/PCeLPLy93V5IgJspqVuS6S/O2d9sLy7y2ni1U4VcYM29HN39+Kd8zzIPAG8yYvXxzd5OghY+3EgG4EcA8xthPpPIRUrUzAIj8CPcCOIt7II1FZHh+gTG2CsBWIprC2zwbwD3SNefw408AeJQ16FKjGrIS/biPQbribrywbmsbPn3981i0Zmt9O5LQSC9UvMZaj7JK1FNCfWajpTaSgx3/+adX6hKzsmzDDky6/OGatxtQX/isR44B8FkAJypuqz/kbqmzAZwA4EIAYIzNBXAngNcAPAjgfMZYibf1JQA3IDJSvw7gAV5+I4ChRLQYwNcBXFSTu/NEucxw/5xVXiu4rpIc6s0b/z57JZ59/W38qY6rRRV5b6lB1wcxauUz0cl1iPW6W5+9yxHXqX3/d85YFrv/Buw+cEZIM8aehl6Cvt9yzeUALteUzwBwmKZ8F4AzXbTUC3fPWoGv3/kKLjl1PM49dmzd+mku+kfz1ntaFC90V0fJWs+G3JN9XgN2vuYrQq194ippzyfYryaSg5eXnLnSwjVbse/QPujVVDTW0fafq3ZAoyBESCPyFgCANzwiX6tZWQnm0OHhk19vm4PPqnd7Wyfum72qZn2qt+SioJ6PIGOQro09uiKI7LN1M0iL4EvNuUy8h4GE9dva8IGrn8T//DV/dt2Qb3z3RGAOAPr1igQowSRsqMbjo7ko1Erdb3MQsE2Kl9z9Ks7/48uYUyODWb0lja5uu1ZzXqeHa1p1e5fzNqrIBrydfxvCdTgP3kl5oXoSAnNAsqIXK7j2zjLGXHQffvrwwmzlKuaUpqJ/wNaYoX0q78gDPhPbcp4O28U0vYPgahAXUStpIiPFVDCB6YzPlUh8seRgqVOL+9Y9f7VdE/1xAGcF9up38v7XyzbswJiL7sPs5dVnmW00G1tgDpCCLfgY3tke6eF/+/SbxrqVoMXD5jBh1EAAwIDW5ip68ofPhK2bOHe0d+LXT76er68GGvsqKRWplTTPpRKJr+ThrVTNoxMqRC2z9exHxOhUMoG9k7VKYk+PO15c5qjpRqPlUwvMAVJuf49RXOkE99SidfHGLTbJgcX/6zdS1m9rw6X3zK2qjav+uRDPLM6nYsgdBJezfq624/0cajtzVZJqvVr3UcYY/jRjWbyoUWG7w6zNwSQ5RP8rmcDewbxB2qOl+rYaTXII+znAnJ+/rbOMNcouZ5VM2owxfPbGF+LfNoN0vfzvZcg2BFM/S9+276y2dVdH7n7zDn79SrexPiAVlXzgseRguTdbu8+9/ja+cddsvLx0E/7vY4fn6ltt1TXJVaI2K7yDjQ7JerL6cRkkhwaE+CjVMdzWWcZl982rvn3lpdsM0oKWrhonun6eWLgO7/vRY3jhzUyQegzTits2iVXAVnNf0ZUt6261EsnBx5XVBmETWre1TXvelg04U2QgZcn6aLFgIvUfc1fj49c923Cr33qjGluMikbLxBwkByQvNlYrWRY6lQRyffLXz6XKbDaHrpAcdP3JWLB6S642nlq03q8vw23rVsxbd3XgwVdX56IjD+pl4K1knk9cWS19VUFvHAQnlbV1lrB1Vyf6t6anAN0E9fa2Nnzq+uc5HXpCvvLHmWgvldHWWUZrczoO4p1sc0jUbY01sdcCQXJA1iBtG8x5h8CO9hJmvLUxVWa1OQjm0GVunLXtx25U9Q90OPw7/8TGHVnVVc2+QTXOoZImNBflzfEE+LmyVgOdQfqLt7yESZc9rElAmODEHz+O037xdGqXPdMkKFLD6IIq38neSuLZ1iLTQKMxmMAcIG0Y7zGI84rN2zVuoH7pM3J1U1PU62NupLFfC+ar+5jvqmCS8FEr6ejd0d6JLSnbj74dXVbWxxbosxrL9/TG+u14ZflmNBVJOq+nTzCHid97CGsVO907W3JIbm6zZjGTB8Hm0IgQBmmPYKG8E5wuRsA2CLpifMj3oLufaj5mq69+5c3WHZXccq3up1K10nFXPoZ3feefzvZX8HgVn7gRLQ1SmWl126spmUqeVNSM72DekLJTbtxRXf6oRrPXBOYAjUG6hu9It/G9fRJgqf/1QC2aNjEQq0G6gQa/TZ3i30Zt7scrQlrT1duZZHb6l/LrJ96I2tC1q5RqjdbSsSk5ZbEgSxdKm9or3hmQJYdq7zNIDg0I8VKSHDTmt3Rvzs1wdCut7sqO6duPGu9Rs6jkBmlDbqcqKclBzEbPTKSdPq6s3lSZ4Zd4Lwt5DJtigeS2VQbiO3427+jAFQ/MR2cd0obXC/LjqHax0EiLJyAwBwBSnIOHWumeWfk2afeJStXVr+cwcbVdjRrA595qjUrc6GvxIbpWekd8/yGvVNU+aiUbVIcHM9yurPrFTD56ShXeyA/un4dfPfE6HpxbnZfaro4Slm2wx+nUCjKzrHbl32C8ITCHD/7sKVz3RLS5euzVUec+bftGxHEOXSU5aO5WXRzmWV37qMzyXOODYgXcoTYr8XQruoWDjx665KVWMlP8myffELUAABu2t2vdgL3SZ+RczOjoq3SSbOuMPJ06PRJT2vD1O2fhuB8+FrdXT6SHXnV0B2+lBsO8VVuwbENksLMFCwnk9eSpVHLoTqh3qPO4qgT1mJCB2m26kxcqJbrtHn0o6/BJvOdNFfD/bnoRX7z1pYzUohtb6oSkm6DkxYzpfmx2CV/PMPEer318MV5dUXkm4MfmR55YlQQk5oVsc3jotbVVtdUAn34KPZ45yBCrAPtHmu8V+mTCTNevrB9flMsMv3xssZUWdbI996YZ3u3bU0B4N5MLxQqYg5eXjgSd/77KqG7UJGr0ytflrOFbKeprKV/sqIZubdBeXvuA0RFBajOTr8nRptL0wjXbcOrPn/a7yIKuWGjJksOVD86vqq0gOTQwYoN0Dd+Rz2otXd+tVrrjxaUYc9F9+M695uR581dv0a6yn1q8HrOWJemF51S4QjN7K5mvqRfDq0ytlI+Wv87Mqozke61Gx+3joeZHr73O+m1ZFZfP4jplkDb2nNSpeMVeIwGwKwXJWkqtDcYbAnOQoQsWqhZ5VwNM+a/DN/88BwDw+2eXaM8/Nn8tpv30KW3Uppr0b/7qrZk69cqTVpuUFVlU9H3moGX9tjatykieA7/0h5cqIKJikirCo/Ozag8ftZKXl5NUx5TjyYVaB192xVxbS4p3O+ZARKOJ6DEimkdEc4noq7x8CBE9RESL+P/B0jUXE9FiIlpARKdI5UcR0Rx+7hribJeIehHRHbx8OhGNqf2tulEgt14pt81BU2aXHNSD/Hidb3c6f1V24te1es+sFemVb52WXvUa+2KzJmf/0jMVRz63euEds3DbC0uz7Ul3ZLIp+zxJnwVELScOOaraR60kFxldWaXjX8cGctGmr83Bq5q7Hf6/K9Q0texhd1QrdQL4T8bYIQCmADifiMYDuAjAI4yxcQAe4b/Bz50F4FAA0wBcS0QiE9d1AM4DMI7/TePl5wLYyBg7AMDVAK6swb1VjFq+Ip+oVBligJiq1MMX+g/Tl8Z7TfijEj1/fQZ/c9GPlscX6lNGuDBzqX6XL3lebTLQ4DPhxe7LHrYoP9hry1HVqgZIN0Et35gsHCpRJ/qi1ksSU6LHWuILt1QvMQrsdsyBMbaKMfYyP94KYB6AkQBOA3ATr3YTgNP58WkAbmeMtTHG3gSwGMBkIhoBYABj7DkWzRI3K9eItu4CMJVqqczzRJ69dv2hMQJaOoiZg6FKPRww1AddVZyD1sCtP1eJGkHXvq/kkIpmVdqZt8qcidaoZ5ca0amdomvd9+jzTisZk359q55FWXzp1pdzt5Nq09cgXeMvvtH3/lDRYLwhn82Bq3uOADAdwJ6MsVVAxEAA7MGrjQQg75m3nJeN5MdqeeoaxlgngM0Ahmr6P4+IZhDRjHXrKlsF2lCI4xzMbynvAK48CE5fq1r3PJ/Vu889GleQ3eCtJOf1saFfryQ9tUrnLc+/Zb6wqlxTPiqj+KV3OdTxpBsfrj3EAViD/Xxvq9Y2h0ZLR+HCbssciKgfgD8D+BpjzJbwX/eGmaXcdk26gLHfMMYmMcYmDR8+3EVybtRDVNG9b3ugmL1OV4ietf5Ik5RV9VIr+UoOyXGex1gwcMJa2QpcCwLXuUr6NNV1XVuR7b+LJYdk684Gm20daDR6vb4qImpGxBj+wBj7Cy9ew1VF4P+FK8RyAKOly0cBWMnLR2nKU9cQUROAgQDM25DVCckm6vmv3dbWiTEX3YffPfNmqlwXDW1bvbtsDvUK7JFJqle+oXqN/UradTtnJqhGz+5Dmp904dFQBch6K9nrFwv5HRxNk576HdRardRok60LjUatj7cSAbgRwDzG2E+kU/cCOIcfnwPgHqn8LO6BNBaR4fkFrnraSkRTeJtnK9eItj4B4FFWxyxUJ171OG7VqBGqScoq3Pd+98ySVHleycFlc6g0b03UpukjTU9Q1Wx4r+vBlZYk1x1pKvvSK1+aJ/utac7KGyNgriNoMtfJ84xY/N991fb2tMrIdY2nkJb2DJOO33fgcKlcvaq23GHuii27VRK/RmNmPq/6GACfBXAiEc3ifx8CcAWAk4loEYCT+W8wxuYCuBPAawAeBHA+Y0yEl34JwA2IjNSvA3iAl98IYCgRLQbwdXDPp3pgV0cJb6zbjm/f/Wr2pEf6DBPENerqJ28QnBjLpo/UlpcJiNxYr35oYab83ldWYuzF92OpJlgrvWEMcOk95uC6amB6rra9qn0wdlhfz/6zZfK7mPKDR/DBnz2VqWPORFpbtVL1ldRr3FW++7fXUr9daZ6aPCUHeQOjlFeXpNtTFzryY95/uN87teHffv8irtJ8C42KBuMN7j2kGWNPw8zSpxquuRzA5ZryGQAO05TvAnCmi5ZaQOTA79tSzJyzRUgP6tOMTYadnhhjsbeKqp/Wps+w0BcnYTPaHCwXA5h61RPacrFD2Rvrt2fOVRJhbII29xH/Xy8D4ZC+LZ415dVs9F+mafWWXVit7GIG2LyVeBu2RIo+DEQkW3TWzAef9lbyjYBUWkwwueyqKJUZxFbSMhOw7fsgt1yNs6J8ZTU5mroKRFx6bzDu0OMipMWHrDMy2j6M848/AAAwaUwc64f2zjIYY/j77FX42h2zAGj0ppombWNA2BRqbXPYwT1O+jRnmaKLJh3yfLrJMzGvFMdcdB/mGNxBv/Whg+NjnyyyJugenY8o7/LM2mXJ/unzuhbyKHV7+gx/+KRgEWhXIuZd1/guJOR25GdctLgT+7zH79w7F//z1zleNADAU8qudI2Iei+eKkWPYw5Wr0GL7vfYccMAAKMG9wEArN2yCwd++wHc+vxbWLR2W1zPgzdYJ4GYORjqVKqX3NEeTWB9NBJTmVW/ao0nJGsdexuf+90L8XEeYcblXXXJ3a/iP26bmerfJ/BM7kEHce2uDrMuxkfvf9NzZjfaOBlkRVol90VqQkFXP02eL0Yep3KbMnNRFzryezSN/98/uwR/mJ6NVlca6hLUeqXfaHEZPY85wDz52gx54qMQg16oZ/42e1Vq9KvisDZXjYW+eFcwk0E6x/JCJqWdGzN03iaMsaoHum6yffDVVakkf64eZPXDgN7N8XG1rrW3PP8W7n1lpZI+Izr2eZ6m+bAcMweL5JDDHqqjROj487yfeBx7XKJuY+s2SPtNGfK4l59xwaJWkvdfcNGuSjwyuip6trNGS33RTJ6x0hXoeczBIjnYxPGi4uYqPio1ACsjOeRUZbhcWae/+bbxWhtck0u1w1y3V/YXb30Zp//yGe9st/KzkwPW0lsxaq7znA3kS0U2Wt27uOPFpfjKbTOd7YuJ1DZRVbsaFKlBKmnFR8pU35trvvOXHPR0yCYLdTJctTlr7zHhsvtec1eqM6rdlEjF7uit9I6ELfukzwquja8WezUV0y6SapuattSuL/7LHHz77kiP6pIcfvTgAv0JT2hz+tdgUN718nLRgXffWZ1zMnO0+PpMIgdzkPq7matydLf+zT/Pwd889govW8aLrs889AlU4iwgrjj7xhfwrRz6+YgGO8G9NWpJVztykzbJIXW9o/0ZS3y3Rq0f2ncjN9lK0OOYg03kTs5lTxaLQq0U/Y4lh+aC0fhmakstue3KR0vlAAAgAElEQVSFpbj1+aUol5k1WnbO8s1YmWN1pb1HrSRjb2fcHv0yZSa9r3al7Kk3l+dBX68YF25+bkl8LNPWyg3zPmolo1rLI1YiD+PVqjOLYlGSrb/XgFZDOxHeWL8df1Tek2vydz2O3gaHBpVBy+3Iz1g2SKuurOncV3ZC+rU6HS3rjmrigXQIkkM3QzacqgPQJjkUlZB8MTCaC5T6qH0M0qZVacmw2hJYtXlnttATNsZXZsw6cfukp2ixTOamFMrZyTBpw6TbfvZ1nVpN3/fmHR3pmA2pO6EOrMZbySQ5yEZ/XfOfP26svkFNXaHG0VHZv4IJ0n27lU1QA1qbU7/l55pSK1kkh/Q5e3829VZX5eysvVqpps1VjZ7HHKQD9UNJjNXZ68RcpUazEpGSekKNgtPToTNiulaxrkHv51Ovu850JoLPBNpsWeGm+zHD9L3L963bkc30WC68c1bqtxzjMZAbvH2+R1dWVvW+5BWw6dl9/MhRmTLd5k3Jc62NOtB1hWuCMp1uUWxvMm2yN1dKraQsvPNo0BphkS0WiO+W3NurQYhz6GaI519m2XUrkzmHgiJRNAkpBmPS1k6Q5wOWmYPvQFm2YQd++dhiMMa8smfqVU12ycEHYrWsa0ZnrNZVThuezdKYL1QDpxwVfxx3TfbLUqunYMsu8bzTbci1Ta3/+Mx3Zcoeem1Npuzafz1Se/3jC9bi9XXZgEbAxaDt9+t6HLYULKbf29qS4NEmWxCcrFaqwpCff2+SyiCYg2/UuAtBcuh2mNVKcQ2d5ECEAlGs+omvVeYNH28lE2TXON/L/v33L+JH/1iA+au34st/SOfd181pOmblXC16EDOojztK2e2tVBkbMOY+8vjavGwOhg5+wlMzZJowMLl0m7ogzDQmjx2CsYY0Ep/73Yt6ohxw3e7MpZUaetMNy+NMXrT4SFWAB5MyfCGX6NLi1AkdXK1UK/tY9T6DtUWPYw6JP372RdhWVYUCoVigjG/zX15egWcXJ1GYnlolLTbtSHLi+zIVEdz2h+lveUWDatVKjnVantTU9ihf1eaQhmkSdqmQ1fPCnmBLUphIkPa2ffq3q5Xc9ZPr1IpSMsgcAynPO1Dxfw/Mt543jYXsjnLJcZukVipa1Epy06771T3XzTs7tPty+CwSKoGQHPJ41tkQJIduBpP+q+9CRKrq3lGBoklH+LTLdV4xpH0A8kkn7//R4xKdfiNFtD+wd7OjZgTdSrlctq8YfSipRYS0aae2vOsycY+2jJzi+XoZpCUKTjl0z3Q7LMtaa7lvRJLRVt/QsQcMy9VXNerDI/YZZLxeHedlgxQsP5vss/cnTk21snjtVkz47j+1dWsVrKZCMIcxnokfXagXE6sUPY45CDAGjPufB4znVBSI0KupkOjPDe9RVY2YXrdr8s/7EfdpyXqu3Pj0m5kynfvdik078YP7zStGH738mq2Rfj+P73o2zsHZjRbqM+8sRxHfRluH1LfPc5bp+urUAzPtqG3YPHLs/WiYg9SPDmOG9UH/Xm6vpVokV2wqkJk5qL9TUkDyo0laZauSnekaHXZ2lFILGl224bifujGHqN0TD442wfyXSaNt1Z1oLNbQA5mD1+pK85oKRGgpypKDSZfs159bp6pp27OejJfeSksElfhm+zyzC+94xUmQa6IUz25bWycWrNkqldsnNr1txWIIR3JP7xk7xNo2kF7Rq7ZHBh2TMyeXs/ejts2cDLOpUMBeA1t5vwlNKvIk4zOBQMZxr656U7mVEHn03Hrue3DA8CRmJuNKLh2v3LzLmpYEAM649tn4WHWlNdFSS3TGBmnCngN6Vd1eiHPoZlS66xYR0Ku5aE2VoG/LoFZyXKcGL7lgmkTumZV2/ezozD8A81xhW6SZXIcFhASwRpM2Oy81neVyHMWug/gQezX7xHAkdTKSoVOtlEdyMJ8ztdJUoGTvc0tXok41ExCR+f2qxaKftVt2YebSTdi4owPHjhuGM44YiZPHR6o5dZ2iPqvlG83SgIpWQ3CejrZaQURINzcVUCAz4/RFg/GGHsgcKtTLFguR5CCSg/m+yFq+78cXrs2U5R1QHRVk98qX9M1qdbBeW7FBWlNWKjN02PZZEP89bk1mIDqHA9HGlR8/HE9+4wSMHzEgPq81SBv60docHHmpmooFL3VczEDcVY0gm9+2Ui7ue/IPHgEALOaZiwsFwsePHMnrmCUHIHG28EE17ruVQqiVWooFEKo3KAfm0M2oRq1ULJC0U5seqgrE14DnwnOvv41bn09LExNGD3Jep3bTUUFUZ55Bn0dyUPHW2/qVomvuIyLs0T8S6z/33jEAIruDVdec455am5JVqc1VeXCfFuwztA+u+dQRUjdRhXmrtjj70Y2dRFWkJ7ipQH4RwXEKkyokB5taKWM/MPcj6FXfj3pJnnFnW5TUy84bq5WKlAmG9cUxVzwaHwe1UjejUrVSgSJ9s0t3m/VGtKuVZPdVG97e3pYpG96vl78oy6v5JJTLXurXx66OksONUvldw29h9JBon40BPKVEqWRnDq7stzKam5K3mpUckgBCMen1l/TfggTd9qOpdhjDIsnOkj2vL28qkldkcbwvhLuqEfnUSjZa9NJQdnykS2xBnlYGUKc5N1YrcemtEsa7QtqJr7FYgwdzIKLfEtFaInpVKvsOEa1Q9pQW5y4mosVEtICITpHKjyKiOfzcNcS/JCLqRUR38PLpRDSmtreYhu392bw+iAgE8jaquvoT5cs3+uVL0rUjD8YfVpmtNW/fOmza0eEQ7/3aycx1qYeqnwmvP3sSfvWZIzG0XyRBuLrKo1ZKpw/Pru4F89Q6DHje9K3Pv6VVo7iEgmadWkm7uHHbJVwoFrLjf8n67djVUUKZMZxz9L6x1GT7TopxKhq7tKFO+HfNWGZss1bJD/OgU1YrUfWT++4oOfwewDRN+dWMsYn8734AIKLxAM4CcCi/5loiEjL5dQDOAzCO/4k2zwWwkTF2AICrAVxZ4b1UDfFqTv3509rzBWkALHlbn7pg7sq0+mCxtEtcLWhzlWXrVD/gvCd1std16Zh1eOnbJ6nUaPsd0rcF0w4bEU+Urg8tz3fY4lArJZKDXz+6sjmavY59SOzdXPQyNicusZWPh6YCpZLNtXWWcPyPH8dXb4922WttLnolNIzVSg5afv3E66nf3/lbdg8HodqpkQYxF2TJoblQ0Kbw3tHe6b2P9W6XW4kx9iSADZ7tnQbgdsZYG2PsTQCLAUwmohEABjDGnmPRE7gZwOnSNTfx47sATCV1edZFcL0cIooH4W+efMOrzYUmVYGlq97NRRy0Z39n27VYafioJHwHbbR6sqiVKiB3aL9e6VxFjjbIc4Usnp1PpttUniciXHb6Yck56Y5905XoYPIwTgzS+nZ6txQzexDratZCcmgqFlIBZYJRPLVoPRiiZ+PTT5H096Re809NrikVM2JX7a6XHISk17uliL69muJ92mX8x22zcOrPn/bKe9Zo20NUY3O4gIhmc7WTSEs4EoAs+y3nZSP5sVqeuoYx1glgM4ChVdBlhVXtofz+2wXHpn4XKtArdpQYxg7ri0n7pjM32jLANhWz4rs+3QewZkvWFqHWscGHD/vecaR2s7Xjvidtu6Q/lvsVSPZcdkgO/L/Y9MenLhAZUT8zZd/kHEv60uWG0pGh32hKXxYbpA2307u5KEVRm2GLgfjX9+xjuTJBJDkkM1jMFCVaCx6Sm2AgtZgMRXxFLdSZebGzPZrwezcX0aeliO1tWbXgyzxQz+ZWLVBqsH1CK2UO1wHYH8BEAKsAXMXLTXFatvgt79guIjqPiGYQ0Yx169blozhu2H9lqwa2FMhtc8i2ydCnpYixSoi9rZmWYsGrn1qsiHxENF0/6vaoAtZJuU4faYp5QKhYgMNGDjBcgXwzhlRVjXNh8mlPyUGXmttEDTnOtzYXvBgiGVbrgH/0dFOxkDLyi3sTnjoEOZ7C3E5BY3NgjOHpxebcYKZ7E6opu5dc/SSHpgKhpamAvr2asL09Kx0UPFVoQBT13UioiDkwxtYwxkqMsTKA6wFM5qeWA5BjyEcBWMnLR2nKU9cQUROAgTCosRhjv2GMTWKMTRo+fHglpDskh/TJXkpgTYEo9ybgDHr/dRuai4XMYK90fLsu8yFN17cu6IiB5dL95pFI7OcTyCvXw0eaXX297DWaGxdxLnKd2OZQYT+AOcWDS7IrFgoZVZqObpu3ku/4bCpQKk5GdLOtrROdZYaCSGsPe8qKmIFIdZ5ZbN8bXTbWD+6TeIOVYsnBsuiztlw5dnWU40VSZKzP1hHP3ceF/MI7XsE/566uJYlVoSLmwG0IAmcAEJ5M9wI4i3sgjUVkeH6BMbYKwFYimsLtCWcDuEe65hx+/AkAj7I6WmZsDTMWGZAEWpXo2ciVz02aTH6ZsZQhW6VDvzWkn4RSG8mhMrWSeQ/u2tocgHw5l9LBXhY9tIfze6LDT+qqOazkfrRpuC03/cxFJ0p9aZ5nqp3kWN76tFiAZHPwmJCrkBzeens7lm3YiVnLNnGa0m1FaiW7jUTuT34FW3fZ92CQV9XyzoRlD8mhXjaHMmPxvZhUzuJ5dHhmVrjsvnm1I7BKODN2EdFtAI4HMIyIlgP4XwDHE9FERON3CYAvAABjbC4R3QngNQCdAM5njIm3+iVEnk+9ATzA/wDgRgC3ENFiRBLDWbW4MRNcK4wjv/9Q/LulWMBfvvxeDOIZT4myaYZ1KJVZnOO9zOyThm5QNxUom5RMM9H5jHm3zcGnDb++y4xhyXpzyoPMPdTom9XZJOT9uPW0uBHdd6Qy2XdoH1x2+mE4aK/+Sh3YJQdLRyMH9XbWS9RKUYWN29tTW5+KfUYAe+xGPAYtq1sXXl4aMYW/vLwcE0cPyoxdgn7iN/Unj3H5Hf7y00fi/D+m9yaRJ3iZOSRBqbVflJjw9rY29GlpUpiDfkEn7lXnyaSDLYFgV8PJHBhjn9IU32ipfzmAyzXlMwAcpinfBeBMFx21gnWcMGCXxOGJCEfukxiSC0Qogedwl9J3qygxFj9YYahTx434qRO/CwVCR0cXSQ4VqpVMDOOLt75kbMc3UlVlpi5vJbm+rGJxGSl9jdaMRTaW48ZpVJlMMp5pbQ7ZsqP3y/pbmFKmqwZpdZIpFihTx9qO5lyTx34E/Vub0NZZRntnGTu5iicz/iS1kj3OIeovnVI9eXgikWAKUnPN0uY6JQ+DdK0lh6Muezg+Hto32uSqQAa1EucOprmif2sTtu5yezJ1B3pehHQOm4MK0wBQkd7u027I1jIHIqebn6ksW8deiUAY1Me+F4TaQntnGds1AVtlxqy587NZOPV1q9EqJrp1+xZGZWaXLEQd0ZZJ/ZaKkNbU0b33zx69b6Zs5eZsskGmkTrVSUbW84uudPagJM4hcwrNHqLDQxe+Pw4GvH/OKgBZ1VwvnoAOsE/Iop30DnHQHgvIrenUSl3lraSOzbe3RxkOTCpn8W5M2ZCba7RRUD3QuJTVDZWPFHkA2CawTsWjo0DAX2ems6OKy7/8h+xKu0lj3KqUObg2OiFy57tXB/3fZ+tTcLjoeeHNxM+gVGZGI51amifxnuwt46LH9fgSA6+ZBlkC0QbBaa6pZGsF0Y76Pgs8cl+GjtkntpgsRa4Jali/XthrYCuG8FXy9lhySNeTmYPt2YsUJ/KKWWaCOgO53F5arSTUs/5qpdnLNzm3QzWlgtEtigCxoDP3XWZROpdGC3SzoQcyBzNc702WHGzzaUma9MqMGfYLjuqs35bNraSTNHTd+UQ/uyZ+QnoFOKxfNi+9+ly+fucr2rZc4vu1jycRr2dc+wy+cIteBZXRVriM5jqbg0tl5KFWius6zsVqJW0/eskwd7+8HXV3O1mtJO5ZlwbmPfsN4XWy/TQb3JKTPng9fjBmaJTHSrWLRcwhTYsOfTh9svOHjsHLkMe6bjMltTeRFly9FgA++otnUntBqFi0ZisOvuRB3D8n6zm008gcbE4akW1k/2/dj/+4fZZyvnGZRY9jDtUYKWWPBNvglz+aMjOsFC2dFQp++nkTCbKh0yU5FCht/P7MlHRAVLT7l98AzpP9crZla9XqpDvJ5mCpV2Zu1iq/Y5XBi7TcsiurPMP96jNHxnSo8GUO7xs3jPed3IsqbRULWU+kgjLg3n/gcByz/7CYXhU6yUFeMBQVemNaFBVXS1PiVmsbCzoVlyvQUW5v085kQWWSHOSfebOyPrEwiqF6cUnWo9703ZtUx+J5CzrlxJfXP/kGNu6we2l1J3oec7Cdk16ubtNwMQBKDk8YOf8MY2Z9tQnFQkFjc8h2aBqo8sfV3mkwHsaV0zumNSkTS6FA3lN1rVZBmWYcj+/Reck+F3JAmGsh4Jt/SVftrMmj43aSxHsJoYdw5lGp6+htn5+Cr510IG83oUHVXadsDrxMfYf9WpuklN3ZvlqKWXpkw7cp1uJH/0gne+zVVPSSHHTGcZlh6nNUJbWH9k2k23ji13R3wQkHZK71gfgedLYb02KLDDFQpDBuGZff3zhuqzr0PObgKTn82zFjMueFK+v1T9nzKsniMmORJLCfGiFtub6oEVH1aiU95I9rV0c0Yo2DGulnoq46iwZdqg61ypuvc5EU0HWxSEpuKEdIu9wb/W0SLMOfUh89E2UyzeYVtPqMdRg1uHfMROTJWe+tlO5LlUwIdv6q81aS+8lsjcr7UZMFpm0OFuag28AoZZA22xz+48QD8L4DE6+xssXmMG7PaEvSvMNStKkzxZQMdjJTnIO4FZcE34jogczBPmEIaDd7p2gVvcKRZvs1aWOXyCBNuOGcSQCAyWOGZPpSMbhPS3ZS0dQ3jTf54xL78JqCvtT7VD/MYg61kqve/zt2rF87OT/nPi3JCk/2VrI149OH7AmjDof4p9SNXCXxINJIDh5qJZWBCHpVbyWCtJ8DE5OawhwcCfF0aiVZZdTMuUPcDdKqEoFezZK3ksWtn+R3xJHap9vyfEYP6ZN6ziXpHWX7cTMqHUSbuvdkSoNh8mSU4252N/Q85uBZTzc+xerANbFc8MeZ8bEIgttveD/M//40nMG3SFTbuOnfJ+PlS07G3ecfg32G9vHzzbYMVAERt2EWh9Vr07+LilrJnjffTCqQnsRN0KmD5MlON2384IzDk/PS5GRXIXpIDvz8K8s3ZxIcyqoRUS8db5FuQ8bw/u7N6OWJSZbuXlNSwreXypmsrBnmINGjG1fNDrVSn17Re0smuuh/hjk0Fb0dAgBlMSaVa11ZpWcsj29Bw4//mVZxHT5yoKRidJKSgpjIdRKeKTleZCfUSA78zm542qxtEF5gjYaexxw8B4pugFaaeE+01dpcNIr3owf3xpC+LZg4epDWLU7HkEyTsdxHqyO/vrpK003M8rVtljQArmfjs3jSuaC6dPTyRyxPTjZGVioz/OQhxwZJ/PL129qwfpvCHCTVyFOL16X6Buyuo2qUtQ6yKkc2SKt66o5SOdOXKg0+OHe1NQhOHQPrtrZh/irzrnRiQlYXHOk4B+PlWqQN0mZvpQKl64oxN391mt4LTjzAqtqzQdyXTnKw2hx0KkTehC131Hv3r1sS6qrgjJB+p8FXZaETbYXo6DIwnzZx7/g4qq+hQyEjLVbbvS8AYOrBe2DN1mzglIp9uduhzmUWyNKmis0qo9plyRzpZg4+q0mW0au7vHvSq87kF0N0/7q9qZ9YuA4vvWX3dbeNFVk18usn3jDSUak2IS05mO0+7Z3llLQEZKW79s6ydV8I9fEec+WjxoheIBkj6iq6pakQMzXRz/gRA1JqVl1/Uf3kWLcWEM+RKP2cTbEyxUKibsurphT3p5McOnPbHNwqxA3b/bYK7mr0OMnBF3qbg3CBjAbB0L4t8Yb2AiMH9UaTtOxjYFpPDHUYqatOl0FaGMd1kCf4dp7ywAT1NlWDW3MxrVaypRX2N/CaUWbA9/+e3fHL1kb62Yl2+Cra0I77k62NO2alqRtS489C7IF79nf2RWQed1Hz6Q7U8bLXgHQ6C6F2GaM4WbRoJAebqkWeTNN5lnQGafE+KfWc1biPVD+KGswH81ZtiZl9sUBYt7UNX7zlpTgxYF5XVp9x9uzr9oy03YWexxw8v1XdSxWr6ET/CXRqRp5sPHx1xRZskTJOmlZwqjFcnZieWJDev4IsKq5T35UkzW0vlR0rp/SdjpBiJIDIWCnTagoCAtzMwTZRDuTJDcuMZVb0eSbYZKK009Or2T307fsj8Drp0uz5CiUHVZ9veod7DWzVbl0qY1i/Xpm03iYMVfTfH52wN6765EQAWc+byWOHpOrqXFltKSXkUyk3ck1QnjgbMbrkjq0M3KLaM+FnDy+Kj0tlhl88uggPzl2Nu16K9iozqZWMqXV8uEODoscxB937m3nJybiQ+5QL6FQZieSQQP2IiZJORA6aF5dIk53HpKETUR9U8rwTzPlazj56THzc1lG29qXeppqmvKmYVmnoPo7vnXYoAP+4geUbs2qeAb3NGk5ds/L2q5SalBPmy6BfhQLRRKbiOx8Zn+7XSJGeyetsDpVIDgfv1T+1ZwEZiBFSq+qJpPbZVCBJ0rGrldSzZ04aFTNuAdPua81F2a1WqJ7Mq2n5jLzGGjmoNw7de0BqQynTY/TbN9tYJQN5b/hymWWYqnnfDT0tefdyaST0POagebeD+7ZkfJpNBmn1elUHKRsPdROpaaiorny2FdF/feBAFIiMq5jekleQqr8/6ZA98fzFU6W+zP0AUTCg6GXzzg7M4mmbZYweEtk15I9Dtwe2mFSOvfKxzLnEeJi9J13ZtY8t1tKbzq2UjU8QUAPFAGDs8H7OfiWCAZg9bkT7rvQlOjz4tfelYg/kMSVDzelkS0In6phsTwLqPeuCQRP30WxdlVHZXUzN/U7Zb2jqHU1/M1K9yEF/OhpkVLJvtryYuOqhhfEe46IJVXIWJJpyK1WSR6tR0POYg2E9mPH317zV7OqAUh8/kT6jqg0iD87oIYk6p0DJpLJsww68rXjK7D+8X6TSMhjHBrQ249H/fD+O3GcQ2jrKKZoLlE6JrOqb1d/N0palE777T/z3n2dn+hNXPCJFKq/blt3b2roS5438Q/KsEShrJmCdhxKgREjLFyh96ebs9yhqEhvBum9eHkPFKphDtq+orQWKR46sagEcqhzexo0Wl0ogO5GmmBRvQyxK1FtrLia5lcR96+ICEskhOacLfBRFjDH8z1+j/cTEN6Zep9r+RN2ojv87UD/7f8xdE9MAAG8p+y3sO7RvfJ3tuato5JxKAj2POfjaHCySg9yGGPxnHjUKT/zXCVEiO2GT0LabXs2MGNSKDx62l+Ijn6gsjvvhY6n88YDI10RGtRIA7De8H/q1NqOtlFYrZaJnFSLV381crWTrS7T5C2k1r/PA8Nmp7MI7XsnQmNprWKmvIuW+yAzvAPqPU02XINeYevAe6TY0q1J5YhFOCWIiFTEee3jEOOjAGMMpP31SKeO0CBWXhu6E3ui/S3JQn4tOwjLtodDcVMjQUonNQdQRRX+bvUqql36nyT7W2XsRjE1nFzTB5TatGsAF3UaDtKE5mRl+QEoS2EjoeczBUJ6ZJDXTiro6kK85ev+h2GdonziK2gR11dRZZpkB6dpNq8S/EFdIflOBUC6ztOSgvHH1LtXtGkX0rM/ELmOEZsMWG2OWJyH1+5QnDnEs15erpyUHfUZcNdmgCeKe+7YUMVbxzNGtfuUxI0sOG7a3x3sgq+34QJ4odfTFMHAHxtLPRYacikLUlSGPTSGtiTZUKbylWMhEhtvyDclQ30f0HUVliyT70va2znScg2Wzn76cIW9rMztR+NAmt59JpS9dp4+Q1rcnv7uph+yhrdPd6HnMwVN0sAXB3fbC0rhMeGz0b+VbiUp9iHFx0iHJyiDjOlpm2WR3DnG4f68mq0E67ou3IQ9adXWsDt6Rg/qkfgvPEdviS/es/vLl9+LWc9+TKrMxmEljEpWO6rueUitxemW1UkrRp7hS6qU3v/gDeUIwSVgmvbl4p50lhhN+/Hhc/uvPHuXuWKUXBpuDXEHqX/+c9QsOVZJRr5SZw8/POgIAcOjeA6O6GbVSOglge2fZGjQpX55x7JDal2noLDOcNnFk/NvmMtuP7xuxvc1/pzWT4CAYlZoGI86ESwZjv6GftDq6MQ0TTuZARL8lorVE9KpUNoSIHiKiRfz/YOncxUS0mIgWENEpUvlRRDSHn7uG+BMhol5EdAcvn05EY2p7i35QJQW9t1J6dUAAPj15H9z5haNxkuD+lAx6MXl///RDM22JcdRZYvG2iXI/gH7AX3b6YTj+oOGRQdpgc5DbYSw9oC89NfLI+cyUfXDrue/JTHrHjhuGJ79xAv7xtffhxIP3wIRRgyI6LdxBN7hHDOyNY3nKaQEbX548dnCm7Avv2w+AolaSPkZdu/LK1dRfFFTm5g7JQpwZ1XFyK2OHJ1JBge+zUCqXsXlnIo0N6pM/VYJ4jxn6FAYYG4q1bUT/dWNKtkGp5+WJeZ+hfXDCQcMz0oFMp6wSHX/pg5nIcu0NQLONpvQdyQGBjDHsOSBLr15yqIQ5mFb66f7i8nJynS23kgq5GZ9cW90BH8nh9wCmKWUXAXiEMTYOwCP8N4hoPICzABzKr7mWiMRS9ToA5wEYx/9Em+cC2MgYOwDA1QCurPRmfOCtVjIYMnUfxOSxQ5IVreRaIlZNstukOql0lssaySGrzxb4+JGj+Efo1qUKA/ovJVuAmJwuO/1wHDtumPY+9xnaBwft1R+//dy7Y3dCn0RqLtgkB51rqai/txJ7AZg/qCSVBNdPK9W+f9qhfNL2UCtJhldTgkLGGPYe2IozjxqFAa1pl8+mgtmjLA/Mj1dIqFyFVTJLDqINlZy9BrTiyH0G489fOhrnHgSy8OIAACAASURBVDvWGrkv+koklCxFcj4j272rHliqtEgSdyhKuZ/UFm2Sg1CJmqKodTBKDga10od5TJHRIG1oT1ajyWvDRvJucjIHxtiTANRdL04DcBM/vgnA6VL57YyxNsbYmwAWA5hMRCMADGCMPcei2fVm5RrR1l0AppJuKVorGFeTbhQMul+1HTFIhPFP9tdW/eNLZZYa/KIfQD/YxECav2qrc9CLdnSblqj0GM/z0zZG5OvLbZsne2kCn8SmRScclNXJykMknd1T9MVtDiD88BPvis+PGNhb65MubAFfPn7/eKezuH3GjB9tmSUOAio6Siy1+50OE0cPsp6PadAMXJPkoAugUzOU7j2wFT//1BH46knjAABH7TuEb0+bvlZduMjqHh1NvvmM5HaArIo0Yh4sQ0PGzGJhVE38uzKN3TyeZLFaSSHgm9MO5vQKZ5X0ebMkwpx1uhuV2hz2ZIytAgD+X3y9IwEsk+ot52Uj+bFanrqGMdYJYDOAumWiMruyuq9VPRIGa1QEwnjY1lnCNY9E0ZYp5qD001lmmQ3ebXl5xLkFa8yJ0VR6XSs4G+JVqaUN22rnlEOl7RoNnPX5i6dq00Z/8t2j4+PDRw5MnTPZWxL1SSI5fHLSaJxw0PD4PCGtAvjQ4Xvh4a+/HwDw39MOxgUnjuP0Jm1lJcvkuTCwjKHfF+9TVG/6m9IvSkSZyOrZ2lzE3JWb4132fvLJCfxyyiw4iAgfmbB36rnrnCl06b9j5mCTHBxRybLBeVdHCZv4jmjCJpO2OSQ0qpNz4jmVlP+ILwZEqnFT+ph5Ss4nwLz4i8eC8h2I52OS9o0G6bKeOZRZbdyfa4FaG6R1T4JZym3XZBsnOo+IZhDRjHXr1umqOGHTQ7ug6hV/fOYEbTsMLN5kB9BvpiKaKWltDtF/7S5i/OR/feDAzLkMLeQebK67TiQHfRu//PSRVoPasQckk59uZQVEOm+dC6FtRSVvTp+2OUTXbNzejn++tiapE5/PqgcLRGmvHGWCiySHNC0rN0XBUU8vXs/HRIWrP49ViamGoO+iD0ar1xMO2gN3vpiszcbvPSDuQmaaNlIyK98Mc4jGZanM9DYvQz//Mml0ppro6vgfPY6rH14IINrSNKYl3aS2XfFbJlvsH93cJCQHvuovM1x6T2w69VaHAsl4M+/nIOhR1WN6yPehjq22Tn/vqnqiUuawhquKwP+L6KflAORRMArASl4+SlOeuoaImgAMRFaNBQBgjP2GMTaJMTZp+PDhuipOGJmDx0BR1RGD+jRr6yxZvwNfvX1m5pyOjs4yi8VfgXgVonUDjP7rtjDM1o0Mr/ZVv0OtBLPk8PNPHYEPv2tERnLQBSQBXNWTQ60nk6aS+a5RaUlCQNzPf/7pFQDAG+uidAixCoYIO9pLWLkpyWibXe0JepP/KgOcsl8k3A7t22JUO6lpJyqFeI8qRJqU1uYiBvVpjiLCOZ39W5viZ0GQ1T1Ccsj2o1Ob6tRKADD1qsdxx4xlUBEvsqSGzj12LK6UVHuif1Fj9ZbkXQhJRnYakGlVn4POQ0swtDjWhEuZi9Zuw83PvRXX0y1ITF/Kr554nfejPy/6lM+/tnIL3pLSccgopSSHtNeYvLDsTlTKHO4FcA4/PgfAPVL5WdwDaSwiw/MLXPW0lYimcHvC2co1oq1PAHiU1TF80CU2x/UMKh25XB9FTViwZiseX6CXbJJJRrI55HBlTRm+FVxy6nj86YtHJ3XhNgy6FryiGx1zkFUUMvr20jOuMgNmLc+m3zDRIUtz6qPo35rkYpJPxfp3Tq9IH6JmaV2xKdnNT1WByHYhnXeU3L+QznTvQ6zobfBZvJr492GSqk3Eboj7lPMpUaRLAyAzSR0tWZuDel/iG1iiSYMenef9aMrUvnTXJluj2lULAjqbg6BZbGIkbHPqe86r69/VkaiKVeik/Q9d8xQ27ujQ1penuB3tJdx7wbH48OEj4n4aAT6urLcBeA7AQUS0nIjOBXAFgJOJaBGAk/lvMMbmArgTwGsAHgRwPmNM3OmXANyAyEj9OoAHePmNAIYS0WIAXwf3fKoX5Jfy78eMxez//QCA7ASn05lmguA07TvVNPx/mQH/d/88tJc03krxKiSfnv/cY8fi3VK8QIHrkMVEecYRIzPXuOgV/dgYjEqLLWXAx659VntOd43t25XJUVVENujOq/ZK2aMsiZfQS3diUyF9XIyVlFRfLqhP/7sfTbtHF4hQKifvKcpBlEy0ohvh1qkbWqoHEaCzOejH5eePG8vP8+cir4wtC7DfPfNm/Ft1E5elvaROVHjLuZMBADdxSUD+XovxfUfqQmGQVt9z3iXo755Zkvr9mSn7xMd58zjJ6qmFa7Zir4GtOOWwvQA0zv4Ozs1+GGOfMpyaqitkjF0O4HJN+QwAh2nKdwE400VHrSC/u/6tTejLcxupw1f/8aRXVtrVl+fHPmvpJvz6ySRvvNoPkBVhv3z8/vGxzyb1gpmVpAkjS2/laqUkAEgv+ajQfTdCBaVfycp0Rv+Fe3Be9aB4xnpVgiI5SJ49JsmhGE+CerVTRH9Utt/wvrF6K0OXp81Bvd+zj9433U4holdMgMVCEpBWIIpjLdZubeP3pumHsvEUWalWnwHgPWOH8vPR77TkoH+5W3d1xOoaDTHaYkHfcePMkd3ypU2FJB5Ifc+uIFIVqi1AVhu6AldVyJ+TkG734x5zyzbsSEmF3YUetxOcDN0mPDYIT6T4d87VrnxeHkSqp44p1cGJUn4fH29fwczEykmXJ8cp6cSSQ/ZDOuXQvVJ1XLTZ/MC1UpimHbExim1jGx1++PF34Yan34xtBTIyNgepXHzEOsMsgFiVY1soWPNSeTB52bNHLku1Q1ESSBETMvWQPRObA2WzifpuTJNhpgbJITmdXdjoGCDB/lySd5B+tiYBVqZJprm5WEjUSsq1eb2CdB51AnlTtMuSlXgOIhuBT2qXrkAPTJ+RHMvjPiM5aK7N2Bw8dam68/KA163OgOyHIE8kvioLxlxuqA56pRWywOVnHIbXvncK+vAIVFcyPwGba66eNjNdpjsytbfHgFZ860OHeBkh5UA6kwFXzjPEDDYH8ew6Os3Pf+ZS+1algJ9dQtgc9h4URQ9fcOIBqettyQzlNlSoEzsB2ofPUhVUVV+2PhFgW7jHqj0GPCxl+9XR/dbb23HPrJXxb/k+motJgkr1WleA4tHKQsKWlM8k7QORB9YnjhqVKrvl+bfi47HD+qXofqe6su4GcOtCAbMPdzktOmTgKznITMQ3t1I6OZ2POoJSzEEbtOSil/+XJYcBrc0xY/BpQ0DnZyBuSauWsTScTsYnt+dJjIUuWbqL9d7Kyy5Kk4HL5mBbIZsMllka7eeLBUo5HzRLAW265+ib6sFXrWRSvwEGxw1QZg/qn/7LxNR5IPpa52viEWSc/stnFBqT46ZiIR67PpKDafMmQC95q33qs/0WMvub3PTskvj4X3g8T5JwMzCHboFpIvEzSCs2B61ayXNyskgOoo1lSu54nS++DSIZmJgwKhlzNm+lpB91dWnXF6eulbxT8sB0L5WkH+ilJiOU+hDjQG23ID0Xk7eSKFM3XMpLb0ETuaxrp1SW7EuFRMolyq6SP8SNny5afA3SifE7+v93Kc22aVdFlabUPiPSZPt+KXOsbvKVGexdXzw69Q02FyhWK2UlB90Wv2kaZajMQY7cdwWuqgxSHhNioZFsEJVtozvQ42wOaUNZcqwOBLNBWn99XMfRv+78mi3p5GRioH1T2VhH/lB9mJAI2kskBw09ToN0BHVTo3Q/TlIA+KsyfJCSHGCmzQffPy3tJyGrNGKbg6rj5ze9YPVW7OwoGSSf6L9Vt+5BcEuxYM1uKugpSc4HTQUpTgBIrdIP3LMfLv1INhmkPr25+jtrtD5gj344nkegi/qzliUuyzp1DCH7XFLR2vx/pNpL6sjHH5mwN15dsRlLN+yI71vO7gtwyYH3o44+V+JK9Z3LzGzJFR9O348i7c9duTnVTtEyv4gY2ELMHBqDOwTJgcNPr5v+bZsQTNCdf2JhOiYiXpVa1EpeXi58lWeTHFytiAErfxiZ3eNstEjnbBHfLjqu/Hg6iCrtyupJiwa3fX5KnH5CbYOBmW0OvEDojvUG6ahQrFz//ZixmTo+jLVXU8GYAkKgSNHeHWKSLBYS9Q9ResfCwX1a9BO25R7i30i/x7HD+uLhr78/Fbzm2646OevSzEQMWq9CLPLx3azOvBKaioQOKUJaRl7d/nAeqLbXgOxeJWom5Q9f83R8bnt7ZzrFvPIdiG9A/G8UyaHnMQdp/XDQXpIeUDXYaQZOdkvNLDym7Ey9b3/4kFSN2DClfDypfaY93pzIlWOzObjaET2u3ZpEsZomSoEh/ZLJduSg5EPSMSc5nbUNew9MZ2bNa5A2QTtxiT5YMg5MgYqm33I7pTLDqe8agUs/Ml5Tx0Ny8GEOhXQeraKsVgIgDyXXLnotmiSIUqXUt5H1VNPQZuivQ/nG1DxPQDRm5Un8IxNGxMfCQ0tcN0GTxLC5kEgOqhpLZ5C2LTR6c/XjDedM0lwXXbhlZzY9uOqwsHzjztR51c06eCt1E8Rzv+HsSaldsNThqxs46gRQkZpGc1pduepW60CSDiCixU9ySKXP0BD86gq7sU8Q/La0vaRpVzSBT09OgoNOPHhP3Pb5KTh4r/5afazIGOu8nYxYrv+A8mqVtAxeTEwsWfWbAhUF1PxYQPodmYyZPryspamAnY6o2Ug6gCI5JFKPrKow9RmTaJmbCkTYJu2PoDala1trcwDQoTC8of2yiSyF5LDfsL5YcsWHsd/wfkm7nAGKVDIDWrNa8uamxOawaUc6uEynvhHj6qozJ2TuLQ6I1NzjQzyP1//dPy9zTs2CkHGdJoU5NIjo0OOYg8CYYem0zFmbg26V7aPn9+s/HeegX5WqYq+c2tvP5pDPCP3zTx2hbQNASud9yIgBSh11okz/Pnr/oWgq6vMD9eslPuh807rcVDXeSjY9u7xqVZMnZu7Zoa/XJV+M+nfT2FIsxNuMmlAspAMemwqF+LkUiFLxHUbJwYMWgn0XM21aDINYJRuE9x/eF8P6JfmF5GbLjGnbEEb4wTzH2ZZd2VV7sVBAZ5lh/bY2fPHWl1PndAvAEmOYevAe+PhRozLfcpKCxfyg3tZEN3eUyqm2TBlrY7USP72jvTPD0LoSPY45mOZK9YXrRDt1EtdvC+iQHPh/eWCqq07x8aqDV56A5MH29ZMPxMNff1+mL+FdNWpwpJI5akx2tzWB48YNw0cm7K2hN+rItqOX36Si97iRcxRZr5fOz1u1Ja2HzkmLqV21rMwSo6lLraQmT1TbNkkOPsysuVjIBLGpKFJarVQgYMzQvhjUpxnfnHYw3jVqUCzxmbqMx66FpAJRSh2UR72mlrVLuq7/d9x+qfNxuu1SGeWynvkWuRFepM7frJlIhceebvzqbA6dJSZ50KX7FGNBp4oV0osuL1JnOb2XudHTjrcr1HYnXfUEJn7vIX3lLkDPYw7xm0m/ePXj1Ul2LcrqT/uOHd+6mGRkY5wxzsGQOz6qkxyfNnFvHLBH2o9akFJmwPgRA9C/tSml7lFhmqRELh41r0zqWk+JSve8RL9DHNtnyj189BdPGwOYcjMHTZnYkP7umSviCSQr3dl/R7RIaiWD0dRHAiwWKOX6qJPwhNG5zFUYRITeLUXMuvQDOIFH1hcNk17SBj+wSJtESZZTwM854Z5ZK7T0ymoldbMn4da6atMulIwR6JE0GjMHjf1KeFeZNmNSUWYsZkSqN9WTC9dH/WpGTeI2nD3XWSprN7NSIdTGYmG6cvOuTJ1Vm3fiIz9/OmUDrBd6HHMQUN+hujLU6SNblK0sq/H+uVNKd6xOHGKAZW0OesnBFLkpPp7OMsOYoX21A/cL79/P2sZGD7HWL1pbv8eu2AltjGLHsKGjxPCjfyww0GJf4WfpypYt3xjFl/z+2SWSgdeuVnKlo2gyWP5V+uTUzXGdAqUmKl2qeNkgbXqXYtIzPRPZCGxC5AGX/p2iVdO26qotrpPVSuo2sf24RLl+Wxseem2NdpVf5AxR5Dj64vv3z9QRgas65qD7xmX7gOo+vIyPC909intp0SwCOkos3s8asDhTFBIaZNz1UrRP2obt7bjk7rmYs2JzXFZP9DjmIHtwyFAnaJ3koHL/oX31O8HZID7QuSsTQ7A6cZjC6OUVempVaph4IoN0NKhME4YQ302Sw1aNHldF317ucBmirCT0hffthy+fkKR5UG0Z6evND1ZW76nP7ONHjlKrqy1nSuRnYcpLpT5O3eQlt2Nk4Mpv3bMsUHo3M93KtUiR4fVXT7xu9GyKJQftWTnK11ABPg4XuskxSw8hvXLfe1DaPVQ8b7H/+aK12zJtRAwxmlSH9++FL2iYg8gvpp/Q0zfKGMPr67bH77xNUREJ6Vb3CP7nQ5HHoTZ3F9LvlTEWb4H7I2mfi6Lhu/8vvjfJCT9+HA/PiwzftmjtWqEHBsEJDw5VrZSeYHUrQdnF79JTxxvTAtjgkxXUJ87Bx621QIStbZ2YtWyTdkUK2JOJAX6upv09mEOBKBMpfPCI/unAPmcreshPSe3jsjMyiYBTcDHzToO3kvoetd5t0qM1qZVGDEq76GqjlImcrqxEwIy37Hma1C0tM21Yr9bX8XF4uPDk7K6FRGlpSDZGAwmttr0NxIIjUhvp6xTIvN+GOgk/zuON/jE3moBVySGR2LJtjePpMXo1FzMqH8ZYZnOq5iLhoxP2xpnSDnmxt5JBZSp/i6YFYS0RJAcO9WPX7TQmMwfTStA3bkBGP2VyNUkOpvQZPioLkapZhdhK0ZSeQZ0UxR6/Mnw9pzqVids3YR9gn7j2GZJ4nsnS3eA+zRl1hU+7si5cTPom1Z+AKy7GtNL771MOSv3WB6e5mUOxQMYxqbZtdGUtCLWSGcL1OA9OkLIJy5CZg0q7GNO2sSWM8Ca1ESBsDize1lWGGoTnkpKFpGNWywFgDAtWp/d3LzOWkQjbOsuZeBIiij2wVGQC57pAcui5zEF5trLB8aMT9tZujFP0UhPkf2lyEjsdbbo+ZanFFGTkM2kLtZJptXLJqUngVp+WYpymOy+IKGMANKnT9Neb25ZF+VGDE0bhFwuSrTOQrxCbCklyOH0cQ3KskxzS3kr6T03d7lWbGbUAtMlGYJ10USA0ezOHyiWHhWvS6h2bfQIAvnT8/th7oCaiGOnJWf2efCY/kTLEJBkAwmMP+PQN0zPndLmVZJhyMZmen+hLnQPUdRdDpCbUGamFB5aKKx6Yn6lXb/Q45iCQybIpPex3jRqo938v6Oun2lWK7zhvSm7aTAM97coqMQejJ4y7L8EUTZLD4L4t8QrHpYKyIdIxpz9G9ePI62lkwgQu9Zm8qGQ7hK7Gv75nXwBRtkyhWlA91YD0GNAZN9N2IfOEYrpGrtOR8hDKolRm2O6MhbDbHORNjnzhqnrhSQca08zYJQcfaZT4Zkt6bybRj2ls62wOMlR+3+mSHBAxy8/cmGZE6mTPGEO7RnIAuB1Fs9AQG4MJdIXNoccxB9NY9ok+9jIwKte+RzFQ+Xx2xv5N3krGlaB7AIngLFvWT9F+NcxBZ3NQ+7RRm0ciE/dk+n6+OnVc0q5hFT60bwuIEkOw7kOW3/WG7To3yixNtjpR39k6nWXmnITFJkg2iAnFFQSXI27SCXPeo7QkqY5hL8mBkrTqLslBxrEHDAOQTU+jwig5GMZigQiPzc/uHa++uzKL1Eo6lWeRyLnPBNA1KTZ6HnMwJFKTB6PRniBdZBq6Lo7u807VJiaMGojnL56q1PGhF846YkVsT8kt6la+WikUspKDqmbKmzTPBDEh+UQCmz504bIpmINOBSA/0s+/b6zmfH7JQcfo125RfNorfEwum0O8h0KOecdV1RZTsUvadlOVfk0GfBlCBWPyRgKSIDgZN/97tP+0cxI2SA7GYUrAa5q9J9T+S2WG9pJecigUyPotCkzad4izTrWoijkQ0RIimkNEs4hoBi8bQkQPEdEi/n+wVP9iIlpMRAuI6BSp/CjezmIiuoZqNUtoYHrs8grHpI6QP2STF4U8CVz3r0dq+ne/+IwnVbGQynUf1UmOfZiZ6eMRH6FN/SqeR7NHII+xDaLMjmjj9067rloN0jlGRHMsObglMHO7UfBUm4U5iIn86P2G4uC9sm64ctPmIDj1d7aeOle4pKjffe7d2nKXt5KPpuK4ccPclTxAUDOs6tVKNpJEgFupbDdIqxJAgRvv1Uk4u8JPF4jIcKNB31Cuvj8xd5hsDj6b/aSShtYJtZAcTmCMTWSMiVSFFwF4hDE2DsAj/DeIaDyAswAcCmAagGuJSMhV1wE4D8A4/jetBnTpYTBI++mHk2OdgRJIJoGJowfhg4ePyJyvRHLQ76DmnvjlezR57YiJ1Camir50evc8kCWHmZecnEngt457VE0cPQh3n39M6lye1axLd+4zCYoV51dumwlAr1YSz8XkoSYnyzMn3nOrU9RJzMUoTYubeAL2GC8myA4KQD77hK2vrEHaPdZkrz6zzSGdGvyqMyfE/WVsDlB/pyG87VxZbVWok/1OC3No8pAcrtFEyNcD9VArnQbgJn58E4DTpfLbGWNtjLE3ASwGMJmIRgAYwBh7jkUj7WbpmprDFOeQWol7eP+YPYoK/H/lwk92V7osbLvY6cpbm/WvWkhMtgEp7qWPRzyDCarNYbAmgFBkyR0xsDWOnBbIs3Xik9xfXZd+AEi/X5shU+5Sb3OI/psmiyP2SXJZ+Ux2Km0CeWwzpjYiGlzqNh89v3tsJnXN59SJNBNHwvuxqX7EWqXTKjmkFyW9W4pxf6prtUty6KxQclDb3cjtU/01WWRFGnIbusAWHfVT5fUMwD+J6CUiOo+X7ckYWwUA/L9wch4JYJl07XJeNpIfq+UZENF5RDSDiGasW5c1/HgRbIhzSAeVVS52N8cfoKF/5be6tywQ5U+Rcc7RYyqiRfZ6UF0mBZo9DNKiL5HArxKoH6kOcXqHOo/+FJO3GBfllaTLW0kHkdYB8NOhR7Rly/LuKWziQ05vJeW3PpmjPx1qKvpUXy7JgT8veUe5bBucgZTKVrWSzGBEN1rJgf88fWKUgFJVtbokB7Vc2DbU97diU5SGY+SgdGZoQZeLOVTiLl8Jqo2QPoYxtpKI9gDwEBHNt9TV3RGzlGcLGfsNgN8AwKRJkyqSZ8VF6vtNq4z019pysqt1TKuw8VKKiP69mvCPC7MfoBwJ+R9Tx+F0TcyFzypve3sS1KOuxAV8vJXEoPdJHmZrw7YjGpC8E93Kt0+Lnrnd9vlKXIWzfWZoQVpXrO4zHbVjX4nL48XX9VDXljpX7KHZiSzVr1Ny8Otbl8wxj0QsM0cXshkL0r+vPzu7wU6cxNKiVioUsvtGANGiSJ2EhTfSN6YdzGlKXyMkX+MTUE4csc8gDOrTjIs/lN7Ma2d71E7vFv84Bxm2qPFaoirJgTG2kv9fC+CvACYDWMNVReD/1/LqywGMli4fBWAlLx+lKa8LEslBUStJv6sRu8UK0fSBDu/fCx/lqbE/dqRWQMJpE5Lyfr30k6Kf3jyppO6TLJColcztiI+wGt9qokR19d79s/lngIReU5TwJyel8yS1NBVwtKEtG9IqORO96b2SdZJDwTHZytigyfOvw79OyWbOlSXAH585IWOrUWGSUpocCxcfm0Pm27DMY7YgRFeAom5PkEydWHJgVtWqnBpcjHOd5CAWL0L6F8/5B2ccDgBYKxIIejLX/q3NmHXpB+LvXaCNe2npXMN9JAdfKbRaVMwciKgvEfUXxwA+AOBVAPcCOIdXOwfAPfz4XgBnEVEvIhqLyPD8Alc9bSWiKdxL6Wzpmroha5BOjn1cQ00Q8RI2FfNq7pooUkOrGNA7EegG9daL5j4rOJZa+ZpsDlxysLqyikm7iiA42fXW4bnja/CrlFelXZLNtMibxOv89QsOemWsM6QvkfH3rxyLU9+V3VNDzvGjU0Nm6TKMX4crq+/ugjJs05itOetWpDDbIHTtt5fKxniKSK2UPD+xKpej34HIsP7tu18FkHwTYlvQwX2aUaAkBU01myUBwK6OcqofGUP6tmhTfcj4iGaM1APVSA57AniaiF4B8AKA+xhjDwK4AsDJRLQIwMn8NxhjcwHcCeA1AA8COJ8xJmbHLwG4AZGR+nUAD1RBlxUmV9JU3ICnTlEHl9EPAF54M8pP8+Crq7Tn5YlUl57ZlxZZVWSqH28Mb2lO8AQfycGkZ04FhDlcb30lFDXtiC98I93nSzlyTJHLQH6bQDXwWRSY6vgGweVp+831290XaWAOjtP3o7unJLV2yephuGlHoqZN0qGkvZhSSe04bYKBtZfKGOChIvPdhXCXRXLYa0CrM9llvW1yAhXbHBhjbwCYoCl/G8DU7BUAY+xyAJdrymcAsKfPrBFMBmn5vfrEOZjgCsBKtefxknVePb5Ib+eor5Pk+HerAHwmpo9p7CNq+65EgabnopLY22Bkd8Frb2cPo98q7g312AK3c4QPAzG9g96SvcU1qQLm5ysC+kwtyAxQ7PPhS2NeuKLt1XvQvSdBS3tn2bhQUOmVs+zKaqXlG5PVuqBNuH+3dZQxoLU5ZjImtY/vnH3vrEhrblJVCrp6Nxede4fXEz0vQlocKC8yr0ht3Ac23h/BTYvPCnmwQXLwgTwhmaQhcR82jVGSPsNN7zc/eLC2PDUhW1QAgNkhQIVJVeZC3ujyroJJ1fKzsybGxz4M2vRc7uaT0tINO7Tnd7QlDgwXf/AQbZ08vEFN0ifDpVbKbMuquW9hmF24ZlssjatQ6RXfRJNikF4s7RchmMPJ4/cEEAVrysb1PY0OAX4PZxt/ziIjsowiJbmVM/CReAAAEdlJREFUfCKl64ket5+DEB1s+muT/t1n9ScmfD/x3z259TasiHxokW/DxPzEALQxx3WW/aNVmFaEPjYdsXI1x5kov72pqpyWWuGzU8Y465i8wQ7dO0kf75PH3xWsaGLOPnYR07vJC3mcXPGxwzPnfRZONjdXgYzkUJYlh2xqdiAZE9MO2wvzvjcNvVuKXp56eRcURoM0/7ZdWWPrjR4rOdgM0ia1hjyATjTkqM/jSeCzQjalYfZZVJQ91Eqiim0yFHnun38jfy5/Ab88Q/y/8StLl48YWFncRSU7tFWL0UPctPpMQCZjfp52TAxcznVkgjrZ6lLb+0AwsFGDe+Mszd7mPipXV9wMAPRVXKDFN6F6BYkYhtuVLMpCpbdll3vTqzx8c8zQPpkNjgRd5XJEZzcLDj2QORhtDkmJKWBMHkymiUlIJOouUjr4rMJM2Tx9JAfZX9o0+R+69wC8d/+huPx0t8mnmqjvtEHaLl34Sg4/+WTG5AXAnFtIwCfJYo0FBy9JxLUxEeBnZzG180O+JaWJEuFFY8uflGeDJhvEN1bNmPL5xtQJWJYc5MSPonw/g5uwTUUmkMcec9y44dryIveu8snMWm/0QOagT58hv1hT8I66c5QO976yAgDw1KL1zro+qyOTnt8np43NPVWgtbmIP35+Cg4bmd35zpcWH+SzOfit5k3BYKadx5J+kmPfnEddAZceHgAGGcbmKYfu6WxH7IFsGhVCh/+B8XsaamRtU5VG6+45IJq0N2zzi//Q4WsnjXPWUXdgKxkkB3FsWoz5wPUkfvu5JJCvtyGoM8rK2v32BqAHMgeBjOQgHU/QbBEK+L0w11aDMnxWTSYVgJdaqcYulj4rWxPIR5WTwyvqaM1G7jI+duRIo+eUTItPqhQTnzjhIP3qr1L4MAcTvd/56KHxsSmPlovfHbVvlAvqUMtCIbtzXbbOYSOzGWpViHe8tc38vUwwRPULHOWRtlp9pjETKBS0Noda5UQ7++h9M+dlQ7ZJO1EsAOu3teGrt8/MnLvk1PG4R0lIWU/0OIO0abr0SWS39yB72oK8GOrhpmpa2fpIBR6CTi6YJh0f+KzWhQ+66dZ8XHMFfvLJifYKDlrkFfFN/zZZW6eazY90qGZiklV1LqnHdPaMI0bi2HHDsEd/8zhX71mnSrn9vKNx2P/+w0qDD96972C84mF0BoBLlWyxAipzkPcEFzaLV5ZtwhIer1FNFgDxCr7w/v203l7yszKpB4VK9Z+vrUmV92oq4Nxjx1ZMWyXoecwhNsCmy320CPsN74fp35qKPfpnDUl52hnYuxmbd3bgR2fqdeYyTBPGSJ4E7yMTzNGSgoH8vEYpfquRHGSJynRPv316CQDgtheW4iKNS6yc1bXdQ9/sA5+AMHXfCfVaNa1Hd8CHsQijqim7LhFZGYMOJ2tUULXawjJPsJcpKaTquTVlv0jaKBYIOzui7+O0Xz4Tn69KcuBs1ycxn2mh9dLSjdryR//r+IrpqhQ9jznw/5ncSp5jwuzjrG9XW4dXETpge119e6MG98HMS042RlADiVqpVpuRV2NzeGheshIyrbhFwM92g6pBZgg+nio+8LE5NDsM6McfZLdxdAV83vEJB+2BKfsNwUXT9LEoefHqd09BPw2j8WEOPhrP/9/eucdIVd1x/PPbXVDZ5SVvF9aloKhQFNhKERBRQfBR0giJhgpiRUnxQWO1Ym1sSk3xUeMzUdJC+kBta23EaiHWFJs+fECKyGtVrFHQVrE+QJKu6Okf9wx72Tvn3jsz987cYX+fZDJ3zz1z5ju/nbm/e37nnN9564P86zHyMbh3MMMpHNpzeP3WmQfHFFz7JpTUc4iYUOFv2tVz2LI7uJMcuGctpkmnG3MwjulKSa38vPuieOEMCM9Ls3jqsMhdt3rXdw0NIbSvYYgtKZRSbPR5nnnkLlxv43cObQnFzOIMfkdNHw270M0eN5jzRgc3fUqaOM6hV7euPHrFBIYcnf9CWiiuKbOF3Ixcf84I57m1W/8NeMkqo+hdn/8myd9z8A8219bUBLap9cpLCSvJIc+u8+AekHZRykB5sXS6nkOOjhegpJzD8P4Nedv3U9+1jo/2fxbax7j+nNLv7r5wzMwqlvNPTuYiF3V35up9jWrseTAW2xZjTn4c4mQodfcuvOewrV/vjBE6TIKkFqcl8Z6FfN/ybXiTo29DV/bsa2PZrJHOOjmcOx06HFjHxHs5orSHDbbnNMTZL9w1IO2i3pGdOU06Yc/Bew7MVkrot5X7wTSEJIVbffl4bpx5Qkl5k+KQu1lP6sJx5gnuKY6FUGzPYfHU4QeP8931JYn/h+xe9e3VycCsw9D0J+m9Z7oOae54b8aPa02AH1cM37VavK42mLI7DnfMdjv7T2zCvHyL2+DQHrzLOSycfOig88MLx/Pqj2aWNN5XLJ3POcTYJrQUuh9Zx7B+9dzpWKAF0Ny3nkVThiXzhiF8nsD0vDSImuXjTBLo+xxJjTnE0eAyX6642H2UkyROWo1yc/aJ7puJOKkhlpx9HDuWzQisVciH6+Lpmh4cZ6/mfJw4yN1z2G1TbR8/oCHvef8NR8eV2zly04kB3lx+HqcN6xtrinMadN6wUkrt1tXW8Ox1Z6TUemHkwkrlSvEblyhnFSfEV+pspRduOutgArR8+BW4Qg2Lpgzj+Tc+YNLw8LGhKP5y/dSSpglDZRIFhrH5B9NDV3NfcPIxPPLi20wMsZ2IxA6/uL5Trp5DbU0NBz43qdxkuMZ0/L9DV+/irBCHWm46nXNw3eRl4OYvNbJ24YgKc4Wd7VLrpT1oK/FHPaDHkYT9DON8HUY19mTDzdNK0gHQ1Kf0AeJKrOgOo8eR4dmETxvWl3/9+NzUdYf1HA588YUzm2tHVl8+nm3v5J9J1BHXlrb+36FrkD3ptTOlkB0lZcKVeC93l51b1n84sPzC0cwd3xS5mrjcuHoyt1zgLWQKW//x/FJvq5C0w0qbd30cXSmDLJjYXGkJsSmHQ3PuB1/rhZV2R+y6lmPi8L4sPD3/Hhc5cl9r194S/tBXnFBZpcm+woSZOWogxw9oCMQoc+HHpGYtZYHGXkdx69eD6ZAL5alrJrG/Lf1NRxZMHMqCieGrQHvZtSFfaY5OnXA4ccfs0ZG9pTeXn1cWLVefOZw9BaRxzyLGwJ59bdzw2OaDZT1CZk7FYc1Vk1jf+p4zxFXK3iyVoNM5h2P71HNsn2DmxVyyvdnjKr/aNWv49xOoNLU1wrolpztXxCbNxacOKcv7RDGnJRs6AK6b7l6bkDWG9+/OD2eNDKzkfmLT7kDdJ6+eVNJ7jWrsGZrAMmuhvygy4xxEZAZwD1AL/NQYs7yc799wRB07ls2IlVNfCfKbKyeE7mVxw4wR3L62NZFZYSMGdi+9kZhMHzmwbO+lpMO8Cc2Bso4RgnL1ulZfPj4yOrH+O2eETpYoF5lwDiJSCzwATAN2AS+JyBpjzLZy6ih0YYrSzqlDw8M8Cyd/idvXtjImItNm1piagdQYipt1S05nf1vhF9KjutYevAA3JbRiPA5hs7NyNDv2lCg3WblNPhV43RjzhjGmDXgUmFVhTUqCdKmt4bFFE1jlyHCqKMUwYmB3xjT1jq7YgRt8aTviXLA7I5noOQCNwNu+v3cB4yukRUmJlioaRH78W6fFnrqoVB9zWoYwp2UIO9/fR2Ov8oxfVRtZcQ75gnCBqeYicgVwBUBTU3DfWUVJirFNvRlbxB2pUl0M65d/NbOSnbDSLsA/HWMw8E7HSsaYFcaYFmNMS79+ye7CpSiKorSTFefwEnCciAwVka7ARcCaCmtSFEXptGQirGSMOSAiVwHr8KayrjTGbK2wLEVRlE5LJpwDgDHmaeDpSutQFEVRshNWUhRFUTKEOgdFURQlgDoHRVEUJYA6B0VRFCWAZGGLw2IQkb1Aq+N0E/BWRBM9gaik/VF1kmgDVG+abUDn1JuUliT0llOL6o2uM8IYE5290hhTlQ9gQ8i592O8fkWpdZJoQ/Wq3jT0JqilZL1l1qJ6I+qEXTv9j8M1rPRRjDpPJlAniTZA9abZBnROvUlpSUJvObWo3uLrHEI1h5U2GGNaCj2XRVRvuqjedFG96ZK03rjtVXPPYUWR57KI6k0X1ZsuqjddktYbq72q7TkoiqIo6VHNPQdFURQlJarGOYjIShF5T0S2+MpOFpF/iMgrIvKkiPSw5V1FZJUtf1lEzvC9Zpwtf11E7pUUdv1OUOt6EWkVkU32kcqelSIyRET+LCLbRWSriFxry48WkWdE5DX73Nv3mqXWhq0ico6vvBz2TVJv6jYuVK+I9LH194nI/R3aypx9I/Rm0b7TRGSjteNGETnT11YW7RumNz37xpnSlIUHcDowFtjiK3sJmGKPLwOW2ePFwCp73B/YCNTYv18EJuBtMPRHYGaGta4HWspg20HAWHvcHXgVOAm4HbjRlt8I3GaPTwJeBo4AhgI7gdoy2jdJvanbuAi99cAkYBFwf4e2smjfML1ZtO8Y4Bh7PArYnXH7hulNzb6p/cNS+hI0c+gF9xPax02GANvs8QPAN3z1nsXbp3oQsMNXfjHwUBa1pv2Pj9D+BDANb5HhIFs2CGi1x0uBpb766+wPqmz2TUJvpWwcpddX71J8F9us2telN+v2teUCfIB345Bp+3bUm7Z9qyas5GAL8DV7PIf23eReBmaJSJ2IDAXG2XONeLvO5dhly7KoNccq2138fhpd3I6ISDPencoLwABjzLsA9jnXZc2353cjFbBviXpzlM3GMfW6yKp9o8iyfS8E/mmM+R/VYV+/3hyp2LfancNlwGIR2YjXPWuz5Svx/rEbgLuBvwMHiLlXdUoUqhVgrjHmy8Bk+7gkTYEi0gD8DlhijPkkrGqeMhNSngoJ6IUy2rgAvc4m8pRlwb5hZNa+IjISuA24MleUp1pm7JtHL6Ro36p2DsaYHcaY6caYccAjeLFkjDEHjDHfNsacYoyZBfQCXsO7CA/2NZF3r+qMaMUYs9s+7wUexguNpYKIdMH7oq42xjxui/8jIoPs+UHAe7bcted32eybkN6y2bhAvS6yal8nWbWviAwGfg/MM8bstMWZta9Db6r2rWrnkBuZF5Ea4GbgQft3NxGpt8fTgAPGmG22q7ZXRL5qu1/z8OJ9mdNqw0x9bXkX4Hy80FQa2gT4GbDdGHOX79QaYL49nk+7rdYAF4nIETYUdhzwYrnsm5Tectm4CL15ybB9Xe1k0r4i0gt4Cm8c6m+5ylm1r0tv6vZNe7AlqQfe3fa7wGd4Hv6bwLV4I/2vAstpH/Btxhvc2Q78CTjW106LNeBO4P7ca7KmFW8GyEZgM7AVuAc7wyYFvZPwus+bgU32cS7QB2+A/DX7fLTvNd+zNmzFN6OjTPZNRG+5bFyk3jeB/wL77HfopIzbN6A3q/bFuzn71Fd3E9A/q/Z16U3bvrpCWlEURQlQ1WElRVEUJR3UOSiKoigB1DkoiqIoAdQ5KIqiKAHUOSiKoigB1DkoSgqIyCIRmVdA/WbxZfFVlEpTV2kBinK4ISJ1xpgHK61DUUpBnYOi5MEmRFuLlxBtDN7ixXnAicBdQAOwB7jUGPOuiKzHy4s1EVgjIt2BfcaYO0XkFLwV8d3wFlddZoz5UETG4eXW2g/8tXyfTlGi0bCSorgZAawwxozGS7m+GLgPmG28HFkrgVt99XsZY6YYY37SoZ1fAN+17bwC3GLLVwHXGGMmpPkhFKUYtOegKG7eNu25bH4F3IS32cozNjNyLV6alBy/7tiAiPTEcxrP2aKfA7/NU/5LYGbyH0FRikOdg6K46ZhbZi+wNeRO/9MC2pY87StKZtCwkqK4aRKRnCO4GHge6JcrE5EuNse+E2PMx8CHIjLZFl0CPGeM+Qj4WEQm2fK5yctXlOLRnoOiuNkOzBeRh/AyZd6Ht8XovTYsVIe3QdPWiHbmAw+KSDfgDWCBLV8ArBSR/bZdRckMmpVVUfJgZyv9wRgzqsJSFKUiaFhJURRFCaA9B0VRFCWA9hwURVGUAOocFEVRlADqHBRFUZQA6hwURVGUAOocFEVRlADqHBRFUZQA/wccKuN04h2XHQAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "sorted_data['inc'].plot()"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 10,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXd8ZGd97/9+plf1skXbqwuu642NCzbG4DjwowQTk1wDCcSEkAsJ/HIvkNzc/JL4hjT4QS4QCOFSktB7scHYYGxjbK+7d9fbm1ZadWlG08/Mc/845zlzRpqRRlqNNKN53q/XvqR95pzR0WjmfJ5vF1JKNBqNRqNx4lrpC9BoNBpN/aHFQaPRaDSz0OKg0Wg0mllocdBoNBrNLLQ4aDQajWYWWhw0Go1GMwstDhqNRqOZhRYHjUaj0cxCi4NGo9FoZuFZ6QtYLF1dXXLz5s0rfRkajUbTUDz55JOjUsru+Y5rWHHYvHkz+/btW+nL0Gg0moZCCHGqmuO0W0mj0Wg0s9DioNFoNJpZaHHQaDQazSy0OGg0Go1mFvOKgxBigxDiZ0KIg0KI/UKI91rrfymEOCuEeMb6d5vjnA8KIY4KIQ4JIV7lWL9SCPG89djHhRDCWvcLIb5qrT8mhNi89L+qRqPRaKqlGsvBAN4vpbwAuBp4txDiQuuxj0opL7P+/QjAeuwO4CLgVuCTQgi3dfyngLuAHda/W631twMTUsrtwEeBvzv/X02j0Wg0i2VecZBSDkopn7K+jwMHgfVznPJa4CtSyoyU8gRwFNgrhFgLtEgpH5Xm+LkvAq9znPMF6/tvADcrq0Kj0Wg0y8+CYg6Wu+dy4DFr6Y+EEM8JIT4nhGi31tYDZxyn9Vtr663vZ66XnCOlNIApoHMh16bRaJaG7z07QCydW+nL0KwwVYuDECICfBP4YyllDNNFtA24DBgE/kkdWuZ0Ocf6XOfMvIa7hBD7hBD7RkZGqr10jUZTJcPxNO/58tP84NnBlb4UzQpTlTgIIbyYwvAfUspvAUgph6SUeSllAfhXYK91eD+wwXF6HzBgrfeVWS85RwjhAVqB8ZnXIaX8jJRyj5RyT3f3vNXfGo1mgWRyBQBSufwKX4lmpakmW0kA/wYclFJ+xLG+1nHY64EXrO+/B9xhZSBtwQw8Py6lHATiQoirred8C/Bdxzlvtb5/I/CAFZfQaDTLSC5vikPWKKzwlWhWmmp6K10L3Ak8L4R4xlr7EPBmIcRlmO6fk8A7AaSU+4UQXwMOYGY6vVtKqbYh7wI+DwSBe6x/YIrPl4QQRzEthjvO79fSaDSLIZc392RaHDTzioOU8mHKxwR+NMc5dwN3l1nfB1xcZj0N3D7ftWg0mtpiWw557VZqdnSFtEajsdFuJY1Ci4NGo7ExCtqtpDHR4qDRaGxyhnIraXFodrQ4aDQam5xlOWS05dD0aHHQaDQ2ynLQ4qDR4tCAnJ1Mce8LuoJVs/QYBR2Q1phocWhAvvzYaf7wP54iX9B1gpqlJavrHDQWWhwakETWoCDNrxrNUmLoVFaNhRaHBiRt9b2Jp7U4aJaWYhGcFodmR4tDA5K2mqNNa3HQLDG6fYZGocWhAUllleWge+5rlhbtVtIotDg0ICntVtLUCNty0G6lpkeLQwNixxwyWhw0S0tOp7JqLLQ4NCDFgLR2K2mWlpyhK6Q1JlocGpBKbqV0Ls8zZyZX4pI0q4RiEZxu2d3saHFoQCplK/3guUFe94lH+MTPjq7EZWlWAVmdyqqx0OLQgKQquJXGExkA/uHHh/j6vjPLfl2axsfQqawaCy0ODUilIrjpjLm+tjXAL46MLvt1aRofVQRXkMW0Vk1zosWhAamUrZTIGIR9blqDXvsYjWYhqFRW0EHpZkeLQ4Nh5Av2B3imWymRMQj7Pfi9bv3B1iyKnMNa0K6l5kaLQ4ORdnxgpzMz3UoGEb+HgMelLQfNonC6knRQurnR4tBgqNYZMDvmoCyHgNdNRouDZhE43UracmhutDg0GMoiCHhdZcQhT9jvJuB12emuGs1CcLqVtGuyudHi0GAoceiO+mfVOdhuJa+btC5i0iwCHXPQKLQ4NBiqxqE74iebL5TEFhJZy63kceuYg2ZRGI7pgjrm0NxocWgwlLuoO+oHSuMOxZhD0a30wW89z8d+emT5L1TTkDitBW05NDeelb4AzcJQlkNPNACYriQlFMqtBJCx3EqPnRhjQ3toBa5U04gYBYnP7SKbL2hxaHK05dBgOGMOUKx1MPIF0rkCYZ9KZS0gpSSVzdtCodHMRy5fIOx3A5DN6/dNM6PFocGYLQ6mWylhpbiG/W78XvPDnTEKJLN5nbmkqZpcXhK2rE9tOTQ32q3UYKg6h56Z4mAVxEX8HoQwj8nkCpbloD/kmurI5QsO16R+3zQzWhwajEpuJSUOYb+HvDQzTqazBtl8QRfEaarGyBdoCfgAbTk0O1ocGoyU5SJyBqSdXyN+j/2hnkhkAb0D1FRPiVtJp7I2NfPGHIQQG4QQPxNCHBRC7BdCvNda7xBC3CeEOGJ9bXec80EhxFEhxCEhxKsc61cKIZ63Hvu4EEJY634hxFet9ceEEJuX/ldtTNK5PO/69yc5OBiz/w/QGTF3d0W3koo5mEVwABNJJQ7actBURy5vJjWAthyanWoC0gbwfinlBcDVwLuFEBcCHwDul1LuAO63/o/12B3ARcCtwCeFEG7ruT4F3AXssP7daq2/HZiQUm4HPgr83RL8bquCR46Ocs8L53jkqDmfIZ3LE/C68LpddIR9nIulgaLloNpnAIxbloMOSGuqJZcvEPIXExo0zcu84iClHJRSPmV9HwcOAuuB1wJfsA77AvA66/vXAl+RUmaklCeAo8BeIcRaoEVK+aiUUgJfnHGOeq5vADcrq6LZ+fmhEQAmk2ZsIZXL25bBxo4Qp8YSQGlA2u8xHx9PaMtBszCMvLQD0tpyaG4WlMpquXsuBx4DeqWUg2AKCNBjHbYecM6o7LfW1lvfz1wvOUdKaQBTQOdCrm01IqXkZ4eGARhPKisgT9ASh82dIU6OJgGzdQZgV0gDTFiCkstL8o62CBpNJXKFAgGvGyG0ODQ7VYuDECICfBP4YyllbK5Dy6zJOdbnOmfmNdwlhNgnhNg3MjIy3yU3PMdGEvRPpACYtMQhlSsULYfOMINTKTJGviQgbcccLMsBtPWgqY5cXuJ1C7tKWtO8VCUOQggvpjD8h5TyW9bykOUqwvo6bK33Axscp/cBA9Z6X5n1knOEEB6gFRifeR1Sys9IKfdIKfd0d3dXc+kNzc8tq2F9W9B2EaWyRbfS5s4QBQn9EykSGQO3S+D3uIoxh6RDHHTcQTMPhYJpYXrdLnwel7YcmpxqspUE8G/AQSnlRxwPfQ94q/X9W4HvOtbvsDKQtmAGnh+3XE9xIcTV1nO+ZcY56rneCDxgxSWamqfPTLKhI8hF61rsmEPGyBO0bv6bOs2eSafGEuYsB58bIYQdc3BaDrqFt2Y+cgVTDLxuF36PSwekm5xq6hyuBe4EnhdCPGOtfQj4MPA1IcTbgdPA7QBSyv1CiK8BBzAznd4tpVR3pncBnweCwD3WPzDF50tCiKOYFsMd5/l7rQqGY2nWtQbpCPt45swkUGo5bOoMA3BqLFnSdE89Pp7QloOmegxrCpztVtLi0NTMKw5SyocpHxMAuLnCOXcDd5dZ3wdcXGY9jSUumiJDsQyXbWijLeRjMplDSknayNMa9ALQGfYR8Xs4NZa023UDjoC0M+agP+iauVGDfjwuF36vW8ccmhxdIV2nSCkZiqXpbfHTEfaSzZtN9JyWgxCCjR0hTo4lyBeKla1Ft1LOfj49/EczH2p+tNfjsiwH/Z5pZnRX1jolljbIGAV6WwK0hcxq6PFElrQjWwlgc1eI05bloNxKXrfAJUrbH2jLQTMfynLwuoQOSGu0ONQrw1blc3fUT7slDpPJnFnn4Cv+2TZ2hDkzkWQolrH78AshSgQEtOWgmR8Vc/CobCXtVmpqtDjUKUOxDAC9LQE6wmaMYTyZNSukPcUb/20vWYNLCM5Opmy3EjBLHLTloJkPJQY6IK0BLQ51y3DctBycbqWJRNayHIo3/kv62vj0nVfic7vojvjt9YCn9E+ri+A082E4Ulm1W0mjA9J1irIceqJ++0M6FEtTkLOtght39XDf+26g0ykO1jGtQS9TqZxuvqeZl5yhUllNcdDWZnOjxaFOGYqlifg9Vq8kiRDw2AmzaLyvPTjreFXzoFCjQjvCPqZSOW05aOZFFcF53DogrdFupbplOJ6mp8W0BNwuQVvQy0NHzH5SV2+dvyeh33IrtYXMeIW2HDSVkFLypV+dYiplpj773C78bm05NDtaHOqU4ViGXmvaG0B7yEcuL9naFaa3JTDHmSaqEK7Dildoy0FTiWMj0/yP77zA9581W515VCqrzlZqarQ41ClDcbMATqEsgKu3VdfJ3I45WOfp9hmaSiSz5sZhcNJMgvB6rN5KOv25qdHiUIeY1dGZEguhI2xaANW4lAA73TXs8+DzuHTjPU1FlMtRTRX0ulxEAh4S2Ty6/2XzosWhDplK5cgaBbqjTstBiUNHVc+h3Eohn5uAx6UtB01FVIHk4JQ5O8TrEbQGveQL0p4Tomk+dLZSHXJm3PyQrm0tZiXd9pI1tAS89ETnjzdA0a0U9Lnxe906uKipiHpvKAvC43LZzR2nUjmiAe+KXZtm5dDiUIc8ecpMWb1sY5u99vLdvbx8d2/Vz6HEIeRza/+xZk5mtlbxuUvFoa99Ja5Ks9Jot1Id8sTJCda1BljfNrueoVpUKmvQZ44N1ZaDphIzxcHjFrQEiuKgaU60ONQZUkqeODnOns3VxRYqoYrgQl7TctCN9zSVmLlx8LgFLZblEEvpmEOzosWhzjgznmI4nuGqLecnDs6AtB75qJmLudxKMW05NC1aHOqMx0+a8YarNp+fo1elsgZ9bsutpC0HTXlmWw4uuz5Gu5WaFx2QrjP2nRynJeBhZ0/0vJ6nGJD24Pe4iKe1e0BTnpnJCl63wOty4xJaHJoZbTnUGcdGptm9tgWXq9LY7uooqXPQloNmDtIzLAevy4XLZcYdtDg0L1oc6oyzEyn6ziNLSXHx+lYuXNvCho6QFZDWMQdNedK5PMLai7hdwt6YtGpxaGq0ONQRuXyBc7E068u05F4oO3uj/Oi919Ma9OL3aMtBU5l0Lk9XxI8QZtM9hRaH5kbHHOqIc1PmMJ/zqW8oR8Crs5U0lckYBcI+N4WQr2SGgxaH5kZbDnXE2UmzbcZSWA5O/F63rnPQVCSdyxPwuumK+PG4i5ZDS9CrU1mbGC0OdcTZCUscltpysOocdIdNTTnSuQJ+j4vOiA+vu3hL0JZDc6PFoY5QlsO6JRYHv9eNlJDLa3HQzCZj5PF73axpCRDyFeeTtwRMcdCbiuZExxzqiIHJFF0Rv12jsFSoPktpI4/Po/cDGpPHT4yzZ1M76VyBlqCXP7llJ2OJrP14a9CLUZCkcnlCPn2raDb0naKOODuZYn1bdS25F4Lqs6RnOjQvX3z0JN995qz9/2fPTPKmTz/Ko8fHSOfy+D0uNnSEuGxDsROwszOrpvnQ4lBHnJ1ILXkwGhyWgw5KNy3/8avT/PuvTtn/PzmWAGB0OkPWKJS1VrU4NDdaHOoEKaVlOdROHHQ6a/OSzBl2wgPA4JQ5EjSWypnZSmXcjbY4JLU4NCNaHOqE0eksGaNQE3FQu0JdCNe8pLJmgWUub82LVuKQNkhry0FTBi0OdcLp8SQA69tDS/7cUb8ZTBydzs5zpGa1ksoaFCQMxUxRGLAy4+Jpg4wVc5iJFofmZl5xEEJ8TggxLIR4wbH2l0KIs0KIZ6x/tzke+6AQ4qgQ4pAQ4lWO9SuFEM9bj31cCLObixDCL4T4qrX+mBBi89L+io2BGg16aV/rkj/3ZRvb8HlcPHhoZMmfu96478BQSZWvxnRZpqx4k3ItKbdSPJ3TloOmLNVYDp8Hbi2z/lEp5WXWvx8BCCEuBO4ALrLO+aQQQr3rPgXcBeyw/qnnfDswIaXcDnwU+LtF/i4Nza+Oj7O1O0xPy9JnK4V8Hq7d1slPDw6t6pz1w0Nxfv+L+/j+swMrfSl1RTZfoGD92VUtjRKHiWSWfEHaXXydRAMevG6hLc4mZV5xkFL+Ahiv8vleC3xFSpmRUp4AjgJ7hRBrgRYp5aPSvDt9EXid45wvWN9/A7hZWRXNgpEv8PiJca7e2lmzn3HzBb2cHk9ydHi6Zj9jpVGukiOr+HdcDKlsMdY0MJkiY+QZnc4AMBI3v/o9sy0Hl0vQEw0wHE8vz4Vq6orziTn8kRDiOcvtpMaWrQfOOI7pt9bWW9/PXC85R0ppAFNA2bukEOIuIcQ+IcS+kZHV4yLZPxBjOmPUWBx6APjpweGa/YyVZti60R0b0eLgJOVIYT47mWI4lrH/r6yCcpYDQHfUX3K8pnlYrDh8CtgGXAYMAv9krZfb8cs51uc6Z/ailJ+RUu6RUu7p7u5e2BXXMb86PgbA1ec5N3ou1rYG2b0myi+PjdbsZ6w0I1ocypJ0WA79Eynbwor4PUXLoUJVfm+L3w5ia5qLRYmDlHJISpmXUhaAfwX2Wg/1Axsch/YBA9Z6X5n1knOEEB6glerdWKuCJ07WLt7gpDvqX9XjQoetm9jpsaSdsqkpupV8HhcDkyk73rCzN8J0xnw/VGrZ0tsSsC2ySvzguQGeOj2xhFesqQcWJQ5WDEHxekBlMn0PuMPKQNqCGXh+XEo5CMSFEFdb8YS3AN91nPNW6/s3Ag/I1Rw1LcO5WJpNHUufwjoTc+jP6r1pqpuYUZB2arCm6Fba2hXm7GSKgSnTctjZW5xTXi6VFaAn6mfKKpSrxN0/PMhnHjy+hFesqQeqSWX9MvAosEsI0S+EeDvw91Za6nPATcCfAEgp9wNfAw4A9wLvllKqd9W7gM9iBqmPAfdY6/8GdAohjgLvAz6wVL9coxBLGbRYaYO1xBz6s3oL4YbjGbum45gOStsoy2FbT4R0rsCBgRgtAQ+9Dku1kuWgrNm54g6xVM4OcCtOj2lxbnTmbbUopXxzmeV/m+P4u4G7y6zvAy4us54Gbp/vOlYzsXTOzimvJX6Pe1U33xuOp9m7pYP7Xxzm2EhipS+nblCWw46eCAA/2T/E1u4w0UDx41+ufQZgC8hwPM3GztnWrZEvkMjmGXGIw1OnJ3jDJ3/Jj//4Bnatic46R9MY6D68C+TgYIxdvVF7CPv5IqUklsrRElgGcVjFloOUkuFYhlsvWsNzZ/0c10FpG2U53HJhL1LCc/2TXL+jm4i/+PGfKyANMFTBclAxrFFHXEIFvIdiaS0ODYxun7EAHj8xzq9/7CF+8Pzggs4z8gVOjZXfyU5nzLYGLcHa67Tf41q1lkMsbZAxCnRH/WzrDuuMJQfKcugM+/mTW3byf353L7933ZZSy6FCKmtP1LQcKmUsKXFIZPMkrOC2WnNmSWkaDy0OC+CLj54E4NFjY2Uf//bT/fw///vhWes/euEcr/jIg4xNz959xawP0nJYDgHv6g1Iq5TMnmiAjR0h+h0dSJsddZMOzrAOoo73XKBMERxAe8iL1y0qZizF0sXWGiruELfWktnVmxnXDGhxqJLheJp7XzgHwL6TxUzbw0NxfvicaUkcGIjxXP/UrDTK4ViaXF7arQucqAHuyxGQ9ntcZiuFwupLBlNVvD1RP+0hn+4H5EBlGgV8pR93p7VaKSAthFUlXcFyiKXKiUPRmtA0LlocquSrj5/BKEh+84o+jgxPM2GNU/zrHxzgQ99+HjAHtUNpRSoUd27lMj5scViOmINHte5efdaDbTm0+GkNeckYBT3cyCKZNXC7BD536cfdaTlUSmUF8zUdssR3MpnlpweG7MecloP6GyhxSGnLoaHR4lAlj58c56J1Lbxpj1nL9+SpCcamM/zy2Jgd8FPB3tSMHVPC+pAMlelRo9xKy5OtpIb+rL6bphLe7miAtqAPgMklHFKTNQr8008ONaSrJJUtEPS6mdmyrDTmUHlueW80YL++X9/Xzzu+uM92HcVSxddjxGrFoQQjkVl977NmQotDlYxOZ1nbGuDSDW143YInTo3z4/1D5AuSbL5AviBty2FmIC5VjeWwDAHp4tCf1Wc5DMfT+D0uWgIe2kKm0E6mlq6b6PNnJ/nnB47arU4aiVQuT9A3++bvFIe5LAdnCw31mqrK6rksh0YUUk0RncpaJaPTGS7tayXgdfOS9a3ct3+IiOPDlTHythtj5odC7aCcQb0fPjdIS9Bjf7iWx620emdJD8cz9LT4EULQZllhS2k5qCyvrFE/8ZrTY0n2nRrnDVf0zXlcKmvMCkaD6Wb0e1xIyZyp2R1hP7G0YdY0WO9l9TWWNhAC2oLeMgHp1fc+aya05VAFhYJkPJGlM2K6K95+3VaGYmme65+iI2yupXMF0taOfKZbKZUzxWLE4Vb6yH2H+MTPjtpmuXMXVyv83tU7S3o8kaUjbObkt4ZqIA5WkkE99Wz60q9O8v6vPztvgkEqly8rDmDGHfwV0liLx5jvzemMMcsqiKVyRPweeqKBMpaDFodGRlsOVTCZypEvSLoi5s3nNy5Zy3Xbu/jJgXNMpXL8zQ8Pks7lydiWQ4WAtMNymEoZGAVJLJ0j7HPjcddep+2A9CqsdYilDTtu0xYyBXtqCd1KarpcPU2ZG5hMI6U5zCfgqhwzSGbLu5UAWgIe5umrZ1vI8bRh1zIULQezgLM76p+VraTdSo2NthyqQL3pOy1xAHN3evueDbZgpHP5ouUwM1spUxpzUFXRg5NpJpO5ZUljhWKh02oMSMdTOVqsm1gt3EpKFOrJclAN9OazBNNzWg6eOeMNUJxBPp0x7OSKouVg9gXrivgcloN2K60GtDhUgRKHLsut5ETdcNO5gm05zHQrJS230uh0hoIVuM7mzX/HR6eXJd4ARcshvQoth6lUsT9VyOfG6xZMrHJxGJw03ZTziX0ymydUwXKIBrxzZiqpY6DUraRqGOJpU5SV5SClLB6T0ZZDI6PdSlWgpmV1OywHhepJky4JSJe3HIyCZDyZxcgXfcSHzsW5eF1rTa571rWu0lRWKU33nLLAhBC0Bn1L61ayRCGbr4+AtJEv2IV/87m6Urk8gQricOOu7nnnNRTdSjn7hp+0s5UM+tqDdEf9pHMFRqezGFYMRFsOjY0WhyoYK+NWUgTs3XjeNu9n+lqT2TzRgId42mA4lsHtyAxJZvPLksYKqzcgnc4VyOVliQXWHvJW7VaaSGS57+AQt1/ZN6sWQFFvlsNwPIOKQ8/rVspWdiu94/qt8/4s1aAvnjbsFFY7lTWVo2Vti+1edfa00uLQ2Gi3UhWMTps39LYysQHbj58rVuTOcitlDTZ3hgEzH39ma4flcisF7Arp1fWhtdOBHSLbtgBx+Nt7DvLfvvEcJ+eYQVBvAenBqWIrlvkSDJK5ym6lamhxZCspUVA3ftNi87CuLQiYljCYcQodkG5stDhUwWg8S0fYVzYXXPlr07l8sQiuTPuMzV1KHDKzxWGZAtJ+R3xkNaFeT2eVeWvQx2QV/ZXOjCf51lNnAeacHpets1TWgcliWvR8Yp+aw3KoBuVWiqUc2UpZg0JBMp0xiAa8bLAmGR4YiAHQ2xrQlkODo8WhCsYSGdtsnknAGXMo0z4jaxQwCtIeAzriEAe1m2tZhhoHcKayrq4Pbbn+VG0hL1PJ+WMOn/z5MQrWVNozc4iDct1k60QcnJZDJWvm1FiCkXiGjFGomMpaDUGvG5cw37vKlZXM5IlnDKQ0379rWgJ43YL9g1MArGkxxWE1NnlsFrQ4VMHIdLZsphIU3UrxtPlBgVJxUKZ1R9hHS8DDcCxt38zUIJRlsxw8qzPmUHQrOcQh6J3XcsjlC3z76X7eeGUfXrfgzMT8bqVcnVRID045LYfyf893fGEfH/jmc8Dsdt0LQQhBxO8pmemQyBp2ympL0IvbJVjXFuTwOTPm0GMNCZqZ1q1pHLQ4VMHY9ByWg7Ubn3L4t51uJWVah3xuelsCDE4VYw6717QAyxdzWLXikFIzMUpjDslsfk6Xy4GBGOlcgZft7GF9W5D+8cozIOotID04Ob84jExneNTqBXU+lgOY6axOayWZyTted/P9u6E9ZFtWa6zxotq11LhocZgHKSWj05k5LAfzQ+fcpTpbFSvLIeQ3g3ZKHKIBDxs6zCDeclkOHrcLj0usut5KynIoiTnYVdKVrYenTk8AcMWmNjZ0hOa2HPLma1ZPAWk1wrOcAEopmU4bFQf9LJRowMO5qVLLYWYigHo/Q3H2tA5KNy5aHOYhkTUDzeXSWKG4G3dmxiSzZSwHr5t1bUEGJlPErIKt9W1KHJYvo9jvca06y0FZbc75BCqzbGqOjKUnT02wtjXA2tagKQ5zBaTrzHIYmEqzycqAKydYGSvWpThfyyHi99j1EB6XIJnNz4r19LWbcTUhoDtqfl502+7GRYvDPIzZ1dHlxcHlEvg8rpKCq7Li4HOzvi3AWCLLUDxNa9DLNds6eeWFvVy0TEVwYBbtrcZU1qDXjc/RBqLYtruyODx9epIrNrYDpktkIpmzUzVnkq2jgHQuX2B0OsPmTvNmXE7sna20gfNKZQXTclBi0x31k8gYs0bcbrSSLiJ+D2GrNkI1ndQ0Hloc5mBsOsPf/PAgABvagxWPC3rdtvvCPcNtM9OtBPDiYJyWgJeeaIDPvGXPsgz6UQQ8rlXXeM/s71Nqfc038OfcVJqzkymu2GSJg+USuf/gEO/4whOzWj/UUyprwsoSUq6bctln0+nS65+vRcZ8RBxWWU9LgETWsDdOHZbLVaWztgS8hC0x0pZD46LFYQ7+7Nsv8OChET502272bumoeFzA67L7+LQFvRUtByUOY4nssgqCE7/XbTcIXC3E0rlZr6eyHCYqpLPa8YaNbYBpOQD8z+/t56cHhzk8FC85vuhWWvlsJZUBpH7ncpaDsoDWtpoCcr4xB1UlDeac7mQmz3A8Q8jnth9TG6howEPIZ67pgHTjosVhDk6PJ7lhZxd33bA2aKAVAAAgAElEQVStYlsFMHdlaofaHvaVikPG6VYqWh8rJg4e16qrc5hK5WZlfLVbczbUrO+ZHBs2Uy4vWGtmjKldr/o7js84L1NHFdIqVVq9h8pdk7IcbtzVDUB7qHxCRbU4M8F6W/wksoY5YCladLd2hH2EfG5LHEwx0gHpxkX3VpqDeCZHNBCd97iAx83ApJnm1x7yMjjpSPlTbiWfh4jfgxAgZXEgzXKzGgPSsXSOnmigZC3sM6ecjVUQh/FklojfY7tb2kOmK0QFcmeKQz3FHNTmIxrw4hKVYg7m++63927iN6/osyv0F4uyDlwCOsNmk71zU6mS110IwfaeCD3RACG/5VbSlkPDosVhDuJpo6oJbQGvy3Y3tId8JHN5pJQIIewPR8hnBkx7on6GYpkVdSutuoB0ymB7d+nfSQhBV6Q4gGYm5uQ4X8nx/+WaTfRGA/zVDw7MFoc6ijmomFbI58bvKf/3VG6l1qCXjZ3nn/CgWmiE/R5bKE6NJblqc6m79ZO/cwU+j8t2K6W05dCwaHGogOpLX404+B3+3I6wDymhfyLFvlPjpLJ5XKKY8rquLchQLLNstQ2zrtXjsvvtrxbKxRwAOiO+WTd5xXgia7ueFB/89QuQUvLhe1+saDnUgziomEPQ2nCUjTmkVXrv0nzEVZpwxO+xrYLheMZOWVWodNa8ldmkA9KNi445VCCVy5MvSCL++W/izkwQNaLyXx86zp989VlOjCYI+Tx2zEIFpVcu5uBeVW6lQsGcqldObDvDPsamy4vDRDJLZ3i2H14IYZ5XSRzqoH1GylHY5ve4ysYc1AYg7F8acVDWQtjvIexzBKdbyqd4u12CoNetYw4NjBaHCqgPV1VuJUd+fbsVS1DdKZ86PVFSgLR+pcXB61pVbqVE1qAgy7cg6Qj7GZvOkDHyXPvhB/j+swP2Y+PT2YpB2o6wb1Ygu57cSspyCHjd1t+zfLaS3+Mqqf04H9TnwFnDAMyK9TgJ+dw6W6mB0eJQgfgCzHKn5aBcFQcHTXEYnErbOd8A66zUwpUSh4DHvarqHOxCrDJV5l0RH6OJLP0TKc5Opnji5Lj92HgyS2eFligdc1gO9WB1pR1upUoxh7jVSnupKBEHx/u5J1recgAI+bU4NDJaHCows/pzLlRnViimDDqzNIIOM/yqLR3s6ImwpfP8skcWy2qzHGJlZjkoOiM+skbBHkCjhvkkswbpXGFOy6GuYw4Ot5LPXb6ocbrKeFm1KLeSGXMoPu/MmIOTsM+j50g3MPOKgxDic0KIYSHEC461DiHEfUKII9bXdsdjHxRCHBVCHBJCvMqxfqUQ4nnrsY8LywkvhPALIb5qrT8mhNi8tL/i4liQW8myHHweV8muSuFcu2hdK/e972Urm8q6iiwHlZVTzrfeGTZvXM+emQTM+QZQrGEoF3OA+heHpDNbyesqm14bT+dKCtfOF2e2UrWWQ9Dn1i27G5hqLIfPA7fOWPsAcL+Ucgdwv/V/hBAXAncAF1nnfFIIod5JnwLuAnZY/9Rzvh2YkFJuBz4K/N1if5mlpOhWqj4gHfC4SuILu3rNGonzbXq2lKy2gLQSunLtIVRbh2csceifSJHLF5hIFAsWy9EZ9jGdMUosrEy+fiqk05bl4Pe4Kor9dMZYUnFosbOV3Lbl4HGJOYvruiN+zk5UboOuqW/mFQcp5S+A8RnLrwW+YH3/BeB1jvWvSCkzUsoTwFFgrxBiLdAipXxUSimBL844Rz3XN4CblVWxkiwmIO33ukuE4NaL1wDn3/RsKQlYO838KpnQpfzvaq6Gky7Lcnj+rDmdLF+Q9E+kGEtYPYHC5YW/wzpPiYiUsq6K4FI5c+ynEKJyzGGJ3Up+j4uOsI+1bUHbcuiO+suOzlVcuqGN46OJOTvjauqXxcYceqWUgwDW1x5rfT1wxnFcv7W23vp+5nrJOVJKA5gCOsv9UCHEXUKIfUKIfSMjI4u89OpYSEBa1TkEvC5CXvN4t0vwigt6AUpS/1YaNSq0HtpALAVqNKvfO/utrALOyWze3kWfHEvY/ZaUCMxEiYYSEWUtuIT5ukm5ssKayuXtTUilOod42rBdQUuBEIJ733s9b3vpZrvAbS6XEsClfWbfqufOTi7ZdWiWj6UOSJfbRsg51uc6Z/ailJ+RUu6RUu7p7u5e5CVWRzxtIER1N/aiW6loOaxvC7JrTRSvW9hFQ/VAcRrc6vAFp5VbqYzl4KyAvnqrWcl7ajRh1z50VAxImzc9FXdQ1oKKaxgrbHWlsgW7kV6lOofpjEF0Cd1KYHZjDVit0b1uQfccaawAL+kzK7NVzEfTWCxWHIYsVxHW12FrvR/Y4DiuDxiw1vvKrJecI4TwAK3MdmMtO/G06bOdy2xWqGylgNdtu5A2d4XxeVx85E2X8ZZrNtfyUheE2mE3QtwhX5B2ELkStlupjOUQ8LrtG+SlfW2EfW5OjiWZSGZxu0TFIUtKVGxxsF4r9VwrHZRO5/L271suhiSlNMWhhuNnO8I++uZoYw9mBtnW7rAd89E0FosVh+8Bb7W+fyvwXcf6HVYG0hbMwPPjluspLoS42oonvGXGOeq53gg8IFfabsdsyVDtbGfnLk59rwaxvObSdezsnb9533KhdtiNMCr0B88N8IqPPGjPDSiHujH6K7Skds4a2NQZ5tRYgvFEjvaQr2Kn3c4K4qAsh5Wukk7l8rZrZ2Zq8n0Hhtg/EDOr+5fQrTSTL/7er/Gem3fMe9xlG9p45szUirviNAunmlTWLwOPAruEEP1CiLcDHwZuEUIcAW6x/o+Ucj/wNeAAcC/wbimleue+C/gsZpD6GHCPtf5vQKcQ4ijwPqzMp+Xgl0dHufbDD9jxBScLCejZbiWvG5dL8Nevu5i3XLNpSa91qWgky+HcVJpcXjJaoQUGzG05QPFG39ceZHNXiFNjScYTmYpprGDueF2isjisdFA6mTVKNiTObKX3fe0Z/vw7Ztb5UmYrzWTXmmiJ264Sl21oY3Q6w+k5RrBq6pN53z1SyjdXeOjmCsffDdxdZn0fcHGZ9TRw+3zXUQsODcU5O5ni8NA0V25qL3ksns4tQBxcJV/vvLo+hQGKAelGqHVQNQwzR146yeTyCAE+dwVxsMa7bugIsa07wo/3D1GQkjWtlf3lLitFc8yOOZgCFKkTcUjlCnbRn8/jstNss0aBeNqw3ThLma20WF66rRO3S/Caf36Yv3rtxbzu8vXzn6SpC5q6QloFM0+OzvZrm5ZDdW4l5aqp5NqoJxopIK3SictZdoq0UcDvcVV0EXVFfPjcLrojfn7rqg24heDkWHLeXW9H2Me4ZbFkbMvB/PvmVtjqSmfzBB0xB5VBNZkqtbDqQRy290T57ruvZU1rgI/ff2SlL0ezAJpcHMwb5ElH0DOdy2PkCwtyK/kdJn69o1xg6QawHJQ4xFKVWzCYwdnKovy2l27hH26/BJdL0Nce4q0vNa26+cTBOQtCuZVUh96VDkirOgdwin1hVj1BNR2Fl4OL17dy9dbOioOXNPXJym8tVhCVI3/Cshwmk1ne8MlfcumGtkW6lerfclCukbl24/XCdMa8xrndSoWyaayKXWui7FpTTAh4903b+e4zA/MmCfS0+Hn6tOmeKYqDVSNSD+Lgmy0Oao65xyUwCrKmMYeF0hH2MZXKkcsX8FZwAWrqi6b+Kym/+8mxBEa+wB/959McH03ws0PDC3MrOeoc6p0uK3un0oS0esKOOaTmcivlyxbAVaIt5OORD7x83vTinqif4XjarI6eUeewEi00skaBP/36sxwdjltuJStbyRKHrFGwi/uu29EF1IdbSaESANQ1auqf+nn3rADK735yNMl3nhng4aOjvGxnNw8eNquvF56tVP9aq9wpc2UA1QvTyq00x+S6dC6/YFGuZufaEw2QzhWIZ4zZqawrYDnsH5ji60/2s70nYlkOxZgDmO/lSevG+56bd7B3S8e8dQjLiSosHJvOzjkDQlM/1P/drIYov/t0xuBLvzrF+rYgd7++mFBVfUC69INaz3jcLtpDXrs1RD1TjDnMYTnkCjURZdWKeiSecbiVrGylFQhI77eGR52dTGEUZDHm4EhNnrTcSjt7o/zhjdsrBulXAtXKpNLYVk390eTiUMzYefbMJK+8qJe+9hAbOswdV0uVloMaur52jvTIeqIz4mc0Xv8f0ngVqazpXL4mWWKqb9BwLDPLrbQSMQclDmomRWBmQDpnxhy8blG2bfxKo9xKOijdODS1WymdyxMNeOwd6qsuMruoXr2lkzPj/QtyKz34pzeu2HS3hdIV8TWE5TBdTbaSUajJ665mIw/H03Yqa8SukF5+cTgwYHaWPW1l1jkb74EpWJPJLG1zVH6vJHZLkgaIdWlMmtxyKLCtO4LHJegI+7hqs9mc7eqtZlPY1uD8FaCKzogfT4NkYXRG/HUTc7j3hXNl22MY+YI9KCaWzjE6neHuHx6Y5dLJ5PI1SSHujphWYDm30nIHpI18gRetaXb91nyEYiqrKmrMM5HM2jPM6w1TtLRbqZFojLtZjcgYZivnKza284bL1+O2muy95tJ1/MMbL+HyDW0rfIW1oduRw7+SxNM5/uDfn+Q/Hzs967Fpx3jJWCrHfQeG+NeHTvD06YmS4zJGoSYpxC1BDz6Pq0Qc7CI4y62UzuWZmiMeslQcH02QMQp0hn12R9hQmVTWyWSOtjmG76wkbqvqfFSLQ8PQ1OKggplffefV/NlvXGCv+zwubt+zoaqOrI1IZ9hHPG2seJW0urH2l5kWplx9fo+LWNpgYNI85ujIdMlxZrbS0r+NhRBWOmsx5jAzIP0PPz7Ef/nsY0v+s2ey33IpvWxXsU19YKbloMShjl2bzqpzTf3T3OJgmMFMIURd+mlrheo3NBzL8NH7Dq+YFaFiCWcnZ4uDshzWtwWJpXL2uMmjw2XEoUbFh91WrUOlxnsDkykGp2o/BvOFszF8HhfXbC3OwFJuJd+MOoe5xnauNJ1lZnNr6pemFodMrtAQLS+WGlUId88Lg3zs/iN8+sFjK3IdKgtpYA5xWNcWxChI22KYLQ61+xv2RP22W8klijdk5VZK5fKksrW1vvIFyT3PD7J3cwfr2op1CzMrpNO5vGk5VBh9Wg90NkgihMak+e6MDmq566xnlOXwwIvmjKavP9m/IvMdVP3C2cnUrH7/qr3HeuuGeMgKyB4fKfbBklKSMWr3N+yJBmy3ks/jsnfptjhk86Ry+ZrOKnjw8DADU2nevHcjvS3FsZwz6xwmUzmy+UJdWw4d2nJoKLQ4NEDh2lKjLId9JyfwugWTyRz3vDC47NehYg4ZozAr/13FHNRuOWOYN+izkykSllWRy0sKsnaV6d1RP5PJHPG0gc/tsiurVbZSOpenIGtb9/Cfj52hK+Ljlgt7S8Zyzow5DMXSAHWbrQRmlfRkKoexwr2pNNXR3OJg1Ka6tt7psiwHoyB55YVr2NIV5iuPn1n263C2xZjpWiq6lYo3xL1WqrGyHlTjxNpZDn772nweN163GZdSdQ9Jy6VUK9fScDzNAy8OcfueDfg8LloCHvv9GppR56DEYSHp18tNZ9iHlNgNApeTzz18gtd/8pFl/7mNTPPdGS2MfIF8QTalWynkc9s3mUs3tHLV5nZOjS3/pC5nW4yzMzKWVAHcekd/oJftNLN1jo6YLiblCqvVHI1eq+L96PC0PTPC6xYlMQfn16XmsePjFCTcahVnmhlU5jXNjDmcm6p/y2ElW2j88tgYz56ZpFDQ40qrpWnFIW3t/prRchBC2NbDpX1tBL1uexe+nMTSObu2ZGbGUjxt4BLQ21K0HK6xpoqpoLTqqlurgPRVmztsV5baoXvdLrtCWolTrSyHJ09NEPS6uXBdi72mrBnlDvW4BH6Pyy6Sa69idOdK0WG30Fj+oPSJ0WkKEqazlavt651/efAYf/ClJ5ft5zXfndGiOHu4+SwHMIPSLmEOYgl43TXPuilHLGWwpiVAyOdmYDJd8th0xiDi95S0xtjUGWJTZ8h2K2Vq7FaK+D3cZNUWqDGkPo+rJCANtbMc9p0a59INrSVdZHtbAvg9LrsGRwjBh267wI7ftNWz5WB1Zl1uyyFfkJwZNzcfMwciOSkUJMPxdMXHV5pvPNnPTw8OLVvMpunFoRlTWQE2doS4aF0rYb+HgNdNxigsu8kdS+doCXpZ1xbk7GSpW0vN01D9rVoCHqIBL2tazAwiKHbVrUURnOLVl6wDKLEcsnmJlNIWhVpkeiUyBgcH4+zZ1FGyvntNlA0doZK1t750M/94+6W8fHePfQOuR5TLa2KZxWFgMmUnDUzOIQ4/OTDEdX/3M7v1ea357EPHue/AUFXHDsfSHB2exihIzpQpGq0FTdt4z76xNKnl8Nevvcj+wKjXIGMUbF/2chBL5awgaznLwZzE5/eY8ZH17eYNsSvi57l+c0Lbclh/L9/dQ8BbTGP1uV1kjYIpppaWJmtgdT1zZpJ8QXLl5vaS9T+8aTt3vWzrrOPfeGUfb7yyb8mvYylRrT2WOyB9wjEjfuacbScDkymrmHB52pD8y4PHuWxDG7dc2DvvsY8eH7O/Pz4yzZaucC0vDdCWQ0PMYKgFbSFfMbjpLRZS1ZKfHhjifz9QHDIfSxsOy2F2zEG1q2gJeFlvZS11RnyMWS0Y0jWOOYBZFf3OG7bZ7iUVkHa+VrVwye07OYEQcMXGUnFwu0TDvmd9HhcRv2fZp8GViMMcwpSqcQzJSaEgGU9kGKnSjfXLo2N2hpqz1qeWNK04ZJo4ID0TtfOudVD6y4+f5h9/cphnzpg7f9Ny8HLB2ijjiSyHh+L2sdMZg4jlUnrny7bx27+2ETAth3jGIJ3LL1vc6E9u2ckfvXwHYAWkHR1joTYxh+f6J9neHWmYNvDV0hbyznmDrgUnRhOo7jhzNUqsdQzJyUQyS0Fiu0jn45fHR7l+RxftIS/HR6fnP2EJaNo7Y6bJA9JO1GtQ6x2TarD3Tz85BKiYg4fbXrIWt0vwrafOMp0xePDwCBPJrD2J7+3XbeHlu03TWxXwjSWyDoFfvr+hLQ6O16oWFteJ0QTbeyJL/rwrTXvIt+yWw8mxBNu6zddyLnFQ7sHl6Bagij5H4pl5Y31nxpOcGU9xzdZOtnZHls1yaN6YQ40zXRoJ23LI1S4LQkpJ/0SS1qCXh46M8sTJcaYzBi0BL10RPzfu7OY7T5/lhbNTPHx0FIDrtnfPep5OexZxxmE5LN8ex+cxA9LOOMNSxxyMfIEzE0ledfGaJX3eeqAt5F32mMPJ0QQXrW/l7ERqzmBzKmemuS7FJunkaIKPP3CEseksn//dq2Y19hy1LAajIJlIZu2WNuXYd2ocgF/b2sn+gRg/t2bc15qmtRyKAemmfQls1GuwlOZ01ijw0fsO260uJpM5Etk877huC0LAj54fREposdwmr79iPediaR4+Osp7Xr6dd1y3hTftmR1gVYVUo9OZFRF4n1XnkK6hW2lgMk0uL9nSWfug43LTHvItWzYQmH2wzkyk2NIZnuXSeubMJKcdxZ/JJXIrTSVzvOafH+ZbT53lwcMjZV1HzrkW87mWnjg5QTTgYWdvlK3dEUbiGbv3WC1p2jujvets0ODeUqKauGWW8Cb31OkJPnb/ETtVT7mUdq6JsqUzzBMnzd2QmtP9igt62dQZ4veu3cL7XrmLP3/1hVw+IxgLxdYfo9PZZQlIz8TrEWRnxBzSS2w5nLBGgW7qDM1zZOPRHvIua53Duak0+YJkY0eI1qDXdisZ+QK/9elHuemffs7f/OAAsHhxePr0BN95+qz9//tfHCKeMXjvzWacqpwbyDn9UIlDxsjz+UdOzOo8vO/kOFduasftEnaW0nK4lppYHKwbi7YcijGHRYhD/0SSd//nUyRnVJ6q1hgqyNw/Ye7Q+tqD7FoT5cBADChaDgGvmwfefyN/8ZoL5/x5ynIYm86uSCFjuZjDUlsOpyxxWI50xeWmLWQOmlquQi7Vc6q3NUBr0Muk9b4cmc6QMcwutp99+ASTyaz9N60m5vDZh47z15aofPbhE3zo28+Tt2IH975wjjUtAW63LN9yAWTnDJVh6xp/fmiEv/z+AV7xkQf5068/S74gmUxmOTw0bY8w3tYdrvicS03T3hm15VBE1TYsJubw8JFRfvjcoH2zV6iuqoeHzDexshz62kPsXtNi1wi0BIrZOO4qJu+FfB5CPjej1ocbltdyCHjcdqtuxVKLw4nRBCGfm+5o/Ra0LRZVCDe5DONVobgr74n6aQt57QrpQasX1d4tpnUaTxv2BqeamMODh0e494VzgBk/SGbznBpLkMyaCRWvuqiXda1BAl5XBcshaxd4jlhCoTYFd169ia8/2c///N4L7DtpjsXds8m8zo2dIX7jkrV2GnotadqA9EpkutQrSiAXc5NTO6Czkyn2ONbVIB+n5RANmO0wdq2J2sctJlXTrHXI2H2FlnOKnwqoqhuIxyWWPCB9aizJps7wqpxOqHo/TSaztouwlijLoSfqpy3oYzJlplGrRoXbe6LAOaYzBqlcaUPFuYilcozEM0gp7c/A/oEYh87FyRgFbr14LS6XYHNnuKTOQjE6naWvPcSZ8STDMSUOSdpCXv76dRcT8rv59IPH+eaTZ/G6BZda8+z9Hjef+O0rzu9FqZKmFYdmb5/hJOBbfBHcqFWQNrOITY0APTORJJk1ODORos+qcr5gbVEcWoILfwt2RfyMJbK0Br3LLu5tVkBVCUJ72LfkqY8nRxPsdrxGq4n2Za6SHo5n8LoF7SEfrY6AtLIcdljpwomMQUpZDlX8Paes4UqxlGF/BvYPxBicStER9nGVVdm+rTtizwB3MjqdoSviI5PLM2JZN6fHk2yyWqN84Nbd7F4T5XMPn2RLV3hFNrHndWcUQpwUQjwvhHhGCLHPWusQQtwnhDhifW13HP9BIcRRIcQhIcSrHOtXWs9zVAjxcbEMW6a0kcfnaGDWzBRTWRd+k1Mm8cx5DCqbQkqz5XX/RJI+q/32hvaQXe3ZshjLIWyO76zliNBKdIS9GIXibrEj5FvS+hCVxrppFWYqgUMclikoPRzL0B3x43IJWoNeMlam2eBkymrLYr4n4xmjWOdQxd9TzSIZmErZQe6nT0/wwMFhbrmgF4/VLHFLV5gzEyl7DrliLJGhK+K355SDaTmovllCCF5/eR/f/6/X8fE3X74Er8TCWYpP1k1SysuklMqr8AHgfinlDuB+6/8IIS4E7gAuAm4FPimEUHL4KeAuYIf179YluK45yeQKNW3Y1kgot9KiLIe4EofSNgCxdA6lu4fOxemfSNni4HIJdvZGEQIivsVYDj7GElnSNRwRWgnVc2dwKo0QpltsKWMOZyZSqzaNFYpdY5erEG44nqbHavuufvZUKsdgLM3a1iBRq0XLdNqoukJaSmknXKjxtT63i8dOjBPPGNzqqE/Z2h0mX5CcHi9tLDk2naUz7KPHaiSZyxc4O5mqqwy1WtwdXwt8wfr+C8DrHOtfkVJmpJQngKPAXiHEWqBFSvmoNIfxftFxTs1o1vnR5fC6BW6XWFRAerSC5RBLGWzuCuPzuHjw8AjJbN52KwFcvL6FzrB/UZZbV8TPeMLMLlnuOhW18x2YTBH0ugn63Lav+nzJ5Qv8+Xeex+sWsxrurRZUzGHZ3EqxjD0Doy2o4h05zk2lWdMSsFu0TGcMkvbwprn/nqlcHsPKqDg4aCZi7LH+XhG/h5du77SP3WpVZj/w4hAPHTGL15JZ00rpjPjpifoZjmUYmEyRL0g2ddTPpuB8P1kS+IkQ4kkhxF3WWq+UchDA+tpjra8HnLMo+6219db3M9drSjqX12msFkIIAh7XIgPS5WMO8UyO9pCPbd0RfvDcIG6XsP2wAO+/ZRf//o69i7rezoiPfEEyFEsvu8CrbJuBKUscvG7bV32+/ONPDvHI0TH+9g2X2O0eVhthnzludXktB0scVKZUMsu5qTRrWwN2c8fxRNZORZ3PTehswXHQshzUlMKX7+4paYyo0pH/149e5C2fe5ypZM5uHNkZ8dET9ZPK5e1sv42ryHK4Vkp5BfDrwLuFEDfMcWy5LaKcY332EwhxlxBinxBi38jI+ZWQZ4yCTmN1EPS5q3Ir/e09B7n/oFnYljUKTKVyRP0e4mnDzlAC03JoCXj4zSvW84oLevnRe67nkr42+/H2sI/da1pmPX81qFYD+wdiduxiuVBupXNTaYI+NyGfe8ncSg8cHObGXd1133r7fBBCmEH9RO0th4yRZyKZo9dK+1SZcRPJLEOxNGvbAoQtt+aIo0p5vs+BSrYAeNGyHG7c1cOvbengzms2lRzbGvTyrhu3cfPuHqQ0xUpZ290Rvz3p8IEXh4H6Knw8L3GQUg5YX4eBbwN7gSHLVYT1ddg6vB/Y4Di9Dxiw1vvKrJf7eZ+RUu6RUu7p7p7dd2chaLdSKX5PdTe5zz9yks8+dAIojnu8eH0rAIOOuIMa5POO67fy2bfuKUlfPV8uWd/Kxo4Qr75kLX9229xFc0uNshxyeUnQ6ybgc5PKnr9bqVCQnBpPsrN3dWYpOWkPeZfFclA3fGU5KHE4NpLAKEjWtAZxuQRhn9tOrID5Yw7OTZCqo+hrD/LVd15jF6s5+e+37uYd12+1r2nUYTm8bGc3LQEPX3+yH5/HZQtZPbBocRBChIUQUfU98ErgBeB7wFutw94KfNf6/nvAHUIIvxBiC2bg+XHL9RQXQlxtZSm9xXFOzUjnCrqvkoOgz23PZK5EOpcnYxR48tQEqWzeNo9VDrYz7mBOcqtNpvTmrjC/+G838f/fcTkv6Wutyc+oRGvQa7d/DvpMt9JSpLIOxtJkjUJd7RxrhdlfqfaWw1BMiYN5w1XxjkePmYNz1lrrkYCHkZhDHOZxK6lgtMeKlwW9bkwrCdAAABTNSURBVML+ud/rqqBxZDpjWw6dET/tYR/vsdpsbOwI1VX25PncHXuBh4UQzwKPAz+UUt4LfBi4RQhxBLjF+j9Syv3A14ADwL3Au6WU6q/wLuCzmEHqY8A953Fdc2LkCzx1emJFMl3qmYB3/piD+lBk8wWeODlu77YutW7QKu6gsjmc1c+rBY/bZf9eASvmkMwamLkUi+ekVSi1WrOUnLSHfPRPJGveY0kN0lEB6Yjfwy0X9tpdf9e0muIQ9nvs97Lf4yoR+28/3W+7URXKclDxAdXSZS5scYhnGJxM4RLQa6295ZrNbO0Oc8HaxblZa8Wit3ZSyuPApWXWx4CbK5xzN3B3mfV9wMWLvZaF8LH7j/Cpnx+jI+wr8YE3O9XsgJ3m9CPHRtluBU13r23B6xa25aCyORZTw9AItIfMBm4hn5mtVJCmYJ7PhLaTqtneKuynNJNbLuzlJwfO8bK//xlffec1XLiuNjfFYuuMoqvmL159Ib84PELGKLDWEoeo32M3u+sM+0o2SX93zyE2dYa4+YLiKE/VgmObNVuhmkrvloAHn8dlu5V6WwJ2LYTP4+K7774Wj6u+PBn1dTXLwO1XbqAgJcPxjHYrOQh45485qCwNn9vFI0dHbd9pT9TPmtaAPfhc9VWqlVtppVFBaZWtBJA+z7jDqbEkPo/LdnWsZn7zSrO4K54xeMTaxdeCoVgat0vQGS7u7Dd0hPh/X7mL7T0ROqz1SMBDwnIldUSK4jCRyHIulrarqRWqAG6r1QSvGnEQQtAdMYs3ByZTrGsLljweDXiXdX57NTTd3XFjZ4hXXWQWqTTqLN5a4Pe4561zUOJww84u9g/E2D8wZftbr9zYzs8PDZPIGLb7aTW6laAYlFZ1DnD+zfdOjCbYVGc+51py4doWPC7BeA0D087qaCe/f8NW7vuTG+zeVWFHIWZH2G/HHF600lTPTaVLprXFLKtxXat5g++Ozu9WMo/zMzKdYWBqtjjUI00nDgDvuH4LoAf9OKkmlVWJw53XbEYAP3x+kC7rg3HnNZuJpw2+88xZ2/20et1K5u8c8BUth4WKQ6EgefLUuN1W4dRYYtW2zCiHEIL2sK+mbTTGEtmK8QBnh56Iw8LtDPvIGAUKBcmhc2aaajZfsMd6gulebQ16bYtBTSecj+6onyHLElnXVv8WYlPeHa/Y2M47b9haUube7ARmBOLKoXytF69r4ddfshYpiyb1FRvbuHh9C1/85Sk7D7xllbqVVNZL0Ou2kxpmzrOYj28+1c9vfupRXvrhB/jy46c5NZZkS9fqz1Ry0hHy1TQoPZbI2q6juYg6Mo2U8KeNvG05QLGLK5ibpJaA1w4yd1URkAZTHE6MJsgaBdvqqGeaUhyEEHzwtgu4fsf51UqsJqqzHKybftDLO28w87aVOAgheMvVmzk0FOehI6YfObrK3UohqwgOFt6X6sf7z9ET9bO1O8wHv/U8GaPQVJYDQHvYW9OU1vFEpiTeUAlnGqqyNFLZPAfPxe3aiIGpYpp2LGXQEvRY3VJd7K4yy6g74ieXN91T2q2kaRiqDUib7Q9cXNLXxu9fv4VXX7LWfvym3WanlPsOmkNQFtOOuxFQAemAM+awgIB0KpvnoSOj3PaStXzp7Xu5bnsXsDonv81FR9hX05jD+HTWrqafC6dbSVkayWyew+fidluMQUcNTyxdtBxe+MtXcfXWTqrBObypEdxKq/PTq1kwAa8ZkJZSVhwyo3ytij/7jdLq5O6on61dYY5bOfurNyA9O1tpITGHh46YqZS3XNiL3+Pm03deyQ+fH+TXtsyurl3NtIVqF3NI5/IksvkFuZV8bpfda+nQuTipXJ5rt3dy7/5zJRlLsXSOXVYlu0pHrQanOKxvAMtBi4MGKAbnM0ahYnHgVCo3b5B5z+Z2jo8m8Lldq7bI0M5W8i0u5vDTg0NEAx72WmIQ9nt4054N85y1+ugI+ZhIZikU5JJnaakA8kLcSkFHgsEzZ8yJcbvXtLC2NcDgVJp9J8fJFyRTyfk/B+VQ4hD0uhc1AXG50eKgASjm68/Rc2oqlZv3Tb1ncwdf29e/al1KUGz8F/F77JhDtQN/skaBnx4c5qZdPXgXsOtcjbSHfRSkWRfTGlram+W4VYNTjeWgrAVnarJqxb2lO8yalgADkyne+5VnyBgF4hljUckW3db7Zl1boCFGwDb3u1NjE6jCPRKrQhz2Wo3HVqtLCWBnb4SPvOlSXnFBr/16TKaqC6w+eHiE8USW112+rpaX2BB0hM3XrhZxB9UUciExh5DDEnzxXJz2kJeWgJd1bUGePjPJ2ckUo9MZpFxcmrayHBohGA1aHDQWyq00VyFcNZbDps4QXRE/0QYwmxeLEII3XNFH0GcWAAa9bsYcXT3n4ltP9dMZ9ulMOYqB/Vqks44vwK0UKeNWMqeymQkCa1sD5AsSn8fFOqvlxmI2PwGvm46wj40djZGyvHptf82CsAOrc7hHqhEHIQTvvGFr01T6gpn+qFqJzMVUMsf9B4f5nas3Nr1LCcyYA8w/T3psOsP3nh3grddsrvp9pToGd1RRg6DEQfXKUqgOuWutnf5Nu7q5pK+Nf/jxoUUXeP6ft11l93Sqd7Q4aADwq5iDUV4ccvkCyWy+qg/F71s1EM1CV8Rvt2Gei/sODpHNF3j95TUfdNgQdNgjQ+cWhx8+P8j/9/0DXL6xncs2VG6W6cy0G0tk8bpFSYFbJZRbKejz2JskgE3WDl/NPn/Npeu4bnsXx0cS9ljQhXLpHNdfb+jtiwZwBKQrWA6qdUYjZFksN11VWg4HBmIEvC4uXre8MyjqlfYqxUEN7XnixLi9dnYyVTK97chQnAv/4sc8bh0znsjQEfZVFfi1LQevu6SlzkbLrXT99i4++TtXcNvFa2kL+finN11aVbO9RkeLgwYoBqQrWQ4xLQ4V6Qz7q4o5HBmOs6Mn2lQut7kI+9z43GYb6/d8+WmeODle9jglvI9ZN34pJXd+9jHe9e9P2sd895kBUrk8n33oOGDGHKrteRT0unGJ0tRkKLqVPG4Xt71kbdP93bRbSQM4U1nLB6S15VCZrqjZI2i+fP3DQ3GutaqhNWqetJefHBji1FiSaMBTdsymEt4nTo6bDfGG4hwfTXB8NMGxkWm2dUe454VBwKwhGZhMzdl0r9x1tAa9RAMe/B4XQoCURbdSs6ItBw1QzFaqFJBW4rBaO62eD51hP0ZB2q9ROaaSOYZiGbuyVmPSEfZxaiwJwAtnpwCYzhjkHS2yVTxnKpXj8HCce184hxDgdgm++WQ/R4biHBtJ8M4btiKBrzx+mrHp6pruKT7x21fw+9dvRQhhV747K5qbEW05aABoC5ofpHOxdNnHteVQmS7rJjKWyNh+9JkcHjY7fO7U4lCCakUCcPBcnOmMwcv+/me882VbueuGbYAZXL5iYxtPnZ7ksePj/Hj/Oa7a1EEk4OFbT51lPJFFCHj7dVs4NjLNFx49RcaornWG4qUOiy7oddMV8TdEoVot0ZaDBoDWkJcL1rbw8JHyk7l0zKEyXdZNaCReObB6eMgUhx29kWW5pkZB3cBv2tVN1ijw+UdOMJbI8lz/lH3M2HSWSze0saUrzP/60UFePBfnVRev4beu2sC5WJqvPHGGq7d00tMS4P2v3EU8nSOdK1RV41COkN9tz4duZrTloLG5YUcXn3vkBMmsQcgxHStj5PmPx07THfXbfYU0RVQVrqrKLceRoWnCPndDNFxbTjojPtwuwXtfsZOfHRrhM78wA8pqpnY6l2c6Y9AV8fOF393LPz9whMdOjPMbL1lLb4uf7/3Rtbhdwu5oe8HaFu7Yu5H/fOx0VdXR5fir115Mb7QxahFqiRYHjc31O7r59C+O89jxcbv99uh0ho/99AgvnovzubftWVAXymZBDXsZjVcWh0Pn4uzojTa9q2Im77huKzft6uGS9a1EAx57PvOp0SRSSjve0BXxsbEzxD/cfmnJ+Zf0za4beP8tOzk3lbYbGy6Um3b1LOq81Yb+pGts9mxux+9x8YsjI4DZ6mHP3/yUL/3qFHdevYmX7+5d4SusT9pCPlyCklGSirHpDO/72jM8cXKcC6ocCtNMbOwMcdPuHlwuYdd/XNLXSjxjMJ7I2mms1aalgmnJfe5tV7GtW7vwzgctDhqbgNfNr23t5MHDpjh85fEzbO0O860/fCl/9dqLVvjq6he3S9AR9pcthPvY/Uf4/rMDvHnvRt7/yp0rcHWNw+Ub2/C6BW976WbAdC2pNNauJs8cWgm0W0lTwi0X9PA/vrufXx0f48nTE/zBy7ZyxcbFtQpoJswq6VK3kpEv8KPnB3nlhWv469ddvEJX1ji868ZtvObSdfg95p715GjSTmldbHBZs3i05aAp4dWXrMPrFnzoW8+TL0hu1P7XquiM+GZVSf/y2Bij01lec6luz10N0YCZMdfXHsIl4NRYghE75qAth+VGi4OmhPawj5t29XB8NEFLwMPlDdQobCUxm++VupW+/+wAUb+HG3fp9twLwedx0dce4uRYkrHpLOEZ3VI1y4MWB80s3nBFH2BmL+nspOrY0RPh9HiSbzzZD5hVvvfuP8crL1qzasel1pJNnSEz5pDILDolVXN+6JiDZhY37e7mxl3dvHnvxpW+lIbhrhu28avj4/z3bz5Hd9TP4XNx4mmDO6/ZtNKX1pBs7gzznWfOEg14qu6RpFla9LZQMwu/x83nf3cv1+3QTeKqxedx8S93XsmOngj/9T+f4tO/OM5Lt3XOOX9AU5mXrG8lnjb41fFxHW9YIbQ4aDRLRMTv4V/fsgeXSzA6neEPb9y+0pfUsLzxyj7eecNW8gXZ9A3wVgrtVtJolpANHSE+97arePTYGNdu71zpy2lYXC7BB2+7gL1bOnSzwhVCi4NGs8RcsbFd14YsETdfoKvyV4q6cSsJIW4VQhwSQhwVQnxgpa9Ho9Fompm6EAchhBv4BPDrwIXAm4UQF67sVWk0Gk3zUhfiAOwFjkopj0sps8BXgNeu8DVpNBpN01Iv4rAeOOP4f7+1ptFoNJoVoF7EoVyTeznrICHuEv+3vbsNtayq4zj+/TF3FIYxTRzjitoElYRjaA6BpPQAJoa96QENcUYFKSjKepNG0osMNFJMg6ZLjpQ90KM0KjhYNuYDhRMV4zDqoERmg46p46QYaL9e7HXpcM65Tvfevc/e+57fBw5n33X2Xfz/f/Y9a+9191lH2ilp5/79+ycQVkTEdOrK4PB34ISBn48H/jG8k+052xttb1y3LuvVREQ0pSuDw0PA2yS9RdJhwAXAtpZjioiYWp34nIPtVyV9BtgOrAK22t7dclgREVNL9sjUfi9IOgg8uoRfPRI4UGMoxwDP1thfnfHVnes01a4P/UFquFzTVL8Tgb8Bb7Z96Hl52718ADuX+HtzXYhjEvE1kOvU1K4P/aWGqd8i+9q/mP278j+HSbq97QAOoc746s51mmrXh/6a0PWcu17DLuf7wmJ27vO00k7bGxNHf6V2y5caLs801W+xufb5ymGu7QCKrsTRR6nd8qWGyzNN9VtUrr29coiIiOb0+cohIiIaksFhiKQTJP1W0h5JuyV9rrQfLeluSXvL8xtL+9mS/ihpV3n+QGlfI+lOSY+Ufq5pM69JqKt25bW7JP2l9LOlrNy74tVZw4E+t0l6eNK5tKHmY3BH+RqBP5fHsW3l1Yo6b7taCQ9gFnhX2T4CeIxqGfGvA1eU9iuAa8v2acBxZXsD8FTZXgO8v2wfBtwHnNt2fn2oXfn5DeVZwC+AC9rOr281LG0fAX4EPNx2bn2rH7AD2Nh2Tq3Vsu0Auv4AfgWcTfWBu9nSNgs8OmZfAf8EDh/z2jeBy9rOp2+1A1ZT3c53ftv59K2GwFrg/vLmOBWDQ831m+rBIdNKr0PSeqoziz8Ab7K9D6A8j7vE/CjwJ9v/HurnKODDwG+ajLdL6qidpO3AM8BB4OcNh9w5NdTwq8B1wMuNB9tBNf393lKmlK6SNG716BUrg8MCJK2lms643PaL/8f+JwPXAp8cap8BfgzcaPuJJmLtmrpqZ/scqrO8w4GRufSVbLk1lHQq8FbbtzUaaEfVdAxeaPsU4KzyuKiJWLsqg8MYklZTHVg/tP3L0vy0pNny+izVGe38/scDtwGbbD8+1N0csNf2Dc1H3r6aa4ftV6hW6J2abwasqYZnAKdL+ivV1NLbJe2YTAbtqusYtP1UeT5I9X+bd08mg27I4DCkXDreDOyxff3AS9uAzWV7M9Vc5vyU0Z3AlbYfGOrraqqFsy5vOu4uqKt2ktYO/CHPAB8CHmk+g/bVVUPb37Z9nO31wJnAY7bf13wG7arxGJyRdEzZXg2cB0zFHV/z8iG4IZLOpLqzaBfwn9L8Jap5y5/yv5UNP277OUlfBq4E9g5080GqO5SepHpTm5/D/Jbt7zaeREtqrJ2AO6imk1YB9wCft/3qJPJoU101tD14ZrweuMP2hsYTaFmNx+BLwO+obohYBfwa+ILt1yaRRxdkcIiIiBGZVoqIiBEZHCIiYkQGh4iIGJHBISIiRmRwiIiIERkcIhog6VOSNi1i//XTsnJq9MNM2wFErDSSZmxvaTuOiOXI4BAxRvng2F1UH546jWrp503AO4DrqVY8fRa42Pa+sjTFg8B7gG2SjgD+ZfsbZZ2jLVTLuD8OXGr7eUmnA1upFsa7f3LZRRxappUiFnYSMGf7ncCLwKeBm4CP2Z5/Y//awP5H2X6v7euG+vk+8MXSzy7gK6X9FuCzts9oMomIpciVQ8TCnhxYb+cHVMswbADuLqs3rwL2Dez/k+EOJB1JNWjcW5q+B/xsTPutwLn1pxCxNBkcIhY2vLbMQWD365zpv7SIvjWm/4jOyLRSxMJOlDQ/EHwC+D2wbr5N0uryPQALsn0AeF7SWaXpIuBe2y8AB8pCcQAX1h9+xNLlyiFiYXuAzZK+Q7Vq503AduDGMi00A9wA7D5EP5uBLZLWAE8Al5T2S4Ctkl4u/UZ0RlZljRhjmpa5jhgn00oRETEiVw4RETEiVw4RETEig0NERIzI4BARESMyOERExIgMDhERMSKDQ0REjPgv03qo4p9lptUAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "sorted_data['inc'][-200:].plot()"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "first_september_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n",
+ " for y in range(1991,\n",
+ " sorted_data.index[-1].year)]"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "year = []\n",
+ "yearly_incidence = []\n",
+ "for week1, week2 in zip(first_september_week[:-1],\n",
+ " first_september_week[1:]):\n",
+ " one_year = sorted_data['inc'][week1:week2-1]\n",
+ " assert abs(len(one_year)-52) < 2\n",
+ " yearly_incidence.append(one_year.sum())\n",
+ " year.append(week2.year)\n",
+ "yearly_incidence = pd.Series(data=yearly_incidence, index=year)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 13,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAD8CAYAAAC2PJlnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHGxJREFUeJzt3X+Q1PWd5/HnC4fMqAsGdDD8UHErxAqarIYpYC+5vTXsAl62BGs1mdONUydVbNTNmq2rUoheWadslW5dXXYpK65UXEWNPzhuLbk9OR3xUsvdEWCIZhUNO5NokEic8YYIpopZB973R39GeibQ0/2dnun+Nq9HVVd/5z3fz2c+/WWYd39+fD+tiMDMzCyLSbVugJmZ5ZeTiJmZZeYkYmZmmTmJmJlZZk4iZmaWmZOImZll5iRiZmaZOYmYmVlmTiJmZpZZU60bUG3nnXdezJ07t9bNMDPLlT179rwfEa2Vlmu4JDJ37ly6urpq3Qwzs1yR9PMs5TycZWZmmTmJmJlZZk4iZmaWmZOImZll5iRiZmaZOYmY1ZHew0f56kM76D1ytNZNMSuLk4hZHVm/rZvdb/ez/qXuWjfFrCwNd5+IWR5dctdWBgaPf/z1Ezv388TO/TQ3TWLfuqtq2DKz0twTKeKhBKuV7bdfydWXz6JlcuG/ZMvkSay4fBbb77iyxi0zK81JpIiHEqxWZkxtYUpzEwODx2lumsTA4HGmNDcxY0rLSc/3Gx6rFx7OwkMJVh/e/3CAGxZdxPULL+TJXfvpK5Egit/wrLvmcxPYSrPhFBG1bkNVtbW1RaV7Z/UePsq659/kxb2/5OhHx2mZPIlll36KO7/y2VO+EzSrhZFveIb4DY+NlaQ9EdFWaTkPZ1H5UIJZrXjuxOqNh7OSSoYSzGrFb3is3jiJJA99/UQvbt3Ky2rYErPS/IbH6onnRMzMzHMiZmY28ZxEzMwss7KSiKS/kLRX0uuSnpLUImm6pE5J3el5WtH5ayX1SNonaVlRfIGk19L31ktSijdLeibFd0qaW1SmI/2Mbkkd1XvpZmY2VqMmEUmzgT8H2iLiMuAMoB1YA2yLiHnAtvQ1kuan718KLAe+K+mMVN2DwGpgXnosT/FVwKGI+DTwHeD+VNd04G5gEbAQuLs4WZmZWW2VO5zVBJwpqQk4C3gXWAFsTN/fCKxMxyuApyNiICLeAnqAhZJmAlMjYkcUZvMfG1FmqK7NwJLUS1kGdEZEf0QcAjo5kXjMzKzGRk0iEfEL4D8D+4GDwAcR8SJwfkQcTOccBGakIrOBd4qqOJBis9PxyPiwMhExCHwAnFuiLjMzqwPlDGdNo9BTuBiYBZwt6U9KFTlJLErEs5YpbuNqSV2Suvr6+ko0zczMqqmc4aw/AN6KiL6I+Aj4e+BfAe+lISrSc286/wBwQVH5ORSGvw6k45HxYWXSkNk5QH+JuoaJiA0R0RYRba2trWW8JDMzq4Zyksh+YLGks9I8xRLgTWALMLRaqgN4Lh1vAdrTiquLKUyg70pDXkckLU713DiizFBd1wIvp3mTF4ClkqalHtHSFDOzBuft7vNh1G1PImKnpM3Aj4BB4BVgA/BbwCZJqygkmuvS+XslbQLeSOffGhHHUnU3A48CZwJb0wPgYeBxST0UeiDtqa5+SfcCu9N590RE/5hesZnlgre7zwdve2INp/fwUf7sqVd44PorvDFhDnm7+9rwtidmiT+hMt+83X2+eBdfaxj+hMrG4O3u88U9EWsYfgfbOIa2u3/2li9yw6KL6PtwoNZNslNwT8Qaht/BNg5/vk9+OIlYQ/EHNplNLK/Ostzwqiuz8ePVWVY19XqTl1ddmdUfD2fZb6i3m7y86sqsfnk4yz5Wrzd59R4+yrrn3+TFvb/k6EfHaZk8iWWXfoo7v/LZXA1reTjO6pmHs2zM6nWJbKOsuvJwnDUiD2fZx+r5j3WeV115OM4amZOIDVOvf6zzfN/A9tuvPOVwnFneOYnYMHn+Y12v6rmHZzZWTiJmE6Bee3hmY+XVWWZm5tVZZmY28ZxEzMwsMycRMzPLzEnEbIzqda8xs4ngJGI2Rr4T3U5nXuJrlpHvRDdzT8Qss3rda8xsIjmJmGXkO9HNykgiki6R9GrR47Ckb0maLqlTUnd6nlZUZq2kHkn7JC0rii+Q9Fr63npJSvFmSc+k+E5Jc4vKdKSf0S2po7ov32xshu5Ef/aWL3LDoovo+3Cg1k0ym1AV3bEu6QzgF8Ai4FagPyLuk7QGmBYRd0iaDzwFLARmAS8Bn4mIY5J2AbcBPwSeB9ZHxFZJtwCfj4hvSGoHromIr0maDnQBbUAAe4AFEXHoVG30Hetm1efPQml8E3XH+hLgpxHxc2AFsDHFNwIr0/EK4OmIGIiIt4AeYKGkmcDUiNgRhcz12IgyQ3VtBpakXsoyoDMi+lPi6ASWV/oizWxsvALNTqXS1VntFHoZAOdHxEGAiDgoaUaKz6bQ0xhyIMU+Sscj40Nl3kl1DUr6ADi3OH6SMmY2zrwCzUZTdk9E0ieAq4H/OtqpJ4lFiXjWMsVtWy2pS1JXX1/fKM0zs3J5BZqNppLhrKuAH0XEe+nr99IQFem5N8UPABcUlZsDvJvic04SH1ZGUhNwDtBfoq5hImJDRLRFRFtra2sFL8nMSql0BZrv3j/9VJJE/h0nhrIAtgBDq6U6gOeK4u1pxdXFwDxgVxr6OiJpcZrvuHFEmaG6rgVeTvMmLwBLJU1Lq7+WppiN4P+8Nl4qWYHmuZPTT1mrsySdRWFu4rcj4oMUOxfYBFwI7Aeui4j+9L07gZuAQeBbEbE1xduAR4Ezga3ANyMiJLUAjwNXUOiBtEfEz1KZm4Bvp6b8ZUQ8Uqqtp+vqrLuefY3v79rPDQsvZN01n6t1c+w0M3LuZIjnTvIj6+osfyhVzlX6n9dLNW089B4+esrPkffvWT74Q6lOU5VOfHq4wcaD794/fXkDxpwr9z+vl2raePPnyJ+enEQaQDn/ebfffuUphxvMquGhr58YCVm38rIatsQmkpNIAyjnP6+HG8xsPHhOZJzU45JbbxZoZtXm1VnjxEtuzSxPsq7O8nBWlXkC28xOJx7OqrJKltzW45CXmVklnESqrJIJbN+zYWZ55+GscTDaklsPeZlZo3BPpELlDEE99PU21q28jPmzprJu5WXDluCCt9e2ieHh0vrXCP9GTiIVqsYQlO/ZsIng4dL61wj/Rl7iW6Zq71L6p4930TqlZdiQ18gei1kW3lG3/tXjv5F38U3GK4l4l9LTUx53Pfbvav2rx38j7+I7zjwE1VjKHYvO43CDf1frXyP9G3l1VgW8S2njKE4OJ9tRIO8r6Py7Wv8a5d/Iw1l2Wil3LLoehxvGSx6H7Kz6PJxldaneljCWu7y6kYYbRpPHITurHx7OsszKeQc72rDRRKskOTTKcMOp5H3IzuqDh7Mss1I7FdfjEsYhXl5dcDoN2dnovIuvTZhy3sHW8ycp+hP4Ck6nITsbP54TsYqVM6/gP1D54A8qs7FyT8QqVm6CaPQ5hUZQi16ZV4M1lrJ6IpI+KWmzpJ9IelPS70qaLqlTUnd6nlZ0/lpJPZL2SVpWFF8g6bX0vfWSlOLNkp5J8Z2S5haV6Ug/o1tSR/Veuo1FOe9gR9uI0k5PXg3WWMqaWJe0EdgeEd+T9AngLODbQH9E3CdpDTAtIu6QNB94ClgIzAJeAj4TEcck7QJuA34IPA+sj4itkm4BPh8R35DUDlwTEV+TNB3oAtqAAPYACyLi0Kna6ol1s/pUz4stbBzvE5E0Ffg94GGAiPiXiPgVsALYmE7bCKxMxyuApyNiICLeAnqAhZJmAlMjYkcUMtdjI8oM1bUZWJJ6KcuAzojoT4mjE1he6Ys0s9rzRyA0pnKGs34b6AMekfSKpO9JOhs4PyIOAqTnGen82cA7ReUPpNjsdDwyPqxMRAwCHwDnlqhrGEmrJXVJ6urr6yvjJZnZRPNii8ZUThJpAr4APBgRVwC/BtaUOF8niUWJeNYyJwIRGyKiLSLaWltbSzTNzGrJq8EaTzmrsw4AByJiZ/p6M4Uk8p6kmRFxMA1V9Radf0FR+TnAuyk+5yTx4jIHJDUB5wD9Kf77I8r8oKxXZmZ1p57v0fGqsWxG7YlExC+BdyRdkkJLgDeALcDQaqkO4Ll0vAVoTyuuLgbmAbvSkNcRSYvTfMeNI8oM1XUt8HKaN3kBWCppWlr9tTTFzMyqyqvGsin3PpFvAt9PK7N+Bvx7Cglok6RVwH7gOoCI2CtpE4VEMwjcGhHHUj03A48CZwJb0wMKk/aPS+qh0ANpT3X1S7oX2J3Ouyci+jO+VjOz3+A9xMbGe2eZ2WnNe4gVeCt4M7MMvGpsbLztiZmd9rxFT3YezjIzMw9nNbJ6+3RAM7MhTiI54KWHZlavPCdSx7z00MzqnXsidcwb1plZvXMSqWNeemhWmucLa89JpM55wzqzU/N8Ye15ia+Z5Y4/4Kr6vMTXzE4blc4Xethr/DiJmFnuVDpf6GGv8eMlvmaWS+VsVeJl8uPPcyJm1rC8Q2/5PCdiZjaCl8mPPw9nmVlD8w6948vDWWZm5uEsMzObeE4iVhe8jt8sn5xErC54Hb9ZPnli3WrK6/jN8s09Easpb3dvlm9OIlZTXsdvlm9lJRFJb0t6TdKrkrpSbLqkTknd6Xla0flrJfVI2idpWVF8QaqnR9J6SUrxZknPpPhOSXOLynSkn9EtqaNaL9zqh7e7N8uvsu4TkfQ20BYR7xfF/groj4j7JK0BpkXEHZLmA08BC4FZwEvAZyLimKRdwG3AD4HngfURsVXSLcDnI+IbktqBayLia5KmA11AGxDAHmBBRBw6VVt9n4iZWeVqcZ/ICmBjOt4IrCyKPx0RAxHxFtADLJQ0E5gaETuikLkeG1FmqK7NwJLUS1kGdEZEf0ocncDyMbTZzMyqqNwkEsCLkvZIWp1i50fEQYD0PCPFZwPvFJU9kGKz0/HI+LAyETEIfACcW6KuYSStltQlqauvr6/Ml2RmZmNV7hLfL0bEu5JmAJ2SflLiXJ0kFiXiWcucCERsADZAYTirRNvMzKyKyuqJRMS76bkXeJbCfMd7aYiK9NybTj8AXFBUfA7wborPOUl8WBlJTcA5QH+JuszMrA6MmkQknS1pytAxsBR4HdgCDK2W6gCeS8dbgPa04upiYB6wKw15HZG0OM133DiizFBd1wIvp3mTF4Clkqal1V9LU8zMzOpAOcNZ5wPPptW4TcCTEfE/Je0GNklaBewHrgOIiL2SNgFvAIPArRFxLNV1M/AocCawNT0AHgYel9RDoQfSnurql3QvsDudd09E9I/h9ZqZNaTew0f5s6de4YHrr5jQ+6y8FbyZWQO469nX+P6u/dyw8ELWXfO5istnXeLrvbPMzHKs1vvPedsTM7Mcq/X+c04iZmY5Vuv95zycZWaWc7X8HHlPrJuZmT9j3czMJp6TiJmZZeYkYmZmmTmJmJlZZk4iZmaWmZOImZll5iRiZmaZOYmYmVlmTiJmZpaZk4iZmWXmJGJmZpk5iZiZWWZOImZmlpmTiJmZZeYkYmZmmTmJmJlZZk4iZmaWWdlJRNIZkl6R9A/p6+mSOiV1p+dpReeuldQjaZ+kZUXxBZJeS99bL0kp3izpmRTfKWluUZmO9DO6JXVU40WbmVl1VNITuQ14s+jrNcC2iJgHbEtfI2k+0A5cCiwHvivpjFTmQWA1MC89lqf4KuBQRHwa+A5wf6prOnA3sAhYCNxdnKzMzKy2ykoikuYAXwG+VxReAWxMxxuBlUXxpyNiICLeAnqAhZJmAlMjYkcUPtj9sRFlhuraDCxJvZRlQGdE9EfEIaCTE4nHzMxqrNyeyF8DtwPHi2LnR8RBgPQ8I8VnA+8UnXcgxWan45HxYWUiYhD4ADi3RF1mZlYHRk0ikv4I6I2IPWXWqZPEokQ8a5niNq6W1CWpq6+vr8xmmpnZWJXTE/kicLWkt4GngS9LegJ4Lw1RkZ570/kHgAuKys8B3k3xOSeJDysjqQk4B+gvUdcwEbEhItoioq21tbWMl2RmZtUwahKJiLURMSci5lKYMH85Iv4E2AIMrZbqAJ5Lx1uA9rTi6mIKE+i70pDXEUmL03zHjSPKDNV1bfoZAbwALJU0LU2oL00xMzOrA01jKHsfsEnSKmA/cB1AROyVtAl4AxgEbo2IY6nMzcCjwJnA1vQAeBh4XFIPhR5Ie6qrX9K9wO503j0R0T+GNpuZWRWp8Ia/cbS1tUVXV1etm2FmliuS9kREW6XlfMe6mVmV9R4+ylcf2kHvkaO1bsq4cxIxM6uy9du62f12P+tf6q51U8bdWOZEzMysyCV3bWVg8MTtdE/s3M8TO/fT3DSJfeuuqmHLxo97ImZmVbL99iu5+vJZtEwu/GltmTyJFZfPYvsdV46p3noeHnMSMTOrkhlTW5jS3MTA4HGamyYxMHicKc1NzJjSMqZ663l4zMNZZmZV9P6HA9yw6CKuX3ghT+7aT98Yeg95GB7zEl8zszrVe/go655/kxf3/pKjHx2nZfIkll36Ke78ymfH3LsZyUt8zcwazHgNj1WTh7PMzOpYNYfHxoOHs8zMzMNZZmY28ZxEzMwsMycRMzPLzEnEzMwycxIxM7PMnETMzCwzJxEzM8vMScTMzDJzEjEzs8ycRMzMLDMnETMzy8xJxMzMMnMSMTOzzEZNIpJaJO2S9GNJeyX9pxSfLqlTUnd6nlZUZq2kHkn7JC0rii+Q9Fr63npJSvFmSc+k+E5Jc4vKdKSf0S2po5ov3szMxqacnsgA8OWI+B3gcmC5pMXAGmBbRMwDtqWvkTQfaAcuBZYD35V0RqrrQWA1MC89lqf4KuBQRHwa+A5wf6prOnA3sAhYCNxdnKzMzKy2Rk0iUfBh+nJyegSwAtiY4huBlel4BfB0RAxExFtAD7BQ0kxgakTsiMKHmDw2osxQXZuBJamXsgzojIj+iDgEdHIi8ZiZWY2VNSci6QxJrwK9FP6o7wTOj4iDAOl5Rjp9NvBOUfEDKTY7HY+MDysTEYPAB8C5JeoyM7M6UFYSiYhjEXE5MIdCr+KyEqfrZFWUiGctc+IHSqsldUnq6uvrK9E0MzOrpopWZ0XEr4AfUBhSei8NUZGee9NpB4ALiorNAd5N8TkniQ8rI6kJOAfoL1HXyHZtiIi2iGhrbW2t5CWZmdkYlLM6q1XSJ9PxmcAfAD8BtgBDq6U6gOfS8RagPa24upjCBPquNOR1RNLiNN9x44gyQ3VdC7yc5k1eAJZKmpYm1JemmJmZ1YGmMs6ZCWxMK6wmAZsi4h8k7QA2SVoF7AeuA4iIvZI2AW8Ag8CtEXEs1XUz8ChwJrA1PQAeBh6X1EOhB9Ke6uqXdC+wO513T0T0j+UFm5lZ9ajwhr9xtLW1RVdXV62bYWaWK5L2RERbpeV8x7qZmWXmJGJmZpk5iZiZWWZOImZmlpmTiJmZZeYkYmZmmTmJmJlZZk4iZmaWmZOImZll5iRiZmaZOYmYmVlmTiJmZpaZk4iZmWXmJGJmZpk5iZiZWWZOImZmlpmTiJmZZeYkYmZWgd7DR/nqQzvoPXK01k2pC04iZmYVWL+tm91v97P+pe5aN6UuNNW6AWZmeXDJXVsZGDz+8ddP7NzPEzv309w0iX3rrqphy2rLPREzszJsv/1Krr58Fi2TC382WyZPYsXls9h+x5U1blltOYmYmZVhxtQWpjQ3MTB4nOamSQwMHmdKcxMzprTUumk15eEsM7Myvf/hADcsuojrF17Ik7v20+fJdRQRpU+QLgAeAz4FHAc2RMTfSJoOPAPMBd4GvhoRh1KZtcAq4Bjw5xHxQoovAB4FzgSeB26LiJDUnH7GAuD/AV+LiLdTmQ7grtScdRGxsVR729raoqurq/wrYGZmSNoTEW2VlitnOGsQ+A8R8VlgMXCrpPnAGmBbRMwDtqWvSd9rBy4FlgPflXRGqutBYDUwLz2Wp/gq4FBEfBr4DnB/qms6cDewCFgI3C1pWqUv0szMxseoSSQiDkbEj9LxEeBNYDawAhjqFWwEVqbjFcDTETEQEW8BPcBCSTOBqRGxIwrdn8dGlBmqazOwRJKAZUBnRPSnXk4nJxKPmZnVWEUT65LmAlcAO4HzI+IgFBINMCOdNht4p6jYgRSbnY5HxoeViYhB4APg3BJ1mZlZHSg7iUj6LeC/Ad+KiMOlTj1JLErEs5YpbttqSV2Suvr6+ko0zczMqqmsJCJpMoUE8v2I+PsUfi8NUZGee1P8AHBBUfE5wLspPuck8WFlJDUB5wD9JeoaJiI2RERbRLS1traW85LMzKwKRk0iaW7iYeDNiPgvRd/aAnSk4w7guaJ4u6RmSRdTmEDflYa8jkhanOq8cUSZobquBV5O8yYvAEslTUsT6ktTzMzM6kA5S3y/BGwHXqOwxBfg2xTmRTYBFwL7gesioj+VuRO4icLKrm9FxNYUb+PEEt+twDfTEt8W4HEK8y39QHtE/CyVuSn9PIC/jIhHRmlvH/DzMl9/PTgPeL/WjRiDPLc/z22HfLc/z22Hxmz/RRFR8VDOqEnExpekrixrs+tFntuf57ZDvtuf57aD21/M256YmVlmTiJmZpaZk0jtbah1A8Yoz+3Pc9sh3+3Pc9vB7f+Y50TMzCwz90TMzCwzJ5FxIOnvJPVKer0o9juSdkh6TdJ/lzQ1xT8h6ZEU/7Gk3y8q8wNJ+yS9mh4zTvLjqt32CyT9L0lvStor6bYUny6pU1J3ep5WVGatpJ7U1mVF8QXpdfVIWp/uD8pL2+v+2ks6N53/oaQHRtRV19d+lLbn4dr/oaQ96RrvkfTlorom9NqPQ/sru/4R4UeVH8DvAV8AXi+K7Qb+TTq+Cbg3Hd8KPJKOZwB7gEnp6x8AbRPc9pnAF9LxFOCfgfnAXwFrUnwNcH86ng/8GGgGLgZ+CpyRvrcL+F0K29dsBa7KUdvzcO3PBr4EfAN4YERd9X7tS7U9D9f+CmBWOr4M+EWtrv04tL+i6++eyDiIiH+kcNNksUuAf0zHncAfp+P5FLbSJyJ6gV8BNVt/HhOza3Ndt30821hKpe2PiF9HxP8Ghn0yUh6u/anaXisZ2v9KRAxtwbQXaFFhl44Jv/bVbH+Wn+0kMnFeB65Ox9dxYk+wHwMrJDWpsE3MAobvF/ZI6lL+x4noFhfT+O3aPO7G2PYh9X7tTyUP1340ebr2fwy8EhED1Pjaw5jbP6Ts6+8kMnFuovCBXnsodDf/JcX/jsIvWhfw18D/pbBdDMANEfE54F+nx9cnqrEa312bx1UV2g75uPanrOIksXq79qXk5tpLupTCh+j96VDoJKdN2BLYKrQfKrz+TiITJCJ+EhFLI2IB8BSF8XciYjAi/iIiLo+IFcAnge70vV+k5yPAk0zQUIvGf9fmcVOltufl2p9KHq79KeXl2kuaAzwL3BgRP03hmlz71J5qtL/i6+8kMkGGVjhImkThM+P/Nn19lqSz0/EfAoMR8UYa3jovxScDf0RhSGy82zkRuzbXddtzdO1PKifX/lT15OLaS/ok8D+AtRHxf4ZOrsW1r2b7M13/SlcB+FHWSomngIPARxTemawCbqOwYuKfgfs4caPnXGAfhYmwlyjspAmF1St7gH+iMPH1N6SVQ+Pc9i9R6H7/E/BqevxbCp80uY1CL2kbML2ozJ0Uelb7KFqJQmGBwOvpew8MveZ6b3vOrv3bFBZxfJh+1+bn6Nr/Rtvzcu0pvBH8ddG5rwIzanHtq9n+LNffd6ybmVlmHs4yM7PMnETMzCwzJxEzM8vMScTMzDJzEjEzs8ycRMzMLDMnETMzy8xJxMzMMvv/r4lGV8I8mfkAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "yearly_incidence.plot(style='*')"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "2020 221186\n",
+ "2023 366227\n",
+ "2021 376290\n",
+ "2024 479258\n",
+ "2002 516689\n",
+ "2018 542312\n",
+ "2017 551041\n",
+ "1996 564901\n",
+ "2019 584066\n",
+ "2015 604382\n",
+ "2000 617597\n",
+ "2001 619041\n",
+ "2012 624573\n",
+ "2005 628464\n",
+ "2006 632833\n",
+ "2022 641397\n",
+ "2011 642368\n",
+ "1993 643387\n",
+ "1995 652478\n",
+ "1994 661409\n",
+ "1998 677775\n",
+ "1997 683434\n",
+ "2014 685769\n",
+ "2013 698332\n",
+ "2007 717352\n",
+ "2008 749478\n",
+ "1999 756456\n",
+ "2003 758363\n",
+ "2004 777388\n",
+ "2016 782114\n",
+ "2010 829911\n",
+ "1992 832939\n",
+ "2009 842373\n",
+ "dtype: int64"
+ ]
+ },
+ "execution_count": 14,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "yearly_incidence.sort_values()"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,