diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb
index 59d72b5b58a3ae26346460dd39e62a39c55243d7..9364c5d20521f65285dc2d79d9b6b022571e0c03 100644
--- a/module3/exo1/analyse-syndrome-grippal.ipynb
+++ b/module3/exo1/analyse-syndrome-grippal.ipynb
@@ -9,7 +9,7 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
@@ -28,13 +28,11 @@
},
{
"cell_type": "code",
- "execution_count": null,
- "metadata": {
- "collapsed": true
- },
+ "execution_count": 5,
+ "metadata": {},
"outputs": [],
"source": [
- "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\""
+ "data_url = \"https://www.sentiweb.fr/datasets/incidence-PAY-3.csv\""
]
},
{
@@ -61,9 +59,976 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 6,
"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",
+ " 202011 | \n",
+ " 3 | \n",
+ " 103527 | \n",
+ " 95314.0 | \n",
+ " 111740.0 | \n",
+ " 157 | \n",
+ " 145.0 | \n",
+ " 169.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1 | \n",
+ " 202010 | \n",
+ " 3 | \n",
+ " 104977 | \n",
+ " 96650.0 | \n",
+ " 113304.0 | \n",
+ " 159 | \n",
+ " 146.0 | \n",
+ " 172.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 2 | \n",
+ " 202009 | \n",
+ " 3 | \n",
+ " 110696 | \n",
+ " 102066.0 | \n",
+ " 119326.0 | \n",
+ " 168 | \n",
+ " 155.0 | \n",
+ " 181.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 3 | \n",
+ " 202008 | \n",
+ " 3 | \n",
+ " 143753 | \n",
+ " 133984.0 | \n",
+ " 153522.0 | \n",
+ " 218 | \n",
+ " 203.0 | \n",
+ " 233.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 4 | \n",
+ " 202007 | \n",
+ " 3 | \n",
+ " 183610 | \n",
+ " 172812.0 | \n",
+ " 194408.0 | \n",
+ " 279 | \n",
+ " 263.0 | \n",
+ " 295.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 5 | \n",
+ " 202006 | \n",
+ " 3 | \n",
+ " 206669 | \n",
+ " 195481.0 | \n",
+ " 217857.0 | \n",
+ " 314 | \n",
+ " 297.0 | \n",
+ " 331.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 6 | \n",
+ " 202005 | \n",
+ " 3 | \n",
+ " 187957 | \n",
+ " 177445.0 | \n",
+ " 198469.0 | \n",
+ " 285 | \n",
+ " 269.0 | \n",
+ " 301.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 7 | \n",
+ " 202004 | \n",
+ " 3 | \n",
+ " 122331 | \n",
+ " 113492.0 | \n",
+ " 131170.0 | \n",
+ " 186 | \n",
+ " 173.0 | \n",
+ " 199.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 8 | \n",
+ " 202003 | \n",
+ " 3 | \n",
+ " 78413 | \n",
+ " 71330.0 | \n",
+ " 85496.0 | \n",
+ " 119 | \n",
+ " 108.0 | \n",
+ " 130.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 9 | \n",
+ " 202002 | \n",
+ " 3 | \n",
+ " 53614 | \n",
+ " 47654.0 | \n",
+ " 59574.0 | \n",
+ " 81 | \n",
+ " 72.0 | \n",
+ " 90.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 10 | \n",
+ " 202001 | \n",
+ " 3 | \n",
+ " 36850 | \n",
+ " 31608.0 | \n",
+ " 42092.0 | \n",
+ " 56 | \n",
+ " 48.0 | \n",
+ " 64.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 11 | \n",
+ " 201952 | \n",
+ " 3 | \n",
+ " 28135 | \n",
+ " 23220.0 | \n",
+ " 33050.0 | \n",
+ " 43 | \n",
+ " 36.0 | \n",
+ " 50.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 12 | \n",
+ " 201951 | \n",
+ " 3 | \n",
+ " 29786 | \n",
+ " 25042.0 | \n",
+ " 34530.0 | \n",
+ " 45 | \n",
+ " 38.0 | \n",
+ " 52.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 13 | \n",
+ " 201950 | \n",
+ " 3 | \n",
+ " 34223 | \n",
+ " 29156.0 | \n",
+ " 39290.0 | \n",
+ " 52 | \n",
+ " 44.0 | \n",
+ " 60.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 14 | \n",
+ " 201949 | \n",
+ " 3 | \n",
+ " 25662 | \n",
+ " 21414.0 | \n",
+ " 29910.0 | \n",
+ " 39 | \n",
+ " 33.0 | \n",
+ " 45.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 15 | \n",
+ " 201948 | \n",
+ " 3 | \n",
+ " 22367 | \n",
+ " 18055.0 | \n",
+ " 26679.0 | \n",
+ " 34 | \n",
+ " 27.0 | \n",
+ " 41.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 16 | \n",
+ " 201947 | \n",
+ " 3 | \n",
+ " 18669 | \n",
+ " 14759.0 | \n",
+ " 22579.0 | \n",
+ " 28 | \n",
+ " 22.0 | \n",
+ " 34.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 17 | \n",
+ " 201946 | \n",
+ " 3 | \n",
+ " 16030 | \n",
+ " 12567.0 | \n",
+ " 19493.0 | \n",
+ " 24 | \n",
+ " 19.0 | \n",
+ " 29.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 18 | \n",
+ " 201945 | \n",
+ " 3 | \n",
+ " 10138 | \n",
+ " 7160.0 | \n",
+ " 13116.0 | \n",
+ " 15 | \n",
+ " 10.0 | \n",
+ " 20.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 19 | \n",
+ " 201944 | \n",
+ " 3 | \n",
+ " 7822 | \n",
+ " 5010.0 | \n",
+ " 10634.0 | \n",
+ " 12 | \n",
+ " 8.0 | \n",
+ " 16.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 20 | \n",
+ " 201943 | \n",
+ " 3 | \n",
+ " 9487 | \n",
+ " 6448.0 | \n",
+ " 12526.0 | \n",
+ " 14 | \n",
+ " 9.0 | \n",
+ " 19.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 21 | \n",
+ " 201942 | \n",
+ " 3 | \n",
+ " 7747 | \n",
+ " 5243.0 | \n",
+ " 10251.0 | \n",
+ " 12 | \n",
+ " 8.0 | \n",
+ " 16.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 22 | \n",
+ " 201941 | \n",
+ " 3 | \n",
+ " 7122 | \n",
+ " 4720.0 | \n",
+ " 9524.0 | \n",
+ " 11 | \n",
+ " 7.0 | \n",
+ " 15.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 23 | \n",
+ " 201940 | \n",
+ " 3 | \n",
+ " 8505 | \n",
+ " 5784.0 | \n",
+ " 11226.0 | \n",
+ " 13 | \n",
+ " 9.0 | \n",
+ " 17.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 24 | \n",
+ " 201939 | \n",
+ " 3 | \n",
+ " 7091 | \n",
+ " 4462.0 | \n",
+ " 9720.0 | \n",
+ " 11 | \n",
+ " 7.0 | \n",
+ " 15.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 25 | \n",
+ " 201938 | \n",
+ " 3 | \n",
+ " 4897 | \n",
+ " 2891.0 | \n",
+ " 6903.0 | \n",
+ " 7 | \n",
+ " 4.0 | \n",
+ " 10.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 26 | \n",
+ " 201937 | \n",
+ " 3 | \n",
+ " 3172 | \n",
+ " 1367.0 | \n",
+ " 4977.0 | \n",
+ " 5 | \n",
+ " 2.0 | \n",
+ " 8.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 27 | \n",
+ " 201936 | \n",
+ " 3 | \n",
+ " 2295 | \n",
+ " 728.0 | \n",
+ " 3862.0 | \n",
+ " 3 | \n",
+ " 1.0 | \n",
+ " 5.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 28 | \n",
+ " 201935 | \n",
+ " 3 | \n",
+ " 1010 | \n",
+ " 2.0 | \n",
+ " 2018.0 | \n",
+ " 2 | \n",
+ " 0.0 | \n",
+ " 4.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 29 | \n",
+ " 201934 | \n",
+ " 3 | \n",
+ " 1672 | \n",
+ " 279.0 | \n",
+ " 3065.0 | \n",
+ " 3 | \n",
+ " 1.0 | \n",
+ " 5.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ "
\n",
+ " \n",
+ " 1816 | \n",
+ " 198521 | \n",
+ " 3 | \n",
+ " 26096 | \n",
+ " 19621.0 | \n",
+ " 32571.0 | \n",
+ " 47 | \n",
+ " 35.0 | \n",
+ " 59.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1817 | \n",
+ " 198520 | \n",
+ " 3 | \n",
+ " 27896 | \n",
+ " 20885.0 | \n",
+ " 34907.0 | \n",
+ " 51 | \n",
+ " 38.0 | \n",
+ " 64.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1818 | \n",
+ " 198519 | \n",
+ " 3 | \n",
+ " 43154 | \n",
+ " 32821.0 | \n",
+ " 53487.0 | \n",
+ " 78 | \n",
+ " 59.0 | \n",
+ " 97.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1819 | \n",
+ " 198518 | \n",
+ " 3 | \n",
+ " 40555 | \n",
+ " 29935.0 | \n",
+ " 51175.0 | \n",
+ " 74 | \n",
+ " 55.0 | \n",
+ " 93.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1820 | \n",
+ " 198517 | \n",
+ " 3 | \n",
+ " 34053 | \n",
+ " 24366.0 | \n",
+ " 43740.0 | \n",
+ " 62 | \n",
+ " 44.0 | \n",
+ " 80.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1821 | \n",
+ " 198516 | \n",
+ " 3 | \n",
+ " 50362 | \n",
+ " 36451.0 | \n",
+ " 64273.0 | \n",
+ " 91 | \n",
+ " 66.0 | \n",
+ " 116.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1822 | \n",
+ " 198515 | \n",
+ " 3 | \n",
+ " 63881 | \n",
+ " 45538.0 | \n",
+ " 82224.0 | \n",
+ " 116 | \n",
+ " 83.0 | \n",
+ " 149.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1823 | \n",
+ " 198514 | \n",
+ " 3 | \n",
+ " 134545 | \n",
+ " 114400.0 | \n",
+ " 154690.0 | \n",
+ " 244 | \n",
+ " 207.0 | \n",
+ " 281.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1824 | \n",
+ " 198513 | \n",
+ " 3 | \n",
+ " 197206 | \n",
+ " 176080.0 | \n",
+ " 218332.0 | \n",
+ " 357 | \n",
+ " 319.0 | \n",
+ " 395.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1825 | \n",
+ " 198512 | \n",
+ " 3 | \n",
+ " 245240 | \n",
+ " 223304.0 | \n",
+ " 267176.0 | \n",
+ " 445 | \n",
+ " 405.0 | \n",
+ " 485.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1826 | \n",
+ " 198511 | \n",
+ " 3 | \n",
+ " 276205 | \n",
+ " 252399.0 | \n",
+ " 300011.0 | \n",
+ " 501 | \n",
+ " 458.0 | \n",
+ " 544.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1827 | \n",
+ " 198510 | \n",
+ " 3 | \n",
+ " 353231 | \n",
+ " 326279.0 | \n",
+ " 380183.0 | \n",
+ " 640 | \n",
+ " 591.0 | \n",
+ " 689.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1828 | \n",
+ " 198509 | \n",
+ " 3 | \n",
+ " 369895 | \n",
+ " 341109.0 | \n",
+ " 398681.0 | \n",
+ " 670 | \n",
+ " 618.0 | \n",
+ " 722.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1829 | \n",
+ " 198508 | \n",
+ " 3 | \n",
+ " 389886 | \n",
+ " 359529.0 | \n",
+ " 420243.0 | \n",
+ " 707 | \n",
+ " 652.0 | \n",
+ " 762.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1830 | \n",
+ " 198507 | \n",
+ " 3 | \n",
+ " 471852 | \n",
+ " 432599.0 | \n",
+ " 511105.0 | \n",
+ " 855 | \n",
+ " 784.0 | \n",
+ " 926.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1831 | \n",
+ " 198506 | \n",
+ " 3 | \n",
+ " 565825 | \n",
+ " 518011.0 | \n",
+ " 613639.0 | \n",
+ " 1026 | \n",
+ " 939.0 | \n",
+ " 1113.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1832 | \n",
+ " 198505 | \n",
+ " 3 | \n",
+ " 637302 | \n",
+ " 592795.0 | \n",
+ " 681809.0 | \n",
+ " 1155 | \n",
+ " 1074.0 | \n",
+ " 1236.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1833 | \n",
+ " 198504 | \n",
+ " 3 | \n",
+ " 424937 | \n",
+ " 390794.0 | \n",
+ " 459080.0 | \n",
+ " 770 | \n",
+ " 708.0 | \n",
+ " 832.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1834 | \n",
+ " 198503 | \n",
+ " 3 | \n",
+ " 213901 | \n",
+ " 174689.0 | \n",
+ " 253113.0 | \n",
+ " 388 | \n",
+ " 317.0 | \n",
+ " 459.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1835 | \n",
+ " 198502 | \n",
+ " 3 | \n",
+ " 97586 | \n",
+ " 80949.0 | \n",
+ " 114223.0 | \n",
+ " 177 | \n",
+ " 147.0 | \n",
+ " 207.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1836 | \n",
+ " 198501 | \n",
+ " 3 | \n",
+ " 85489 | \n",
+ " 65918.0 | \n",
+ " 105060.0 | \n",
+ " 155 | \n",
+ " 120.0 | \n",
+ " 190.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1837 | \n",
+ " 198452 | \n",
+ " 3 | \n",
+ " 84830 | \n",
+ " 60602.0 | \n",
+ " 109058.0 | \n",
+ " 154 | \n",
+ " 110.0 | \n",
+ " 198.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1838 | \n",
+ " 198451 | \n",
+ " 3 | \n",
+ " 101726 | \n",
+ " 80242.0 | \n",
+ " 123210.0 | \n",
+ " 185 | \n",
+ " 146.0 | \n",
+ " 224.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1839 | \n",
+ " 198450 | \n",
+ " 3 | \n",
+ " 123680 | \n",
+ " 101401.0 | \n",
+ " 145959.0 | \n",
+ " 225 | \n",
+ " 184.0 | \n",
+ " 266.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1840 | \n",
+ " 198449 | \n",
+ " 3 | \n",
+ " 101073 | \n",
+ " 81684.0 | \n",
+ " 120462.0 | \n",
+ " 184 | \n",
+ " 149.0 | \n",
+ " 219.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1841 | \n",
+ " 198448 | \n",
+ " 3 | \n",
+ " 78620 | \n",
+ " 60634.0 | \n",
+ " 96606.0 | \n",
+ " 143 | \n",
+ " 110.0 | \n",
+ " 176.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1842 | \n",
+ " 198447 | \n",
+ " 3 | \n",
+ " 72029 | \n",
+ " 54274.0 | \n",
+ " 89784.0 | \n",
+ " 131 | \n",
+ " 99.0 | \n",
+ " 163.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1843 | \n",
+ " 198446 | \n",
+ " 3 | \n",
+ " 87330 | \n",
+ " 67686.0 | \n",
+ " 106974.0 | \n",
+ " 159 | \n",
+ " 123.0 | \n",
+ " 195.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1844 | \n",
+ " 198445 | \n",
+ " 3 | \n",
+ " 135223 | \n",
+ " 101414.0 | \n",
+ " 169032.0 | \n",
+ " 246 | \n",
+ " 184.0 | \n",
+ " 308.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1845 | \n",
+ " 198444 | \n",
+ " 3 | \n",
+ " 68422 | \n",
+ " 20056.0 | \n",
+ " 116788.0 | \n",
+ " 125 | \n",
+ " 37.0 | \n",
+ " 213.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "
1846 rows × 10 columns
\n",
+ "
"
+ ],
+ "text/plain": [
+ " week indicator inc inc_low inc_up inc100 inc100_low \\\n",
+ "0 202011 3 103527 95314.0 111740.0 157 145.0 \n",
+ "1 202010 3 104977 96650.0 113304.0 159 146.0 \n",
+ "2 202009 3 110696 102066.0 119326.0 168 155.0 \n",
+ "3 202008 3 143753 133984.0 153522.0 218 203.0 \n",
+ "4 202007 3 183610 172812.0 194408.0 279 263.0 \n",
+ "5 202006 3 206669 195481.0 217857.0 314 297.0 \n",
+ "6 202005 3 187957 177445.0 198469.0 285 269.0 \n",
+ "7 202004 3 122331 113492.0 131170.0 186 173.0 \n",
+ "8 202003 3 78413 71330.0 85496.0 119 108.0 \n",
+ "9 202002 3 53614 47654.0 59574.0 81 72.0 \n",
+ "10 202001 3 36850 31608.0 42092.0 56 48.0 \n",
+ "11 201952 3 28135 23220.0 33050.0 43 36.0 \n",
+ "12 201951 3 29786 25042.0 34530.0 45 38.0 \n",
+ "13 201950 3 34223 29156.0 39290.0 52 44.0 \n",
+ "14 201949 3 25662 21414.0 29910.0 39 33.0 \n",
+ "15 201948 3 22367 18055.0 26679.0 34 27.0 \n",
+ "16 201947 3 18669 14759.0 22579.0 28 22.0 \n",
+ "17 201946 3 16030 12567.0 19493.0 24 19.0 \n",
+ "18 201945 3 10138 7160.0 13116.0 15 10.0 \n",
+ "19 201944 3 7822 5010.0 10634.0 12 8.0 \n",
+ "20 201943 3 9487 6448.0 12526.0 14 9.0 \n",
+ "21 201942 3 7747 5243.0 10251.0 12 8.0 \n",
+ "22 201941 3 7122 4720.0 9524.0 11 7.0 \n",
+ "23 201940 3 8505 5784.0 11226.0 13 9.0 \n",
+ "24 201939 3 7091 4462.0 9720.0 11 7.0 \n",
+ "25 201938 3 4897 2891.0 6903.0 7 4.0 \n",
+ "26 201937 3 3172 1367.0 4977.0 5 2.0 \n",
+ "27 201936 3 2295 728.0 3862.0 3 1.0 \n",
+ "28 201935 3 1010 2.0 2018.0 2 0.0 \n",
+ "29 201934 3 1672 279.0 3065.0 3 1.0 \n",
+ "... ... ... ... ... ... ... ... \n",
+ "1816 198521 3 26096 19621.0 32571.0 47 35.0 \n",
+ "1817 198520 3 27896 20885.0 34907.0 51 38.0 \n",
+ "1818 198519 3 43154 32821.0 53487.0 78 59.0 \n",
+ "1819 198518 3 40555 29935.0 51175.0 74 55.0 \n",
+ "1820 198517 3 34053 24366.0 43740.0 62 44.0 \n",
+ "1821 198516 3 50362 36451.0 64273.0 91 66.0 \n",
+ "1822 198515 3 63881 45538.0 82224.0 116 83.0 \n",
+ "1823 198514 3 134545 114400.0 154690.0 244 207.0 \n",
+ "1824 198513 3 197206 176080.0 218332.0 357 319.0 \n",
+ "1825 198512 3 245240 223304.0 267176.0 445 405.0 \n",
+ "1826 198511 3 276205 252399.0 300011.0 501 458.0 \n",
+ "1827 198510 3 353231 326279.0 380183.0 640 591.0 \n",
+ "1828 198509 3 369895 341109.0 398681.0 670 618.0 \n",
+ "1829 198508 3 389886 359529.0 420243.0 707 652.0 \n",
+ "1830 198507 3 471852 432599.0 511105.0 855 784.0 \n",
+ "1831 198506 3 565825 518011.0 613639.0 1026 939.0 \n",
+ "1832 198505 3 637302 592795.0 681809.0 1155 1074.0 \n",
+ "1833 198504 3 424937 390794.0 459080.0 770 708.0 \n",
+ "1834 198503 3 213901 174689.0 253113.0 388 317.0 \n",
+ "1835 198502 3 97586 80949.0 114223.0 177 147.0 \n",
+ "1836 198501 3 85489 65918.0 105060.0 155 120.0 \n",
+ "1837 198452 3 84830 60602.0 109058.0 154 110.0 \n",
+ "1838 198451 3 101726 80242.0 123210.0 185 146.0 \n",
+ "1839 198450 3 123680 101401.0 145959.0 225 184.0 \n",
+ "1840 198449 3 101073 81684.0 120462.0 184 149.0 \n",
+ "1841 198448 3 78620 60634.0 96606.0 143 110.0 \n",
+ "1842 198447 3 72029 54274.0 89784.0 131 99.0 \n",
+ "1843 198446 3 87330 67686.0 106974.0 159 123.0 \n",
+ "1844 198445 3 135223 101414.0 169032.0 246 184.0 \n",
+ "1845 198444 3 68422 20056.0 116788.0 125 37.0 \n",
+ "\n",
+ " inc100_up geo_insee geo_name \n",
+ "0 169.0 FR France \n",
+ "1 172.0 FR France \n",
+ "2 181.0 FR France \n",
+ "3 233.0 FR France \n",
+ "4 295.0 FR France \n",
+ "5 331.0 FR France \n",
+ "6 301.0 FR France \n",
+ "7 199.0 FR France \n",
+ "8 130.0 FR France \n",
+ "9 90.0 FR France \n",
+ "10 64.0 FR France \n",
+ "11 50.0 FR France \n",
+ "12 52.0 FR France \n",
+ "13 60.0 FR France \n",
+ "14 45.0 FR France \n",
+ "15 41.0 FR France \n",
+ "16 34.0 FR France \n",
+ "17 29.0 FR France \n",
+ "18 20.0 FR France \n",
+ "19 16.0 FR France \n",
+ "20 19.0 FR France \n",
+ "21 16.0 FR France \n",
+ "22 15.0 FR France \n",
+ "23 17.0 FR France \n",
+ "24 15.0 FR France \n",
+ "25 10.0 FR France \n",
+ "26 8.0 FR France \n",
+ "27 5.0 FR France \n",
+ "28 4.0 FR France \n",
+ "29 5.0 FR France \n",
+ "... ... ... ... \n",
+ "1816 59.0 FR France \n",
+ "1817 64.0 FR France \n",
+ "1818 97.0 FR France \n",
+ "1819 93.0 FR France \n",
+ "1820 80.0 FR France \n",
+ "1821 116.0 FR France \n",
+ "1822 149.0 FR France \n",
+ "1823 281.0 FR France \n",
+ "1824 395.0 FR France \n",
+ "1825 485.0 FR France \n",
+ "1826 544.0 FR France \n",
+ "1827 689.0 FR France \n",
+ "1828 722.0 FR France \n",
+ "1829 762.0 FR France \n",
+ "1830 926.0 FR France \n",
+ "1831 1113.0 FR France \n",
+ "1832 1236.0 FR France \n",
+ "1833 832.0 FR France \n",
+ "1834 459.0 FR France \n",
+ "1835 207.0 FR France \n",
+ "1836 190.0 FR France \n",
+ "1837 198.0 FR France \n",
+ "1838 224.0 FR France \n",
+ "1839 266.0 FR France \n",
+ "1840 219.0 FR France \n",
+ "1841 176.0 FR France \n",
+ "1842 163.0 FR France \n",
+ "1843 195.0 FR France \n",
+ "1844 308.0 FR France \n",
+ "1845 213.0 FR France \n",
+ "\n",
+ "[1846 rows x 10 columns]"
+ ]
+ },
+ "execution_count": 6,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"raw_data = pd.read_csv(data_url, skiprows=1)\n",
"raw_data"
@@ -78,9 +1043,73 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 7,
"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",
+ " 1609 | \n",
+ " 198919 | \n",
+ " 3 | \n",
+ " 0 | \n",
+ " NaN | \n",
+ " NaN | \n",
+ " 0 | \n",
+ " NaN | \n",
+ " NaN | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "
"
+ ],
+ "text/plain": [
+ " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n",
+ "1609 198919 3 0 NaN NaN 0 NaN NaN \n",
+ "\n",
+ " geo_insee geo_name \n",
+ "1609 FR France "
+ ]
+ },
+ "execution_count": 7,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"raw_data[raw_data.isnull().any(axis=1)]"
]
@@ -94,9 +1123,976 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 8,
"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",
+ " 202011 | \n",
+ " 3 | \n",
+ " 103527 | \n",
+ " 95314.0 | \n",
+ " 111740.0 | \n",
+ " 157 | \n",
+ " 145.0 | \n",
+ " 169.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1 | \n",
+ " 202010 | \n",
+ " 3 | \n",
+ " 104977 | \n",
+ " 96650.0 | \n",
+ " 113304.0 | \n",
+ " 159 | \n",
+ " 146.0 | \n",
+ " 172.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 2 | \n",
+ " 202009 | \n",
+ " 3 | \n",
+ " 110696 | \n",
+ " 102066.0 | \n",
+ " 119326.0 | \n",
+ " 168 | \n",
+ " 155.0 | \n",
+ " 181.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 3 | \n",
+ " 202008 | \n",
+ " 3 | \n",
+ " 143753 | \n",
+ " 133984.0 | \n",
+ " 153522.0 | \n",
+ " 218 | \n",
+ " 203.0 | \n",
+ " 233.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 4 | \n",
+ " 202007 | \n",
+ " 3 | \n",
+ " 183610 | \n",
+ " 172812.0 | \n",
+ " 194408.0 | \n",
+ " 279 | \n",
+ " 263.0 | \n",
+ " 295.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 5 | \n",
+ " 202006 | \n",
+ " 3 | \n",
+ " 206669 | \n",
+ " 195481.0 | \n",
+ " 217857.0 | \n",
+ " 314 | \n",
+ " 297.0 | \n",
+ " 331.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 6 | \n",
+ " 202005 | \n",
+ " 3 | \n",
+ " 187957 | \n",
+ " 177445.0 | \n",
+ " 198469.0 | \n",
+ " 285 | \n",
+ " 269.0 | \n",
+ " 301.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 7 | \n",
+ " 202004 | \n",
+ " 3 | \n",
+ " 122331 | \n",
+ " 113492.0 | \n",
+ " 131170.0 | \n",
+ " 186 | \n",
+ " 173.0 | \n",
+ " 199.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 8 | \n",
+ " 202003 | \n",
+ " 3 | \n",
+ " 78413 | \n",
+ " 71330.0 | \n",
+ " 85496.0 | \n",
+ " 119 | \n",
+ " 108.0 | \n",
+ " 130.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 9 | \n",
+ " 202002 | \n",
+ " 3 | \n",
+ " 53614 | \n",
+ " 47654.0 | \n",
+ " 59574.0 | \n",
+ " 81 | \n",
+ " 72.0 | \n",
+ " 90.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 10 | \n",
+ " 202001 | \n",
+ " 3 | \n",
+ " 36850 | \n",
+ " 31608.0 | \n",
+ " 42092.0 | \n",
+ " 56 | \n",
+ " 48.0 | \n",
+ " 64.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 11 | \n",
+ " 201952 | \n",
+ " 3 | \n",
+ " 28135 | \n",
+ " 23220.0 | \n",
+ " 33050.0 | \n",
+ " 43 | \n",
+ " 36.0 | \n",
+ " 50.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 12 | \n",
+ " 201951 | \n",
+ " 3 | \n",
+ " 29786 | \n",
+ " 25042.0 | \n",
+ " 34530.0 | \n",
+ " 45 | \n",
+ " 38.0 | \n",
+ " 52.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 13 | \n",
+ " 201950 | \n",
+ " 3 | \n",
+ " 34223 | \n",
+ " 29156.0 | \n",
+ " 39290.0 | \n",
+ " 52 | \n",
+ " 44.0 | \n",
+ " 60.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 14 | \n",
+ " 201949 | \n",
+ " 3 | \n",
+ " 25662 | \n",
+ " 21414.0 | \n",
+ " 29910.0 | \n",
+ " 39 | \n",
+ " 33.0 | \n",
+ " 45.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 15 | \n",
+ " 201948 | \n",
+ " 3 | \n",
+ " 22367 | \n",
+ " 18055.0 | \n",
+ " 26679.0 | \n",
+ " 34 | \n",
+ " 27.0 | \n",
+ " 41.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 16 | \n",
+ " 201947 | \n",
+ " 3 | \n",
+ " 18669 | \n",
+ " 14759.0 | \n",
+ " 22579.0 | \n",
+ " 28 | \n",
+ " 22.0 | \n",
+ " 34.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 17 | \n",
+ " 201946 | \n",
+ " 3 | \n",
+ " 16030 | \n",
+ " 12567.0 | \n",
+ " 19493.0 | \n",
+ " 24 | \n",
+ " 19.0 | \n",
+ " 29.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 18 | \n",
+ " 201945 | \n",
+ " 3 | \n",
+ " 10138 | \n",
+ " 7160.0 | \n",
+ " 13116.0 | \n",
+ " 15 | \n",
+ " 10.0 | \n",
+ " 20.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 19 | \n",
+ " 201944 | \n",
+ " 3 | \n",
+ " 7822 | \n",
+ " 5010.0 | \n",
+ " 10634.0 | \n",
+ " 12 | \n",
+ " 8.0 | \n",
+ " 16.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 20 | \n",
+ " 201943 | \n",
+ " 3 | \n",
+ " 9487 | \n",
+ " 6448.0 | \n",
+ " 12526.0 | \n",
+ " 14 | \n",
+ " 9.0 | \n",
+ " 19.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 21 | \n",
+ " 201942 | \n",
+ " 3 | \n",
+ " 7747 | \n",
+ " 5243.0 | \n",
+ " 10251.0 | \n",
+ " 12 | \n",
+ " 8.0 | \n",
+ " 16.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 22 | \n",
+ " 201941 | \n",
+ " 3 | \n",
+ " 7122 | \n",
+ " 4720.0 | \n",
+ " 9524.0 | \n",
+ " 11 | \n",
+ " 7.0 | \n",
+ " 15.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 23 | \n",
+ " 201940 | \n",
+ " 3 | \n",
+ " 8505 | \n",
+ " 5784.0 | \n",
+ " 11226.0 | \n",
+ " 13 | \n",
+ " 9.0 | \n",
+ " 17.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 24 | \n",
+ " 201939 | \n",
+ " 3 | \n",
+ " 7091 | \n",
+ " 4462.0 | \n",
+ " 9720.0 | \n",
+ " 11 | \n",
+ " 7.0 | \n",
+ " 15.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 25 | \n",
+ " 201938 | \n",
+ " 3 | \n",
+ " 4897 | \n",
+ " 2891.0 | \n",
+ " 6903.0 | \n",
+ " 7 | \n",
+ " 4.0 | \n",
+ " 10.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 26 | \n",
+ " 201937 | \n",
+ " 3 | \n",
+ " 3172 | \n",
+ " 1367.0 | \n",
+ " 4977.0 | \n",
+ " 5 | \n",
+ " 2.0 | \n",
+ " 8.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 27 | \n",
+ " 201936 | \n",
+ " 3 | \n",
+ " 2295 | \n",
+ " 728.0 | \n",
+ " 3862.0 | \n",
+ " 3 | \n",
+ " 1.0 | \n",
+ " 5.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 28 | \n",
+ " 201935 | \n",
+ " 3 | \n",
+ " 1010 | \n",
+ " 2.0 | \n",
+ " 2018.0 | \n",
+ " 2 | \n",
+ " 0.0 | \n",
+ " 4.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 29 | \n",
+ " 201934 | \n",
+ " 3 | \n",
+ " 1672 | \n",
+ " 279.0 | \n",
+ " 3065.0 | \n",
+ " 3 | \n",
+ " 1.0 | \n",
+ " 5.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ " ... | \n",
+ "
\n",
+ " \n",
+ " 1816 | \n",
+ " 198521 | \n",
+ " 3 | \n",
+ " 26096 | \n",
+ " 19621.0 | \n",
+ " 32571.0 | \n",
+ " 47 | \n",
+ " 35.0 | \n",
+ " 59.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1817 | \n",
+ " 198520 | \n",
+ " 3 | \n",
+ " 27896 | \n",
+ " 20885.0 | \n",
+ " 34907.0 | \n",
+ " 51 | \n",
+ " 38.0 | \n",
+ " 64.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1818 | \n",
+ " 198519 | \n",
+ " 3 | \n",
+ " 43154 | \n",
+ " 32821.0 | \n",
+ " 53487.0 | \n",
+ " 78 | \n",
+ " 59.0 | \n",
+ " 97.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1819 | \n",
+ " 198518 | \n",
+ " 3 | \n",
+ " 40555 | \n",
+ " 29935.0 | \n",
+ " 51175.0 | \n",
+ " 74 | \n",
+ " 55.0 | \n",
+ " 93.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1820 | \n",
+ " 198517 | \n",
+ " 3 | \n",
+ " 34053 | \n",
+ " 24366.0 | \n",
+ " 43740.0 | \n",
+ " 62 | \n",
+ " 44.0 | \n",
+ " 80.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1821 | \n",
+ " 198516 | \n",
+ " 3 | \n",
+ " 50362 | \n",
+ " 36451.0 | \n",
+ " 64273.0 | \n",
+ " 91 | \n",
+ " 66.0 | \n",
+ " 116.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1822 | \n",
+ " 198515 | \n",
+ " 3 | \n",
+ " 63881 | \n",
+ " 45538.0 | \n",
+ " 82224.0 | \n",
+ " 116 | \n",
+ " 83.0 | \n",
+ " 149.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1823 | \n",
+ " 198514 | \n",
+ " 3 | \n",
+ " 134545 | \n",
+ " 114400.0 | \n",
+ " 154690.0 | \n",
+ " 244 | \n",
+ " 207.0 | \n",
+ " 281.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1824 | \n",
+ " 198513 | \n",
+ " 3 | \n",
+ " 197206 | \n",
+ " 176080.0 | \n",
+ " 218332.0 | \n",
+ " 357 | \n",
+ " 319.0 | \n",
+ " 395.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1825 | \n",
+ " 198512 | \n",
+ " 3 | \n",
+ " 245240 | \n",
+ " 223304.0 | \n",
+ " 267176.0 | \n",
+ " 445 | \n",
+ " 405.0 | \n",
+ " 485.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1826 | \n",
+ " 198511 | \n",
+ " 3 | \n",
+ " 276205 | \n",
+ " 252399.0 | \n",
+ " 300011.0 | \n",
+ " 501 | \n",
+ " 458.0 | \n",
+ " 544.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1827 | \n",
+ " 198510 | \n",
+ " 3 | \n",
+ " 353231 | \n",
+ " 326279.0 | \n",
+ " 380183.0 | \n",
+ " 640 | \n",
+ " 591.0 | \n",
+ " 689.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1828 | \n",
+ " 198509 | \n",
+ " 3 | \n",
+ " 369895 | \n",
+ " 341109.0 | \n",
+ " 398681.0 | \n",
+ " 670 | \n",
+ " 618.0 | \n",
+ " 722.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1829 | \n",
+ " 198508 | \n",
+ " 3 | \n",
+ " 389886 | \n",
+ " 359529.0 | \n",
+ " 420243.0 | \n",
+ " 707 | \n",
+ " 652.0 | \n",
+ " 762.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1830 | \n",
+ " 198507 | \n",
+ " 3 | \n",
+ " 471852 | \n",
+ " 432599.0 | \n",
+ " 511105.0 | \n",
+ " 855 | \n",
+ " 784.0 | \n",
+ " 926.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1831 | \n",
+ " 198506 | \n",
+ " 3 | \n",
+ " 565825 | \n",
+ " 518011.0 | \n",
+ " 613639.0 | \n",
+ " 1026 | \n",
+ " 939.0 | \n",
+ " 1113.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1832 | \n",
+ " 198505 | \n",
+ " 3 | \n",
+ " 637302 | \n",
+ " 592795.0 | \n",
+ " 681809.0 | \n",
+ " 1155 | \n",
+ " 1074.0 | \n",
+ " 1236.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1833 | \n",
+ " 198504 | \n",
+ " 3 | \n",
+ " 424937 | \n",
+ " 390794.0 | \n",
+ " 459080.0 | \n",
+ " 770 | \n",
+ " 708.0 | \n",
+ " 832.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1834 | \n",
+ " 198503 | \n",
+ " 3 | \n",
+ " 213901 | \n",
+ " 174689.0 | \n",
+ " 253113.0 | \n",
+ " 388 | \n",
+ " 317.0 | \n",
+ " 459.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1835 | \n",
+ " 198502 | \n",
+ " 3 | \n",
+ " 97586 | \n",
+ " 80949.0 | \n",
+ " 114223.0 | \n",
+ " 177 | \n",
+ " 147.0 | \n",
+ " 207.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1836 | \n",
+ " 198501 | \n",
+ " 3 | \n",
+ " 85489 | \n",
+ " 65918.0 | \n",
+ " 105060.0 | \n",
+ " 155 | \n",
+ " 120.0 | \n",
+ " 190.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1837 | \n",
+ " 198452 | \n",
+ " 3 | \n",
+ " 84830 | \n",
+ " 60602.0 | \n",
+ " 109058.0 | \n",
+ " 154 | \n",
+ " 110.0 | \n",
+ " 198.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1838 | \n",
+ " 198451 | \n",
+ " 3 | \n",
+ " 101726 | \n",
+ " 80242.0 | \n",
+ " 123210.0 | \n",
+ " 185 | \n",
+ " 146.0 | \n",
+ " 224.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1839 | \n",
+ " 198450 | \n",
+ " 3 | \n",
+ " 123680 | \n",
+ " 101401.0 | \n",
+ " 145959.0 | \n",
+ " 225 | \n",
+ " 184.0 | \n",
+ " 266.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1840 | \n",
+ " 198449 | \n",
+ " 3 | \n",
+ " 101073 | \n",
+ " 81684.0 | \n",
+ " 120462.0 | \n",
+ " 184 | \n",
+ " 149.0 | \n",
+ " 219.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1841 | \n",
+ " 198448 | \n",
+ " 3 | \n",
+ " 78620 | \n",
+ " 60634.0 | \n",
+ " 96606.0 | \n",
+ " 143 | \n",
+ " 110.0 | \n",
+ " 176.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1842 | \n",
+ " 198447 | \n",
+ " 3 | \n",
+ " 72029 | \n",
+ " 54274.0 | \n",
+ " 89784.0 | \n",
+ " 131 | \n",
+ " 99.0 | \n",
+ " 163.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1843 | \n",
+ " 198446 | \n",
+ " 3 | \n",
+ " 87330 | \n",
+ " 67686.0 | \n",
+ " 106974.0 | \n",
+ " 159 | \n",
+ " 123.0 | \n",
+ " 195.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1844 | \n",
+ " 198445 | \n",
+ " 3 | \n",
+ " 135223 | \n",
+ " 101414.0 | \n",
+ " 169032.0 | \n",
+ " 246 | \n",
+ " 184.0 | \n",
+ " 308.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ " 1845 | \n",
+ " 198444 | \n",
+ " 3 | \n",
+ " 68422 | \n",
+ " 20056.0 | \n",
+ " 116788.0 | \n",
+ " 125 | \n",
+ " 37.0 | \n",
+ " 213.0 | \n",
+ " FR | \n",
+ " France | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "
1845 rows × 10 columns
\n",
+ "
"
+ ],
+ "text/plain": [
+ " week indicator inc inc_low inc_up inc100 inc100_low \\\n",
+ "0 202011 3 103527 95314.0 111740.0 157 145.0 \n",
+ "1 202010 3 104977 96650.0 113304.0 159 146.0 \n",
+ "2 202009 3 110696 102066.0 119326.0 168 155.0 \n",
+ "3 202008 3 143753 133984.0 153522.0 218 203.0 \n",
+ "4 202007 3 183610 172812.0 194408.0 279 263.0 \n",
+ "5 202006 3 206669 195481.0 217857.0 314 297.0 \n",
+ "6 202005 3 187957 177445.0 198469.0 285 269.0 \n",
+ "7 202004 3 122331 113492.0 131170.0 186 173.0 \n",
+ "8 202003 3 78413 71330.0 85496.0 119 108.0 \n",
+ "9 202002 3 53614 47654.0 59574.0 81 72.0 \n",
+ "10 202001 3 36850 31608.0 42092.0 56 48.0 \n",
+ "11 201952 3 28135 23220.0 33050.0 43 36.0 \n",
+ "12 201951 3 29786 25042.0 34530.0 45 38.0 \n",
+ "13 201950 3 34223 29156.0 39290.0 52 44.0 \n",
+ "14 201949 3 25662 21414.0 29910.0 39 33.0 \n",
+ "15 201948 3 22367 18055.0 26679.0 34 27.0 \n",
+ "16 201947 3 18669 14759.0 22579.0 28 22.0 \n",
+ "17 201946 3 16030 12567.0 19493.0 24 19.0 \n",
+ "18 201945 3 10138 7160.0 13116.0 15 10.0 \n",
+ "19 201944 3 7822 5010.0 10634.0 12 8.0 \n",
+ "20 201943 3 9487 6448.0 12526.0 14 9.0 \n",
+ "21 201942 3 7747 5243.0 10251.0 12 8.0 \n",
+ "22 201941 3 7122 4720.0 9524.0 11 7.0 \n",
+ "23 201940 3 8505 5784.0 11226.0 13 9.0 \n",
+ "24 201939 3 7091 4462.0 9720.0 11 7.0 \n",
+ "25 201938 3 4897 2891.0 6903.0 7 4.0 \n",
+ "26 201937 3 3172 1367.0 4977.0 5 2.0 \n",
+ "27 201936 3 2295 728.0 3862.0 3 1.0 \n",
+ "28 201935 3 1010 2.0 2018.0 2 0.0 \n",
+ "29 201934 3 1672 279.0 3065.0 3 1.0 \n",
+ "... ... ... ... ... ... ... ... \n",
+ "1816 198521 3 26096 19621.0 32571.0 47 35.0 \n",
+ "1817 198520 3 27896 20885.0 34907.0 51 38.0 \n",
+ "1818 198519 3 43154 32821.0 53487.0 78 59.0 \n",
+ "1819 198518 3 40555 29935.0 51175.0 74 55.0 \n",
+ "1820 198517 3 34053 24366.0 43740.0 62 44.0 \n",
+ "1821 198516 3 50362 36451.0 64273.0 91 66.0 \n",
+ "1822 198515 3 63881 45538.0 82224.0 116 83.0 \n",
+ "1823 198514 3 134545 114400.0 154690.0 244 207.0 \n",
+ "1824 198513 3 197206 176080.0 218332.0 357 319.0 \n",
+ "1825 198512 3 245240 223304.0 267176.0 445 405.0 \n",
+ "1826 198511 3 276205 252399.0 300011.0 501 458.0 \n",
+ "1827 198510 3 353231 326279.0 380183.0 640 591.0 \n",
+ "1828 198509 3 369895 341109.0 398681.0 670 618.0 \n",
+ "1829 198508 3 389886 359529.0 420243.0 707 652.0 \n",
+ "1830 198507 3 471852 432599.0 511105.0 855 784.0 \n",
+ "1831 198506 3 565825 518011.0 613639.0 1026 939.0 \n",
+ "1832 198505 3 637302 592795.0 681809.0 1155 1074.0 \n",
+ "1833 198504 3 424937 390794.0 459080.0 770 708.0 \n",
+ "1834 198503 3 213901 174689.0 253113.0 388 317.0 \n",
+ "1835 198502 3 97586 80949.0 114223.0 177 147.0 \n",
+ "1836 198501 3 85489 65918.0 105060.0 155 120.0 \n",
+ "1837 198452 3 84830 60602.0 109058.0 154 110.0 \n",
+ "1838 198451 3 101726 80242.0 123210.0 185 146.0 \n",
+ "1839 198450 3 123680 101401.0 145959.0 225 184.0 \n",
+ "1840 198449 3 101073 81684.0 120462.0 184 149.0 \n",
+ "1841 198448 3 78620 60634.0 96606.0 143 110.0 \n",
+ "1842 198447 3 72029 54274.0 89784.0 131 99.0 \n",
+ "1843 198446 3 87330 67686.0 106974.0 159 123.0 \n",
+ "1844 198445 3 135223 101414.0 169032.0 246 184.0 \n",
+ "1845 198444 3 68422 20056.0 116788.0 125 37.0 \n",
+ "\n",
+ " inc100_up geo_insee geo_name \n",
+ "0 169.0 FR France \n",
+ "1 172.0 FR France \n",
+ "2 181.0 FR France \n",
+ "3 233.0 FR France \n",
+ "4 295.0 FR France \n",
+ "5 331.0 FR France \n",
+ "6 301.0 FR France \n",
+ "7 199.0 FR France \n",
+ "8 130.0 FR France \n",
+ "9 90.0 FR France \n",
+ "10 64.0 FR France \n",
+ "11 50.0 FR France \n",
+ "12 52.0 FR France \n",
+ "13 60.0 FR France \n",
+ "14 45.0 FR France \n",
+ "15 41.0 FR France \n",
+ "16 34.0 FR France \n",
+ "17 29.0 FR France \n",
+ "18 20.0 FR France \n",
+ "19 16.0 FR France \n",
+ "20 19.0 FR France \n",
+ "21 16.0 FR France \n",
+ "22 15.0 FR France \n",
+ "23 17.0 FR France \n",
+ "24 15.0 FR France \n",
+ "25 10.0 FR France \n",
+ "26 8.0 FR France \n",
+ "27 5.0 FR France \n",
+ "28 4.0 FR France \n",
+ "29 5.0 FR France \n",
+ "... ... ... ... \n",
+ "1816 59.0 FR France \n",
+ "1817 64.0 FR France \n",
+ "1818 97.0 FR France \n",
+ "1819 93.0 FR France \n",
+ "1820 80.0 FR France \n",
+ "1821 116.0 FR France \n",
+ "1822 149.0 FR France \n",
+ "1823 281.0 FR France \n",
+ "1824 395.0 FR France \n",
+ "1825 485.0 FR France \n",
+ "1826 544.0 FR France \n",
+ "1827 689.0 FR France \n",
+ "1828 722.0 FR France \n",
+ "1829 762.0 FR France \n",
+ "1830 926.0 FR France \n",
+ "1831 1113.0 FR France \n",
+ "1832 1236.0 FR France \n",
+ "1833 832.0 FR France \n",
+ "1834 459.0 FR France \n",
+ "1835 207.0 FR France \n",
+ "1836 190.0 FR France \n",
+ "1837 198.0 FR France \n",
+ "1838 224.0 FR France \n",
+ "1839 266.0 FR France \n",
+ "1840 219.0 FR France \n",
+ "1841 176.0 FR France \n",
+ "1842 163.0 FR France \n",
+ "1843 195.0 FR France \n",
+ "1844 308.0 FR France \n",
+ "1845 213.0 FR France \n",
+ "\n",
+ "[1845 rows x 10 columns]"
+ ]
+ },
+ "execution_count": 8,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"data = raw_data.dropna().copy()\n",
"data"
@@ -122,7 +2118,7 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
@@ -152,10 +2148,8 @@
},
{
"cell_type": "code",
- "execution_count": null,
- "metadata": {
- "collapsed": true
- },
+ "execution_count": 10,
+ "metadata": {},
"outputs": [],
"source": [
"sorted_data = data.set_index('period').sort_index()"
@@ -179,9 +2173,17 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 11,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n"
+ ]
+ }
+ ],
"source": [
"periods = sorted_data.index\n",
"for p1, p2 in zip(periods[:-1], periods[1:]):\n",
@@ -199,9 +2201,32 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 12,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 12,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXm8HUWZN/59znb33OwhZCEBwhI2gRBABJR9xgVmxBFmFByZQf3xU19nXGAGZV4QB1xHdEQZZdURAVFQZAlh35OwhiwkIQnZc5Ob3Nwldznn1PtHV/XpPqe6u/p0ny33+X4+93P6VtfydHV1PfUs9RQJIcBgMBgMRlQkak0Ag8FgMPYNMENhMBgMRixghsJgMBiMWMAMhcFgMBixgBkKg8FgMGIBMxQGg8FgxAJmKAwGg8GIBcxQGAwGgxELmKEwGAwGIxakak1ANTFx4kQxa9asWpPBYDAYDYUlS5bsEEJMCso3qhjKrFmzsHjx4lqTwWAwGA0FIlpvko9VXgwGg8GIBcxQGAwGgxELmKEwGAwGIxYwQ2EwGAxGLGCGwmAwGIxYEMhQiOhWItpOREsdaeOJaAERrZK/4xz3riKi1US0kojOdaQfT0RvyXs3ERHJ9CYi+p1Mf5mIZjnKXCrbWEVElzrSZ8u8q2TZTPSuYDAYDEYUmEgotwM4ryjtSgALhRBzACyU/4OI5gK4CMARsszPiCgpy9wM4HIAc+SfqvMyALuEEAcD+BGAG2Vd4wFcA+BEAPMBXONgXDcC+JFsf5esg8FgMBg1RCBDEUI8A6C7KPl8AHfI6zsAXOBIv1sIMSSEWAtgNYD5RDQVwBghxIvCOnP4zqIyqq77AJwppZdzASwQQnQLIXYBWADgPHnvDJm3uP19Hrm8wD2LNiCby9eaFAaDwXChXBvKFCHEFgCQv5Nl+jQAGxz5Nsq0afK6ON1VRgiRBdADYIJPXRMA7JZ5i+va5/G/r7yHr//+Tdz5otE+IwaDwaga4jbKkyZN+KSXU8avrlKCiC4nosVEtLirq8srW8Ogq3cIANCzd6TGlDAYDIYb5TKUbVKNBfm7XaZvBDDDkW86gM0yfbom3VWGiFIAOmGp2Lzq2gFgrMxbXFcJhBC3CCHmCSHmTZoUGIqm7pHLW6qudFLHVxkMBqN2KJehPAhAeV1dCuABR/pF0nNrNizj+ytSLdZLRCdJG8glRWVUXRcCeELaWR4FcA4RjZPG+HMAPCrvPSnzFre/zyObt4SxZII9vhkMRn0hMDgkEf0WwAcBTCSijbA8r24AcA8RXQbgPQCfAAAhxNtEdA+AZQCyAK4QQuRkVV+A5THWAuBh+QcAvwJwFxGthiWZXCTr6iai6wAskvmuFUIo54BvALibiL4N4DVZx6hANmcxlFSCJRQGg1FfCGQoQoiLPW6d6ZH/egDXa9IXAzhSkz4IyZA0924FcKsm/V1YrsSjDjkpoaRY5cVgMOoMrDdpMGSlDYUlFAaDUW9ghtJgUNtPZKABBoPBqBswQ2EwGAxGLGCG0mBQmi7PjTcMBoNRIzBDaTDYmi7BLIXBYNQXmKE0GEgGCmB2wmAw6g3MUBoMSkLJ55mlMBiM+gIzlAYD+3YxGIx6BTOUBoNyF2b5hMFg1BuYoTQo2CbPYDDqDcxQGgzEbsMMBqNOwQylwWB7ebGIwmAw6gzMUBoMHHGFwWDUK5ihNBh4XyODwahXMENpMNj7UBqAozy+bBtmXfkQtvYM1poUBoNRBTBDaTA0ktvwb195DwDw1qaeGlPCYDCqAWYoDYZGUnmxvYfBGF1ghtJosN2GG4CjMBiMUQVmKA2GgttwjQkJAXZxZjBGB5ihNBgaS43UUMQyGIyIYIbSYOApmsFg1CuYoTAqDlZ4MRijA8xQGBVDY6nnGAxGVDBDYVQMbItnMEYXmKEwGAwGIxYwQ2lQNIIrLqu8GIzRBWYojIqjAXgfg8GIAcxQGhTUAMv/+qeQwWDECWYoDAaDwYgFzFD2QdyzaAPmX/94HdlZ6oUOBoNRSURiKET0FSJ6m4iWEtFviaiZiMYT0QIiWiV/xznyX0VEq4loJRGd60g/nojekvduIqnPIaImIvqdTH+ZiGY5ylwq21hFRJdGeY5GhB+z+Prv38T23qGa2y4aQCvHYDBiRNkMhYimAfgSgHlCiCMBJAFcBOBKAAuFEHMALJT/g4jmyvtHADgPwM+IKCmruxnA5QDmyL/zZPplAHYJIQ4G8CMAN8q6xgO4BsCJAOYDuMbJuBgWWC5gMBjVRFSVVwpACxGlALQC2AzgfAB3yPt3ALhAXp8P4G4hxJAQYi2A1QDmE9FUAGOEEC8Ka9l9Z1EZVdd9AM6U0su5ABYIIbqFELsALECBCY0K+Bnl1a16UXnVCRkMBqPCKJuhCCE2Afg+gPcAbAHQI4R4DMAUIcQWmWcLgMmyyDQAGxxVbJRp0+R1cbqrjBAiC6AHwASfukYN/JiFYjX5Wqu82M+LwRhViKLyGgdLgpgNYH8AbUT0Kb8imjThk15umWI6LyeixUS0uKury4e8fQ98CBeDwagmoqi8zgKwVgjRJYQYAXA/gPcD2CbVWJC/22X+jQBmOMpPh6Ui2yivi9NdZaRarRNAt09dJRBC3CKEmCeEmDdp0qQyH7UxUS+qpjohg8FgVBhRGMp7AE4iolZp1zgTwHIADwJQXleXAnhAXj8I4CLpuTUblvH9FakW6yWik2Q9lxSVUXVdCOAJaWd5FMA5RDROSkrnyLR9Ho3kOdVItDIYjOhIlVtQCPEyEd0H4FUAWQCvAbgFQDuAe4joMlhM5xMy/9tEdA+AZTL/FUKInKzuCwBuB9AC4GH5BwC/AnAXEa2GJZlcJOvqJqLrACyS+a4VQnSX+yyNBBOpg4gAIZCvFxGFwWCMCpTNUABACHENLPddJ4ZgSSu6/NcDuF6TvhjAkZr0QUiGpLl3K4BbQ5I8qlAv/KRe6GAwGJUF75RvMIRRI9V6HmeVF4MxusAMZR9GvexDYTAYowPMUPZB1Ms+FAV2X2YwRgeYoezL4I2NDAajimCG0qDw02bZoVdqzVEkomrevv3nZbhn8YbgjAwGo6aI5OXFqB385mhLMhC1V3nFJKD88rm1AIC/mzcjICeDwaglWELZh8FGeQaDUU0wQ2lQmPCKemEn9UIHg8GoLJihNChM7CO1FlDYJM9gjC4wQ9kXUWfnoTAYjNEBZigNioZSeTFjYzBGBZihNChMpuhaz+N+p0oyGIx9D8xQGgxhmARHG2YwGNUEM5RGhcERwMxOGAxGNcEMpcEgin5989ZYQmGFF4MxusAMpcEQhkfUWuPFEhKDMbrADKVBYRTLi2d0BoNRRTBDaTCYbGhUUX5rHRySVV4MxugCM5QGRSPslGcwGKMLzFAaDZJJNNbGxlpTwGAwqgFmKA0Gk7lZ2VBqvQ+F9zUyGKMLzFAaFI2wU57BYIwuMENpMKi9JWbMoj44Sq2dAxgMRnXADGUfhNI01frERtZ4MRijC8xQGgxKMmkkL696oYPBYFQWzFAaDKLkwi9vrY3y1ZNRNu3eW7W2GAyGHsxQQmD19l786Y3NtSbDGKNFMnjg9U045YYn8MKaHbUmhcEY1WCGEgJn/fAZfPG3r9WUhoLKyxtKMqi127BCpcl4df0uAMDKrb2VbYjBYPiCGUqDwiSScK35SbUUXuox2QmAwagtmKE0GMxieY1O8AmRDEZtwQylwSDChF6pD41XxV0D6uU5GYzRjkgMhYjGEtF9RLSCiJYT0clENJ6IFhDRKvk7zpH/KiJaTUQriehcR/rxRPSWvHcTyaUmETUR0e9k+stENMtR5lLZxioiujTKc+xzqJPQK9USlZTUxgIKg1FbRJVQfgzgESHEYQCOAbAcwJUAFgoh5gBYKP8HEc0FcBGAIwCcB+BnRJSU9dwM4HIAc+TfeTL9MgC7hBAHA/gRgBtlXeMBXAPgRADzAVzjZFyVRq1PQgQMQ69UnIr6gHodzE8YjNqibIZCRGMAnAbgVwAghBgWQuwGcD6AO2S2OwBcIK/PB3C3EGJICLEWwGoA84loKoAxQogXhTVT31lURtV1H4AzpfRyLoAFQohuIcQuAAtQYEIVR67WW9BhqvKqPZ1A5emwa2cRhcGoKaJIKAcC6AJwGxG9RkS/JKI2AFOEEFsAQP5OlvmnAdjgKL9Rpk2T18XprjJCiCyAHgATfOqqCnI1nKhNJuf6Cb1S3Qme2QmDUVtEYSgpAMcBuFkIcSyAfkj1lgd037vwSS+3jLtRosuJaDERLe7q6vIhzxy1lFCE/ds4Sq8oVDSCezSDwbAQhaFsBLBRCPGy/P8+WAxmm1RjQf5ud+Sf4Sg/HcBmmT5dk+4qQ0QpAJ0Aun3qKoEQ4hYhxDwhxLxJkyaV8ZilaByVV+XpMEIEOsJEVWaNF4NRW5TNUIQQWwFsIKJDZdKZAJYBeBCA8rq6FMAD8vpBABdJz63ZsIzvr0i1WC8RnSTtI5cUlVF1XQjgCWlneRTAOUQ0Thrjz5FpVUE+X62WShGGSdSan6gJPkpMsTDnvlRbxcZgMNxIRSz/RQC/IaIMgHcB/CMsJnUPEV0G4D0AnwAAIcTbRHQPLKaTBXCFECIn6/kCgNsBtAB4WP4BlsH/LiJaDUsyuUjW1U1E1wFYJPNdK4TojvgsxqilDcUEduiVOpCkgGiSkonrs81QmJ8wGDVFJIYihHgdwDzNrTM98l8P4HpN+mIAR2rSByEZkuberQBuDUNvXKitDUUdsGUw0VaaGENEs6GY1C9VXhHaYTAY0cE75ctATRlKGJVXzb28LEShI4y6jCUUNzZ0D2DWlQ/h6XficUZhMILADKUM1IPKyz/asMpTezqBiDaUEM4HbENx49X3rCjM9y7eEJCTwYgHzFDKQC1tE7bbsA8JcUgGccBmbBWmw34do4SfDAxn8eTK7YH5OFgmo9pghhICCfl91oPbsAlqzVAUopBhZJQfZTaUq+5/C/942yKs6eqrNSkMhgvMUEIgIVd8td0pL39NzpSvE5VXFM5mVNT28hodLOXdrn4AQN9g1ih/rUeBEAKPL9uGbK6G/vaMqoAZSggkpIhSWwnFIPSKfWJjpWkJoEPKDJG8vEK1V1kIIfDNPy7Fss17KtySPwo2soB8FafEDE+s2I5/unMxbn5qTa1JYVQYzFBCoJ5UXo0QHLLg4hyhjjpyj+7qG8JdL63HJbe+UqUW9SjYyAyfvMbDdfPuvQCArXsGa0sIo+JghhICtsqrDtyGGyeSVzTU0zM0Ja3TFoZGcgE5Kwwyk/zqRQM4lLVUXZnUvjHdrNi6B7OufAgvvbuz1qTUHfaNN1wlJGxVUj1Ncz6oF5VXFBuKgdpd1R/HBNo3lMXOviF9O7JDB7O1ZSjqMU0XNrW2pSmG0pRKBuRsDCxaawXl+POb2vCBoxrMUEKA6kDlZXIEcCF8fX0wvmg2lOpubDzzB0/h+G8/rr2nXvtIrrb9qp5zJOvPbetlX86+JqG0ZqwAI/1DNZZU6xD7xhuuEpJ1YJQvTLCNE9Y9mg0lPjpMsG2PXjoBam+TUlBkZE0llBqTrfZtpRL1weCioq3JkrT6h8y87EYTmKGEQD3YUMKg1lSaeiP5IUzZSk+cte5PBcXYgtzX68WGsq8hmbCmzXrRANQTmKGEgO3lVQ/7UPxUXvYO9foY8JFsKHXyDEB8E8je4RyeW7Wj7PLKLbye+sYP1YqYwKg9mKGEQCEsfI0JQdDHWR/7UOJAqH0oFd+IEk81V93/Jj71q5exbkd/WeULtjyz/FEn8seXbUNXr7cqMAj7qqDEDLIUzFBCoC4klArmjhtxrEzrSa0QF4NeI3e67xkcKau8qerV3q8SYRwMZXP4pzsX41O/fDk4cwBq7W0WF/ZVBhkHmKGEQNL+kL2Xhi+u2Yk/vLaxYjSECr1SJ99vpImkTp4BiG9CjBpxQS1sgphtHAxdlV27szxpykVIHeDXL63HX/342VqTsc8i6omNowpkMxTvPBf/z0sAgL85dnpFaTGxodSLyivaeSj1g7gYtClD8ALZKs0qBl+J4dnrYYFz9R+XxlZXHTxO3YEllBCoi30oIYZx7VUMMcTyClG48mHy42kgabAw8UPYcRhH/0cZSwXVW/1gJEKgyjoSuOoOzFBCoJ52yhuFXqk9mQCqd2Jjxd2G45JQEtHGkek4jGPiiyMeWxAd2/cMRjL6l4OB4fI3JdbLd1WPYIYSAkpVYbqhrCIwaLpedsrHcXJkmK5ulO/cVnmVa0ORX201vLzCxI4rl5D531mIE67XRyeIG2q3/u6B4bLrsLcWM2cpATOUELBXlnV+YmO9oVrRhiv9gcem8lJG+QrbUOLQzKg2ovRtHMcYxIVZE1oBAO+W6bINFPqiHp6n3sAMJQTqaae836qfIhp96wmhbCiVI8OqPzajfLRxRBElnDAwD/TjjXra2NictsKmRGGQdfD51y2YoYRAXexDCdF2rT/g0Od2REXEZoLojFtCKdvLK/TJoREiFUi1WiQbSvlFY4etDo6wOdmWUJixlIAZSgjYxlCDJUrFJ1E/t+E62ylf6eCQcdmMgvorru5M2KFTyi1v/QaNwzjojdNTsPZehwVEWRTWz1PUH5ihhIBaGZoY5Ss1mZuoIOpN5RXJbdVkA2cM7QDB/RVXd0ZdbBQYqFk70SIVlF9WoZ5UXoghDlq9fFf1CGYoIZCUvWUyoCo16MKs2Gu9lIon9Ip53qhdHsxQ6mMiIUO34TjIjeOZqY42bpgyYz/Ui+Rfj2CGEgJhjKmVXsX4feimE061EEXVEcpmFJGDBjUVV29GnWBNJQ81TKPQHefkWR+j0UKUb4O9vLzBDCUEKAxDqVBE4jDqnXoZ8NUKvVJpCaV+GLRpvuj01pPdIw7EEZaoToZBXYIZSgiEicFUD5NPPdAAVDP0SoWN8jF3Z9T6gorbEkqk1XjZRStaV1SwDaUyYIYSAolQRvlK2VDMDa21Hvd2+5EICaPyioZqGeWjeqWZFovD/hHHM8cRMSEuxOERGAej3lcRmaEQUZKIXiOiP8v/xxPRAiJaJX/HOfJeRUSriWglEZ3rSD+eiN6S924iqVsioiYi+p1Mf5mIZjnKXCrbWEVEl0Z9DhMkQ7gN14OXV70M+Grp8COv+ANtEvH2Z9QxEvR+4yA3jmdWbuy6gVCrMRrHPhRGKeKQUL4MYLnj/ysBLBRCzAGwUP4PIpoL4CIARwA4D8DPiCgpy9wM4HIAc+TfeTL9MgC7hBAHA/gRgBtlXeMBXAPgRADzAVzjZFyVQpiT8uph0NWagjgElGqqvKr9zsptz7SUHTalrFbCteUH8uYnNZOioxnlYyRkH0MkhkJE0wF8GMAvHcnnA7hDXt8B4AJH+t1CiCEhxFoAqwHMJ6KpAMYIIV4U1hd2Z1EZVdd9AM6U0su5ABYIIbqFELsALECBCVUMNkMxGFEVC89iIG5XMzSHH+IIfR7mILHoKi+zdqKiWufVxCKhxECkn09bte0RFHFTKcA2FD9ElVD+C8DXATjX7FOEEFsAQP5OlunTAGxw5Nso06bJ6+J0VxkhRBZAD4AJPnVVFGoc1VLlFQb/8adlGMqWH6Y7LlRPQim/HaD6Xl7l21BM3YbN7W3VgG4RVCvSeKd8ZVA2QyGijwDYLoRYYlpEkyZ80sst426U6HIiWkxEi7u6uowI9YJqwMQoXyn1iX0+hU8ecnTPkyuiPXM0xKByMSgcV2SAQKN8pNrDtxcVcdQep1Feh6pLKDG0W08SylMrt+OpldtrTYaNKBLKKQA+RkTrANwN4Awi+jWAbVKNBfmrnnYjgBmO8tMBbJbp0zXprjJElALQCaDbp64SCCFuEULME0LMmzRpUnlPaldm/Zi5DUdrypMEE6t8ncBWRUWRUKr4oNU2ykdFUN/EsQEvzmfWVVXtLo1D3VgP2geFz9y2CJ+5bVGtybBRNkMRQlwlhJguhJgFy9j+hBDiUwAeBKC8ri4F8IC8fhDARdJzazYs4/srUi3WS0QnSfvIJUVlVF0XyjYEgEcBnENE46Qx/hyZVlGoD9jEPlLLiMR1FOkCQNSd8uZ5Kh96JVr9YduLSkcsoVeiV+FbV62YdLS9OfWlSqwnVGIfyg0AziaiVQDOlv9DCPE2gHsALAPwCIArhBBKwf8FWIb91QDWAHhYpv8KwAQiWg3gXyA9xoQQ3QCuA7BI/l0r0yoKNYDMdspXSOVlYOh28pNymcvgSA5/eG1jPJviqmVDiTj9OV+ZVt8ftw0lYjQF49ArFVbvZHN5bO0Z9LzvV0XVJRSYu/57gRmJN1JxVCKEeArAU/J6J4AzPfJdD+B6TfpiAEdq0gcBfMKjrlsB3FouzeVAjSOTj6zSg86v/jiC8d3w8Arc/sI6TGxvwqlzIqoKI6CaZ8o7J5m8AJJF3RjfGiFarDXTPokl9IpBFdf+eRnufHE93vjWOehsTWvo8K6rVhJKNJUXcxQv8E75EFArvZrulDeYJOLQeO3oGwIAdPdHOXu72kb5CA0VtaVb1cf9TsutztRNOg4GaCLdPPb2NgBA/3DWow7vslW3R8TgwFEIusmMpRjMUEJADZ+hkWBdRaVsKGHtBeUyl4yM1T+Si64aiKQ2C9NOZJWXW0Lxux8HotYXVDyW0CsGedRzqJMoS+sQrl/XvZrZUKKUZUbiBWYoIaDGUf+QfiXmzlvZQec7ecYgomRS1tAYzkYPm1ytTWRxGuX1k1+0+hWiSlTVNMqHCYTqpWn1WwTxTvl9C8xQQkCNIy/RvlgHX0kaKo20LaGUz1BisMlX9eN1G+V19+tMQgno2Tg2NpqUzfswjOCy1Z2d4zlgi728vMAMJQzkCNo7rN997lRzVf6ALe97bi+v8sSVOCSUeNx5w0goUVVIDglFy1AiVe/bXqhypkb5Kkso5Zx4Was9HZEkFPXLDKUEzFBCQI0fL6O80524YgdsGRpko0Lpw+OwBVV6H0o5eXVwSSgamusv2nBQ/d62i7jaAAqSuZc7vR/jrLY9Io44d+zl5Q1mKCEQtA/FxVAqNuiqM5iVXBPHLvdodZgj6gSdD5BQ4p78ypZQQhartcrLbxFUbQnF3ocSQ5+wl1cpmKGEgBpAnhJKnai84mjZjsoaaVdidDrCrCTj9fLSSCgOqTMO5lK2Ud64/jgYurnKy1NCiVh/JRDJbbieYq/UGZihhEBBQtHrs3I554RUWRr8M0Vvp3BIV/S6Ku42bLdTdjMl5fWr6XjebxwBCoHgfo3HhhKcRzESr+fxs6VVe24uSM2VHZOjFcxQQkCNwXqQUEyHdbkexAmbocRgvCy7BrNJMY4giIB7hR00+cXxfsuuwrBcHBvwTN6/yuLJUHz2oVR7tZ8X7t/y6mCW4gVmKCGghpGRDSXihzIwnEXP3hFPGiqt8kpQHLrm6qhc7BwGDe0eGMZnb19kRwJwwm1DKa1LxLRgiCvcfuDGxjhCrxjkUQspr7HiO1adUmE1JuoA5meCKG7S+zqYoYSAHXrFY/e42ygfra0Pfu8pHPN/H/OkwQ9xfJhxGOUVKm6HCbHq/PVL6/HEiu247fm1JfeC9qHkXAzHgK4AlG9DMStYLXWlqQ1F6+iAePs0CPkA5mcE3ofiCWYoZaAaXl7be0tX0E74GzoLKDtOZAxG+ThsG2b8xFulUoxBGTanKZUsrcfJMDRl41Z5Rd/YGFB/PgYJ0UjlKNsLcPPSh7NxXns3tmR9N467bgFWb+8LJsgHBUeF6BIKq75KwQwlBAJtKFVwGy5M0t71x9F0HDuK49gzY7SxLu9uzw/qSOSmVOnQD5rc3Cqv4La8oFxXK+02HMcIDPOcXnuv/Ni0s5/99jw99vY2dPcPRz6dULUQx4mNQXXc+MgKfOa2V8pupxHBDCUECgdseXh5GahEhBB4Yc2OiuqLnSv1ciWURIxuXhVfIYeIajyUVRKKjqH4v78gt+KwiOo2HERCYWNj+QgzToP6RPfZBEUnsO/JX5OziHxpCKEeDaIlqI6bn1qDp1aWdwT3u119+KsfP4vdA+VH+64FmKGEQBgJxWvgP7x0K/7+f17G/77yXiQafFVeMfCqRAwh4f0VSGHr8Mkj3L9+ULHJ0kEMReuR5GgzhkgIlY7lZVcfoZkw799LwvBTiTmT/PpDRW4wOTrCFwH2HhPEoTYLwk+fXI3lW/bg8eX1c168CZihhEDQKslE5aU8t97c0BONFr/VXBwqLyWgBMxGPXtH8PQ7/quwaBJKcOHCvBmcV72jpEZ0c3sclZaNS0KJGkHAdCIruFNHpzVMe1516O661Yzedav3FcQIVm/vxfOrd3jej+MUywKDLLuKQBS8LBvLTsMMJQSCDtgK2scAAG1N1iGZXhGLA2kImZ/K3Ili75QPaPAPr27EZ257BX2akP7xuA0b5AnRjgqerDu7Izj0ij5vWIgYJrVierT35W8cdjATeAWm9jNiu2woPoSaSihn/fAZ/MMvX/a8r5hbHO7wlZzs49gHVgswQwmBOCSU4roC2yyqx2QTn9tbqbwBaXpmx96RPISwzqAvoSMiDYChhBJiglb2Lz1DcdQZYED265fvPboCd764zvN+VD2+aTFT47FvWyGKerVj78kKkPr83l9Kvq/bNe7eYaDUlnHsQ6mkhFJw3KhcG5VALGfKjxrIl5vLCwghSkLD5wxXW2GgO9scCPDyiqFdU08kNUF77c2x6iifDjOjvKIlqC6B51bvBBAsoQS5uPr1y38/uQYAcMnJs7woKWkvDEy95+JQzYSh0Wszr5Iqgm0o3nWr4xT2DJYn2RfaiC6hVMOGkpBLfb8W6jGmGEsoIeB8fTqGYbKxMawCqhzG5Bzn5aq8bJE7IJ86Ilh7EJfhxOcHI4ZiG1r9Ocqjb2/T7pAvrqf4WiGuWF6RJRTDiSwWe0GIvF7Po8aw7raTNBOV19TO5hAUeSMeG4pZHbrn2tk3hDteWOdDR7ANpR7tKyyhhIDz5eeEKOk8tw1F/7JtocZwLNRq0JieG6GeeVjDUKoVvl6RGHR2y9aevY4yGobhiiasoSUmo3zUFa5dKkitGoPKK0xZr/4vBI/0r9+vP1QdibKnPMwVAAAgAElEQVR36rrbi6YGDCfl5IVAsmhh939+9zqeXbUDpxw8EQdPbi8pY6JyzhX1XbmH6cUJllBCIEhCMTkCOOw4Lm7HTAXkElEiIag5pc7wOyq44jYU+RskzSWTheGuE2aC96Ho84ZF2BWuV/lAlZf89RPc1u/s9/WKCvPqvJ7HT+VlKvXFIW1Z5YPbCkLYnfK6fO9s6wUA9GucWYCChsDvgw9aANUCLKGEgPOl6bxNnGleq7Wwk0hxPSar/jgGV8HQ7Z8vm/O2odhlKyyhFFRe/rlTDruJfnJztus/+UXp46h6fFMGnfeZyBU+9P2nkBfAuhs+rK8jBhuKv1HerK1YYnC56olDwjTMr2HovdIW1OthE0rYoY+84bLZCoFE1NVjDGAJJQScH3JOM4G6wtcHeIKFnRT8aCm9Fx326jbgq1FMVKvyMlxJK3zkJ8/i7qINn2FWpEEupU5DfJBKqyEklIDiJpJb0Mo/nJeXPt1PQnGODl86Q343XjDtO9865K+pfVP33GoBlvUQH+3QRz5tVOeE2HBghhICznc2ohkITiZjslfFBGWpvGLQeJmOT0XfSNZH5WVY2dJNe3Dl/W+FpsNedQYxFIeOWSdBViuWV/SNjWb51DP4qSMVdPuIgHCLE293eqv9ICbt91zqfZk+u/cmS0VTlAVBOClHl089jxcdZHB8hEvFHkPkhjjADCUEXC6OmhcYRkIxRYnKy2iFFX21UtgAFiShWB0xolN5hYix5U2HQR6haAlQeSWDVF4OCUVT3nQTXhAKwSy969gzOIKfLFylbUfYv/40qOp176YYQx4LgjArX6/nUc3rHSHMVtkF9Z0ZLV7vJxaVV4hgpM78Tij6gsas3924Y8vFAWYoZUInqjo/Dk+Pl4j6VxP9bXVtKHIVrOmPMOqFKCoXUxVS0mVDKb0fZCPJuYygESYkA3vAdx9ZgR8seAcPL91Sci9s6BWdOrIYXvuIwjymt5eXlFA098xtKO66guA1Ucej8oouoSh4MT7bhuJXtsiGUg9ghhICQgh7UtINBJdRPshAaYhSo3wwnHnKdSU0jQOV9VF5+dhiS+B9OFNwaZXHb3Ml4DbK609khO/9uPahqKJ+E81I1rq3Z68upI37N6gdE5WXV55w0Yb16eq9aPsczj4NnjxNg0NWVEIxWMAUbzHwgjdDUfV40+Ear6zyajwIFCYl3cA2UYmEjZZazIBMJvo4dvAWVv3++Qpuw7qM5h+vV7+YdJepG6dzD4OOsbv2EWnKx7UPxWQfQ0vGOgBsrzakjRoD/rBtKD72LQVPhhJY0tFe4CSuqd8loQTXbapqDJJQ4jkPxaz9ciQUk2OiXdHNWUJpPAgBpOVeBq2E4phUg+J9hTV6O2kIKm86tHb2DeGB1zf51hFEpx16xUflZSKVeX04RsxRZgli1s7bOv5XrX0oJnsq0tLek/XxnjNtx8SG4s3Q/cuaMNmsPeb9mbhfW6oO0wWZrt+AAuOMJGEavD8ng/brQq/nIQO34X3KhkJEM4joSSJaTkRvE9GXZfp4IlpARKvk7zhHmauIaDURrSSicx3pxxPRW/LeTSR7k4iaiOh3Mv1lIprlKHOpbGMVEV1a7nOEgUBB5aVTr+QMXnD0fSjBMG3iC79+FV+++3Vs7Rn0rCNoQleT1bBW5WUuoUTyuoHZ6jUotIp7tewvgUbSwRv0i5+q0pTZq/vDuXzge9S9P5M2XEw6YBGlu+12ZPFux64josrLj7mZwkTCVCpLP1qsex5uwwYSirNovcT1iiKhZAH8qxDicAAnAbiCiOYCuBLAQiHEHAAL5f+Q9y4CcASA8wD8jIjUwd43A7gcwBz5d55MvwzALiHEwQB+BOBGWdd4ANcAOBHAfADXOBlXpeCUUII8Viql8jI5hc/0Y9m02wpFolN3mBoebbdhrZcXZB3BtHjaUEwEFOFfh0KQhBG04otbQjHpF12WsEZ5wFtKUbr6ciUUEybrF8vL1MvL6RVl8vxez2MioXz13jdw4nce97xvomLdHBDmR8HLvGWfh+LHjAxserm8wE+fWIU9gyOe9cSJshmKEGKLEOJVed0LYDmAaQDOB3CHzHYHgAvk9fkA7hZCDAkh1gJYDWA+EU0FMEYI8aKwRsqdRWVUXfcBOFNKL+cCWCCE6BZC7AKwAAUmVDHkBZBJettQnGmeH2jIDVpebsOmei0/k7z6MHWLYVPDr1/oFVMpBwg2ovrVo/IEMZSgDzBoT0Rs+1BCSG56OuSvodsw4G0jUROXt1Henz6X3bAMlZcpk3bbJPxpArzHQoGheFdy35KN2LbHO4ioyb6n7/xleSG/r+Tlv7HRz5/CpO+fWLEd33/sHVz3p2XeFcWIWGwoUhV1LICXAUwRQmwBLKYDYLLMNg3ABkexjTJtmrwuTneVEUJkAfQAmOBTl462y4loMREt7uoq73xnBSGEfXSsbiC4Vrie4r+qy6zNEhuKpq0SOs2qtvPp1Ct2OJMAQpWuWjchmU70fu2YHKssAu7b+QIlEP8FgYnO2kT1YK8JfLL6LgQCW5C0uCQUD4aS8GcowX2qb89Fh63yKr1vEqHbWQfgvbvcCa8FnVJVx+Gl5/f+hkYKNJoySl0bfs/q3tjo9f1Y5dVJsZVGZIZCRO0Afg/g/wgh9vhl1aQJn/Ryy7gThbhFCDFPCDFv0qRJPuQFIy+ErfLS2VBcRnnPCdLbJ1/bZtF4MproDSv3kyBMPhogyMvLQpQVpYmrrq3yCqGe0X2Azn7QGXXdEozH+zVYKRjZlnyOYLbbDmjKedtrL4qKHuDlch21TwGnhOJf3o95mbjkO6Fb8AkhCmFgotjsDN7fjPGt9nVxPmHwzKqM14ZTK09p/mIkE95ORJVAJIZCRGlYzOQ3Qoj7ZfI2qcaC/N0u0zcCmOEoPh3AZpk+XZPuKkNEKQCdALp96qoo8qLgNqyNNmwioRiIy7r8Cupfv/LOO37bUPyM2cZqJAOVVxSjvLPaIE+wYFqddZXed6YFSyhebRhIKML9q4PfiX2i6NcLTnq9jO7JAAklSF3pntT0eQpGeX8JxT98fYE+EztkkEraZEx69Unh1EfvspPHNHm2ZeLIoL7vAZ+jwk2ku5SPir4SiOLlRQB+BWC5EOKHjlsPArhUXl8K4AFH+kXSc2s2LOP7K1It1ktEJ8k6Lykqo+q6EMAT0s7yKIBziGicNMafI9MqCuGUUMq0oagxWu4GLXui91N5GRtt9W1Y7Xjfc8JE5RWJoRiomVRqsFHev66g1bJLvVOmEdvKY543iA5fOPJ5MRS14PCSMIP2RJqoJP1c5U03i7qYgQlD0TyPc4xGYSgmDiu9DiN4MblZA+aomh4YLt2HVKg3uO8L++aqs/MxSvj6UwB8GsBbRPS6TPs3ADcAuIeILgPwHoBPAIAQ4m0iugfAMlgeYlcIIVRvfQHA7QBaADws/wCLYd1FRKthSSYXybq6ieg6AItkvmuFEN0RnsUIQhT2BwSd2BgUvr7cSKUFJuBDp1HN/vYYUzr9VF4FaSqYFpPIAkHqgXA2FF09hWvdh26yIjSTULwn2OI8fh54wdJDMDMOMsrHsQ/FnKF4txVk3/Jq04m9jsnZ5PPzYrKqrF/X9DnC0hfTYrJfTT2vLlJCcR6LlvJpjRNlMxQhxHPwth2e6VHmegDXa9IXAzhSkz4IyZA0924FcKspvXHAaUPxO2Ark0x465Nz4RiKp1HeT+VlOHhUPt0HWrjnzw18VV4hdsoHTUZWHo+y+dK8+roctOnsRgE2FBP1jAnzNFEFmjDqILhVK/o8yYAVbLDbcHDerG03LL2/eXdhD5TfmDaZhN1tahiKI+qAketxgBrQN3SOj3rNyBtUllnT1efZhskCVj1DtRgK75QPgbwAMqlglVc6SZ4fcGgJpZihmKi8DGUUP9uDavfx5duNdgTrVV7qt7wJoLhskJopjMoryAZWtoQSk1He3nfh647t347zdlAQQm+VVxgpyCuP9/3vPboysHwxHfUiofguCHx2yjvvebkNK/q9TnR00gF4L2TUM0Q9R8YUzFBCIFBCEQJEQCqZCBwopjrNsozyQn9dkq+IJiecSUvW7/Kmz0dCgc9E4lWPX3rgJstQ6pnS+0HGUhN7jslCwSQWlOpPvfRouBgxoFediuyt8vJv10SPr8Z6mE2SpXU42jEKJVP6PE4JxeQ9BTkqmDLA4rbc9/TlVV+Y2lACpUOWUOoPThuKbsBm8wKpBCGZoMBNXgbfhJXPwyjvt0oTKEhSfs34SigO+gdHvJmfWgH5nYdipPLyaMIkRLfdJwGdGqSyct7XTSbOicwktE7QIU9+jEE9i/YkzKJfTwQwSKAgoXj1Xd7VZ7r7jua8mKwdbdiP2DA2lODFmO55lQtuKkGRdtsHMdnisn4qL6+Fp0ruH876qFcNGIotoVQHzFBCIC8EUgE2lAQREuSj8lIMxVBCKTHKy1/fFZYA3jdjbGDdqgbdYHQOYsWcdPCTUPxUHcXw1OEHTGjudoJWwM5r7wkHCJZQvJoyUYsJg36xgyFGUHmZGLLDGOV1NTjfm9eYL5wBFEFCMWDmIuB51WbDlnTSSDXpZUMJ6z5ewlBMTnZ1OG547UVxLbY86rHfK0so9Ye8EEj7hK/P2RKKt0rKPtfBUEQpHtMmBkEB4ThPwSefvKWjxUm+86TDYviHXpG0hlADFcNZrXe8r2CprbgNXdahbEG9EGxD0bfldM8NkmL8VTxS5eUj+QXBhN6ErfLyH69edbgCFAaoAYOoNgkOCZhJDjq12HDOer8tmWSgOzRgZlcKkgyK6QKKmXCw9OFlR3H3vTZL1fafKDBDCQEBR3BID6N8IkFIJRKBZ8qbqIGc+W0ahD69GPbGOJ88fgZ+8x3M1qgezpbm8ZOAihG0uvWrx8Su5CyfTJD2mZzMQBsJwUD6+NhPn7OvgwJemkgoI1rGJusJmKJN3mFhp3ywhKjrf+fkGMRQdDTMmdweWB5wjwPPXf0uplP6POr9NqeThiovfZ+YbJB0MZ2i5zbZ9e983nU7+7V5jGwoysuLjfL1g427BrBw+TbLhuLj5ZUXVnj7RCJ4oPnaQHwGiskkLYT/DvmSujTfjZM8vxP/1ErQ9zwUg7FsovLyWsGa2JWcdHjp0IddKi+NCs9gItnj2H8QJKH421Ck95xG3VHOMbhezJYCVV6Fax25Rqt1D4byzT8uxartfWhvSvmWd9ahq0eXrlV5ZR0qrwhG+ZzBmMzm84U9az4qL+99KIXrB1/XBwFxh/7X1zNiaL+KC1E2No4a/NWPn0WvnCjSAUcApxKEJOlXwEBhEPgNaD9vI5MwIwKO0wl9RRRFt//k6aeei03lFUFCsW0oAe2oPkslSMvkhrN5tGaSGBjOebiFO11BTaQufwnFr4qsD6POGk4SJvsUbBuYwSStZSguKcifjuLd+ne9tB4A0JRKoG8oyGsqj1SCkM0Lb+k/QCKzGUomaRTO3Uvl5ZTm/JhoJpnASC5X0m9G6ru8wLSxLdi0ey8mj2nW5jE5ZriwB6g6YAnFAL2OVaeRUd7Hy0ut6P0Ygl+ICFXMl6EIYSSh+OnyTb1q1L2RiCovEwnFq0/tM+UDGIr6AJMJ8jTKt8qjd7UHqOULoSxMpC6vbjOxg/ltbBwJ4dChbGlBq3qvY4Jdh5JppiVnP3l6eck2/AIdKnq9kM0J2znE81kctOj6bdghoZis2L0WUiYqLz96RwxsKLm8QFtTEplkwjX/uPMUrr2epyChVIelMEMJCWWg1g0kZZRPJchTJ20iofj5sBdccb1pFHAeIeovyQDBRnm/SMLq49K6tyrJwWAwl6OLLqTLvIZeXulkQtueJaFYQrveKJ/3ff8l+YMkKhOjvIZ5hImy4HcgHFCYaJwHQrnpcE6epfdN7DRBDEWl+z1WXgg0BTEUl52ltC2nhGKk8grYS+ZHcy5fYCjlenkliNCSSWKvR4BIoz1A7DZc30j7hIN+bvUOdA8MI51MBPr1+62m/XSsaoz7SQ1CwOHl5ZnNd8J3hyHxfhbbU8znfBiTxZGZl5G+rHm0Yeu+1/sZzhUkFJ2dIiecJ3b6NuVLjy0Z+izYFX3Fkp8QwnjV6ZzUgpwedMdAA8GrcZMwIvaiI6vfpKe863xDr+QFmlLq3Xi1U3hI3QJnOKQNxURC8d6HkrfHSnEWl8rMUy1qMZTWTNK1IdMJo30oalOp4SIkKpihhETKw9AGAFt6BjE4Yg0kL520yXkMJp4zngZqWW/SROcloZtsTFRezhWcVuUlk0x8/r1DhQevwpxqQH83aUsVmE6Rh9E2hxbJUHQMLueaJMqXulTT/qpEyVCK8jj7YOMuvVThzJvxUdFa6davl1rF+V6CjPK6STzvc7+tqK+DQtEUwh4FSw666MoqrSmdCFTPWXR5eHm5VNI+9Hr0vXNu8HuWZILQkk567pY32Rel+nbYdCd1RDBDCQlbh148SByDLJ0kT520iUeS26PFfS8olpdKV0H/fCUU2/ag35RYOCujvI9P3TVRD3l61BgY5Z1qPX+1ieWskPZw6x7O5tGcSiJB+gnYUiGZq7y81XjS7uSnSlShV4ryOOleuGI7/OCWqPylJS+GEmQjCZrEFb1ElsrJWUd7s9snyN8ob6DyCmBuQ9kcUgnr/Xt9P85n8GQoRpKBn8rLwIYirNM0LZWXh4TiKOqt8rLa2tHnfaRxnGCGEhIpj/NQBh1iaTqZ8FndBKtnXAwlpFE+X8xQPFvx1+XnRWGF5WUPctKp0zcXvLx8iFDlNRIO4GaoQSt+K7/PKlcaqZMeNq7hbB6ZVAKppJ7hWDYytQ/Js5lAWnzjn0l4ec+F2aiWy+eRTnl7JTrThzzUUW5pWXPfwXB0NhJVf6s0hDvp72hOu/KaSyhek6djPGrGk3q/iYR3JAvX5lYPhm9yrko2522/UguJdFIvKQPWt5MkxKDystK7eofw2nveMfniAjOUkPCSUJzxrvxUXmrl5MtQfHSsrtW4pg71/auzwn1hMyfNLVEIueKpvgv4gMMY5YddqhX9MwbZUAD/s7PzUi+dSia00sFQNo+mVAKpBOltKD6rTm17AW60vgxF0lc8mXgxd307sBcFQTv7vQzmIwESinMhpbORqAVRi3R2cEoAk9qtUw0/8/5ZFo0Bi6wgLy+3DaWUFvV+kwlvlaWzH7zUUQPDOTSng2hxjhU9nc0pb1tOLq+M8ilPlZfJxsaRXB4dzSnc/o8nGIVjigpmKCGRSpClEil6gc6VjZ/KS31QxhJKsVHeuRrXqSBkWspWeXm3UzgCWG9Qt1eEnhKKlU4U/TyUEQ81gEm8IiGAKfLI1fUeu4pVvgSRJ8NQK9hkgjxsKMLRr7r6rURFS1CAUH/vOYu+7r5hbVkT5PP+0bEtGvwZSpC9oF96ILVmklqVl3LlVc4OxXnmHTAOV3zoYIsWP+nSQOXljNbg5TacSSWsWHtGKi99nsGRHMa3ZuS1lw0k72lDUf83pRP+DCVBaEknPFVeJt5m2ZxAayaJDx462fb8rCSYoYREWk44zvnz+dU78IEbnyzk8VF5qQHrZ5D1E6n9XIqd+ZMUrPJS0E1SIzmBlrT3ngygEGyvPZPSG2SF+9eXBg9XSndwSG8134Q2axL3C/dtrfosxwqvndSZVALppP5Dzwl/N1w1KY+Tk00cKq/eoayRGkZfR95XTZTPC3szbi4vtAsH5ySvc0FXk93YlrSWKalxrhhKsQSQTiZsGoc8VDuAYighvLx0RvlcHk2ppAze6iWhOPta0x+5PEZyAuPaFEPxUhUK28Gj+FtXjKoplfSJSyaQJEJrJoWBEb19yzkEPRcM+bytpq0GmKGExIcOnYwEuTfGPfr2Vvv6nz4w23JLDRiweeE9QTpXpcX1OCUffbA+K81E5aVWLNqVZT5vi/Ve/vhqddrZmtYHMbRtKOFUXmEllLyAHb5jwMNnH7AmiFQygXRCz/A37d6LTNJaMGh3yuecRvnS+tVkqSZPT6O8MFd5AUB3f2E8FJfx69uRnECzXBTo8qk+75DGcR1DcEp8usdRDHxsa0Y/jmyVV6mEMpITSCXJljx0ixK7HoeE4vVtOfvGyyhvLRjIe5EUIKGo5x0vGYqXfSOXLyzIivtFMapmHwlFOcVYRvnyHVacY7YaYIYSAh85eir2H9tSElww4RAlTz1kEtLJhPbjemdbL9btHLD/9xpMa3b0e+YZyeftPSZ+EXH9VDMKSR+Gks0HSygqCurY1nSsXl5e+3C85s68sHYVA0DfkPcqd0R6aemCQ27ott7L3Ys2eG5MzQdKKFbbbZK5aTeM5gt7d4JUXuNaLaP1TscCo5huv0l4JJe336FOxVNgKGlJf6mL8tJNewq0a+pQE+q4Nr2EYhvlFUNx2DbUxssmW0Lxk9jyDhtKsL3Hy204k0ygJZ3EYDanXdC5GIrukK4ihuInoShmXuyyq6S+5rS3hJLLWy7uLWnvjY3ub8ObyaoxWw0wQwmBw/brAICSWF0phzTQJFdAugn2yt+/6frfazB9849L7etS92T/VaeaOMwkFOtXy1BywTaUfjl5j2vNIJsXJfSoPjJRebkYisc+HE/GJID9OlsAAFt2e+/NyOYs8T+VLLWR9A05w+voVSLZvLD3IWknIzkhdrZYE/RejarC5RkXIKFMkTGcNjmeqVh94jWhAdZ7bfGTUOR7V9JdsadX8epbL6FkkUwQ2jIpD0lXSijplKTX7ZabTBCICJlUwjc0Sz7vf/y2qk93rTCUzaMpnUBzJul5zoiTqekWBEoCVmpNHRO0Np/m7UVOsT1VfU9NqYTvsdbJhLWxcWBEz/zczjvaarB3pLC3qhpghmKA0w+ZBAD4/OkHAbAma+fklipaAXjZUKaOtSa9U+dMBKCXUIoHTrE9IedY+ejLW7+FjY0+q2CfUwGVi2w6SZ7G0oKEYn1cxSs6Vb+ZhKKXSoxWYXnLk6WjOYWd/cPaPIqeVJKkSrJYr239/70Lj/Y8fsBp5NZ1yertfQCAidJ7Sbe3Qz1DwsORwaY1LzB36hgAwPItBSlB9dMEuUL2m4SHc3lfbyTVvtoPUjw5FqsPdf3fP5RDazqJpnTScxwBTgnFuXAoqGOakglP12Urb962oXhH1g1mKEpCAaA1drvo09ShVF4TfFReWSmFKkm1uF/UdxEkoSRl6BVP5ueMju3xbewdztnzRTXADMUA3/nbo/D01z5oM45ihuGUUI6bOc479IoADprUhiP27wQA3Pni+pIsjyy17DFfO/dQGV21VAQfK1fA/Rr1jvqAgzY25vLCHuh6lZcVtyqVSHhKKGpVv5/0alq7w+1hpfrIyIbiUjU4GYqbZh2Unri9KeV5GJGqN61sJEXvR00uSqXpdcSzn8rrH29fBKCgAu3T0KImkOZ0EiM575392Vwe49oyaEknXc+k+uDEA8cDCJZQ7MWHphnV5x1NehuK6pOTD5wAQD+W9g7n0NqURFMqoaUlW8xQik7FVAbjprReTQxY4ycvEGhDUZN9JpnQSg7D2Tya0kmblgENvU4mo1tIFVR83iqvYsmvWEJR33FzOumpvlNeXq0+zM8p7XmNo8GRnP281QAzFANMG9uCAya02f93NKdcunq1Ojhu5ljb6KdbrfUNZdHWlMIyueJ84PVNrvv5vMAXfvMqAGD2xDbJuJwqEut6Uoc1ge/sL939arskBnx8zg9B750jA11q1EMKaqJ734xxAEpdXMtVebmjxjo/mtJyQihPpQTamlK2o4AOlsqLtBvKBh2HL1k2ML19Ksh1FQDGtFgTiY65qcml1SfEC2BNZqkkoa0piX7nJCf7oy3jbUgHCmoXP/WoqqtglHdPWmqCVittrcprJIfWTApjmtN6iczBQIvpHZHvA7A8noKOug3q+4Ghgn1jUCPtKBtKs88k3esIa69z/Tcxyg873rHOpb4goXg776jFQIsP83O+Ly9hd2A4Z0tk1QAzlDLQ3pRyDbz7X90IALjzshMBeKu8dg8MY1xrBh8+aj8AwJwpHa773QOFCfno6Z0WY3IMasUElITz1qaekjbUyrkprSYsfw8tQK/yyuYEkomEVj2k0NU3hFSCMHWspesv/rhsCSWsUT7vnnTUhPfaht2YdeVDrh2/akLOpCyG4muUz1l6aZ3UpSaX5nQCbZmk1ltscCRXWHX6GNRVnn//w9KSe0pqUROSrm8HhrMYzlqTbUsmiQEHYyou72XIHsrmkRcFZqFTiagJvN3Dy6vAULy91hat7UYqQRjbmkbfUNZzZ7/qE+dCRqkgAYtZeDEUVabdQ5JSUGN6QntGyyyGsjk0pRJ2RGk9Qyn0tW6yVwbyglFe754MwHZBHyphKHkrplzS24aiNmG22LTqxqO/xycgbSjMUOob7U0pPLWyyxYzt+0ZstMB6yPOC/dqBwB29A1jQnsGnzxhJmaOb4XTbN7dP4x5337c/l+tlJ0fqLp/9PROdDSnsEbq7BWEEPbpbmo158VQ5l+/0L72UnmlkyoUv36wbukZxOSOJvu5S3Z1y4/F0in7MxXnKtz5IQ9n87Zx+t7FGwAAC5Zts++r50slCO1NSV+Vl9r3kNK4BavVXks6ifbmlFZdNZTN25Ovn/2jsP9AlKhElHSq9s3oIgxcLZ0y0skE9u9swXvdBc/A7XKszRjfCgDalTgA/PLZdwEAY6QHl96GYqXZXl5Fk6OiXdVRzPy29OzF1j2DWLW9z7bVLFrb7cqjmJLa7LnHEckgm88jKVVemVTCcx+KOgVzUkcTiLzVfP/zzLt2vr2aib5vKGut+pWEoqlHMZTOFr3nYrGE4qfyyiQT1iFbRe94JC+QTni7p6t6m1IJWx2pk/76h7L2d+7VJ4NslK9/vPjuTgDWxHbrc2sBuO0oB0ywPnany2XP3hFs2r3XNuZNbM+49hd86VJP6VYAAByUSURBVLevudpokqsbNaidIUU6mlMY15opCTPyxsYe/GDBOwDg6bIIlOr2dQylf8gaiMVqNwUhBJZu6sF+nc32B1q8qTCbs1wfc3nh6a+vsMJheM4WqbwmSxXfFhli3akTVnnTyQTaMin0eQQ5BKyPqzlteXl52VCa00m0N5XWI4TFHNoyKa1BfdG6wkR64uwJ9gTb1etWS37uriUAYEt1XX2lYePvf3WTbNMy8DuPFd7Wa+U/cJKlgt3jEWrm+49Z46DgqVf6Dn8lx66SYorDldgG6PaM638FZx+t2NILAPjKPa8X1WHl2V86pDidJgaGc3bE4aZUwlY7FqNnwHrGzpY0Wn2i726W42NMc7pkgu3uH8a2PUM4ZEq7vSjQ9V3v4Aha0klLHaXps9Xb+5BOEmaMa/Vkbu9KW+J+nc3IpBIl/WrthyJ7Q6kOQ1LlpWjVLXB+t3iDLa15BffcyyqvxkHvYBbX/nkZAOCyU2fb6fNmWQZTtRrN5wWO+b+PASi4gI5vy7g+rhVbe111N6eTVph1Oaidg/+I/TsxpiVVMoh2OVRmapLQ6YGLw5N0azyjdg0MY0JbRrpzln4063cO4J1tfTh77n72Cqj44xrJ522Pp90D3jG2Hlm6BWu6CjQ53W2Hsnm0NaXsiQeAy2tFGTiVUV734QHAPYs34KV3u9HRnJbBH/Wr8ZZ0Eh3NKfQW1ZN1GIbTyUSJmvC259fa1wdPbsdPLj7O97mVp9/bm/do7wNq34tb7bl9zxA6mlK2TW97r38U2aaU5dXUN1RKx5/esKTZyR0Wc/Py8lLvsNgJpLCpMY2/PW46AOD9B0105VES49TOFkzuaMIy+bzZXB4DwzlbOpo5oQ0rt+7RSrLquN7OlrRnbCv1/gp7N9x5lMPInCnt2H9sqTu2Qu9gFh3NKRmPr/Tb2bR7L6Z2tqAlk0RzKqllKLvk9zRtbIsMw1RkrxuxmIWXhKIWL02phP0dezEMQKngS++v29GP/mGWUOoeKpidM6KBiu0DWC6FHU0pvCcn7sXrCzr/k6THzPi2DLodRvXi8NLWhq9CpFH1Ud108bHYr7MZHU3pknOxtzkOSVKeYDqxfd2OAdf/z63e4fr/5qfWYGA4h3FtGUxoy5TQtntgGB/8/lMAgGNnji14zTjVVrm8vcIG3MzOiVXbevHVe937c7r7HYbRnLWhTXnVAIWQ7EBBukon/Y3yX7/PamP9zn6kE1QilSmVW0umVEIZGM7iD1JqaE4ntWqMYhQ2WurpOWy/MSAq9YxzYkxzGpmUW+15+wvr0DuUtQMrKhWYE85J+chplnp0z96s675zIjxymuWe7OXlNdGWUNzPop7tF586HqccPAFEwPRxLa486l2ObU1j/7Et9jj4hVRPqQnzmOmd2LZnCDv6SseJWkyNaUmjNaOfxC+59RUAwBfPmGPtLi/Ko+oY25qxv1Uds39zUw+mdjZ7empu7Rm01Xe6dgDgL29tAWA5M+jsqcpQnkyQ1oYyMJxDNi/Q2ZK21cnFErMav58//SCMaU6VzAUA8JnbrD7ReYNWCsxQysBXzj4ECQJWbi3YMJIOlRcRYeaEVntX/NodhXyfPukAAMD4tiZ09w9jwbJtnuFCxrakbbWWmhAmyom1o7l0VbJ1T4GhqAlYJ32sk4xO7XMACgxra88gbnxkhUVjawaTxzSVrIJ/8/J79vWM8a1olvsD7n6lkK4GsZpgejxW6uf/9/P2xPTbfz5J0lxob3DEMk46F647HQxOTQpjW9MWQ/GxoQDAmq5+jG3NoNdhQF64fBv++NomdDSlpMorjb0jOVsd8ePHV+HrclNqUzqBdKp0kig2zqqJ4GdPrZZ0DuNDkgm3N6Vw+NQxmDa2Be92lTKUpPREu/jEma4JyckoMqkEJnU04RFH2B8AWLF1D+Z+61H7f8VQeh0Syg8eeweHffMRAMDfzZuOsS3WWCmWytQCQdlrnOMLKNgIO5rTILI2NxYz0K17BkFk2TX2G9OMFVt7IYTA9x5dCQCY2tks+8eq6+o/vlXSH0rlN6Y5LQ+cKn3Hr0jbTV7u0yqe6Hfvtb6DzhZLQm3NJEsm4a09g1i+ZQ+OnTkOqUTp5uQ1XX3Y0D1g2/TGtqSxSzOuH5M2vramlCXhlzCUrDwvXs8cd+8tjGklwRXTqr7rGeNb0NGcLrHXAgXb2MkHTSi5VykwQykDnS1pzDtgPH7+9Bo77ZAij61ZE9ps1dJ6yVgWX32WHT9rYnsGIzmBf75zMa7WeAMBlhTzXvcA8nlhD5gxUvLoaE6XSA7qGNdUgmw7zi+eedcVa+wXT6/B9x5dif07m/GXL5+Kb19wJIDChO88W3x8WwaTO5pL7ABOve/+nc32rnxnWJkX37WkngMnWqqZLR5HzKpJ68CJbXZ4baUKfPndndi0ey/am1Iu9YQzNI1a8Y5rzUhDqrDDqADWJHz9Q8vs///rk+/DpI4mCAHcu9jyzrvsjsVYtb3PthUUSxcpRyyk5pRl1C3+gIsnhgPGW8/97KodWLujH5+7a4ktjXzutAMBAEfu34nF67pdjOKB1zchlxf4xnmHoV2ucJXkoBYVV3/4cACWfWb5lj2u93vtn5bZk+k3PzIXAOSEU5iEf/rkavv6a+cehskdTWjNJF1OHr989l1bnWu5sBM2dBfewUguj8//+lVZf8rut4Gi1fC2nkFMbG9COpnACbPHo6t3yDUJK/XwnCntAEolto27BvDVe98AYLljt2RKbSjfe3SFff1382agJW1FPlaefE+u3I6v/M6qYz/JDMZoJuGT/tNyVDn+gHGY0J7B6u199rvZO5zDmT94Gpt7Bu06JrRnXIubwZEcfvz4Kvv/1nQSnS1p7C6SztfvHJDSh+USXiylOG1G7U2Wza5Ymnp7s+XhObG9CZ2tesY2OJLD+w+agLPnTim5Vyk0NEMhovOIaCURrSaiK6vZ9j+cNNO+nja2BafJ3fQKsyZaEsqS9d1Y3z2AWRNabfUPUDCqApb7byaVsCcahXOO2A8bd+3Faxt22xOC+ngTZHmNzbryIby9uQdCCGzdM4jDp47BK/9+lq0XB4AXHCqt/3zY+viUkVTtablHelApiWhiewYnHzQBkzqa0DuYdemk1+3oRypBePnfzrQZ5OyJbThUMlUhhD3ZHDKlA50taSxcUfDMAoD7lmzErCsfAgCccdhkPPHVD9q63h9Jx4JP3vISAMvg+/NPHYe2TBIzx7dirWNVr3Ty49syOO0QS39/6neftFUCD76xGf/zrGXf+Nq5h+KCY6fZz/xvf3gLT60snHqo3k9HkSHUadTs6hvCfp3NeGbVDvz86TU2sy12kOhsTeOLZxyMZILwoe8/hZcd3k+HyhA+J8wej809gzj1u09i464BXH7nYnz5bsuorZirU+X16gZLdapWyArK0A+4436p92GpvPQS4qSOJiQShIMnt2NNV4GhfPuh5QAs21RLOon9x7Zg464Co/72nwtMWknDbU0p9BVJD8+t3mGrX2dPtBY5L0unlmNmjLXfxceO2R8AcODEdlf56yUdgCXZtRYxlO17BvHfTxYWdjMntNoLA6U++5ffFRwF1J6aMS0p3LN4I+Z9ewEA4PUNu+08R07rxBmHTcaq7X32RL1pd+HZJ0qaJ7Y3uVR0d7/yHn70uDV2j505FokEYVJ7k2tBdu/iDVi2ZQ9OOnCCbXB37jF59b1d+OubngUAdLZkkEwQpna2lNh73tjYY9M6c3wr1hUx4j2DI1i1vc/elFotNCxDIaIkgP8G8FcA5gK4mIjmVqv98983zb7WGb0OmmR9GB+/+UU89OaWkgnnyGmd9vWq7X0YzuZx0OR2l/FZGW4/fvML+Fe5SpsgJ71jZ46z8334pudwyzPv4qmVXZg2tsV2aXz4y6cCsAx6z67qcqkjbvj40QCAk2ZbA+4nT6zGonXd2C7VGn/4/05BR3Ma0yTjudoRX+zltd047ZBJrontsP06sHXPIH76xCrX4P/w0VNx5uGT8craXdjaM4hBGZdowbLCqvo/PnqEq29GcgLbHOqVc46YgvOOnIq3rz0PZxw2Gau7+tA7OIIl67tx/V+sCWd8W8benwMAi6XX1Qurd9pps6QhW01iAPCZ2xbZ1+ceYe0PUqqCnX1DWL+z3/aYAizV1ZL1u9DdP4wbHl6BE65/HEII7OofxozxLXjl38+085504AStF49SISnvtY279uK/Hl9lq0oAYO7+ljqyJZ3E4Ege9yzegDc3FCYRAPjc6YUFyIU3v4A/vLYREzusd9+cTuAIWcf0ca14Y2MPNu/eay8MpoxpwnXnF/q9OZXEs6t2YHAk51L1XHj8dBBZXk1ruvptZvBrqfY8dc5EW73X2ZK2VXibdu/Fb15ej02792KVlHxmS2ahNu+u2lZwRCEifPSY/bFwxTYMjuTw4pqd2Nk35GKQqWQCB01qx9JNPbZq882Nhb1YT/zr6fJ5W2S/DmAom9Ou3tU73tE3jIHhLC747+etNhKE2RPbbFoXSQnyZ08VmNYF8tuf3NGEzbv3QgjLLV4t1qaMacI9nzvZyjOmwFDyeYGvSVve/Nnj7fhzv1u0wa77G/cV7IkHyUXnxI4mlwPPpt17cdPCVZgzuR3TxrbgoEnt2N47ZEtcubzA6d99EgBw1PTCN1ENpIKz1C3mA1gthHgXAIjobgDnA1jmWypGPP4vp+OsHz6Nv5s3veTeeUfuh3+55w37fzVBKEzuaMa15x+Bbz3wtp12zPSxeO4bZ9gqiyljmjGpw73CUR/vxfNnYPmWPbjrJSt8ixrMSl0DAIdPHYPjZo7F/a9twv2vbbJjkn3340fj4MnWB9PZmsa5R0zBo29vwyd+/iIAa1WqJJjzjtwPRzw7Br9/dSN+LzdwAsBnP1DwagMs8b9n7wi+/9g79gT80Jc+gOZ0EnOnjsH9r26yVQpO3HXZfMyU6jkA+P4njsFX730DF/+PJZ3c8LdH4aL5BWnwnLlTcPsL63DSdxba0tFR0zptRvu50w/EL55+F3//y5cxucNt/zlr7mQAhZMCnfjjFafYUoFaDHzsp8+X5Dv/fdMwlM27wubM+/bj2Nk/jEtOPsAlGao+duKWTx+Pw6XtStkPAEtiU3jl38+0N9/Nm2UtHL7umGhmSzXiN849DL942jJuL16/y3b+uHj+DPzn3x5t51dqx/ff8ISdds1Hj8BfHzXV/v9VuVn0sG8+YqtBAeBS6YAytbMZz63egU/e8hIy8ryYw6eOwV1yMy8AHDtjHG59fq0teSr88pJ5AICZ41td6U9+9YOu/8+ZOwV/emOzbd9xQjnCfPSY/XHXS+txxDWPuu63pJM4UL636eOsdi6U4xmwJv+X/63A7A/brwNLZH85bU5//tIHABSkKaf0BwBvfOscdMoo0FM6mzEwnMPsq/6CL505B0PZPE6cPR6/+swJtuPI1M4W7BoYKemTE2aNt73qrvvzMgwMZfHosq02892/sxmT5YItQcAz73ThjB88hc+eMhvfesBa3L1f2kYU47nq/rdw9twp2LR7r81EjztgHKqJhpVQAEwDsMHx/0aZVjUcPLkdr33zbPzzqQeW3GvNpLDgK6fZ///w795XkueSk2fhU1J19t0Lj8ah+3VgXFvGnswB90d308XH2tdEhOsuOBLPfO1Dronr7x2TLwCc6BB5n36nCwBw/Cz3ILvp4mPxyXkzHM/VYTsZNKeT+NZHSgW/C49zM9FPnzSrJM/h+1kT50elOqMYF50wA6fOcasKP3zUVOzf2WyvdD8wx+2GetwB4zBtbAv6h3PoG8pifFsGf/riB+z7XzvnUDsqtGImHz56Ktb+51/bwQWnjW3BP54yy5YAAbiORz3IoY5UuOajc7Huhg/jyGmduPb8I3GtY3WvVo/zZ493lZnc4WZcZx0+GedIKQiwdPVHTXOvIK89/wgXU/rAwe7nH9taOIc9kSAsufosHFN0tOsnHO8SAE4pqgNAyXGwzrGlpNGF/3o6DpPv0Em3cpmeNtbt0fX3J7rHHgBcftqBOEvq8JMJwhvXnIP3HzQBv/mnE0tUdx89Zn9ccvIBJXU89pXT8B8fs/r7hFnjbPugwoePmopXv3m2/X/xfQC453Mnu04sPGHW+JI8937+ZPt5Z00oHQP/cvYhNjMBgI8eXRjXNy20bCdfl7YvhY8cPdVeACi8/q2z0ZxOorM1jb+WUTN+sOAde9/a5acdiEccc4ei6d2uflz9x6UQAL56ziH42nmHASh8I39+cwu+fPfr+O4jKzFtbAuWXH2WvSm1WqCgHcz1CiL6BIBzhRD/JP//NID5QogvFuW7HMDlADBz5szj168vDchYSby4ZicGhrM483C9YSwvgzRGiQiazeXx+PLtOOOwyfZGNoVcXuC259diYnsTfv70GvzzqQfi48eXSlSA5SH04pqduOB901xuuqqNlFyZOs8FKcaG7gE0p5OY2J4pOXL0hdU78PDSrcgLgS+eMQeTpf5e9zy3v7AO82ePx9HTS8/B7hvK4sU1O/Hkyu246IQZJXl69o5ga88gXlizA2Oa057Pm8sLPPTWFpxx2GTXJADADn2yclsvOppT9qq3GI8s3Yr/ePBt/PCTx+DkAyeUPLMQAm9u7MFPnliN7/zNkfaqs5jeF1bvQO9QFn977LSS6NUDw1ms3zmAtzb14OjpnfYEo+uXd7v6tH02OJLDpt17sXD5Npx2yCTPOn78+Cr86PF38J2/OaqEQXT1DmHRum789InVmNrZjG//zZGY2ulmKht3DeDBNzbj/QdNxOBIDsfOHGszclMs3dSD8W0Z/PLZtfjoMVNd6l3AimrQMzCCe5dsxGlzJmnVOlt69mLPXsuTb/bENtt2opDLCzyydCsmtGewfmc/Dp86pqTfcnmBHX1DaE4lsa13ELMmtJV8X4A15l96dye6+obw+dMO0o7pbXsG8fzqHZg/e7xrLK3b0Y+7XlqPtkwSpx0yCRPbm3DAhFbXOMrnBb732EocNKkdY1vSGN+ewXFFfbJqWy+eXbUDa7r6cNCkdpw9d4qtWo0DRLRECDEvMF8DM5STAfyHEOJc+f9VACCE+E+vMvPmzROLFy+uEoUMBoOxb8CUoTSyymsRgDlENJuIMgAuAvBgjWliMBiMUYuGNcoLIbJE9P8DeBRAEsCtQoi3A4oxGAwGo0JoWIYCAEKIvwD4S63pYDAYDEZjq7wYDAaDUUdghsJgMBiMWMAMhcFgMBixgBkKg8FgMGIBMxQGg8FgxIKG3dhYDoioF8BKnyydAHp87gPATADv+dw3qSMoTxx1AI1DaxCdcbUz2voVYFrLycO0luJQIURHQB7YkTJHwx+AxQH3bzGooyuGOnzzxFFHI9EaRGc90dpI/cq0VvR5RhutvnOn+mOVlxt/MsizO+C+SR1BeeKoA2gcWoPojKud0davANNaTh6mtUyMNpXXYmEQj6bSdVQLjUJro9Cp0Ej0Mq2VwWij1bSO0Sah3FIndVQLjUJro9Cp0Ej0Mq2VwWij1aiOUSWhMBgMBqNyGG0SCoPBYDAqhFHPUIjoViLaTkRLHWnHENGLRPQWEf2JiMbI9DQR3SHTl6szWOS9p4hoJRG9Lv8m15jWDBHdJtPfIKIPOsocL9NXE9FNVHwqVH3RWo1+nUFET8p3+jYRfVmmjyeiBUS0Sv6Oc5S5SvbfSiI615Fe0b6NmdaK9m1YWologszfR0Q/Laqrrvo1gNZ669eziWiJ7L8lRHSGo654+9XEFWxf/gNwGoDjACx1pC0CcLq8/iyA6+T13wO4W163AlgHYJb8/ykA8+qI1isA3CavJwNYAiAh/38FwMkACMDDAP6qjmmtRr9OBXCcvO4A8A6AuQC+C+BKmX4lgBvl9VwAbwBoAjAbwBoAyWr0bcy0VrRvy6C1DcAHAHwewE+L6qq3fvWjtd769VgA+8vrIwFsqlS/jnoJRQjxDIDuouRDATwjrxcA+LjKDqCNiFIAWgAMA9hTDTqB0LTOBbBQltsOy3VwHhFNBTBGCPGisEbUnQAuqEda46bJC0KILUKIV+V1L4DlAKYBOB/AHTLbHSj00/mwFhZDQoi1AFYDmF+Nvo2L1jhpiotWIUS/EOI5AIPOeuqxX71orQbKoPU1IcRmmf42gGYiaqpEv456huKBpQA+Jq8/AWCGvL4PQD+ALbB2nn5fCOGcNG+TIu43K6FGCknrGwDOJ6IUEc0GcLy8Nw3ARkf5jTKtHmlVqFq/EtEsWCu6lwFMEUJsAayPGJb0BFj9tcFRTPVhVfs2Iq0KVelbQ1q9UI/9GoR67dePA3hNCDGECvQrMxQ9PgvgCiJaAkukHJbp8wHkAOwPS33wr0R0oLz3D0KIowCcKv8+XWNab4U1QBYD+C8ALwDIwhJti1EtV7+wtAJV7FciagfwewD/RwjhJ3l69WHV+jYGWoEq9W0IWj2r0KTVul/9UJf9SkRHALgRwOdUkiZbpH5lhqKBEGKFEOIcIcTxAH4LS+8MWDaUR4QQI1I18zykakYIsUn+9gL4X1RPraClVQiRFUJ8RQjxPiHE+QDGAlgFa+Ke7qhiOoDNxfXWCa1V61ciSsP6OH8jhLhfJm+TagGldtku0zfCLUGpPqxK38ZEa1X6NiStXqjHfvVEPfYrEU0H8AcAlwgh1HwWe78yQ9FAeWUQUQLA1QB+Lm+9B+AMstAG4CQAK6SqZqIskwbwEVjqnZrRSkStkkYQ0dkAskKIZVIU7iWik6QofgmAB+qR1mr1q+yHXwFYLoT4oePWgwAuldeXotBPDwK4SOqhZwOYA+CVavRtXLRWo2/LoFWLOu1Xr3rqrl+JaCyAhwBcJYR4XmWuSL9GsejvC3+wVspbAIzA4tiXAfgyLM+JdwDcgMIG0HYA98IybC0D8DVR8PhYAuBNee/HkJ40NaR1FqzIyssBPA7gAEc982AN8jUAfqrK1ButVezXD8AS9d8E8Lr8+2sAE2A5C6ySv+MdZf5d9t9KODxjKt23cdFajb4tk9Z1sJw5+uS4mVvH/VpCaz32K6zFW78j7+sAJleiX3mnPIPBYDBiAau8GAwGgxELmKEwGAwGIxYwQ2EwGAxGLGCGwmAwGIxYwAyFwWAwGLGAGQqDUScgos8T0SUh8s8iRzRnBqPWSNWaAAaDYW2IE0L8PDgng1G/YIbCYMQEGajvEViB+o6FtYHzEgCHA/ghrI2xOwB8RgixhYieghW37BQADxJRB4A+IcT3ieh9sCIJtMLadPZZIcQuIjoeVuyzAQDPVe/pGIxgsMqLwYgXhwK4RQhxNKyjDa4A8BMAFworhtmtAK535B8rhDhdCPGDonruBPANWc9bAK6R6bcB+JIQ4uRKPgSDUQ5YQmH8v/buGBfCKArD8PsFjUaltQJLYAGWICJiG5ZBQyMRGtFqRDWdDYhOYwGI0F7FfyeSySAmJxnF+5SnOLl/9eXc/DlXtZ7a176kC+CA4VGj277FfIFhJc3Y5WSDJCsMQTPqpTPgakr9HNiq/wRpNgaKVGtyl9EbcP/DRPH+h96Z0l/6N7zykmqtJRmHxzZwB6yOa0mW+rsU32qtvQLPSTZ7aRcYtdZegNckG72+U398aXZOKFKtB2AvyQnD1tcj4AY47FdWiwyPiN3/0mcPOE6yDDwC+72+D5wm+eh9pX/DbcNSkf6X13VrbX3OR5HmwisvSVIJJxRJUgknFElSCQNFklTCQJEklTBQJEklDBRJUgkDRZJU4hOemcP7BIlCPQAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
"source": [
"sorted_data['inc'].plot()"
]
@@ -215,9 +2240,32 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 13,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 13,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXl8nOV57/29ZkajXbJWb7Kx8QIYCCY4BkK2hgST9DTQBlrzSYPb0uM0IT3J27w9CT19XxJSTkNOmrQkDQkJNEDTAFkhKYQ6ECAsAQxhsw1YBi+yjSRbu0azX+eP537kkTSSRtJoGen6fj7zmdE9z33PPY9m5vdcy33doqoYhmEYRi4EZnsChmEYRuFgomEYhmHkjImGYRiGkTMmGoZhGEbOmGgYhmEYOWOiYRiGYeSMiYZhGIaRMyYahmEYRs6YaBiGYRg5E5rtCeSb+vp6XbVq1WxPwzAMo6B49tlnj6lqw3jHzTvRWLVqFTt37pztaRiGYRQUInIgl+PMPWUYhmHkjImGYRiGkTM5i4aIBEXkdyLyC/d3rYjsEJG97r4m49hrRKRZRF4VkS0Z7eeIyEvuuRtFRFx7sYjc5dqfEpFVGX22udfYKyLb8vGmDcMwjMkxEUvjU8CejL8/BzyoquuAB93fiMgGYCtwOnAx8E0RCbo+NwHbgXXudrFrvwroVNW1wNeAG9xYtcC1wLnAZuDaTHEyDMMwZpacRENEmoDfB76b0XwJcJt7fBtwaUb7naoaU9U3gGZgs4gsBapU9Un1NvG4fVgff6wfARc6K2QLsENVO1S1E9jBCaExDMMwZphcLY1/Bv4nkM5oW6yqRwHcfaNrXw4cyjiuxbUtd4+Htw/po6pJoBuoG2OsIYjIdhHZKSI729vbc3xLhmEYxkQZVzRE5L8Bbar6bI5jSpY2HaN9sn1ONKjerKqbVHVTQ8O4acaGYRjGJMnF0rgA+JCI7AfuBN4rIv8OtDqXE+6+zR3fAqzI6N8EHHHtTVnah/QRkRBQDXSMMZZhGAXIb/a288ax/tmehjEFxhUNVb1GVZtUdRVegPshVf1T4F7Az2baBtzjHt8LbHUZUavxAt5POxdWr4ic5+IVVw7r4491mXsNBR4ALhKRGhcAv8i1GYZRgHz6zue5+dF9sz0NYwpMZUX4l4C7ReQq4CBwOYCq7hKRu4HdQBK4WlVTrs/Hge8BpcD97gZwC3CHiDTjWRhb3VgdIvJF4Bl33HWq2jGFORuGMUuk0kpHJE4knhr/YGPOMiHRUNWHgYfd4+PAhaMcdz1wfZb2ncAZWdqjONHJ8tytwK0TmadhGHOPrkgcVYgmTDQKGVsRbhjGjNAZiQMQTaTHOdKYy5hoGIYxI3T0JwCzNAodEw3DMGaEjn5naSTN0ihkTDQMw5gRfNGImaVR0JhoGIYxI5yIaZhoFDImGoZhzAiD7ikLhBc0JhqGYcwInYMxDbM0ChkTDcMwZoTj/eaemg+YaBiGMSNkrtPwqgQZhYiJhmEYM4If0wCIWdptwWKiYRjGjNDZHycU8HY7iFkwvGAx0TAMY9qJJlL0x1Msrirx/rZgeMFiomEYxrTjxzOWLXKiYcHwgsVEwzCMacePZyytLgVsrUYhY6JhGMa00+mKFS5b5IuGWRqFiomGYRjTzvH+GGDuqfnAuKIhIiUi8rSIvCAiu0TkC6798yJyWESed7cPZvS5RkSaReRVEdmS0X6OiLzknrvRbfuK2xr2Ltf+lIisyuizTUT2uts2DMMoOLoinqWxZDAQbu6pQiWXnftiwHtVtU9EioDHRMTfpvVrqvqVzINFZAPedq2nA8uAX4nIerfl603AduC3wH3AxXhbvl4FdKrqWhHZCtwA/ImI1ALXApsABZ4VkXtVtXNqb9swjJkk5rKlasrDgFkahcy4loZ69Lk/i9xtrOWclwB3qmpMVd8AmoHNIrIUqFLVJ9VbDno7cGlGn9vc4x8BFzorZAuwQ1U7nFDswBMawzAKiLizLKpKigATjUImp5iGiARF5HmgDe9H/Cn31CdF5EURuVVEalzbcuBQRvcW17bcPR7ePqSPqiaBbqBujLGGz2+7iOwUkZ3t7e25vCXDMGaQeDKNCJQXBwFb3FfI5CQaqppS1Y1AE57VcAaeq2kNsBE4CvyTO1yyDTFG+2T7ZM7vZlXdpKqbGhoaxnwvhmHMPLFUmnAwQEmRJxq2uK9wmVD2lKp2AQ8DF6tqqxOTNPAdYLM7rAVYkdGtCTji2puytA/pIyIhoBroGGMswzAKiHhymGiYe6pgySV7qkFEFrnHpcD7gFdcjMLnD4GX3eN7ga0uI2o1sA54WlWPAr0icp6LV1wJ3JPRx8+Mugx4yMU9HgAuEpEa5/66yLUZhlFAxJNpwqEAJSHvJ8cW9xUuuWRPLQVuE5Egnsjcraq/EJE7RGQjnrtoP/AxAFXdJSJ3A7uBJHC1y5wC+DjwPaAUL2vKz8K6BbhDRJrxLIytbqwOEfki8Iw77jpV7ZjC+zUMYxZIpDzRCAUDhAJilkYBM65oqOqLwNlZ2j86Rp/rgeuztO8EzsjSHgUuH2WsW4Fbx5unYRhzF9/SACgpCpqlUcDYinDDMKaduAuEA5QUBSwQXsCYaBiGMe3Ek2mKnGgUh4LmnipgTDQMw5h2YkPcUwFbp1HAmGgYhjHt+IFw8GMaZmkUKiYahmFMO/FkmuJM0bCYRsFiomEYxrQTT52IaZQUBSx7qoAx0TAMY9rxV4SDBcILHRMNwzCmnfiwQLiJRuFiomEYxrSTSOkJ0QjZ4r5CxkTDMIxpJzPltrgoOLgpk1F4mGgYhjHtxJOpoSvCzdIoWEw0DMOYduK2TmPeYKJhGMa0k5k9VRIKkkwryZRZG4WIiYZhGNNKKq2klSHZUwDRpIlGIWKiYRjGtBJ34pDpngLbva9QMdEwDGNa8UUjc0U4mGgUKrls91oiIk+LyAsisktEvuDaa0Vkh4jsdfc1GX2uEZFmEXlVRLZktJ8jIi+55250277itoa9y7U/JSKrMvpsc6+xV0S2YRhGQRFLeeIwmHIb8iyNREpnbU7G5MnF0ogB71XVs4CNwMUich7wOeBBVV0HPOj+RkQ24G3XejpwMfBNt1UswE3Adrx9w9e55wGuAjpVdS3wNeAGN1YtcC1wLrAZuDZTnAzDmPv4lkaxszR8iyNhgfCCZFzRUI8+92eRuylwCXCba78NuNQ9vgS4U1VjqvoG0AxsFpGlQJWqPqmqCtw+rI8/1o+AC50VsgXYoaodqtoJ7OCE0BiGUQD4FoVvaYSCApwQE6OwyCmmISJBEXkeaMP7EX8KWKyqRwHcfaM7fDlwKKN7i2tb7h4Pbx/SR1WTQDdQN8ZYhpF3frO3nZsf3Tfb05h3DA+Eh83SKGhyEg1VTanqRqAJz2o4Y4zDJdsQY7RPts+JFxTZLiI7RWRne3v7GFMzjNH56XOH+fYjr8/2NOYdwwPh/n0ybTGNQmRC2VOq2gU8jOcianUuJ9x9mzusBViR0a0JOOLam7K0D+kjIiGgGugYY6zh87pZVTep6qaGhoaJvCXDGCQST9EfT872NOYd8WGB8CLnnkqYe6ogySV7qkFEFrnHpcD7gFeAewE/m2kbcI97fC+w1WVErcYLeD/tXFi9InKei1dcOayPP9ZlwEMu7vEAcJGI1LgA+EWuzTDyTiSRIppI20rlPBPz3VNBP6bh3cftPBckoRyOWQrc5jKgAsDdqvoLEXkSuFtErgIOApcDqOouEbkb2A0kgatV1U/I/jjwPaAUuN/dAG4B7hCRZjwLY6sbq0NEvgg84467TlU7pvKGDWM0BpyVEUmkqAraEqZ8MXpMw9xThci4oqGqLwJnZ2k/Dlw4Sp/rgeuztO8ERsRDVDWKE50sz90K3DrePA1jqvTHvGubSCxFVUnRLM9m/uCLg79HeFHIuafM0ihI7HLKMBwDboWyxTXyy2iBcBONwsREwzAcEScW/TETjXwyPBBu7qnCxkTDMByRuLM0YlYTKZ8Mj2n4i/vM0ihMTDQMA1BVBpxoRMw9lVfiw7KnzD1V2JhoGAZe+qe/2KzP3FN5JT6sjEiRuacKGhMNw4BBKwNOuKmM/DDc0rAyIoWNiYZhMFQoLBCeX4bHNGxFeGFjomEYDBUNszTySzyVIhgQggFPLPx7szQKExMNw2Coe8osjfwST6YHXVIAIkI4GBiMdRiFhYmGYTB0QZ8t7ssviZQOuqR8ioJiNb4KFBONec7Dr7Zx/X/unu1pzHmGBMJtnUZeiSXThEPBIW1FoYC5pwoUE415TCqtfOHnu/nOb96gpTMy29OZ0/hxjGBAzNLIM/FkerDulE8oYO6pQsVEYx7zqz2tvHGsH4BHXzs2y7OZ2/gL+mrLw7YiPM/EU+nBzCmfcFDM0ihQTDTmMTc/+jpNNaUsqy7hkdfaxu+wgPGLFTZUFJulkWfiydSQQDh47imLaRQmJhrzlAf3tPLsgU7++ztP5t2nNPJ483G7shsD37poqCy2mEaeSaR0sBy6T1EwYCvCCxQTjXlINJHi8z/fxdrGCq7YvJJ3r2+gL5bkuQOdsz21OctAPIkI1JWHrYxInhmecgueaNjOfYVJLtu9rhCRX4vIHhHZJSKfcu2fF5HDIvK8u30wo881ItIsIq+KyJaM9nNE5CX33I1u21fc1rB3ufanRGRVRp9tIrLX3bZhjMu/Pb6fQx0DXHfJ6YRDAd6+tg6AZw+aaIxGJJ6itChIeXHIChbmmXhyZEyjyGIaBUsu270mgc+o6nMiUgk8KyI73HNfU9WvZB4sIhvwtms9HVgG/EpE1rstX28CtgO/Be4DLsbb8vUqoFNV14rIVuAG4E9EpBa4FtgEqHvte1XVfv3G4IVDXaxpKOfta+oBqCopIhwM0Bu1H8PRiCRSlIWDlBUH6bcV4XkllkpTHR66E6LnnjLRKETGtTRU9aiqPuce9wJ7gOVjdLkEuFNVY6r6BtAMbBaRpUCVqj6pqgrcDlya0ec29/hHwIXOCtkC7FDVDicUO/CExhiDtt4oi6tKhrSVFQdtpfMYDMRTlIaDVIRDxJNp+0HLI9ndU2IxjQJlQjEN5zY6G3jKNX1SRF4UkVtFpMa1LQcOZXRrcW3L3ePh7UP6qGoS6AbqxhjLGIO23tgI0SgPhyyVdAz6Y0nKwyHKij3j24Lh+SORShPOGgg3YS5EchYNEakAfgx8WlV78FxNa4CNwFHgn/xDs3TXMdon2ydzbttFZKeI7Gxvbx/zfcx3VJW2nhiNlcVD2svN0hiTgYRnaZSHvZXLlnabP0YLhJtoFCY5iYaIFOEJxvdV9ScAqtqqqilVTQPfATa7w1uAFRndm4Ajrr0pS/uQPiISAqqBjjHGGoKq3qyqm1R1U0NDQy5vad7SPZAgnkrTMEI0QvZDOAaRuBfTKPctDTtXeWPUQHjS3FOFSC7ZUwLcAuxR1a9mtC/NOOwPgZfd43uBrS4jajWwDnhaVY8CvSJynhvzSuCejD5+ZtRlwEMu7vEAcJGI1Dj310WuzRiFtt4YAI1Z3VP2QzgaXvZUiPJiz9LoM/dU3si2IrwoGCCRNkujEMkle+oC4KPASyLyvGv7O+AKEdmI5y7aD3wMQFV3icjdwG68zKurXeYUwMeB7wGleFlT97v2W4A7RKQZz8LY6sbqEJEvAs+4465T1Y7JvdWFQVuPE41hlkZZOMixvthsTKkgiMSTXvZU2I9pmMDmC889NbRgYdjcUwXLuKKhqo+RPbZw3xh9rgeuz9K+EzgjS3sUuHyUsW4Fbh1vnoZHW28UGCkaFeaeGhPfPVXh3FOWdps/4sl09hXh5p4qSGxF+DxjNPeUl3JrP4SjMRBPURYOUeYC4RbTyA/ptBJPpSkZVho9ZIv7ChYTjXlGW09syBWzT3mxxTRGQ1UH3VN+INxKieQHv1RISdGw/TSsjEjBYqIxz2jtjY5wTYEXCI8l01ZZNAuxZJq0Qmk4SFWJt3K5eyAxy7OaH0Rd9eDh+2mEQwGStrivIDHRmGe098RGuKaAwSto89WPxN+AqSwcpDQcpLQoSEdffJZnNT+IJb2LlOIiqz01XzDRmGe0jWppuEVr5nYZgR+/8OMZteVhOvpNNPJBLOHcU8NjGoEAybSSTpu1UWiYaMwz2npjNFaOtDTKbNHaqEQTQ/3udRVhjpto5IVo0rmnika6pwBbq1GAmGjMI/piSSLxFI1VIy2NClu0Niox/4ctZJZGvvEtjeLQ8EC4l4JrRQsLDxONeUS7S7dtqBgpGrZobXTiw/zuJhr5wxfkkhExDe9vS8woPEw05hHH3YrvuorwiOds0droDAZrncuktizM8X5bPZ8PTpzbkSm3gKXdFiAmGvMI3w9fn9XSsED4aIwQjYow0USaARPYKTNayq25pwoXE415xHGXJlpbPtLSOJFya6IxnFhiaEyjzp0/szamzugpty4QnjRLo9Aw0ZhHdLgfuTFFwyyNEYywNMo9S83iGlNnMKYxinsqadlTBYeJxjzieH+c8nBwRMkGgLIi3z1lLpfhxIf53WsHLQ0TjanipzOPZmnErWhhwWGiMY/o6I9TlyWeARAICGVh270vG8NdKL57ylaFT53hrj8ff/tXWxVeeJhozCM6+uNZXVM+ZeGQZU9l4cQ6De/rUOOLhlkaU8YX5OEpt6GAi2mYaBQcJhrziGN98cGr5GxU2D7hWfF/2PxVylUlIYqCYu6pPDB4brPsEQ6WcjsRkqk0L7Z0zfY0ctrudYWI/FpE9ojILhH5lGuvFZEdIrLX3ddk9LlGRJpF5FUR2ZLRfo6IvOSeu9Ft+4rbGvYu1/6UiKzK6LPNvcZeEdmGMSod/bFxLQ0rIzISf9Wy/8MmItSUhek00Zgy0USKUEAIDRMN3z1llW5z557nj/ChbzzOoY7IrM4jF0sjCXxGVU8DzgOuFpENwOeAB1V1HfCg+xv33FbgdOBi4Jsi4js0bwK24+0bvs49D3AV0Kmqa4GvATe4sWqBa4Fzgc3AtZniZJxAVceMaQCU20ZMWYmnRv6w1ZZb/al8EEumsyZmDKbcmqWRM7uP9gCwr71vVucxrmio6lFVfc497gX2AMuBS4Db3GG3AZe6x5cAd6pqTFXfAJqBzSKyFKhS1SdVVYHbh/Xxx/oRcKGzQrYAO1S1Q1U7gR2cEBojg95YkkRKx3RPlduWr1mJJdIjFp/VVYQHU5iNyRNLpkacW7CYxmTY2+aJRSFYGoM4t9HZwFPAYlU9Cp6wAI3usOXAoYxuLa5tuXs8vH1IH1VNAt1A3RhjGcPoGGNhn0952Hbvy0YsmaZ42NVwbXmxBcLzQDSLIMMJ91Tc3FM5s8+JxsFCEQ0RqQB+DHxaVXvGOjRLm47RPtk+mXPbLiI7RWRne3v7GFObv/irl2uz1J3yMfdUdmLJ1IhAbV15mON9cTyj2Jgs47mnrGBhbvTFkhzuGgAKRDREpAhPML6vqj9xza3O5YS7b3PtLcCKjO5NwBHX3pSlfUgfEQkB1UDHGGMNQVVvVtVNqrqpoaEhl7c07/BLiNSXjx7T8FJuzdIYjmdpDP0qNNWU0htLmrUxRWKJ1GBWWiYW05gYvpURDgY42DEwq3PJJXtKgFuAPar61Yyn7gX8bKZtwD0Z7VtdRtRqvID3086F1Ssi57kxrxzWxx/rMuAhF/d4ALhIRGpcAPwi12YMw/9xG9/SSNrV8zDiyZEulDWNFQC8fqx/NqY0b8jm+oPMlFv7LOZCsxONc0+u5VBHZFa/w7lYGhcAHwXeKyLPu9sHgS8B7xeRvcD73d+o6i7gbmA38EvgalX1fSIfB76LFxzfB9zv2m8B6kSkGfgbXCaWqnYAXwSecbfrXJsxDD/TZ6xAeFk4RFqtsuhwYsn0iBXLa+qdaMxypkqhE01kD4QPVrm1goU5sbetj6Kg8M519fTFknRFErM2l9B4B6jqY2SPLQBcOEqf64Hrs7TvBM7I0h4FLh9lrFuBW8eb50Knoz9O2Sh1p3x8v30smd1lsFDJluGzvKaUcCjAvnazNKZCLJmmqrRoRLu5pyZGc1svq+vLWe0uZg52RAYrF8w09ssxTzjeN/bCPjhRWylmV3dDiCVGxjSCAWF1XblZGlMklsX1B5lVbs3qzYXmtj7WNlawsrYMmN1guInGPKEnmmRR2cgrukz8L6+JxlBiyfSI7CmANY3lZmlMkdg47qm4fRZzojOSoLGyhKaaUsBEw8gD/bHk4D7go+H77e2LOpR4lpgGwMn1FRzsiNj5mgKjpdyKCKGAmHsqRyLxJGXhIOXFIeorwrO6wM9EY57QH08O7gM+GicsDVurkUksmRrhngLP0killYMdZm1MltFWhIPnojLRGJ94Mk0ipYNbNi+pLqGtd/aqFZhozBP6Y6nB3flGww9++wX6DI/R/O5rGrygo7moJo9XoiV7ckZRUCyTLwf8IqO+J6GmLDyr64dMNOYJfbEkFcWjZ07BCfeUxTSGki3lFuDkhgoCAj/c2WIuqkkSTaZG7KXhEw6ZpZEL/h445e77XVMWpitiomFMkUguMY0ic09lY7RVyxXFIf7ug6fxqz2tXP0fz9miyAmSSiuJlI5haZho5EIkNtzSKJpVS2PcdRrG3CedVvrj47unfBeMXTUPJZ7K7p4C+Mt3nkxnJM6//nof7X0xGitLZnh2hUs8mX1/cJ+QuadyIjLc0igP0xNNkkylR+xTMhOYpTEPiLh9mM09NXHGuxoGuGBtPQCvvWlrNiZCNDF0G93hFAUDdgGTA/1ZYhoAXQOzsyrcRGMe4Jc7zzkQbu6pQca7GgY4ZXElAK+8OVZxZ2M4J/YHzy7IJaGgfRZzIOIqU5f7ouEW8c7WzpImGvOAQdEYd52GZU8Nx//RGu1qGKCuopj6ijCvtfbO1LTmBeOd25KiAFH7LI7LoKXhPAm1ztLonKX6UyYa8wB/j4ycYxoWfBzEvxoeyz0FsH5xJa+2mntqIox3bkvDQQYSZmmMx2BMw10U+pUfZisYbqIxD+gbdE+NE9NwbgKzNE7gu6fGK+B4ypJK9rb2krZaSTnjxzRGS7ktCQUHjzFGx/cklLrFfX6NudlKuzXRmAf4HypbET5xcnFPgRfXiMRTtHTO7gY4hcR4lkZJkYlGLviWhr8i3A+Ed5hoGJNleHbFaIQCgohlT2Xi+9THE431S7xg+KsW18gZ36IdLcmg2GIaORGJe+uI/MrApeEgJUUBC4Qbk8ePaYxnaYgIxaGAiUYGg1fDY+xDArDO7eT3qmVQ5cyge2q0mIZZGjkRiScpDw89hzVl4bkbCBeRW0WkTURezmj7vIgcHraTn//cNSLSLCKvisiWjPZzROQl99yNbstX3Lawd7n2p0RkVUafbSKy19387WCNYfTnGNMAz1VgufEnGEy5HcfSqCwpYll1CXvbLBieK7Fx0pnNPZUb/bHUCC9CTVl4Tlsa3wMuztL+NVXd6G73AYjIBmArcLrr800R8X/JbgK24+0Zvi5jzKuATlVdC3wNuMGNVQtcC5wLbAaudfuEG8PI1T0FOEvDvqg+/rnIZSfDNY0Vg3s1G+OTU8ptMm3lWcYhEk+OuCCsKS+ic67GNFT1USDXfbkvAe5U1ZiqvoG3F/hmEVkKVKnqk+p9Qm4HLs3oc5t7/CPgQmeFbAF2qGqHqnYCO8guXgue/liS0qIgwcBou/KeoLgoYNlTGcRytDQA1jZWsK+9zzKociSXxX3+inxjdPrjo1gac9U9NQafFJEXnfvKtwCWA4cyjmlxbcvd4+HtQ/qoahLoBurGGMsYRl8OZdF9wkGLaWSS6zoN8EQjmkhzuMsyqHJhvDIifgpp1CzfMYnERloateXhuWtpjMJNwBpgI3AU+CfXnu1SV8don2yfIYjIdhHZKSI729vbx5r3vKQ/h7LoPsVWumEIsXF+2DJZ6/bXaLZ9w3NiPEH2kw8srjE2/fEUpUVDLwoXlYXpHkiQnIWFupMSDVVtVdWUqqaB7+DFHMCzBlZkHNoEHHHtTVnah/QRkRBQjecOG22sbPO5WVU3qeqmhoaGybylgsbzeeZmaRQXmaWRib86fqzaUz5rXQbVPotr5IQvBqPFi0qsrE1OZItp1JYVoQrds1C0cFKi4WIUPn8I+JlV9wJbXUbUaryA99OqehToFZHzXLziSuCejD5+ZtRlwEMu7vEAcJGI1Dj310WuzRhGXyw5bt0pH0u5HcrgWoIc3FN1FcXUlBVZMDxHBuLeVq+jxdr8WIeVEhmbSLaYhl+0cBZcVOP+0ojID4D3APUi0oKX0fQeEdmI5y7aD3wMQFV3icjdwG4gCVytqv4n4uN4mVilwP3uBnALcIeINONZGFvdWB0i8kXgGXfcdaqaa0B+QdEfS1FfEc7p2OJQcNZKKs9FJhIIB8/aaG7r4/HmY6yuL2fZotLpnF5B4/3YjS7GpeaeyolIbOQ6jbryYgCO98VZ2ziz8xlXNFT1iizNt4xx/PXA9VnadwJnZGmPApePMtatwK3jzXGh0x9LclJdWU7HhkOBQT++kZFym+NmNmsbK7jzmUN85LtPceGpjdzyZ2+bzukVNNmukDMpGRQNs3xHI51WIokUZcPcz/WV3kVie19sxudkK8LnAf3x5LirwX2KQwGrcptBLJkmHAwQyCFdGWDjikWowpqGcn6z9xg9UbPaRmMgkRzMkMqGX8jQ3FOjE02mUGWEpdFQ4Vkax3pNNIxJkG3F6GgUh4IWeMwgnhx9q9dsXH7OCn57zYV8+bKziKfSPLSnbRpnV9iM554qMffUuPglgoZbGjVlYYIBMUvDmDiq6iyNHFNuLXtqCLFkKqfMKZ9AQFhSXcLZKxaxpKqE+146Oo2zK2wG4qnBuEU2TDTGJxL3N1gbeh4DAaGuPMyx3pkPhJtoFDiRuDNfJ+CesnUaJ4gl0jnHMzIJBISLz1jCw6+1D9b+MoYykBjP0rCU2/EYtDSyeBIaKovN0jAmjl93KucV4ZZyO4RYMj1uhdvROH9NHfFkmn222C8ruQbCLaYxOpHBunIjP6P1FcW0W0zDmCgntnrNfUV43IrEDRJNpCYU08hkZa2XsXaow8qRi0BCAAAgAElEQVSKZGMgnhq17hRYym0uDG71muX73VBZzDGzNIyJMlgWfQKL+8D2CfeJxHOv2zWcFU40DnZE8jmleUMknswxEG6fxdGIjFHB2heNmb4ANNEocPyUz8qSopyOP7Hlq31RwXPvjfXDNhYVxSFqy8Mc6jTRyMZ42VPBgFAUFCtYOAaDnoQsolFfUUwipTNeSsREo8DpGfCuRKpKc6095X2JLfjoEYmlcrbSsrGippRDZmmMIJVWYsn0mOs0wCuPPhA30RgNP2aZ7Tw2VHprNWY6rmGiUeD0OkujKldLI+hbGvZFBWdp5BgPysaK2jITjSz4we3xrLiSsFVdHoseZ0Vkuyj0SweZaBgToifqLI1cRcOlOdqWrx6R+BQtjdoyDncNkLKNmYYQGbxCHvvclhQFLKYxBt0DCcrCwawFNRt9S2OGg+EmGgWOb2lUlEwsEG4xDY/xgrXjsbK2jERKebMnmsdZFT6+y6lsnHRmc0+NTVckQXVp9gvChooSwCwNY4L0DHh1p3LZ6hVOlAA30fD87tFEOucSLNlYUeMyqI6biyqTnN1TRUELhI9B18DoolFVGiIcDHCsb2ZXhZtoFDi90QRVOVoZkGFpWG78iRINU4ppeKXRLYNqKP76gvEC4aVFQVunMQbdkQSLyrKLhohQXxE2S8OYGD3RRM7ptnBiFzWzNE78sE3F0li2qJSAYMHwYfgup7FqT4EXY7OYxuh0DyRYVDr6XjmLq0t45c2eGV2rYaJR4PQMJHNOt4UT7ikLhGcsjJyCpVEUDLBsUSkHzD01hFwFucQsjTHpGoiP6p4C+ONNK9h1pIdfzWC15XFFQ0RuFZE2EXk5o61WRHaIyF53X5Px3DUi0iwir4rIloz2c0TkJffcjW7bV9zWsHe59qdEZFVGn23uNfaKiL8lrJFBb2xiloafPWWWRn4sDYBTl1Sx60h3PqY0b4iMsb4gExONsekawz0FcPk5TZxcX87/eeCVGbsQzMXS+B5w8bC2zwEPquo64EH3NyKyAW+71tNdn2+KiP+puQnYjrdv+LqMMa8COlV1LfA14AY3Vi3e1rLnApuBazPFyfDoGUhOLqZhwcdBS2Mq2VMAb2mq5vVj/YOZbEZG9tS4MQ1zT41GNJEilkxTPYZohIIB/nbLKbzW2sfvfeVhbn9y/7S7qsYVDVV9FG/v7kwuAW5zj28DLs1ov1NVY6r6BtAMbBaRpUCVqj6p3ju6fVgff6wfARc6K2QLsENVO1S1E9jBSPFa8PROMKZh2VMniOSY4TMeZzZVowq7jvTkY1rzgkiOomHZU6PjlwcZK6YB8IEzl/Jvf/Y2llSX8OCeNpwTZ9qYrF2+WFWPAqjqURHxtzZfDvw247gW15Zwj4e3+30OubGSItIN1GW2Z+lj4G3A1BOdaEzDFvf5RAYrBE/NPXXm8moAXj7czXkn1015XvMBP+U2F/eUrdPITlfEE42xYho+v3dqI+85pWFQrKeTfAfCs0mcjtE+2T5DX1Rku4jsFJGd7e3tOU10PjCQSJFK6ySzp+yL2j/GXgUTob6imGXVJbzYYnENn0g8STAg425wVVIUJGal+rPSFfHWX4wV08hERKZ8AZQLkxWNVudywt37ofsWYEXGcU3AEdfelKV9SB8RCQHVeO6w0cYagarerKqbVHVTQ0PDJN9S4TFYrHBC7inbLc0nMsGy8mNxZlM1Lx020fCJxFOUFQXHdZWUWGLGqHQN5G5pzCSTFY17AT+baRtwT0b7VpcRtRov4P20c2X1ish5Ll5x5bA+/liXAQ+5uMcDwEUiUuMC4Be5NsPhl0WfiHsqFAwQDIh9SYF+3+8+hZRbn7c0LeKNY/0zXqZ6rhJNpCjJwYIrcTE2c1GNpNu5p3K1NGaKcX9tROQHwHuAehFpwcto+hJwt4hcBRwELgdQ1V0icjewG0gCV6uq/2n4OF4mVilwv7sB3ALcISLNeBbGVjdWh4h8EXjGHXedqg4PyC9oeie4l4ZPWThIn+1rnbMLJRdOW1oJQHNbH+ecZEl+4+2l4TO4EZO5S0cwGAgvGzsQPtOMKxqqesUoT104yvHXA9dnad8JnJGlPYoTnSzP3QrcOt4cFyon3FMTc680VMzOhvRzjf6Y98OWj2yTJVVeOZFWK1wIeKIx3mpwgNKwJ9iWdjuSroE4wYBQPsWYW76xFeEFzER37fOpryjm2CxsSD/XGJhiWfRMllZ7FUePdptogHduc7I0zD01Kl2RBItKi6Y9hXaimGgUMIN7aUwgpgFQXxmelQ3p5xpT3YApk0VlRRSHAmZpOLyS8+N/Lv0gb48tjBxB10BizIV9s4WJRgEz0V37fOrKizneP7PllOciU92AKRMRYUl1iVkajkg8Ne4aDYCacs9f32GfxxH0DHiWxlzDRKOA6RlIEg4GBtNoc6W+opiuSIJEamH7kftjU9uAaThLqkpoNdEAvDVEuZzbWhONURlrA6bZxESjgPHKoocm7POsr/S+qMdnePOWuUauGT65sqS6hKM9A3kbr5DJ9dz66aSdJhoj6BqIz7nMKTDRKGh6o0mqJnElUl/h7S280OMaXkwjfytol1SX0Nods9XNeIHt0qLxz21xKEhFcYiOiInGcMzSMPJOz4BnaUwUXzQWetptJJbKazrjkqoS4qn0gne1qCqReHIwnXY8asvDZmkMoy+WpDeaZHFVyWxPZQQmGgVMZyROzSTM1wbf0ljgabe5Zvjkip92++YCz6DqiyVJa+4JGjXlYToilj2VyeFOz83ZVFM6yzMZiYlGAXO8L05dxcRFw49pzPSG9HMJ72o4NaVd+4bjXxW+ucCD4Ye7vB+8ZYty+8GrLSuio39hX8AM53CXtxPkchMNI1+oKsf7Y4OupolQFg5RFg4u6JhGPJUmmdY8WxreF3yhWxr+VXKuP3g15WE6+83SyKTFLA0j30TiKaKJ9GDK4kSpryhe0KLh76WRz+yp+oowATFLw7c0mnK0NOrKwws+DjScw50DhEMB6ssnflE43ZhoFCh+umzdpEVjYa8K9/fSyNfiPvAqCDdWlphodA4QDgZytoJrysMMJFJWSiSDls4Bli8qJRCYWyVEwESjYDnufMCTcU/5/Y71Ltyru0gey6JnsmxRCQeOR/I6ZqHR0jXAskUlOf/g1bpkjk5Lux2kpWtgTrqmwESjYPEtjcm6p+oWuHvKT/HMdx78WSsW8eLhrgW92v5w58CEArhWSmQkhzsjLM/RvTfTmGgUKL6lMZnsKYCGijAdkTip9MJciHake2IZPrnytlW1RBNpdh3pyeu4hcThroEJ/eBZKZGhRBMpjvXFzdIw8otfcLBukoGyhspiVOH4ArU2jnR5cYdl1fn9Ym5yGzDt3L8w9wuLJlK098ZYvqgs5z6+aJh7yqNlgtlnM82URENE9ovISyLyvIjsdG21IrJDRPa6+5qM468RkWYReVVEtmS0n+PGaRaRG92WsLhtY+9y7U+JyKqpzHc+cbwvTnk4mFMl0Ww0ujUFbQt0gd/hrgFqy8OTPn+j0VhVwsraMnbu78zruIWCX+V3Ij94fkzDLA2PweyzmtyFdybJh6Xxe6q6UVU3ub8/BzyoquuAB93fiMgGvK1cTwcuBr4pIv439iZgO96e4uvc8wBXAZ2quhb4GnBDHuY7LzjeF6N2kq4pOLEQbaHu/3DUBWung00n1bDzQMeCrEF1xP3gTcQ9VVVaRECsaKFPS6db2LeAYhqXALe5x7cBl2a036mqMVV9A2gGNovIUqBKVZ9U71t2+7A+/lg/Ai6UubaN1SxxvD8+adcUwOIqr29rz8K0NI50RQcX4+WbTatqOdYXZ/8CzKKaTPmLYEBYVBa2ooWO/cf6CYcCc7LuFExdNBT4LxF5VkS2u7bFqnoUwN03uvblwKGMvi2ubbl7PLx9SB9VTQLdQN0U5zwvON4Xp34KlkZ9RTEiC9fSODLBYO1EeEtTNQB7ji68YHhLZwQRr+LvRKi1BX6DvPJmL+sXVxCcg2s0AKa6sukCVT0iIo3ADhF5ZYxjs50BHaN9rD5DB/YEazvAypUrx57xPOF4f4wzlldNun9RMEBdeTFtvQtPNHqiCXpjyWlzT61pqABgX1vftIw/l3mhpZu1DRUUBSd2PdpYWczhroX3WczGnqM9/N4pjeMfOEtMydJQ1SPuvg34KbAZaHUuJ9x9mzu8BViR0b0JOOLam7K0D+kjIiGgGhiRlqKqN6vqJlXd1NDQMJW3VBCoKh39ceomubDPZ3FV8YJ0Tx2ZYEG9iVIaDrJ8USnN7QtLNJKpNDv3d3DuybUT7ruusYJ9bX0LMg6USXtvjGN9cU5dOvkLwulm0qIhIuUiUuk/Bi4CXgbuBba5w7YB97jH9wJbXUbUaryA99POhdUrIue5eMWVw/r4Y10GPKQL/VMF9ESTJFI66RIiPourShake+qon247jYHGNY0V7FtgorHrSA/98RTnrp64B3nd4kr6YkmOLPASLK+86bk0T1taOcszGZ2puKcWAz91cekQ8B+q+ksReQa4W0SuAg4ClwOo6i4RuRvYDSSBq1XVLzbzceB7QClwv7sB3ALcISLNeBbG1inMd97gr62Y7MI+n8VVxbzY0p2PKRUUg6W7pykQDrCmoZxn3uggndY5WT9oOnjqjeMAk7I01i/2fiRfa+2ds1lDM4EfBzttydy1NCYtGqr6OnBWlvbjwIWj9LkeuD5L+07gjCztUZzoGCc4NlhCZGruqcbKEo73x0ik0hP2QRcyR7oGCAWEhsrpqyC6pqGCgUSKN3ui02rRzAW6BxK09kR56vUOTq4vp7Fy4rGi9Yu9ONDe1t457c+fbl452suSqpLB0ipzkfyV+DRmjP3H+gE4qXZqi38WV5Wg6u0VPl3pp3ORI10DLKkumdbslLWNLhje3jfvReO6n+/mx8+1IAJb37Zi/A5ZWFQWpqGymNdaF5ZLbzi7j/Zw6hx2TYGVESlIXmvtpTgUYMWURWNhrtVobu9j5RTP3Xj4GVTN8zyDSlV55LU2Tl1SyVuaFvGHZzeN32kU1i+uYG9rbx5nV1j0x5Lsa+/jtDkcBAcTjYJkb1sfaxqmnse9EFeF90YT7D7Sw6ZVE/e7T4T6ijBVJaF5Hwx/tbWXY31xrnrHau65+gI2r578eV3XWMnetj7SC6SI5n0vHeXz9+4azBj71Z5WEimd8+45E40CZG9r76APeCo0OkujbQGJxrMHOkkrnDuFH7dcEBHWNFbw2pvzWzQe23sMgAvW1k95rPWLK4nEU4MViOczqspXd7zG957Yz7//9gAA9zx/hGXVJYNFL+cqJhoFRm80wZHuKOsWT93vWVdeTDAgC8o99fQbHYQCwtkrF037a527uo7nDnbSE52/+18/3nyMkxvK8xK32bDMc8s8e2D+F3vcdaSH5rY+FpUVcf19e3hwTyuPvtbOH2xcNuez7Uw0Coy9zke+Pg+iEQwIS6pKONS5cGokPf1GB2c2VVOWx21eR+P9GxpJppVHXm2f9teaDeLJNE+90cE78mBlALxleTUNlcX88uU38zLeXOae5w9TFBTu2n4+1aVFXHXbTpJp5UNnLZvtqY2LiUaB4QcK8+GeAjh1SeWCqZEUTaR4oaVrSn73ibBxRQ115WF+tad1Rl5vpvnN3nYi8RTvXp+fKgyBgHDx6Ut4+NV2Im4P9/lIMpXm3heO8O71jZyypJIHP/MePnXhOq48/yQ2zPEgOJhoFBx7W/soKQrkrdb+aUur2NfeTzSRGv/gAud3B7tIpJTN0xwE9wkGhPee2sivX2mbl9u//nBnC/UVYd6VJ9EA+MAZSxhIpOatdaaq/K+fvkxrT2wwPbmiOMT/8/71XHfJGRRCEW8TjQJib2svj+5tz0vmlM+GZVWk0jrvU0MBnth3jGBAZszSAHjfhsX0RJM8se/4jL3mTNDRH+fBV1q5dOPyvC4M3by6lpqyIv7zpaN5G3Mu8ZX/epW7dh7ir9+7lvdtWDzb05kUJhoFwn/tepOL/vlRDnUM8BcXrM7buH5O+O4FsKf1Y83HeEtTNZUlRTP2mu9e30B9RTE3P7pvxl5zOkmllfteOsoXf7GbREq5fNPkFvONRigY4A/PbuL+l9+cdxcyP/vdYf711/u4YvMK/ub962d7OpPGRKMA6Ikm+PufvcypS6p4/HPv5cPnTH4B1XBOqi2jLBxk9zyPa/REE7zY0p23oG2ulBQF+e/vXM3jzcf53cHCzwr60v17+MT3n+OnvzvMO9bWc8qS/K9e/sTvraEkFOArD7ya97Fni1/tbuV//vhFzl1dWzBuqNEw0ZjDqCpP7jvOZ+5+gWN9Mb70R2dSm+eaNIGAcOqSyoIUDVUllswtFvPU6x2k0srb18ysaAB85LyTqC4t4qaHC9va+I+nDvKd37zBR887iV1f2MIdV22eltepryhm+7vW8Mtdb/LEvmPT8hozyd07D7H9jp2ctqSSm/70nIKv81bYs58n9EYT/OS5Fr78y1e465mDg3sl/9vj+7niO7/l16+08akL13PWiulZW7BhWRV7jvYU1F4Gqsrf/+xlzr5uBz96tmXUubf1RPn/fvYy3/3N65QUBXjrSdO/PmM4FcUhPvzWJh5+rb1gEw72H+vnul/s4p3r6rn2DzZQXhya1qvlv3znak5uKOd//OB3HC3gxX6d/XGu+/luzl1dxw+2n5f3i77ZwERjlnmzO8plNz3J39z9At96ZB+f/fFLvOOGh3ii+Rhff2gvF6yt4/lrL+JT71s3bXM4q2kRvdEkf/ztJwtiYZWq8vWHmvn+UwepKQvz//7wBd76xR388bef5A1XzNHnf9+3hzt+e8CtJ2igOBSclTm/c3098WSanfvn/vnNJJpIsetIN5/98YsUBQL8n8vOIjQDV8rlxSFu/ug5DMRTfPrO56f99aaLbz26j/54ki9ccvqMrA2aCebHuyhQeqIJLv/2E3T0xbll2ybetb6BPUd7+MT3n+MjtzyFKnzu4tOoKJ7ef9MfvbWJ7oEEtzz2Bn/xvWf4xV+/Y8rFEKeLA8f7+dsfvsjT+zu4dOMyvnL5Wdy18xC7j/Rw/8tvcum/Ps63/vQczl9Tx4stXfzs+SN84j1r+LMLVlFZPHMB8OFsXlVLUVB4rPkY71g38y6yyRBNpLjkG4/zqlsbdMOHz5zw3t9TYW1jJX+75RQ+//PdPHugg3NOmrmst8kSiSd54VA355xUw4Hj/dz2xH4uOWtZXhbjzhWkkFwSubBp0ybduXPnbE8jJz77oxf54bOHuPtj5w8poLe3tZcP3/QE7zmlkRuvOHvG5nPgeD//7euPUVMWJpVWSooCXLF5JX9+weo5s8n9lbc+ze8OdPK5D57K1retHDKvg8cj/MVtz7D/WD9XvXM1O3a10j2Q4OG/fc+MZkyNxh9/60kGEil+/tfvmNV5ROJJnj/URUlRkGgixZGuKC8f7iaVVpYtKuUPzlpKU00Z//CL3Xz3sTe47pLTeevKGs5YXj0rcz3/Hx/ivJNr+fZHN834649GOq3s2NPKv//2AHXlYS45ezmvt/dz86P7aO2JsaK2lGO9ccqLg/zk4xewsm5uXoRlIiLPquq4J7kgLA0RuRj4FyAIfFdVvzRdr/XGsX4eaz5GMpVmbWMFZyyrJq1KZUkR4ZBnlu860s1je4/xZxesmpS7Q1W55/kj3LXzEH/17jUjKq6uW1zJbz77XsrDM+tKOamunK/+8UY+++MX2XRSDR39cf7hP/dQXhziis0r8/payVSaXUd6qCgJ0VRTmtN5fK21l0dfa+cz71/PR849acTzK+vK+Mkn3s7V33+Obz/yOmsayvnqn2ycE4IBXlG/f37wNTr747O2yc7r7X1sv+PZEemsZeEg4VCArkiCLz/wCosrS2jtjfLR807iyvNXzcpcvXmF+Oh5J/GvDzfz8uHuWRGu4bT1RPmbu1/gseZjLF9UynMHOvnZ80cAOKupmk+/bz3ff+oAK2rK+NqfbBysJj1fmPOWhogEgdeA9wMtwDPAFaq6O9vxk7U0uiMJvrrjVb7/1EGSWUozBwSWVpfSUFnMCy1dqML7Nyzmmx95K0XBAMf7YvzuYBfH+mJsOX3JkB+FeDLN84e6WNdYQWckzhd+vptHXmvnLU3V3P2x8ykpmh0/+3ioKh/6xuP0RhM8+Jn3ZLU2kqk0T+/vYElVCavry8cMjqbTys4DnTyx7xh3P3NocD/omrIi/uitTWxcsYjV9eWsri+nPItL7nM/fpGf/u4wT15z4ZgBxWQqzb72ftYvrphTqY3PHujgwzc9yUfPO4lPv28dLx/p4eyVi6jKELV0WjncNUBjVTHFoSAHj0e45bHXeWLfcc5asYi3rqzhtKWVbFyxaELvbSCe4ru/eZ2bHtlHSVGQz3/odCqLQxQXBWis9P53wYBwuGuAH+48xNGuKHUVYT753rWz7otv743xgX95lJ5oko+962Teua6Bc06qIRgQWjojJFNKY1XxtM8zmkjxH08d5F8e3EssmeLvf38DW9+2gv54iucPdbG2sYJl1SVz6jM3EXK1NApBNM4HPq+qW9zf1wCo6j9mO36yonG8L8aFX32E3z9zKX/17jWUhoO8fLibfe39hALC8b4YhzoHONw5wNkrF1FXEeZ/3/cKZy6v5m2ravnB0wcZcJkxFcUhPrRxGSfXl7O3tY9f7WnleH+cUEAQgZJQkE+9bx3b3r5qzqff3f/SUT7+/ee48Yqz+dBZyzjSNcAT+47z7IEO+mMpnjvYSUunl91ycn05f7vlFFY3lHO0K8rR7igpVRaVFrF+cSU3/PIVHnqlDRGvNPkVm1eSSis7drfyX7tbSWWIdWNlMTVlYSKJJJFYikg8xUAixRWbV/CPf/SW2TodU0JVufbeXdz+5IHBtpPry/nT807i9if30xdLEkum6Y0mWV1fzkUbFvNvj+9HUd62qpZdR3roHvAq5p5zUg3vXt9AR3+cdJbvcEDEqz5bXcqR7gFuengfR7ujbDl9Mf//H5xecPtwt/VG+V8/fZkdu706XhuWVnH6sip+9FwL/tuvLA7RWFXM4qoSFleV8NaTanjXunqa2/o41BHhaE+UFw91U1se5gNnLmHP0R7e7I5RXBSgtTtKICCcc1IN5cUh0mlFVTm5oYJEKs0Pnj7Ib/YeI5ZMuwyy0wd3Z5wvzCfRuAy4WFX/0v39UeBcVf1ktuOnEtPojSYm5Mr4yXMtfOPXzbze3s8HzljCVe9YTXEoyLce2ccjr7XTF0tSXVrE29fU8cEzl/Ly4W4SKeWv3nPypPZRng1SaeX9X3uE19v7qS0P0+HSgatLi1hUVsTyRaVcsXkl3QMJbn9y/5jbdRYFhc994DQue2sT1WVDz3M0kWL/8X7eaO/n9WP9vN7eT280QXlxiLJwkLJwkIriIq44d0XBnLvR+NXuVl483M3q+jL+4Rd7ON4f56ymajYsqyYYgNX1FXzviTc41DHAB89cwrV/cDqLq0oGrZBHXmvn6w/tpbUnRmVxiFBw5JVtIqX0xU4U/XtLUzX/64Once7JdTP5VvNOR3+ch19t4ysPvEprb4wrzz+J05dV09Ybpa0nRltvlNaeGEe7BgYtWZ9QQDhtaRUtnRE6IwmCAWFxZTGxZJqGymIGEikOHM9e8bmhspjfP3MpW05fwnkn1xasNTEW80k0Lge2DBONzar61xnHbAe2A6xcufKcAwcOZB1rOlBVOiOJEe4SVaWjP05tebjgP2CHuwa49/kj7Gvv49Qllbx9TT2nLqkcUfc/mUqzY3crKVWWVpeytLqEUFBo64nxYks3G1csGtwzwfB4szvKq629vGtd/ZDPSSSepLmtj7c0ZV9XkkoriVR6VNemqvJmT5T23hgVxSFW1ZXP+X0aJkI0kWIgnho1NqSq7DrSw3MHOzl1SRUnN5RTUxYmGBCiCc+ddNrSKqpLh168dPbHSaTSBAKCKrz6Zi/RRIp3rW8YjGnOV+aTaMyIe8owDGMhk6toFIJ0PgOsE5HVIhIGtgL3zvKcDMMwFiRzPuVWVZMi8kngAbyU21tVddcsT8swDGNBMudFA0BV7wPum+15GIZhLHQKwT1lGIZhzBFMNAzDMIycMdEwDMMwcsZEwzAMw8gZEw3DMAwjZ+b84r6JIiK9wPDNhauB7jy+zFwfrx7I1z6Zc/295ns8n3ydw0J4v3P58wdz/xzOh/NXD5SrasO4o6nqvLoBO7O03Zzn15jr4404B3NobnN6vHyfw0J4v3P581cI53A+nL+JvOZCcU/9fIGNl0/m+nudy+cOCuP92jmcW+Plm7zObz66p3ZqDvVT5jN2DqaOncPJY+duaszG+ZvIa85HS+Pm2Z7AHMDOwdSxczh57NxNjdk4fzm/5ryzNAzDMIzpYz5aGoZhGMY0YaJRAIjIChH5tYjsEZFdIvIp114rIjtEZK+7r3Htde74PhH5RsY4lSLyfMbtmIj882y9r5kkX+fQPXeFiLwkIi+KyC9FpH423tNMkedz9yfuvO0SkS/PxvuZaSZx/t4vIs+6z9izIvLejLHOce3NInKjzMYOb/lM7bLb9NyApcBb3eNK4DVgA/Bl4HOu/XPADe5xOfAO4K+Ab4wx7rPAu2b7/RXSOcSrDN0G1Lu/v4y3Sdisv8cCOHd1wEGgwf19G3DhbL+/OXj+zgaWucdnAIczxnoaOB8Q4H7gAzP9fszSKABU9aiqPuce9wJ7gOXAJXhfPNz9pe6YflV9DIhmGQ4AEVkHNAK/mcapzxnyeA7F3crdVV4VcGT638HskcdzdzLwmqq2u79/BXx4mqc/60zi/P1OVf3P1C6gRESKRWQpUKWqT6qnILf7fWYSE40CQ0RW4V2JPAUsVtWj4H0w8UQgV64A7nIfvgXFVM6hqiaAjwMv4YnFBuCWaSQXUu8AAAOBSURBVJzunGKKn79m4FQRWSUiIbwfvBXTN9u5xyTO34eB36lqDE9oWjKea3FtM4qJRgEhIhXAj4FPq2rPFIfbCvxg6rMqLKZ6DkWkCE80zgaWAS8C1+R1knOUqZ47Ve3EO3d34Vm4+4FkPuc4l5no+ROR04EbgI/5TVkOm/GLPhONAsH9WP0Y+L6q/sQ1tzqTFXffluNYZwEhVX12WiY7R8nTOdwIoKr7nJV2N/D2aZrynCFfnz9V/bmqnquq5+PViNs7XXOeS0z0/IlIE/BT4EpV3eeaW4CmjGGbmAXXqIlGAeB857cAe1T1qxlP3Qtsc4+3AffkOOQVLDArI4/n8DCwQUT8wm7vx/NRz1vy+fkTkUZ3XwN8Avhufmc795jo+RORRcB/Ateo6uP+wc6F1Ssi57kxryT373z+mO3MAruNf8PLRFE8V8jz7vZBvGyUB/Gu1h4EajP67Ac6gD68K5QNGc+9Dpw62++rUM8hXlbQHjfWz4G62X5/BXTufgDsdrets/3e5uL5A/4e6M849nmg0T23CXgZ2Ad8A7dAeyZvtiLcMAzDyBlzTxmGYRg5Y6JhGIZh5IyJhmEYhpEzJhqGYRhGzphoGIZhGDljomEYM4yI/JWIXDmB41eJyMvTOSfDyJXQbE/AMBYSIhJS1W/N9jwMY7KYaBjGBHFF536JV3TubLxS11cCpwFfBSqAY8CfqepREXkYeAK4ALhXRCqBPlX9iohsBL4FlOEt2PoLVe0UkXOAW4EI8NjMvTvDGBtzTxnG5DgFuFlV3wL0AFcDXwcuU1X/B//6jOMXqeq7VfWfho1zO/BZN85LwLWu/d+A/6FejSbDmDOYpWEYk+OQnqgL9O/A3+FtmLPDbaYWBI5mHH/X8AFEpBpPTB5xTbcBP8zSfgfwgfy/BcOYOCYahjE5htff6QV2jWEZ9E9gbMkyvmHMCcw9ZRiTY6WI+AJxBfBboMFvE5Eitx/CqKhqN9ApIu90TR8FHlHVLqBbRN7h2j+S/+kbxuQwS8MwJsceYJuIfBuvSunXgQeAG517KQT8M952nWOxDfiWiJThVR/+c9f+58CtIhJx4xrGnMCq3BrGBHHZU79Q1TNmeSqGMeOYe8owDMPIGbM0DMMwjJwxS8MwDMPIGRMNwzAMI2dMNAzDMIycMdEwDMMwcsZEwzAMw8gZEw3DMAwjZ/4vMPhA39ojqPoAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
"source": [
"sorted_data['inc'][-200:].plot()"
]
@@ -252,10 +2300,8 @@
},
{
"cell_type": "code",
- "execution_count": null,
- "metadata": {
- "collapsed": true
- },
+ "execution_count": 14,
+ "metadata": {},
"outputs": [],
"source": [
"first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n",
@@ -274,7 +2320,7 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 15,
"metadata": {},
"outputs": [],
"source": [
@@ -298,9 +2344,32 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 16,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 16,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAD8CAYAAAC7IukgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+w1fV95/HnC0GwiSgQMPxQYSPJBE0Xyx10x+5uxC2QHxOwNQ2rVWbqDIniju1mVqTJjqnSmZhpYsu6ISE1FbWK1MSRTaF4jTq1XQJcgr+Q0Hu7UkSIXOaiYmegXnjvH9/Pqd97cjmcc/lyftz7esycOd/7OZ/P53zul8N9n8+P7/ejiMDMzKxIwxrdADMzG3wcXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwjm4mJlZ4YY3ugH19JGPfCSmTp3a6GaYmbWU7du3H4qI8bWUGVLBZerUqXR0dDS6GWZmLUXSP9daxsNiZmZWOAcXMzMrnIOLmZkVrqrgImmPpFckvSipI6V9Q9KbKe1FSZ/N5V8uqUvSbknzcumzUj1dklZKUkofKenxlL5F0tRcmcWSOtNjcS59Wsrbmcqeffqnw8zMilBLz+XqiJgZEW25tPtS2syI2AAgaQawCLgUmA98V9JZKf8qYAkwPT3mp/SbgcMRcQlwH3BvqmsscBdwBTAbuEvSmFTm3vT+04HDqQ4zM2sCZ2JYbAGwNiKORcTrQBcwW9JEYHREbI5sh7KHgIW5MmvS8RPANalXMw9oj4ieiDgMtAPz02tzUl5S2VJdLe/gu0f53e9v5uCRo41uipnZgFQbXAJ4WtJ2SUty6bdJelnSD3M9isnAG7k8+1La5HRcnt6nTET0Au8A4yrUNQ54O+Utr6vlrfxpJ9v29LDymc5GN8XMbECqvc7lqojYL2kC0C7pF2RDXPeQBZ57gG8Dvw+on/JRIZ0BlKlUVx8pGC4BuOiii/rL0jQ+8fWNHOs98W8/P7JlL49s2cvI4cPYveIzDWyZmVltquq5RMT+9HwQeBKYHRFvRcTxiDgB/IBsTgSyXsSFueJTgP0pfUo/6X3KSBoOnAf0VKjrEHB+ylteV3nbV0dEW0S0jR9f0wWmdffCHVfzhZmTGDUi+2cZNWIYC2ZO4oVlVze4ZWZmtTllcJH0IUnnlo6BucCraQ6l5Frg1XS8HliUVoBNI5u43xoRB4Ajkq5McyY3AU/lypRWgl0HPJvmZTYBcyWNScNuc4FN6bXnUl5S2VJdLWvC6FGcO3I4x3pPMHL4MI71nuDckcOZcO6oRjfNzKwm1QyLXQA8mVYNDwcejYi/lfSwpJlkw1F7gC8DRMROSeuA14BeYGlEHE913QI8CJwDbEwPgAeAhyV1kfVYFqW6eiTdA2xL+e6OiJ50vAxYK2kFsCPV0fIOvXeMG664mOtnX8SjW/fS7Ul9M2tByjoBQ0NbW1v43mJmZrWRtL3sMpRT8hX6ZmZWOAcXMzMrnIOLmZkVzsHFzMwK5+BiZmaFc3BpUb7/mJk1MweXFuX7j5lZM6v23mLWJHz/MTNrBe65tBjff8zMWoGDS4vx/cfMrBV4WKwF+f5jZtbsfG+xKhx89yi3PbaD+6+/3D0EMxtyfG+xM8Qrs8zMauNhsQq8MsvMbGDcc6nAK7PMzAbGwaUCr8wyMxuYqoKLpD2SXpH0oqSOlDZWUrukzvQ8Jpd/uaQuSbslzculz0r1dElambY7Jm2J/HhK3yJpaq7M4vQenZIW59KnpbydqezZp386flVpZdaTt17FDVdcTPd7x87E25iZDSpVrRaTtAdoi4hDubRvAT0R8U1JdwJjImKZpBnAY8BsYBLwDPDxiDguaStwO/AzYAOwMiI2SroV+PWI+IqkRcC1EfElSWOBDqCNbDvl7cCsiDictlL+cUSslfQ94KWIWFXp9/BOlGZmtav3arEFwJp0vAZYmEtfGxHHIuJ1oAuYLWkiMDoiNkcW0R4qK1Oq6wngmtSrmQe0R0RPRBwG2oH56bU5KW/5+5uZWYNVG1wCeFrSdklLUtoFEXEAID1PSOmTgTdyZfeltMnpuDy9T5mI6AXeAcZVqGsc8HbKW16XmZk1WLVLka+KiP2SJgDtkn5RIa/6SYsK6QMpU6muvo3JguESgIsuuqi/LGZmVrCqei4RsT89HwSeJJtPeSsNdZGeD6bs+4ALc8WnAPtT+pR+0vuUkTQcOA/oqVDXIeD8lLe8rvK2r46ItohoGz9+fDW/rpmZnaZTBhdJH5J0bukYmAu8CqwHSqu3FgNPpeP1wKK0AmwaMB3YmobOjki6Ms2Z3FRWplTXdcCzaV5mEzBX0pi0Gm0usCm99lzKW/7+ZmbWYNUMi10APJlWDQ8HHo2Iv5W0DVgn6WZgL/BFgIjYmVZyvQb0Aksj4niq6xbgQeAcYGN6ADwAPCypi6zHsijV1SPpHmBbynd3RPSk42XAWkkrgB2pDjMzawK+caWZmVXkG1eamVlTcHAxM7PCObiYmVnhHFzMzKxwDi5mZlY4BxczMyucg4uZmRXOwcXMzArn4GJmZoVzcDEzs8I5uJiZWeEcXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhas6uEg6S9IOST9JP39D0puSXkyPz+byLpfUJWm3pHm59FmSXkmvrVTaO1nSSEmPp/QtkqbmyiyW1Jkei3Pp01LezlT27NM7FWZmVpRaei63A7vK0u6LiJnpsQFA0gxgEXApMB/4rqSzUv5VwBJgenrMT+k3A4cj4hLgPuDeVNdY4C7gCmA2cJekManMven9pwOHUx1mZtYEqgoukqYAnwP+oorsC4C1EXEsIl4HuoDZkiYCoyNic0QE8BCwMFdmTTp+Argm9WrmAe0R0RMRh4F2YH56bU7KSypbqsvMzBqs2p7LnwF3ACfK0m+T9LKkH+Z6FJOBN3J59qW0yem4PL1PmYjoBd4BxlWoaxzwdspbXpeZmTXYKYOLpM8DByNie9lLq4CPATOBA8C3S0X6qSYqpA+kTKW6+pC0RFKHpI7u7u7+spiZWcGq6blcBXxB0h5gLTBH0iMR8VZEHI+IE8APyOZEIOtFXJgrPwXYn9Kn9JPep4yk4cB5QE+Fug4B56e85XX1ERGrI6ItItrGjx9fxa9rZman65TBJSKWR8SUiJhKNlH/bET8XppDKbkWeDUdrwcWpRVg08gm7rdGxAHgiKQr05zJTcBTuTKllWDXpfcIYBMwV9KYNOw2F9iUXnsu5SWVLdVlZmYNNvzUWU7qW5Jmkg1H7QG+DBAROyWtA14DeoGlEXE8lbkFeBA4B9iYHgAPAA9L6iLrsSxKdfVIugfYlvLdHRE96XgZsFbSCmBHqsPMzJqAsk7A0NDW1hYdHR2NboaZWUuRtD0i2mop4yv0zcyscA4uZmZWOAcXMzMrnIOLmZkVzsHFzMwK5+BiZmaFc3AxM7PCObiYmVnhHFzMbMg7+O5Rfvf7mzl45GijmzJoOLiY2ZC38qedbNvTw8pnOhvdlEHjdO4tZmbW0j7x9Y0c6/1gm6pHtuzlkS17GTl8GLtXfKaBLWt97rmY2ZD1wh1X84WZkxg1IvtTOGrEMBbMnMQLy65ucMtan4OLmTWdes2BTBg9inNHDudY7wlGDh/Gsd4TnDtyOBPOHXVG33cocHAxs6ZTzzmQQ+8d44YrLubJW6/ihisupvu9Y2f8PYcC33LfzJpG+RxIiedAGsu33DezluY5kMHDwcXMmobnQAaPqoOLpLMk7ZD0k/TzWEntkjrT85hc3uWSuiTtljQvlz5L0ivptZWSlNJHSno8pW+RNDVXZnF6j05Ji3Pp01LezlT27NM7FWbWDDwHMjhUPeci6b8DbcDoiPi8pG8BPRHxTUl3AmMiYpmkGcBjwGxgEvAM8PGIOC5pK3A78DNgA7AyIjZKuhX49Yj4iqRFwLUR8SVJY4GO9L4BbAdmRcRhSeuAH0fEWknfA16KiFWVfgfPuZiZ1e6MzblImgJ8DviLXPICYE06XgMszKWvjYhjEfE60AXMljSRLDBtjiyiPVRWplTXE8A1qVczD2iPiJ6IOAy0A/PTa3NS3vL3NzOzBqt2WOzPgDuA/DKOCyLiAEB6npDSJwNv5PLtS2mT03F5ep8yEdELvAOMq1DXOODtlLe8rj4kLZHUIamju7u7yl/XzMxOxymDi6TPAwcjYnuVdaqftKiQPpAylerqmxixOiLaIqJt/Pjx/WUxM7OCVdNzuQr4gqQ9wFpgjqRHgLfSUBfp+WDKvw+4MFd+CrA/pU/pJ71PGUnDgfOAngp1HQLOT3nL6zKrie+Ia1a8UwaXiFgeEVMiYiqwCHg2In4PWA+UVm8tBp5Kx+uBRWkF2DRgOrA1DZ0dkXRlmjO5qaxMqa7r0nsEsAmYK2lMWo02F9iUXnsu5S1/f7Oa+I64ZsU7nbsifxNYJ+lmYC/wRYCI2JlWcr0G9AJLI+J4KnML8CBwDrAxPQAeAB6W1EXWY1mU6uqRdA+wLeW7OyJ60vEyYK2kFcCOVIdZ1XxHXLMzx7d/sSHr4LtHWbFhF0/v/CVH3z/BqBHDmHfpR/na5z7pi/bMcnz7F7Ma1HI1eKvOy7Rqu631ObjYkFbt1eCtOi/Tqu221udhMbMKWvUuva3abmtOHhYzK1ir3qW3Vdttg4eDi1kFrXqX3lZttw0ep7MU2WxIKM3LXD/7Ih7dupfuFpkcb9V22+DgORczM6vIcy5mZtYUHFzMzKxwDi5mZgXyhasZBxczswL5wtWMV4tZSzr47lFue2wH919/uZfXWlPwjVD7cs/FWpK/HVqz8YWrfbnnYi3F3w4/4N5bc/GFq32552Itxd8OP+DeW/Op9kaoQ4F7LtZ0Kn0j97dD996a2fdv/OA6wxULL2tgSxrvlD0XSaMkbZX0kqSdkv44pX9D0puSXkyPz+bKLJfUJWm3pHm59FmSXkmvrUzbHZO2RH48pW+RNDVXZrGkzvRYnEuflvJ2prJnF3NKbKCKWoJ5qm/kQ/3boXtv1gqq6bkcA+ZExHuSRgB/L6m0PfF9EfGn+cySZpBtU3wpMAl4RtLH01bHq4AlwM+ADcB8sq2ObwYOR8QlkhYB9wJfkjQWuAtoAwLYLml9RBxOee6LiLWSvpfqWDXwU2GnKx8UVlz7qZrLV/uNfKh/O3TvzVrBKXsukXkv/TgiPSrdkGwBsDYijkXE60AXMFvSRGB0RGyO7IZmDwELc2XWpOMngGtSr2Ye0B4RPSmgtAPz02tzUl5S2VJdVmef+PpGpt75NzyyZS8RWVCYeuff8Imvbzx14Rx/I6/eUO+9WfOras5F0lnAduAS4H9HxBZJnwFuk3QT0AF8NQWAyWQ9k5J9Ke39dFyeTnp+AyAieiW9A4zLp5eVGQe8HRG9/dRldfbCHVefdC/6WvgbefWGeu/Nml9Vq8Ui4nhEzASmkPVCLiMbgvoYMBM4AHw7ZVd/VVRIH0iZSnX1IWmJpA5JHd3d3f1lsdNUZFDwN/LBz7dHGRpqWi0WEW9Leh6Yn59rkfQD4Cfpx33AhbliU4D9KX1KP+n5MvskDQfOA3pS+qfLyjwPHALOlzQ89V7ydZW3eTWwGrJb7tfy+1r1ito7xN/IB7/TnZuz1nDK/VwkjQfeT4HlHOBpssn07RFxIOX5Q+CKiFgk6VLgUWA22YT+T4HpEXFc0jbgvwFbyCb0/1dEbJC0FPhURHwlTej/dkT8bprQ3w78RmrOz4FZEdEj6a+BH+Um9F+OiO9W+l28n4tZ45Qv2CjxEurmN5D9XKrpuUwE1qR5l2HAuoj4iaSHJc0kG47aA3wZICJ2SloHvAb0AkvTSjGAW4AHgXPIVomVZnwfAB6W1EXWY1mU6uqRdA+wLeW7OyJ60vEyYK2kFcCOVIeZNami5uasNZwyuETEy8Dl/aTfWKHMnwB/0k96B/ArYx0RcRT44knq+iHww37S/x9Z78jMWoAXbAwtvkLfzOqmqLk5a36nnHMZTDznYmZWu4HMufjGlWZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwjm4mJlZ4RxczMyscA4uZmZWOAcXszrzfiY2FDi4mNVZfj8Ts8HKN640q5Py/Uwe2bKXR7bs9X4mNii552JWJy/ccTVfmDmJUSOy/3ajRgxjwcxJvLDs6ga3zKx4Di5mdeL9TGwocXCxQa3ZJs9L+5k8eetV3HDFxXS/d6zRTTI7I04ZXCSNkrRV0kuSdkr645Q+VlK7pM70PCZXZrmkLkm7Jc3Lpc+S9Ep6baUkpfSRkh5P6VskTc2VWZzeo1PS4lz6tJS3M5U9u5hTYoNJs02ef//GNlYsvIwZk0azYuFlfP/GmrbIMGsZp9wsLAWAD0XEe5JGAH8P3A78NtATEd+UdCcwJiKWSZoBPEa2BfEk4Bng4xFxXNLWVPZnwAZgZURslHQr8OsR8RVJi4BrI+JLksYCHUAbEMB2YFZEHJa0DvhxRKyV9D3gpYhYVel38WZhQ0f55HmJJ8/NandGNguLzHvpxxHpEcACYE1KXwMsTMcLgLURcSwiXge6gNmSJgKjI2JzZBHtobIypbqeAK5JQW0e0B4RPRFxGGgH5qfX5qS85e9v5snzBmm2YUhrnKrmXCSdJelF4CDZH/stwAURcQAgPU9I2ScDb+SK70tpk9NxeXqfMhHRC7wDjKtQ1zjg7ZS3vK7yti+R1CGpo7u7u5pf1wYBT543RrMNQ1rjVHWdS0QcB2ZKOh94UtJlFbKrvyoqpA+kTKW6+iZGrAZWQzYs1l8eG5xKk+fXz76IR7fupdvfps8YX8Nj5Wq6iDIi3pb0PDAfeEvSxIg4kIa8DqZs+4ALc8WmAPtT+pR+0vNl9kkaDpwH9KT0T5eVeR44BJwvaXjqveTrMgPoM1m+YmGl70N2ul6442pWbNjF0zt/ydH3TzBqxDDmXfpRvva5Tza6adYg1awWG596LEg6B/gvwC+A9UBp9dZi4Kl0vB5YlFaATQOmA1vT0NkRSVemOZObysqU6roOeDbNy2wC5koak1ajzQU2pdeeS3nL39/M6szDkFaump7LRGCNpLPIgtG6iPiJpM3AOkk3A3uBLwJExM60kus1oBdYmobVAG4BHgTOATamB8ADwMOSush6LItSXT2S7gG2pXx3R0RPOl4GrJW0AtiR6jCzBvEwpOWdcinyYHImlyIffPcotz22g/uvv9zf1szqpFX/37Vau8/IUmSrjlfJmNVfq/6/a9V218I9l9Pki/XM6q9V/9+1arvdc2kAX6xnVn+t+v+uVds9EA4up6mWVTK+etlaWTN9flt1dVqrtnsgHFwKUO2dbofCOKsNXs32+W3VO0y3artr5TmXOmjVcVYz8OfXPOfStIbSOKsNPv78Nq9mGqos5+BSB4NhnLWZP8R2Zg2Gz+9g1WxDlXk13VvMBq7Vr17Of4hXXPupRjfH6qzVP7+DTSvcKNRzLlaRx9vNms/Bd4+e9EahZ6JH6TkXK5zH282aTysMVXpYzCpqhQ+x2VDU7EOVDi52Ss3+ITYbipp9vyLPuQxxrXZ3VjOrP8+5WM2aeSmjmbUuD4sNUa2wlLHVuBdo9oFqtjm+UNJzknZJ2inp9pT+DUlvSnoxPT6bK7NcUpek3ZLm5dJnSXolvbYybXdM2hL58ZS+RdLUXJnFkjrTY3EufVrK25nKnl3MKRkavAqseO4Fmn2gmp5LL/DViPi5pHOB7ZLa02v3RcSf5jNLmkG2TfGlwCTgGUkfT1sdrwKWAD8DNgDzybY6vhk4HBGXSFoE3At8SdJY4C6gDYj03usj4nDKc19ErJX0vVTHqoGfiqHFq8CK416g2a86Zc8lIg5ExM/T8RFgFzC5QpEFwNqIOBYRrwNdwGxJE4HREbE5slUEDwELc2XWpOMngGtSr2Ye0B4RPSmgtAPz02tzUl5S2VJdVqWhcnfWM829QLNfVdOcSxquuhzYAlwF3CbpJqCDrHdzmCzw/CxXbF9Kez8dl6eTnt8AiIheSe8A4/LpZWXGAW9HRG8/dVmVmn0pY6twL9DsV1W9WkzSh4EfAX8QEe+SDUF9DJgJHAC+XcraT/GokD6QMpXqKm/3Ekkdkjq6u7v7y2J22twLNOurqp6LpBFkgeWvIuLHABHxVu71HwA/ST/uAy7MFZ8C7E/pU/pJz5fZJ2k4cB7Qk9I/XVbmeeAQcL6k4an3kq+rj4hYDayG7DqXan5fs1q5F2jWVzWrxQQ8AOyKiO/k0ifmsl0LvJqO1wOL0gqwacB0YGtEHACOSLoy1XkT8FSuTGkl2HXAs2leZhMwV9IYSWOAucCm9NpzKS+pbKkuMzNrsGp6LlcBNwKvSHoxpf0R8F8lzSQbjtoDfBkgInZKWge8RrbSbGlaKQZwC/AgcA7ZKrGNKf0B4GFJXWQ9lkWprh5J9wDbUr67I6InHS8D1kpaAexIdZiZWRPw7V/MzKwi3/7FzMyagoOLmVkTavWtxR1czKzl/5ANRq1+OyHfuNLM+vwhW3HtpxrdnCFtsNxOyBP6ZkNY+R+yklb7QzaYHHz3KCs27OLpnb/k6PsnGDViGPMu/Shf+9wnB3TXhyLu1u0JfTOrie+L1nyKvp1Qo4bXPCxm1oTqtTeM74vWnIrYWrzRw2sOLmZNqJ5zIEX8IbNiFXE7oRfuuPqkw2v14OBi1kQa8W3T90UbnBrdK/Wci1kT8RxI8YbyMutG3q3bPRezJtLob5uD0VBeZt3IXqmDi1mT8RxIMRo9oT3U+ToXMxuUir5eZCjzdS5mZomHGBvLw2JmNmh5iLFxPCxmZmYVeVjMzMyawimDi6QLJT0naZeknZJuT+ljJbVL6kzPY3JllkvqkrRb0rxc+ixJr6TXVkpSSh8p6fGUvkXS1FyZxek9OiUtzqVPS3k7U9mzizklZmZ2uqrpufQCX42ITwJXAkslzQDuBH4aEdOBn6afSa8tAi4F5gPflXRWqmsVsASYnh7zU/rNwOGIuAS4D7g31TUWuAu4ApgN3JULYvcC96X3P5zqMDOzJnDK4BIRByLi5+n4CLALmAwsANakbGuAhel4AbA2Io5FxOtAFzBb0kRgdERsjmyi56GyMqW6ngCuSb2aeUB7RPRExGGgHZifXpuT8pa/v5mZNVhNcy5puOpyYAtwQUQcgCwAARNStsnAG7li+1La5HRcnt6nTET0Au8A4yrUNQ54O+Utr8vMzBqs6uAi6cPAj4A/iIh3K2XtJy0qpA+kTKW6+jZGWiKpQ1JHd3d3f1nMzKxgVQUXSSPIAstfRcSPU/JbaaiL9Hwwpe8DLswVnwLsT+lT+knvU0bScOA8oKdCXYeA81Pe8rr6iIjVEdEWEW3jx4+v5tc1M7PTVM1qMQEPALsi4ju5l9YDpdVbi4GncumL0gqwaWQT91vT0NkRSVemOm8qK1Oq6zrg2TQvswmYK2lMmsifC2xKrz2X8pa/v5mZNVg1V+hfBdwIvCLpxZT2R8A3gXWSbgb2Al8EiIidktYBr5GtNFsaEcdTuVuAB4FzgI3pAVnwelhSF1mPZVGqq0fSPcC2lO/uiOhJx8uAtZJWADtSHWZm1gR8hb6ZmVXkK/St6Q3ljZvMhhIHF6ur/MZNZjZ4+a7IVhfeuMlsaHHPxerCe8ObDS0OLlYX3rjJbGjxsJjVjTduMhs6vBTZzMwq8lJkMzNrCg4uZmZWOAeXQcwXLJpZozi4DGK+YNHMGsWrxQYhX7BoZo3mnssg5AsWzYrnYebaOLgMQr5g0ax4HmaujYfFBilfsGhWDA8zD4wvojQzq+Dgu0dZsWEXT+/8JUffP8GoEcOYd+lH+drnPjlkRgPOyEWUkn4o6aCkV3Np35D0pqQX0+OzudeWS+qStFvSvFz6LEmvpNdWpq2OSdshP57St0iamiuzWFJneizOpU9LeTtT2bNr+aXNzKrlYeaBqWbO5UFgfj/p90XEzPTYACBpBtkWxZemMt+VdFbKvwpYAkxPj1KdNwOHI+IS4D7g3lTXWOAu4ApgNnCXpDGpzL3p/acDh1MdZmZnRGmY+clbr+KGKy6m+71jjW5S0zvlnEtE/F2+N3EKC4C1EXEMeF1SFzBb0h5gdERsBpD0ELAQ2JjKfCOVfwK4P/Vq5gHtEdGTyrQD8yWtBeYA16cya1L5VVW20cysJt+/8YMRoRULL2tgS1rH6awWu03Sy2nYrNSjmAy8kcuzL6VNTsfl6X3KREQv8A4wrkJd44C3U97yuszMrAkMNLisAj4GzAQOAN9O6eonb1RIH0iZSnX9CklLJHVI6uju7j5ZNjMzK9CAgktEvBURxyPiBPADsjkRyHoRF+ayTgH2p/Qp/aT3KSNpOHAe0FOhrkPA+SlveV39tXV1RLRFRNv48eNr/VXNzGwABhRcJE3M/XgtUFpJth5YlFaATSObuN8aEQeAI5KuTPMpNwFP5cqUVoJdBzwb2froTcBcSWPSsNtcYFN67bmUl1S2VJeZmTWBU07oS3oM+DTwEUn7yFZwfVrSTLLhqD3AlwEiYqekdcBrQC+wNCKOp6puIVt5dg7ZRP7GlP4A8HCa/O8hW21GRPRIugfYlvLdXZrcB5YBayWtAHakOszMrEn4IkozM6toIBdRDqngIqkb+Od+XvoI2VxOq3G768vtrq9WbTe0bttP1u6LI6KmSeshFVxORlJHrVG5Gbjd9eV211erthtat+1Fttt3RTYzs8I5uJiZWeEcXDKrG92AAXK768vtrq9WbTe0btsLa7fnXMzMrHDuuZiZWeEGZXA5yR40/17S5rSnzP+RNDqlj5C0JqXvkrQ8V+b5tC9Nad+aCU3U7rMl/WVKf0nSp3Nl+t07pwXaXe/zfaGk59K/+05Jt6f0sZLa035B7bkbs9a8X1ELtLtu57zWdksal/K/J+n+srrq/Rkvsu3NfM5/S9L2dG63S5qTq6u2cx4Rg+4B/CfgN4BXc2nbgP+cjn8fuCcdX0+2TQDAr5HdcWBq+vl5oK1J270U+Mt0PAHYDgxLP28F/gPZTT43Ap9pkXbX+3xPBH4jHZ8L/CMwA/jPDIkDAAADj0lEQVQWcGdKvxO4Nx3PAF4CRgLTgH8Czqr3OS+43XU75wNo94eA3wS+AtxfVle9P+NFtr2Zz/nlwKR0fBnw5kDP+aDsuUTE35HdSibvE8DfpeN24HdK2YEPKbsR5jnAvwLv1qOd5Wps9wzgp6ncQeBtoE3Zfd9GR8TmyD4Rpb1zmrrdZ7J9JxMRByLi5+n4CLCLbPuGBWT7BJGeS+fv3/YriojXgdJ+RXU950W1+0y172RqbXdE/EtE/D1wNF9Pgz7jhbS93gbQ7h0RUboR8E5glLJ7RdZ8zgdlcDmJV4EvpOMv8sEdl58A/oVs64C9wJ/GB/cwA/jL1HX9n2e6630SJ2v3S8ACScOV3SR0Vnqt0t459VRru0sacr6VbYh3ObAFuCCym62SnkvDFgPZr+iMOs12l9T9nFfZ7pNp6Gf8NNte0grn/HeAHZFt/ljzOR9KweX3gaWStpN1D/81pc8GjgOTyIYMvirp36XXboiITwH/MT1urG+TgZO3+4dk/8AdwJ8B/5fsZqE17XdzBtXabmjQ+Zb0YeBHwB9ERKVeayF7DBWlgHZDA855De0+aRX9pNXlM15A26EFzrmkS8m2k/9yKamfbBXP+ZAJLhHxi4iYGxGzgMfIxp0hm3P524h4Pw3T/ANpmCYi3kzPR4BHacxQQr/tjojeiPjDiJgZEQuA84FOKu+d08ztbsj5ljSC7D/dX0XEj1PyW2kYoDQEczClD2S/omZud93PeY3tPpmGfMYLanvTn3NJU4AngZsiovR3suZzPmSCS2lFhqRhwNeB76WX9gJzlPkQcCXwizRs85FUZgTweT7Yt6bh7Zb0a6m9SPotoDciXovKe+c0bbsbcb7T+XkA2BUR38m9lN9jKL9f0ED2K2radtf7nA+g3f1qxGe8qLY3+zmXdD7wN8DyiPiHUuYBnfNKs/2t+iD7pnwAeJ8s4t4M3E62UuIfgW/ywQWkHwb+mmzy6jXgf8QHqz22Ay+n1/6ctMKmSdo9FdhNNkH3DNldS0v1tJF9YP8JuL9Uppnb3aDz/ZtkXfuXgRfT47PAOLJFB53peWyuzNfSed1NbrVMPc95Ue2u9zkfYLv3kC0WeS99tmY06DNeSNub/ZyTfRH8l1zeF4EJAznnvkLfzMwKN2SGxczMrH4cXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwv1/q4MdzlEclD8AAAAASUVORK5CYII=\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
"source": [
"yearly_incidence.plot(style='*')"
]
@@ -314,9 +2383,54 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 17,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "2014 1600941\n",
+ "1991 1659249\n",
+ "1995 1840410\n",
+ "2012 2175217\n",
+ "2003 2234584\n",
+ "2019 2254386\n",
+ "2006 2307352\n",
+ "2017 2321583\n",
+ "2001 2529279\n",
+ "1992 2574578\n",
+ "1993 2703886\n",
+ "2018 2705325\n",
+ "1988 2765617\n",
+ "2007 2780164\n",
+ "1987 2855570\n",
+ "2016 2856393\n",
+ "2011 2857040\n",
+ "2008 2973918\n",
+ "1998 3034904\n",
+ "2002 3125418\n",
+ "2009 3444020\n",
+ "1994 3514763\n",
+ "1996 3539413\n",
+ "2004 3567744\n",
+ "1997 3620066\n",
+ "2015 3654892\n",
+ "2000 3826372\n",
+ "2005 3835025\n",
+ "1999 3908112\n",
+ "2010 4111392\n",
+ "2013 4182691\n",
+ "1986 5115251\n",
+ "1990 5235827\n",
+ "1989 5466192\n",
+ "dtype: int64"
+ ]
+ },
+ "execution_count": 17,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"yearly_incidence.sort_values()"
]
@@ -331,9 +2445,32 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 18,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 18,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEKCAYAAAAGvn7fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGalJREFUeJzt3X2UJXV95/H3h5kBhmkYjAONDsr4QAhIqzgXXWQ13WhcdNCcGE5QQcVIGo0PRCdnM8v6sLrLOj5MsphgkklUiAodw8PZyBjUE2hQNEgPqA2OEBdmlSEMAjLSMAuMfPePX7XctP1wq+69XcXPz+ucPn3vrbpVn/rdut9b9auqexURmJlZPvaqO4CZmfWWC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLzNJ+TnzVqlWxZs2aWYc9+OCDrFixop+zr6zJ2aDZ+Zytuibnc7bqyubbunXrPRFxUFczjYi+/a1duzbmctVVV805rG5NzhbR7HzOVl2T8zlbdWXzARPRZe11V4yZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWWmVGGX9B5JN0u6SdJFkvbtVzAzM6um48IuaTXwbqAVEUcDS4DX9SuYmZlVU7YrZimwXNJSYD/gzt5HMjOzbihK/OappLOAc4DdwFcj4tRZxhkFRgEGBwfXjo2NzTqtqakpBgYGqmTuuyZng97nm9yxq2fTGlwOO3d3Pv7Q6pU9m/dCftVe115yturK5hsZGdkaEa1u5tlxYZf0JOAS4BTgfuAfgIsj4vNzPafVasXExMSsw8bHxxkeHi6bd1E0ORv0Pt+aDVt6Nq31Q3vYNNn5N1Vs37iuZ/NeyK/a69pLzlZd2XySui7sZbpiXg7cHhE/iYhHgUuBF3czczMz670yhf1HwH+QtJ8kAS8DtvUnlpmZVdVxYY+I64CLgRuAyeK5m/uUy8zMKir1tb0R8UHgg33KYmZmPeArT83MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy03Fhl3SEpO+0/f1M0h/1M5yZmZXX8U/jRcQtwPMBJC0BdgCX9SmXmZlVVLUr5mXA/4mI/9vLMGZm1j1FRPknSZ8BboiIv5hl2CgwCjA4OLh2bGxs1mlMTU0xMDBQet6LocnZoPf5Jnfs6tm0BpfDzt2djz+0emXP5r2Q9nbr5TKXMd/yNnm9c7bqyuYbGRnZGhGtbuZZurBL2hu4E3hOROycb9xWqxUTExOzDhsfH2d4eLjUvBdLk7NB7/Ot2bClZ9NaP7SHTZMd9/CxfeO6ns17Ie3t1stlLmO+5W3yeuds1ZXNJ6nrwl6lK+aVpK31eYu6mZnVo0phfz1wUa+DmJlZb5Qq7JL2A34LuLQ/cczMrFudd4YCEfEQ8OQ+ZTEzsx7wladmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmSn703gHSrpY0g8kbZN0XL+CmZlZNaV+Gg84F7giIk6WtDewXx8ymZlZFzou7JIOAF4KnA4QEY8Aj/QnlpmZVaWI6GxE6fnAZuD7wPOArcBZEfHgjPFGgVGAwcHBtWNjY7NOb2pqioGBgerJ+6jJ2aD3+SZ37OrZtAaXw87dnY8/tHplz+a9kPZ26+UylzHf8jZ5vXO26srmGxkZ2RoRrW7mWaawt4B/AY6PiOsknQv8LCLeP9dzWq1WTExMzDpsfHyc4eHh8okXQZOzQe/zrdmwpWfTWj+0h02Tnffwbd+4rmfzXkh7u/VymcuYb3mbvN45W3Vl80nqurCXOXh6B3BHRFxX3L8YeEE3Mzczs97ruLBHxF3AjyUdUTz0MlK3jJmZNUjZs2LeBXyhOCPmNuAtvY9kZmbdKFXYI+I7QFd9P2Zm1l++8tTMLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsM6V+QUnSduAB4OfAnm5/SdvMzHqv7G+eAoxExD09T2JmZj3hrhgzs8woIjofWbod+CkQwF9HxOZZxhkFRgEGBwfXjo2NzTqtqakpBgYGqmTuq8kduxhcDjt3L/68h1av7Gi8Xrfd5I5dPZtW2bbrdJl7ob3dernMZcy3vE19T4CzdaNsvpGRka3ddnOXLexPjYg7JR0MfA14V0RcM9f4rVYrJiYmZh02Pj7O8PBwybj9t2bDFtYP7WHTZJVequ5s37iuo/F63XZrNmzp2bTKtl2ny9wL7e3Wy2UuY77lbep7ApytG2XzSeq6sJfqiomIO4v/dwOXAS/sZuZmZtZ7HRd2SSsk7T99G3gFcFO/gpmZWTVl+hsGgcskTT/vwoi4oi+pzMysso4Le0TcBjyvj1nMzKwHfLqjmVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZKF3ZJSyTdKOnyfgQyM7PuVNliPwvY1usgZmbWG6UKu6RDgXXA3/YnjpmZdUsR0fnI0sXAR4D9gT+OiJNmGWcUGAUYHBxcOzY2Nuu0pqamGBgYqJK5ryZ37GJwOezcXXeSuTU5X9lsQ6tX9i/MDO3r3OSOXYs233bzLW9T3xPgbN0om29kZGRrRLS6mefSTkeUdBJwd0RslTQ813gRsRnYDNBqtWJ4ePZRx8fHmWtYnU7fsIX1Q3vYNNlx0yy6Jucrm237qcP9CzND+zp3+oYtizbfdvMtb1PfE+Bs3agjX5mumOOB10jaDowBJ0j6fF9SmZlZZR0X9oj4LxFxaESsAV4HXBkRp/UtmZmZVeLz2M3MMlOpozYixoHxniYxM7Oe8Ba7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8tMx4Vd0r6Svi3pu5JulvShfgYzM7NqyvyC0sPACRExJWkZ8A1J/xQR/9KnbGZmVkHHhT0iApgq7i4r/qIfoczMrDqlet3hyNISYCvwbOC8iPiTWcYZBUYBBgcH146Njc06rampKQYGBuac1+SOXR3n6rXB5bBzd22zX1CT85XNNrR6Zf/CzNC+ztW1fs23vAu9J+r0RM3WhNe5bNuNjIxsjYhWN/MvVdh/8STpQOAy4F0RcdNc47VarZiYmJh12Pj4OMPDw3POY82GLaVz9cr6oT1smqz0O9+Losn5ymbbvnFdH9P8e+3rXF3r13zLu9B7ok5P1GxNeJ3Ltp2krgt7pbNiIuJ+YBw4sZuZm5lZ75U5K+agYksdScuBlwM/6FcwMzOrpsz+/FOAC4p+9r2AL0bE5f2JZWZmVZU5K+Z7wDF9zGJmZj3gK0/NzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMlPmN0+fJukqSdsk3SzprH4GMzOzasr85ukeYH1E3CBpf2CrpK9FxPf7lM3MzCroeIs9Iv4tIm4obj8AbANW9yuYmZlVo4go/yRpDXANcHRE/GzGsFFgFGBwcHDt2NjYrNOYmppiYGBgznlM7thVOlevDC6Hnbtrm/2CmpzP2arrd76h1SsrP3eh92ud5stWVx1pb+uybTcyMrI1IlrdzL90YZc0AFwNnBMRl843bqvViomJiVmHjY+PMzw8POdz12zYUipXL60f2sOmyTK9VIuryfmcrbp+59u+cV3l5y70fq3TfNnqqiPtbV227SR1XdhLnRUjaRlwCfCFhYq6mZnVo8xZMQI+DWyLiD/tXyQzM+tGmS3244E3AidI+k7x96o+5TIzs4o67tCLiG8A6mMWMzPrAV95amaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpaZMr95+hlJd0u6qZ+BzMysO2W22M8HTuxTDjMz65GOC3tEXAPc18csZmbWA4qIzkeW1gCXR8TR84wzCowCDA4Orh0bG5t1vKmpKQYGBuac1+SOXR3n6rXB5bBzd22zX1CT8zlbdf3ON7R6ZeXnLvR+rdN82eqqI+1tXbbtRkZGtkZEq5v597ywt2u1WjExMTHrsPHxcYaHh+d87poNWzrO1Wvrh/awaXJpbfNfSJPzOVt1/c63feO6ys9d6P1ap/my1VVH2tu6bNtJ6rqw+6wYM7PMuLCbmWWmzOmOFwHfAo6QdIekt/YvlpmZVdVxh15EvL6fQczMrDfcFWNmlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWmVKFXdKJkm6R9ENJG/oVyszMqivzm6dLgPOAVwJHAa+XdFS/gpmZWTVltthfCPwwIm6LiEeAMeC3+xPLzMyqUkR0NqJ0MnBiRJxR3H8j8KKIeOeM8UaB0eLuEcAtc0xyFXBPldCLoMnZoNn5nK26JudzturK5jssIg7qZoZLS4yrWR77pU+FiNgMbF5wYtJERLRKzH/RNDkbNDufs1XX5HzOVl0d+cp0xdwBPK3t/qHAnb2NY2Zm3SpT2K8HDpf0DEl7A68D/rE/sczMrKqOu2IiYo+kdwJfAZYAn4mIm7uY94LdNTVqcjZodj5nq67J+ZytukXP1/HBUzMze2LwladmZplxYTczy4wLu5lZZp6QhV3Sakmr684xG0nPlPQeSSfUnWWmJmeDZudztuqanK/J2aB6vidUYZe0RtLVwBXAxyW9pO5M7ST9R+BrpO/SeZukt9cc6ReanA2anc/ZqmtyviZngy7zRUSj/4B9226/FvhEcfvNwD8AQ8V91ZDtBOAZ0/MHPgCcVtx/EfAlYLiOfE3O1vR8zpZnviZn63W+Rm6xSzpA0l9JuhX4hKTDikG/A/youD0G/BA4Y/ppi5jvKEnfA/4b8FlJJ0Rq7aOAQwAi4jrgm8BbFjNfk7M1PZ+z5Zmvydn6la+RhR04EdiXtGCPAB+QtJy0W/JqgIh4GLgYeElx/7F+hZF0qKQD2h46BbgkIl5K+oB5g6TDgQun8xUuA46WtE+/8jU5W9PzOVue+ZqcbbHy1VbYlSyV9FZJX5d0lqRnFYOfDTwSEXuAPwN+CpwGfBV4iqRfK8a7FfixpOP6lPFISV8GvgF8WNL01xT/P2C/4vYXgbuAdaRP1Ce37WHcR/p2y+f9KmVrej5nyzNfk7Mtdr7aCnuxq/GbwJuAjwH7AH9TDL4LuLv4ZPoxaWGeRWqA7/P41wIvA+4tHu8JSSva7j4fuCMi1gBXAp8oHr8PeFjS/hFxH/CvwFOLHN8E3luMtzfwc2B77tmans/Z8szX5Gx15lu0wi7pOEkflXR6cV/AkcAVEfGliPgYcJikFwM7SJ9gRxZP3wYMFI/9BfAqSa8mfSgMAt/tMtuTJJ0v6Xpgo6SDinxDwLWSFBH/CNwvaR1pT2H/YjjF/YOBx0h7GAdL+hvgImBPRNydY7am53O26pqcr8nZmpJvUQq7pOcAfwk8APyepPcW814NPFAsNMD5wBtIhXoP8OLi8RtIR4wfiohrgA3A6cDxwH+PiMfaplHFS4v5vYp0UOJs4ADSl50dUuxdAFxQ5Pt2sSyvBIiIbxXTWBoR24AzgZuB/xkRb6E7Tc7W9HzOlme+JmdrRr65Tpep+kfasj6DtNuxtHjsT4Gzitst4JPAycDLga+0PfdppF0VSIX8RtKvMB0D/G/gKW3jlj4dqWjYM4GrSd05q4rHvwi8u7j9DGBjMfxYUn/YkrZl+0kxndWkPYl3Ap8FPgWs6KLdGput6fmcza+r2+7f//V0i13S80kHOH8b+CDwvmLQDtJvpkL65LkW+F3gn4FDJD1X0rJI/ek7JL0kIq4kfd3lR4FLgYsi4t+m5xVFy5R0EvAa4EPAcaS+fUhn20zvHfwY+Drwyoi4nvSJO1LMcwq4Djg2InYAbyR1Bd0FvC8iHiwbqG1P49VNyzaD266axrUbuO26yfZEaLsyP433SyS9EDgc+GpE/IS0NX5rRJwu6QXAOZJawDjwnyTtFxEPSfou8HukczQvBP4A+KSk3cAkcHsxi78CLoyIXSUyKSJC0rGk3ZyvA1sinR7568BtEXGlpNtJV6++AtgK/I6kVRFxj6R/BR6U9HTgz4HTJB1M+tWoe0m7TkTEBDBRod1apL2aB4CPA3cDz6w7m9uuWrYnQru57fJru/mU2mJXskzSmyTdSOrYPxCYLrw/B7YXW983kHYtjgMe4vFTeAAeJe2CHELaKr+J1L9+NXBPRNwBaau8YlF/KfAZ0lHllwMfKUZ5DLhV0vKIuL3I91zSi3Un6XzS6eVYQmqfS4qMpwJrgc1R8RxXSSslfbaY5u3AuRFxt6S9SJ/kdWZbUrTdb5J2BRvTdsV6NyDpfBrWdsU8Q9IwzVzn9pG0oqFtd0DD225A0r6SLqBhbbegTvprgBXAi4vbBxbBPjnLeGeRLoNdXdw/mdSffhjpKwCuLh7fl9QNs6rtuccAe3eSZ8Y89wPexuNb/suAPwLeUQx/EvC9YvqnkPq71hTDTiqWZVVxexJYSerf/3J7HmCvLrJdRLpibIDUtXRm2zjTxyHeCfyPxcrW9rqeQVrZ1pMO8DSl7aazXVqsVwc1rO32B7aQfkkM4D1NaLcZ+b4M/HVx/2PA2+puO9J74s2k9/8lTWu7tnxXAn9fPNaY9a7TvwW32CWdDdwGbJE0GBH3k/qF7iz6xl+jxy8Q+hbpAOj0hUbXkg6iPhQRFwA/lfQ50kHRW4Bf9CFFxI0R8chCeWZkOwS4HBgGPkc6QPFa0l7CnmK6PyUdeH03qe/rYB4/jfIa0rn0j0TE5cCnSVeznkc6Yv1oW75Sn6ozsv0d8PYi263AEZI2FltRv690wdUVpD2Yvmcr8q0gvblOIF0/8ArScY9jSVtKdbZde7bNpLMFXku6huE36m67wnLStRfPkrSKtM4vKaZZS7vNkm9v0rr2VFIXx9GSPlJX20laRjrGdjLw8Yj43WLQMW3TrK3tZuT7WERMb3FPAkfV2XaldfAJNkzavfhb4D3FY8eSitYdRfALgU3FsHOAD7c9/3rgmOL2PqRTgI7txacSaeV9Udv900lbJm8Gvt32+FOBO4vb7yBdtvuk4vlfAp7eNu6qPmV7E+lI968Df1/8vR74X6Rz+RctW9v0Dmy7/Z9Jb6ZT6267WbL9MemUsWc2qO3eTOprfT/wVtKBtOvrbrdZ8r2PtMezqgltR9oDO3XGY6cA1zWh7ebI9/QiQ+3rXcfL0cGCTp+acwowXtxeRtqaWlncP4y0tX4saRfwYtKW1j+RPqn26Uv41Mcl+MVvt76Ax7t77iWdMzo97tcoCi1p9+mrxTh/skjZjgG+Mb3ito23jHRw+YTi/jn9zjYj5wGk4xs7gQ8X9+8FButqu1my3VXMdwVFN19dbdf2er6F1M32WuALxWP31N1uc+QbKx5rP124lvWO1EVxK7CpmP8HivpxH3BwA9a59nxXkb6Y69C617vSy1FigZ9MulDoOcX9pTOGnw+cPL0CkboezqRPRX2OlfkCHj9f/nPAR4vbv0ba43h62wtzNG1fCbxI2d7R/lhx+5Ci7Z672NnaMvwh6XzbzaR+7W8WbzjV2XYzsp1HOq3s2U1oO9JXRi8h9aFeTdoyvgl4f93r3Cz5/pl0htkLGtJ2XyHtgT2NtBV8FmnDsCnrXHu+L5Au/T+8CW3X6d900emIpE8BP4uIDcX9vUjnXb4DeA5wSpTsJ+8VSYeS+rTeFRG3Kn2h2GiRazXwnejNVWXdZHt7RNxWPHYMqVtqXZHtD+vI1k7pOoQzSW+yI0kr66HU2HZt2Y4mvdn+nHSW1UnU1HaSBkjdHPuQ2uk3SBeenE3aUj6cGtttlnyHk45P/BbpmNfLSO1Xy3qn4rTn4vbzSO/Ta0mX1Ne+zs3IdzTpSvdzSd80W9t6V0bZ89g3A+cWBxmOJK3Ex5NelLPrKuqFYyjOgZd0Bqn//2xSF9IPIp1+WXe2HxXZbietHHtIW/E31pit3b2kg4Dvi4i/k3QacHND8t1P6ie+ifS6LqO+tttDOnviUdKW+s9J6/8k8N4GtNtc+R6W9BpSwa9tvZsumoX7Sced3h8RFzag7Wbme4C08boN+K/Uu951rOwW++tIB0ofJn3j2JURcUufspUi6VrSwbXtpHNIPxQR36s1VGFGtruADQ1qt5WkLbg3kL7/fjNwXkQ8Ou8TF8Es2T4dEZvqTfXLigtPpvuy76o7z0xFvpOBz0Y666TuPPuQfnPhjaQ96r8EPhXpa7prN0u+zRHxZ/WmKqfjwi7puaTzOS8mHSzq2VfldqvYg/ggaUv485GuWmuEJmcDkLSU1P3yMClfk17XxmaDdFEX8FiU2TpaRE3OJ+lM0mm1n2va6wrNz7eQUlvsZmbWfE39aTwzM6vIhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlpn/D0QBdzhJVkBDAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
"source": [
"yearly_incidence.hist(xrot=20)"
]
@@ -364,7 +2501,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.6.1"
+ "version": "3.6.4"
}
},
"nbformat": 4,