From 15134480d0a35a27808c7b9c21dcde9da46a818f Mon Sep 17 00:00:00 2001 From: 3177f1da89e77b25849339303ddf1357 <3177f1da89e77b25849339303ddf1357@app-learninglab.inria.fr> Date: Mon, 1 Feb 2021 12:53:32 +0000 Subject: [PATCH] no commit message --- module3/exo1/analyse-syndrome-grippal.ipynb | 2212 ++++++++++++++++++- 1 file changed, 2174 insertions(+), 38 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5..a9aadfa 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": 1, "metadata": {}, "outputs": [], "source": [ @@ -28,13 +28,11 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 2, + "metadata": {}, "outputs": [], "source": [ - "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" + "data_url = \"./incidence-PAY-3.csv\"" ] }, { @@ -61,9 +59,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020210332145417480.025428.03226.038.0FRFrance
120210231784114356.021326.02722.032.0FRFrance
220210132180917786.025832.03327.039.0FRFrance
320205332122016498.025942.03225.039.0FRFrance
420205231642812285.020571.02519.031.0FRFrance
520205132161917370.025868.03327.039.0FRFrance
620205031684513220.020470.02620.032.0FRFrance
72020493129399923.015955.02015.025.0FRFrance
820204831380410641.016967.02116.026.0FRFrance
920204731908515285.022885.02923.035.0FRFrance
1020204632480120503.029099.03831.045.0FRFrance
1120204534251636857.048175.06556.074.0FRFrance
1220204434456738521.050613.06859.077.0FRFrance
1320204334373737523.049951.06657.075.0FRFrance
1420204233514529812.040478.05345.061.0FRFrance
1520204132787723206.032548.04235.049.0FRFrance
1620204032044316381.024505.03125.037.0FRFrance
1720203931981015900.023720.03024.036.0FRFrance
1820203832556221142.029982.03932.046.0FRFrance
1920203731848514649.022321.02822.034.0FRFrance
202020363103907646.013134.01612.020.0FRFrance
21202035399186842.012994.01510.020.0FRFrance
22202034360843090.09078.094.014.0FRFrance
23202033361063411.08801.095.013.0FRFrance
24202032359183330.08506.095.013.0FRFrance
25202031343512269.06433.074.010.0FRFrance
26202030381795442.010916.0128.016.0FRFrance
27202029386875860.011514.0139.017.0FRFrance
28202028383405701.010979.0139.017.0FRFrance
29202027340662406.05726.063.09.0FRFrance
.................................
186119852132609619621.032571.04735.059.0FRFrance
186219852032789620885.034907.05138.064.0FRFrance
186319851934315432821.053487.07859.097.0FRFrance
186419851834055529935.051175.07455.093.0FRFrance
186519851733405324366.043740.06244.080.0FRFrance
186619851635036236451.064273.09166.0116.0FRFrance
186719851536388145538.082224.011683.0149.0FRFrance
18681985143134545114400.0154690.0244207.0281.0FRFrance
18691985133197206176080.0218332.0357319.0395.0FRFrance
18701985123245240223304.0267176.0445405.0485.0FRFrance
18711985113276205252399.0300011.0501458.0544.0FRFrance
18721985103353231326279.0380183.0640591.0689.0FRFrance
18731985093369895341109.0398681.0670618.0722.0FRFrance
18741985083389886359529.0420243.0707652.0762.0FRFrance
18751985073471852432599.0511105.0855784.0926.0FRFrance
18761985063565825518011.0613639.01026939.01113.0FRFrance
18771985053637302592795.0681809.011551074.01236.0FRFrance
18781985043424937390794.0459080.0770708.0832.0FRFrance
18791985033213901174689.0253113.0388317.0459.0FRFrance
188019850239758680949.0114223.0177147.0207.0FRFrance
188119850138548965918.0105060.0155120.0190.0FRFrance
188219845238483060602.0109058.0154110.0198.0FRFrance
1883198451310172680242.0123210.0185146.0224.0FRFrance
18841984503123680101401.0145959.0225184.0266.0FRFrance
1885198449310107381684.0120462.0184149.0219.0FRFrance
188619844837862060634.096606.0143110.0176.0FRFrance
188719844737202954274.089784.013199.0163.0FRFrance
188819844638733067686.0106974.0159123.0195.0FRFrance
18891984453135223101414.0169032.0246184.0308.0FRFrance
189019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1891 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202103 3 21454 17480.0 25428.0 32 26.0 \n", + "1 202102 3 17841 14356.0 21326.0 27 22.0 \n", + "2 202101 3 21809 17786.0 25832.0 33 27.0 \n", + "3 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "4 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "5 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "6 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "7 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "8 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "9 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "10 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "11 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "12 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "13 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "14 202042 3 35145 29812.0 40478.0 53 45.0 \n", + "15 202041 3 27877 23206.0 32548.0 42 35.0 \n", + "16 202040 3 20443 16381.0 24505.0 31 25.0 \n", + "17 202039 3 19810 15900.0 23720.0 30 24.0 \n", + "18 202038 3 25562 21142.0 29982.0 39 32.0 \n", + "19 202037 3 18485 14649.0 22321.0 28 22.0 \n", + "20 202036 3 10390 7646.0 13134.0 16 12.0 \n", + "21 202035 3 9918 6842.0 12994.0 15 10.0 \n", + "22 202034 3 6084 3090.0 9078.0 9 4.0 \n", + "23 202033 3 6106 3411.0 8801.0 9 5.0 \n", + "24 202032 3 5918 3330.0 8506.0 9 5.0 \n", + "25 202031 3 4351 2269.0 6433.0 7 4.0 \n", + "26 202030 3 8179 5442.0 10916.0 12 8.0 \n", + "27 202029 3 8687 5860.0 11514.0 13 9.0 \n", + "28 202028 3 8340 5701.0 10979.0 13 9.0 \n", + "29 202027 3 4066 2406.0 5726.0 6 3.0 \n", + "... ... ... ... ... ... ... ... \n", + "1861 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1862 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1863 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1864 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1865 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1866 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1867 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1868 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1869 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1870 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1871 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1872 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1873 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1874 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1875 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1876 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1877 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1878 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1879 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1880 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1881 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1882 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1883 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1884 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1885 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1886 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1887 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1888 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1889 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1890 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 38.0 FR France \n", + "1 32.0 FR France \n", + "2 39.0 FR France \n", + "3 39.0 FR France \n", + "4 31.0 FR France \n", + "5 39.0 FR France \n", + "6 32.0 FR France \n", + "7 25.0 FR France \n", + "8 26.0 FR France \n", + "9 35.0 FR France \n", + "10 45.0 FR France \n", + "11 74.0 FR France \n", + "12 77.0 FR France \n", + "13 75.0 FR France \n", + "14 61.0 FR France \n", + "15 49.0 FR France \n", + "16 37.0 FR France \n", + "17 36.0 FR France \n", + "18 46.0 FR France \n", + "19 34.0 FR France \n", + "20 20.0 FR France \n", + "21 20.0 FR France \n", + "22 14.0 FR France \n", + "23 13.0 FR France \n", + "24 13.0 FR France \n", + "25 10.0 FR France \n", + "26 16.0 FR France \n", + "27 17.0 FR France \n", + "28 17.0 FR France \n", + "29 9.0 FR France \n", + "... ... ... ... \n", + "1861 59.0 FR France \n", + "1862 64.0 FR France \n", + "1863 97.0 FR France \n", + "1864 93.0 FR France \n", + "1865 80.0 FR France \n", + "1866 116.0 FR France \n", + "1867 149.0 FR France \n", + "1868 281.0 FR France \n", + "1869 395.0 FR France \n", + "1870 485.0 FR France \n", + "1871 544.0 FR France \n", + "1872 689.0 FR France \n", + "1873 722.0 FR France \n", + "1874 762.0 FR France \n", + "1875 926.0 FR France \n", + "1876 1113.0 FR France \n", + "1877 1236.0 FR France \n", + "1878 832.0 FR France \n", + "1879 459.0 FR France \n", + "1880 207.0 FR France \n", + "1881 190.0 FR France \n", + "1882 198.0 FR France \n", + "1883 224.0 FR France \n", + "1884 266.0 FR France \n", + "1885 219.0 FR France \n", + "1886 176.0 FR France \n", + "1887 163.0 FR France \n", + "1888 195.0 FR France \n", + "1889 308.0 FR France \n", + "1890 213.0 FR France \n", + "\n", + "[1891 rows x 10 columns]" + ] + }, + "execution_count": 3, + "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": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
165419891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1654 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1654 FR France " + ] + }, + "execution_count": 4, + "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": 5, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020210332145417480.025428.03226.038.0FRFrance
120210231784114356.021326.02722.032.0FRFrance
220210132180917786.025832.03327.039.0FRFrance
320205332122016498.025942.03225.039.0FRFrance
420205231642812285.020571.02519.031.0FRFrance
520205132161917370.025868.03327.039.0FRFrance
620205031684513220.020470.02620.032.0FRFrance
72020493129399923.015955.02015.025.0FRFrance
820204831380410641.016967.02116.026.0FRFrance
920204731908515285.022885.02923.035.0FRFrance
1020204632480120503.029099.03831.045.0FRFrance
1120204534251636857.048175.06556.074.0FRFrance
1220204434456738521.050613.06859.077.0FRFrance
1320204334373737523.049951.06657.075.0FRFrance
1420204233514529812.040478.05345.061.0FRFrance
1520204132787723206.032548.04235.049.0FRFrance
1620204032044316381.024505.03125.037.0FRFrance
1720203931981015900.023720.03024.036.0FRFrance
1820203832556221142.029982.03932.046.0FRFrance
1920203731848514649.022321.02822.034.0FRFrance
202020363103907646.013134.01612.020.0FRFrance
21202035399186842.012994.01510.020.0FRFrance
22202034360843090.09078.094.014.0FRFrance
23202033361063411.08801.095.013.0FRFrance
24202032359183330.08506.095.013.0FRFrance
25202031343512269.06433.074.010.0FRFrance
26202030381795442.010916.0128.016.0FRFrance
27202029386875860.011514.0139.017.0FRFrance
28202028383405701.010979.0139.017.0FRFrance
29202027340662406.05726.063.09.0FRFrance
.................................
186119852132609619621.032571.04735.059.0FRFrance
186219852032789620885.034907.05138.064.0FRFrance
186319851934315432821.053487.07859.097.0FRFrance
186419851834055529935.051175.07455.093.0FRFrance
186519851733405324366.043740.06244.080.0FRFrance
186619851635036236451.064273.09166.0116.0FRFrance
186719851536388145538.082224.011683.0149.0FRFrance
18681985143134545114400.0154690.0244207.0281.0FRFrance
18691985133197206176080.0218332.0357319.0395.0FRFrance
18701985123245240223304.0267176.0445405.0485.0FRFrance
18711985113276205252399.0300011.0501458.0544.0FRFrance
18721985103353231326279.0380183.0640591.0689.0FRFrance
18731985093369895341109.0398681.0670618.0722.0FRFrance
18741985083389886359529.0420243.0707652.0762.0FRFrance
18751985073471852432599.0511105.0855784.0926.0FRFrance
18761985063565825518011.0613639.01026939.01113.0FRFrance
18771985053637302592795.0681809.011551074.01236.0FRFrance
18781985043424937390794.0459080.0770708.0832.0FRFrance
18791985033213901174689.0253113.0388317.0459.0FRFrance
188019850239758680949.0114223.0177147.0207.0FRFrance
188119850138548965918.0105060.0155120.0190.0FRFrance
188219845238483060602.0109058.0154110.0198.0FRFrance
1883198451310172680242.0123210.0185146.0224.0FRFrance
18841984503123680101401.0145959.0225184.0266.0FRFrance
1885198449310107381684.0120462.0184149.0219.0FRFrance
188619844837862060634.096606.0143110.0176.0FRFrance
188719844737202954274.089784.013199.0163.0FRFrance
188819844638733067686.0106974.0159123.0195.0FRFrance
18891984453135223101414.0169032.0246184.0308.0FRFrance
189019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1890 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202103 3 21454 17480.0 25428.0 32 26.0 \n", + "1 202102 3 17841 14356.0 21326.0 27 22.0 \n", + "2 202101 3 21809 17786.0 25832.0 33 27.0 \n", + "3 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "4 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "5 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "6 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "7 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "8 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "9 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "10 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "11 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "12 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "13 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "14 202042 3 35145 29812.0 40478.0 53 45.0 \n", + "15 202041 3 27877 23206.0 32548.0 42 35.0 \n", + "16 202040 3 20443 16381.0 24505.0 31 25.0 \n", + "17 202039 3 19810 15900.0 23720.0 30 24.0 \n", + "18 202038 3 25562 21142.0 29982.0 39 32.0 \n", + "19 202037 3 18485 14649.0 22321.0 28 22.0 \n", + "20 202036 3 10390 7646.0 13134.0 16 12.0 \n", + "21 202035 3 9918 6842.0 12994.0 15 10.0 \n", + "22 202034 3 6084 3090.0 9078.0 9 4.0 \n", + "23 202033 3 6106 3411.0 8801.0 9 5.0 \n", + "24 202032 3 5918 3330.0 8506.0 9 5.0 \n", + "25 202031 3 4351 2269.0 6433.0 7 4.0 \n", + "26 202030 3 8179 5442.0 10916.0 12 8.0 \n", + "27 202029 3 8687 5860.0 11514.0 13 9.0 \n", + "28 202028 3 8340 5701.0 10979.0 13 9.0 \n", + "29 202027 3 4066 2406.0 5726.0 6 3.0 \n", + "... ... ... ... ... ... ... ... \n", + "1861 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1862 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1863 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1864 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1865 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1866 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1867 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1868 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1869 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1870 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1871 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1872 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1873 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1874 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1875 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1876 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1877 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1878 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1879 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1880 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1881 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1882 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1883 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1884 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1885 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1886 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1887 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1888 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1889 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1890 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 38.0 FR France \n", + "1 32.0 FR France \n", + "2 39.0 FR France \n", + "3 39.0 FR France \n", + "4 31.0 FR France \n", + "5 39.0 FR France \n", + "6 32.0 FR France \n", + "7 25.0 FR France \n", + "8 26.0 FR France \n", + "9 35.0 FR France \n", + "10 45.0 FR France \n", + "11 74.0 FR France \n", + "12 77.0 FR France \n", + "13 75.0 FR France \n", + "14 61.0 FR France \n", + "15 49.0 FR France \n", + "16 37.0 FR France \n", + "17 36.0 FR France \n", + "18 46.0 FR France \n", + "19 34.0 FR France \n", + "20 20.0 FR France \n", + "21 20.0 FR France \n", + "22 14.0 FR France \n", + "23 13.0 FR France \n", + "24 13.0 FR France \n", + "25 10.0 FR France \n", + "26 16.0 FR France \n", + "27 17.0 FR France \n", + "28 17.0 FR France \n", + "29 9.0 FR France \n", + "... ... ... ... \n", + "1861 59.0 FR France \n", + "1862 64.0 FR France \n", + "1863 97.0 FR France \n", + "1864 93.0 FR France \n", + "1865 80.0 FR France \n", + "1866 116.0 FR France \n", + "1867 149.0 FR France \n", + "1868 281.0 FR France \n", + "1869 395.0 FR France \n", + "1870 485.0 FR France \n", + "1871 544.0 FR France \n", + "1872 689.0 FR France \n", + "1873 722.0 FR France \n", + "1874 762.0 FR France \n", + "1875 926.0 FR France \n", + "1876 1113.0 FR France \n", + "1877 1236.0 FR France \n", + "1878 832.0 FR France \n", + "1879 459.0 FR France \n", + "1880 207.0 FR France \n", + "1881 190.0 FR France \n", + "1882 198.0 FR France \n", + "1883 224.0 FR France \n", + "1884 266.0 FR France \n", + "1885 219.0 FR France \n", + "1886 176.0 FR France \n", + "1887 163.0 FR France \n", + "1888 195.0 FR France \n", + "1889 308.0 FR France \n", + "1890 213.0 FR France \n", + "\n", + "[1890 rows x 10 columns]" + ] + }, + "execution_count": 5, + "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": 6, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "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": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXm8XEWZ//95ertrbvYESIAk7AFEIIYgowgMm84ILvyGcSEqI+qXcZz5Oio4+sVRUXR0dBwVRUUBlV0lshrDvoawhJCEkJCQfb/JvTf35i7dXb8/TtXpc7qrTp3uPr2c3Of9et1X9z1LnafrnKqnnqXqkBACDMMwDFMtiUYLwDAMwxwYsEJhGIZhIoEVCsMwDBMJrFAYhmGYSGCFwjAMw0QCKxSGYRgmElihMAzDMJHACoVhGIaJBFYoDMMwTCSkGi1APZk0aZKYMWNGo8VgGIaJFS+88MIuIcRk23GjSqHMmDEDS5YsabQYDMMwsYKI1oc5jl1eDMMwTCSwQmEYhmEigRUKwzAMEwmsUBiGYZhIYIXCMAzDRIJVoRDRjUS0g4he9WybQEQLiWi1/Bzv2Xc1Ea0holVEdL5n+6lEtEzu+xERkdzeQkS3y+3PEdEMzznz5TVWE9F8z/aZ8tjV8txM9VXBMAzDVEMYC+U3AC4o2nYVgEVCiKMALJL/g4hmA7gUwPHynJ8SUVKecz2AKwAcJf9UmZcD2COEOBLADwB8R5Y1AcA1AE4DMBfANR7F9R0AP5DX3yPLYBiGYRqIVaEIIR4H0F20+SIAN8nvNwG42LP9NiHEkBBiHYA1AOYS0cEAuoQQzwjnncM3F52jyroLwDnSejkfwEIhRLcQYg+AhQAukPvOlscWX/+A589Lt6BnYKTRYjAMw5RQaQxlqhBiKwDIzyly+zQAGz3HbZLbpsnvxdt95wghsgB6AEwMKGsigL3y2OKyDmje3NWPz976Ev719pcaLQrDMEwJUQflSbNNBGyv5JygskoFIrqCiJYQ0ZKdO3eaDosF+4YcHbqtd6jBkjAMw5RSqULZLt1YkJ875PZNAA71HDcdwBa5fbpmu+8cIkoBGAvHxWYqaxeAcfLY4rJKEELcIISYI4SYM3mydSmapiaXd/RmOqnTqQzDMI2lUoWyAIDKupoP4B7P9ktl5tZMOMH3xdIt1kdE82QM5LKic1RZHwTwsIyzPATgPCIaL4Px5wF4SO57RB5bfP0DmqxUKMkEKxSGYZoP6+KQRHQrgHcBmEREm+BkXl0H4A4iuhzABgCXAIAQYjkR3QFgBYAsgCuFEDlZ1GfgZIy1AXhA/gHArwDcQkRr4Fgml8qyuonoGwCel8d9XQihkgO+BOA2IvomgJdkGQc82VweAJBihcIwTBNiVShCiH807DrHcPy1AK7VbF8C4ATN9kFIhaTZdyOAGzXb18JJJR5VKJdXKsHzURmGaT64Z4oRI0qhcAyFYZgmhBVKjMhLhZIgVigMwzQfrFDihNQjxhxphmGYBsIKJUYou8RJdGMYhmkuWKHECGJXF8MwTQwrlBhRsFAaKgbDMIwWVigxgtwYCmsUhmGaD1YoMYKkjcIWCsMwzQgrlBjCCoVhmGaEFUqMUK4udnkxDNOMsEKJIWyhMAzTjLBCiRFKkbA+YRimGWGFEkdYozAM04SwQokRwv1kjcIwTPPBCiWGNHsMZfX2Psy46j4sebPbfjDDMAcMrFBihFrDK9/kGuWJ1bsAAPe+srXBkjAMU09YocSI5lYjBXjJMYYZnbBCiSFxUSwMw4wuWKHECZU2HBONwsvsM8zoghVKDOGOmmGYZoQVSozgdGGGYZoZVihM5HBMnmFGJ6xQYgR7uhiGaWZYocSIuCmUmInLMEyVsEJhIod4IgrDjEpYocQIHvEzDNPMsEJhakbcXHQMw1QHK5QYEZf5J+zxYpjRCSuUGBEPdcIwzGiFFcoBxo1PrsPZ33u00WIA4ImYDDPaqEqhENG/EdFyInqViG4lolYimkBEC4lotfwc7zn+aiJaQ0SriOh8z/ZTiWiZ3PcjkmlCRNRCRLfL7c8R0QzPOfPlNVYT0fxqfkdcCOPx+vq9K7B2V3/thQmAPV4MMzqpWKEQ0TQA/wJgjhDiBABJAJcCuArAIiHEUQAWyf9BRLPl/uMBXADgp0SUlMVdD+AKAEfJvwvk9ssB7BFCHAngBwC+I8uaAOAaAKcBmAvgGq/iYpqDmIR8GIaJiGpdXikAbUSUAtAOYAuAiwDcJPffBOBi+f0iALcJIYaEEOsArAEwl4gOBtAlhHhGOFHnm4vOUWXdBeAcab2cD2ChEKJbCLEHwEIUlNABTEx6aI7KM8yopGKFIoTYDOB7ADYA2AqgRwjxFwBThRBb5TFbAUyRp0wDsNFTxCa5bZr8Xrzdd44QIgugB8DEgLIYhmGYBlGNy2s8HAtiJoBDAHQQ0UeCTtFsEwHbKz2nWM4riGgJES3ZuXNngHjNT9xcSDETl2GYKqnG5fW3ANYJIXYKIUYA/AHA2wFsl24syM8d8vhNAA71nD8djotsk/xevN13jnSrjQXQHVBWCUKIG4QQc4QQcyZPnlzhT20O4tJBs8OLYUYn1SiUDQDmEVG7jGucA2AlgAUAVNbVfAD3yO8LAFwqM7dmwgm+L5ZusT4imifLuazoHFXWBwE8LOMsDwE4j4jGS0vpPLmNYRiGaRCpSk8UQjxHRHcBeBFAFsBLAG4A0AngDiK6HI7SuUQev5yI7gCwQh5/pRAiJ4v7DIDfAGgD8ID8A4BfAbiFiNbAsUwulWV1E9E3ADwvj/u6EKK70t8SF2Ln8oqZvAzDVEfFCgUAhBDXwEnf9TIEx1rRHX8tgGs125cAOEGzfRBSIWn23QjgxjJFZhiGYWoEz5SPETzznGGYZoYVSgwJo1YauZAkT0NhmNEJK5QYUY6O4PgFwzD1hhVKjAianGM6trE0hxQMw9QHVigxJEw3nW+kyyuCmSj7h3P4yC+fw5od+yKQiGGYesAKJUaUExeJu8vr6Td24ck1u3DtfSsaLQrDMCFhhXKA0gwZYVEoNeIIP8PEBlYoByiNtFBYBzDM6IQVygFK3F1eDMPED1YoMaKstOEDxOXFMEx8YIUSI8pREg11eUVQBisjhokfrFAOUJqhP67GSipnzg3DMM0BK5QYUc6ovaHzUCLUAhzgZ5j4wArlAIVdRgzD1BtWKDGiLCXRBAqFlRrDjC5YocSI8vRJvJdeaeRqyQzDVAYrlAOUA6c/5iAKw8QFVigxopxReyOD8gzDjE5YoRygxF2dxF1+hhmNsEKJEWXFUA6QHpnThhkmPrBCOUBpiqVXqjm38eIzDFMmrFDiRFxeARzlxMboimIYpsawQokRcVnLq5lkYBimfrBCOUBp7DwUhmFGI6xQYkRZy9fH3jqI/Q9gmFEHK5QYEkZZNMM8lKpWG5anBmV59ewfQd/gSMXXYBgmWlihlMGKLb148NVtDbt+XNKGo3wPfNAyLif9519w4tf+Etm1GIapjlSjBYgT7/7REwCAN697T4MlYRiGaT7YQokRsYuhNIMMDMPUDVYoMULFJMJ4lDjLi2GYesMKJYaEC8rXXg4bVc2Uj0wKhmHqRVUKhYjGEdFdRPQaEa0kotOJaAIRLSSi1fJzvOf4q4loDRGtIqLzPdtPJaJlct+PSEZ1iaiFiG6X258johmec+bLa6wmovnV/I64UJ7LK96vAA6T5cUwTHNRrYXyPwAeFEIcC+AkACsBXAVgkRDiKACL5P8gotkALgVwPIALAPyUiJKynOsBXAHgKPl3gdx+OYA9QogjAfwAwHdkWRMAXAPgNABzAVzjVVzMgTPCZ4XCMPGhYoVCRF0A3gngVwAghBgWQuwFcBGAm+RhNwG4WH6/CMBtQoghIcQ6AGsAzCWigwF0CSGeEc6w+uaic1RZdwE4R1ov5wNYKIToFkLsAbAQBSV0wBKXtOGCDE0gBMMwdaMaC2UWgJ0Afk1ELxHRL4moA8BUIcRWAJCfU+Tx0wBs9Jy/SW6bJr8Xb/edI4TIAugBMDGgLEbCnTnDMPWmGoWSAnAKgOuFECcD6Id0bxnQOS9EwPZKz/FflOgKIlpCREt27twZIF4MKENJNFKdRKHLmmH5fYZhyqMahbIJwCYhxHPy/7vgKJjt0o0F+bnDc/yhnvOnA9git0/XbPedQ0QpAGMBdAeUVYIQ4gYhxBwhxJzJkydX8DO1ZUZSTtnXLefYA6Q/DpopzzBMc1GxQhFCbAOwkYiOkZvOAbACwAIAKutqPoB75PcFAC6VmVsz4QTfF0u3WB8RzZPxkcuKzlFlfRDAwzLO8hCA84hovAzGnye31YVmSMm10cgRfgyqh2GYGlDt0iufBfA7IsoAWAvg43CU1B1EdDmADQAuAQAhxHIiugOO0skCuFIIkZPlfAbAbwC0AXhA/gFOwP8WIloDxzK5VJbVTUTfAPC8PO7rQojuKn9LaHJ5gWSi/iPncqyOfL52ctiIwoI7UCwshhlNVKVQhBAvA5ij2XWO4fhrAVyr2b4EwAma7YOQCkmz70YAN5Yjb1TkYmCiNEMMIhIJ2OPFMLGBZ8pXQK5RMZRygvIN1CdRXLrx6jB+PPb6Tsy46j5s2jPQaFGYUQorlApolIUSt042CqXGBkp47ljiZNK/uGFvgyVhRiusUCogHweXV0Pzhht47QOQkVwef1m+zWqhsvJlGg0rlAponMsr/LGNfGNjM8RvDiR+/PAaXHHLC3h0Vbh5VDyplWkUrFAqIBYWSqMFQHPIcCCwac9+AMDOfUOBx0X5psxqeOS1HRjK5uwHMgccrFAqoGEWSjnHNtJCiWKmvFDvfmlsJzmUzeH/3v4ytuzd3zAZ3CqIgYZevK4bH//N8/j+X15vtChMA2CFUgGNThsO41I6UEIojR5zP/LaTvzhpc24ZsHyhsmgpjw10o0Zlm29gwCAzQ1UwEzjYIVSAQ3L8iorbbjxnU8zyFAtLWmniQyONM6Fo5afsdVmo5UvAAxnnRm1Lcn4dy23PLseM666D72DI40WJTbE/643gEZbKGHWt2roPJSY6ZH9wzns6BvU7mtNOa/sGco2bukB5fIKW6+NrH9XoaTj37X87tn1AIBN3WxthSX+d70BxMH10FiXV/PXj5d/uOEZzL12kXaf6hiHGmihKGzPXTPE5FUwPnMAWCjtGWcwMTCcbbAk8SH+d70B5Bo0WC3vFcC1kyO0DNWcW0f5X9nUY9ynOsZGWiiKREiN0UiFrqz3RAPWuoua9oyzMlX/cOMHE3GBFUoFNNrlFabDaOg8lAgv3QyjbqCx9zzvZrwFH9ckVXXAoBaAjYNHollghVIBjXrAyhl5NkUbqEKGZnGbqXtdrUJ5eeNeY5zGhrJM4hBDaXSad01ojkcxFrBCqYBGWyhhiPv7UNTy+43unlTnPFLl+wAu/slTuPCHT1R0rhuUt9RsFJ35xu4BPLVmV9XlNMWApkoORN1Ya6p9H8qoJA5LrzTFO4CraJDN0h8pOXK56iXa3T9c0XlKUYQdx1TzeJ75X48gL4A3r3tPRedzHzy6YQulAoIslN8/twEvb6zNaq9x0Scu1bi8mmSIq1xe2QZapW58uw6LQ8bA+A7FK5v24qT//At2W5arCUOzuF/jACuUMlBBuiCF8uU/LsPFP3mqXiIZaezikNGV0WifvKrGRro51byj0BZKDWWxUZgz09hO+IbH16Jn/wierMJ9x9ZW+bBCKQN3CYyGzZSvzbFRE8m1m2ZQ2DwWinWQ0AQ9YLMsO9aalhNSRxqf7j2aYIVSBmq03KgYSjk0g4TVuAqaxc2g9EgjR9zucxdSqTWDrCYRNnYP1GUpk1a1ZE4Eqx7HoLk3DaxQykCNFBv3xsbw123sPJTqr62quNGDbvVTGtmnJCjcfIgwS/LUGpuH8h3ffQSX/vzZmsuhlszpGahcedmUI1MKK5QyCNuwm4ImELGaamqWKm50LAAodNJxiKEodIOf/iFnCZMVW3trfv0ZkzoAAOu7B6ouKxbtvUlghVIGCdf10Jjrx+eNjVGUEWyi1KujdzvxRk4WlJ82y7gZ5k0EiTBcx+VrMimna4viMTlQMt/qASuUMiCLy6sZRrOKuAflbWXU6/cpxdZQl5daAiRGJoru/sRtpF/I1o6X3I2EFUoZJCzB0WaaQd8Mjbcql1eV+yu6pk7gSJRjdYW4A5k6zENRVCxzgJlUz+ah5I8iuaOJmnXTwwqlDNx5KIbGVuvU0rJesFVDOepybctvrcWoUT+qrv56UYka3kBpfCeqO72emXtRPh7NMDiLC6xQysA2D6UZHrxmmFgWxbXdiY2GcXdNLBTttuh+S7UF2Oo1yhhKpc9ykAj1fCSjiH1R2Pk/jAsrlDKwzQeotcsrLhMbXRmqaM1KaZs6yVo0cl2HHU1Qt7pC3L4xrIXSDDI3OIYSZeyrGdpSXGCFUgYJiy+75golxDGq/42739caQ6nB79MVmRfVd0zVyhrW4otyHkqliysXltqvjXIOS5QTUtlCCQ8rlDJw56E0cVDenYzVyOXr65DlVQt0HUe51kHYcsuhMLkyXDlRVF3FLq8AnVbXjjmSa5W3hhrDCqUsEpalV2q9JEuY4pvJQqlpllctLBRdklcU8aAqiwh7Lwvxs+qu51wzepdXQ2IokZTVBI0pJlStUIgoSUQvEdG98v8JRLSQiFbLz/GeY68mojVEtIqIzvdsP5WIlsl9PyI5zCaiFiK6XW5/johmeM6ZL6+xmojmV/s7wv1W5zMOacMNDcpHEci2yF+vRl6udaAto8r6cOMBMcjyCnK61bNjjsJV2QwJLnEjCgvlcwBWev6/CsAiIcRRABbJ/0FEswFcCuB4ABcA+CkRJeU51wO4AsBR8u8Cuf1yAHuEEEcC+AGA78iyJgC4BsBpAOYCuMaruGqFbfn62sdQ7OVHOUqtluicDrUpuxityyuS0X5154eVodwlWoKodEXtQnJVY2Mo0aYNR1fWgU5VCoWIpgN4D4BfejZfBOAm+f0mABd7tt8mhBgSQqwDsAbAXCI6GECXEOIZ4QwFbi46R5V1F4BzpPVyPoCFQohuIcQeAAtRUEI1wzaxscq3xIYmTGM50GMoUY0aveUEZSZVty5ZRLLW8XqVWhNBpzXEQmkC999ooloL5YcAvgjA25VOFUJsBQD5OUVunwZgo+e4TXLbNPm9eLvvHCFEFkAPgIkBZdUU23spsjXWKOFiKI6Q//HHVxv33hZXlspRdWwK8kb1y7x1qitTBOwLS/UWSnmdYxT3vdIigtxzcRvpN1M8Mi5UrFCI6O8A7BBCvBD2FM02EbC90nP8FyW6goiWENGSnTt3hhLUhG1xyFqPZMopfWA4h3W7+2smSxiqSrW17Y9Id3vvWa1cXtVaDOWeHYnLqwYWSj1jEVHEUBSNjqE8/Np2PP56dX1XvajGQjkDwHuJ6E0AtwE4m4h+C2C7dGNBfu6Qx28CcKjn/OkAtsjt0zXbfecQUQrAWADdAWWVIIS4QQgxRwgxZ/LkyZX9Uolt+fp6rUIc+Hw3wYqzdXF5RWSjeDvfwCyvqlxelZ/rPd/+m6N7vULFCqXoU7evHkQaQ2mwifKJ3yzBZTcubqgMYalYoQghrhZCTBdCzIATbH9YCPERAAsAqKyr+QDukd8XALhUZm7NhBN8XyzdYn1ENE/GRy4rOkeV9UF5DQHgIQDnEdF4GYw/T26rKQ3P8oqJLzea5Uqky8u09EpEVeGTVadQdMeVSbUdfCN8+LVIJKhvDEXJUfk1o0xyGC2kalDmdQDuIKLLAWwAcAkACCGWE9EdAFYAyAK4Ugih3s/5GQC/AdAG4AH5BwC/AnALEa2BY5lcKsvqJqJvAHheHvd1IUR3DX6Lj0YvvaIIa6BUaqy8uGEP2tJJHHdwV4UlOFQXyLbsr7xo43V0HZ4by6nC9KtWVlHyJZhILJSAZ3nTngFMH9+u3+m6mjR1Wcf3CEWZlMJB+fBEolCEEI8CeFR+3w3gHMNx1wK4VrN9CYATNNsHIRWSZt+NAG6sVOZqMCqUJoihRLFA4Pt/+jQA4M3r3lPR+fVof1E1cm85uhIbGY9QNCRt2HDRZ9fuxqU3PIsf/MNJeN/J00v2Bym/Rqw2HE0MJYJCRgk8U74MlPncqLW8iuXQ0QzvFS9QeX3YXBVRNXJ/DEUXlI/AfVd1EeUFmKOJoei3r5Sv7315w17t/uCgfLVShSeKuAdFGJMaLbBCqYChEb3t3kyrDQMFF10cUVVpThuO3kLR3b4oOpOogvI2Cm8YrO56gPl3F+6LKbZlVn71XW04OjiGEh5WKGWg2kP/UFa7v5lmyjeSKEf1xt8TVVDeMzaold8/ssUhQ5YTxejcdC21PWFSKAFl1rNjzkfo82ILJTysUMpAdTj9w3qFUu/1pZqd6haHtLi8Ki/ah++e1cpCqfp8Oeq3FBRlDMWUAq8GTUlDzxGk/Oo5nyOKNdh4La/yYYVSBuq5GhzJaffX/hXA8jOgkTSBgVKXeSj1CspH4j6q8rko9/RazkNRsiQSzW2hlLu6QBDs8goPK5QyUM+VSXE0egJUsxBJZo36NBRWi6B8UNpwNUQ3sTHkcTVVKBaXV0AMpa4Wivxs5CTP0QgrlDJws7wattqwksN8TMrji2i0tVJVbdiyvKop23cZj4WidXlFcI2Ilq+3oWSNRGZDGWrQZDBQAs9vRAwlmhTq6ssYLbBCKQP1XJkUR61dXmFIm5zbdSQSl5elLL8iqPyCvrRh7f4oRrhVFuBaHuEOjEJm44rayuVlyf5ouIUSgbVWWLG48e06LjS+94kT8rkyurxqPbExhOsj5Rk6NirjK5oXPCnXiSnbSP+90usAepdlNBlr1Voo4VAZabWc2Gh3eQWVWbVYoYnCWiuUwQolLKxQysBmodQrbTguI6Zq5CyMMIP3A9W5v2ydRTO4TIRFubrHRfh8mGS2KpQAGeproQS7p8OV4Xw2geMhNrBCKQP1kNbaQhFCYEfvYGXnNvDFWq4MUbq8jPvLc3ld/YdlePqNXaXlWCydaILydbJQIhxR2y0U/XmBM+WrFaoMogjKq/vGFkp4WKGUQcFC0SfpZ3PRPHi/fGId5n5rEdbu3KcVIOgqB8qzb/sd3v22EWQ+L3Dr4g340C+eK93nmylfWlA0QfnqKKyca7lOhCNqU8Zi2LThWsWjwqLkj8IdeqC0qXrACqUM1INlUhxRLQ75xBpnJL2+e6Dsc5vh2bdZF+HKCG7M/vkjwVcazDrzhnR9oC0o38jX6ZYrQ5QjalMJhSwvSwxFa+2Fu/Y37l2BS294JtzBBiKxUFQZ7PMKDSuUMlAdl/md8tE8eGnZ8+WKFJco+VJKU4ymohjZudaYIShfejkj+4cdhdKSSpbs8ymmZp2HUuZx0UzmC3Z5GWfKB0gbti5/9eQ6PLu2urdR5CNQrmEC+/uHczj1GwvxyGs7zAeNIlihlIF6NqtJG968dz/W7Qp+Na9yJxSXFzRpzHOU+63RKw9HEiw3FFJOP7FfrmzQki593L1KRL84ZPjrmKhaKbkur+ByouhE3UsailBLstjThjUF1HGwE4X7L4zFt767H7v7h/HtB1aWXf7gSA4X/fhJvLxRv3JzHGGFUgauyytEUN7U+M+47mGc9b1HA6+jUn8ryVBpBgulMFKuJiDqL8t8FftvHso6vWAqUfq4+6u4WS2UsC4v57OWMtvKDsrOK1eu6p4f4fusrAy7HGrQVsllVmztxdJNPfjaguWViNeUsEKpgDBpw9W0aTXbPVsU/A/zgHv31HoeysIV27XrmkWZ5WXCH/sIPnokpxRKaYX4BwEaOaL4LdUqFKtydYh2drjB1RjSCtIrlPJkCLL4czLRYsSwimUhhlLeNb1EWZ86Cq8baIJRYESwQimDcpZeqWaUmJRPmin4H1SybSmRqHh541588uYl+Pq9KwJkqbx8m/IsJ8tL1WNSp1A8/ZHW5RVBb1KvNzZGYRkWl1WM6uRNqxEH1Ve59WBSFgBw1wsbcfUfluEXT6wNvFZ1MRR7GW7ySAXlq3fKhG3PcYAVShmoW2saOfkVSuXXUf7p4ge5mR6t3v0jAIANu0sz0QqNrAp3g6WhljMPRXVMGo+XNVvMvY8B1t7+4Rw+9ItnsWZHn0HW6gibvupaDzV8h4tSzqb9KtNRV9fl9o3Lt/Qa9+0dcJ6/Pf3D2v1RzJQPE4dR9VFJx59wl8c3HxO3BDNWKGVQCMob3thomdMQmgoW3nP3+b5XOzI2n59wR1elxwT50cNfO7gMn4Vi6UBf2dQDAEhqfIDecgL9/gG/5Zm1u/D0G7vxzfv0gdmqLZSwx0UYQzFdVA2mTJZIUNyv3E738dd3Wo+xPR+1jqEEWVE2TANHL3GbVMkKpQxU52m0UDwuKuuSHgENz5SdFebZivL5C+oc3Jc5BbSnWsZSvPVrm/9zjQx66ibj2WIotjXFAHtgtuqZ8q5OszxTdYihqLie6RrqmakmhjJtXBsAYFx7xniMLT4YxZycMC6vEctk5rte2ISNlvlkQfVS6+WcooYVShnY0ob9FkpwWWEmQRYfEcaVFGUMJeg3uAFFnYUSsC/8tYPP9e4uTl4wYQvKB70PJfB+KteFYXf1ijWsy0t9Vt8JmX6v1eWVN+8P66ZMulmOlY/+XWutCvefkjDo3hcSPkq70r7BEfz7nUsx/9eLteeGeSMkWygHMOrWhnnBlu1BCDP6r2QOhjB8r4TA3xDidbO1DMp7Cdtp6OZOWJdtCeFGsmXr+DLSKqiUsKeEkfXRVTuwZe/+ENfUlzGcUxZKsEKxWShhnv+gLC/XKjTsr1dQftiNz5U+W+tlfNEU5wnjGo6ZgcIKpRxsFoq3AdiygwIblOX6gdTL5RUwLA+b5hqENShfgYWiy/Kyv2Cr0EGaOlmyTfIrw3LVnl/0aTsy6Bof+/XzOO8Hj1tLMq42HGCBAJ4Yi/a5CFcPqjarybBz6yySQY35GGWx6azfzVJxj21La88N405ll9cBjXNzc3mh7Vy8bizbgxDO5aU/JrSFUqW5HJguGdAY3G0hLj8wnMW8by3CU2v8KwG7lzZZaZ40etOKAAAgAElEQVQdYS0Ubdqw8H7X/Bah/14OtmvYCHsf1XVsz96+oWyYqxquodqASQalgEvPt7kXFUpBB1oolgypKCyUMHGYkQALRc3RyqT03WyQNafwKtU4pBCzQikD7/3UNVrvzbcpjOJ1urzYGksQvhF3+af7COqoRwJGooXr2yVYvX0ftvUO4roHXis6N7gMn/skZEXpFUpwfYVxYxZcXvrripAdqYnwqw07B1STeVR8TbMswRaKbq8vMy/EYKWq0bmwX8dGmNRjVddJjZGqVmgwLVMTxkXpfbZjoE9YoZSD937qOjFvO7aNmivxD4cZoUT5zAU+6PIHamUK2QEC5mydsJ2nVxbbcbq0YfvikPrvXtwBQBULWQoh8IvH12JbT+l7cMLeUyVfFArFNvI3JqYEBO3DxlByZSgU84DDPvoPW3aYLC9dUH5YLfmj0zZeGQNkCGvVNQusUMrA34npFEqhIdv8+mEejuJr2EaHxVT7/AWN/N0JXbrrFn2GobRjCO4QfMo9oKq9VahzS/hHzbrzw1gowWnDttRkwPG3X3v/SvzTzc9rZAxXo+o6plTWclwm5qVX1H79ebmAjtxflwHXlvczOHAf7p32VVkoUo4w81B0EzkL1oteVvf3hagLILrXY9QSVihl4L2dOgvDlzZcjYXiZlD5jwnjQorymQu2UMwur8LLjezCmDrjQmPWnxc2KO/9Dba0YV3LDpOGbXNRhnH1KHfc9t4hfSFlMGzQsOV4kGzL11eSNuyt3jBLtIRZvdusxP2fleD+1iDXb0DasLJQTMovTJwnzGCkmWCFUga+GIpmFJgrYzQROLGRVHlFCsU+oClSOtU9gWFiKLqnvDA6tF/D1Bnb5tyIkMrbW4e2oLx+VO39HtyJ2twvQWWo6wwOmxfbtHUoBQvFsJJDhOuS2eoi6FwguH2oMirJhCy+VhRuosC04aw5KK/26Z47oPDchp3YGIeML1YoZSCEKEy60sZQ8p7vwTc/zOjL7PIKklH/vZhXNu3FC+uDX2IU3OhlDCVAhiiWvTBVk8/lFRjgLXy3BeUrdXmp64dZ4cA8YdCpzxGNdrS9vbL4OiaFUk7nalR8+eD9QRZK2KWJwigUG6o+q0k9DjOpVbVjXZhE7TMpP9c9GDhRuVSeZqZihUJEhxLRI0S0koiWE9Hn5PYJRLSQiFbLz/Gec64mojVEtIqIzvdsP5WIlsl9PyJpIxJRCxHdLrc/R0QzPOfMl9dYTUTzK/0d5SAQ/K4SX1DecvODG4tpjZ8QLgDrEQ7v/fFT+MD1wa9ZDWqMhUXxdDLIDqGMBmD6pWFWGw4KyttcXn6XVum1wgTlbZ1WmBnibnaUrj5DVqP6rWpkXEw5HbQ9KK/fHzix0TdPy3xt9dyEcQubKMzoDz4uiDAxyxHXCintSpUr1tQOwrjUynGjNwPVWChZAJ8XQhwHYB6AK4loNoCrACwSQhwFYJH8H3LfpQCOB3ABgJ8SkXon6/UArgBwlPy7QG6/HMAeIcSRAH4A4DuyrAkArgFwGoC5AK7xKq6aIYCM+64SXecT3jwNM1O4uNEWXF5BJor2a0WEi6GUHuNaFyEagGnegvrXVE/+BAlz+ba0S9vy9TaFEyRj2GsA5lcVOOfYR7KOfM6nKShfnoVikcWiGG1ZXoHPVk4praCBAgLlGAmQQ/H5O5Zi3rcWGfeHmYdy+5KNAPQKLmuxtMJYT2Es5J19Q/j5Y280xTyVihWKEGKrEOJF+b0PwEoA0wBcBOAmedhNAC6W3y8CcJsQYkgIsQ7AGgBziehgAF1CiGeEUyM3F52jyroLwDnSejkfwEIhRLcQYg+AhSgooZqRF8JNAdTFULxKphqF4r2e7v9Al1eEicOB+fchZkRXkslWWoZdtsAVbi1Woz1t2NugTbLYLBT7sdkQLkQb6jiThVLOCNe2jIyxkwx4RsPGAwppw2b5bJ2x6/IKOOzuFzdhW29pmrbC5nYFgE179stjSw8qKMbggUjgWl4h5kF9/s6l+PYDrwUu918vIomhSFfUyQCeAzBVCLEVcJQOgCnysGkANnpO2yS3TZPfi7f7zhFCZAH0AJgYUJZOtiuIaAkRLdm5074cdhB5AaSlhaIzY/MRKRQyHBPmAQ8bQwlDYKPPmeehqC1hrm9baNOYvopwdW2PkRS+661OfVk+WS0ddZhMnSALJWx9qjoxBuXLcUEaDrUp+iDXXdiMpcLI3lyxQRMogereU6IoJ7Cvqw+rhaKUc5C1FiKRYd+g826YoWxpQke9qVqhEFEngLsB/KsQIkhF6ryeImB7pef4NwpxgxBijhBizuTJkwPEs5MXoqBQNA+7bwRmeQjDrD9lCsoHzjL2fa9OowQ1xiD/tgjRUBS20Zvxt3otlKARnsVlYJtbVFZQ3jhJM9gKAjz1qTVR9OUWo4qIIsvLlo1mVq7mxSPLzfIKesZsnfxIgBxhqXZgpNq4cTFZFSsKcnd6bqVJDpWybFtKvx5UpVCIKA1HmfxOCPEHuXm7dGNBfu6Q2zcBONRz+nQAW+T26ZrtvnOIKAVgLIDugLJqigCQTprXGSpnccgwSsE0D6Ve6YNhGn3QSDSUy8vgIrEFLP0uL7NytsZQLBZKGIvPGpT3WTn6Y7IBZk7Y18wWll4J7sDCYDoybNqwbq8vrT7A1RkmyysoKcS7P0xTsbnvguqtPZM0HuMOiiyDpiGDi7K4XJOcyYBEoXpTTZYXAfgVgJVCiP/27FoAYL78Ph/APZ7tl8rMrZlwgu+LpVusj4jmyTIvKzpHlfVBAA/LOMtDAM4jovEyGH+e3FZThBBIqaC8ptGWE5QPdHEIQ4NSo8NA6yC4Ay2HIBmDgq/Vjuy8201Kze/yMpfvPd1mgeg69VAWinXwoC9PV0bQGzBtqOsYYyhlBeUt96UCy9IWr3KOKS1Lex1LZz9iWWZfd2yJLJYUaSGEuwCk7hDVfmwWSpCrKszzp+K6USy5Uy2pKs49A8BHASwjopflti8DuA7AHUR0OYANAC4BACHEciK6A8AKOBliVwohVE1+BsBvALQBeED+AY7CuoWI1sCxTC6VZXUT0TcAqHUqvi6ECJ5UEQHeGErQiATQd4Q+90qI0X/xMW7KZliXV5UKJagxBi69ohRfGS6v4o5UtY1qOx5bg/QpFE053rJNl7G59nwrI5s6wKBBgvq0VKfabZInirRhtd20P8hyDeMSzoacyxXkWgMQ2NEXM5LLozWdNO43iTEwnAtMUrApX9tSOcXnmgxxcl/H3XgqVihCiCdhnrNzjuGcawFcq9m+BMAJmu2DkApJs+9GADeGlTcK8kIgE+Dyst38sFkurp+6OCgf4txqlYiwdLAKt0FrjwluSF5MWUG2FWfDWoPefbq6KceCMVmGdvem/np+OfPG/eHThi11FkmWV7BlEJS5tHbnPqss3u3Bz588xtAZD8gVB8IMakxlWFOkLa/8HrEoFK9BMTCcRXumtDsOY92qTjjIo1AveKZ8SIQQEAKuy8s0IgmcSR+yEzTN0C1n5A9UFpT3Fr1qW5/xONXYdetG2dbh8mKcIGeoAxdReqz2MEuDtFko3kZqVAZyR5igvNlFY7f4bPh+q+75LMMla9qdtyitwuzvUv70ciHMaUufDjrGuX4+UI79I+EVim51Au+55oy2wnl6iyxcUB4A+of0bq8wLi/13DWDy4sVSkjUvVSzrXWjgZwQbtBeFyj2rRwaJkuqJMsruDFH8y7xQhlX/2GZ8TjVSHT1EGbZb4VpCRdblpffQgkX1LRNttPesxANupygvN1VpKvP4i96bFlU3v22zse2QoE5uSDcvTfXZWlZOlwLxSCIcnmFscpM17EN4GwuvKxlgqb3udmvWcOt+BhrZhsrlPigbqV6+5q2c8oLdya97t6GHSGaOlP1X14YOh6DRVMOYd3sqrHoHuLyYijqHL3ytKWvesuwHaf7bd7GrvNlZ0M0aPtoP4yFEjKzIACr+86zTZ/R5pVTfw1bMDzM5FuTfI5c4WIoQYOrfF6497I6l5f/s1TW4PtqTWDwWijD+rdo5kLcE2UYxz5teDShHph0wNIr2bxAJuUE97QuMc8Nr6ixWDpHtenEaWONZdsI605THbHW5eVaWCHKsXTSpjL8qw3blXPxOQpvRlRQYBUI6mSNl5fXLXw3luEmJ2jOL/o0YVNc3rrWZ7SVXtMkp9llpfaHl1VXvlNW0MTGvPEY9UymkxRqLpTJ5eW+5jqErLrfa4uheDcPGCwU2yABKATl2UKJEQWFEuDS8gTtrQ06hMuruAxbIFo9+LaF84IIa9WEcXmFccGZXF62lxt5f35wXQa7JbxzAGzrs5lkGbC8oz1MGcFp5IHFa69jU466gYBN+XqPqWY5EUdW/Xb/zHDz+Up8nRxDI87O1lQSQthlMXXEQanxxedpY1a2GEoIl5dfaRkUikaeRsEKJSRuDCVgHkouL5BOBQftFcEBR32jtQWY1RY3jbACCzishRLk8nLny5Th8irWKNalVyyKQuFzeWlE9XasukGCP5NHf43vL3wdQJXL17sdcem+sJ20LcvQe7ptHpXpUkFyemWw3Xujy8tnxQdMWA3orIdyTufc6k46DBTFqMwLExNNMgS350IMxW7hLF63W3uMP8tQL0chKF9Bg48YVigh2Ng9gL+u3A4ASCcCLBBPDMWaURRiddmS1YY9vW7wO+krJ7RCURZKXpR0dOq/MPEY9+2ORdttI2Hv1lzAyMzW6NVoFtA3SJsbKYwVFmYeSpB7J+z8EVuQ2OdOquC3AoX7ZXNZFT+fn/iNM2XsoK7WUOebZCwcV3q8Qrkxg2axezFbKMFzXfwxFJ2M8hkOkVhywxNrrcfEIShfzcTGUcO7/+cJ9Em3RioghpLLe9f60oxYLA1eYQpI+wPRus7N+XTfgmh6g2Bg/Kbw/biDu4zHFQezM6mCGrMF1H3lVBjc9Wc0mctXxyUTpC1rOJdHgpzfXYlVOehRSLq39jnnFb4bR/4BP2LE4BYsxpb04d2vixvYEhgA+4KHanvxbP2HX3NWYGpJywGXJV086BrOPmmhaOpNXbvVE88MmLdoHJwFTdL07k8mSNumRjzn5/Oi5PnwPk/nzT7IcA398V7cGErAEi71gi2UEPR5fOSFGIq+wWYCXF5hVyMeMZjKPteJTqHILidhcXmFjTkY+seSMopHRm4MKMTI2taYzQrH8z1EPCqZIG2DHM7m0SZ7G9tkVZ0o3uycMC8Dq8ZCsennvKXzyQfcs+JzbIMRWydrWv5FdfI2CyWZCA6oKz0SFJRXLi9bvVUaQ1H7M8mEwWMRHJ9Tlz1ySif2GeJwYSwUtbkZLBRWKGWiVvY0urxSwRaM7rvpuNIGFWzhuBaKsWRI2czBRF+wO2gtL8++4uPUf2W5vIp+j80XH3pGv0qmMHRQw9k8MqkEUgnS+uxti0t6g6lh5sxUouTDzoC2PV/eTVv3lr4HxJuFaJPTmJ0nt+uC/kAh7d649Equ0EmHsVB0bsqChRJ8reJresnnRYh5KIVssqAYiqkMta2jJRVyHopBflUXIV2jtYQVSpm4L9gqjm8IgfW7BzwTH4M7p6DG4vpuAzr7oFF5wrK2T5DrzT9zPChts3BccedRzgu21KgqQXp3gEmEMKN+7750St9BKYWSTJBxpnwqIG42ZEk7dmQNHgio65gIevmWF//E2uABje7FUl43mDHLS243Wpa5YAslaxhAKFQdZ1IJSwzFfB11T9pCxlB0z7m/Lgzn5Qqy6qrDazFoMwjlts6WpDuzv+QYz2m2RAZ2ecUQd3HIopt7/7JtAICn33CyNXQjJ3+OvX1EWuryCu6YvO6C4uN9xwXMh7EtRVI4z9tY9C6vMPF9t56KzCp16TCj/jArN2eS+g5qOFewUEyB6lRAKri3Q6tmqf2gt30qud7c1a+/gHt9Ebh4qXfbvsFSF0uYjDabS8u1UAz71cq6Ju+MqoeWlN6NVCyHzs1TEpS3uJ9tbdXmnjO5vHz31LCyBgB0ZMwWin9io14OVQfs8oohKcNaXTv7/CM+3c31p3UGjUj1I3zbaMUbgA7CO/oKitOEWb4eAEayRWWoxSGrsVAswV9fXQRaKM5nOpnQ+9uzebSkko7P3pJooRPFu/S4eQXd4M4L8Fu0xc+O+n/Z5h7tud7rBMXwvNv65Fv+/DIU9ut+ivddJcOGJde92X9Bc0RsQf1MKmEZ0NgVSmvAJGOvVa17zm0ZXN5j0gbl50uB1ip457OjJWW0UPwTeIPluOmZ9Vi/O3jQUWtYoZRJyvAym+LHxRZDCdNYgq6he7jUNpVNYrpCUHqpahhEwSMeb2MpdXn5ywpCdaQlMRTXbaY/r9yZ8i0GF8pQNo9MMmFUOF6ForP4wri8bMHw4nOLn42wacM2C8V7r3UvdbItzGibGKmOUeMZnZWigs+2iZMZg4uyIKtSKOZ7EjQPxadQdC6vbHBdeGVNJxPa9pgNsOIBmflFQGs6aZwpH8ZS8g5G1DvuGwUrlDIxNdjiex3k29Wd70V1OsVtxeryUhaKJSrvGzkVXcQ/ord31ECwy8s2T2PY8GY9e1C+8L17YNhYvs8nrxsxZ3NuDMWUuZd2XV4a+T331Pp6X5gVirdjHCha1ylssDUngteS88qqVyheC6X0moOec3TPt7Jg1DLs3mMmj2kBAMw//XBXVh3qWXJclOYBjapH3e9QykJl71Uy/8hrMRhjFzaXV04Evk0xJ5z97Zmku5hlMb64qdHlZVf09YIVSgi8LqTOVqexBMUeUgkKNMUB26Qt2ZkGuKPCuLzCZBSVWChqRJ9MBFso3nkoWXNdhE3ZLFZKtgl0avuEjgw27B4wlq9Od4K8+nsSFEPJ5URgZp/q0Ma3pwPWJQvuvIqP6e4fLtoX3hoLCsp776euI/bu111FKbpkgrTnq0uqYLiasQ4AXa0pXHjCQXjvW6f5ji0pQxbbkk4GZgmqetQ9o0Oycw5SKMMBLkagkL03vj1t7OzdLC+DNTWSyxcyzQztNUGEtnQSA8NZrRIPs9pwNp/HYRPand/V4MA8K5QQeEMS6s1u3gfk+39ZhW/etxKAEwg0je59y3yEsVCKg7MWl4R6+AoK0N55FV/DdRGlk4FKbyibd69TPEnOe5YtjuIqlKJrqfNsM+WnjGlBz/7SeIBbjsctYbonLakEkkl9lpffQtFbOADQ1Za2BpqBAAvFc0z3voJC8cYtlLwm/PEem0Ip7SRtaa6qkx3XltZ2XEpOFQwvXngznUy4bcmkGNUz3mJwQRZkNQeiXQslY46hDHmUhG6goSyU8R0Z3+RVv6xOuW1p/bOVzQttf6FwXF6EtoyjPPULrXoHE1oxkM0Jt86D3k9fD1ihhMD7sJw+ayIA/42+/tE33O+3XTEPqSRpG5y/gZlvvBqhFnfG3tGt6X0sgD0o70v5LZJTNeK2THCDHhjOYXx72pGrqIywKb2AeRJnYVl7/Xmq3PHtmUCF4v6edFLb6FXmVCqh7xRyeXhiKKXlq/prSydDvdHRGEPx3M/dHgulWKYhQ+eWzwtk8yKwE3WD1emEtuPpG/RO0iy9hvLzj20PVijKMvAeM5JzsuWSASnY3jJsMRT13AS1M3fCqiG7r7gsL+q3TmjPYDCb01oP3t+rS9nN5goKxfRsJRPkyqnL9PIvlml2vXW2lLoZGwErlBCo+/juEw/CVLkWkbdf8HbgMyd1IGNwF33hrqXud9NAc2A4i17ZsIs7KO9ISp8W6nxaZ8oH+FzVg9+aSmIkV7pOl2LfUBZj26RCKWmQ5bu8TC5Es8vL+RzXng5UKEq21rTTQXl/z0gujz0DI3hi9S7zxMZ83l1uJ8jl1Z5JGhVwmCyvkXxhxr7X5VXcGQ4asquUlehmNmndPM62zha9QvjIr55zv+v68gGPhZLNi5LnU11TKTX/wpvOfB71bNqyvFqsCkW5SkvlUL9NPZ+6OvMqZt19c62x9gyEwXpQ96Ytk9TuH8k71i9gSF0WTlBeWRe6TC/veaa2OJLLo0MqFJ3lWU9YoZTB8YeMhXSn+zoXr0JpTSeRSpb64/N5gb0DhY7PZKF86/6VhWOKHqBhS5Cv1OWlx6vsSiyUnL9TMDXqgaEsxrVnnPLyeqUEhLFQ9K4L2+KQSlONa0+jd3DEaB2o83UjRW9mTdIUQ8l7X0lQWr5ynYxvz4Ragjwoy2vSmAyIgC17C5k6xZ2dyUJRiqqwVpbZQulqTVldI7qlV7ydLFDaySorS+fyykrFnCjDQsmLINeYp06L6khdd5y0oHX3JbSF0iGVkqbe1b1sS6cwnMv7k2byzkz7lgALJS+c9b1UW9NlenllMw1CB0dyrvJkCyUGnHn0ZADAp888AknNCMt7E9MyBbW44/COPjoySWODXrvTccFM7MiUdHDZXGHEE2pio+H3eBtHqcurYKF4//eSzwv0D+cwrk3fYMMuMQMEWCjyf5M+Ug1tvBxBDhgCp6p8nS9bKYN/OftIpJLmeShBFspDy51VqCd0ZIzrMfkUimkGec7J0Dp8Qjte29bn2w44SgAwj0AL7izzQEDVRUdLyhdDUMybNQGA05nrXV7O71P3vVi5uRZKutSfn80LpBPkxlBMnaMqw7ZES7E7zbdPLvip3EBahWJJkPHGUABoA/NKKY1pTUEI/cChxTIvKEnBLq8wmaH7R3Ku8uQYSgz41vtPxGNfeBeSCccHXDxHo7jTzSQTJameahHBb1x8AvqHc7j5mfUl19k/nMPTb+xGKkGYO3NCyUOczQu3MfdrOi83y8vi8vKWWzLKVEFROdLVZsDI8w+f2AEAWFc0g9v34iGry8s5oGTuhQrKW9JL1cjM9JIr1+KSjdYrm1Ksh0/sQDJRes/U9VNu1lzp/mfWOisjTOh0FIrumGy+MDfDlIqdzeeRTiYwbXyb796q0ffJh40HYO4wCqNls8tLKbPOlpTWRTOuLYOjp3YiZVg9V933sW7nVfx8qvhbqT/fSaFNuGnNtvk4QR2xVxYAJcpxSGbutQW4krydt97l5dyDCe0BCiWrFLR6tkoHUjaF4rVQdHJ669iUiDM4knfbgcmCrResUEIwbVyb23kSETozKfQP+S0OALj6wmMBOOt9FY9E1fHqWKDUnH/fT58C4LibWlKlgdPhbN7N598zUBo3yBeN7kw+fZ9CMbi8WgMCmqrDmzm5A8kElYzM/TOuK7NQVBm23Hs1MjNZB6oOWtOlDVs14NZ0EhnNPVPHtwTUhWJcWwZ5oe8URnKF+IgpS0vNWWhLp9Dv6eyUvJ0WH7k34A7oM4KGc3kQOR2gruPZP5JDWzqJBJEhbdi5tooj9hYt3+JmeWmC8o7CJLcuTb/Duz6Wt8xSWbKeGElpO8kkE4WRv+ae9A0V2k+Qy6tgoZizydx5N5rZ90FB+eFsHq3pRCGGorFQvNcNiuF1tKSQTBCGc/p6rResUCqgszWFfZ4HUnUAn3zHLADQurz2yMl34+WIB0BJMFm5Oj595hFoTZdOdhrK5jBrcicyqQSWbdpbIpe6ZGu6dMTkxdthlVoo9oawo28IADC5M4OOTLKkIXhH+nYLpRBc9b7DezjndEBCAP/xx2U48sv3a88b2+bUp1fB+48rtlAKAg26CiWB9kyqZEJhLi8wkhMYIzvzoImeKrX4qTWlb97bN5R1O6biOTvq927asx8pOcnNK4dac0spTlMKq+o0VXBWH5R3rKCWdFLboe8fzqEtkwRB33m9vt15Pg8e6yiUnv36+TLtcsQ+6EvNdRSmGrGbRtIqgN7ZYo4J5OSofKLBHeVYKEn3GdZ11N61zALThkNYKLoMq5GAgYxPzmTClVMXQ/FZKJoqU3K2peUglC2U+NGWTuKJ1bsA+B9GFXDsbEmht2itJDW3YEJHBj/8h7c62zwzvD91yxLf8cUWylV3v4Jd+4Yxrj2NWZM63FiLQgiBe17eDKDwEOtG3Gt37sO/31nINitusEoZtAZYOdvlSrVTu1rRplEoWd9ILfgB9yqCwvpMzueUMU7H9bvnNpR05tmc40ZSsQWjhaJ82a6CLMjjbYydLSnfe2+AQmNWHUbQRM9tPU6dfPJm/30UQuCh5dswti2NBOnr84FXt2HV9j4s3dSDg8e1Ysve/a7cO6XyVhPXTCP7ry1YDqDgAtS+8Ckr0JJMoCWpTxtWFoouqQQAfv3UmwCAQ8a1AQBuW7zRf36RBaPSkIVwUppTCZKZXmbXXe9+55yDulpcmYpZJQdeXYYYXv9QFm2ZgstLpwyUddVqmEOyfziHllTBejAplJRHSY5oLJSWVOlzpxjK5tCSSmJMi9mN7d2my1ZTg4+2dBKZVMI3QNzRN4gf/vV1bbm1ghVKBazd1Y+tPYNYs6MPX7jrFQDAebOnuvsPk4FVb2f64HJnNeIJHRlMkCMrlR46MJx1g7uAs45WsYVy2/NO4929b1g792LJ+j34qZwPox5iXQf48ka/ZVOsUFQsQjVWXcfy6uZeAMBBY1vRli5dets7Git2i3jJ54XbOQAFC0A1CuXeU/hSfmXWkBqRF1sXisGiUWRWY6G0SIVS3PDUaG9Mqz7H/2ePFeYf/fv5xwAoZDgpHlm1A4MjeSzf0ouOlpSrILw8sXqn+31iRwYjOeHKrazBWZM7ARQ63GLUKtdqRK1bYPLGp9ahbyiLlnRCO/IfGM6iPZOS1prZdaKe3ztf2OTbrixfZcHs6ndkdxV3JgUiQkvKvFx77+AI0klyLTqdHNcseBVAIRuuuLNfubUXR08Z495z3TPYN5hFOknobElp28kbO/dh+vg2d3BW7FYDgDd392NqV6s7T8lbTiEZxOyCHMrm0ZJOuKtv6AZFdywp1HGf5ne4Vnam1EL51ZPr8MO/rsb9y7aWXrxGsEKpgv6hHP74kmMVvOctB7vb58wYj77BLLbIlxht2D2Au2A3umYAABsPSURBVGTjm9hZqlCK4yFJItdCKY5BXHzyIRjTmip5uHo8ZbgdoKahFC/rUdwRK6tJvftbNyL+68rtOGFaFw7qakWbpvMZyeVdGYrdIl6++9AqX8eiRprKspra5VcovqyhnJM1pBSKrjGu2taH6x54DUChTryjUeU+ak07iqnYbaaupxp8sQtRlX3m0ZPRmk7iwhMOwjQ5eld4Ff9bpo/FSo8C9f4Whds5uQrFeYaOP8R5HbNOIXk5dIJz/eL68CrulpQ+y3BwJI+2TNKxOkfMAwFlgRSjBiNTu1oxtasFK7c611TPqroHh09sd91nxfQNjqCrNe3GJXSWgYpnfu29xwPwWzHqvURHTOlEazqJCR0ZbN5bumBi3+AIxrSmkUro54xt3juImZM6XHeUTo49/SM4ZFyrG+/RrUSgfofOQhkccawgpfh0K0ArkgnS7n/8dcdTorNQdu9T/Yu5DUYNK5QKOPmwcQCAvZ7OwhsbUQ/8+m7HLfWVe15197VnUq5C2SM79x8sfN1XfjLhDV46SiWZIHzyHTNx9rFTMaY1XdJh7PB0NMrfrouhbOj2r3v1pbuXud/7h7L4jz86sh45xRkR79rnfxi//5dVWLa5BydNHwciQls6UdLYRnICkzsdZWCadPjq5h53hK86YfXgq4apXF4Kr2tjOJtH2tMYdTGUz932kvtdN4pUFkl7JiXjYlnXVdSzfwR/XroFADCmNV1yrk6u9kyplaOuCziTXtft3GdcWXnGxPaSzulb9ztK6+CxrWhJJbBzX6lC8Sr9I6Ql01cUI9i0p3DfncGK3n3Slk6iPZMsqU8hJ+H9n3cdgc6WFE6cNhZnHTPZd4wajIxrT+PgsW3YK///3kOrABQUygnTxrou42J692fR1ZZ2Y17Fg5W1O/e5gzPlBvQ+F0PZPIZzhayn8e1p32BL8eqWXhzU1Wp07+3oHcSUrlajQsnm8m6Gn7rHw574mJK7y7D2H+A8Y2Pb0sikEmhJJUpcruoZ+NSZsxyXrMZC+fq9K9zf6QwUCnKqONFeze+vFaxQKuC/PvgWAMBSj/tIKQnAGYEBwHq5aOGqbY6L6MITDvId+/LGvXh541489Oo2X/lEnuBlNo+B4RxyeYGJspMe01oao/G+gU9NPNvTXzoyeVPK9LfHTSn9XbLhA8C08U4n7x0RCyHwvw+vAQBMH+/8xmWbe/DkGn/n0D+Udc83Pcx/979Put+/8p7jABSsJzXiPGisX6Hs7i/Isnf/CMa1pd2UTZ2fOOVZdnl8UZ28sH4PbnrmTQCO1dgpy1HzWf5zwXJcKyeZjrHEUJas7wbgpFxv6RnEGzv3AQB+++x6/PPvHaV26yfnYeakTvQOZkusRKX4fzn/bW7nVGxBEBEOGdeGB1/d5lNIi9d145ivPOj+r+pMjWaFELjoJ0/h0799EQCw6PNnojXtWCi6uVLtmSQmdbb4BihqX14UXKHtmaQvuQMAtvc65xzU1YqDulpdq0i5xpTV+4B0wXhdfYrewRF0taaMqbTe57xVk8WlnjelUMa0pkvayp7+YSzduBcnHTrOSaAp6uxf3dyD3f3DmDqm1U3TL36OF8jBxvNv7kEm5TxnXuvAXVVAPnc6i1ApFEfOVMlLz9SzOn18O7ra9ApFceiEdmSSBVfm9t5B182+cc9+7NC8obMWsEKpgCMmd2JSZwb/7bEslBIBgKljHDN4o7QGcnmBtx8xET/98CkAnIbQlk7ituc34uKfPIUPnDodAHDB8Y7CSVIhW2dbz2CJy6BLury8o57tMiicSSZwwiFjAQDX3r8Sy7cUXsr02VtfwuOv78TFbz0Ev5z/NsydOQHHHdzl7u/1WBPKwvAqFO935YJRnaEa2W/eux+7+4ddq0P3qlkvF731EFc5qY5WxYumj/e7j9btKoyy9w4MY1x7xnUp3FvkJ35o+TY31nPCtC5MGuM07FueXY+RXB4fuP5pvLKpB5lUAmNaUm5WkVJM3i5Gzbvwr3VVmsn2lumO5Xr/K1vxp5c24yt/Klimb5sxHrNlXS9Zv8fdvmvfEP68dAvmzpyAI6d0lszTGNOawsfePgOA48ba0D3gulkB4DO/fcH9/ovL5qAl5bg+1Gj3zd0DWL6l1z3miMmdmDGpHUIU1jHL5QVmXHUfBkfyaE0ncej4Np9FAwDfe8h51pVFqIs57egdRCaVwNi2NObMGI8dfUOulQI4lgkA95krnr90z8ub8eiqnT4LZb/HJTuUzeFDv3CWh5k3a0JJ0H3XviHM+/YiAAXl1dWWLomh/O1/PwbA8TRM6MhgrRwAKJnUYGdqVwvGtWeQIKdsxWvbevGLJ9a5/491lU7htypPgEpg0A14evaPuAqnqzVdorRWbHXu26SODMa2pUtcV/uHc0glCP/0NzMxtasVXW0ptw3duaSQMPHnpVsw91uL8KrlBW1REGuFQkQXENEqIlpDRFfV8br4J5kiDDhpvsotAjjZXh2ZJH7++Frs6R/Grn3DOOPISSDPWwmLR14dmSQmdjoPVzJBOPtYx4K475Ut7mhTXWOsfAiP+PL9eGnDHuTzAtt6B3HcwV149svn+ILZa3Y4jaVvsODCOePISQCA6ePasHJrL7b2OD5m5Z665u9nY1x7Gukk+d4At0Y2vC+/+1i8U64e8KkzZ/l+zxnXPQzAUTSzJnXg4ZU7fL/z1sUbMOOq+wAA7z9lGv7n0pNdeX/33AZs6xl0F9vMJBM4b/ZUnDTd6Yi8r8B9YvUuZJLOe0zGt6exdONePPa6M+IdHMnhU7c4He2YlhTu/ew7MEkqyHtf2eq6YNQ1iMi1dJTSmFIUvxnTmsKi17bjJ4+swUgu7+uk/mHOoQCAj58xA1PGtGDltl786+0v+85PJRM45fBxIAI+dcsLuPHJdXhy9S7M+eZfARQC/gXXnHAHE0qWQTnq/b93FLL0vAtJHjG5w/3NKnivG5keNWUMALiWlHpGAMfCmT6+HX2DWddVtG5XP258yulAp8h71d5SGjt7bl03OjJJEBFmSLevuifnHz/VHST96mNvA+DEILx87janzt7Ysc9VFl7X2z/e8Kz7/RsXneBOOlQDnUdeKzxrykIe05rC0o17MeOq+zAwnMXr2/vcOjvhEMdtt3xLr9vhr9u1z1dGMkGY0JHxuX4vuf4ZrJSd/e1XzMPkzlafHHv6h/FFmazzthnj5e8oPC89AyP4yC+fw+BI3lXQh4xrw6aiWM+yTY4COP6QsThsQrvr8VC8uqUH2bzA6Uc4C9YePXUMVmztlenfKRRz9NQxJduiJrYKhYiSAH4C4EIAswH8IxHNrtf1P33mEe734mwkoBBof+d3HwFQGDEplLUBAL95+k3MmtzpjtxmTup0A5s/engNzv3B4wCckQoAHHtQ4cF430+fxvcXrsJjr+/EtHGtrjvtF5fNAeCMch57facvzfh9JzvvpDhNLrVx0Y+fws6+IWzvG8Q7j56Mj58xE0SEg8a24qZn1rud0pI398jzp7tlKevim/et9LlyvnThMTj3+KlYumkvNnYPYHAkh3xe4EGPe+/qC4/z1d9jr+/ERs/I+KipY3DDZXPwpyvPwLj2NJbKuTeqTlXG2ufPczKs5t+4GAB8S5eoevbGuH7++Fr3+wXSDdkllfWegWGs2dGHnz9WOKZn/wj6BrN4dXMv/uuhVfjYrxe7o9FzZ0/Ft99/IgBnoHHarIlYvK5ggXhpSSVdF8rX713hW4xx/tsPB+Cs8gwAP3/sDTwuXUIzZef8f8460j3+op88hT8v3eI+C9PGtbmd6CHj2nDr4g3oH8q6z+ExU8dgwT+fAaAQY/vjS5shhMBfVxYyDN9/yjQ3sL/glS3YN5TFM28U5taogcTYthQ27RnAUDaHHb2DuPeVLXh54173ejMmOTIrJaHSqp1z05g1uQMPLt+GXF7gL8u3+Z6dHX1DOHhsK9ozSZ819+KGgov5qKlj0JZJYqIMuufyAne/WMiIOlTGV7o87Wxj936cJ9sSAMw+pMuVU13nawsKMYm3zXDax6TOFjejbP3uftf6u/KsI3DarImu9aviW95VMJSL8D9lrAMA7lm62XUTK4UyqTPjc1Fv6xnE9xe+jhkT23HYxHbMnNSBDd0DruXaP5TFJT97BkDB8jvz6MkYHMnj+sfewA//+jrSScInzpgJwHF1qvhcLSlVY/FhLoA1Qoi1AEBEtwG4CMCKwLMi5Lr3n4ir/rAMJxzSVbLvrk+fjg/+7Bn34ZsnRxGKez/7N7j4J0+5DfDE6WPx4dMOwymHjcdsWd77T5nuWxpfPfxnHDkJXzj/GDfm8ZNHnGO8o5JzZ09FgoCfP7YWP39srWvxPPC5d7hrU11y6qH40t3LsKNvCG+71hkpf3Te4W4ZV11wHK78/YuY+61FPtm9ClQpuVsXb8CtizcAAL53yUmYMqYVsw/uwkhO4B1SAXi57Yp5bjnJBOFDpx2G3z+3wZ1PsfDf3ukmBhARzjl2Ku5+cRO29z7tuhPu/PTpAID3vvUQ172krB/F//v72e41ivnuB96CS+Y4ynGWHN2rRurl7UdMxMFjW7FVdopPrdmNM//rUQDAh047zJ1/BDjzd7zukb85chK+K2NuJr7ynuNcJX3iNMdt9oeXNuMP0rV1mnxlwqfPnIXvPOgE6Zdu3IvP3urEZz515ixXOQOFzL3jr3nI3XbnZ053labq5B5avh2f+e2Lbn2ePmsiWlJJTBvndMZf/dOr+KrHbff0VWe7cYu3zZiA3z67wRe/AZw6BZwEAy83fWKu7/8zj56MXz/1Jo4omrAKOIOh1nQSZx07xfdcKe7+zNvd76kk4dbFG3GrZ07Mum+/2/UGHHtQoW2e/8OCMrnnSke5Hj7Bue9qMKJY8pVz3Wdm175hvLatD5+/Y6mbnXb9h09xByPK7frdB1fhuw8WrN9fXDYH6WQCJ04bi2Wbe/C1BcuxcmsvnlvX7R6jnr+2TAobugcw71uL8H/OOgLXyvcrnTbTufezJnUilxe4ZsFyvP2IiXhsVSH+pKzGOYc7CvBHi1bLcyfgq393HL54wTHufas1sbVQAEwD4J1ZtUluqxuXzj0ML371XLfBe5kzYwI+fNphAJzGVZxOevjEDtz3L+8A4Iwu//28Y0BErjIBgC+efwzeeqjTwXx03uGuPxYArjzrSDxz9dk49fDx7jZ1PcUphxX2PfzaDqQS5HacgOOa+8u/vRPHeEzhYzzWz7tPPMh1Nyk+9c5Zvv/ViNWLCvifdewUUGk/jo/OOxzziursM2cegXSSsHxLL46Y3IGZkzp8+y966yEAnCAoAHzt72fjJFk3Xa1pN1HCy+prL8QFJxTSua/5+9luPAJwUr1Vx3Po+HY3fqH4zgdOxJvXvQfTx7fjiS+e5VO2imI3woeLjvnGxSf47puyZhTTxrX53KeTx7SUWLzqvTNEhJe+eq77TCguOfVQ3/8fk6NSL10el+yYlpSbkv3g8m1YsbUX586eit/902nObzqos+R8oNBxAcC7ji5N6rjg+IPw/73NkSWVTODFr56LuTMn4O7PnO7GChRfec9srWW/6psX4Cw5+Pns2UeW7P/t5af5nvmzjvHL8cULjvG5lt9x1KSSMn78oZPdZ+eoqaW/9XuXnOQbgHzh/KMBAHe/uAnLZBziwhMP9l3nC3IekleOc+XctC9d4CzJ9Jun33SVyamHj8fiL5/jurHfeqjTzrb1DuL/3bMcQ9k8/vmsI/FlmbByjmxTv39uA/759y/hzhc24ZTDxmHF18935Rjbnsb5xxfmw/3nRceDiOqmTACAbGstNStEdAmA84UQ/yT//yiAuUKIzxYddwWAKwDgsMMOO3X9+tJFGWvFcDaPP728GW+bMaGkg1SoyWSVkssLLFyxHWcdO9md0Ogt+6an12NMawq/fXY9rvn7411/azF3LtmIiZ0ZvOvoKb4RN+CkSCYThJGcMJrNe/qH0bN/BIdNaC85//k3u3HzM+vRkUnic397FKaOaS05RpXx+8Ub8P5TpuHgsW0l+zftGcA9L2/BwHAWnz37qJKG8tq2XmyQQehPvnOW604oZtHK7Zh9SFfJNdTM55XbejG+PeNTBF7+66HXsGjlDvzkw6e4abpe9g/ncP2jazCmNY1PFilgxbNrd2Pl1l6858SDMaXIHTqUzWH97gGs2taHKWNatAMWwPHFb9673zcIAZyEgaFsHsu39ODZtd34x7mH+bIQFVt79uOc7z+GsW1p3P8v73AnE6rfsL67H9fetxKb9+7Ht993YokcPQMjuGPJRpwwbSyEEHjrYePKepbzeYGXN+1FayqJWxdvwL/+7VFuJqN7jf0jWL29Dy+s34OPzDvcjcMoBkdyWLuz3303u3J1efnriu0Y257Gmh37cOSUTteVpcjlBbb3DqKzNYWtewdx1JTOkuczlxdYu3MfHlq+DUdPHYPzZAJNMY+scuI47zp6stvR5/ICv3hiLbb1DOLEaWMxd+YETB7TUvL8/urJdRBC4IjJnehoSWHuTL+cr27uwVNrdmHz3v2YOakD7zt5WomiBgqvtY4SInpBCDHHelyMFcrpAL4mhDhf/n81AAghvm06Z86cOWLJkiWm3QzDMIyGsAolzi6v5wEcRUQziSgD4FIACxosE8MwzKgltkF5IUSWiP4ZwEMAkgBuFEIsb7BYDMMwo5bYKhQAEELcD6A0TYRhGIapO3F2eTEMwzBNBCsUhmEYJhJYoTAMwzCRwAqFYRiGiQRWKAzDMEwkxHZiYyUQUR+AVQGHjAUQtMbzYQA2BOwPU4Ztf1Rl2GRtBjnjUp8sZ7RlxEVOIB7tCKi9nMcIIezLFQshRs0fgCWW/TdY9u8McQ1bGYH7IywjUNZmkDMu9clyjk45w8g6WuS09Z3qj11efv5s2b/Xsj9MGbb9UZVhk7UZ5IxLfbKc0ZYRFzmBeLQjoD5yWhltLq8lIsR6NLU6v57EQdY4yAiwnFETFzmB+MhaaznDlj/aLJQbGnx+PYmDrHGQEWA5oyYucgLxkbXWcoYqf1RZKAzDMEztGG0WCsMwDFMjRr1CIaIbiWgHEb3q2XYSET1DRMuI6M9E1CW3p4noJrl9pXoHi9z3KBGtIqKX5V/pK+3qJ2eGiH4tty8lond5zjlVbl9DRD8i0r1TsSnkrFl9EtGhRPSIvIfLiehzcvsEIlpIRKvl53jPOVfLOltFROd7ttesPiOWs2nqk4gmyuP3EdGPi8qq9fMZpazNVKfnEtELsu5eIKKzPWXVtE59hEkFO5D/ALwTwCkAXvVsex7AmfL7JwB8Q37/EIDb5Pd2AG8CmCH/fxTAnCaR80oAv5bfpwB4AUBC/r8YwOkACMADAC5sUjlrVp8ADgZwivw+BsDrAGYD+C6Aq+T2qwB8R36fDWApgBYAMwG8ASBZ6/qMWM5mqs8OAH8D4NMAflxUVq2fzyhlbaY6PRnAIfL7CQA216tOvX+j3kIRQjwOoLto8zEAHpffFwL4gDocQAcRpQC0ARgG0NuEcs4GsEietwNOSuEcIjoYQJcQ4hnhPGk3A7i42eSMUh6DjFuFEC/K730AVgKYBuAiADfJw25CoW4ugjOQGBJCrAOwBsDcWtdnVHJGJU9Ucgoh+oUQTwIY9JZTp+czEllrTQVyviSE2CK3LwfQSkQt9ahTL6NeoRh4FcB75fdLABwqv98FoB/AVjizUr8nhPB2nr+Wpu9Xa2pW2uVcCuAiIkoR0UwAp8p90wBs8py/SW5rNjkVNa9PIpoBZ3T3HICpQoitgNOg4VhNgFNHGz2nqXqrW31WKaeiWerTRF2fzyplVTRjnX4AwEtCiCHUuU5Zoej5BIAriegFOObmsNw+F0AOwCFwXAqfJ6JZct+HhRAnAniH/PtoA+W8Ec6DswTADwE8DSALx+Qtph5pfuXKCdShPomoE8DdAP5VCBFkaZrqrS71GYGcQHPVp7EIzbaaPJ8RyAo0YZ0S0fEAvgPgU2qT5rCatXlWKBqEEK8JIc4TQpwK4FY4vmjAiaE8KIQYkS6apyBdNEKIzfKzD8DvUR9Xg1ZOIURWCPFvQoi3CiEuAjAOwGo4nfd0TxHTAWwpLrcJ5Kx5fRJRGk5D/Z0Q4g9y83bpIlDulx1y+yb4LSdVbzWvz4jkbLb6NFGX5zMiWZuuToloOoA/ArhMCKH6rLq2eVYoGlS2BhElAHwFwM/krg0AziaHDgDzALwmXTaT5DlpAH8Hx83TEDmJqF3KByI6F0BWCLFCmsh9RDRPmueXAbin2eSsdX3K3/4rACuFEP/t2bUAwHz5fT4KdbMAwKXSJz0TwFEAFte6PqOSswnrU0s9ns+oZG22OiWicQDuA3C1EOIpdXDd23zUUf64/cEZMW8FMAJHm18O4HNwsipeB3AdChNAOwHcCSfotQLAF0QhE+QFAK/Iff8DmV3TIDlnwFlVeSWAvwI43FPOHDgP/hsAfqzOaSY5a12fcLJ2hCz/Zfn3bgAT4SQJrJafEzzn/Iess1XwZMnUsj6jkrNJ6/NNOMkb++RzMrtOz2cksjZbncIZqPV7jn0ZwJR61Kn3j2fKMwzDMJHALi+GYRgmElihMAzDMJHACoVhGIaJBFYoDMMwTCSwQmEYhmEigRUKwzQJRPRpIrqsjONnkGdVZ4ZpNKlGC8AwjDNRTgjxM/uRDNO8sEJhmIiQi/g9CGcRv5PhTOS8DMBxAP4bzsTYXQA+JoTYSkSPwlm/7AwAC4hoDIB9QojvEdFb4awo0A5nQtonhBB7iOhUOGugDQB4sn6/jmHssMuLYaLlGAA3CCHeAufVBlcC+F8AHxTOWmY3ArjWc/w4IcSZQojvF5VzM4AvyXKWAbhGbv81gH8RQpxeyx/BMJXAFgrDRMtGUVhL6bcAvgznhUcL5ermSThL0yhuLy6AiMbCUTSPyU03AbhTs/0WABdG/xMYpjJYoTBMtBSvZdQHYHmARdFfRtmkKZ9hmgZ2eTFMtBxGREp5/COAZwFMVtuIKC3fWWFECNEDYA8RvUNu+iiAx4QQewH0ENHfyO0fjl58hqkctlAYJlpWAphPRD+HsyLs/wJ4CMCPpMsqBedlYsst5cwH8DMiagewFsDH5faPA7iRiAZkuQzTNPBqwwwTETLL614hxAkNFoVhGgK7vBiGYZhIYAuFYRiGiQS2UBiGYZhIYIXCMAzDRAIrFIZhGCYSWKEwDMMwkcAKhWEYhokEVigMwzBMJPz/g+rkWj3YNfcAAAAASUVORK5CYII=\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": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXl8nNV1979nRqN9lyV5x8Y2izFgsGNMSEoJCSZJG0gKiWkCbl9aEkKadPmkDe3bkqVuQ9OUlKYhIcFhSd4EQhZIwhJDVhww2GDwhrG8YMuWbdnatxnNzHn/eO7II3kkjaQZaTxzvp/PfGbmzHPv3HlsPb/nLPdeUVUMwzAMIxl8Uz0AwzAM4/TBRMMwDMNIGhMNwzAMI2lMNAzDMIykMdEwDMMwksZEwzAMw0gaEw3DMAwjaUw0DMMwjKQx0TAMwzCSJm+qB5Bqpk2bpvPmzZvqYRiGYZxWbN68+biq1o52XNaJxrx589i0adNUD8MwDOO0QkTeTOY4C08ZhmEYSWOiYRiGYSRN0qIhIn4ReUVEfubeV4vIehHZ7Z6r4o69XUQaRGSXiKyKsy8Tka3us7tFRJy9QEQedvaNIjIvrs0a9x27RWRNKn60YRiGMT7G4ml8CtgZ9/4zwLOqugh41r1HRBYDq4HzgKuBr4mI37W5B7gFWOQeVzv7zUCrqi4E7gLudH1VA3cAlwArgDvixckwDMOYXJISDRGZDbwX+Fac+RrgAff6AeDaOPv3VTWoqvuABmCFiMwAylX1efU28XhwSJtYX48CVzovZBWwXlVbVLUVWM9JoTEMwzAmmWQ9ja8Afw9E42z1qtoE4J7rnH0WcDDuuEZnm+VeD7UPaqOqYaAdqBmhr0GIyC0isklENjU3Nyf5kwzDMIyxMqpoiMgfAcdUdXOSfUoCm45gH2+bkwbVe1V1uaour60dtczYMAzDGCfJeBqXAe8Tkf3A94F3iMh3gKMu5IR7PuaObwTmxLWfDRx29tkJ7IPaiEgeUAG0jNCXYWQ82w618/KB1qkehmGklFFFQ1VvV9XZqjoPL8H9S1X9CPA4EKtmWgM85l4/Dqx2FVHz8RLeL7oQVqeIrHT5ipuGtIn1dZ37DgWeBq4SkSqXAL/K2Qwj41n785386892TPUwDCOlTGRG+BeBR0TkZuAAcD2Aqm4XkUeAHUAYuE1VI67NrcD9QBHwpHsA3Ac8JCINeB7GatdXi4h8AXjJHfd5VW2ZwJgNY9I40R3EJ4kirIZx+jIm0VDVXwO/dq9PAFcOc9xaYG0C+yZgSQJ7H050Eny2Dlg3lnEaRibQ0h2itCDrVuoxchz7H20YaSAaVVp7+snz2aILRnZh/6MNIw109oWJRJW+cGT0gw3jNMJEwzDSQEtPCIC+fhMNI7sw0TCMNNDSHQSgrz+KVwhoGNmBiYZhpIGW7v6B18FwdIQjDeP0wkTDMNJAa3do4HWw30TDyB5MNAwjDcRyGoAlw42swkTDMNJAS5ynYclwI5sw0TCMNDBYNCw8ZWQPJhqGkQZazdMwshQTDcNIAy09IfJ83rpTJhpGNmGiYRhpoKU7RH15IQB9VnJrZBEmGoaRBlq6Q8ysdKJhnoaRRZhoGEaK6Y9E6ewLM6OiCDDRMLILEw3DSDGtbo7GzEpPNGxyn5FNmGgYRoqJldsOhKdscp+RRYwqGiJSKCIvisirIrJdRD7n7J8VkUMissU93hPX5nYRaRCRXSKyKs6+TES2us/udtu+4raGfdjZN4rIvLg2a0Rkt3uswTAynFa37tT0cstpGNlHMpswBYF3qGqXiASA50Qktk3rXar6n/EHi8hivO1azwNmAs+IyFluy9d7gFuAF4AngKvxtny9GWhV1YUishq4E/iQiFQDdwDLAQU2i8jjqto6sZ9tGOkj6DyLyuJ8wCb3GdnFqJ6GenS5twH3GGmt52uA76tqUFX3AQ3AChGZAZSr6vPqrRX9IHBtXJsH3OtHgSudF7IKWK+qLU4o1uMJjWFkLCFXYluc78fvE/M0jKwiqZyGiPhFZAtwDO8ivtF99AkReU1E1olIlbPNAg7GNW90tlnu9VD7oDaqGgbagZoR+ho6vltEZJOIbGpubk7mJxlG2ghFPNHIz/NRmOczT8PIKpISDVWNqOpSYDae17AEL9S0AFgKNAFfdodLoi5GsI+3Tfz47lXV5aq6vLa2dsTfYhjpJuZp5Pt9FAb8lgg3sooxVU+pahvwa+BqVT3qxCQKfBNY4Q5rBObENZsNHHb22Qnsg9qISB5QAbSM0JdhZCwx0QjkeaJhJbdGNpFM9VStiFS610XAO4HXXY4ixvuBbe7148BqVxE1H1gEvKiqTUCniKx0+YqbgMfi2sQqo64DfunyHk8DV4lIlQt/XeVshpGx9EdOehoFAZ95GkZWkUz11AzgARHx44nMI6r6MxF5SESW4oWL9gMfBVDV7SLyCLADCAO3ucopgFuB+4EivKqpWBXWfcBDItKA52Gsdn21iMgXgJfccZ9X1ZYJ/F7DSDux7V29nIafoCXCjSxiVNFQ1deAixLYbxyhzVpgbQL7JmBJAnsfcP0wfa0D1o02TsPIFGKJ8II8H4UBS4Qb2YXNCDeMFDOQ04glws3TMLIIEw3DSDGhcBS/T/D7xKqnjKzDRMMwUkx/JEq+3/vTsvCUkW2YaBhGigmFo+TnOdHIs/CUkV2YaBhGiglFTopGQcBvnoaRVZhoGEaKCYYHh6es5NbIJkw0DCPFDApPWSLcyDJMNAwjxQxKhOf56Y8okehIC0MbxumDiYZhpJjBnob3bMlwI1sw0TCMFBOfCC8M+AETDSN7MNEwjBQTCkcJ+L1V/Qc8jbBVUBnZgYmGYaSYUETJz/M8DPM0jGzDRMMwUkworuS2IM9Ew8guTDQMI8WEwhEKBib3xRLhFp4ysgMTDcNIMYMS4c7TsAl+RrZgomEYKSZxItxEw8gOktnutVBEXhSRV0Vku4h8ztmrRWS9iOx2z1VxbW4XkQYR2SUiq+Lsy0Rkq/vsbrftK25r2IedfaOIzItrs8Z9x24RWYNhZDj9ET259pTzNEJhm9xnZAfJeBpB4B2qeiGwFLhaRFYCnwGeVdVFwLPuPSKyGG+71vOAq4Gvua1iAe4BbsHbN3yR+xzgZqBVVRcCdwF3ur6qgTuAS4AVwB3x4mQYmYiXCPf+y+fneR5HbN9wwzjdGVU01KPLvQ24hwLXAA84+wPAte71NcD3VTWoqvuABmCFiMwAylX1eVVV4MEhbWJ9PQpc6byQVcB6VW1R1VZgPSeFxjAykvgZ4Xk+79lEw8gWksppiIhfRLYAx/Au4huBelVtAnDPde7wWcDBuOaNzjbLvR5qH9RGVcNAO1AzQl/GBGk41snnfrqdqK2JlFJUdVAiPJBnomFkF0mJhqpGVHUpMBvPa1gywuGSqIsR7ONtc/ILRW4RkU0isqm5uXmEoRkxntl5jG9v2M+J7tBUDyWrCDlxyHeJ8FhCvD9i4mxkB2OqnlLVNuDXeCGioy7khHs+5g5rBObENZsNHHb22Qnsg9qISB5QAbSM0NfQcd2rqstVdXltbe1YflLO0hOKuOfwFI8ku4iJQ8zTiE3yM0/DyBaSqZ6qFZFK97oIeCfwOvA4EKtmWgM85l4/Dqx2FVHz8RLeL7oQVqeIrHT5ipuGtIn1dR3wS5f3eBq4SkSqXAL8KmczJkivE4vuoJWCppJQOOZpuPCUiYaRZeQlccwM4AFXAeUDHlHVn4nI88AjInIzcAC4HkBVt4vII8AOIAzcpqqxK9OtwP1AEfCkewDcBzwkIg14HsZq11eLiHwBeMkd93lVbZnIDzY8Yp5Gt3kaKWVANFypbZ6Fp4wsY1TRUNXXgIsS2E8AVw7TZi2wNoF9E3BKPkRV+3Cik+CzdcC60cZpjI3emGgETTRSyUnRcJ6Gq54K2Sq3RpZgM8JzlJiHEfM4jNQQinjnM5YA9/mEPJ8QjppoGNmBiUaO0mOeRlqIzfyOLVgIXl7DwlNGtmCikaNYeCo9DJTcDhINsfCUkTWYaOQoJxPhFp5KJSerp/wDNs/TMNEwsgMTjRylt9/maaSDoYlwMNEwsgsTjRylx+ZppIWYOMQS4QCBPCFsOQ0jSzDRyFF6gpbTSAfBYTyNkHkaRpZgopGDqCo9A+Ep8zRSSUwc4qun8i08ZWQRJho5SCgSJeJWt7UZ4aklUSI8zy9WcmtkDSYaOUhvnHfRYzmNlBITjUBeXE7DPA0jizDRyEHiQ1JdltNIKf2RwQsWgomGkV2YaOQgscopv0+s5DbFJCq5zbcZ4UYWYaKRg8Q8jeqSfJvcl2KGmxFunoaRLZho5CAx0agtLaDHwlMpJRg+NTyV5/fZMiJG1mCikYPEEuG1ZQV0hyK2T3gKCYWjBPyCt8+Yh5XcGtmEiUYOEvM0ppUWACeXFDEmTn8kOsjLAC88FTZhNrKEZLZ7nSMivxKRnSKyXUQ+5eyfFZFDIrLFPd4T1+Z2EWkQkV0isirOvkxEtrrP7nbbvuK2hn3Y2TeKyLy4NmtEZLd7rMGYMLHk97SyfMDmaqSSUDg6KJ8BrnrKwlNGlpDMdq9h4O9U9WURKQM2i8h699ldqvqf8QeLyGK87VrPA2YCz4jIWW7L13uAW4AXgCeAq/G2fL0ZaFXVhSKyGrgT+JCIVAN3AMsBdd/9uKq2Tuxn5zbxOQ1wczXKpnJE2UNC0cjzEbLqKSNLGNXTUNUmVX3Zve4EdgKzRmhyDfB9VQ2q6j6gAVghIjOAclV9XlUVeBC4Nq7NA+71o8CVzgtZBaxX1RYnFOvxhMaYAD1xOQ2wuRqpJBRJIBo+q54ysocx5TRc2OgiYKMzfUJEXhORdSJS5WyzgINxzRqdbZZ7PdQ+qI2qhoF2oGaEvowJ0BsKI+KV3IKtP5VKvER4gvCUiYaRJSQtGiJSCvwQ+GtV7cALNS0AlgJNwJdjhyZoriPYx9smfmy3iMgmEdnU3Nw84u8wPJEoCvgpLfCik5bTSB2hRInwPJ8tjW5kDUmJhogE8ATju6r6IwBVPaqqEVWNAt8EVrjDG4E5cc1nA4edfXYC+6A2IpIHVAAtI/Q1CFW9V1WXq+ry2traZH5STtPTH6E430+JEw1bfyp1hMLRQSvcwsml0b2orGGc3iRTPSXAfcBOVf2vOPuMuMPeD2xzrx8HVruKqPnAIuBFVW0COkVkpevzJuCxuDaxyqjrgF+6vMfTwFUiUuXCX1c5mzEBeoJhivL9FOd7K7HanhqpI1EiPN9tyGRlt0Y2kEz11GXAjcBWEdnibP8I3CAiS/HCRfuBjwKo6nYReQTYgVd5dZurnAK4FbgfKMKrmnrS2e8DHhKRBjwPY7Xrq0VEvgC85I77vKq2jO+nGjF6QhGKA3kWnkoDoUiUwsBg0chz4ar+yKn5DsM43RhVNFT1ORLnFp4Yoc1aYG0C+yZgSQJ7H3D9MH2tA9aNNk4jeXr7I87TcOEpS4SnjFA4Slnh4D+rmFD0hxXyp2JUhpE67LYnB+kJRSgp8JOf5yPgFwtPpZBgOEJhnn+QLRae6o9aBZVx+mOikYN41VPe3XBxfp6JRgoJhk8NTwXiwlOGcbpjopGD9IbCA0nwiqIA7b39Uzyi7KGvP0LBEE8jLz48ZRinOSYaOUh3KDIgGlUl+ZzoDk3xiLKHYDhKwSmehheeCpmnYWQBJho5SG/IS4QD1JTk02KikTKC/VEKA0NzGhaeMrIHE40cpK8/MnBhqzbRSBmqSl84knByH2Czwo2swEQjx4hElXBUByp8alx4ymYrT5z+iKLKqaLh3lt4ysgGTDRyjNi2o7G4e3VJPqFw1PYKTwHBsHcOh4anAj5XcmuiYWQBJho5RuzCFrsbjq1029JlIaqJ0tfvBHkYT8NEw8gGTDRyjKDzNGLrI9WUOtHoMdGYKCcFeYinYYlwI4sw0cgxggN3w7FEuLcRU0t3cMrGlC0Eh4T+YsRKbvstEW5kASYaOcYp4aliz9M4YeGpCdPXn9jTsJJbI5sw0cgxBu6GY6IRC09Z2e2EGc7TyDPRMLIIE40c4+SFzbsbLsn3Fi400Zg4weES4bHwlC0jYmQBJho5xtDwlIgMzNUwJsZwJbcD4Slb5dbIAkw0coyh1VNgs8JTxbAltwMLFppojIUtB9uI2m6HGUcy273OEZFfichOEdkuIp9y9moRWS8iu91zVVyb20WkQUR2iciqOPsyEdnqPrvbbfuK2xr2YWffKCLz4tqscd+xW0TWYEyIRCEUE43UMGzJ7cA8DbsAJsu2Q+1c+78b+NWuY1M9FGMIyXgaYeDvVPVcYCVwm4gsBj4DPKuqi4Bn3XvcZ6uB84Crga+JSOyv6B7gFrx9wxe5zwFuBlpVdSFwF3Cn66sauAO4BFgB3BEvTsbYSXRhs0ULU0PMiztlu1efrXI7VnYc7gCg4VjXFI/EGMqooqGqTar6snvdCewEZgHXAA+4wx4ArnWvrwG+r6pBVd0HNAArRGQGUK6qz6u30NGDQ9rE+noUuNJ5IauA9araoqqtwHpOCo0xDkLhRJ5GgYlGCggOU3Jrk/vGTkOzJxYHW3umeCTGUMaU03Bho4uAjUC9qjaBJyxAnTtsFnAwrlmjs81yr4faB7VR1TDQDtSM0JcxThKVhdaU5tMVDA94Icb4GK7k1u8T/D6xVW7HwO6jnQAcaOmd4pEYQ0laNESkFPgh8Neq2jHSoQlsOoJ9vG3ix3aLiGwSkU3Nzc0jDM0YuLD5T94Nx9afsgl+E2MgPDXE0wCv7NY8jeTZ7cJSB1vM08g0khINEQngCcZ3VfVHznzUhZxwz7GMVSMwJ675bOCws89OYB/URkTygAqgZYS+BqGq96rqclVdXltbm8xPylkGchpxd8NzqooB2H+ie0rGlC309UcQOTkvI56A32c5jSTpCYU51NZLvt9HY2sPEaugyiiSqZ4S4D5gp6r+V9xHjwOxaqY1wGNx9tWuImo+XsL7RRfC6hSRla7Pm4a0ifV1HfBLl/d4GrhKRKpcAvwqZzPGSax6KjZ3AGBBXQkAe5tNNCZCMBylIM+HKwocRMDvM08jSfY2d6MKl5xZTX9EOdLRN9VDMuLIS+KYy4Abga0issXZ/hH4IvCIiNwMHACuB1DV7SLyCLADr/LqNlWNBctvBe4HioAn3QM8UXpIRBrwPIzVrq8WEfkC8JI77vOq2jLO32rgXdjy/T58vpMXtunlhRTn+000JkiwP3JKEjxGwC82IzxJYhVT7zinjt/tPs6BEz3Mqiya4lEZMUYVDVV9jsS5BYArh2mzFlibwL4JWJLA3ocTnQSfrQPWjTZOIzlC7m44HhFh/rQS9jRbeeNECIajp5Tbxgj4fTYjPEl2H+vE7xPevsgLNR9s6eHSBTVTPCojhs0IzzGC4cgp1T0AC2pL2XvcRGMi9I3gaeT7fTa5L0kajnVxRk0xZ9QU4/cJBywZnlGYaOQYsfDUUM6sLaGxtXdgeW9j7AQTeHExAn6fLSOSJK09/dSVFRDw+5hRUWiikWGYaOQYwXB0YIXbeBbUlqJqFVQTwQtPJfY08qzkNml6QmFK8r3I+dzqYhONDMNEI8fwkrWJPQ2APcdMNMZL3zDnFqzkdiz0BCMU5XviO72ikOZO21UykzDRyDFCkcQhlDOnlQKw15Lh48bz4hL/SeX7fTYjPEm64zyNqmJbFy3TMNHIMYL90YTJ2qJ8P3Oqi/j51ibaeuyPdDwEw5GEs8EBAnkWnkqWnmCE4oLYHvb59PZHLNeWQZho5BjBcGTQXhrx3PFH57G3uZsPfeMFekP2RzpW+vqH9zRscl9yqOopngZAq93IZAwmGjnGSBU+71xcz5euv4BdRzvZ9KbNoRwrwfDwJbd5Ph8hC0+NSjAcJaoMeBpVxQEAWrv7p3JYRhwmGjnGSHF3gLctnAbAriOdkzWkrMEL/Q2T08gTQraK8Kj0OA93wNMoMU8j0zDRyDFGuhsGqCktYFppPm8cNdEYK339kWFLbgvz/AOr4BrD0x0MA1CcH/M0TDQyDRONHCPRMiJDOau+jF1HrYpqrIwU+isI+Af2EDeGZ8DTKIh5GrHwlIlGpmCikWOMdGGLcfb0MnYf7SRqS1InjaqOeG6LAn6rAEqC7tBwnoblNDIFE40cI9gfHbZ6KsbZ9WX0hCI0ttquackSm7iXaLY9ePuGm2iMTk/QO0fFLqcR8PsoK8izuRoZhIlGDuHdDY+c0wA4a3oZALssr5E0sdDTcJ5GYcBPOKqErex2RIZ6GuAlw23uUOZgopFDhKNKVIe/sMU4q94TDUuGJ8/JHRGH9zQA+iwZPiI9TjRiOQ3wym5bLDyVMZho5BCh2P7gI5TcApQW5DG7qojXrew2aWI7IhaOkNMAbNLkKHQHYyW3gz0NS4RnDsls97pORI6JyLY422dF5JCIbHGP98R9druINIjILhFZFWdfJiJb3Wd3uy1fcdvCPuzsG0VkXlybNSKy2z1i28Ea4yRW8jlaeApgUV0pu83TSJrRPI2Y3fIaIxPzNIrjPI3q4nwruc0gkvE07geuTmC/S1WXuscTACKyGG+r1vNcm6+JSOyv6B7gFrw9wxfF9Xkz0KqqC4G7gDtdX9XAHcAlwArgDrdPuDFOYhe20RLhAAvrStl7vJuIVVAlRTI5DTj5b2AkJuZpFMWJb2WxeRqZxKhXD1X9Ld6+3clwDfB9VQ2q6j6gAVghIjOAclV9XlUVeBC4Nq7NA+71o8CVzgtZBaxX1RZVbQXWk1i8jCQJjnJhi2dhXSmhcJTGVtvLIBliXtzwk/tcTsPmaoxIb3+EooAff9we9tUlAbpDERPcDGEiOY1PiMhrLnwV8wBmAQfjjml0tlnu9VD7oDaqGgbagZoR+jLGyVjCUwvrvKXSG47ZJL9kCLqw07DzNFyMvtfCUyPSHQxTUjD4/2elm6vRZsnwjGC8onEPsABYCjQBX3Z2SXCsjmAfb5tBiMgtIrJJRDY1NzePNO6cZiARnoynUetVUJloJEdwlHNbaDmNpOgJndyAKUa1rT+VUYxLNFT1qKpGVDUKfBMv5wCeNzAn7tDZwGFnn53APqiNiOQBFXjhsOH6SjSee1V1uaour62tHc9PyglOJmtH/2evKA4wrbTARCNJ+gY8jeHXnvKOs/DUSHQHTy6LHiM2K9wm+GUG4xINl6OI8X4gVln1OLDaVUTNx0t4v6iqTUCniKx0+YqbgMfi2sQqo64DfunyHk8DV4lIlQt/XeVsxjgZS3gKYGFdCQ22k19SxNZMKs4fZZ6GeRoj0hOKnHIOT64/ZeGpTCBvtANE5HvAHwLTRKQRr6LpD0VkKV64aD/wUQBV3S4ijwA7gDBwm6rG/kpuxavEKgKedA+A+4CHRKQBz8NY7fpqEZEvAC+54z6vqrbJwwQYS/UUeHmNx7YcZtuhdvw+4dwZ5ekc3mlNT/9oomE5jWToDoUpLRh8WaopKQCgpdv2Cs8ERhUNVb0hgfm+EY5fC6xNYN8ELElg7wOuH6avdcC60cZoJMdYqqcAFtaW0tkX5o+/+hyzKov43d9fgZteYwyh180vGBqPjzFQcmuiMSI9wQh1ZQWDbNUl+fgEmjtNNDIBmxGeQ4yWrB3KhXMqAVhQW0pjay/bDnWkbWynOyfDU4nvw06GpyynMRLxW73G8PuE6pICmrssp5EJmGjkECeXEUkup3HR3Cqev/0dPPqxS8nzCT/f2pTO4Z3W9Ia8vdfj5xfEY9VTydETigxs9RrPtNJ88zQyBBONHGKgeipJTwNgRkURlcX5vHXhNJ7c1oRXo2AMpbf/1ARuPAG/jzyfWE5jFBJVTwHUlhVwvMtEIxMw0cghYuGpZBPh8bxnyXTePNHDjiYLUSWiJxSheBQPrtB27xuRSNTbyCpRiK+2tMA8jQzBRCOHGGtOI55LF9QAsP2wiUYiehNMShtKYcBHny2FMSw9CfbSiFFbVkBzV9A83QzARCOH6OuPIAL5/rH/s8+sLMIn0Nhia1EloicUHjYJHqPQtnwdkYFigoQ5jQJC4SidwfBkD8sYgolGDtEdjFCSnzeustmA38eMiiIOmGgkJNHyF0Mx0RiZbicIw+U0wMpuMwETjRzCuxtOrnIqEXOrizlo+4YnZLREOMT2CbecxnCMNKs+JhrHTTSmHBONHKI7FBm0jeZYmVNtnsZw9IQig/aASERhnnkaI9EVjOU0Tv0/Oq3UeRpWQTXlmGjkED3BiXsazZ1B27I0AckkwovyTTRGoqPXW1uqoihwymcWnsocTDRyiESzbcfCnOpiANuYKQHJhP4K8vz0WnhqWNqdaFQWnyoalUUB/D6xuRoZgIlGDtE7zGzbZImJxkETjVPwVmcdrXrKZ2tPjUBMNMoTeBo+n9is8AzBRCOH6E6w7PRYmFPlicaBEyYa8UTdpLRRcxpWPTUibT39+ATKhsm7TSst4LitPzXlmGjkEF5OY/zhqWml+RQF/FZBNYTeUZZFj1EU8NMXtvDUcLT1hqgoCuAbZv2u2jKbFZ4JmGjkEN2hCCUT8DRExCqoEjDaBkwxCgM+KyIYgfbe8MB+4ImYXl7I/hPdAzPHjanBRCOH6AmFKZ5AyS14FVQWnhpMTAiKkpkRHo7YUhjD0NYTSpjPiHH98jl09oVZ99y+AdvT249wx2Pb+PIvdtl5nSRGFQ0RWScix0RkW5ytWkTWi8hu91wV99ntItIgIrtEZFWcfZmIbHWf3e22fcVtDfuws28UkXlxbda479gtIrEtYY1xEApH6Y/ohDwNgHOml9PQ3GWx+Th6+odfMymewoAfVQhFLESViPbefipHEI1lZ1TxrsX1fOM3e2ntDrFpfwsffWgzD73wJv/zy4aBeR5GeknG07gfuHqI7TPAs6q6CHjWvUdEFuNt13qea/M1EYn9Jd0D3IK3b/iiuD5vBlpVdSFwF3Cn66sab2vZS4AVwB3x4mSMjZOLwU3M0zh/dgWRqNpqt3H0DHgao4sG2EZMw9HW05+w3DaeT686m+5QmE89vIW1T+ykvryAz13jbQja0m1J8slgVNFQ1d/i7d0dzzXAA+71A8C1cfbvq2pQVfdh3A6mAAAgAElEQVQBDcAKEZkBlKvq8+r5kA8OaRPr61HgSueFrALWq2qLqrYC6zlVvIwkiV3YSiZQcgtwwewKALY2tk94TNlCLDw1+tLosd37zEtLxGieBsBZ9WX8+wfO57dvNPPKgTb+5p1nMbuqCMAqqyaJ8d521qtqE4CqNolInbPPAl6IO67R2frd66H2WJuDrq+wiLQDNfH2BG2MMZIqT2N6eSHTSgvYeshEI0bSnkae7d43HJGo0tHXn3A2+FA+9Ja5hMJRNu5r4bpls9nZ1AmYpzFZTOwKciqJauV0BPt42wz+UpFb8EJfzJ07d/RR5iDdweQqfEZDRDh/Vrl5GnGMtA9EPDFRsfDUqXT29aMKFSNUT8Vz46XzuPHSeQBUl3ptTths8UlhvNVTR13ICfd8zNkbgTlxx80GDjv77AT2QW1EJA+owAuHDdfXKajqvaq6XFWX19bWjvMnZTfdKfI0AM6fXcnuY51W+uhIvnrKwlPD0dbjlhBJwtMYSk2JEw3zNCaF8YrG40CsmmkN8FicfbWriJqPl/B+0YWyOkVkpctX3DSkTayv64BfurzH08BVIlLlEuBXOZsxDnqCqclpAFwwq4Kowg7bxQ+Im9yXxCq38ccbJ2kbYd2p0SgM+CnJ93PCchqTwqi3nSLyPeAPgWki0ohX0fRF4BERuRk4AFwPoKrbReQRYAcQBm5T1dhfyK14lVhFwJPuAXAf8JCINOB5GKtdXy0i8gXgJXfc51V1aELeSJJUehrnziwH4I2jXSyfVz3h/k53ks1pFAQspzEcIy1WmAw1pQWc6Lbw1GQw6hVEVW8Y5qMrhzl+LbA2gX0TsCSBvQ8nOgk+WwesG22MxuikqnoKoK6sABE40tE34b6ygd5QBJ+Mvvd6kZXcDktbj+clJJMIT0RNab4lwicJmxGeI5xc6mLinkbA76O2tIAj7bYGFZxc4Xa0bXQtpzE87QN7aSSXCB9KTUm+ldxOEiYaOUJPMLkKn2SZXlHIkQ4LBwD09odHDU3BySW/O/r60z2k045YInzcnkZJgVVPTRImGjlCdyhCvt9HwJ+af/Lp5YXmaTh6klxyvrIogAiWsE1Ae28/Jfl+8kcJ8Q1HtQtP2fpT6cdEI0fwFitMjZcBztNot5wGJLc/OECe30dFUYDWHhONobT1JDexbzhqSvIJR5WOXisDTzcmGjlCdzAyoa1ehzK9opCOvrDN1cDtiJhk2K+62BK2iWjvDSU9sS8R00q9PcSPWwVV2jHRyBGS2cN6LEwvLwQwbwPv3CaT0wCoKsk3TyMBnqcx/puaajfBzwQ5/Zho5AjdociE99KIZ3qFiUYMLzyV3LmtKs63nEYCmtr7Bm5ExkONLSUyaZho5Ai9ofCE99KIZ8DTsLkadPT2U57kXXKNeRqnEI5EOdLRx2y3B/14qClx4SkT5LRjopEjdAcjKZmjESPmaTTluKfR7y54syqLkjq+qiSf1u5+q/KJ40hHH5GoMqsquXOYCAtPTR4mGjlCTyicktngMYrz8ygvzONojnsaR9r7iCpJi0Z1SYBQJEq37RU+wKFWr3R79gREIz/Pq0xr7rTwVLox0cgRusdQ4ZMs0ysKc97TONTmXfCSvUuuchVCrXZHPECjE41khXc46soKONaZ2/8fJwMTjRyhJxhOaXgKYHpFUc57GofbxnbBq7ZlvE8hJrwzJyga9eWFHDNPI+2YaOQA0ajS0x9JaSIcYFZlIW+e6Mnp+HwstJLsBa+qxDyNoTS29lBbVjCwh/p4qSsr4JgtbZN2TDRygA63K1r5BGbcJuLC2ZW09/azp7k7pf2eThxq62VaaX7SF7waS9iewqG23gmHpgDqygs51tmX0zcxk4GJRg5waIwhlGSJ7aWx+c3c3eZkrBe8AU/Dym4HaGztnVASPEZdWQH9EaW1xxaETCcmGjnA4TYv7zDRmPFQFtSWUFUcYNP+1pT2ezpxqLV3TKWiZQV55PnEchqOaFRpauubULltjHo3d8iS4ellQqIhIvtFZKuIbBGRTc5WLSLrRWS3e66KO/52EWkQkV0isirOvsz10yAid7stYXHbxj7s7BtFZN5ExpurNLnVaGdUjn/GbSJEhGVnVLPpzdwUDVUds6chIm6uhokGQHNXkFAkOqGJfTHqyr0Jfkctr5FWUuFpXKGqS1V1uXv/GeBZVV0EPOveIyKL8bZyPQ+4GviaiMQCwfcAt+DtKb7IfQ5wM9CqqguBu4A7UzDenONQWy/5fh/T3KzZVLJ8XhX7jndzPAeXbzjeFSIYjo457GeLFp6ksbUHgNkp8ILry5ynkeMVfekmHeGpa4AH3OsHgGvj7N9X1aCq7gMagBUiMgMoV9Xn1ctgPTikTayvR4ErZbTt0YxTONzWx4zKQny+1J+6t8zzHMlcDFGdnKMxtrvkqhJbHj3GvuOeaMypTp2nYWW36WWioqHAL0Rks4jc4mz1qtoE4J7rnH0WcDCubaOzzXKvh9oHtVHVMNAO1ExwzDlHU1svMytSm8+Icd7MCnwCO5s60tJ/JnNonJPSakoKLKfheL2pg4I8H/NqJi4ahQE/5YV55mmkmYmKxmWqejHwbuA2EfmDEY5NdJurI9hHajO4Y5FbRGSTiGxqbm4ebcw5x+G23pTnM2IUBvzMqS6mobkrLf1nMq81tpHv93FmbcmY2tWWFXjLj0StNHTnkQ7Onl5GXop2lKyzCX5pZ0L/Uqp62D0fA34MrACOupAT7vmYO7wRmBPXfDZw2NlnJ7APaiMieUAFcEp9p6req6rLVXV5bW3tRH5S1hEe44J642FBbSl7juWeaLywr4WlcyrHPCltUX0pPaHIQHgrV1FVdjZ1cu708pT1WV9ekPOrFKSbcYuGiJSISFnsNXAVsA14HFjjDlsDPOZePw6sdhVR8/ES3i+6EFaniKx0+YqbhrSJ9XUd8Eu1mTtj4mhnkKimvtw2ngW1Jew73k0kh+6cu4Jhth1q55Izq8fc9qz6MgB2H+tM9bBOK5o7g7R0hzhnRlnK+qwrO308jV1HOrn1O5u59N+f5XM/3U7vabKI5UQ8jXrgORF5FXgR+LmqPgV8EXiXiOwG3uXeo6rbgUeAHcBTwG2qGjtLtwLfwkuO7wGedPb7gBoRaQD+FleJZSRPU4rW9RmJhXWlBMPRgRh/LrD5zVYiUWXF/HGIRp13kXzjaO55Z/HscHmwc2ekztOoK/eWEjkd7i3/+SfbeK7hOGdPL+PbG/az+psvTPWQkmLcK9ip6l7gwgT2E8CVw7RZC6xNYN8ELElg7wOuH+8YjbjF4CrSk9MALzwFsKe5i7kpSGhmMqFwlNePdPDC3hPk+YRlZ1SN3mgIFcUB6soK2J3jovH6Ec/TSmV4anp5IaFIlBPdoYF9wzORQ229vLi/hb9711n81ZWL+Oovd/Ofv3iD413BjB432IzwrCc2G3xGmnMa4IlGtnP/7/fxvq9u4Ou/2cP5syvGvXLwWfVlOR+e2tnUwcyKQiqKU7cm2sI67//iG0cz+9w+vsVL216z1CsUXTrHu/nI9HGDiUbW03Csi5qSfEpTuD/4UKpK8qkpyachB5Lhv97VzMyKQlbMq+aGFXPH3c+i+lJ2H+3K2QoqVWVrY3tKQ1MA5ziv5fWmzL74PrblEBfPrRzwzM+a7sTuSGaPG0w0sp6X9rewfN7YQyhjZUFtadZ7Gr2hCJv2t/Ke82fw8Ecv5YPL54zeaBjOqi+jtz93Kqhea2zjk997he5gGPBCU3uPd/OHZ6e22rG2rIBppQW8fiRz5w1tO9TO60c6B7wMgNrSAqqKA+w6DUKWJhpZTFN7LwdaelgxP/3zIRfUlfBGlt85b3qzhVAkymWLpk24r7PqvTvLXAlRfeO3e3n81cP82xM7AXj81cP4fcJ7zp+R8u86d0YZOzPY0/jOC29SGPBxbZxoiAhn1ZdZeMqYWl7c501puWQcFT5jZcX8atp7+3ntUHvav2uqeK7hOAG/sGLexM/nWfVl+H2SE8uvdPb188yOo1QVB/juxgP8+JVGHt9ymLcvmkZNGpK+50z3Lr7hSDTlfU+U9t5+frLlENdcOOuUXM7Z08t440hnxld+mWhkMS/ua6G0IC/lceNEXHF2HX6f8MyOo2n/rqliQ8NxLppbRUkK8kNlhQFWnlnNU9uPZPxFYqI8te0IwXCUez6yjAtmV/A3D7/KobZe3nfhzLR83znTywmGo+w/0ZOW/ifCDzc30tcf5cZLzzjls7Pqy+gMhmlqz+zJiSYaWcyL+7x8hj8NCxUOpbI4n+VnVPHMzuwUjQMneth2qIPLz0pdDP7qJTPY29zN7iwvIPjJlkOcUVPMJfOr+eGtb+UL1y7hvRfMYNV509PyfbHJglOd1whHojyxtYl+5/GoKt954U0umlvJklkVpxx/9nRv3LsyPERlopGlnOgKsvtYF29JQSglWd61uJ7Xj3RysCXz7vAmyqMvNyIC779o1ugHJ8mq8+oRgSe2NqWsz0zjW7/by4aGE3xw+RxEhIDfx40rz+B///TilHhsiVhYV4rfJ2w/PLWi8YPNjXz8uy/zwO/3A/D7PSfYe7ybG1ee6mXAyUmfm/Zn9k6YJhpZyu/3nADgrQsmb1Hgd55bD3jhiGwiGlV+uLmRty2cltKZ9XVlhSw/o4ontjZlZYjqsS2H+Nef7+TdS6bzscsXTNr3FuT5Pa9mcyN9/VOzNIeqcv+G/QB87dd76Ozr56Hn36SqODBs8r+iOMBVi+u5f8N+mhMshdIdDPPszqNEo4qq8vyeE3z9N3t4bvdxwBObyRAcE40sZUPDccoK8jg/gRucLuZNK2HFvGrWbdhHKJx5Scjx8Nzu49z59OscauvlumWzR28wRq5fNoc3jnaxPstyQa8caOXTj77GJfOr+crqpZMSIo3nr96xiGOdQb734oFJ/d4Yz+89wa6jndy48gxaukN8+FsbWb/zKB98y5wRF7j8zLvPIRiO8t/PvnHKZ//82DZufmATP37lEPc9t48bvvkCX3zydW5+4CWe3n6EWx7azP/9yba0rwFnopGlbNhznJULalK25HSyfPyKBTS19/HjVxpHPzjD+cGmg3zkvo184zd7mV1VlJYY/AcunsWZtSX8x9O7MrLaZzxsO9TOXz64mfryAu75yDIK8sa2CnAquHRBDZfMr+Zrv97DrimYMPfA7/dTXZLPP733XG5YMZcTXSHee/4Mbnn7mSO2O7O2lA9fMpfvvXiQhrhy7N++0cyPXj5EYcDHvz+5ky//4g2uOLuWZ//uciqKAnz0oc1Eoso9H1mWdoE20TiNCIaTc7UPtvRwsKWXyyYxNBXj8rNqWTKrnHt+vee0Drm8sPcEt/9oK5ctrOHVf7mK3/39FWNeAj0Z8vw+/n7V2TQc6+I7L7yZ8v4nm417T/DBbzxPvl/49p+toLokf8rGcvt7ziXYH+Hd//1bvv6bPZP2vce7gjy78xjXLZtNYcDPv3/gfDZ85h3cfcNFSZUYf/LKRRQH/HzxydcBaO0OcfuPtnJmbQkP/p9LON4Vwiew9v3ns6C2lP/98MUsrCvl6x9ZxvxpY9vbZTyYaJwm3L9hH+f9y9N85Zk3hnU/+/oj/OvPdgxMoHpbCiahjRUR4caVZ7D/RE/GV4EMR1cwzN898ipzq4v52oeXUVEcIJ27DK86bzqXn1XL2id28vKB03feRiSq/NNPtlFbVsCPb7tsYB2oqWLpnEp+8+kreOe59Xzp6V1sm6Q5RD955RDhqHL9OMOZNaUFfPyKhTyz8xj3b9jHpx7eQnNnkP/64FJWzK/mix84n69++OKB/Npb5lXzzN9ezqWTdJNoonEa8MTWJj73sx3UlxfylWd2s/Rzv+A9//07Nu49Mei4b2/Yz7ee28eT244wt7p4YCHByeZti7yy1A0NJ0Y5MrOIRpXXj3Tw+Z9u53B7L1+6/gIqilK3mN5wiAj/vXopMyqK+Ph3Xp6y5O1EefzVQzQc6+Ifrj6H+vL0rao8FqpK8vnSdRdSXZLP3z/6WtrPrary6OZGLpxTyaL68e8T8ueXzeOc6WV89qc7+O0bzdzxvsUsnVMJwOoVc7ni7LpRekgfJhoZTHtvP7f9v5f5+HdfZumcSp7528v5+kcu5gMXz6InFOYj923kB5u8bddPdAX52q8auPKcOrb8y7v46Sfelta745GYVVnEmdNK2NBwfEq+fzyoKrd+dzNXf+V3PLKpkT976zyWnTF55cqVxfnc+ScXcKSjjx++fHrkg8KRKC/sPUFHXz/NnUHuWr+bxTPKuTpN8y/GS0VxgH+9dgk7mjr4k3t+n7aS8J+8coi/fHAzrx/pnHDRRGHAzxOffDuP3XYZ9964jD+dwOKYqSZ9S5+eZvT1R3h0cyNXnFOX1q1Rx8JXf7mbp7Yd4VNXLuJjly+gKN/P1UtmcPWSGbT3eILy6Udf4+UDbbxxtJOe/gi3v+ccKounLo4c460La/jRy4foj0QJTHIyPp5wJMrLB9rw+yAShaMdfWw71E5HX5i6sgKuXjKdc2eU8+jmRp7efpSP/sGZXHluPcvHsU/GRFl5ZjUXzK7gvt/t44a3zMU3yRVHI7Fpfwvf/N1ewhHlg2+ZQ1tPiHXP7WfX0U6qigPk+X109PZz55+vyKhxx1h13nTuW7Ocv3l4C9d9/ff84KNvTWrvl6e3H+Gu9W9wpKOPmpJ8Ll1Qw7/80Xnk5/lOOe6vH97CrMoiPnDxLD6Qgvk8Pp9wofMuMgk5HZKVInI18N+AH/iWqn5xuGOXL1+umzZtGvN3HG7r5fIv/YoPX3IGn33feeMfbAJieyHn+YVZlUVJTWrqCoa59N+e5fKza/nqn16c8Jj+SJTPPr6d7248QF1ZAX/zrrMmtFx3KnlqWxMf+87L/OBjl07qBMN4WrtDfOJ7L58SJsv3+ygvyqOlO0RUoa6sgNaeEBfPreJ7f7lySi96P331MH/1vVe460MX8v6LUl/iO1a6g2E+99PtPLKpcUAcYnMI5lQX8dE/WMDT249wrCPIXR9ayuKZ6V+yZiK8fqSD1fe+QGlBHt/9i0s4o2b4xPG+49289+7fMbOyiEvmV3O0o49ndh7jxpVn8McXzuSpbUfo7Y8QCkd5evsRFtSW8PBHL01LwcRkICKbVXX5qMdlumiIiB94A2/r2EbgJeAGVd2R6PjxigbAp3/wKj997TAb/uEdCascVJUtB9sI+H0snlE+6sVlx+EOfr/nOI9ubhzYpaww4ON9F87kkvk1zK8t4cxpJQk9g29v2MfnfrqDH3/8rVw0d/i7XlVl97Eu5k8rmdI7+qG09/Rz0Rd+wdsX1fKl6y/gjSNdLKovHRTrVlWa2vsoLwpQWuBdxB98fj9PbG1iYV0pb5lXzdn1ZayYXz2m0uFIVHl080H+8xdv0N7bzz+951zOqCnG7xOmlRawoLaU/DwfbT0hfvjyIXYd6aA4P4+PXb6A6Wnc4TAZwpEof/zVDexs6uBPL5nLVYvrWXlmDYUBP8e7grT19FNXXkBZQV5aw4/hSJQntx3hzqe8OSq3Xr6Av3rHInw+b3mamZVFzK8pyUivYjRea2zjpnUvAnDDirk0tfVyyZk1XH5WLXl+4YnXmth6qINXDrRyojvEU3/9dmZUeNGHf3tiJ/f+di/g/S2XFQbI9/uYUVHI3TdclNZtldNNNonGpcBnVXWVe387gKr+e6LjJyIaDcc6edddv+Xjf7iAT686h9buEM/vPcHGvSdo6+1n99GugX2Np5UW8MkrF/LWBdNoau+lqa2PYDhCSUEe582s4Lsb3+TB570SynNnlHPjyjMoK8zj93uO85NXDtMbl5CrLA5QX1ZIXzhCdzBCbyhMdyjCxXMr+dHHLxvXb8kE7t+wj3/9+U7CrtqrqjjAJ69cxI9ePjSwj0RLd4jqknz+dMVcvrPxTdp6+lkxr5oDLT0c6fAWbjtzWgnXLJ1Fe28/4WjiuQyzq4qYV1NCR1+Yb/1uL68f6eTiuZV8/polCdf5yWR6QmHufPJ1HnrhTaLqbdV71XnT+d6LBwi6SZNFAT/15QXUlRVSV17A4pleLqGpvY+9x7tp7uhj66F2Igrvv2gmh9v62NPcRVHAz4muEL39ES6eW0VNaT5RVSJRZWZlEeWFAX6w6SDP7DxKR1+Yc6aX8bn3ncclZ05++XY62X+8m798cBMNzd4mZce7QoM+ry8voKQgj39+72KuOOdk0jkcibL2iZ3MqCjkxpXzKMo/Pb2KRGSTaFwHXK2qf+He3whcoqqfSHT8REQD4GMPbeap7UeoKg7Q1tuPKhTn+6ktK6CyOJ8PLZ9DUb6PR15q5Pm9I1cH/cXb5vOXf3DmKZUk/ZEoB1t62He8m33Hu9l7vJvjnUGK8v0U5+dRnO+nJN/PH184c0IVGJnAloNt/GL7ERbPLOcrz+ym4VgXc6uLefuiaUSiyln1ZTz26mFePdjGRXMr+eIHLuDs6WWoKs2dQV7a38rdz+5m19FOivP9FOSd6nFE1SsaiDGnuoh/uPoc3nv+jCkrBkgFXcEwL+1r4UtP72JHUwfvvWAG7zy3jubOIMc6ghztDHKso4+jHX2nrOgqAovqSukOntzoaUZFIcFw1Asz+Xy8cayTRH/+pQV5vHvJdK48t553La6f9Nnck0U0qgTDUQoDPl5rbGfb4XZ6ghHeurCG82aeXjcaqSCbRON6YNUQ0Vihqn8Vd8wtwC0Ac+fOXfbmm+OfJNXR188PNzey43AHc6qLuWxhDRfMrjwl9KOqbGg4wfGuIDMqCplZWURhwE97b4iXD7Qxu6qIty6Y/HkSmUxnXz+/33OCK86uG5RIjESVLQdbWTon8Yq8qkpff3TEu7qW7hCNrT0UBfycUVNySqLydCYSVU50B6krGz50duBED7/Z3czc6mLOmV5GdUk+Ab9v4NzOrio+5ealo6+fvlAEn08QYP+Jbo60B7n87Nq0bg9sZCbZJBqTFp4yDMPIVZIVjdPhduwlYJGIzBeRfGA18PgUj8kwDCMnyXgfVFXDIvIJ4Gm8ktt1qrp9iodlGIaRk2S8aACo6hPAE1M9DsMwjFzndAhPGYZhGBmCiYZhGIaRNCYahmEYRtKYaBiGYRhJY6JhGIZhJE3GT+4bKyLSCeyaQBcVQCq3+Eplf6keW4xpQCo2v8jkc5euPlN17mJk+jm085c5/aXy3E0DSlS1dtQjVTWrHsCmCba/N8XjSVl/qR5bqs7Z6XDu0jjGlJy70+Uc2vnLnP5See7G0peFp07lpxncX6rHlmoy+dyls89Ukunn0M5fZvU36WRjeGqTJrF+inESO2fjx87dxLDzN35See7G0lc2ehr3TvUATkPsnI0fO3cTw87f+EnluUu6r6zzNAzDMIz0kY2ehmEYhpEmTDSyEBGZIyK/EpGdIrJdRD7l7NUisl5EdrvnKmevccd3ichXh/R1g4hsFZHXROQpEcnqnaVSfO4+5M7bdhH5j6n4PZPNOM7fu0Rks/s/tllE3hHX1zJnbxCRu+V03oYxCVJ87taKyEER6Ur5QFNZTmaPzHgAM4CL3esy4A1gMfAfwGec/TPAne51CfA24GPAV+P6yQOOAdPc+//A2xBryn/jaXDuaoADQK17/wBw5VT/vgw8fxcBM93rJcChuL5eBC4FBHgSePdU/77T6NytdP11pXqc5mlkIarapKovu9edwE5gFnAN3sUL93ytO6ZbVZ8D+oZ0Je5R4u7yyoHD6f8FU0cKz92ZwBuq2uzePwP8SZqHP+WM4/y9oqqx/1PbgUIRKRCRGUC5qj6v3lXwwVibbCVV58599oKqNqVjnCYaWY6IzMO7I9kI1Mf+I7nnupHaqmo/cCuwFU8sFgP3pXG4GcVEzh3QAJwjIvNEJA/vD31O+kabeYzj/P0J8IqqBvEulo1xnzU6W04wwXOXVkw0shgRKQV+CPy1qnaMo30ATzQuAmYCrwG3p3SQGcpEz52qtuKdu4eB3wH7gXAqx5jJjPX8ich5wJ3AR2OmBIflRKlnCs5dWjHRyFLcBf+HwHdV9UfOfNS5/bjnY6N0sxRAVfe4EMEjwFvTNOSMIUXnDlX9qapeoqqX4q2HtjtdY84kxnr+RGQ28GPgJlXd48yNwOy4bmeT5aFRSNm5SysmGlmIyz/cB+xU1f+K++hxYI17vQZ4bJSuDgGLRSS2iNm78OKsWUsKzx0iUueeq4CPA99K7Wgzj7GePxGpBH4O3K6qG2IHuzBMp4isdH3eRBLn/HQmVecu7Ux1xYA9Uv/Aq+ZRvHDSFvd4D15Fz7N4d7zPAtVxbfYDLUAX3l3eYmf/GJ5QvIa3bk7NVP++0+jcfQ/Y4R6rp/q3ZeL5A/4v0B137Bagzn22HNgG7AG+ipuMnK2PFJ+7/3D/F6Pu+bOpGqfNCDcMwzCSxsJThmEYRtKYaBiGYRhJY6JhGIZhJI2JhmEYhpE0JhqGYRhG0phoGMYkIyIfE5GbxnD8PBHZls4xGUay5E31AAwjlxCRPFX9+lSPwzDGi4mGYYwRt5jcU3iLyV2Et4T1TcC5wH8BpcBx4M9UtUlEfg38HrgMeFxEyvCWrP5PEVkKfB0oxpvE9n9UtVVElgHrgB7gucn7dYYxMhaeMozxcTZwr6peAHQAtwH/A1ynqrEL/tq44ytV9XJV/fKQfh4E/sH1sxW4w9m/DXxSvXWrDCNjME/DMMbHQT253s93gH/E2whnvdtgzg/E72fw8NAORKQCT0x+40wPAD9IYH8IeHfqf4JhjB0TDcMYH0PX3+kEto/gGXSPoW9J0L9hZAQWnjKM8TFXRGICcQPwAlAbs4lIwO1zMCyq2g60isjbnelG4Deq2ga0i8jbnP3DqR++YYwP8zQMY3zsBNaIyDfwVh/9H+Bp4G4XXsoDvoK3DedIrAG+LiLFwF7gz539z+ZnT5wAAABdSURBVIF1ItLj+jWMjMBWuTWMMeKqp36mqkumeCiGMelYeMowDMNIGvM0DMMwjKQxT8MwDMNIGhMNwzAMI2lMNAzDMIykMdEwDMMwksZEwzAMw0gaEw3DMAwjaf4/1DtGyXXii+sAAAAASUVORK5CYII=\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": 11, + "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": 12, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2344,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": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+QldWd5/H3B0FxElEg6PBDhY3GCppZHLrQLfdHxB0gPyrgjJmwOkrVWEXijy1nNrWiE7fMKFMVU5M4y7ohMWNG1FF0TCzdRIa0UWucWQSaoFE0TjMrgygjWI2KUwVj43f/eL53fLhpbt9ubve9t/vzqrrVT5/nnHNPPzT9vefH8xxFBGZmZkdrTLMbYGZmI4MDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNcTYZjdgOH3sYx+LmTNnNrsZZmZtZcuWLW9FxJT+8o2qgDJz5ky6urqa3Qwzs7Yi6R/ryechLzMzawgHFDMzawgHFDMza4i6AoqkHZJekPScpK5M+7qk1zPtOUmfLeW/UdJ2Sa9IWlhKn5v1bJe0SpIy/ThJD2b6RkkzS2WWSerO17JS+qzM251ljz36y2FmZoM1kB7KhRExJyI6Smm3Z9qciHgcQNJsYClwNrAI+I6kYzL/amA5cGa+FmX6lcC+iDgDuB24LeuaBNwMnAfMA26WNDHL3JbvfyawL+swM7MmGYohr8XA2og4GBGvAtuBeZKmAhMiYkMUu3rdAywplVmTxw8DF2XvZSHQGRE9EbEP6AQW5bn5mZcsW6mrre159wC/+70N7Nl/oNlNMTMbkHoDSgA/lbRF0vJS+rWSfiHpB6Wew3TgtVKeXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdbW1VT/rZvOOHlY90d3sppiZDUi996FcEBFvSDoZ6JT0S4rhq1spgs2twLeA3wfUR/mokc4gytSq6zAZAJcDnHbaaX1laQln3bSOg70f/Ov3923cyX0bd3Lc2DG8svIzTWyZmVl96uqhRMQb+XUP8AgwLyLejIhDEfEB8H2KOQ4oegunlorPAN7I9Bl9pB9WRtJY4ESgp0ZdbwEnZd7quqrbfmdEdEREx5Qp/d7o2TTPXH8hX5gzjfHjin+S8ePGsHjONJ5ZcWGTW2ZmVp9+A4qkj0g6oXIMLABezDmRiouBF/P4MWBprtyaRTH5vikidgP7JZ2fcyBXAI+WylRWcF0CPJnzLOuBBZIm5pDaAmB9nnsq85JlK3W1pZMnjOeE48ZysPcDjhs7hoO9H3DCcWM5+YTxzW6amVld6hnyOgV4JFf4jgXuj4i/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDURcSjrugq4GzgeWJcvgLuAeyVtp+iZLM26eiTdCmzOfLdERE8erwDWSloJbM062tpb7x3ksvNO59J5p3H/pp3s9cS8mbURFR/2R4eOjo7ws7zMzAZG0paqW0b65DvlzcysIRxQzMysIRxQzMysIRxQzMysIRxQzMysIRxQ2pCf92VmrcgBpQ35eV9m1opG1Z7y7c7P+zKzVuYeShvx877MrJU5oLQRP+/LzFqZh7zajJ/3ZWatys/yqsOedw9w7QNbuePSc90bMLNRx8/yaiCvqjIz65+HvGrwqiozs/q5h1KDV1WZmdXPAaUGr6oyM6tfXQFF0g5JL0h6TlJXpk2S1CmpO79OLOW/UdJ2Sa9IWlhKn5v1bJe0KrcCJrcLfjDTN0qaWSqzLN+jW9KyUvqszNudZY89+svxqyqrqh65+gIuO+909r53cCjexsys7dW1ykvSDqAjIt4qpX0T6ImIb0i6AZgYESskzQYeAOYB04AngE9ExCFJm4DrgGeBx4FVEbFO0tXAb0TEVyQtBS6OiC9JmgR0AR0UWw1vAeZGxL7cZvhHEbFW0neB5yNida2fwzs2mpkN3HCs8loMrMnjNcCSUvraiDgYEa8C24F5kqYCEyJiQxRR7J6qMpW6HgYuyt7LQqAzInoiYh/QCSzKc/Mzb/X7m5lZE9QbUAL4qaQtkpZn2ikRsRsgv56c6dOB10pld2Xa9DyuTj+sTET0Au8Ak2vUNRl4O/NW12VmZk1Q77LhCyLiDUknA52Sflkjr/pIixrpgylTq67DG1MEwOUAp512Wl9ZzMysAerqoUTEG/l1D/AIxfzImzmMRX7dk9l3AaeWis8A3sj0GX2kH1ZG0ljgRKCnRl1vASdl3uq6qtt+Z0R0RETHlClT6vlxzcxsEPoNKJI+IumEyjGwAHgReAyorLpaBjyax48BS3Pl1izgTGBTDovtl3R+zoFcUVWmUtclwJM5z7IeWCBpYq4iWwCsz3NPZd7q9zczsyaoZ8jrFOCRXOE7Frg/Iv5a0mbgIUlXAjuBLwJExLZcgfUS0AtcExGHsq6rgLuB44F1+QK4C7hX0naKnsnSrKtH0q3A5sx3S0T05PEKYK2klcDWrMPMzJrED4c0M7Oa/HBIMzMbVg4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEHUHFEnHSNoq6cf5/dclvS7puXx9tpT3RknbJb0iaWEpfa6kF/Lcqtxbntx//sFM3yhpZqnMMknd+VpWSp+Vebuz7LFHdynMzOxoDKSHch3wclXa7RExJ1+PA0iaTbEn/NnAIuA7ko7J/KuB5cCZ+VqU6VcC+yLiDOB24LasaxJwM3AeMA+4WdLELHNbvv+ZwL6sw8zMmqSugCJpBvA54M/ryL4YWBsRByPiVWA7ME/SVGBCRGyIYiP7e4AlpTJr8vhh4KLsvSwEOiOiJyL2AZ3Aojw3P/OSZSt1mZlZE9TbQ/kz4Hrgg6r0ayX9QtIPSj2H6cBrpTy7Mm16HlenH1YmInqBd4DJNeqaDLydeavrMjOzJug3oEj6PLAnIrZUnVoNfByYA+wGvlUp0kc1USN9MGVq1XUYScsldUnq2rt3b19ZzMysAerpoVwAfEHSDmAtMF/SfRHxZkQciogPgO9TzHFA0Vs4tVR+BvBGps/oI/2wMpLGAicCPTXqegs4KfNW13WYiLgzIjoiomPKlCl1/LhmZjYY/QaUiLgxImZExEyKyfYnI+L3ck6k4mLgxTx+DFiaK7dmUUy+b4qI3cB+SefnHMgVwKOlMpUVXJfkewSwHlggaWIOqS0A1ue5pzIvWbZSl5mZNcHY/rMc0TclzaEYatoBfBkgIrZJegh4CegFromIQ1nmKuBu4HhgXb4A7gLulbSdomeyNOvqkXQrsDnz3RIRPXm8AlgraSWwNeswM7MmUfFhf3To6OiIrq6uZjfDzKytSNoSER395fOd8mZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGY2Kux59wC/+70N7Nl/oNlNGbEcUMxsVFj1s2427+hh1RPdzW7KiHU0z/IyM2t5Z920joO9H27ldN/Gndy3cSfHjR3DKys/08SWjTzuoZjZiPbM9RfyhTnTGD+u+HM3ftwYFs+ZxjMrLmxyy0YeBxQza5rhmNc4ecJ4TjhuLAd7P+C4sWM42PsBJxw3lpNPGD9k7zlaOaCYWdMM17zGW+8d5LLzTueRqy/gsvNOZ+97B4f0/UYrP77ezIZd9bxGhec1WpMfX29mLcvzGiOTA4qZDTvPa4xMdQcUScdI2irpx/n9JEmdkrrz68RS3hslbZf0iqSFpfS5kl7Ic6tyb3ly//kHM32jpJmlMsvyPbolLSulz8q83Vn22KO7FGY2nDyvMfLUPYci6b8BHcCEiPi8pG8CPRHxDUk3ABMjYoWk2cADwDxgGvAE8ImIOCRpE3Ad8CzwOLAqItZJuhr4jYj4iqSlwMUR8SVJk4CufN8AtgBzI2Jf7lv/o4hYK+m7wPMRsbrWz+A5FDOzgWvoHIqkGcDngD8vJS8G1uTxGmBJKX1tRByMiFeB7cA8SVMpgtGGKKLYPVVlKnU9DFyUvZeFQGdE9ETEPqATWJTn5mfe6vc3M7MmqHfI68+A64HysoxTImI3QH49OdOnA6+V8u3KtOl5XJ1+WJmI6AXeASbXqGsy8Hbmra7rMJKWS+qS1LV37946f1wzMxuofgOKpM8DeyJiS511qo+0qJE+mDK16jo8MeLOiOiIiI4pU6b0lcXMzBqgnh7KBcAXJO0A1gLzJd0HvJnDWOTXPZl/F3BqqfwM4I1Mn9FH+mFlJI0FTgR6atT1FnBS5q2uy6wmP3XWbGj0G1Ai4saImBERM4GlwJMR8XvAY0Bl1dUy4NE8fgxYmiu3ZgFnAptyWGy/pPNzDuSKqjKVui7J9whgPbBA0sRcRbYAWJ/nnsq81e9vVpOfOms2NI7macPfAB6SdCWwE/giQERsyxVYLwG9wDURcSjLXAXcDRwPrMsXwF3AvZK2U/RMlmZdPZJuBTZnvlsioiePVwBrJa0EtmYdZkfkp86aDS0/esVGjT3vHmDl4y/z023/xIH3P2D8uDEsPPvX+drnPukb6sxq8KNXzKoM5O7sdpxnacc228jigGKjSr13Z7fjPEs7ttlGFg95mZW041Nw27HN1l485GU2CO34FNx2bLONTA4oZiXt+BTcdmyzjUxHs2zYbESqzLNcOu807t+0k71tMMndjm22kcdzKGZmVpPnUMzMbFg5oJiZWUM4oJiZDYBvID0yBxQzswHwDaRH5lVe1lb2vHuAax/Yyh2XnutlsTas/HDR/rmHYm3Fnw6tWXwDaf/cQ7G2MNo/Hbpn1ny+gbR/7qFYWxjtnw7dM2sN9T5cdLRyD8VaRq1P4aP10+Fo75m1mu9d/uG9fSuXnNPElrSmfnsoksZL2iTpeUnbJP1xpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEtiA9WoZZT9fQofjZ8OR3vPzNpLPT2Ug8D8iHhP0jjgbyVVtu69PSL+tJxZ0myKLXzPBqYBT0j6RG4DvBpYDjwLPA4sotgG+EpgX0ScIWkpcBvwJUmTgJuBDiCALZIei4h9mef2iFgr6btZx+rBXwobrHIgWHnxpwZcvt5P4aPx0+Fo7ZlZe+q3hxKF9/Lbcfmq9QCwxcDaiDgYEa8C24F5kqYCEyJiQxQPELsHWFIqsyaPHwYuyt7LQqAzInoyiHQCi/Lc/MxLlq3UZcPkrJvWMfOGn3Dfxp1EFIFg5g0/4ayb1vVfuMSfwmsbjT0za091zaFIOgbYApwB/O+I2CjpM8C1kq4AuoCv5h/96RQ9kIpdmfZ+Hlenk19fA4iIXknvAJPL6VVlJgNvR0RvH3XZMHnm+guPuEf7QPhTeG2jsWdm7amuVV4RcSgi5gAzKHob51AML30cmAPsBr6V2dVXFTXSB1OmVl2HkbRcUpekrr179/aVxQapkYHAn8JHHj+iZPQZ0CqviHhb0tPAovLciaTvAz/Ob3cBp5aKzQDeyPQZfaSXy+ySNBY4EejJ9E9XlXkaeAs4SdLY7KWU66pu853AnVA8vn4gP6/1r1H7cPhT+MhztHNr1n763Q9F0hTg/QwmxwM/pZgQ3xIRuzPPHwLnRcRSSWcD9wPzKCblfwacGRGHJG0G/iuwkWJS/n9FxOOSrgE+FRFfyUn5346I381J+S3Ab2Zzfg7MjYgeSX8F/LA0Kf+LiPhOrZ/F+6GYDT3vcT/y1LsfSj09lKnAmpxHGQM8FBE/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDW5wgvgKuBu4HiK1V2V2du7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsysyRo1t2btp9+AEhG/AM7tI/3yGmX+BPiTPtK7gF8Zz4iIA8AXj1DXD4Af9JH+/yh6QWbWQrzIYvTynfJm1nDe43508p7yZmZWk/eUNzOzYeWAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYjaEvCeIjSYOKGZDqLwniNlI54dDmg2B6j1B7tu4k/s27vSeIDaiuYdiNgSeuf5CvjBnGuPHFf/Fxo8bw+I503hmxYVNbpnZ0HFAMRsC3hPERiMHFBtRWmkSvLInyCNXX8Bl553O3vcONrtJZkOq34AiabykTZKel7RN0h9n+iRJnZK68+vEUpkbJW2X9IqkhaX0uZJeyHOrJCnTj5P0YKZvlDSzVGZZvke3pGWl9FmZtzvLHtuYS2LtrJUmwb93eQcrl5zD7GkTWLnkHL53eb/bSZi1tX432Mo/+h+JiPckjQP+FrgO+G2gJyK+IekGYGJErJA0G3iAYnveacATwCci4pCkTVn2WeBxYFVErJN0NfAbEfEVSUuBiyPiS5ImAV1AB8Xe9VuAuRGxL/et/1FErJX0XeD5iFhd62fxBlsjV/UkeIUnwc2OXsM22IrCe/ntuHwFsBhYk+lrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIgiit1TVaZS18PARRnIFgKdEdETEfuATmBRnpufeavf30YhT4IPr1YaWrTWUdcciqRjJD0H7KH4A78ROCUidgPk15Mz+3TgtVLxXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdVW3fbmkLklde/furefHtTbkSfDh1UpDi9Y66roPJSIOAXMknQQ8IumcGtnVVxU10gdTplZdhydG3AncCcWQV195bGSoTIJfOu807t+0k73+9Nxwvr/GahnQjY0R8bakp4FFwJuSpkbE7hzO2pPZdgGnlorNAN7I9Bl9pJfL7JI0FjgR6Mn0T1eVeRp4CzhJ0tjspZTrslGqPOm9ckmtzzw2WM9cfyErH3+Zn277Jw68/wHjx41h4dm/ztc+98lmN81aQD2rvKZkzwRJxwP/Gfgl8BhQWXW1DHg0jx8DlubKrVnAmcCmHBbbL+n8nAO5oqpMpa5LgCdznmU9sEDSxFxFtgBYn+eeyrzV729mQ8RDi1ZLPT2UqcAaScdQBKCHIuLHkjYAD0m6EtgJfBEgIrblCqyXgF7gmhwyA7gKuBs4HliXL4C7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsxsiHlo0Y6k32XDI8lQLhve8+4Brn1gK3dceq4/rZk1UDv+32rHNtfSsGXDVh+vejEbGu34f6sd29wI7qEcJd9QZzY02vH/Vju2uR7uoQwT31BnNjTa8f9WO7a5kRxQjtJAVr347mJrB63ye9qOK8rasc2N5IDSAPU+VXa0jqtae2ml39N2fGJzO7a5UTyHMgxG6riqjSz+PbUj8RxKCxnt46rWHvx72npaZfixXg4ow6Cdx1Xb7RfaBq+df09HqlYafqzHgJ7lZYPXrncXl3+hV178qWY3x4ZYu/6ejjTt+hBOz6FYnzyebtY8e949cMSHcDajx+g5FDsqHk83a552HX70kJf1qV1/oc1GinYcfnRAsSNqx19os5GiHff38RzKKDbSnohqZkPDcyjWr3Zbkmhmrc1DXqNQuy5JbDXu4Zkdrp4tgE+V9JSklyVtk3Rdpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEsy8nkFV2O4h2d2uHp6KL3AVyPi55JOALZI6sxzt0fEn5YzS5pNsYXv2cA04AlJn8htgFcDy4FngceBRRTbAF8J7IuIMyQtBW4DviRpEnAz0AFEvvdjEbEv89weEWslfTfrWD34SzF6eAXX0XEPz6xv/fZQImJ3RPw8j/cDLwPTaxRZDKyNiIMR8SqwHZgnaSowISI2RLES4B5gSanMmjx+GLgoey8Lgc6I6Mkg0gksynPzMy9ZtlKX1WE0PxH1aLmHZ9a3Ac2h5FDUucBG4ALgWklXAF0UvZh9FMHm2VKxXZn2fh5Xp5NfXwOIiF5J7wCTy+lVZSYDb0dEbx91WR3acUliq3APz6xvda/ykvRR4IfAH0TEuxTDSx8H5gC7gW9VsvZRPGqkD6ZMrbqq271cUpekrr179/aVxWzA3MMz+1V19VAkjaMIJn8ZET8CiIg3S+e/D/w4v90FnFoqPgN4I9Nn9JFeLrNL0ljgRKAn0z9dVeZp4C3gJEljs5dSruswEXEncCcU96HU8/Oa9cc9PLNfVc8qLwF3AS9HxLdL6VNL2S4GXszjx4CluXJrFnAmsCkidgP7JZ2fdV4BPFoqU1nBdQnwZM6zrAcWSJooaSKwAFif557KvGTZSl1mZtYE9fRQLgAuB16Q9Fym/RHwXyTNoRhq2gF8GSAitkl6CHiJYoXYNbnCC+Aq4G7geIrVXesy/S7gXknbKXomS7OuHkm3Apsz3y0R0ZPHK4C1klYCW7MOMzNrEj96xczMavKjV8zMbFg5oJiZNcFI3F7bAcVsFBmJf8Ta1Uh8dI8fDmk2ipT/iK28+FPNbs6oNJIf3eNJebNRoPqPWMVI+CPWbhq5X/xwPfHak/Jm9q/8/LHW0chH97TasJmHvMyabDg+Zfr5Y63laLfXbtVhMwcUsyYbrnmNo/0jZo1ztI/ueeb6C484bNZMDihmTTLcnzL9/LGRo1V7nJ5DMWsSz2s0zmhcDt2KT7x2D8WsSVr1U2Y7Go3LoVuxx+mAYtZEntc4Oq06OT1a+T4UM2tbjbynw47M96GY2YjnYcPW4iEvM2trHjZsHR7yMjMbwRpx46yHvMzMbFgfz1LPnvKnSnpK0suStkm6LtMnSeqU1J1fJ5bK3Chpu6RXJC0spc+V9EKeW5V7y5P7zz+Y6RslzSyVWZbv0S1pWSl9VubtzrLHNuaSmJm1v7NuWsfMG37CfRt3ElGsgJt5w08466Z1/RcepHp6KL3AVyPik8D5wDWSZgM3AD+LiDOBn+X35LmlwNnAIuA7ko7JulYDy4Ez87Uo068E9kXEGcDtwG1Z1yTgZuA8YB5wcylw3Qbcnu+/L+swMzOac+NsvwElInZHxM/zeD/wMjAdWAysyWxrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIhi4uaeqjKVuh4GLsrey0KgMyJ6ImIf0AksynPzM2/1+5uZjXrNWAE3oFVeORR1LrAROCUidkMRdCSdnNmmA8+Wiu3KtPfzuDq9Uua1rKtX0jvA5HJ6VZnJwNsR0dtHXWZmxvCvgKs7oEj6KPBD4A8i4t2c/ugzax9pUSN9MGVq1XV4Y6TlFMNsnHbaaX1lMTMbkYb78Sx1rfKSNI4imPxlRPwok9/MYSzy655M3wWcWio+A3gj02f0kX5YGUljgROBnhp1vQWclHmr6zpMRNwZER0R0TFlypR6flwzMxuEelZ5CbgLeDkivl069RhQWXW1DHi0lL40V27Noph835TDY/slnZ91XlFVplLXJcCTOc+yHlggaWJOxi8A1ue5pzJv9fubmVkT1DPkdQFwOfCCpOcy7Y+AbwAPSboS2Al8ESAitkl6CHiJYoXYNRFxKMtdBdwNHA+syxcUAeteSdspeiZLs64eSbcCmzPfLRHRk8crgLWSVgJbsw4zM2sS3ylvZmY1+U55aymjcQMks9HGAcWGxXA+/sHMmsNPG7Yh5Q2QzEYP91BsSHnfdLPRwwHFhpQ3QDIbPTzkZUPOGyCZjQ5eNmxmZjV52bCZmQ0rBxQzM2sIB5QRyjcSmtlwc0AZoXwjoZkNN6/yGmF8I6GZNYt7KCOMbyQ0GzwPFR8dB5QRxjcSmg2eh4qPjoe8RiDfSGg2MB4qbgzf2Ghmo96edw+w8vGX+em2f+LA+x8wftwYFp7963ztc590754G3tgo6QeS9kh6sZT2dUmvS3ouX58tnbtR0nZJr0haWEqfK+mFPLcqtwEmtwp+MNM3SppZKrNMUne+lpXSZ2Xe7ix7bD0XxcysLx4qbox65lDuBhb1kX57RMzJ1+MAkmZTbN97dpb5jqRjMv9qYDnFHvNnluq8EtgXEWcAtwO3ZV2TgJuB84B5wM25rzyZ5/aIOBPYl3WYmQ1aZaj4kasv4LLzTmfveweb3aS20+8cSkT8TbnX0I/FwNqIOAi8mnvEz5O0A5gQERsAJN0DLKHYU34x8PUs/zBwR/ZeFgKdlT3kJXUCiyStBeYDl2aZNVl+dZ1tNDP7Fd+7/MMRnZVLzmliS9rX0azyulbSL3JIrNJzmA68VsqzK9Om53F1+mFlIqIXeAeYXKOuycDbmbe6LjMza5LBBpTVwMeBOcBu4FuZrj7yRo30wZSpVdevkLRcUpekrr179x4pm5mZHaVBBZSIeDMiDkXEB8D3KeY4oOgtnFrKOgN4I9Nn9JF+WBlJY4ETgZ4adb0FnJR5q+vqq613RkRHRHRMmTJloD+qmZnVaVABRdLU0rcXA5UVYI8BS3Pl1iyKyfdNEbEb2C/p/JwfuQJ4tFSmsoLrEuDJKNYyrwcWSJqYQ2oLgPV57qnMS5at1GVmZk3S76S8pAeATwMfk7SLYuXVpyXNoRhq2gF8GSAitkl6CHgJ6AWuiYhDWdVVFCvGjqeYjF+X6XcB9+YEfg/FKjEiokfSrcDmzHdLZYIeWAGslbQS2Jp1mJlZE/nGRjMzq6neGxtHVUCRtBf4xz5OfYxibqaduM1Dr93aC27zcGm3Nh9te0+PiH4noUdVQDkSSV31RN9W4jYPvXZrL7jNw6Xd2jxc7fXThs3MrCEcUMzMrCEcUAp3NrsBg+A2D712ay+4zcOl3do8LO31HIqZmTWEeyhmZtYQIzKgHGEPl38raUPuyfJ/JE3I9HGS1mT6y5JuLJV5Ovd1qez7cnKLtPlYSX+R6c9L+nSpTJ/7zrR4m4flOks6VdJT+e+8TdJ1mT5JUmfur9NZetjpgPf3afE2t+R1ljQ5878n6Y6qulryOvfT5iG/zoNo729J2pLXcouk+aW6GneNI2LEvYD/CPwm8GIpbTPwn/L494Fb8/hSikfuA/waxZ3/M/P7p4GOFmzzNcBf5PHJwBZgTH6ruEO0AAADpklEQVS/Cfh3FA/RXAd8pg3aPCzXGZgK/GYenwD8PTAb+CZwQ6bfANyWx7OB54HjgFnAPwDHDOd1bnCbW/U6fwT498BXgDuq6mrV61yrzUN+nQfR3nOBaXl8DvD6UFzjEdlDiYi/oXiMS9lZwN/kcSfwO5XswEdUPGzyeOBfgHeHo51lA2zzbOBnWW4P8DbQoeIZaxMiYkMUvymVfWdats1D1ba+RMTuiPh5Hu8HXqbY+mAxxb465NfKNfvX/X0i4lWgsr/PsF3nRrV5KNrWqDZHxD9HxN8CB8r1tPJ1PlKbh8sg2rs1IioP0d0GjFfxzMWGXuMRGVCO4EXgC3n8RT58kvHDwD9TPIZ/J/Cn8eEzwwD+Irut/2Moh4+O4Ehtfh5YLGmsiodwzs1ztfadGS4DbXPFsF5nFZvGnQtsBE6J4gGm5NfKEMVg9vcZMkfZ5opWvM5H0srXuT/Ddp0H0d7fAbZGsRFiQ6/xaAoovw9cI2kLRRfxXzJ9HnAImEYxRPBVSf8mz10WEZ8C/kO+Lh/eJh+xzT+g+IfvAv4M+L8UD+Mc0F4xQ2SgbYZhvs6SPgr8EPiDiKjVG23InjyN0IA2Q+te5yNW0Udaq1znWobtOg+0vZLOpthC/cuVpD6yDfoaj5qAEhG/jIgFETEXeIBibBmKOZS/joj3cyjm78ihmIh4Pb/uB+5n+IcO+mxzRPRGxB9GxJyIWAycBHRTe9+ZVm3zsF5nSeMo/gP+ZUT8KJPfzK5/ZZhlT6YPZn+fVm1zK1/nI2nl63xEw3WdB9peSTOAR4ArIqLy96+h13jUBJTKSgtJY4CbgO/mqZ3AfBU+ApwP/DKHZj6WZcYBn+fDfV+a2mZJv5ZtRdJvAb0R8VLU3nemJds8nNc5r8ldwMsR8e3SqfKePOX9dQazv09LtrnFr3OfWvw6H6meYbnOA22vpJOAnwA3RsTfVTI3/BoPdja/lV8Un4x3A+9TROArgesoVkL8PfANPryp86PAX1FMVL0E/Pf4cBXHFuAXee5/kqtlWqDNM4FXKCbinqB4Emilng6KX+B/AO6olGnVNg/ndaZYlRP5Xs/l67PAZIoFA935dVKpzNfyWr5CafXLcF3nRrW5Da7zDooFHu/l79LsNrjOv9Lm4brOA20vxYe7fy7lfQ44udHX2HfKm5lZQ4yaIS8zMxtaDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQ/x/rFzgxQKYDrwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2383,55 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2042389\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": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2446,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEKCAYAAAAGvn7fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGbhJREFUeJzt3X2UJXV95/H3h5kBhmkYjAONDkr7QAjIqDgXXWQ13WhcdNCcGE5QQcVIGo0PRCdnM8v6sLrLOj5MsphgkklUiAodw8PZyBjUE2hQNEgPqA2OEBdmlSEMAjLSMAuMfPePX7XctP1wq27dvsXPz+ucPn3vrbpVn/rdut9b9auqexURmJlZPvbqdwAzM6uXC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLzNJeTnzVqlUxNDQ067AHH3yQFStW9HL2lTU5GzQ7n7NV1+R8zlZd2Xxbt269JyIO6mqmEdGzv7Vr18ZcrrrqqjmH9VuTs0U0O5+zVdfkfM5WXdl8wER0WXvdFWNmlhkXdjOzzLiwm5llxoXdzCwzLuxmZpkpVdglvUfSzZJuknSRpH17FczMzKrpuLBLWg28G2hFxNHAEuB1vQpmZmbVlO2KWQosl7QU2A+4s/5IZmbWDUWJ3zyVdBZwDrAb+GpEnDrLOKPAKMDg4ODasbGxWac1NTXFwMBAlcw91+RsUH++yR27apvW4HLYubvz8desXlnbvBfyq/a61snZqiubb2RkZGtEtLqZZ8eFXdKTgEuAU4D7gX8ALo6Iz8/1nFarFRMTE7MOGx8fZ3h4uGzeRdHkbFB/vqENW2qb1vo1e9g02fk3VWzfuK62eS/kV+11rZOzVVc2n6SuC3uZrpiXA7dHxE8i4lHgUuDF3czczMzqV6aw/wj4D5L2kyTgZcC23sQyM7OqOi7sEXEdcDFwAzBZPHdzj3KZmVlFpb62NyI+CHywR1nMzKwGvvLUzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMdF3ZJR0j6TtvfzyT9US/DmZlZeR3/NF5E3AI8H0DSEmAHcFmPcpmZWUVVu2JeBvyfiPi/dYYxM7PuKSLKP0n6DHBDRPzFLMNGgVGAwcHBtWNjY7NOY2pqioGBgdLzXgxNzgb155vcsau2aQ0uh527Ox9/zeqVtc17Ie3tVucylzHf8jZ5vXO26srmGxkZ2RoRrW7mWbqwS9obuBN4TkTsnG/cVqsVExMTsw4bHx9neHi41LwXS5OzQf35hjZsqW1a69fsYdNkxz18bN+4rrZ5L6S93epc5jLmW94mr3fOVl3ZfJK6LuxVumJeSdpan7eom5lZf1Qp7K8HLqo7iJmZ1aNUYZe0H/BbwKW9iWNmZt3qvDMUiIiHgCf3KIuZmdXAV56amWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZabsT+MdKOliST+QtE3Scb0KZmZm1ZT6aTzgXOCKiDhZ0t7Afj3IZGZmXei4sEs6AHgpcDpARDwCPNKbWGZmVpUiorMRpecDm4HvA88DtgJnRcSDM8YbBUYBBgcH146Njc06vampKQYGBqon75HJHbsYXA47dy/+vNesXtnReHW33eSOXbVNq2zbdbrMdWhvtzqXuYz5lrep7wlwtm6UzTcyMrI1IlrdzLNMYW8B/wIcHxHXSToX+FlEvH+u57RarZiYmJh12Pj4OMPDw+UT99jQhi2sX7OHTZNle6m6t33juo7Gq7vthjZsqW1aZduu02WuQ3u71bnMZcy3vE19T4CzdaNsPkldF/YyB0/vAO6IiOuK+xcDL+hm5mZmVr+OC3tE3AX8WNIRxUMvI3XLmJlZg5Ttb3gX8IXijJjbgLfUH8nMzLpRqrBHxHeArvp+zMyst3zlqZlZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmSv2CkqTtwAPAz4E93f6StpmZ1a/sb54CjETEPbUnMTOzWrgrxswsM4qIzkeWbgd+CgTw1xGxeZZxRoFRgMHBwbVjY2OzTmtqaoqBgYEqmXtqcscuBpfDzt2LP+81q1d2NF7dbTe5Y1dt0yrbdp0ucx3a263OZS5jvuVt6nsCnK0bZfONjIxs7babu2xhf2pE3CnpYOBrwLsi4pq5xm+1WjExMTHrsPHxcYaHh0vG7b2hDVtYv2YPmyar9FJ1Z/vGdR2NV3fbDW3YUtu0yrZdp8tch/Z2q3OZy5hveZv6ngBn60bZfJK6LuylumIi4s7i/93AZcALu5m5mZnVr+PCLmmFpP2nbwOvAG7qVTAzM6umTH/DIHCZpOnnXRgRV/QklZmZVdZxYY+I24Dn9TCLmZnVwKc7mpllxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWWmdGGXtETSjZIu70UgMzPrTpUt9rOAbXUHMTOzepQq7JIOBdYBf9ubOGZm1i1FROcjSxcDHwH2B/44Ik6aZZxRYBRgcHBw7djY2KzTmpqaYmBgoErmnprcsYvB5bBzd7+TzK3J+cpmW7N6Ze/CzNC+zk3u2LVo82033/I29T0BztaNsvlGRka2RkSrm3ku7XRESScBd0fEVknDc40XEZuBzQCtViuGh2cfdXx8nLmG9dPpG7awfs0eNk123DSLrsn5ymbbfupw78LM0L7Onb5hy6LNt918y9vU9wQ4Wzf6ka9MV8zxwGskbQfGgBMkfb4nqczMrLKOC3tE/JeIODQihoDXAVdGxGk9S2ZmZpX4PHYzs8xU6qiNiHFgvNYkZmZWC2+xm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMdFzYJe0r6duSvivpZkkf6mUwMzOrpswvKD0MnBARU5KWAd+Q9E8R8S89ymZmZhV0XNgjIoCp4u6y4i96EcrMzKpTqtcdjiwtAbYCzwbOi4g/mWWcUWAUYHBwcO3Y2Nis05qammJgYGDOeU3u2NVxrroNLoedu/s2+wU1OV/ZbGtWr+xdmBna17l+rV/zLe9C74l+eqJma8LrXLbtRkZGtkZEq5v5lyrsv3iSdCBwGfCuiLhprvFarVZMTEzMOmx8fJzh4eE55zG0YUvpXHVZv2YPmyYr/c73omhyvrLZtm9c18M0/177Otev9Wu+5V3oPdFPT9RsTXidy7adpK4Le6WzYiLifmAcOLGbmZuZWf3KnBVzULGljqTlwMuBH/QqmJmZVVNmf/4pwAVFP/tewBcj4vLexDIzs6rKnBXzPeCYHmYxM7Ma+MpTM7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8yU+c3Tp0m6StI2STdLOquXwczMrJoyv3m6B1gfETdI2h/YKulrEfH9HmUzM7MKOt5ij4h/i4gbitsPANuA1b0KZmZm1Sgiyj9JGgKuAY6OiJ/NGDYKjAIMDg6uHRsbm3UaU1NTDAwMzDmPyR27Sueqy+By2Lm7b7NfUJPzOVt1vc63ZvXKys9d6P3aT/Nl61cdaW/rsm03MjKyNSJa3cy/dGGXNABcDZwTEZfON26r1YqJiYlZh42PjzM8PDznc4c2bCmVq07r1+xh02SZXqrF1eR8zlZdr/Nt37iu8nMXer/203zZ+lVH2tu6bNtJ6rqwlzorRtIy4BLgCwsVdTMz648yZ8UI+DSwLSL+tHeRzMysG2W22I8H3gicIOk7xd+repTLzMwq6rhDLyK+AaiHWczMrAa+8tTMLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsM2V+8/Qzku6WdFMvA5mZWXfKbLGfD5zYoxxmZlaTjgt7RFwD3NfDLGZmVgNFROcjS0PA5RFx9DzjjAKjAIODg2vHxsZmHW9qaoqBgYE55zW5Y1fHueo2uBx27u7b7BfU5HzOVl2v861ZvbLycxd6v/bTfNn6VUfa27ps242MjGyNiFY386+9sLdrtVoxMTEx67Dx8XGGh4fnfO7Qhi0d56rb+jV72DS5tG/zX0iT8zlbdb3Ot33jusrPXej92k/zZetXHWlv67JtJ6nrwu6zYszMMuPCbmaWmTKnO14EfAs4QtIdkt7au1hmZlZVxx16EfH6XgYxM7N6uCvGzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDOlCrukEyXdIumHkjb0KpSZmVVX5jdPlwDnAa8EjgJeL+moXgUzM7NqymyxvxD4YUTcFhGPAGPAb/cmlpmZVaWI6GxE6WTgxIg4o7j/RuBFEfHOGeONAqPF3SOAW+aY5CrgniqhF0GTs0Gz8zlbdU3O52zVlc13WEQc1M0Ml5YYV7M89kufChGxGdi84MSkiYholZj/omlyNmh2Pmerrsn5nK26fuQr0xVzB/C0tvuHAnfWG8fMzLpVprBfDxwu6RmS9gZeB/xjb2KZmVlVHXfFRMQeSe8EvgIsAT4TETd3Me8Fu2v6qMnZoNn5nK26JudztuoWPV/HB0/NzOyJwVeempllxoXdzCwzLuxmZpl5QhZ2Saslre53jtlIeqak90g6od9ZZmpyNmh2Pmerrsn5mpwNqud7QhV2SUOSrgauAD4u6SX9ztRO0n8Evkb6Lp23SXp7nyP9QpOzQbPzOVt1Tc7X5GzQZb6IaPQfsG/b7dcCnyhuvxn4B2BNcV99yHYC8Izp+QMfAE4r7r8I+BIw3I98Tc7W9HzOlme+JmerO18jt9glHSDpryTdCnxC0mHFoN8BflTcHgN+CJwx/bRFzHeUpO8B/w34rKQTIrX2UcAhABFxHfBN4C2Lma/J2Zqez9nyzNfkbL3K18jCDpwI7EtasEeAD0haTtoteTVARDwMXAy8pLj/WK/CSDpU0gFtD50CXBIRLyV9wLxB0uHAhdP5CpcBR0vap1f5mpyt6fmcLc98Tc62WPn6VtiVLJX0Vklfl3SWpGcVg58NPBIRe4A/A34KnAZ8FXiKpF8rxrsV+LGk43qU8UhJXwa+AXxY0vTXFP8/YL/i9heBu4B1pE/UJ7ftYdxH+nbL5/0qZWt6PmfLM1+Tsy12vr4V9mJX4zeBNwEfA/YB/qYYfBdwd/HJ9GPSwjyL1ADf5/GvBV4G3Fs8XgtJK9ruPh+4IyKGgCuBTxSP3wc8LGn/iLgP+FfgqUWObwLvLcbbG/g5sD33bE3P52x55mtytn7mW7TCLuk4SR+VdHpxX8CRwBUR8aWI+BhwmKQXAztIn2BHFk/fBgwUj/0F8CpJryZ9KAwC3+0y25MknS/pemCjpIOKfGuAayUpIv4RuF/SOtKewv7FcIr7BwOPkfYwDpb0N8BFwJ6IuDvHbE3P52zVNTlfk7M1Jd+iFHZJzwH+EngA+D1J7y3mvRp4oFhogPOBN5AK9R7gxcXjN5COGD8UEdcAG4DTgeOB/x4Rj7VNo4qXFvN7FemgxNnAAaQvOzuk2LsAuKDI9+1iWV4JEBHfKqaxNCK2AWcCNwP/MyLeQneanK3p+Zwtz3xNztaMfHOdLlP1j7RlfQZpt2Np8difAmcVt1vAJ4GTgZcDX2l77tNIuyqQCvmNpF9hOgb438BT2sYtfTpS0bBnAleTunNWFY9/EXh3cfsZwMZi+LGk/rAlbcv2k2I6q0l7Eu8EPgt8CljRRbs1NlvT8zmbX1e33b//q3WLXdLzSQc4fxv4IPC+YtAO0m+mQvrkuRb4XeCfgUMkPVfSskj96TskvSQiriR93eVHgUuBiyLi36bnFUXLlHQS8BrgQ8BxpL59SGfbTO8d/Bj4OvDKiLie9Ik7UsxzCrgOODYidgBvJHUF3QW8LyIeLBuobU/j1U3LNoPbrprGtRu47brJ9kRouzI/jfdLJL0QOBz4akT8hLQ1fmtEnC7pBcA5klrAOPCfJO0XEQ9J+i7we6RzNC8E/gD4pKTdwCRwezGLvwIujIhdJTIpIkLSsaTdnK8DWyKdHvnrwG0RcaWk20lXr74C2Ar8jqRVEXGPpH8FHpT0dODPgdMkHUz61ah7SbtORMQEMFGh3VqkvZoHgI8DdwPP7Hc2t121bE+EdnPb5dd28ym1xa5kmaQ3SbqR1LF/IDBdeH8ObC+2vm8g7VocBzzE46fwADxK2gU5hLRVfhOpf/1q4J6IuAPSVnnFov5S4DOko8ovBz5SjPIYcKuk5RFxe5HvuaQX607S+aTTy7GE1D6XFBlPBdYCm6PiOa6SVkr6bDHN24FzI+JuSXuRPsn7mW1J0Xa/SdoVbEzbFevdgKTzaVjbFfMMScM0c53bR9KKhrbdAQ1vuwFJ+0q6gIa13YI66a8BVgAvLm4fWAT75CzjnUW6DHZ1cf9kUn/6YaSvALi6eHxfUjfMqrbnHgPs3UmeGfPcD3gbj2/5LwP+CHhHMfxJwPeK6Z9C6u8aKoadVCzLquL2JLCS1L//5fY8wF5dZLuIdMXYAKlr6cy2caaPQ7wT+B+Lla3tdT2DtLKtJx3gaUrbTWe7tFivDmpY2+0PbCH9khjAe5rQbjPyfRn46+L+x4C39bvtSO+JN5Pe/5c0re3a8l0J/H3xWGPWu07/Ftxil3Q2cBuwRdJgRNxP6he6s+gbf40ev0DoW6QDoNMXGl1LOoj6UERcAPxU0udIB0VvAX7RhxQRN0bEIwvlmZHtEOByYBj4HOkAxWtJewl7iun+lHTg9d2kvq+Defw0ymtI59I/EhGXA58mXc16HumI9aNt+Up9qs7I9nfA24tstwJHSNpYbEX9vtIFV1eQ9mB6nq3It4L05jqBdP3AK0jHPY4lbSn1s+3as20mnS3wWtI1DL/R77YrLCdde/EsSatI6/ySYpp9abdZ8u1NWteeSuriOFrSR/rVdpKWkY6xnQx8PCJ+txh0TNs0+9Z2M/J9LCKmt7gngaP62XaldfAJNkzavfhb4D3FY8eSitYdRfALgU3FsHOAD7c9/3rgmOL2PqRTgI6t41OJtPK+qO3+6aQtkzcD3257/KnAncXtd5Au231S8fwvAU9vG3dVj7K9iXSk+9eBvy/+Xg/8L9K5/IuWrW16B7bd/s+kN9Op/W67WbL9MemUsWc2qO3eTOprfT/wVtKBtOv73W6z5HsfaY9nVRPajrQHduqMx04BrmtC282R7+lFhr6vdx0vRwcLOn1qzinAeHF7GWlramVx/zDS1vqxpF3Ai0lbWv9E+qTapyfhUx+X4Be/3foCHu/uuZd0zuj0uF+jKLSk3aevFuP8ySJlOwb4xvSK2zbeMtLB5ROK++f0OtuMnAeQjm/sBD5c3L8XGOxX282S7a5ivisouvn61XZtr+dbSN1srwW+UDx2T7/bbY58Y8Vj7acL92W9I3VR3ApsKub/gaJ+3Acc3IB1rj3fVaQv5jq03+td6eUoscBPJl0o9Jzi/tIZw88HTp5egUhdD2fSo6I+x8p8AY+fL/854KPF7V8j7XE8ve2FOZq2rwRepGzvaH+suH1I0XbPXexsbRn+kHS+7WZSv/Y3izec+tl2M7KdRzqt7NlNaDvSV0YvIfWhXk3aMr4JeH+/17lZ8v0z6QyzFzSk7b5C2gN7Gmkr+CzShmFT1rn2fF8gXfp/eBPartO/6aLTEUmfAn4WERuK+3uRzrt8B/Ac4JQo2U9eF0mHkvq03hURtyp9odhokWs18J2o56qybrK9PSJuKx47htQtta7I9of9yNZO6TqEM0lvsiNJK+uh9LHt2rIdTXqz/TnpLKuT6FPbSRogdXPsQ2qn3yBdeHI2aUv5cPrYbrPkO5x0fOK3SMe8XkZqv76sdypOey5uP4/0Pr2WdEl939e5GfmOJl3pfi7pm2b7tt6VUfY89s3AucVBhiNJK/HxpBfl7H4V9cIxFOfASzqD1P9/NqkL6QeRTr/sd7YfFdluJ60ce0hb8Tf2MVu7e0kHAd8XEX8n6TTg5obku5/UT3wT6XVdRv/abg/p7IlHSVvqPyet/5PAexvQbnPle1jSa0gFv2/r3XTRLNxPOu70/oi4sAFtNzPfA6SN123Af6W/613Hym6xv450oPRh0jeOXRkRt/QoWymSriUdXNtOOof0QxHxvb6GKszIdhewoUHttpK0BfcG0vffbwbOi4hH533iIpgl26cjYlN/U/2y4sKT6b7su/qdZ6Yi38nAZyOdddLvPPuQfnPhjaQ96r8EPhXpa7r7bpZ8myPiz/qbqpyOC7uk55LO57yYdLCotq/K7VaxB/FB0pbw5yNdtdYITc4GIGkpqfvlYVK+Jr2ujc0G6aIu4LEos3W0iJqcT9KZpNNqP9e01xWan28hpbbYzcys+Zr603hmZlaRC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDP/H+KofDj+oV4qAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -341,9 +2479,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +2500,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4, -- 2.18.1