diff --git a/module3/exo1/analyse-syndrome-grippal .ipynb b/module3/exo1/analyse-syndrome-grippal .ipynb new file mode 100644 index 0000000000000000000000000000000000000000..7fec51502cbc3200b3d0ffc6bbba1fe85e197f3d --- /dev/null +++ b/module3/exo1/analyse-syndrome-grippal .ipynb @@ -0,0 +1,6 @@ +{ + "cells": [], + "metadata": {}, + "nbformat": 4, + "nbformat_minor": 4 +} diff --git a/module3/exo1/influenza-like-illness-analysis.ipynb b/module3/exo1/influenza-like-illness-analysis.ipynb index 87092fc69cd90ff457ea56284b789cb0de199a41..4df4548385c8d13d09f9f821853645cde5cfd974 100644 --- a/module3/exo1/influenza-like-illness-analysis.ipynb +++ b/module3/exo1/influenza-like-illness-analysis.ipynb @@ -9,10 +9,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 1, + "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", @@ -30,10 +28,8 @@ }, { "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\"" @@ -63,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
020204133442127989.040853.05242.062.0FRFrance
120204032139217158.025626.03226.038.0FRFrance
220203931982315907.023739.03024.036.0FRFrance
320203832556221142.029982.03932.046.0FRFrance
420203731848514649.022321.02822.034.0FRFrance
52020363103907646.013134.01612.020.0FRFrance
6202035399186842.012994.01510.020.0FRFrance
7202034360843090.09078.094.014.0FRFrance
8202033361063411.08801.095.013.0FRFrance
9202032359183330.08506.095.013.0FRFrance
10202031343512269.06433.074.010.0FRFrance
11202030381795442.010916.0128.016.0FRFrance
12202029386875860.011514.0139.017.0FRFrance
13202028383405701.010979.0139.017.0FRFrance
14202027340662406.05726.063.09.0FRFrance
15202026340392389.05689.063.09.0FRFrance
16202025328531488.04218.042.06.0FRFrance
17202024330581690.04426.053.07.0FRFrance
18202023341682468.05868.063.09.0FRFrance
19202022335801947.05213.053.07.0FRFrance
20202021361144026.08202.096.012.0FRFrance
21202020393156775.011855.01410.018.0FRFrance
222020193116798722.014636.01814.022.0FRFrance
2320201831639812851.019945.02520.030.0FRFrance
2420201731808214454.021710.02721.033.0FRFrance
2520201632416519893.028437.03731.043.0FRFrance
2620201534104935377.046721.06253.071.0FRFrance
2720201437166664531.078801.010998.0120.0FRFrance
28202013310774299187.0116297.0164151.0177.0FRFrance
29202012310728398610.0115956.0163150.0176.0FRFrance
.................................
184619852132609619621.032571.04735.059.0FRFrance
184719852032789620885.034907.05138.064.0FRFrance
184819851934315432821.053487.07859.097.0FRFrance
184919851834055529935.051175.07455.093.0FRFrance
185019851733405324366.043740.06244.080.0FRFrance
185119851635036236451.064273.09166.0116.0FRFrance
185219851536388145538.082224.011683.0149.0FRFrance
18531985143134545114400.0154690.0244207.0281.0FRFrance
18541985133197206176080.0218332.0357319.0395.0FRFrance
18551985123245240223304.0267176.0445405.0485.0FRFrance
18561985113276205252399.0300011.0501458.0544.0FRFrance
18571985103353231326279.0380183.0640591.0689.0FRFrance
18581985093369895341109.0398681.0670618.0722.0FRFrance
18591985083389886359529.0420243.0707652.0762.0FRFrance
18601985073471852432599.0511105.0855784.0926.0FRFrance
18611985063565825518011.0613639.01026939.01113.0FRFrance
18621985053637302592795.0681809.011551074.01236.0FRFrance
18631985043424937390794.0459080.0770708.0832.0FRFrance
18641985033213901174689.0253113.0388317.0459.0FRFrance
186519850239758680949.0114223.0177147.0207.0FRFrance
186619850138548965918.0105060.0155120.0190.0FRFrance
186719845238483060602.0109058.0154110.0198.0FRFrance
1868198451310172680242.0123210.0185146.0224.0FRFrance
18691984503123680101401.0145959.0225184.0266.0FRFrance
1870198449310107381684.0120462.0184149.0219.0FRFrance
187119844837862060634.096606.0143110.0176.0FRFrance
187219844737202954274.089784.013199.0163.0FRFrance
187319844638733067686.0106974.0159123.0195.0FRFrance
18741984453135223101414.0169032.0246184.0308.0FRFrance
187519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1876 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202041 3 34421 27989.0 40853.0 52 42.0 \n", + "1 202040 3 21392 17158.0 25626.0 32 26.0 \n", + "2 202039 3 19823 15907.0 23739.0 30 24.0 \n", + "3 202038 3 25562 21142.0 29982.0 39 32.0 \n", + "4 202037 3 18485 14649.0 22321.0 28 22.0 \n", + "5 202036 3 10390 7646.0 13134.0 16 12.0 \n", + "6 202035 3 9918 6842.0 12994.0 15 10.0 \n", + "7 202034 3 6084 3090.0 9078.0 9 4.0 \n", + "8 202033 3 6106 3411.0 8801.0 9 5.0 \n", + "9 202032 3 5918 3330.0 8506.0 9 5.0 \n", + "10 202031 3 4351 2269.0 6433.0 7 4.0 \n", + "11 202030 3 8179 5442.0 10916.0 12 8.0 \n", + "12 202029 3 8687 5860.0 11514.0 13 9.0 \n", + "13 202028 3 8340 5701.0 10979.0 13 9.0 \n", + "14 202027 3 4066 2406.0 5726.0 6 3.0 \n", + "15 202026 3 4039 2389.0 5689.0 6 3.0 \n", + "16 202025 3 2853 1488.0 4218.0 4 2.0 \n", + "17 202024 3 3058 1690.0 4426.0 5 3.0 \n", + "18 202023 3 4168 2468.0 5868.0 6 3.0 \n", + "19 202022 3 3580 1947.0 5213.0 5 3.0 \n", + "20 202021 3 6114 4026.0 8202.0 9 6.0 \n", + "21 202020 3 9315 6775.0 11855.0 14 10.0 \n", + "22 202019 3 11679 8722.0 14636.0 18 14.0 \n", + "23 202018 3 16398 12851.0 19945.0 25 20.0 \n", + "24 202017 3 18082 14454.0 21710.0 27 21.0 \n", + "25 202016 3 24165 19893.0 28437.0 37 31.0 \n", + "26 202015 3 41049 35377.0 46721.0 62 53.0 \n", + "27 202014 3 71666 64531.0 78801.0 109 98.0 \n", + "28 202013 3 107742 99187.0 116297.0 164 151.0 \n", + "29 202012 3 107283 98610.0 115956.0 163 150.0 \n", + "... ... ... ... ... ... ... ... \n", + "1846 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1847 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1848 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1849 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1850 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1851 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1852 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1853 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1854 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1855 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1856 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1857 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1858 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1859 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1860 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1861 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1862 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1863 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1864 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1865 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1866 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1867 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1868 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1869 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1870 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1871 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1872 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1873 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1874 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1875 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 62.0 FR France \n", + "1 38.0 FR France \n", + "2 36.0 FR France \n", + "3 46.0 FR France \n", + "4 34.0 FR France \n", + "5 20.0 FR France \n", + "6 20.0 FR France \n", + "7 14.0 FR France \n", + "8 13.0 FR France \n", + "9 13.0 FR France \n", + "10 10.0 FR France \n", + "11 16.0 FR France \n", + "12 17.0 FR France \n", + "13 17.0 FR France \n", + "14 9.0 FR France \n", + "15 9.0 FR France \n", + "16 6.0 FR France \n", + "17 7.0 FR France \n", + "18 9.0 FR France \n", + "19 7.0 FR France \n", + "20 12.0 FR France \n", + "21 18.0 FR France \n", + "22 22.0 FR France \n", + "23 30.0 FR France \n", + "24 33.0 FR France \n", + "25 43.0 FR France \n", + "26 71.0 FR France \n", + "27 120.0 FR France \n", + "28 177.0 FR France \n", + "29 176.0 FR France \n", + "... ... ... ... \n", + "1846 59.0 FR France \n", + "1847 64.0 FR France \n", + "1848 97.0 FR France \n", + "1849 93.0 FR France \n", + "1850 80.0 FR France \n", + "1851 116.0 FR France \n", + "1852 149.0 FR France \n", + "1853 281.0 FR France \n", + "1854 395.0 FR France \n", + "1855 485.0 FR France \n", + "1856 544.0 FR France \n", + "1857 689.0 FR France \n", + "1858 722.0 FR France \n", + "1859 762.0 FR France \n", + "1860 926.0 FR France \n", + "1861 1113.0 FR France \n", + "1862 1236.0 FR France \n", + "1863 832.0 FR France \n", + "1864 459.0 FR France \n", + "1865 207.0 FR France \n", + "1866 190.0 FR France \n", + "1867 198.0 FR France \n", + "1868 224.0 FR France \n", + "1869 266.0 FR France \n", + "1870 219.0 FR France \n", + "1871 176.0 FR France \n", + "1872 163.0 FR France \n", + "1873 195.0 FR France \n", + "1874 308.0 FR France \n", + "1875 213.0 FR France \n", + "\n", + "[1876 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" @@ -80,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
163919891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1639 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1639 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -96,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
020204133442127989.040853.05242.062.0FRFrance
120204032139217158.025626.03226.038.0FRFrance
220203931982315907.023739.03024.036.0FRFrance
320203832556221142.029982.03932.046.0FRFrance
420203731848514649.022321.02822.034.0FRFrance
52020363103907646.013134.01612.020.0FRFrance
6202035399186842.012994.01510.020.0FRFrance
7202034360843090.09078.094.014.0FRFrance
8202033361063411.08801.095.013.0FRFrance
9202032359183330.08506.095.013.0FRFrance
10202031343512269.06433.074.010.0FRFrance
11202030381795442.010916.0128.016.0FRFrance
12202029386875860.011514.0139.017.0FRFrance
13202028383405701.010979.0139.017.0FRFrance
14202027340662406.05726.063.09.0FRFrance
15202026340392389.05689.063.09.0FRFrance
16202025328531488.04218.042.06.0FRFrance
17202024330581690.04426.053.07.0FRFrance
18202023341682468.05868.063.09.0FRFrance
19202022335801947.05213.053.07.0FRFrance
20202021361144026.08202.096.012.0FRFrance
21202020393156775.011855.01410.018.0FRFrance
222020193116798722.014636.01814.022.0FRFrance
2320201831639812851.019945.02520.030.0FRFrance
2420201731808214454.021710.02721.033.0FRFrance
2520201632416519893.028437.03731.043.0FRFrance
2620201534104935377.046721.06253.071.0FRFrance
2720201437166664531.078801.010998.0120.0FRFrance
28202013310774299187.0116297.0164151.0177.0FRFrance
29202012310728398610.0115956.0163150.0176.0FRFrance
.................................
184619852132609619621.032571.04735.059.0FRFrance
184719852032789620885.034907.05138.064.0FRFrance
184819851934315432821.053487.07859.097.0FRFrance
184919851834055529935.051175.07455.093.0FRFrance
185019851733405324366.043740.06244.080.0FRFrance
185119851635036236451.064273.09166.0116.0FRFrance
185219851536388145538.082224.011683.0149.0FRFrance
18531985143134545114400.0154690.0244207.0281.0FRFrance
18541985133197206176080.0218332.0357319.0395.0FRFrance
18551985123245240223304.0267176.0445405.0485.0FRFrance
18561985113276205252399.0300011.0501458.0544.0FRFrance
18571985103353231326279.0380183.0640591.0689.0FRFrance
18581985093369895341109.0398681.0670618.0722.0FRFrance
18591985083389886359529.0420243.0707652.0762.0FRFrance
18601985073471852432599.0511105.0855784.0926.0FRFrance
18611985063565825518011.0613639.01026939.01113.0FRFrance
18621985053637302592795.0681809.011551074.01236.0FRFrance
18631985043424937390794.0459080.0770708.0832.0FRFrance
18641985033213901174689.0253113.0388317.0459.0FRFrance
186519850239758680949.0114223.0177147.0207.0FRFrance
186619850138548965918.0105060.0155120.0190.0FRFrance
186719845238483060602.0109058.0154110.0198.0FRFrance
1868198451310172680242.0123210.0185146.0224.0FRFrance
18691984503123680101401.0145959.0225184.0266.0FRFrance
1870198449310107381684.0120462.0184149.0219.0FRFrance
187119844837862060634.096606.0143110.0176.0FRFrance
187219844737202954274.089784.013199.0163.0FRFrance
187319844638733067686.0106974.0159123.0195.0FRFrance
18741984453135223101414.0169032.0246184.0308.0FRFrance
187519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1875 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202041 3 34421 27989.0 40853.0 52 42.0 \n", + "1 202040 3 21392 17158.0 25626.0 32 26.0 \n", + "2 202039 3 19823 15907.0 23739.0 30 24.0 \n", + "3 202038 3 25562 21142.0 29982.0 39 32.0 \n", + "4 202037 3 18485 14649.0 22321.0 28 22.0 \n", + "5 202036 3 10390 7646.0 13134.0 16 12.0 \n", + "6 202035 3 9918 6842.0 12994.0 15 10.0 \n", + "7 202034 3 6084 3090.0 9078.0 9 4.0 \n", + "8 202033 3 6106 3411.0 8801.0 9 5.0 \n", + "9 202032 3 5918 3330.0 8506.0 9 5.0 \n", + "10 202031 3 4351 2269.0 6433.0 7 4.0 \n", + "11 202030 3 8179 5442.0 10916.0 12 8.0 \n", + "12 202029 3 8687 5860.0 11514.0 13 9.0 \n", + "13 202028 3 8340 5701.0 10979.0 13 9.0 \n", + "14 202027 3 4066 2406.0 5726.0 6 3.0 \n", + "15 202026 3 4039 2389.0 5689.0 6 3.0 \n", + "16 202025 3 2853 1488.0 4218.0 4 2.0 \n", + "17 202024 3 3058 1690.0 4426.0 5 3.0 \n", + "18 202023 3 4168 2468.0 5868.0 6 3.0 \n", + "19 202022 3 3580 1947.0 5213.0 5 3.0 \n", + "20 202021 3 6114 4026.0 8202.0 9 6.0 \n", + "21 202020 3 9315 6775.0 11855.0 14 10.0 \n", + "22 202019 3 11679 8722.0 14636.0 18 14.0 \n", + "23 202018 3 16398 12851.0 19945.0 25 20.0 \n", + "24 202017 3 18082 14454.0 21710.0 27 21.0 \n", + "25 202016 3 24165 19893.0 28437.0 37 31.0 \n", + "26 202015 3 41049 35377.0 46721.0 62 53.0 \n", + "27 202014 3 71666 64531.0 78801.0 109 98.0 \n", + "28 202013 3 107742 99187.0 116297.0 164 151.0 \n", + "29 202012 3 107283 98610.0 115956.0 163 150.0 \n", + "... ... ... ... ... ... ... ... \n", + "1846 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1847 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1848 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1849 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1850 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1851 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1852 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1853 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1854 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1855 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1856 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1857 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1858 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1859 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1860 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1861 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1862 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1863 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1864 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1865 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1866 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1867 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1868 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1869 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1870 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1871 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1872 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1873 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1874 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1875 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 62.0 FR France \n", + "1 38.0 FR France \n", + "2 36.0 FR France \n", + "3 46.0 FR France \n", + "4 34.0 FR France \n", + "5 20.0 FR France \n", + "6 20.0 FR France \n", + "7 14.0 FR France \n", + "8 13.0 FR France \n", + "9 13.0 FR France \n", + "10 10.0 FR France \n", + "11 16.0 FR France \n", + "12 17.0 FR France \n", + "13 17.0 FR France \n", + "14 9.0 FR France \n", + "15 9.0 FR France \n", + "16 6.0 FR France \n", + "17 7.0 FR France \n", + "18 9.0 FR France \n", + "19 7.0 FR France \n", + "20 12.0 FR France \n", + "21 18.0 FR France \n", + "22 22.0 FR France \n", + "23 30.0 FR France \n", + "24 33.0 FR France \n", + "25 43.0 FR France \n", + "26 71.0 FR France \n", + "27 120.0 FR France \n", + "28 177.0 FR France \n", + "29 176.0 FR France \n", + "... ... ... ... \n", + "1846 59.0 FR France \n", + "1847 64.0 FR France \n", + "1848 97.0 FR France \n", + "1849 93.0 FR France \n", + "1850 80.0 FR France \n", + "1851 116.0 FR France \n", + "1852 149.0 FR France \n", + "1853 281.0 FR France \n", + "1854 395.0 FR France \n", + "1855 485.0 FR France \n", + "1856 544.0 FR France \n", + "1857 689.0 FR France \n", + "1858 722.0 FR France \n", + "1859 762.0 FR France \n", + "1860 926.0 FR France \n", + "1861 1113.0 FR France \n", + "1862 1236.0 FR France \n", + "1863 832.0 FR France \n", + "1864 459.0 FR France \n", + "1865 207.0 FR France \n", + "1866 190.0 FR France \n", + "1867 198.0 FR France \n", + "1868 224.0 FR France \n", + "1869 266.0 FR France \n", + "1870 219.0 FR France \n", + "1871 176.0 FR France \n", + "1872 163.0 FR France \n", + "1873 195.0 FR France \n", + "1874 308.0 FR France \n", + "1875 213.0 FR France \n", + "\n", + "[1875 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -123,10 +2117,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 6, + "metadata": {}, "outputs": [], "source": [ "def convert_week(year_and_week_int):\n", @@ -154,10 +2146,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -180,9 +2170,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", @@ -200,9 +2198,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+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXm8XEWZ939Pr3dP7s1GSAIJEJaACxACbogii+II4zKDzkhG8EUdZnTUGQdcXnxhUHRcxhVFRRAXRGUEBcSwKjtJCIRsJGQnl+QmN7n70re73j9O1ek63VXnnO4+vdzk+X4++dzOWarq1KlTTz1LVZEQAgzDMAxTKbF6F4BhGIY5OGCBwjAMw0QCCxSGYRgmEligMAzDMJHAAoVhGIaJBBYoDMMwTCSwQGEYhmEigQUKwzAMEwksUBiGYZhISNS7ALVk+vTpYv78+fUuBsMwzKRixYoVe4UQM4KuO6QEyvz587F8+fJ6F4NhGGZSQUTbwlzHJi+GYRgmEligMAzDMJHAAoVhGIaJBBYoDMMwTCSwQGEYhmEiIVCgENFNRLSHiF7QjnUR0TIi2ij/dmrnriKiTUS0gYjO046fSkSr5blvExHJ42ki+rU8/hQRzdfuWSrz2EhES7XjC+S1G+W9qcqrgmEYhqmEMBrKzQDOLzh2JYAHhBALATwg/w8iWgTgYgAnynu+T0Rxec8NAC4HsFD+U2leBmC/EOIYAN8E8BWZVheAqwGcDmAJgKs1wfUVAN+U+e+XaTAMwzB1JFCgCCH+AqC34PCFAG6Rv28BcJF2/DYhxJgQYguATQCWENFsAB1CiCeEs+fwzwruUWn9FsDZUns5D8AyIUSvEGI/gGUAzpfn3iqvLcz/oGfFtl6s6+6vdzEYhmGKKHdi4ywhRDcACCG6iWimPD4HwJPadTvlsYz8XXhc3bNDpjVBRH0ApunHC+6ZBuCAEGLCkNZBz3tueAIAsPX6C+pcEoZhGC9RO+XJcEz4HC/nHr+0igtEdDkRLSei5T09PbbLGIZhmAopV6DslmYsyL975PGdAOZp180FsEsen2s47rmHiBIApsAxsdnS2gtgqry2MK0ihBA3CiEWCyEWz5gRuBQNwzAMUyblCpS7AKioq6UA7tSOXywjtxbAcb4/Lc1jA0R0hvSBXFJwj0rrvQAelH6W+wCcS0Sd0hl/LoD75LmH5LWF+TMMwzB1ItCHQkS/AnAWgOlEtBNO5NX1AG4nossAbAfwPgAQQqwhotsBrAUwAeAKIURWJvUxOBFjzQDulf8A4CcAbiWiTXA0k4tlWr1EdC2AZ+R11wghVHDAfwK4jYj+C8CzMg2GYRimjpAz4D80WLx4sZjsqw3Pv/JuAOyUZximdhDRCiHE4qDreKY8wzAMEwksUBiGYZhIYIHCMAzDRAILFIZhGCYSWKAwDMMwkcAChWEYhokEFigMwzBMJLBAYRiGYSKBBQrDMAwTCSxQGIZhmEhggcIwDMNEAgsUhmEYJhJYoDAMwzCRwAKFYRiGiQQWKJOIybTVwJW/ex7Hff7e4AsZhjloCNxgi2HK4bZndtS7CAzD1BjWUCYRk0hBYRjmEIQFCsMwDBMJLFAmEaygMAzTyLBAYRiGYSKBBcokYjJFeTEMc+jBAmUSweKEYZhGhgUKwzAMEwksUCYRbPFiGKaRYYHCMAzDRAILlEmEYC8KwzANDAuUSQSbvBiGaWRYoDAMwzCRwAKFYRiGiQQWKAcZuw6M4OjP3oM1u/rqXRSGYQ4xKhIoRPRJIlpDRC8Q0a+IqImIuohoGRFtlH87teuvIqJNRLSBiM7Tjp9KRKvluW8TEcnjaSL6tTz+FBHN1+5ZKvPYSERLK3mOyUIYH8oD6/cgmxP45VPbq18ghmEYjbIFChHNAfBxAIuFECcBiAO4GMCVAB4QQiwE8ID8P4hokTx/IoDzAXyfiOIyuRsAXA5gofx3vjx+GYD9QohjAHwTwFdkWl0ArgZwOoAlAK7WBRfDs+oZhqk9lZq8EgCaiSgBoAXALgAXArhFnr8FwEXy94UAbhNCjAkhtgDYBGAJEc0G0CGEeEI4i1X9rOAeldZvAZwttZfzACwTQvQKIfYDWIa8EDpoCRM2TDUoB8MwjImyBYoQ4mUAXwOwHUA3gD4hxJ8BzBJCdMtrugHMlLfMAaBv47dTHpsjfxce99wjhJgA0Adgmk9aBzWlhA1ziDHDMLWmEpNXJxwNYgGAwwG0EtE/+t1iOCZ8jpd7T2E5Lyei5US0vKenx6d4BwfEKgrDMHWiEpPX2wBsEUL0CCEyAO4A8HoAu6UZC/LvHnn9TgDztPvnwjGR7ZS/C4977pFmtSkAen3SKkIIcaMQYrEQYvGMGTPKfNTGgJUOhmEamUoEynYAZxBRi/RrnA1gHYC7AKioq6UA7pS/7wJwsYzcWgDH+f60NIsNENEZMp1LCu5Rab0XwIPSz3IfgHOJqFNqSufKYwzDMEydSJR7oxDiKSL6LYCVACYAPAvgRgBtAG4nosvgCJ33yevXENHtANbK668QQmRlch8DcDOAZgD3yn8A8BMAtxLRJjiaycUyrV4iuhbAM/K6a4QQveU+y2ShtA22WJ9hGKa2lC1QAEAIcTWc8F2dMTjaiun66wBcZzi+HMBJhuOjkALJcO4mADeVWORJTRgRQRznxTBMneCZ8gzDMEwksECZRHAoMMMwjQwLlIMMFTbMwodhmFrDAmUyEUJIKA8KCxSGYWoNCxSmqpQWmeZlza4+LL3paYxP5CIsEcMw1YIFyiRiMm4BXImmdOXvVuORF3uwrrs/ugIxDFM1WKBMIkpay6tBhE8UpeDlZBhmcsAC5SCj0TrfSkxeDMNMLligTCImY9c8GcvMMEx5sEBhqgorKAxz6MACZRJRivmoUTrySnw56l5eToZhJgcsUCYRk3Etr0oEm7q30fxCDMOYYYHCMAzDRAILlElEaWHDjUGjmN4Yhqk+LFAONhpsLa9GmQ/DMEz1YYEyiQjTOTeauyEKHwrDMJMDFihMValEJrA8YZjJBQuUycQk7GGjmCnPUV4MMzlggTKJKG1H+caQPo1RCoZhagELlIMMarDhPPtBGObQgQXKJGJSds4VOeUn4wMzzKELCxSmqkRhemu02f8Mw5hhgTKJKKlzbpDBPSsZDHPowAJlEhGmc3b3lK9qScLDG2wxzKEDC5SDDO58GYapFyxQJhGNonWUQiWO9TC37h0cw2gmW3YeDMNEBwuUElj/Sj/uXd1d72KEolEipKpt8lr8X/fj0pufiSAXhmEqJVHvAkwmzv+fvwIAtl5/QV3yDyMkGs3kVQu59vhL+6qfCcMwgbCGMoloEKWjJKLYsXEyPjfDHIqwQGGqSwTCgAUKw0wOWKAcpDRKH1zRasON8hAMw4SiIoFCRFOJ6LdEtJ6I1hHR64ioi4iWEdFG+bdTu/4qItpERBuI6Dzt+KlEtFqe+zbJBamIKE1Ev5bHnyKi+do9S2UeG4loaSXPcTDRaLPKoxAKjbLQJcMw/lSqoXwLwJ+EEMcDeA2AdQCuBPCAEGIhgAfk/0FEiwBcDOBEAOcD+D4RxWU6NwC4HMBC+e98efwyAPuFEMcA+CaAr8i0ugBcDeB0AEsAXK0LroOVyThij0IY2J67USLZGIZxKFugEFEHgDMB/AQAhBDjQogDAC4EcIu87BYAF8nfFwK4TQgxJoTYAmATgCVENBtAhxDiCeH0ED8ruEel9VsAZ0vt5TwAy4QQvUKI/QCWIS+EGDSO8KlmOXIN8owMwzhUoqEcBaAHwE+J6Fki+jERtQKYJYToBgD5d6a8fg6AHdr9O+WxOfJ34XHPPUKICQB9AKb5pHVQE2a070ZGVbswIYlix0abUMo1itRkGAZAZQIlAeAUADcIIU4GMARp3rJgMu4Ln+Pl3uPNlOhyIlpORMt7enp8iheeeplawmTbaH1sZTPl/e9lgcIwjUUlAmUngJ1CiKfk/38LR8DslmYsyL97tOvnaffPBbBLHp9rOO65h4gSAKYA6PVJqwghxI1CiMVCiMUzZswo4zGLyTawraXR+thqOuUb7VkZ5lCnbIEihHgFwA4iOk4eOhvAWgB3AVBRV0sB3Cl/3wXgYhm5tQCO8/1paRYbIKIzpH/kkoJ7VFrvBfCg9LPcB+BcIuqUzvhz5bGakK2XhhLRNZMNNnkxzOSg0qVX/hXAL4goBWAzgA/BEVK3E9FlALYDeB8ACCHWENHtcITOBIArhBBqVb+PAbgZQDOAe+U/wHH430pEm+BoJhfLtHqJ6FoAahGna4QQvRU+S2gaW0NprLJFo6GYaeDXwDCHJBUJFCHEKgCLDafOtlx/HYDrDMeXAzjJcHwUUiAZzt0E4KZSyhsV9RIoYYRF3pHdGL1tNGHD5jRYQ2GYxoJnypdBLleffEN1nw3Wx1bS5wfdKur0Hhqd+Vfeja/+aX29i8EcgrBAKYN6+VB0bKP2RptVXlFphH8arKHY+f7DL9W7CMwhCAuUMqifycv823ZNIxCF6Y2d8g6Pb9qLvuFMvYvBMFZYoJTBZOjIGqWE0ZTD5kOJJPFJwch4Fh/48VO49BbeTIxpXFiglEH9oryE4ZfligbpbCOJ8uK1vJCRjrv13f11Lkk49g6O4ZmtNQu8ZBoEFihl0Nhhw/UuQXQEPUqtXsP4RA6fun0VdvQO1yZDA9lsYy2pE8R7bngc7/vBE/UuBlNjWKCUQWP4UCaHU76SLlA9Y72d8k9s3oc7Vr6Mz/7v6prkZyKTdTSUoEduFK1t2776CV+mfrBAKYNGmClvNXm5kVGN0bFU0+RVK4HSlHA+k9FMNuDK6pHJKeHq/8wNIk8OKn63YifmX3k39g2O1bsoDQ8LlDLINbLJq94FKCCK8li1sQgfdmhsAnstHUZT0tm2ZzRTv4kvEyE1lMkQMDLZuH25s7D5i7sH61ySxocFShnUTUMJETbcaEPUai69EuWjvvM7j2Lxf91vPJdONoCGonwoQSavGpSlFBrFBFcJrWlnQZGhsYk6l6TxYYFSBo3glLeuwFvjcgRRiektcKZ8hE+7Ze+Q9ZzaVnk8W0cNRUZ5UcAOz42moTRYccqiOeVoqEPjLFCCYIFSAupjrptTPswGW8L7t95UN2y48rRD5S/rfSJbv0pVeQcJlEZ574oGK05ZJGJOpTda3TYiLFBKICa/5saI8rJd01itvqK1vAICDGr1pGrttkrf+6Y9A76akB9hBwoN9vobTmOqhEYJdGlkWKCUgByoTIqPpFGKGMlHWOcoL1dDqXBV0Ld94y94y9cerqgMtTID3rO6G30jlS/z0ijtsBIClEJGgwVKCZCrodQn/1BbAFe/GCUxWZzyvvnLfCbq6DvLuRqKfxmiKOL2fcP451+sxL/d9mzFaU2GwRcTHSxQSiAewuR156qX8fimvVUvS5Bf4WBQz4OfoUYaiswmW0cfijvJswYTG0cnnGi2HftHKk6r3vLk/rW7ccaXHsDYROURevV+lskAC5QSCGPy+sRtq/CBHz9VlfxDOeWrknP5HFRO+TpqKCrnoBF/FEWM0rRb74HNtXevxSv9o3i5AuFIQZEQjAsLlBKot1Nex+qoDjmSrRWR7NhY59WGXQ2lngLF1TyDLowit+g60Hp/Kk0JJ+R3bIJ3Y6sFLFBKoO5hw2EmNjYYkUR52TSUGo1+1Ui9vqPtcAOFSH0WkWiX5kTW7uqvSURiU4STUifLN1dPWKCUQCzWOBqKjUZr9FEUp+5O+RrnZyxDyLyjKKIaOEWRlulT+dML3XjHt/+KPzzfHUEO/qhlcyqJWGODV3hYoJSAa/JqgF7b2smGDC+tFZWMQhtlzkUDvO7Q7zMKDUV1oJFoEIYkntzs7JOyp3+08vQDOKKrBQCwvY5bDxxKsEApASVQ6rU4ZKjl6xug89OJREOp81L9Qcvo14KwbS4SgRKhE9pUHrWETVpqD9VEaSjVDF9n8rBAKQEV/dIIGoqNRjDP6BwU81Dc/Oof5RXdhXZcDaXypIxpuOKqBvUZiWxkm1doWKCUQL2jvPQReXAn2yASpYoz5Q8pk1fIMkTZNKN4bpOGEqWPJiyNtiTRwQoLlBKo99IrpWTbKN9PNBpKfU1euQYweYV91mjCtKNLyyhQULvFFiPVthrlo2pgWKCUgLIt13PVWUVQKG39S+hQSTmCPuBaayh17U9qqKFE2nEakqrlPEH1zVbySLUUgJMdFiglEJO1VTcNxfof7fBB2OjrvQVwvWd7O2UIeV0EdZKLUICaBFykUWQhiSYEuv7toNFhgVIC8TovDlkKjaKeV3XplcqTluk3hiYURRmiKWu4SZThUjL5UGS0ZA3rNUpBy9hhgVICQfNQqt2J6+kHLb3SKERRnmpHeQV1FI1QpaF9KFUU4OXQKJ1wJdoF1dl3OplggVICbsOyfCW1HXH5H2+Upl9dZ2g0Txm84GL9azNs24qirGGXyg+DXxq1qNW8MKg8rUYbrDUiFQsUIooT0bNE9Ef5/y4iWkZEG+XfTu3aq4hoExFtIKLztOOnEtFqee7bJHViIkoT0a/l8aeIaL52z1KZx0YiWlrpc4QhKGy42uHEYVJvtCZf0VpeAeejqu7AGfnRZFMRYTuzaOaORBfYYSq2GzZcww46SkHL2IlCQ/kEgHXa/68E8IAQYiGAB+T/QUSLAFwM4EQA5wP4PhGpqbI3ALgcwEL573x5/DIA+4UQxwD4JoCvyLS6AFwN4HQASwBcrQuuauHOlLc0zmqPZD0z5QOuaZTBVCUO7SBtKzqTl2ZKNCTaCCPTsCWIpOOUPsJqmc+ohjMFo4jQUqVtBE210alIoBDRXAAXAPixdvhCALfI37cAuEg7fpsQYkwIsQXAJgBLiGg2gA4hxBPC+XJ/VnCPSuu3AM6W2st5AJYJIXqFEPsBLENeCFUNNbKy7YtRz/0yFI0QkeShmk75iD7woFWcG6IfqaFTPq+hRDGqNznlK0629HJE8G02wOfd8FSqofwPgM8A0OOeZgkhugFA/p0pj88BsEO7bqc8Nkf+LjzuuUcIMQGgD8A0n7SqSr1NXnqvErSW1yMv9uDF3QNVLk8wEXX5VUzb2+mZOsBGENLhnfJR+D28fyvBb1RfS0HNPpTaULZAIaJ3AtgjhFgR9hbDMeFzvNx7vJkSXU5Ey4loeU9PT6iC2ojH/BeHrOWikWFy+vAty6tejiCqGXUUlQnCY/IqIf9aUsvl66s8r7EuK04cDFFe2/cN46ZHt9S1DEFUoqG8AcC7iGgrgNsAvJWIfg5gtzRjQf7dI6/fCWCedv9cALvk8bmG4557iCgBYAqAXp+0ihBC3CiEWCyEWDxjxozynlTibrBlaVjVNnmFac8N0Pd5qMiHEnRvZD4ULUlDmo1g6gi/lld0ZqpqRejVcktd13x3EDjlP/DjJ3HNH9dWtLdLtSlboAghrhJCzBVCzIfjbH9QCPGPAO4CoKKulgK4U/6+C8DFMnJrARzn+9PSLDZARGdI/8glBfeotN4r8xAA7gNwLhF1Smf8ufJYTbCHDdcuysseSdsAvZ9GRVFeQU758pMuyCfA5NUAdRq2bUXjQ4kwLZ80alGrKv8oVgivt4biCpL6N0criSqkeT2A24noMgDbAbwPAIQQa4jodgBrAUwAuEIIofbl/BiAmwE0A7hX/gOAnwC4lYg2wdFMLpZp9RLRtQCekdddI4TorcKzeAhqnLVchdi+YGKecgeCo5ks7n6+G+8+ZU7Fo8lqmmBqNbGx3iNToLZRXnkBak9raGwCYxM5dLWmAspTfCy/9Ep55SsF9SyVvENeyys8kQgUIcTDAB6Wv/cBONty3XUArjMcXw7gJMPxUUiBZDh3E4Cbyi1zOahO3Lb0StXnoYQxeUVQhG/e/yJ++MhmdDQncc6iWRWlFc1M+epqhN6wYUP+EfckQoiSBXUtl14Js5bXOd94BLv6RrH1+gv8y2NcHVLlU/0eWuUQTTh1Y0iURggSscEz5UvA1VByZolS031SbKP2CBpb7+C483dorOK0qqqhRJA2EBzlFfVrLa9vq53JS+Xll9SuvnDb95o+lfyIv/rfi+sPilDQ1o0Qgr7esEApAfUixyYsAqWGa3nZr8n/LtdYlUw4zWI8imX6K/GhBJ2vxjwUw/moR9LlpFZbp7zKs3raJVCbRVbVs0QzD6UxevJG3jGWBUoJqNc4ODZhPF9bH0ppx0shFXeaRcYiOEshyg2fwh4vlaCZ8lF3JOWkF/aOamqEUaUVFC0ZJSqLSj7NKCPFoqBRBJsJFigloBrUUJ0ESpgoL4+GUqZDPSU1lEwEQ8hoIoWqa/PSX5vpFUb9/ZYlUGqqoQSbvEpNSyc/cbIWAkU55cvPKxeBUIoSi8W9IWCBUgYjmfr4UMLNQ/EfbYchGXcE0XgUGkolI8PA6KuInPKeiSjVy8fNoozkaho2HKGt3i/IoRYafRTCK6/l1FeiRBlgUC1YoJRAKU75qu+NUsVIj3jAvi+lUFWnfERV4PWhVN8pX1WTV4Rhw9FMBrSbEGsx4o8iryhCj6OkpsE/JcICpQRUZ2PbU17vgKvxzr3ah/Uil7LnkESwD7dbnIoS8b+3Gj4U03uL3odS+j21Xb4+OsxBDupvDTQUN89KTF6N5UNpkGIYYYFSAkGqr246aYh958vEnXgWQVrRpGFOJaoPvNbL10exrpQ17QijmaLRLu0+lFqMtCPRUArSqheuqbCBJQoLlBJQrzHM8vVVaXz+pn7nuJZvuWHDMardPIEw2IoRVX8U4EKJfh5KGa6p0BMbS0/amle1fCj5Tn5y+FAazinfIN+lCRYoJRDkTPRoKFWIxAjTjKJoa2o12Gp1KNGlEc2HFbSWV7TLmZTn/1L3BA0SIilrxSnkMUfNSYFSE6d85QEAtRSAfqjcG2WgZ4IFSgm4GkooH0plL31P/6hvw7Huh1JRrg7RLtddfhru6DLgfBgmsjn8yy9XYu2u/qJznr7GOKIOn48Nvazl+VCcv4F+sQjKGqXPwGjycvOpOPlAotAu8kEKERQoAmoxIbRcWKCUQoDtNyqT15a9Q1jypQfwo79u9mYvzL9t15Rr81KdVpTmk7LurfC8zqaeQfzx+W588teris4F74cS7ai/nLahmlawhlJy0sUECPJS8AtyqMnERvX3IAgbVjRKOUywQCmBvA/FPETwOuXLz+fl/SMAgIc3lL4hWJThxJGYqypPwtoZlPJhjcm5Q+lkcZMPXMsr4mU7ypvYGO6eam3bWy7mlQfs53Qu+t5j+Np9GyrKPwpzVaOYvBQcNnyQEGSPjWoeSkJOLCx0/ofpLMJoMUG4Tvkolk2polAqJW21/lo6Udzkg+osapNXOXUS9paoy1pxWsb0w/k1Vu04gO8+tCmSAlRSL2HNZmd//WH86unt5WcUQJTBEtWCBUoJ5DWUYIFi+1j2D41jza4+33zUVsN+H1yYRlXuTHeK0ikfhY5iM++VkMTYhLP1TjoRLzoXZPKKcjmTstNzfSgBl1VpMmK5+fktvVKbxSEr1y7CRIoJIfBSzxCuumN1WXl8+vbn8PMnt4W6lsOGDxKC4ufDTGx8/4+exAXfftS3cSqBUqShlNiOxsv8YvNRXo2ioZgTKaV8al0ytayMjmctL8OLizpSrSynfA2Xry8liaBn8QsbLifopFSiWL4+H5Vmv6ZSM9TvVu7E53//QqhrG8X0ZoIFSgnkN9gq3+S1/pUBAMCAZYFJAEjGYjI9bwsWnt/BnWzZGop0/QZ9I919I1ixzX+jzEqaflCnUlqUl3NxPObvQzERTShuZeZQN8orwC0fjZkyfBrl1J1qV34j7bBF6B/N4I6VOwPTqbYPJRPFVg8BqHfbKBt9mWCBUgKBGkoIp7xayXdg1C5QlFnDFp7sW0btd7mrBYc1ef3wkc346M9X+penElODm4btfCkdn/M3bmjxOc9AwH5vJXgmT1biQwmcKV962kV5Bdjq9XcaNDL3MyH63RrWrPOZ3zyPT93+HNZ1F4eD63lUZPIKUaZyrQHl0MDyhAVKKaj2FMaHYmvAqj/wG2XYRkSeyXHWj938uxQopFN+NJPFyHi2vExKwFaKUp5P1aUyJ3rPFV9nutePbE7gw7c8g+VbzRpbpT4UdU8jTGycCBDAnrRMF6hOPsQ3AMC3jb3S7+wcOZKxXRMsvIIIsyVyFFs9hIVNXgcZ1pnyJYzcbELJScf569durJ2sblop0/zhruUVcPtETgR+SFH7H8pNe8W2/QDyEWw6UcxD2d0/ivvX7cG//PJZ4/mKfSgh74miqwnqsHRTapAmYdKYwpiQ9Ps27x30zQOw108Uu0+G2VOlpgKlgVUUFigloBqWbR5KKSM3P4FjOxemGUUTNqzu908gG0agRDJTPliA+9E3nMFPHt0CwKah+Ps39Ee01YlrJgzh26pk+fpaRHkFJaELlKBnMZu8nL9+34BHo/NpYmHro7KZ8qpM9muCzNMPrt+NVTsOlF+IkOWoN4l6F2Ay4dpSw6zlFdDx+IcEW50GgdcIy+9yCGq4mWwOOeE8i6mjBqobMx82ad2+bSpnkPZQuLy9IVDMdZaHM0VWT0UpzKecLQwCNRStPoMWujSGDYfIR9d8bAM4c6qF+QfnFUQYjSrIh3LpzcsBAFuvv6Dscig4bPggIQofiunaQlw13a8stuOVdlxa2kHahXoGPy0linBNaxoh09aFSDzA5GVKNIx2kddQzATtuRJE2FvCOP+7+0bw0Po9pRdCUorJyxzkEKw16IMzv28lrE8pisUh/b4nv28gaj8j+1AOElQHK4TZjlnKWl7lmbxCGb3yv8psd2F8OED+eX0FSnlFCJVG+LkZ+euCnPJBUV6296reWZAt3y8NP9QtpYQN2/L5ux8+gQ/d/IzVFl9K2w00eZk0lBCdvH7Kz9+YT7Oy82Hu9VOUMhP2DHb1jZSfuQH2oRwkiIBG7nXK+6flp8aH2YI1jGml3JFMWLvzhHxIP/txGC3psU178bovP2AfydnMeyEfT3+OWIAPJcjkZcsz3zkGDwbKs3g5NwVusBWiI1ZrxQ2Nm0PXg8rn9W8ECBTT/bnidArxCC0/DSWgQqJJpz4DAAAgAElEQVScKe/rlPf5nqPSUPKmwkiSqwosUEpAf4+m0ZUuRGyjr/ykwdJNXqV+E5W3O/8UotJQrrt7Hbr7RvFSjzeax7+LDv9h6UItYfShaJ29Ibcw2kXQjOxKBX3Y1Yb1Z7EJlLa04zrtG8n45hVUlnDXGjQUN5Q33EAkjIZiuySMMAhOO3iAlZFmQJMrMbIIMJk/Lw55kKC3SZPtWJ/ZHhg27DOqt9mlvYeDryn3Gwqj4gP5Z/R1SIYoQxl+Y5l0uAfU69MYNuyJ4jLkE8L/4Zq8LGUIo+UAwD2ru42rD4R9lXraGctKCWpy7ZjlfPAKBZomXpYPRf71W8rEo+37aCjyr3UF8CiivNRfXx+Kc85kUg1asaJUX2cjb7DFUV4loTVyg0DwaChBdmhfDcWsooTplCJZekP5ikJqKL4mrxLKE2Z0H+Z4IXpnYhJegcvXh3DKBzluw2oo//wLZ+WBwmigvMkrvA/FJuiVULW9t3zzsz2r9tu49pl/fYXZG71UH4rV76iepcpRXkoLMQ1YgiLAStU4GlhBYQ2lFLw+lOJGEmZUVUrYcJHJq8Qylktop7z8UCqN8rLO4QjoDPSjfh1GKUurBJ63PKqqgjAaSjU7BL1Z2kbGqtOzvbegQUCwANbSMmp8qqzh3lnhmnY67jJFAQIlipnyfmkooWEyqQaZvEotG4cNHyQI5G2kZh9K+Alf/lFeMr+CNLy2fnsZK8WiIBURxuQV5mMJjlyyHA9pFgmauFjKEuzlRnlV6kPJR3kFXKf9tnVkyixjFSgBxQsSjp5FUo0+qeARv3crCP/yAGaLQdi8AgmRhtL2TEEfQSavUsvWyCYvFiglIIRAQq5WaxoRhXHKhzlva2Bh/CNh29rYRBa/fGq72WQRwmkKhDN5lTKaKrrU3ds++PowS9nYyhPkdA8jDNy1tiw9vvfdlSFQQg4V9PLZVsBVCy7b6ixoEBBUH94IreL7VbH8JpfrZfOLiFSDEduzZHKqLdvzCiKMxp7x0VCCTHalmrwOSqc8Ec0jooeIaB0RrSGiT8jjXUS0jIg2yr+d2j1XEdEmItpAROdpx08lotXy3LdJGoqJKE1Ev5bHnyKi+do9S2UeG4loabnPUQoC+d0UTS+1lHBKv0ZWkckrZMfzrfs34rP/uxr3vNBtyN/710a4iY3B5bFNCsyvKWYTsMGaQ+G5csKCw5irgj7yiic2qntKUFFs7yVSk5dxPpY2k950v6vNhXtnYTpQ2zXKLFttH4rS0k3bIwQty1KqCauB5UlFGsoEgE8LIU4AcAaAK4hoEYArATwghFgI4AH5f8hzFwM4EcD5AL5PRGr7vBsAXA5gofx3vjx+GYD9QohjAHwTwFdkWl0ArgZwOoAlAK7WBVe1EAJIxtVeJYYPSWs4NoERZrXhvMnLe9y7iKFt2G5N1sPewTEAwKDPMvpBH6HqkHxNXiFav1snBfnFAnau1I/6aigeH0rAQCDgvK1OAqOd9N9ldAhhw4b1sgY55W0aTNCIPChsWFco/DQYP0ER5lsC4FaIdX29bLAwuOIXK3HWfz9kPZ/3w9jT+POaVwCYt0cI0lCClq8pxO+buvWJrdjRO1xaghFStkARQnQLIVbK3wMA1gGYA+BCALfIy24BcJH8fSGA24QQY0KILQA2AVhCRLMBdAghnhDO1/qzgntUWr8FcLbUXs4DsEwI0SuE2A9gGfJCqGrkhHAFSuDExpDmIls+JkKZvHxzLb7fZKKxaUiFqA7BbwRWymjKPnM7OG3/pdD164rPe004AffbylKShlK6RAm9KoD22xY2rKwyE1YfipDXmcVXkPbg6dx9NL6wq0WEGZTY0lITDv1ClO9e3Y2t++ydcJjQ4/vXOUvZmJb28QsqAMrRUMzXj4xn8YU71+Dvf/hESelFSSQ+FGmKOhnAUwBmCSG6AUfoAJgpL5sDYId22055bI78XXjcc48QYgJAH4BpPmmZynY5ES0nouU9PT3lPaDE0VDso+ZSPoIwJprCTkQYrikuY2kdj8khHiaqBQg3sTFU50n+dvAwzxpWQAdpIKb3Gsa0FrjRVKUCxR0AhJsZDtg1kLxTvrxnCVoENRuyvv2qIRvy3VLANWE0lCDyJuDgNExO+aDdHMP6RPK+TfN5Jcj3Do6HSq8aVCxQiKgNwO8A/JsQwrxtmrzUcEz4HC/3Hu9BIW4UQiwWQiyeMWOGT/GCEUK4PhSTiq03anvYsP8cAMD+wXk6JUsfHo2GEi4N9QyVChRVhMI6c00NIT64sJs1BWkgpvcaZMIBSpvgZw+oKL/TU+jPMp41L/kR5EPJ+bSNwjyCnPKmJ3JNXr5RXuEnCftdo7SwigRKQEeukzTYvIIGmqW+d1u9qe8x3OrM1aEigUJESTjC5BdCiDvk4d3SjAX5Vy1ruhPAPO32uQB2yeNzDcc99xBRAsAUAL0+aVWVnMg3GOPGQQEjNx3fKC+bsNBuqTQWXX0kphFvmNVVAV1D8evMw5epcJQZNAEuF3IU650JX7qGEmZCqcrDpj8EdcK2vPP5qvdlvaSogOOWBQvdQY1tdrksh/1Z/AW0PlgyazD2c+41WtF8NZSAeShRRnn5lffoGa0AgqO8jCts6O04RIy07bt0zXuVj0vKppIoLwLwEwDrhBDf0E7dBWCp/L0UwJ3a8Ytl5NYCOM73p6VZbICIzpBpXlJwj0rrvQAelH6W+wCcS0Sd0hl/rjxWVXJCIOX6UII0FP+0wsyU93PKB80MDsTnujCrwTrngyc2htJQ3Lk93nTUnWHChsOGYZvq3WM6MwjHUMIg0Cnv3wkDQZF/vskb07Y55dUg2u6Utw829PNAiGhHQxbu4qN+TnmPhhLcyWYtz5opQUMJWn3Z73MYGnO0QVM70IVEkKncvpWx+XpvPnWUJJJKll55A4APAlhNRKvksc8CuB7A7UR0GYDtAN4HAEKINUR0O4C1cCLErhBCqNr7GICbATQDuFf+AxyBdSsRbYKjmVws0+olomsBPCOvu0YIYd7MO0KECB82bBv9Bdl8C9MpzD/wGmuq5jxMi9mpooXdxthXoJQwXCr8INS9YZz1oXf/M1zmNXmV3kH6lTF/3pyejq9AkX9LWRzSVqZ4gMlLdYq2vLzmO0MHmhPmi1X6IbSGsBMbg+ahjMoOOoxAzuRySMfixSfkvX4a++DYhDWfIFO4fk//6ATam5K+5bQOSCazQBFCPAp7mzvbcs91AK4zHF8O4CTD8VFIgWQ4dxOAm8KWNwpyImhio7+q70krZFSSjn7YrqGEa1RuB+WztlVQuGPeh1L6s5gofKYgU0PYqLrSnPKGJXVCvNdS5qHYpL4+ks3lhMfBG6Q15K8LLlOQH89tm1YfSoCAzvmfX7Or37d8hedCaSimgUBOYDQTXkOZyAqkDT1iKWt5BfmUgvqNTXsGMWdqs285bX1HPX0nCp4pXwICcE1eNtU1f94/rTAaSqFwCDNxMmz/7TrlfaK8Hly/x1dAhZnYWIpTvrBOlJAI4y8KK6DNPpT876APvtx5KN78gjWUwv01wibvGXRYblJyKtApb8kjSGhN+NTX0Fh+3lPYpVfK9aGMTuTNR2EFigl3YOO3OrIKNAioD1M71d/T8Jh9Xli+PKVruLWCBUoJ5IRAMuEfNqyWBg82gQR3gEVX6B+yNVTIN9uiy/zmoQDA8m37rWmoEVHhhxhmuXcdNWIu/FDyM6rN9+VCdjqeEbNxHoq/ScKroQTnYTwfok48k/kK69Q3de26MCavmL+ZSD1vGB+KScD61VeQw969LsBM5JenQt/YKszg3bZJVtBSREIIt7xGk5cmuE11rtfhcIjNuGzffmT7rlQAC5SQCCECZ8pnhcjPU7F9LAFRKUCwc9CWP+A0/sOnNAEA/n7xPOM1gL9pTM9n1OIkzOWE21kUOn89TuwKfCiqDiv3oei/DSNILV+T+S7M1s6u38HaCfuXwclH63gKyxHWlBli0OGGDQfsh2JrI56VgAMESuHZjMfZbky+OA8/35I8Zfqe9P1eKtJQLCtXmMoXNCAJ2phv2LKLpo6tHI3gQ2GBEhL1EoN8KEpDCV4QLthMVNhwhOEaUzlb0wm0pOJob7K7yIR2fXH++d8p01oSKDDPZAvNM+FGoQrTPBQlwAF7x6hHTlXilB/3ROEY5qF4Rv3mPEqZ2GjrqD2CzRLxVtLKyFYfisrDXzjacvIKWP/zheXVn9HvWcJqKK5Z1PBi1Cq/6UQslJZjX33ZX0MJGnAELXSply2MhsI+lIMA9QrzM+XNDSOsyct3dCZUngVmD8+I3FJO4XQYgZsg+ozs9I8iYREoHhu3j3kmnMmruCxe4WK+z6OhVOCU15cXN9VHGA0laN9wr4ZivmbCr05DDj7D7Mmjmq5tzoPb/mz1HuAT8Ju3o3fa/hMb9U7Y593Kc34aSnMqHqod2gSKSjtMQIZRwHqCLYrP6+kOVWTyYg1l0qBeetBaXn4mMQBub+uroVRo8lKOdr/m5UauGG265jx19NFQ4YcYZjFFE57Z0WGeNcQ1TnnMZVOMac5bW7RQPk9zHlfesdqav3NfsFDyq9MgrcFN22OOslwj07I65V1/QJkj8qxeBu/58YJINhtBZqLCspj2Q1EDheZkPJzJK8CnZDU1Vaih6Pf8ZvmOovMKdVkYH1y94C2AQ7Bz/zBeeLkPQPBqw8pEFGRrL2dio/5f3322KcRWsUqwGaOe8sesozat8Rb6ULyaVHAjN80lCDNvI2w+QWGsuoZi9qGUZos3EUZr8/pyCgRKSF9UUEQRkH/nQRMbbVkGjchtAx8hBM7++iMAgMOnNKHfZ6XrbEAnXJiXSRiopWeak/FQpqRoNJTia0rxsXT3jVrLF+RTtAUV1BIWKCF4x7f+6jZ+38UhRd7kZesA3BV6Q9iFi9A/1BCdrF/fp8xpQSYv26jHzzzjTcteBhfDLphBCwwWHg/jQ0nGyTjqHp/IIRknZLLC7EMJIdwUcdNMURRqOSHqtOB5VGcXJM+CHOb6NbaOWvWrYTrQoDBZ3TGuDzyaknHsH84Y0y9Mw888rNqeqRxjcg5KOqyGYmvr7hwTW1md86l4zNjZZ3w0Nifd/LEOi99T9ylaNdwG0FDY5BUCfSTlO1M+p5m8Aj5G25alQL5hFC4f7jHdWE1eDoRiH4znOtU4jQJFK4u107GPpsOYzEzYfCi2gZfHhxLC5BWPkdmHks2hORkvKkO+XLqGYkhfO2haywnwaij2aCG7yUu1iSATYhhToauhWNb6cjVkWx4BwlGvD1370+8L6uSz2iAgTACLqZ2OZZXJy9zRF2LdUyWkhpJKxIztQy+/39I+x85qs4el623dKlDqr6GwQCkRFeVls5W6a33ZlgYPYfKyaS9hfAaOUz7YK6/uDlqu3WYW0T++4rDhYO1CxxTlFWakHTbKS51LxszRPmMTObTKKdImQZ8Vec3D1IHqo/Awe4jY6kSv68L2E3bSWjaXn7ho9cXJ4u7qG7Gk4S+8PEIrwGdgMyemEzF/gZJVEVpx32f3W6B0vGSnvL/AsPpQsrpAMWgoIU1erekERjLZwLk9tnLYovZqCQuUEkm4uwgWn1u140DelGS4IJfTQmH9PhJ5b+HHGmaUAght61xrFr7qsz5Qs5sB7J2f/r8wAkVdYXPmVupDUR9o0hI+Oj7hr6Fkczn3vZuy0Z36NpNXmLL6+Q3crWyNd+bxBIYEmAp395vt9SrvME5ocwepDzbMAQ9NSfNovjCPVEDIrxKa44Y5NbpTPshiAJjbuj5pMUhDcbQpg4ANmNionqEtnUA2J4yLeoYx77KGMglJWD7W53ceAAA8ubkXqXjMOFoIY47QzxWOMD0dtU+nFI+ZFlQpRJkKTCq4vWPLH9e1GL95KIEF0ZatMDu/7SPtYC1GL0MqHjOOQscmckglYkjEyDJPAL7BFl4NxVaG4Hev12NhOdU9Az6ObHVdfosF/w7QlpbK29qBauUcM3TkenNQfgzA2+GlE/FQfq9UPBZSQzGYvGTZmpJxa114NKgy13FTbSYZj1lmyvunobJok1ry6HjAXCj2oRw8KKd8YQPtG8k7GBNxMs5CDhsKqT7owk7SY/Kyjj6dXePCLiJoXjW5uCyF+PlQ9PvDhA2r++1RXub7wvtQlIZiHkGOT+SQTsQQj5FdQ1Hv3eLUd68N4XC3mSb8tD79ng2vDBjvd8oqECNHU7KOyoW/QAnShvSimTQDmzlUr4OmpNP1BAUopAP8H37ryYUJG/bMQQpcJcFcBl2bCgwbNppUnWPK7DqcKX4v3m/KXA5dIP7XH9eaL6oyLFBKxLYOkv6Sk5ZRVViBokbqhdd4nN22UbsQ7ijZrzP325NCCOGaeGxqtN9o2rPmWAgVxRWgFo0jqpny6UTc2vEoDcXoQ8kJVzM1FUU3eVkDCHSBYlnyxG9xSP09dFt8H6qs8RghTmSNjlJl0cttKocQ5jakl8UkUHTtwtZhpxNxt7zG58jm0/DXUHLyWUzvVYYNSx+K2f+VrwNTW9fbS9AKB6l4zOJTMgcmKHSTF2CeJBum79Dr98ePbjFeU21YoIRAH+zbzAn6yCQZJ6Md1BbxUoi7y1yRySu4UWVzAnGiwJ391N1BEzSty3P4mbxCbCalY9ZQ8r9tH7J+eP+wfR9tvYMzaihZKVB8BgJ+Jq9RzawTaiVh63yH/PHCTkW/X8+vqKzCESixmL0sqh1mssI4oPAuj1J8v56/qSNX9zen4laT0qyOtFte23MAwT4Uv22o1TfYJP1jZv+XXr7iC9SGVzbtFfBGeZmEcEabn+a39llr2imnac5MmImxeh285bjKtjsvFxYoIdBDQVtTcmRVaI7SfifjMeMo1Kb+F5K1mLy8TnnzvTmR30fDry93/TQWm667JlnAZC+g8rBhdb9tFGZdQkQIdDQlkIrHsG3fsDV91Y8lE2aBMT6RQyouNRSL7ytv8ipOX3VK09tSIR3uwdfsG/IKyAmP38I+SS+XE4iR0lAsAsVnUqpTPm1Ubrh/ODPhfhN+GkpzMl6wTppz/DvvPxldrY5AsYdQO5q2TcgXlt9UDuW/yQsUk4aiCRRDGkqwd7YkrQulqvK5g82CbCayOaSTan6avb6UycuUTzhfjnN8wfRWxGP16dpZoIRADwVtThWr6jc/tgUf+ukz7v9tJq8wtnZA26ynoO3Zwmp1cjJslOAf5aVGdrZ5F8rGbXP0qU4tlSh2dPut5WQi446YzfVj60+EcDqc6W0pjw+rEJVWOh6zOG+zSLk+FIPvK6s5un18KB3NyVCTAU2dH+A1HfYOFgiUMKqevC4RI8QswtHJP99hjRm0nUyAE3lkPIu2pgQSMfJEcRWWtSUV96Sv6j4Zj7lmWav5RpruEjHy3yrYZ5O38WwORE6Isi0v7zpuBoGSUQIlhdFMzjekV0UKFs0hyglXqJnnoTjH2pUPxaCh6EUPWqC0JRX3HXRUExYoIdA/qjccM905pr3hnz25zf39g388FQmLyUsfDflNbOyV5puidZACFjFU98RDOOVV4zM6VbMi/3FYWq/aQ3tqc9Jg8irOxw/XxGeJ7AryF3U0J30FivqIm1LmyKLNPUNIxGSUl8Vpqkbkfnb4Fkv6gL+JMH9N/njv0JjxGQrTKiSTzSEphaPffBelaZs6ni17h9zfpiSGx7NoScaRSsSMAsntYFNxd3KhfjwZJ9cXaV0FQfmCLEJeoc5ZfWPxmBtIY15AMv/8JqGkOvfO1pS83q7RNVkESiYrXKFmXprF+as0FNO+8p4owQCTV2s6YSxnLWCBEgLV0M4/8TDMnuJsz6m/1LjWeZ94eIcTnmp4od9/eJP72750isCz250Q5MKOI8yaUk4nq0xePh1Pzm4qyAaMqID8vg2dLSnflXFLMXmVOg/FqR4KFig5Fe1TrE31DWcwkRO4e3U3EhYfSzYHbRXp4vTVx9uSTIQyZ9muUWVrSyeKTV565JRPZ5HJOp2ov8krhza5xEdhxzM4NoHtvXnzoU1DaU45AsU0cNJHyqaJjfqAxyYrHE3LboYEnG9FpWk0eclgCzU4MpmSvEEDdpNXV0vKmoYuQIFiwTSRzeW/J4vPEtAESoBT3s+HEo8RmpJx3zZSTViglMCJh3cYZyHrk9nSyRgScbMD746VL7u/bR/J01t6PZ2jJzooK3zXElPXxyh4HoqfMzOTE0i7AsXcMAflVqVTW4o1lDCzej35TRQ75b0RX7Y7BYiAKc1J9PsIFPWBNyXjRTbsIW1Do4Rf2LDPiFp1ylNbkhgZD57pbPvY1TUzO9LYdcAbyWVbH6sQZ12yGGIWDSWbE8jmhBtRVJjWUIgtaIfHJ9CSSiCdiFkHJADQnEp4zGt5DSWGuDJ5+ZgIVfizTQB7ZuRbLALpRMztzE0d9ViA1l+ooZgCIlT7b5H5FNaJY/Kym90KTV4mDcUzl8VSHyPjzgTddCLGGkoj8+ZjnYiJf37LMSAixKhAQ9EFSiLuOOV9Zq0eNb3V6uDbuGfQ83/d5JTJ5oLDLaXJC/DvzFX+5g4hh6S0X9s+5mFl8mpJejplIJwmpZMx+HP0ft8W5ZXJCiRjhPamhO+EPyUUm5PxIm1K1cPnLzjBmbthifKyOVwB4K8v9gAAprWlMJETFrNIsMlLHT9mRhvWdXvnmkxkBRbObPOU2cR41lmg1KahqDzam5IAik1ehfZ70/sb1jUUP4GS9Gow6j3EpY9Hv9aURiIes74T/VkKfyuc+UVxV3MwddRBK02re7panfry03KU1lds8sp/t0aTqhvlZfeheELTLZ/USCaLpqR6L+xDaVi+9O5X4ZH/OMvtqJMFM671RtLRlLAKlJZUHB9+4wJs3juEx1/ah/5R76g6mxP4/O9fAAB87h0nAPCOiCakhtKcjFtHko5T3gkb9uvK1YcyZijnRFY6RC2aFpDXUGa2N2F3/5hn9BfW3g9Is4U778bslPdbaiKZiKEtnSgSap7rcnkNpdAnpOp3bmcL4jEyvjddoJiE229W7AQATJORS4OGd6Mvdmhd502L0hnJZIv8LvO6WgAEaShZx+QVM89DGXcFillDUe/x1CM7AZg7r5FMFi2pONKJOEYNHVfeKZ/wOuU1H4oKo/VbMj5G5Gvy0jte4zyUrNfkZdJQdCFjNnnlTbuF1xfmrQRCobaU1TUUn4mPSiCNGNpymND00UwWzakYayiNzpypzThyWqv7//amhKfTUI3l0jcsABG5S6Hr5HLCcWam88tT7+n3Ol6//ucN7m81qhrL6E5Dx5zR2ZJ0HfeFKEd1IhbzdfyrBmqe6eyEySZjdk1raGwCzck4Fs5yRs0DYxntfr3xW4sAwPnY1PfhnSke/AFlso6zvCWV8DXVqGdIJ4t9JKqDaErGZMSauT6UU9Uv2qpLmkUGDdqSMus1Jc2TK/W0O5qd0bC+v3hWliERI98InkxWIJmwz0NR5XAFSoEJZ0TO0m6R7c8kQIfHs+4W0ybNUJlkmlOFYcNKQ4m5YbS2jk+ZGf00FCUgprWmjJrDWMYRrq5AMVyjl980D8U1efn4UPK+L7tTvsnHsqC+wfZ0AvEYWYSWrqHY64NNXpOQtnTC02k8t9PZfOuqdxwPAEYN5YC08Xe2JN1jhZPx/vTCKwCAeV3NbgemN4zRTBbpZAxHz2zDaplnIWqmdDJhHm0DTiehRkFWgRKTviCLUOoZHMP09hRaU8WORF2YBi29krFoM/p+4COZHM7++sP40E+fLrjXEbBt6TgyWWH1TUxkHV9LWmqWepnGXIESR0sqXrQFqxCOCUuNHv0ijjqlWeTOVbuKzqkBSGdLymjvB4CXpd9Edfb6CHxwLO+38JvYODyeRXMyYTV5qbzVe7OZvNT5wiSEENi0ZxDNyQSmNidxwLCniRKMram4pwNW7zoRI7eDtQnH0YwzdyMZN5vVgLz/q6vVHNLraig+Jq9BfSBk0lCUyavN7kMZL9RQinwo/k559Y2nk3HrZmAeDcXy+kcySqCwU35S0ZJK4BFpN9dRZpG2dKLISazCQLtaU/jjv74RALBPm2vwg0dewmYZrjmWyTdA9UE+8mIPfr9qF3b0juCkOVOwee9QkXNu5fb92LhnELEYOR+ipeN6y9cedjtOc9iwMzp0JpWZ03ilbxSHdTQZP1bPzoMBKopuqtKvVWWf2ZHGgeFxvNQzhIc2eOtcCRT1Idu0lEwuh2Qs5k720oukyt2cjBcNFID8x66c2LY9RAC46X/z/heLzj20YQ/iMTL6cQCnPn/51HYAwNxOJ5Jwh4y2yuUEegbGMLMjjaakfY7Buu5+rOvuR/9oxpmHYtRQvOaVwvevBgYdzVKAFrSh369yAkse3rAH+4bGsfrlPo8mBShbfgwdTUmMZnJuHuq5E3HKaygW4dg/msGU5iSak3GjIHDK4LSHmXLWvSnAoDmZ96GMGjrqgdEJEDkaqqmzHxnPuoEfgM2H4hxz24hhLxs3bNgndDmtzLeGdqwGJDGC0cwIeH0oPA9lErG2ux99Ixk8t+MAvnzPOgDAa+dNdc8f0dWCHftHcEDTQB7duBcAML0t7ZpGdA3l+nvXu79jRJpAcRrn0pvyo/POliSyOVHkN3j39x9370/5BAZs1WaVm4SOcromY8WmO6dMWWzcM4hZHU1G+7QuhAZG7dFXQL7TBLw+FNUJzWxvwv4hs3lPmeaUQDH5LgCn01LRd4D3gx/VZlObfDFjBQ7Xwvr61dPb3d8XvGq2Mf99g2N4dvsBZHMCrek4egbHiq55WYvqUrPIlY9t//A4JnICM9vTaGtKoH/E/Jw/+stmAMCq7QeQjMWMI+7frHD2LM875QtNXtKM1OaUoXC0/OJuJ2hk39A4npda8jNb93uuGRqbQGsq4ZrueuX7U4KnJZlwndQ2YdE/kkFHUxItKfv2veqbUb4rvQ3mcgLruxePfp8AABweSURBVAewcFabq20V+iwBR6C0pe1+z017BjG3s9nVqEwCZeu+YTQn4+53PV4w6Mhkc/moSZNAUTtLJmJos5gRn3jJ6T/mdbVYA1AcH0rcjb4LszBr1LBAqYCh8Qn8UH7E7zl1rnv8tPldyOaEG7G1d3AMX/yDs/pnV2vKbXi9lo4yHiPN5OU04COkQ/ZvT57jdgZ6w9IbT1MiVhQ4oCj8IHoGiju33qFxdLWm5Ein+CNbuX0/eofGcc6iWXkNxSNQnHynNCfR57PNKwD8rRSC09vSno5DlX1WR9q693gm62geamRocsw/uXkfbn58KzqaksYJbqOahtJq0lDk+XbL6POqO1YDcMyU8Rjhn14/3zVZKXRB9+q5U7G+u7/oY+8bybcFVU7VMe2R72hmexNmtKWN7wzIB1jMmpJGazpuFLDfeXCTTMs8qh/W/BJAcYffIdveFW85Gqcv6AKQnyGuGBqbQGs64frX1nX3A8i31/amhNuebSsn949OoKM5gaaUWUMRQuDIaU4apx/VVVTW3uFxDIxN4JiZbThsShNiBLy8v3hRzf7RjGwbZoGy68AIFkxvy2s5hpF/79A4ZnWkXb+Tnk5ORv25SzYZNP5RuVIDERX5ZxW/lYEfR05rNQ7S+kczeH5nH5LxmLMVsci381U7DuDTtz/nu6hoVLBAKYOz5MJrurqu+0ZUQ1frS93w8EvuufnTWl17vRIoz+044EnfUcG9GsqR01pw0pwOfP19r3E7LL3h6aPWllRcBgYUN1590pr6vy5kbn1yG/pGMuhsSWFaWxr7CkbTLx8YwQd+9BQA4OR5nUaHp3L8Bi2Jcueq/LycJQs6PZqIKvvM9ibr/eMTOSQT5GvyuvjGJ91yq9029UmnSgg1p+LO6FCv09EM/vB8N4D8iN6m9SmB2pp2IvB0gaG/p6NmtGJoPOsKCYXui1CDCZXXbVILmtmRxoz2NPYMmDfGuluW9d0nz0V7U7JowKG/58XznSiuQtOI6+iWPoPCzk1N8Lz8TUfjM+c7PsNCk1fvcAZTW5I4XE4C7h0ax2gmi2vkkuptTQnM7WxGjPJ+w0JcDSWZwPhErsgf9LuVL2PbvmGkEvkBhd4G+12fZQrJeAztTUnXj6mz5uV+zOpIW1dJeKV/FLPa027gzUjBXiVCCNz13C4Mj2eRihcLFCWAVPsxaSj9IxnXpNaWThgHUKcc2Yl5Xc2Y2Z42aig/+esWAMCD6/cUBTxs2jOI363c6WuujQoWKGVw7YUnAXBMXwo1kxZwQlBjBGzf5/hE1MhobmezO9Lpak1hc88gHnmxB5v3eueemDSU/tEJdLakEIuRpqHkP5BXtN33xqVv4cBwpuhD3Cr9NPqocqss54HhcXxBhi13taYwoy1d1PF98a417u/ZU/M+lIc1/4bqhOZ2tuDASMaoeo9msvjEbavc/89sb/LMDledw2FTvAJFT+vASAZTm1NudM3gmL/dWDnNlanx8Zf24ldP7wCRMyBoSzmdlzK3/b+71uJa1Qm6Gor5o9wr/WHxmDM6fHSTY6J4YN1uXPDtRwEAn33H8Vgw3YkW3Nwz5Ll/vxQo3/+HU1xfnCqHMpEdf1g7Zk9pxks9Qx5BPzaRxYXffdT9/yfOXlgUgfWrp3fg+C/8CQDwqjlTcPhUp7MvFMKq3pUGURiJqDS41nQ8rxkW1Pue/lHM6mjC9PY0YgS8uHsAK7fnzWLJuDMazwngic37CqsSgNQcmpNoTsU85VI8s6XXrSPTxEUlPFRH3dFcbEradWAEG3YP4OQjOo0a/brufuzuH8OsjiY3nQMjXquCCsjZMzCGZKJ4wcwhbb4WYPYZ9WkCpaMpiUGDBtI7NI6FM9vRUTBQUKjVMS59w3w3JFuV499/85xbB9WGBUoZzOtqwZHTWvDf9+XDfI+QWgngjOJmT2nGNqkNvNI/iqOmt+Khfz/LvWZaawoPbejB0puedsNm/+n18wE4S7moUXd3nyMoBqRqDuSjgHqHzALl5HmdSCViWP1yH/5Ji4xatnY3Lr91BQDgqc+djZs/dBoAuGYpXXhMaU5iRnuxeUUXnMl4zC3Trdp6Zr9+xrHTL5jeimxOuJ2tju47ufcTb0JXawoDoxMYn8hheHwCtz7hpDdHdnyK/dpI/sBwBp2tSbcMjxQ47X+ulenvFs/FdOkXuPXJbRBC4AM/egrruvvR1ZJCIh5z/SSqk9WXrlGLIQbtmHiaHPn/ZvlOPLh+Ny67Zbl77rwTD8Pxh3WACFi+tdc9nsnmXMF1/omH5bcOkCNdAuHoGa1ob0rilCMdX93pX3rAvX/1zj63Y5sztdkddOgDjpse2+L+/sbfvQbt6QSmtabw0p68YPveQ5vw3/dtQIyAo2c45qqd+/Pvad/gmBtwkIjH3OXWC02Nu/tHMasjjbZ0AifM7sD6VwbcTq4tXdypFZpwrrpjNUYzOXQ0JdCcUhFv+Tw2vDKAXy932th/v/fVRUurbNw94PoTZ3U4A5L2dPFqCq+//kEAwClHdGJaWwqbevIDu+37hvH2b/1VppFGSyqBllQcewfybXnbviF898GN7v+VUNDb6PZep347W1Mygqu4/RwYzmCqJvgKI+cmsjms2dWPaa0pTGlOYnBsokizfHHPAKa3pfHZd5yAKS1e35XCVPdRM6kFChGdT0QbiGgTEV1Zy7wve+MC9/eFrz0ccztbPOdndaRx56pd6BvJYHvvME4/qsvtKID8BkNA3gz1+qOnAXB2XDxqeisWTG/FXc85IagDoxOuIJkunZD/52fLcdOjWzCayWK3FDz/+8+vx/sWz3VHJ3+VwQAA8PFfPev+7mhKus7Mn8jNeJR5Kp2I4azjZmBmexp9IxmPqWRX3wi6WlNY+YVzAOQ1iPdKH9L+oXE8sH4PAODN0jT4wLrd3rq7+Rmc882/AABu/8jrcMLsDtev9MzWXvziye2u9nf+SYdhWmsKr5k7BUB+4cJMNofeoXFMbU65nd9Nj23BZtkp7OgddieJvvX4mfjqe1+DGdJv8NPHtnqEjfoAC537urktk82hqzWFh9bvwXcf3IjRTNYTsXPKEU5H/6aFM/CaeVOxY/8wLr05L0wAx/49oz2NBdNb8fVlL+Lzv1+NA8PjWPi5e91rVISeyhMA7l+32+0YZ8v6nsgJVyjrDv1jpd+iQzrv3Y3UNM1u4ax2EBGOmdmGl7ROVA2QulrTmNKcRHs64RH8Z3w5L8SAfGixruWMZrLYP5zBLFl386e3Yuu+IfzgEcfXqAYxAHDthScCgGfrgR29w26ggz4pcVjTgs77n7+4v99zylzXPKcGP/q7nSMj5jqaE3hg/R4s/Nw9EEK4W3YDwElzOnDWcTPx/M4DblvXrQYqjeltaezTFu38yK0rcP86p63f8A+nuN+TKsf4RA7vv9ExD5++oMvx0WnPIYTAVXc8j8dfymtpc6Y2Y580ESrWSz/TghmtOGKaU5btBds1rNp+AGcc1QUiwsKZ7QCAldv2u/1AMk7uJnHVZNIKFCKKA/gegLcDWATg/US0qFb5X/K6+e5vNfLVUaOUM770AHqHxt1RtGL+9LwA+p/7NyIZJyw6vAOA0/hiMcJp8zvx9JZezL/ybvQMjLkfjm4GuuaPa3HWfz/sRu8sOrwDRIQfL10MAJjRnsaz2/ejZ2DMNRt89T2vBgAcd5jT8P68djceebEHu6WW8/sr3oCpLSnMltrB/5MBBRPZHFZu24+zj5/pCgAAOHpGK57e0ovbl+/AX6WpZ/aUJpx17AxMb0vjrxv3om84g5HxLCayOVfgHDOzDUukY1d1lNf8Ya1nxJqMx7DiC+fgG3//WgCOvymXE24nvG9oHLEY4Uy5PM5//u55AMCzml9KjRxnaO/pC3fmTXdvO2EWAKcTVmlu3D2AHzyS9331jWSwZ2AMG3YP4Gt/fhF/851H3aihM47qwm2Xv8699owFXW4ElAnlEP/5k9vx0Z+vcI+/89VOlJgyI978+Fa88HIfxiZyrjDUJ9i+6asP4ao7Vns6/cvPPBqAY14dz+bwsAxv7xvOoL0pgU++7Vj32uZUHMu37cfwuHfE++5T5oCIMKezGWu7+7FK1qUyCf3HeccBcLS2GOXNdy/1DOLC7z4GIL8p3YJprdi2bxj3y0GFEoyAI3wBx4+yd3AMT23ehzW78vU2PpHD0TOc512xzTGZ6R3tR958FGIxwtypzre0c/8IsjmBO5/LzwNS7z7vAxPoGRzDu7RyHjmtFQumt0AIpxMGnHYIOKP6M45yBnoz2vNrrO0bHHM7+otPm4e3v2o2Ugln0nHPoPMd/XntK25U4KyOJqQTMfzq6e3uQOHF3YP41dPOd6tW4VB9idIuRjNZfOwXThu54FWzcdT0NlnXec3yi3etwcsHRtz+48TDOzC9LYV7X+h2zV3KTF9tqq8DVY8lADYJITYDABHdBuBCADXbTPlHlyzG//nZcrdT1Pnx0sU4++uPuJ34afO913zugkXo7ht1NYhjZrZjbmcLln3yTLfT+MDpR+L25Tvde9SIOZWI4f5PnYm3fcMZqb3SP+qavJTm096UxLUXnYQv/P4F/O33H8fRM1qRSsTw/tPm4e9Om+em89E3H40fPPKSJyxZ2c/f+erZ+NFfNuNXT2/3hMeqzlsxtSWFFdv24zO/fd499pfPvAVEhAXTW3D36m7cvbrbc89r5k7BDz54qvv/tx4/E50tSWzYPQABgVkdaTx51dnu+aOmt+Ko6a245o9r8Yfn8x3Gp891Osgv/s0ivPXrj+CZrfsx/8q7PXld+gZHm+zUzHWKGz94Ks498TAAeTPPRd97rOi6MxfOwPGHtbudyMY9g3jtNcsAAH9/2jzXWQ04HY/uuzr+sHZ89wMnu/+fpgm2Jzf3yrIl8a2LnWuUaeKlniG88zuOb0QJveltafzHece52oR6L3OmNuOxK9/qpvuquY7GpO/T86lzjsXHz17o/n/tLkcLXPR/78P/fWd+LHaxbB8z2p3BwEXfewyHS4H/zlfPxhVvOQaAI+zndDbj1ie3eUyeAPBBOeB683Ez8N2HnMiyo6a3usvHAE6gyRFdLfjuQ5vcaxTHzGzD+xbPw8z2NJJxwqd/8xw+LTtHhRKOyjdw3T3rcJ0M4weALV9+h/v7hMPasWytI9SWXJfXtO75+JtkWeQ39+OnPHms+r/nuCP7kfEs1nb346O3rnAHE1/8m0X4xzOOdK/P5gR+/uR2/PzJ/Pei6vbtJx2GHz/q7J3Ulk7gT2vyAQmfOsd5FuVref31D+JT5xyLhzfswY7eEUxvS+HIaa3oanXyveqO5zGacebI3Pz4Vqeu5XdJRFiyoAv3rM6n/6Zja7OD46TVUADMAbBD+/9OeaxmnLNoFlZ+4RycJzsknaNntOHqv8l/pGefMNNzfkpzErd8aIn7/++83+lMFs5qdzun186biivffrx7/fuXHOFef8zMdrz0pXfg0+fkR5wfO+toTx6L5XpMgNM5jU/kcFqB8PvMecfhP2W0DuCEvyrTT1Myji+8s1jpO2fRLM//P/Zmb77HH9bumm1U51PIl979KncrAMD5CK6Ro6gXdw/i7SfN9uzpQkT425Od16uW9//dx17njuiOmtGGD2tmSMWWL78Dr5LmsliM8OlzjnXNc4AjyBRHSYGi89X3vBpbr78A87pacPfH34SPvPmoomuOm9Xh+f8Fr/bOR/ncBSfgGGmGAIB/fWtxnaz8wjmeRUbVe1coZz7g1On7lxzh2eb5zGOne64//rB2FFI48Ln+Pa9yf6sIrIf//Sy3HpQQA4Bd0qRaWEcfOdP77gFHw1WawWnzu/DDD56KNx4zHfd98kzPdUSEn192etH973rN4bj/U2/GrI4mEBG++K4TPeent6Xxl/94i+uMJyKceLj3HfzoksWe9nOGNCfr/OLDp+OE2c59JxzWUXT+q+99tcdMpN7bn9a84pqpPvi6+Z5rPl/wvVz02sNxqWyX6vt8dNNeV5hMaU7ioX8/C6dLLWjR7Cnuvd9Y9iJWbj+As4+fiV9/xNGA25uSOGF2B/YPZ/Bvv17lBrYs++SZOPHw/L1/f1q+r/jP848v8kVWC6rH5JcoIKL3AThPCPFh+f8PAlgihPjXgusuB3A5ABxxxBGnbtu2rSitapHNCdyzuhsLprfipDlTjNeMjDuzioM2xPJjxbZezOpoKvLjAMBdz+1CNpfDL57cjtcdPQ2ffNux7kqvOtv2DeHPa3bjopPnuOYVxUQ258yal6q6zRa7dlc/FkxvLXoeIQRueOQlbNo9iGNmteHSNyxwO4NCVu/swzNbe/GPZxzpGfWrdDbsHsAdK1/G4VOasPT18z35ZHPOsiCrX+7D3sExXPbGBR6/lZ7On154BUsWdHm0BcAxsyRihLXd/ZjWlvIIPZ0bHn4Jv1m+A9/7h1PcTklnNJPFb5bvwPbeYXz2HScY3+/G3QN45MUenHnsDBw7yysAhBDY3juMnftH0N036hGChdc9vaUXp83vKnqvo5ks+kcy+P2ql3HyEZ1FWrLiq39aj+8//BKuf/ercLE2aAGciK0H1+/BD/+yGacc0Ymr37WoyHy7fGsvVmzbjzcunA4hYG3rNvYPjWPf0DjW7OrDtn3D+Je3HFP0LPsGx/Dk5l6MTWTxN685vOi99o1k3GjKGe3pojYshMB9a15BV2sam3sGcfzsDs9kZMBpP7v7R9HelMCuA6NYOLOtqBy5nMDmvUN4YN1uHDalCRe+1jyGfWZrL3oGxnDeiYd5Bgq3P7MDG3YPYPaUJpx34mHoaE66wlfx2Ka9eHD9HrzluJkYHMvgbSfM8nxzvUPjuPeFbqzZ1Y+jprdi8fyuomcB1FYGwZvthYGIVgghFgdeN4kFyusAfFEIcZ78/1UAIIT4su2exYsXi+XLl9tOMwzDMAbCCpTJbPJ6BsBCIlpARCkAFwO4q85lYhiGOWSZtE55IcQEEf0LgPsAxAHcJIRYE3AbwzAMUyUmrUABACHEPQDuqXc5GIZhmMlt8mIYhmEaCBYoDMMwTCSwQGEYhmEigQUKwzAMEwksUBiGYZhImLQTG8uBiAYAbPC5ZAoA+6p+DkcA2O5zPiiNMHlEkUal5YyqHEHXBJWzVuU4WMoJ8LuvRhqH+nd/nBCieD2fQoQQh8w/AMsDzt8YIo2eStIImUcUaVRUzho+i285G6VOJ0s5+d03Zp1OlnLargnqO9U/Nnl5+UOIaw4EnA9KI0weUaRRaTmjKkfQNUHlrFU5DpZyAvzuq5EGf/chONRMXstFiPVoqp1GLeByRstkKScweco6WcoJTJ6yVqucYdM91DSUGxskjVrA5YyWyVJOYPKUdbKUE5g8Za1WOUOle0hpKAzDMEz1ONQ0FIZhGKZKHPIChYhuIqI9RPSCduw1RPQEEa0moj8QUYc8niSiW+TxdWoPFnnuYSLaQESr5L+ZpvxqVM4UEf1UHn+OiM7S7jlVHt9ERN+mKHbfqU45q12f84joIfke1xDRJ+TxLiJaRkQb5d9O7Z6rZL1tIKLztOPVrtMoy1q1ei21nEQ0TV4/SETfLUiroeo0oKyNVKfnENEKWXcriOitWlpVrVMAh1bYsCVE7kwApwB4QTv2DIA3y9+XArhW/v4AgNvk7xYAWwHMl/9/GMDiBinnFQB+Kn/PBLACQEz+/2kArwNAAO4F8PYGLWe163M2gFPk73YALwJYBOCrAK6Ux68E8BX5exGA5wCkASwA8BKAeI3qNMqyVq1eyyhnK4A3AvgogO8WpNVodepX1kaq05MBHC5/nwTg5VrVqRAcNgwhxF8A9BYcPg7AX+TvZQDeoy4H0EpECQDNAMYB9DdgORcBeEDetwdOKOFiIpoNoEMI8YRwWtjPAFzUaOWMsjw2hBDdQoiV8vcAgHUA5gC4EMAt8rJbkK+fC+EMJsaEEFsAbAKwpEZ1GklZoyxTFOUUQgwJIR4FMKqn04h1aitrtSmjnM8KIXbJ42sANBFRuhZ1CrDJy8YLAN4lf78PwDz5+7cAhgB0w5mN+jUhhN55/lSqvF+oijoZvpzPAbiQiBJEtADAqfLcHAA7tft3ymONVk5FTeqTiObDGdk9BWCWEKIbcD5mOJoT4NTTDu02VXc1rdMKy6qoer2GLKeNRqzTIBqxTt8D4FkhxBhqVKcsUMxcCuAKIloBR80cl8eXAMgCOByOKeHTRHSUPPcPQohXAXiT/PfBOpbzJjgNZjmA/wHwOIAJOKpuIbUI8yu1nECN6pOI2gD8DsC/CSH8tE1b3dWsTiMoK1CDei2hnNYkDMfqXad+NFydEtGJAL4C4CPqkOGyyOuUBYoBIcR6IcS5QohTAfwKjg0acHwofxJCZKSJ5jFIE40Q4mX5dwDAL1EbE4OxnEKICSHEJ4UQrxVCXAhgKoCNcDrvuVoScwHsKky3AcpZk/okoiScj/QXQog75OHd0jygTC975PGd8GpPqu5qUqcRlbXq9VpiOW00Yp1aabQ6JaK5AP4XwCVCCNV31aROWaAYUFEaRBQD8HkAP5CntgN4Kzm0AjgDwHppspku70kCeCccM09dyklELbJ8IKJzAEwIIdZK1XiAiM6QavklAO5stHLWoj7l8/8EwDohxDe0U3cBWCp/L0W+fu4CcLG0Ry8AsBDA07Wo06jKWu16LaOcRhq0Tm3pNFSdEtFUAHcDuEoI8Zi6uGbfftRe/sn2D86IuRtABo4UvwzAJ+BEU7wI4HrkJ4C2AfgNHGfXWgD/IfIRICsAPC/PfQsyqqZO5ZwPZ1XldQDuB3Ckls5iOA3+JQDfVfc0UjlrVJ9vhKPyPw9glfz3DgDT4AQKbJR/u7R7PifrbQO0CJka1GkkZa12vZZZzq1wgjgGZXtZ1MB1WlTWRqtTOAO2Ie3aVQBm1qJOhRA8U55hGIaJBjZ5MQzDMJHAAoVhGIaJBBYoDMMwTCSwQGEYhmEigQUKwzAMEwksUBimQSCijxLRJSVcP5+0VZ0Zpt4k6l0AhmGcCXJCiB8EX8kwjQsLFIaJCLl435/gLN53MpyJnJcAOAHAN+BMjN0L4J+EEN1E9DCc9cveAOAuImoHMCiE+BoRvRbOigItcCaiXSqE2E9Ep8JZA20YwKO1ezqGCYZNXgwTLccBuFEI8Wo4WxtcAeA7AN4rnLXMbgJwnXb9VCHEm4UQXy9I52cA/lOmsxrA1fL4TwF8XAjxumo+BMOUA2soDBMtO0R+DaX/397do0QUBFEUPhc1MTEydQUuQRfgEkRkmG24DE00EUQTMTURo8ncwGDmFnQQTdvg9STiDw4FY3C+sIKiX3Sp5lF9BRwxPHR037earzCsppm7/twgyQZD0Ex66QK4+aJ+CezVf4K0GANFqvV5l9ErMP1honj7Q+980V/6N7zykmptJZmHxz7wAGzOa0nW+lsV32qtzYDnJLu9dAhMWmsvwCzJTq8f1B9fWpwTilTrERglOWPYBHsC3AHH/cpqleExsekvfUbAaZJ14AkY9/oYOE/y3vtK/4bbhqUi/S+v29ba9pKPIi2FV16SpBJOKJKkEk4okqQSBookqYSBIkkqYaBIkkoYKJKkEgaKJKnEB0IbTWG6KkNRAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -216,9 +2237,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+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXmYnFWZ9/+5q6qr9056zdYJIQuSEDRIDCDuyCI6AyrMxNcRxsHBQZxXZ+Y3/sTXeRnxx7ziq6KMAzMoDAG5BAYdQQQxbCqyJUFICCGkQ7bO1p3e16quqvv3x3Oe7uru6urq7urq7vT9ua66qvrUc06feq6q5/vcyzm3qCqGYRiGkQmBqZ6AYRiGMXMw0TAMwzAyxkTDMAzDyBgTDcMwDCNjTDQMwzCMjDHRMAzDMDLGRMMwDMPIGBMNwzAMI2NMNAzDMIyMCU31BLJNVVWVLl26dKqnYRiGMaPYunXrcVWtHu24E040li5dypYtW6Z6GoZhGDMKEdmfyXHmnjIMwzAyxkTDMAzDyJiMRUNEgiLyRxF5xP1dISKbRGS3ey5POvY6EakTkV0icmFS+5kist29d4uIiGvPF5H7XfuLIrI0qc+V7n/sFpErs/GhDcMwjPExFkvjS8DOpL+/CjypqiuBJ93fiMhqYANwGnARcKuIBF2f24CrgZXucZFrvwpoUdUVwM3ATW6sCuB64CxgPXB9sjgZhmEYuSUj0RCRWuCjwI+Tmi8BNrrXG4FLk9rvU9WIqu4F6oD1IrIAKFPV59Ur4nH3kD7+WA8C5zkr5EJgk6o2q2oLsIkBoTEMwzByTKaWxveBrwCJpLZ5qnoEwD3XuPZFwMGk4+pd2yL3emj7oD6qGgPagMo0Yw1CRK4WkS0isqWxsTHDj2QYhmGMlVFFQ0Q+BjSo6tYMx5QUbZqmfbx9BhpUb1fVdaq6rrp61DRjwzAMY5xkYmmcC/ypiOwD7gM+JCI/AY45lxPuucEdXw8sTupfCxx27bUp2gf1EZEQMAdoTjOWYcxqXj7QwmuH2qZ6GsYsZFTRUNXrVLVWVZfiBbifUtW/AB4G/GymK4GH3OuHgQ0uI+pkvID3S86F1SEiZ7t4xRVD+vhjXeb+hwKPAxeISLkLgF/g2gxjVvO/H3qN//v4rqmehjELmciK8G8BD4jIVcAB4HIAVd0hIg8ArwMx4FpVjbs+1wB3AYXAY+4BcAdwj4jU4VkYG9xYzSLyTWCzO+4GVW2ewJwN44SgqTNKYV5w9AMNI8uMSTRU9RngGfe6CThvhONuBG5M0b4FWJOivRcnOineuxO4cyzzNIwTGVWluStKVUn+VE/FmIXYinDDmGF0R+NEYgl6++KjH2wYWcZEwzBmGM1dUQB6YyYaRu4x0TCMGUZLtxONvsQoRxpG9jHRMIwZRr+lYe4pYwow0TCMGYZvaUTM0jCmABMNw5hhNHV6ohGNJ4gnhm2QYBiTiomGYcwwfEsDIGLBcCPHmGgYxgyjuauv/7UFw41cY6JhGDOMlq4BS8OC4UauMdEwjBlGs4mGMYWYaBjGDKO5O0oo4FUNMPeUkWtMNAxjhtHSFWVeWQFgq8KN3GOiYRgziERCaemOsnCuEw1zTxk5xkTDMGYQ7b19JBQWzCkEbIGfkXtMNAxjBtHkguAL53qiYZaGkWtMNAxjBtHSLxoW0zCmhlFFQ0QKROQlEXlVRHaIyDdc+z+LyCERecU9Lk7qc52I1InILhG5MKn9TBHZ7t67xZV9xZWGvd+1vygiS5P6XCkiu93jSgxjFtPS7S3sm+8Hws09ZeSYTCr3RYAPqWqniOQBz4qIX6b1ZlX9TvLBIrIar1zracBC4AkROcWVfL0NuBp4AXgUuAiv5OtVQIuqrhCRDcBNwJ+LSAVwPbAOUGCriDysqi0T+9iGMTPxtw0pLw4D5p4ycs+oloZ6dLo/89wj3S5plwD3qWpEVfcCdcB6EVkAlKnq86qqwN3ApUl9NrrXDwLnOSvkQmCTqjY7odiEJzSGMSuJxjzLoqwgDzBLw8g9GcU0RCQoIq8ADXgX8RfdW18UkW0icqeIlLu2RcDBpO71rm2Rez20fVAfVY0BbUBlmrGGzu9qEdkiIlsaGxsz+UiGMSPxRaO0wHMSmKVh5JqMRENV46q6FqjFsxrW4LmalgNrgSPAd93hkmqINO3j7ZM8v9tVdZ2qrquurk77WQxjJhONe6KRHwoQDgUsEG7knDFlT6lqK/AMcJGqHnNikgB+BKx3h9UDi5O61QKHXXttivZBfUQkBMwBmtOMZRizEt/SyAsFKAgFbJ2GkXMyyZ6qFpG57nUh8GHgDRej8Pk48Jp7/TCwwWVEnQysBF5S1SNAh4ic7eIVVwAPJfXxM6MuA55ycY/HgQtEpNy5vy5wbYYxK4k40QgHAxTkBa2ehpFzMsmeWgBsFJEgnsg8oKqPiMg9IrIWz120D/g8gKruEJEHgNeBGHCty5wCuAa4CyjEy5rys7DuAO4RkTo8C2ODG6tZRL4JbHbH3aCqzRP4vIYxo+mLDxYNC4QbuWZU0VDVbcAZKdo/k6bPjcCNKdq3AGtStPcCl48w1p3AnaPN0zBmA9FYgrygEAgIBXkBC4QbOcdWhBvGDCIaSxAOej9bz9Iw0TByi4mGYcwgovEEeSEnGiFzTxm5x0TDMGYQyZZGfp6l3Bq5x0TDMGYQ0XiCcCjZPWWWhpFbTDQMYwYRjQ0WjYjFNIwcY6JhGDOIQYHwkGVPGbnHRMMwZhDJ7ikvpmHuKSO3mGgYxgxisKVhKbdG7jHRMIwZRN+wQHgcb8cdw8gNJhqGMYMYHAgPkFDoi5toGLnDRMMwZhCRWIK8pBXhYHXCjdxiomEYM4jBgXAnGhbXMHKIiYZhzCCisQT5SSm3gNXUMHKKiYZhzCCGBsLBLA0jt5hoGMYMYuiKcMC2EjFyiomGYcwgooMC4d6zBcKNXJJJudcCEXlJRF4VkR0i8g3XXiEim0Rkt3suT+pznYjUicguEbkwqf1MEdnu3rvFlX3FlYa937W/KCJLk/pc6f7HbhG5EsOYxQwKhIc8S6PPVoUbOSQTSyMCfEhV3wGsBS4SkbOBrwJPqupK4En3NyKyGq9c62nARcCtrlQswG3A1Xh1w1e69wGuAlpUdQVwM3CTG6sCuB44C1gPXJ8sToYxm0gklL649q8IzwsK4AmJYeSKUUVDPTrdn3nuocAlwEbXvhG41L2+BLhPVSOquheoA9aLyAKgTFWfV28J691D+vhjPQic56yQC4FNqtqsqi3AJgaExjBmFX0JVx885IuG92yL+4xcklFMQ0SCIvIK0IB3EX8RmKeqRwDcc407fBFwMKl7vWtb5F4PbR/UR1VjQBtQmWYsYxrz2qE2bvr1G7a9RZaJOjdU/jDRMEvDyB0ZiYaqxlV1LVCLZzWsSXO4pBoiTft4+wz8Q5GrRWSLiGxpbGxMMzUjF/z6taPc9sweIuZrzyq+aOQNcU+ZaBi5ZEzZU6raCjyD5yI65lxOuOcGd1g9sDipWy1w2LXXpmgf1EdEQsAcoDnNWEPndbuqrlPVddXV1WP5SMYk0B2ND3o2soMfuzD3lDGVZJI9VS0ic93rQuDDwBvAw4CfzXQl8JB7/TCwwWVEnYwX8H7JubA6RORsF6+4Ykgff6zLgKdc3ONx4AIRKXcB8AtcmzGN6emLAdAViU3xTE4sfEtjIBBu7ikj94QyOGYBsNFlQAWAB1T1ERF5HnhARK4CDgCXA6jqDhF5AHgdiAHXqqp/y3kNcBdQCDzmHgB3APeISB2ehbHBjdUsIt8ENrvjblDV5ol8YGPy8S2MrqiJRjbpF42QuaeMqWNU0VDVbcAZKdqbgPNG6HMjcGOK9i3AsHiIqvbiRCfFe3cCd442T2P60BWJD3o2ssMw95R7jlrsyMghtiLcyDq+e6rbLI2sMtQ95T/HEhbTMHKHiYaRdfrdUxbTyCrD3VMupmGWhpFDTDSMrNMTNffUZDDUPRUMCCIW0zByi4mGkXUGUm7N0sgmQ91T4FkbUUu5NXKIiYaRdXyx6DRLI6v0DbE0wBMQszSMXGKiYWQdszQmh8iQFeHeayFmomHkEBMNI6uoKj19FtOYDIbuPQXmnjJyj4mGkVV6+xL4+xSapZFdhgbCwRMNc08ZucREw8gqyULRaSm3WSV1IFxMNIycYqJhZJXkTQptw8Ls4otDnlkaxhRiopEh2+vb+MqDr1rQcRSShcIW92WXkVJubZdbI5eYaGTIvzy6kwe21PPHg61TPZVpje+eCgbELI0sM1BPY6DMjLmnjFxjopEB2+vbeP6tJgB+u8uKPKXDXw1eURw2SyPLROIJwqEAXmUBD3NPGbnGRCMDbv/9W5Tkh1i9oIzfvmmikQ7fuqguybet0bNMNJYgPzj4J5sXDNAXM/eUkTtMNEbhtUNt/GrbYT591hIuPn0+2w+1cbwzMtXTmrb4QlFdmk+3rdPIKn3xxKAgOHhB8ahZGkYOMdFIQyKhfP0Xr1FRHOYLH1zB+0+pAeD3u83aGAnfPVXlLA1VuwvOFtFYYlAQHCAcFGIJEw0jd2RS7nWxiDwtIjtFZIeIfMm1/7OIHBKRV9zj4qQ+14lInYjsEpELk9rPFJHt7r1bXNlXXGnY+137iyKyNKnPlSKy2z2uJIf8ctthXjnYynUfWcWcwjxOW1hGWUGIrftbcjmNGYXvnqoqDZNQb7GfkR2iscSghX0AoYC5p4zckkm51xjwD6r6soiUAltFZJN772ZV/U7ywSKyGq9c62nAQuAJETnFlXy9DbgaeAF4FLgIr+TrVUCLqq4QkQ3ATcCfi0gFcD2wDlD3vx9W1ZxctV892EZROMgn3rkIgEBAqCgO09FrvvqR8LcQqS7JBzx3VWE4OJVTOmGIxoeLRl7IAuFGbhnV0lDVI6r6snvdAewEFqXpcglwn6pGVHUvUAesF5EFQJmqPq+ez+Ju4NKkPhvd6weB85wVciGwSVWbnVBswhOanNDQ0cu8soJB2SpF4ZBlBaWhOxojGBDmFoW9vy2ukTVSuafygmIxDSOnjCmm4dxGZwAvuqYvisg2EblTRMpd2yLgYFK3ete2yL0e2j6oj6rGgDagMs1YOaGhI0J1af6gtpL8kG3El4auSJyicJCSfM+6sK1Eskc0rsMC4bY1upFrMhYNESkBfgZ8WVXb8VxNy4G1wBHgu/6hKbprmvbx9kme29UiskVEtjQ2Zi9I3dDeS80Q0SjKD1oqaRp6op5oFIU9z6dtWpg9orH4sJTbUFCI2YpwI4dkJBoikocnGPeq6s8BVPWYqsZVNQH8CFjvDq8HFid1rwUOu/baFO2D+ohICJgDNKcZaxCqeruqrlPVddXV1Zl8pIxo6IhQU1owqK0439xT6ejui1MUDlHsLI0uWxWeNVIFwr2t0c3SMHJHJtlTAtwB7FTV7yW1L0g67OPAa+71w8AGlxF1MrASeElVjwAdInK2G/MK4KGkPn5m1GXAUy7u8ThwgYiUO/fXBa5t0umMxOiOxqkpG2xpFIeD5p5KQ080RmHegKVhAps9UgXCzT1l5JpMsqfOBT4DbBeRV1zb14BPichaPHfRPuDzAKq6Q0QeAF7Hy7y61mVOAVwD3AUU4mVNPeba7wDuEZE6PAtjgxurWUS+CWx2x92gqs3j+6hjo6G9F2C4eyocMvdUGgZiGiYa2SZ1INw2LDRyy6iioarPkjq28GiaPjcCN6Zo3wKsSdHeC1w+wlh3AneONs9s09DhrfqeVzbYPVWSH6I7GkdVB2VVGR7dfXHKCkIUuTRb27Qwe0RjKVaEBwPEE0oioQQC9n00Jh9bET4CvmikCoTHE9pfr9kYTE80RnE4RLFvaZhVljUisQQFQxf3uR1v+2xVuJEjTDRGYMA9NdzSAHO7jES3y57KDwUIiJ2nbBKJJSjIG7xQ0ndXmYvKyBUmGiPQ0BEhHApQVjjYgzcQ4DW3Syp6onEKw0FEhDmFebT19E31lE4Yevvi5A9zTzlLwyxfI0eYaIyAv0ZjaNyipD+V1O6gU9EVjfXHM8qLwzR3Rad4RicOkViC/Lzh24gAlkFl5AwTjRHw1mjkD2u3VNKRSSSU3r4Ehe4cVRaHaeo00cgGsXiCeEIpCA12T+UFvJ+wrdUwcoWJxgg0dESGZU4BtmgtDX5yQIG7G64wSyNr9LpzO9zS8CxhWxVu5AoTjRFItYUIMJAVZJbGMCIxT0j9u+GK4nwTjSwRcbsH5w+1NILmnjJyi4lGCnr74rT3xoZtVghQbO6pEYkOuRuuLA7T0h0lkbC74Iky1Irz8UXD3FNGrjDRSIF/d1xZMrKlYYvWhuNf2PL7LQ2vEJNlUE2c3hEsDUu5NXKNiUYK+kWjODzsPT8zyLb8Ho7vnvL3R6os8c5fk7moJsyAII+wuM8sDSNHmGik4Hintxrcv+glkx8KEAyIbfmdAr+0q39hq3Cia3GNiRMZKRBuMQ0jx5hopMC/yFUUD3dPiYjtdDsCQ++Gy4t80YhM2ZxOFHz31LCUW3NPGTnGRCMFA6Ix3NIAq6kxEr57yve7+5Zac5fFNCbKSJZGf0zDVoQbOcJEIwVNXVHygkJZQepNgIvzbXv0VAzNnhpwT5mlMVFGTLkNWUzDyC0mGilo7oxSURwecetzc0+lZqh7Kj/k1dWwQPjEGSnlNuRWhPdZWrORI0w0UtDUFUkZz/Ax91RqUmX42Krw7DBqyq25pzJGVfnjgRa84qDGWMmk3OtiEXlaRHaKyA4R+ZJrrxCRTSKy2z2XJ/W5TkTqRGSXiFyY1H6miGx3793iyr7iSsPe79pfFJGlSX2udP9jt4hcSQ5o6oqmTLf18ar3maUxlFQuFBON7DBSyq25p8bOC2818/Fbn2PzvpapnsqMJBNLIwb8g6quAs4GrhWR1cBXgSdVdSXwpPsb994G4DTgIuBWEfGvIrcBV+PVDV/p3ge4CmhR1RXAzcBNbqwK4HrgLGA9cH2yOE0WzV3REYPg4O0/ZSm3w0l1Yas00cgKA4Fw20Zkorx+pB2AuobOKZ7JzGRU0VDVI6r6snvdAewEFgGXABvdYRuBS93rS4D7VDWiqnuBOmC9iCwAylT1efXswruH9PHHehA4z1khFwKbVLVZVVuATQwIzaThxzRGwtxTqYkOWREOZmlkiwH31BBLo3+XW3O1ZEpdQwcAB1u6p3gmM5MxxTSc2+gM4EVgnqoeAU9YgBp32CLgYFK3ete2yL0e2j6oj6rGgDagMs1Yk0YkFqcjEqMqxcI+HwuEpyZVWmhFSZimrqj5jyfIaO6pmFkaGeNbGAeaTTTGQ8aiISIlwM+AL6tqe7pDU7Rpmvbx9kme29UiskVEtjQ2NqaZ2uikW9jnU5wfoqcvTtwyVgbRv41IcLB7KhpL2LYrEyQS86r2Dc3oM/fU2FBV3jzmicZBE41xkZFoiEgenmDcq6o/d83HnMsJ99zg2uuBxUnda4HDrr02RfugPiISAuYAzWnGGoSq3q6q61R1XXV1dSYfaUT8okFp3VNhf9NCuxAmE4klyAsKgcDAhW1xeREA+5vsBzoRIn2JYVYGQMida3NPZcbxzihtPX2EgwGzNMZJJtlTAtwB7FTV7yW99TDgZzNdCTyU1L7BZUSdjBfwfsm5sDpE5Gw35hVD+vhjXQY85eIejwMXiEi5C4Bf4NomjYEdbtNkT/mFmMxFNQjvwjY4ULu8pgSAPY0WdJwIkVh8WBAcvG1twsGAWRoZ4rumzlpWQWt3n+3APA5SL3kezLnAZ4DtIvKKa/sa8C3gARG5CjgAXA6gqjtE5AHgdbzMq2tV1b+6XgPcBRQCj7kHeKJ0j4jU4VkYG9xYzSLyTWCzO+4GVW0e52fNiNG2EIGBnW794KTh4btQkjmpsoiAwJ7Grima1YnBSJYGeDvd2jqNzPCD4B86tYbf7z7OweZu5iyaM8WzmlmMKhqq+iypYwsA543Q50bgxhTtW4A1Kdp7caKT4r07gTtHm2e28FcvV6WJaYSDnmhE7Ic6iGhs+IUtPxSktryIt8zSmBCRWIKCFJYGeHGNmMXXMmJ3Qycl+SHetbQC8OIaa0w0xoStCB9CU2eEUEAoKxxZT/0Lox/4NTwisURKF8ry6mKzNCZIb99wK84nLxiwyn0ZUtfQyfKaEpZUerE2i2uMHRONIbT39lFWmDfivlMwkFJqlsZgIrH4oMwpn2XVJew93mllXydAJIUV5xM291TGtHT3UVOaT1lBHnOL8kw0xoGJxhC6I/H+mMVI+MHeqP1QB+FZGsO/UsurS+jtS3C4rWcKZnViEInFR3RPhSwQnjHd0RjF7ve9pKLIRGMcmGgMoTMSoyQ/fajH3FOpGSlYu6y6GLBg+EToTRMIzwuKFWHKkK5InEKXMj+/rIDGDtu2f6yYaAyhKxqjeBTR8GtgR/rs7i6ZaHx4yi14lgZgwfAJ4GWmjRwIN0sjM5ItjfIi2+JmPJhoDKEzEh9VNAYsDfuhJpMq5RagqiTM3KI8HnrlsK0MHyde9tQIMY2QiUYmJBJKdzROkft9lxeHae3usy1uxoiJxhC6IwN3IiPhZwiZe2owkb5EvxWWjIhwwyVr2H6ojU//+EW7wI0DL3sqnaVhF77R6HHrqgYsjTyi8YSVORgjJhpD6IqM7p7y76YtED6YdBk+f/qOhXzjT0/j1YOt7DicbusyIxUjJRmAt5WIpdyOjl+iOdnSAGgxF9WYMNEYwtgC4fZDTSad3x3gPSuqAHjzaEeupnTCkG5FeDgUsBuYDOiJDrU0nGh0m2iMBRONJFSVrmic4vz07qmwiUZK0t0NAyyuKKIgL8AbJhpjQlXpTZNymx8K2ncxA/y94opc9lRFcR6ABcPHiIlGEpFYgnhCR8+eCvrZU+YLTSbVNiLJBAPCyppS3jxmojEW+uKK6vBaGj4FeQH7LmaAvyu1f1PoWxqt3bZp4Vgw0UjCr8bnb30+EiJCfihAxPzIg/BiGumttLfNL2WXicaY8BMuRjq3hXnB/iCvMTJ+wNu3NHzRMEtjbJhoJOGbr6NZGuDd9dk6jQFicc9KS5U9lczb5pXS2BGxH+oY8F1PI6XcFuQFbcflDOh2N4X+jg9lhXkEBFotpjEmTDSS8NcQlIwS0wAv7db8yAOMVI50KKfMLwVgl8U1MmagPnjq72VBXoBeu4EZla7+QLh3UxgMCHMK82g20RgTJhpJDPg8R7c0wsGArdNIIlPRONWJhsU1MidV7fVkCvKC9MbitkhtFLr7U24HxLe8OExLl8U0xoKJRhKd/eZrBu6pvIBZGklE+y9s6a20mtJ85hTmWQbVGPDdoCNbGkFULZtvNPrdz0m/7/KisKXcjpFMyr3eKSINIvJaUts/i8ghEXnFPS5Oeu86EakTkV0icmFS+5kist29d4sr+YorC3u/a39RRJYm9blSRHa7h18OdtLwv1SjrdMA7wdsufEDDARr03+lRISVNSXsabB9qDKl1z+3aSwNsL3QRqM7GkNkcGzI9p8aO5lYGncBF6Vov1lV17rHowAishqvVOtprs+tIuLfHt0GXI1XM3xl0phXAS2qugK4GbjJjVUBXA+cBawHrnd1wieN/uypTGIaIbM0kvHPxWiBcIAVNSXU2eaFGTNgaYyccgsD4mKkpisSpzgcGlQrp7woz1Jux8iov3BV/R1e3e5MuAS4T1UjqroXqAPWi8gCoExVn1fP8Xo3cGlSn43u9YPAec4KuRDYpKrNqtoCbCK1eGUNf5uBzCwNy41PZjQXSjIrakpo7oraHV6G+FbcSIv7CkJWsz4Tevpiw2rlVBSHae6OWjxoDEwkpvFFEdnm3Fe+BbAIOJh0TL1rW+ReD20f1EdVY0AbUJlmrEmjawwxjbBZGoPI1D0FsLzG2yq9zlxUGdE7iqVR6C6EtlYjPV0pdrAuLw4TjSXs3I2B8YrGbcByYC1wBPiua09VI1XTtI+3zyBE5GoR2SIiWxobG9PNOy2dEa9caSYuFtu6YTDRDLOnAFZUm2iMhdEW9/W7pyymkZbuaIzCIdZaeZFtJTJWxiUaqnpMVeOqmgB+hBdzAM8aWJx0aC1w2LXXpmgf1EdEQsAcPHfYSGOlms/tqrpOVddVV1eP5yMB/g63o7tXwAtKRs2H3E8kw+wpgEVzCynMC5poZMioMQ1zT2WEZ2kMFQ1/p1uLa2TKuETDxSh8Pg74mVUPAxtcRtTJeAHvl1T1CNAhIme7eMUVwENJffzMqMuAp1zc43HgAhEpd+6vC1zbpJFJ1T4fC4QPZizuqUBAWFZdbMHwDOlfXzBCnRdfqE000tMdjQ1zPfvbo9sCv8wZ9QopIj8FPgBUiUg9XkbTB0RkLZ67aB/weQBV3SEiDwCvAzHgWlX1v8nX4GViFQKPuQfAHcA9IlKHZ2FscGM1i8g3gc3uuBtUNdOA/LjoisRG3XfKx9xTgxlL9hR4wfAt+1p4+UALZQUhVtSUTub0ZjTdfYP3TBpKoYlGRnRF4ywqHyy8lb5odFmt8EwZ9Qqpqp9K0XxHmuNvBG5M0b4FWJOivRe4fISx7gTuHG2O2SKV+ToSlj01mNFcKENZUV3CQ68c5hO3Psep80v59ZffN5nTm9H0ROPD1hckYzGNzOiODLc0qkrzAWjsMNHIFFsRnkRnBlX7fPJDAauWlsRowdqhrF0yF/AsjjeOdvCWuapGpDsapzAvOGh9QTIFZmlkRFc0PqyUc2l+iPxQgOOd5p7KFBONJLqjo1ft8/FjGpbf7THa/khDee/Kap776oe4+6+8HIrHXjs6aXOb6XRH4yPGM8DcU5nSHY31l3r1ERGqSvLN0hgDJhpJdEXiGa3RAC/4qOoVyDEy37AwmYVzC1k4t5AzlszlsdeOTNbUZjy9ffH+tRip8C2NHnNPjUg0lqAvrsMsDYDq0nyOd5poZIqJRhJeffDMYxqA7XTr6A+EB8f+lbp4zQJeO9TOwebubE/rhKA7GqMob+SbGf+7aJbGyAxkoA0/j9WlZmmMBRMNh6q6dRqZWRpWJ3wwkViccCgwot89HecsrwRgx+G2bE/rhKA7mt7SCASEcChge0+lobu/at/w81hVYpbHcnkjAAAgAElEQVTGWDDRcETjCWIZ1Af38e/ubKdbj0hfgoIxuKaSWVxRBMDB5p5sTumEoWeUmAZ4cQ3b5XZkBmpppLY0mrqixCyxJSNMNBxj2RYdBrKEzNLwGIuVNpQ5hXmUFYQ4YO6plIwWCAcv7bYnapbGSAzU0kgd01C1BX6ZYqLhKMwL8t3L38G5K6oyOt5iGoPJ5MKWjiWVRRxsMdFIRU9fnMJREjT86n1GarrSxTRKvAV+FtfIjPHdGp6AFIaDfPLM2tEPdPippeYS8BjLFiypWFxexC4rAZsSb6O99Pd3BaGgBcLT4FsaqW5sqt0CP1urkRlmaYyTcNDcU8l0RyZoaVQUUd/SQyJhKcxD8ay4USyNcNBWhKehvcfbkHBOYd6w96pKbFX4WDDRGCe+pWGBcI+uaOb7dqWitqKIaCxBg/1wh9EzSvYUQEEoYDUh0tDmRGNu0ciiYRlUmWGiMU4spjGYnmg8ZWZKpizxM6gsrjGIaMzL6isaZcv5gryg7YWWhlYnGqUFw0WjOD9EcTholkaGmGiME8ueGkxXNDbqhS0di8sLATjQZKKRjG89jGpp5AXMPZWGtu4oZQUhgoHU64iqbIFfxphojBOzNAbTHYlTlOFq+lQsKi9ExCyNofRE02+L7lNo2VNpae3pY64ruJSKalvglzEmGuOkf0W43d15q+knGNPIDwWZX1ZgazWGMFoBJp+CvKCt00hDW09fyniGz7w5BdQ1dFqMMgNMNMZJ/4pwW0VKJJYgoUzI0gBvZbi5pwbjb38xunvKUm7T0drdlzJzyufyM2tp6Ijw05cOAN6N0C9fPcw//eI1/u3pulxNc0YwqmiIyJ0i0iAiryW1VYjIJhHZ7Z7Lk967TkTqRGSXiFyY1H6miGx3793iyr7iSsPe79pfFJGlSX2udP9jt4j4JWGnBX6JTbM0vNXgwIQsDYBV80vZeaSduKXd9tPTN/L6gmTy8wL02l3yiLT1pBeN959SzdnLKvjXp3bTFYnxzJuN/O1P/8i9L+7nu7/ZZangSWRiadwFXDSk7avAk6q6EnjS/Y2IrMYr13qa63OriPjf9tuAq/Hqhq9MGvMqoEVVVwA3Aze5sSrwSsueBawHrk8Wp6nGYhoDpNsMbiycXjuXrmicvcetIJNPpue2MC9INJawi9sItHZH07qnRISvXHQqxzuj/OODr/KtR99gaWUR131kFQkdyL4yMhANVf0dXu3uZC4BNrrXG4FLk9rvU9WIqu4F6oD1IrIAKFPV59WrWnT3kD7+WA8C5zkr5EJgk6o2q2oLsInh4jVlhAJCQCx7CgYubBNZEQ7w9to5AGyrt91ufXpcTKMwzdbokFS9z25ihpFIqBfTKBw5EA7wziXlfP2jq3h0+1F2HevgHy88lXlzCgBosiB5P+P9lc9T1SMAqnpERGpc+yLghaTj6l1bn3s9tN3vc9CNFRORNqAyuT1FnylHRMgPBS1wRvK+PhOzNJZXl1CYF2T7oTY+8c7Mt3Q5kck4ptFfUyNBmiShWUlnNEZCU68GH8rn3ruMeEJ581gnF58+n+f3NAHQ1BVl5WRPdIaQ7b2nUiVBa5r28fYZ/E9FrsZzfbFkyZLRZ5klwq7k62ynO5JZWuhoBAPCaQvL2G6WRj8Zu6fCVvJ1JNq63RYiadxTyXz+/cv7X1e4zQybbF+qfsabPXXMuZxwzw2uvR5YnHRcLXDYtdemaB/UR0RCwBw8d9hIYw1DVW9X1XWquq66unqcH2nseHXC7UeaLUsD4PTaOew43G61DRw9Y8ieAhONVLQ60ZibgaUxlMpib4uRpi5zT/mMVzQeBvxspiuBh5LaN7iMqJPxAt4vOVdWh4ic7eIVVwzp4491GfCUi3s8DlwgIuUuAH6Ba5s25OcFLHuKgbUEE41pgBfX6OmLs6exa8JjnQj0Z0+Nstre36HA9p8aTmuPZyWkW9w3EuXOOjFLY4BRf+Ui8lPgA0CViNTjZTR9C3hARK4CDgCXA6jqDhF5AHgdiAHXqqr/Lb4GLxOrEHjMPQDuAO4RkTo8C2ODG6tZRL4JbHbH3aCqQwPyU0pxOESHSzedzaQrcDNWVi0oA+DNYx28bX7phMeb6XRH44SDAUKj1F4vyBuIaRiDSbdZ4WiEggHKi/LM0khiVNFQ1U+N8NZ5Ixx/I3BjivYtwJoU7b040Unx3p3AnaPNcaqw2sIe6UppjpUFZd4eVMfaeyc81olATzQ2qmsKvJRbwDYtTIHvnsokEJ6KypJ8mrvM0vCxFeEToKokbKJBUobPBDYs9CkrDFGYF+RIm4kGZF4R0Y9pmHtqOG1pamlkQkVx2Ao0JWGiMQGqSvI53mFfpu5onMK84Ig7iI4FEWH+nAKOmqUBQHff6LU0AMrcBbG91xahDaW1O0pBXqBfWMdKVUnY1mkkYaIxASpL8unpi/e7Z2YrXZFYVjKnfOaXFXDULA3A1SnJ4NxWFFtq6EhksrAvHRXFYXNPJWGiMQGqXA73bLc2uqMT2xZ9KPPnmGj4dEdjFI2yGhzorxXR0j27v4upGG2zwtGoLM6npbvP0sAdJhoToMoVpG+c5aZrV2Ri26IPZf6cAo6199o+SmRW6hU8t155UZjmLnNPDaW1py/jhX2p8G8Om02QARONCVFttYWBzIO1mTK/rIBYQmkyl0B/vCgTKorzaLFzNoy2CVoaFW6Bn7moPEw0JoAVpPfoisaysrDPZ77bJM5cVGMTZM/SsAvbUI609TC/rGDc/SttK5FBmGhMgEqLaQCZB2szxf+BWwYVtPf09WdGjUZlSdhcKENo7+2jvTdGratBPx4qXZLBbL859DHRmAB5wQBzi/Jm/ZdpoqVeh7Kg39LoydqYM5G2nj46IjEWzs3sLrm8KGzuqSEcavG+Q4smIhol5p5KxkRjgtiqcG+X22xmT1WW5BMMyKy3NPoveHOLMjq+ojhMS3fUEgiS8M9hbXlm5zAVcwvzCAWExo7Z/Tv3MdGYILYq3LM0JrotejLBgFBTmj/rV4Ufah3bXXJ5UZiE2gK/ZOpbvJrzi+aO39IIBISqknwaTDQAE40J41kas9dsjSfUFf7JnqUBA2m3s5nDvmhkeMHrX+BnbpR+DrX2kB8K9KfNjpd5ZSYaPiYaE8TbSmT2fpn6t0XPoqUBsHBuIfuOd2d1zJnGodYewmO44JU70bC4xgD1LT0sKi/Eq8gwfqpLC2iY5TcxPiYaE6SqJExHJDZri99MdAfRkVhbO5dDrT2z+od6qKWHRXMzv+D5WT4WsB3gUGvPhFxTPmZpDGCiMUH8tRqz1SXg+90XZuGHmcy6peUAbNnfktVxZxL1Y7zg9VsalnbbT31Lz4SC4D41pQU0d0WJWnlnE42JUu22Epmtd8SH+0Vj/IunUnHawjnkhwJs2Td7RcO3NDKloshiGsl0R2M0d0UntEbDZ16ZbRnkMyHREJF9IrJdRF4RkS2urUJENonIbvdcnnT8dSJSJyK7ROTCpPYz3Th1InKLKwmLKxt7v2t/UUSWTmS+k8E8txBttpqufobTgjnZtTTCoQBrF89ly/5pVawxZ/T2xTneGRnT+oLCcJCCvIDFNBz+DU02RKPGicZsT86A7FgaH1TVtaq6zv39VeBJVV0JPOn+RkRW45VyPQ24CLhVRPyUm9uAq/Fqiq907wNcBbSo6grgZuCmLMw3q/hfptlqaRxq7aGiOJzRpnpjZd3ScnYcbp+VW8+PNXPKp8I2LeznYMv4zmEqakrdzWH77Lw5TGYy3FOXABvd643ApUnt96lqRFX3AnXAehFZAJSp6vOqqsDdQ/r4Yz0InCcTTYPIMpXF3kK0Y7P0y3S4tSfrrimfdUsriCeUVw60Tsr405mxrtHwKXcL/AzYd7wLgCUVWYhp+O6pjtl5c5jMREVDgd+IyFYRudq1zVPVIwDuuca1LwIOJvWtd22L3Ouh7YP6qGoMaAMqJzjnrBIMCNUl+bPWbD3S2svCLLumfN5ROxeA14+0T8r405lD47xLrigOW0zD8caRDiqKw/1xx4kw228Ok5locv25qnpYRGqATSLyRppjU1kImqY9XZ/BA3uCdTXAkiVL0s94EphXls+xWRrTONzawznLJ0fHK4rDVBSH2dPYOSnjT2derW+jND805qy0mtIC3jjaOEmzmlnsPNrOqgWlE16jAd7NYVVJmAazNCZmaajqYffcAPw3sB445lxOuOcGd3g9sDipey1w2LXXpmgf1EdEQsAcYFhkVFVvV9V1qrquurp6Ih9pXNSUzc6FP+29Y9tQbzwsry5mT0PXpI0/XXlxbxPvOrlizHXXV84robEjQlv37I5rxBPKrqMdrJpflrUx55UVmKXBBERDRIpFpNR/DVwAvAY8DFzpDrsSeMi9fhjY4DKiTsYLeL/kXFgdInK2i1dcMaSPP9ZlwFMu7jGtmFc2O91ThydpjUYyy6tLZp2l0dDRy1uNXZx1csWY+54yrwSANxs6sj2tGcXe411EYglOXZA90agptQV+MDFLYx7wrIi8CrwE/EpVfw18CzhfRHYD57u/UdUdwAPA68CvgWtV1V9GfQ3wY7zg+B7gMdd+B1ApInXA3+MysaYb80oLaOnuIxKbXavCj7R6QjmZorGipoSmruisSiN9aa9nTJ+1bOxuv5U1pQC8eWx2i8ZOFwdbtaA0a2POVo/CUMYd01DVt4B3pGhvAs4boc+NwI0p2rcAa1K09wKXj3eOuaJ/rUZ7hMVZyNSYKfSvBp+kQDh4lgbAnsZO1hWP/c57JtETjbOnsZMX32qmOBxkzcKx3yUvmltIUTjI7mOzyzobyhtH2wkFhBU1JVkbc35ZAU1dUXr74hRkWIL3RMRWhGeB/rUasyxIdri1h1BAspKdMhLJonGic9szdXzsX5/lJy/u58ylFYSCY/95BgLCypoSszSOdLC8uoT8UPYu7r4A1TWc+N/FdJhoZAF/4c9sC5LVNXRSW1445mDtWFhUXkh+KDArfqjPvNnISZVFrDupnD9ft3j0DiOwcl4pb85iSyORULYfasuqawpglYuPTNcU8FyFe000ssC8WbjFgKqyeV8z65ZOrssoGBBOripmT+OJnUHV2h1l+6E2PnFGLf/1N+/mo29fMO6xTplXwvHOyKyJA73wVhN//8Ar/THFLftbaOyI8IG31YzSc2wsqSiiMC/IG0empxX3Dw+8yl/fvWXS/4+JRhYoLwqTF5xdC3/qGjpp6e5j/TgyfMbK8poSdh2dnj/UbPH8niZU4T0rJ77mZeU87w579yywzgBufWYPP3/5EDdv2g3Aw68eoiAvwPmr52X1/wQDwinzS3nj6PSzNFSV5/Y0kR+a/Eu6iUYWCASEmllWpOVFP8MnB6Jx1skVHGrtOaHjGs/WHackP8Tb3Sr4ibDauVFmw2aPjR0Rnt3dyNyiPG7/3R4e236EX207wvmr51Ocn93CYACr5pey80h7zlxBmXKotYej7b28a5ItfzDRyBoL5hRwsGX2VJp7aW8z88rys7Kvz2ict8q7Y3zi9WOT/r+mij/UHefsZRXkjSP4PZR5ZQW8vXYOj792NAszm948su0wCYW7Prue5dUlXHPvy7R09/Gn71g4Kf/v1PmltHT3Tbv1Gltd3Rm/Ds1kYqKRJU5dUMobRzqm3R3IZKCqvLS3mfUnV2Zli4bRWDS3kNULynhyZ8PoB89Adh5pZ19TN+8/JXu7GXxkzQJerW/rT4s+EVFVfvHHQ6xeUMbaxXN55H++h69/dBWXrl2Y1XOZjB8M3znFwfBILM5j24+QSHjXm837minJD3FqFlfAj4SJRpZYtaCMjkiM+pYT90fqc7DZM4XX5+CuxufDq+exZX/zCVnK9L+21BMOBvjY27N3d/yRNfMB+PUJbG384MndvFrfxob1XqZZfijI5967jO9vOIPwJPn2/YvyjsNTKxobn9vHNfe+zINbvb1et+xr4Ywlcyc1k9HHRCNLrJ7m6XjZ5A97jgNwzvKqnP3P81fNI6Hwmx0n1kUwGkvwi1cO8eHVNf3lWrPB0qpiTp1fyqPbj2RtzOnET17Yz/ef2M1lZ9bymbNPytn/nVOUx+mL5vDAloPE4lNT+jWeUDY+tx+Am594k2Ptvew61pGTeAaYaGSNt80vRQRen+I7kFzwbN1x5pXls7y6OGf/c82iMk6dX8qPfv9Wv0k+03n6jQb+5dGdNHdFufzM8a/LGInLzqxl6/4WnnMif6LwXN1xrn94Bx86tYZvfeL0nLhIk/nbD61gf1M3v3jl8OgHTwJP7DzGodYerjjnJI609fLB7zyDKrx3ZW5u4kw0skRROMTJVcVT7uucbBIJ5fk9TZy7oiqnP1YR4doPrmBPYxe/PgGsjf/8w14+e9dm7npuHytrSiblB/8XZ5/EwjkF3PTYGydMrG3r/mauufdlllcX84MNa8e1an6inL96HqsXlPGDJ9/krSnI6Nv43D4WzS3kf39sNZeuXciahXO493NnccaS3LiLTTSyyKoFZTPWPZXpZos7j7bT3BXl3By6pnwuPn0By6qKufWZupz/72zy9BsN3PDI61yweh7b/vkCfvN375uUi19BXpC/O/8UXq1v42cvH8r6+Lnm6V0NfOpHL1JelMcdV76L0oK8KZmHiPD1j62iuTPKBTf/jnte2J+z/32gqZvn9jSx4V2LCQUDfH/DGTzwN+dw7orc/R5NNLLI6gVl1Lf00N47s2oZ/PCp3Zx+/W/4j9/uGdH10xmJ8Y1f7uDmTW8C5PRL6hMMCP/jrCW8dqh9xmYFtXRF+ccHX+XU+WV8f8NaygryJtVi+8Q7a1l3Ujlf/8X2GW0FR2Jxvv7fr7Gsqpiff+HcKd8Y9N3Lq3jmHz/IOcsr+eYjr+dsDdHPXq5HBD55Zu3oB08SJhpZZO1ib2HWn//HCzz9xsxID31g80G+85s3qS7N5/889gZrb/gNf/rDZ3n14OC63Lc+Xcd//mEfT+xsYPWCMubPmbzCS+l470ovlfIPdTPLTx+LJ9h5pJ1/eug1Wrv7+O7l76AonP3FZ0MJBoRb/+KdlBXk8YV7X56y4O1EuX/zQQ619vC/PrqKiiwmDEyE6tJ8vvtn76AgFOCrP9tGNDa55zaRUB7cWs97VlRNajmC0TDRyCLvXl7Jtz/5diKxOJ+/Zyvb69umekojcrwzwuc2buYrP9vGuSsqefr/+QA/2LCWS89YRFNnlD/7j+f51TYv8+ZQaw93PLuXS9cu5OV/Op/7P3/2lM37lHklVJXkzyjRSCSUv/zPzXzkB7/nkW1H+MIHV7B6HNuej5ea0gL+v0vXsPd4F4/OkBTcaCzB83ua6IrEONLWw78+Vcf6kyt4zxRYuOmoKS3g+j85jc37WvjUj17gaNvk7Apx/+YDfPauzRxq7eGyKbQyAORECZD5rFu3TrdsmfxNu9LR0hXlo7f8nrgqReEQffEEn3hnLV/4wPJpsw//dT/fxs+2HuLvzj+Fz567dNC8mjojfP6erWzZ38JfvnspW/e3sOtYB0/9w/upLZ/6eiFfvu+PPFt3nM3/68M5z5xJJhpL8PKBFvKCAWLxBEfbe9le30Z3X5z5ZQVcfPoCVtSUcOeze7nhkdf50nkred8pVbxzSXnO551IKB+++bcUhYP88ovvmdLzloyql1hx5x/2EhDhz9Yt5lhHL3f8fi9vHe+iujSfeEKJ9MW596/P7rfmpxuPbDvMVx7cRlVJPg98/pyMLPFf/PEQtz2zh2MdvVSX5PO+U6r52sWrhq21+PnL9fz9A69yUmUR715eyfV/ctqkXEdEZKuqrhvtuMm3j7OAiFwE/AAIAj9W1W9N8ZTSUl4c5t8+/U7++u6t/b7XW57cjaryDxe8Lav/K5FQXj/STn4oQG15EYXh0b9MTZ0Rfv7yIT55Zi3XfGD5sPcrS/K596/P4qs/285dz+1j4ZwCvv3Jt08LwQAvnvKLVw6z61hHTlbApqKhvZdr7n25f/sGn/xQgOL8EC3dUb636U3mleXT1BnlvFNr+PKHV07ZxToQED73nmV87b+388TOhqxv5jce2nr6+PovXuOXrx6mqiSfhCq/cVvFLK8u5l8+fjoPvXKI7mic729Y219bZTrysbcvpLa8iE//6AU+/eMX+MnnzmJBmuJkrx9u5ysPbmNZdTEfPX0B9S2eNZ8fCvDu5VU8sfMYkVicSCzBr7Yd4exlFfzkqrOmJFtsKNPe0hCRIPAmXunYemAz8ClVfT3V8dPB0kjF5+/ZwvN7mnjuuvMoSbGRmqqydX8LJQUhTqkpJZBmZaeqVy/gD3VN/NfWg7zltg0vDge55IxFvGtpOSdXlXByVTFzCodnmNzy5G6+t+lNnvj797GiZuSaA6rKm8c6WVFTkpOVpplyuLWHd3/rKT6yZj7fvHQNO4+0s2pBGVUlA8WgVJVDrT1UFIcpCodo6Ojlrj/sY9Prxzh1QRnvWlrOKfNKedfSijF9tr54gp++dIDvP7Gbnmic//0nq5k/p4BwMEBVST7LqovJCwY43hnhZ1vr2dPYSWlBHl/4wHIqSyavWFUm9PbF+cgPfs/+pi6uOGcp562qYf3JFeSHgjS099IRiTGvrCDl9zObRGMJfvnqYb79+Bsc74zydx9eyefeuwzw9jRbUlHESZVF08YaGgsvvtXEVRu3UBgOctmZtRxt6+Xdyyv7Y3GPbDvMziMdvLSvid6+BL/+0nv7vxfX/XwbP33pIABF4SDF+SHCwQCLKwr54f9456Dv92SQqaUxE0TjHOCfVfVC9/d1AKr6f1IdP11F45WDrVz6b3/gqx85lb95/3KOd0Z4bk8Tm/c2097bx+uH2/u3sp5fVsDfnb+SM08q51BrL0fbeojGEpQW5LF6YRm3/+6t/u0D3lE7h0+fdRL5eQF+9+ZxHtl2mEhSQK6yOExVST49fXG6ozG6o3G6o3Hef0o1G/9q/ZSci2zwb0/X8Z3f7ML/+laX5nPtB5Zz3+aDNHRESKjS2t1HTWk+l6+rZeNz++mKxli/tII9jV0c7/Q2nHvbvFI++vYFtHRHiafIHBNgcUURJ1UW09wV4T9+9xZvNXZxzrJKvnHJaZwyL7uFfiabjt4+vvnI6/zX1npUvRoR7z+lmvs2H6Av7n3+onCQeWUF1JTmU1NWwNsXzeHDq+dxsLmb/U1dHGuP8Gp9K3nBAJesXcj+pm72NXVRkBfkeEeEvniCM08qZ25RmIQq8YSypKKIgrwg920+yNNvNNAZifGO2jl845I109blNF7ePNbB5zZuob6lm/KiME1Dtr5ZMMcT5hsuWcM5ywe2wo/E4nzzkddZUV3Cp85aktWqg5lwIonGZcBFqvo59/dngLNU9Yupjp+uogHw6R+/wB/qmqgoDvfvoVSSH6KyJExlcZgN65cgwE9fOsDLB1pHHEcEvvjBFVxxztJhpVYjsTgHm3vYe7yLvcc72Xu8i+OdUYrDQQrDIYrDQYryQ3zijEUsrcrdiu7J4MW3mvjd7kZOnV/G/318Fweau1leXczZyypJqBc0v3/zQd442sE5yyq58eNrWFZdgqpyrD3CC281ccuTu3nreBfF4WDK/YriCaW9N9b/97KqYr528SrOW1UzI++Efdp7+3h+TxPfeXwXuxs6+eQ7a3nPykoa2iMca4/Q0NFLQ3uEo+29HGgevHtzQOCUeaW0dvdxtL0XEVhQVkAklqCiOIwycknUOYV5XHz6fD68ah4ffFtNWot6JhNPKH3xBPmhAC8faGHnkQ56+7ybtZXT9EbjRBKNy4ELh4jGelX926RjrgauBliyZMmZ+/fnbrHNWGjuivLzl+t542gHJ1cVc+6KKtYsLBvmp1RVntnVSHtvHwvmFLJgTgEFeUFauqP88UALy6tLJr1i3kyjtTvKln0tfOBt1YPOZzSWYPuhNs5YPDflBSqRUCKxRNpYUGNHhCNtPRSFgyytLJ4WfuVs0RdP0NbTl9b1saexk+fqjrOsuoSVNSWUF4f7g/8vH2hlWXXxsP5tPX1E+uIEAoIAexq7aOr0qullEnczcs+JJBonhHvKMAxjOpOpaMyEW6bNwEoROVlEwsAG4OEpnpNhGMasZNqn3KpqTES+CDyOl3J7p6rumOJpGYZhzEqmvWgAqOqjwKNTPQ/DMIzZzkxwTxmGYRjTBBMNwzAMI2NMNAzDMIyMMdEwDMMwMsZEwzAMw8iYab+4b6yISAewK0vDVQHZLNwwB8hWkY1sjjUZ4/lk6xzOhM+b7TGn8/dvJoxn5y9zqoBiVa0e9UhVPaEewJbpOJYb7/bpONZkjJftczgTPu8kzHHafv9myHh2/ibhXJl7Krf8cpqONRnjZZuZ8HntHE6v8bLNdP+8OTl/J6J7aotmsH9Krseardg5HD927iaGnb/MGcu5OhEtjdun6VizFTuH48fO3cSw85c5GZ+rE87SMAzDMCaPE9HSMAzDMCaJWSUaIrJYRJ4WkZ0iskNEvuTaK0Rkk4jsds/lrr3SHd8pIj9MGqdURF5JehwXke9P1efKJdk6h+69T4nIdhHZJiK/FpGqqfhMuSLL5+7P3XnbISLfnorPk2vGcf7OF5Gt7ju2VUQ+lDTWma69TkRukZlchjHXZDPla7o/gAXAO93rUuBNYDXwbeCrrv2rwE3udTHwHuBvgB+mGXcr8L6p/nwz6Rzi7bDcAFS5v7+NV2xryj/jDDh3lcABoNr9vRE4b6o/3zQ8f2cAC93rNcChpLFeAs7BKwP/GPCRqf58M+UxqywNVT2iqi+71x3ATmARcAneDw/3fKk7pktVnwV6RxpTRFYCNcDvJ3Hq04YsnkNxj2J3l1cGHJ78TzB1ZPHcLQPeVNVG9/cTwCcnefpTzjjO3x9V1f9O7QAKRCRfRBYAZar6vHoKcrffxxidWSUayYjIUrw7kReBeap6BLwvJp4IZPajURUAAAOdSURBVMqngPvdl29WMZFzqKp9wDXAdjyxWA3cMYnTnVZM8PtXB5wqIktFJIR3wVs8ebOdfozj/H0S+KOqRvCEpj7pvXrXZmTArBQNESkBfgZ8WVXbJzjcBuCnE5/VzGKi51BE8vBE4wxgIbANuC6rk5ymTPTcqWoL3rm7H8/C3QfEsjnH6cxYz5+InAbcBHzeb0px2Ky76Rsvs0403MXqZ8C9qvpz13zMmay454YMx3oHEFLVrZMy2WlKls7hWgBV3eOstAeAd0/SlKcN2fr+qeovVfUsVT0Hb6+13ZM15+nEWM+fiNQC/w1coap7XHM9UJs0bC0nuGs0m8wq0XC+8zuAnar6vaS3HgaudK+vBB7KcMhPMcusjCyew0PAahHxN0g7H89HfcKSze+fiNS453LgC8CPszvb6cdYz5+IzAV+BVynqn/wD3YurA4ROduNeQWZ/+aNqY7E5/KBl4mieK6QV9zjYrxslCfx7taeBCqS+uwDmoFOvDuU1UnvvQWcOtWfa6aeQ7ysoJ1urF8ClVP9+WbQufsp8Lp7bJjqzzYdzx/wdaAr6dhXgBr33jrgNWAP8EPcQmd7jP6wFeGGYRhGxswq95RhGIYxMUw0DMMwjIwx0TAMwzAyxkTDMAzDyBgTDcMwDCNjTDQMI8eIyN+IyBVjOH6piLw2mXMyjEwJTfUEDGM2ISIhVf33qZ6HYYwXEw3DGCNus7xf422WdwbeFt1XAKuA7wElwHHgL1X1iIg8AzwHnAs8LCKlQKeqfkdE1gL/DhThLTT7K1VtEZEzgTuBbuDZ3H06w0iPuacMY3y8DbhdVd8OtAPXAv8KXKaq/gX/xqTj56rq+1X1u0PGuRv4f90424HrXft/Av9Tvb2lDGPaYJaGYYyPgzqwn9FPgK/hFfrZ5IrABYEjScffP3QAEZmDJya/dU0bgf9K0X4P8JHsfwTDGDsmGoYxPobuv9MB7EhjGXSNYWxJMb5hTAvMPWUY42OJiPgC8SngBaDabxORPFfHYURUtQ1oEZH3uqbPAL9V1VagTUTe49o/nf3pG8b4MEvDMMbHTuBKEfkPvN1V/xV4HLjFuZdCwPfxyoym40rg30WkCG/X5M+69s8Cd4pItxvXMKYFtsutYYwRlz31iKqumeKpGEbOMfeUYRiGkTFmaRiGYRgZY5aGYRiGkTEmGoZhGEbGmGgYhmEYGWOiYRiGYWSMiYZhGIaRMSYahmEYRsb8/wWoFjir5DRyAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -251,10 +2295,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", @@ -273,10 +2315,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 12, + "metadata": {}, "outputs": [], "source": [ "year = []\n", @@ -299,9 +2339,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": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAD8CAYAAAC7IukgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+w1fV95/HnC0GwiSgQMPxQYSPJBE0Xyx10x+5uxC2QHxOwNQ2rVWbqDIniju1mVqTJjqnSmZhpYsu6ISE1FbWK1MSRTaF4jTq1XQJcgr+Q0Hu7UkSIXOaiYmegXnjvH9/Pqd97cjmcc/lyftz7esycOd/7OZ/P53zul8N9n8+P7/ejiMDMzKxIwxrdADMzG3wcXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwjm4mJlZ4YY3ugH19JGPfCSmTp3a6GaYmbWU7du3H4qI8bWUGVLBZerUqXR0dDS6GWZmLUXSP9daxsNiZmZWOAcXMzMrnIOLmZkVrqrgImmPpFckvSipI6V9Q9KbKe1FSZ/N5V8uqUvSbknzcumzUj1dklZKUkofKenxlL5F0tRcmcWSOtNjcS59Wsrbmcqeffqnw8zMilBLz+XqiJgZEW25tPtS2syI2AAgaQawCLgUmA98V9JZKf8qYAkwPT3mp/SbgcMRcQlwH3BvqmsscBdwBTAbuEvSmFTm3vT+04HDqQ4zM2sCZ2JYbAGwNiKORcTrQBcwW9JEYHREbI5sh7KHgIW5MmvS8RPANalXMw9oj4ieiDgMtAPz02tzUl5S2VJdLe/gu0f53e9v5uCRo41uipnZgFQbXAJ4WtJ2SUty6bdJelnSD3M9isnAG7k8+1La5HRcnt6nTET0Au8A4yrUNQ54O+Utr6vlrfxpJ9v29LDymc5GN8XMbECqvc7lqojYL2kC0C7pF2RDXPeQBZ57gG8Dvw+on/JRIZ0BlKlUVx8pGC4BuOiii/rL0jQ+8fWNHOs98W8/P7JlL49s2cvI4cPYveIzDWyZmVltquq5RMT+9HwQeBKYHRFvRcTxiDgB/IBsTgSyXsSFueJTgP0pfUo/6X3KSBoOnAf0VKjrEHB+ylteV3nbV0dEW0S0jR9f0wWmdffCHVfzhZmTGDUi+2cZNWIYC2ZO4oVlVze4ZWZmtTllcJH0IUnnlo6BucCraQ6l5Frg1XS8HliUVoBNI5u43xoRB4Ajkq5McyY3AU/lypRWgl0HPJvmZTYBcyWNScNuc4FN6bXnUl5S2VJdLWvC6FGcO3I4x3pPMHL4MI71nuDckcOZcO6oRjfNzKwm1QyLXQA8mVYNDwcejYi/lfSwpJlkw1F7gC8DRMROSeuA14BeYGlEHE913QI8CJwDbEwPgAeAhyV1kfVYFqW6eiTdA2xL+e6OiJ50vAxYK2kFsCPV0fIOvXeMG664mOtnX8SjW/fS7Ul9M2tByjoBQ0NbW1v43mJmZrWRtL3sMpRT8hX6ZmZWOAcXMzMrnIOLmZkVzsHFzMwK5+BiZmaFc3BpUb7/mJk1MweXFuX7j5lZM6v23mLWJHz/MTNrBe65tBjff8zMWoGDS4vx/cfMrBV4WKwF+f5jZtbsfG+xKhx89yi3PbaD+6+/3D0EMxtyfG+xM8Qrs8zMauNhsQq8MsvMbGDcc6nAK7PMzAbGwaUCr8wyMxuYqoKLpD2SXpH0oqSOlDZWUrukzvQ8Jpd/uaQuSbslzculz0r1dElambY7Jm2J/HhK3yJpaq7M4vQenZIW59KnpbydqezZp386flVpZdaTt17FDVdcTPd7x87E25iZDSpVrRaTtAdoi4hDubRvAT0R8U1JdwJjImKZpBnAY8BsYBLwDPDxiDguaStwO/AzYAOwMiI2SroV+PWI+IqkRcC1EfElSWOBDqCNbDvl7cCsiDictlL+cUSslfQ94KWIWFXp9/BOlGZmtav3arEFwJp0vAZYmEtfGxHHIuJ1oAuYLWkiMDoiNkcW0R4qK1Oq6wngmtSrmQe0R0RPRBwG2oH56bU5KW/5+5uZWYNVG1wCeFrSdklLUtoFEXEAID1PSOmTgTdyZfeltMnpuDy9T5mI6AXeAcZVqGsc8HbKW16XmZk1WLVLka+KiP2SJgDtkn5RIa/6SYsK6QMpU6muvo3JguESgIsuuqi/LGZmVrCqei4RsT89HwSeJJtPeSsNdZGeD6bs+4ALc8WnAPtT+pR+0vuUkTQcOA/oqVDXIeD8lLe8rvK2r46ItohoGz9+fDW/rpmZnaZTBhdJH5J0bukYmAu8CqwHSqu3FgNPpeP1wKK0AmwaMB3YmobOjki6Ms2Z3FRWplTXdcCzaV5mEzBX0pi0Gm0usCm99lzKW/7+ZmbWYNUMi10APJlWDQ8HHo2Iv5W0DVgn6WZgL/BFgIjYmVZyvQb0Aksj4niq6xbgQeAcYGN6ADwAPCypi6zHsijV1SPpHmBbynd3RPSk42XAWkkrgB2pDjMzawK+caWZmVXkG1eamVlTcHAxM7PCObiYmVnhHFzMzKxwDi5mZlY4BxczMyucg4uZmRXOwcXMzArn4GJmZoVzcDEzs8I5uJiZWeEcXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhas6uEg6S9IOST9JP39D0puSXkyPz+byLpfUJWm3pHm59FmSXkmvrVTaO1nSSEmPp/QtkqbmyiyW1Jkei3Pp01LezlT27NM7FWZmVpRaei63A7vK0u6LiJnpsQFA0gxgEXApMB/4rqSzUv5VwBJgenrMT+k3A4cj4hLgPuDeVNdY4C7gCmA2cJekManMven9pwOHUx1mZtYEqgoukqYAnwP+oorsC4C1EXEsIl4HuoDZkiYCoyNic0QE8BCwMFdmTTp+Argm9WrmAe0R0RMRh4F2YH56bU7KSypbqsvMzBqs2p7LnwF3ACfK0m+T9LKkH+Z6FJOBN3J59qW0yem4PL1PmYjoBd4BxlWoaxzwdspbXpeZmTXYKYOLpM8DByNie9lLq4CPATOBA8C3S0X6qSYqpA+kTKW6+pC0RFKHpI7u7u7+spiZWcGq6blcBXxB0h5gLTBH0iMR8VZEHI+IE8APyOZEIOtFXJgrPwXYn9Kn9JPep4yk4cB5QE+Fug4B56e85XX1ERGrI6ItItrGjx9fxa9rZman65TBJSKWR8SUiJhKNlH/bET8XppDKbkWeDUdrwcWpRVg08gm7rdGxAHgiKQr05zJTcBTuTKllWDXpfcIYBMwV9KYNOw2F9iUXnsu5SWVLdVlZmYNNvzUWU7qW5Jmkg1H7QG+DBAROyWtA14DeoGlEXE8lbkFeBA4B9iYHgAPAA9L6iLrsSxKdfVIugfYlvLdHRE96XgZsFbSCmBHqsPMzJqAsk7A0NDW1hYdHR2NboaZWUuRtD0i2mop4yv0zcyscA4uZmZWOAcXMzMrnIOLmZkVzsHFzMwK5+BiZmaFc3AxM7PCObiYmVnhHFzMbMg7+O5Rfvf7mzl45GijmzJoOLiY2ZC38qedbNvTw8pnOhvdlEHjdO4tZmbW0j7x9Y0c6/1gm6pHtuzlkS17GTl8GLtXfKaBLWt97rmY2ZD1wh1X84WZkxg1IvtTOGrEMBbMnMQLy65ucMtan4OLmTWdes2BTBg9inNHDudY7wlGDh/Gsd4TnDtyOBPOHXVG33cocHAxs6ZTzzmQQ+8d44YrLubJW6/ihisupvu9Y2f8PYcC33LfzJpG+RxIiedAGsu33DezluY5kMHDwcXMmobnQAaPqoOLpLMk7ZD0k/TzWEntkjrT85hc3uWSuiTtljQvlz5L0ivptZWSlNJHSno8pW+RNDVXZnF6j05Ji3Pp01LezlT27NM7FWbWDDwHMjhUPeci6b8DbcDoiPi8pG8BPRHxTUl3AmMiYpmkGcBjwGxgEvAM8PGIOC5pK3A78DNgA7AyIjZKuhX49Yj4iqRFwLUR8SVJY4GO9L4BbAdmRcRhSeuAH0fEWknfA16KiFWVfgfPuZiZ1e6MzblImgJ8DviLXPICYE06XgMszKWvjYhjEfE60AXMljSRLDBtjiyiPVRWplTXE8A1qVczD2iPiJ6IOAy0A/PTa3NS3vL3NzOzBqt2WOzPgDuA/DKOCyLiAEB6npDSJwNv5PLtS2mT03F5ep8yEdELvAOMq1DXOODtlLe8rj4kLZHUIamju7u7yl/XzMxOxymDi6TPAwcjYnuVdaqftKiQPpAylerqmxixOiLaIqJt/Pjx/WUxM7OCVdNzuQr4gqQ9wFpgjqRHgLfSUBfp+WDKvw+4MFd+CrA/pU/pJ71PGUnDgfOAngp1HQLOT3nL6zKrie+Ia1a8UwaXiFgeEVMiYiqwCHg2In4PWA+UVm8tBp5Kx+uBRWkF2DRgOrA1DZ0dkXRlmjO5qaxMqa7r0nsEsAmYK2lMWo02F9iUXnsu5S1/f7Oa+I64ZsU7nbsifxNYJ+lmYC/wRYCI2JlWcr0G9AJLI+J4KnML8CBwDrAxPQAeAB6W1EXWY1mU6uqRdA+wLeW7OyJ60vEyYK2kFcCOVIdZ1XxHXLMzx7d/sSHr4LtHWbFhF0/v/CVH3z/BqBHDmHfpR/na5z7pi/bMcnz7F7Ma1HI1eKvOy7Rqu631ObjYkFbt1eCtOi/Tqu221udhMbMKWvUuva3abmtOHhYzK1ir3qW3Vdttg4eDi1kFrXqX3lZttw0ep7MU2WxIKM3LXD/7Ih7dupfuFpkcb9V22+DgORczM6vIcy5mZtYUHFzMzKxwDi5mZgXyhasZBxczswL5wtWMV4tZSzr47lFue2wH919/uZfXWlPwjVD7cs/FWpK/HVqz8YWrfbnnYi3F3w4/4N5bc/GFq32552Itxd8OP+DeW/Op9kaoQ4F7LtZ0Kn0j97dD996a2fdv/OA6wxULL2tgSxrvlD0XSaMkbZX0kqSdkv44pX9D0puSXkyPz+bKLJfUJWm3pHm59FmSXkmvrUzbHZO2RH48pW+RNDVXZrGkzvRYnEuflvJ2prJnF3NKbKCKWoJ5qm/kQ/3boXtv1gqq6bkcA+ZExHuSRgB/L6m0PfF9EfGn+cySZpBtU3wpMAl4RtLH01bHq4AlwM+ADcB8sq2ObwYOR8QlkhYB9wJfkjQWuAtoAwLYLml9RBxOee6LiLWSvpfqWDXwU2GnKx8UVlz7qZrLV/uNfKh/O3TvzVrBKXsukXkv/TgiPSrdkGwBsDYijkXE60AXMFvSRGB0RGyO7IZmDwELc2XWpOMngGtSr2Ye0B4RPSmgtAPz02tzUl5S2VJdVmef+PpGpt75NzyyZS8RWVCYeuff8Imvbzx14Rx/I6/eUO+9WfOras5F0lnAduAS4H9HxBZJnwFuk3QT0AF8NQWAyWQ9k5J9Ke39dFyeTnp+AyAieiW9A4zLp5eVGQe8HRG9/dRldfbCHVefdC/6WvgbefWGeu/Nml9Vq8Ui4nhEzASmkPVCLiMbgvoYMBM4AHw7ZVd/VVRIH0iZSnX1IWmJpA5JHd3d3f1lsdNUZFDwN/LBz7dHGRpqWi0WEW9Leh6Yn59rkfQD4Cfpx33AhbliU4D9KX1KP+n5MvskDQfOA3pS+qfLyjwPHALOlzQ89V7ydZW3eTWwGrJb7tfy+1r1ito7xN/IB7/TnZuz1nDK/VwkjQfeT4HlHOBpssn07RFxIOX5Q+CKiFgk6VLgUWA22YT+T4HpEXFc0jbgvwFbyCb0/1dEbJC0FPhURHwlTej/dkT8bprQ3w78RmrOz4FZEdEj6a+BH+Um9F+OiO9W+l28n4tZ45Qv2CjxEurmN5D9XKrpuUwE1qR5l2HAuoj4iaSHJc0kG47aA3wZICJ2SloHvAb0AkvTSjGAW4AHgXPIVomVZnwfAB6W1EXWY1mU6uqRdA+wLeW7OyJ60vEyYK2kFcCOVIeZNami5uasNZwyuETEy8Dl/aTfWKHMnwB/0k96B/ArYx0RcRT44knq+iHww37S/x9Z78jMWoAXbAwtvkLfzOqmqLk5a36nnHMZTDznYmZWu4HMufjGlWZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwjm4mJlZ4RxczMyscA4uZmZWOAcXszrzfiY2FDi4mNVZfj8Ts8HKN640q5Py/Uwe2bKXR7bs9X4mNii552JWJy/ccTVfmDmJUSOy/3ajRgxjwcxJvLDs6ga3zKx4Di5mdeL9TGwocXCxQa3ZJs9L+5k8eetV3HDFxXS/d6zRTTI7I04ZXCSNkrRV0kuSdkr645Q+VlK7pM70PCZXZrmkLkm7Jc3Lpc+S9Ep6baUkpfSRkh5P6VskTc2VWZzeo1PS4lz6tJS3M5U9u5hTYoNJs02ef//GNlYsvIwZk0azYuFlfP/GmrbIMGsZp9wsLAWAD0XEe5JGAH8P3A78NtATEd+UdCcwJiKWSZoBPEa2BfEk4Bng4xFxXNLWVPZnwAZgZURslHQr8OsR8RVJi4BrI+JLksYCHUAbEMB2YFZEHJa0DvhxRKyV9D3gpYhYVel38WZhQ0f55HmJJ8/NandGNguLzHvpxxHpEcACYE1KXwMsTMcLgLURcSwiXge6gNmSJgKjI2JzZBHtobIypbqeAK5JQW0e0B4RPRFxGGgH5qfX5qS85e9v5snzBmm2YUhrnKrmXCSdJelF4CDZH/stwAURcQAgPU9I2ScDb+SK70tpk9NxeXqfMhHRC7wDjKtQ1zjg7ZS3vK7yti+R1CGpo7u7u5pf1wYBT543RrMNQ1rjVHWdS0QcB2ZKOh94UtJlFbKrvyoqpA+kTKW6+iZGrAZWQzYs1l8eG5xKk+fXz76IR7fupdvfps8YX8Nj5Wq6iDIi3pb0PDAfeEvSxIg4kIa8DqZs+4ALc8WmAPtT+pR+0vNl9kkaDpwH9KT0T5eVeR44BJwvaXjqveTrMgPoM1m+YmGl70N2ul6442pWbNjF0zt/ydH3TzBqxDDmXfpRvva5Tza6adYg1awWG596LEg6B/gvwC+A9UBp9dZi4Kl0vB5YlFaATQOmA1vT0NkRSVemOZObysqU6roOeDbNy2wC5koak1ajzQU2pdeeS3nL39/M6szDkFaump7LRGCNpLPIgtG6iPiJpM3AOkk3A3uBLwJExM60kus1oBdYmobVAG4BHgTOATamB8ADwMOSush6LItSXT2S7gG2pXx3R0RPOl4GrJW0AtiR6jCzBvEwpOWdcinyYHImlyIffPcotz22g/uvv9zf1szqpFX/37Vau8/IUmSrjlfJmNVfq/6/a9V218I9l9Pki/XM6q9V/9+1arvdc2kAX6xnVn+t+v+uVds9EA4up6mWVTK+etlaWTN9flt1dVqrtnsgHFwKUO2dbofCOKsNXs32+W3VO0y3artr5TmXOmjVcVYz8OfXPOfStIbSOKsNPv78Nq9mGqos5+BSB4NhnLWZP8R2Zg2Gz+9g1WxDlXk13VvMBq7Vr17Of4hXXPupRjfH6qzVP7+DTSvcKNRzLlaRx9vNms/Bd4+e9EahZ6JH6TkXK5zH282aTysMVXpYzCpqhQ+x2VDU7EOVDi52Ss3+ITYbipp9vyLPuQxxrXZ3VjOrP8+5WM2aeSmjmbUuD4sNUa2wlLHVuBdo9oFqtjm+UNJzknZJ2inp9pT+DUlvSnoxPT6bK7NcUpek3ZLm5dJnSXolvbYybXdM2hL58ZS+RdLUXJnFkjrTY3EufVrK25nKnl3MKRkavAqseO4Fmn2gmp5LL/DViPi5pHOB7ZLa02v3RcSf5jNLmkG2TfGlwCTgGUkfT1sdrwKWAD8DNgDzybY6vhk4HBGXSFoE3At8SdJY4C6gDYj03usj4nDKc19ErJX0vVTHqoGfiqHFq8CK416g2a86Zc8lIg5ExM/T8RFgFzC5QpEFwNqIOBYRrwNdwGxJE4HREbE5slUEDwELc2XWpOMngGtSr2Ye0B4RPSmgtAPz02tzUl5S2VJdVqWhcnfWM829QLNfVdOcSxquuhzYAlwF3CbpJqCDrHdzmCzw/CxXbF9Kez8dl6eTnt8AiIheSe8A4/LpZWXGAW9HRG8/dVmVmn0pY6twL9DsV1W9WkzSh4EfAX8QEe+SDUF9DJgJHAC+XcraT/GokD6QMpXqKm/3Ekkdkjq6u7v7y2J22twLNOurqp6LpBFkgeWvIuLHABHxVu71HwA/ST/uAy7MFZ8C7E/pU/pJz5fZJ2k4cB7Qk9I/XVbmeeAQcL6k4an3kq+rj4hYDayG7DqXan5fs1q5F2jWVzWrxQQ8AOyKiO/k0ifmsl0LvJqO1wOL0gqwacB0YGtEHACOSLoy1XkT8FSuTGkl2HXAs2leZhMwV9IYSWOAucCm9NpzKS+pbKkuMzNrsGp6LlcBNwKvSHoxpf0R8F8lzSQbjtoDfBkgInZKWge8RrbSbGlaKQZwC/AgcA7ZKrGNKf0B4GFJXWQ9lkWprh5J9wDbUr67I6InHS8D1kpaAexIdZiZWRPw7V/MzKwi3/7FzMyagoOLmVkTavWtxR1czKzl/5ANRq1+OyHfuNLM+vwhW3HtpxrdnCFtsNxOyBP6ZkNY+R+yklb7QzaYHHz3KCs27OLpnb/k6PsnGDViGPMu/Shf+9wnB3TXhyLu1u0JfTOrie+L1nyKvp1Qo4bXPCxm1oTqtTeM74vWnIrYWrzRw2sOLmZNqJ5zIEX8IbNiFXE7oRfuuPqkw2v14OBi1kQa8W3T90UbnBrdK/Wci1kT8RxI8YbyMutG3q3bPRezJtLob5uD0VBeZt3IXqmDi1mT8RxIMRo9oT3U+ToXMxuUir5eZCjzdS5mZomHGBvLw2JmNmh5iLFxPCxmZmYVeVjMzMyawimDi6QLJT0naZeknZJuT+ljJbVL6kzPY3JllkvqkrRb0rxc+ixJr6TXVkpSSh8p6fGUvkXS1FyZxek9OiUtzqVPS3k7U9mzizklZmZ2uqrpufQCX42ITwJXAkslzQDuBH4aEdOBn6afSa8tAi4F5gPflXRWqmsVsASYnh7zU/rNwOGIuAS4D7g31TUWuAu4ApgN3JULYvcC96X3P5zqMDOzJnDK4BIRByLi5+n4CLALmAwsANakbGuAhel4AbA2Io5FxOtAFzBb0kRgdERsjmyi56GyMqW6ngCuSb2aeUB7RPRExGGgHZifXpuT8pa/v5mZNVhNcy5puOpyYAtwQUQcgCwAARNStsnAG7li+1La5HRcnt6nTET0Au8A4yrUNQ54O+Utr8vMzBqs6uAi6cPAj4A/iIh3K2XtJy0qpA+kTKW6+jZGWiKpQ1JHd3d3f1nMzKxgVQUXSSPIAstfRcSPU/JbaaiL9Hwwpe8DLswVnwLsT+lT+knvU0bScOA8oKdCXYeA81Pe8rr6iIjVEdEWEW3jx4+v5tc1M7PTVM1qMQEPALsi4ju5l9YDpdVbi4GncumL0gqwaWQT91vT0NkRSVemOm8qK1Oq6zrg2TQvswmYK2lMmsifC2xKrz2X8pa/v5mZNVg1V+hfBdwIvCLpxZT2R8A3gXWSbgb2Al8EiIidktYBr5GtNFsaEcdTuVuAB4FzgI3pAVnwelhSF1mPZVGqq0fSPcC2lO/uiOhJx8uAtZJWADtSHWZm1gR8hb6ZmVXkK/St6Q3ljZvMhhIHF6ur/MZNZjZ4+a7IVhfeuMlsaHHPxerCe8ObDS0OLlYX3rjJbGjxsJjVjTduMhs6vBTZzMwq8lJkMzNrCg4uZmZWOAeXQcwXLJpZozi4DGK+YNHMGsWrxQYhX7BoZo3mnssg5AsWzYrnYebaOLgMQr5g0ax4HmaujYfFBilfsGhWDA8zD4wvojQzq+Dgu0dZsWEXT+/8JUffP8GoEcOYd+lH+drnPjlkRgPOyEWUkn4o6aCkV3Np35D0pqQX0+OzudeWS+qStFvSvFz6LEmvpNdWpq2OSdshP57St0iamiuzWFJneizOpU9LeTtT2bNr+aXNzKrlYeaBqWbO5UFgfj/p90XEzPTYACBpBtkWxZemMt+VdFbKvwpYAkxPj1KdNwOHI+IS4D7g3lTXWOAu4ApgNnCXpDGpzL3p/acDh1MdZmZnRGmY+clbr+KGKy6m+71jjW5S0zvlnEtE/F2+N3EKC4C1EXEMeF1SFzBb0h5gdERsBpD0ELAQ2JjKfCOVfwK4P/Vq5gHtEdGTyrQD8yWtBeYA16cya1L5VVW20cysJt+/8YMRoRULL2tgS1rH6awWu03Sy2nYrNSjmAy8kcuzL6VNTsfl6X3KREQv8A4wrkJd44C3U97yuszMrAkMNLisAj4GzAQOAN9O6eonb1RIH0iZSnX9CklLJHVI6uju7j5ZNjMzK9CAgktEvBURxyPiBPADsjkRyHoRF+ayTgH2p/Qp/aT3KSNpOHAe0FOhrkPA+SlveV39tXV1RLRFRNv48eNr/VXNzGwABhRcJE3M/XgtUFpJth5YlFaATSObuN8aEQeAI5KuTPMpNwFP5cqUVoJdBzwb2froTcBcSWPSsNtcYFN67bmUl1S2VJeZmTWBU07oS3oM+DTwEUn7yFZwfVrSTLLhqD3AlwEiYqekdcBrQC+wNCKOp6puIVt5dg7ZRP7GlP4A8HCa/O8hW21GRPRIugfYlvLdXZrcB5YBayWtAHakOszMrEn4IkozM6toIBdRDqngIqkb+Od+XvoI2VxOq3G768vtrq9WbTe0bttP1u6LI6KmSeshFVxORlJHrVG5Gbjd9eV211erthtat+1Fttt3RTYzs8I5uJiZWeEcXDKrG92AAXK768vtrq9WbTe0btsLa7fnXMzMrHDuuZiZWeEGZXA5yR40/17S5rSnzP+RNDqlj5C0JqXvkrQ8V+b5tC9Nad+aCU3U7rMl/WVKf0nSp3Nl+t07pwXaXe/zfaGk59K/+05Jt6f0sZLa035B7bkbs9a8X1ELtLtu57zWdksal/K/J+n+srrq/Rkvsu3NfM5/S9L2dG63S5qTq6u2cx4Rg+4B/CfgN4BXc2nbgP+cjn8fuCcdX0+2TQDAr5HdcWBq+vl5oK1J270U+Mt0PAHYDgxLP28F/gPZTT43Ap9pkXbX+3xPBH4jHZ8L/CMwA/jPDIkDAAADj0lEQVQWcGdKvxO4Nx3PAF4CRgLTgH8Czqr3OS+43XU75wNo94eA3wS+AtxfVle9P+NFtr2Zz/nlwKR0fBnw5kDP+aDsuUTE35HdSibvE8DfpeN24HdK2YEPKbsR5jnAvwLv1qOd5Wps9wzgp6ncQeBtoE3Zfd9GR8TmyD4Rpb1zmrrdZ7J9JxMRByLi5+n4CLCLbPuGBWT7BJGeS+fv3/YriojXgdJ+RXU950W1+0y172RqbXdE/EtE/D1wNF9Pgz7jhbS93gbQ7h0RUboR8E5glLJ7RdZ8zgdlcDmJV4EvpOMv8sEdl58A/oVs64C9wJ/GB/cwA/jL1HX9n2e6630SJ2v3S8ACScOV3SR0Vnqt0t459VRru0sacr6VbYh3ObAFuCCym62SnkvDFgPZr+iMOs12l9T9nFfZ7pNp6Gf8NNte0grn/HeAHZFt/ljzOR9KweX3gaWStpN1D/81pc8GjgOTyIYMvirp36XXboiITwH/MT1urG+TgZO3+4dk/8AdwJ8B/5fsZqE17XdzBtXabmjQ+Zb0YeBHwB9ERKVeayF7DBWlgHZDA855De0+aRX9pNXlM15A26EFzrmkS8m2k/9yKamfbBXP+ZAJLhHxi4iYGxGzgMfIxp0hm3P524h4Pw3T/ANpmCYi3kzPR4BHacxQQr/tjojeiPjDiJgZEQuA84FOKu+d08ztbsj5ljSC7D/dX0XEj1PyW2kYoDQEczClD2S/omZud93PeY3tPpmGfMYLanvTn3NJU4AngZsiovR3suZzPmSCS2lFhqRhwNeB76WX9gJzlPkQcCXwizRs85FUZgTweT7Yt6bh7Zb0a6m9SPotoDciXovKe+c0bbsbcb7T+XkA2BUR38m9lN9jKL9f0ED2K2radtf7nA+g3f1qxGe8qLY3+zmXdD7wN8DyiPiHUuYBnfNKs/2t+iD7pnwAeJ8s4t4M3E62UuIfgW/ywQWkHwb+mmzy6jXgf8QHqz22Ay+n1/6ctMKmSdo9FdhNNkH3DNldS0v1tJF9YP8JuL9Uppnb3aDz/ZtkXfuXgRfT47PAOLJFB53peWyuzNfSed1NbrVMPc95Ue2u9zkfYLv3kC0WeS99tmY06DNeSNub/ZyTfRH8l1zeF4EJAznnvkLfzMwKN2SGxczMrH4cXMzMrHAOLmZmVjgHFzMzK5yDi5mZFc7BxczMCufgYmZmhXNwMTOzwv1/q4MdzlEclD8AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -315,9 +2378,54 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2012 2175217\n", + "2003 2234584\n", + "2019 2254386\n", + "2006 2307352\n", + "2017 2321583\n", + "2001 2529279\n", + "1992 2574578\n", + "1993 2703886\n", + "2018 2705325\n", + "1988 2765617\n", + "2007 2780164\n", + "1987 2855570\n", + "2016 2856393\n", + "2011 2857040\n", + "2008 2973918\n", + "1998 3034904\n", + "2002 3125418\n", + "2009 3444020\n", + "1994 3514763\n", + "1996 3539413\n", + "2004 3567744\n", + "1997 3620066\n", + "2015 3654892\n", + "2000 3826372\n", + "2005 3835025\n", + "1999 3908112\n", + "2010 4111392\n", + "2013 4182691\n", + "1986 5115251\n", + "1990 5235827\n", + "1989 5466192\n", + "dtype: int64" + ] + }, + "execution_count": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -332,9 +2440,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+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGalJREFUeJzt3X2UJXV95/H3h5kBhmkYjAONDsr4QAhIqzgXXWQ13WhcdNCcGE5QQcVIGo0PRCdnM8v6sLrLOj5MsphgkklUiAodw8PZyBjUE2hQNEgPqA2OEBdmlSEMAjLSMAuMfPePX7XctP1wq+69XcXPz+ucPn3vrbpVn/rdut9b9auqexURmJlZPvaqO4CZmfWWC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLzNJ+TnzVqlWxZs2aWYc9+OCDrFixop+zr6zJ2aDZ+Zytuibnc7bqyubbunXrPRFxUFczjYi+/a1duzbmctVVV805rG5NzhbR7HzOVl2T8zlbdWXzARPRZe11V4yZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWWmVGGX9B5JN0u6SdJFkvbtVzAzM6um48IuaTXwbqAVEUcDS4DX9SuYmZlVU7YrZimwXNJSYD/gzt5HMjOzbihK/OappLOAc4DdwFcj4tRZxhkFRgEGBwfXjo2NzTqtqakpBgYGqmTuuyZng97nm9yxq2fTGlwOO3d3Pv7Q6pU9m/dCftVe115yturK5hsZGdkaEa1u5tlxYZf0JOAS4BTgfuAfgIsj4vNzPafVasXExMSsw8bHxxkeHi6bd1E0ORv0Pt+aDVt6Nq31Q3vYNNn5N1Vs37iuZ/NeyK/a69pLzlZd2XySui7sZbpiXg7cHhE/iYhHgUuBF3czczMz670yhf1HwH+QtJ8kAS8DtvUnlpmZVdVxYY+I64CLgRuAyeK5m/uUy8zMKir1tb0R8UHgg33KYmZmPeArT83MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy03Fhl3SEpO+0/f1M0h/1M5yZmZXX8U/jRcQtwPMBJC0BdgCX9SmXmZlVVLUr5mXA/4mI/9vLMGZm1j1FRPknSZ8BboiIv5hl2CgwCjA4OLh2bGxs1mlMTU0xMDBQet6LocnZoPf5Jnfs6tm0BpfDzt2djz+0emXP5r2Q9nbr5TKXMd/yNnm9c7bqyuYbGRnZGhGtbuZZurBL2hu4E3hOROycb9xWqxUTExOzDhsfH2d4eLjUvBdLk7NB7/Ot2bClZ9NaP7SHTZMd9/CxfeO6ns17Ie3t1stlLmO+5W3yeuds1ZXNJ6nrwl6lK+aVpK31eYu6mZnVo0phfz1wUa+DmJlZb5Qq7JL2A34LuLQ/cczMrFudd4YCEfEQ8OQ+ZTEzsx7wladmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmSn703gHSrpY0g8kbZN0XL+CmZlZNaV+Gg84F7giIk6WtDewXx8ymZlZFzou7JIOAF4KnA4QEY8Aj/QnlpmZVaWI6GxE6fnAZuD7wPOArcBZEfHgjPFGgVGAwcHBtWNjY7NOb2pqioGBgerJ+6jJ2aD3+SZ37OrZtAaXw87dnY8/tHplz+a9kPZ26+UylzHf8jZ5vXO26srmGxkZ2RoRrW7mWaawt4B/AY6PiOsknQv8LCLeP9dzWq1WTExMzDpsfHyc4eHh8okXQZOzQe/zrdmwpWfTWj+0h02Tnffwbd+4rmfzXkh7u/VymcuYb3mbvN45W3Vl80nqurCXOXh6B3BHRFxX3L8YeEE3Mzczs97ruLBHxF3AjyUdUTz0MlK3jJmZNUjZs2LeBXyhOCPmNuAtvY9kZmbdKFXYI+I7QFd9P2Zm1l++8tTMLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsM6V+QUnSduAB4OfAnm5/SdvMzHqv7G+eAoxExD09T2JmZj3hrhgzs8woIjofWbod+CkQwF9HxOZZxhkFRgEGBwfXjo2NzTqtqakpBgYGqmTuq8kduxhcDjt3L/68h1av7Gi8Xrfd5I5dPZtW2bbrdJl7ob3dernMZcy3vE19T4CzdaNsvpGRka3ddnOXLexPjYg7JR0MfA14V0RcM9f4rVYrJiYmZh02Pj7O8PBwybj9t2bDFtYP7WHTZJVequ5s37iuo/F63XZrNmzp2bTKtl2ny9wL7e3Wy2UuY77lbep7ApytG2XzSeq6sJfqiomIO4v/dwOXAS/sZuZmZtZ7HRd2SSsk7T99G3gFcFO/gpmZWTVl+hsGgcskTT/vwoi4oi+pzMysso4Le0TcBjyvj1nMzKwHfLqjmVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZKF3ZJSyTdKOnyfgQyM7PuVNliPwvY1usgZmbWG6UKu6RDgXXA3/YnjpmZdUsR0fnI0sXAR4D9gT+OiJNmGWcUGAUYHBxcOzY2Nuu0pqamGBgYqJK5ryZ37GJwOezcXXeSuTU5X9lsQ6tX9i/MDO3r3OSOXYs233bzLW9T3xPgbN0om29kZGRrRLS6mefSTkeUdBJwd0RslTQ813gRsRnYDNBqtWJ4ePZRx8fHmWtYnU7fsIX1Q3vYNNlx0yy6Jucrm237qcP9CzND+zp3+oYtizbfdvMtb1PfE+Bs3agjX5mumOOB10jaDowBJ0j6fF9SmZlZZR0X9oj4LxFxaESsAV4HXBkRp/UtmZmZVeLz2M3MMlOpozYixoHxniYxM7Oe8Ba7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8tMx4Vd0r6Svi3pu5JulvShfgYzM7NqyvyC0sPACRExJWkZ8A1J/xQR/9KnbGZmVkHHhT0iApgq7i4r/qIfoczMrDqlet3hyNISYCvwbOC8iPiTWcYZBUYBBgcH146Njc06rampKQYGBuac1+SOXR3n6rXB5bBzd22zX1CT85XNNrR6Zf/CzNC+ztW1fs23vAu9J+r0RM3WhNe5bNuNjIxsjYhWN/MvVdh/8STpQOAy4F0RcdNc47VarZiYmJh12Pj4OMPDw3POY82GLaVz9cr6oT1smqz0O9+Losn5ymbbvnFdH9P8e+3rXF3r13zLu9B7ok5P1GxNeJ3Ltp2krgt7pbNiIuJ+YBw4sZuZm5lZ75U5K+agYksdScuBlwM/6FcwMzOrpsz+/FOAC4p+9r2AL0bE5f2JZWZmVZU5K+Z7wDF9zGJmZj3gK0/NzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMlPmN0+fJukqSdsk3SzprH4GMzOzasr85ukeYH1E3CBpf2CrpK9FxPf7lM3MzCroeIs9Iv4tIm4obj8AbANW9yuYmZlVo4go/yRpDXANcHRE/GzGsFFgFGBwcHDt2NjYrNOYmppiYGBgznlM7thVOlevDC6Hnbtrm/2CmpzP2arrd76h1SsrP3eh92ud5stWVx1pb+uybTcyMrI1IlrdzL90YZc0AFwNnBMRl843bqvViomJiVmHjY+PMzw8POdz12zYUipXL60f2sOmyTK9VIuryfmcrbp+59u+cV3l5y70fq3TfNnqqiPtbV227SR1XdhLnRUjaRlwCfCFhYq6mZnVo8xZMQI+DWyLiD/tXyQzM+tGmS3244E3AidI+k7x96o+5TIzs4o67tCLiG8A6mMWMzPrAV95amaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpaZMr95+hlJd0u6qZ+BzMysO2W22M8HTuxTDjMz65GOC3tEXAPc18csZmbWA4qIzkeW1gCXR8TR84wzCowCDA4Orh0bG5t1vKmpKQYGBuac1+SOXR3n6rXB5bBzd22zX1CT8zlbdf3ON7R6ZeXnLvR+rdN82eqqI+1tXbbtRkZGtkZEq5v597ywt2u1WjExMTHrsPHxcYaHh+d87poNWzrO1Wvrh/awaXJpbfNfSJPzOVt1/c63feO6ys9d6P1ap/my1VVH2tu6bNtJ6rqw+6wYM7PMuLCbmWWmzOmOFwHfAo6QdIekt/YvlpmZVdVxh15EvL6fQczMrDfcFWNmlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWmVKFXdKJkm6R9ENJG/oVyszMqivzm6dLgPOAVwJHAa+XdFS/gpmZWTVltthfCPwwIm6LiEeAMeC3+xPLzMyqUkR0NqJ0MnBiRJxR3H8j8KKIeOeM8UaB0eLuEcAtc0xyFXBPldCLoMnZoNn5nK26JudzturK5jssIg7qZoZLS4yrWR77pU+FiNgMbF5wYtJERLRKzH/RNDkbNDufs1XX5HzOVl0d+cp0xdwBPK3t/qHAnb2NY2Zm3SpT2K8HDpf0DEl7A68D/rE/sczMrKqOu2IiYo+kdwJfAZYAn4mIm7uY94LdNTVqcjZodj5nq67J+ZytukXP1/HBUzMze2LwladmZplxYTczy4wLu5lZZp6QhV3Sakmr684xG0nPlPQeSSfUnWWmJmeDZudztuqanK/J2aB6vidUYZe0RtLVwBXAxyW9pO5M7ST9R+BrpO/SeZukt9cc6ReanA2anc/ZqmtyviZngy7zRUSj/4B9226/FvhEcfvNwD8AQ8V91ZDtBOAZ0/MHPgCcVtx/EfAlYLiOfE3O1vR8zpZnviZn63W+Rm6xSzpA0l9JuhX4hKTDikG/A/youD0G/BA4Y/ppi5jvKEnfA/4b8FlJJ0Rq7aOAQwAi4jrgm8BbFjNfk7M1PZ+z5Zmvydn6la+RhR04EdiXtGCPAB+QtJy0W/JqgIh4GLgYeElx/7F+hZF0qKQD2h46BbgkIl5K+oB5g6TDgQun8xUuA46WtE+/8jU5W9PzOVue+ZqcbbHy1VbYlSyV9FZJX5d0lqRnFYOfDTwSEXuAPwN+CpwGfBV4iqRfK8a7FfixpOP6lPFISV8GvgF8WNL01xT/P2C/4vYXgbuAdaRP1Ce37WHcR/p2y+f9KmVrej5nyzNfk7Mtdr7aCnuxq/GbwJuAjwH7AH9TDL4LuLv4ZPoxaWGeRWqA7/P41wIvA+4tHu8JSSva7j4fuCMi1gBXAp8oHr8PeFjS/hFxH/CvwFOLHN8E3luMtzfwc2B77tmans/Z8szX5Gx15lu0wi7pOEkflXR6cV/AkcAVEfGliPgYcJikFwM7SJ9gRxZP3wYMFI/9BfAqSa8mfSgMAt/tMtuTJJ0v6Xpgo6SDinxDwLWSFBH/CNwvaR1pT2H/YjjF/YOBx0h7GAdL+hvgImBPRNydY7am53O26pqcr8nZmpJvUQq7pOcAfwk8APyepPcW814NPFAsNMD5wBtIhXoP8OLi8RtIR4wfiohrgA3A6cDxwH+PiMfaplHFS4v5vYp0UOJs4ADSl50dUuxdAFxQ5Pt2sSyvBIiIbxXTWBoR24AzgZuB/xkRb6E7Tc7W9HzOlme+JmdrRr65Tpep+kfasj6DtNuxtHjsT4Gzitst4JPAycDLga+0PfdppF0VSIX8RtKvMB0D/G/gKW3jlj4dqWjYM4GrSd05q4rHvwi8u7j9DGBjMfxYUn/YkrZl+0kxndWkPYl3Ap8FPgWs6KLdGput6fmcza+r2+7f//V0i13S80kHOH8b+CDwvmLQDtJvpkL65LkW+F3gn4FDJD1X0rJI/ek7JL0kIq4kfd3lR4FLgYsi4t+m5xVFy5R0EvAa4EPAcaS+fUhn20zvHfwY+Drwyoi4nvSJO1LMcwq4Djg2InYAbyR1Bd0FvC8iHiwbqG1P49VNyzaD266axrUbuO26yfZEaLsyP433SyS9EDgc+GpE/IS0NX5rRJwu6QXAOZJawDjwnyTtFxEPSfou8HukczQvBP4A+KSk3cAkcHsxi78CLoyIXSUyKSJC0rGk3ZyvA1sinR7568BtEXGlpNtJV6++AtgK/I6kVRFxj6R/BR6U9HTgz4HTJB1M+tWoe0m7TkTEBDBRod1apL2aB4CPA3cDz6w7m9uuWrYnQru57fJru/mU2mJXskzSmyTdSOrYPxCYLrw/B7YXW983kHYtjgMe4vFTeAAeJe2CHELaKr+J1L9+NXBPRNwBaau8YlF/KfAZ0lHllwMfKUZ5DLhV0vKIuL3I91zSi3Un6XzS6eVYQmqfS4qMpwJrgc1R8RxXSSslfbaY5u3AuRFxt6S9SJ/kdWZbUrTdb5J2BRvTdsV6NyDpfBrWdsU8Q9IwzVzn9pG0oqFtd0DD225A0r6SLqBhbbegTvprgBXAi4vbBxbBPjnLeGeRLoNdXdw/mdSffhjpKwCuLh7fl9QNs6rtuccAe3eSZ8Y89wPexuNb/suAPwLeUQx/EvC9YvqnkPq71hTDTiqWZVVxexJYSerf/3J7HmCvLrJdRLpibIDUtXRm2zjTxyHeCfyPxcrW9rqeQVrZ1pMO8DSl7aazXVqsVwc1rO32B7aQfkkM4D1NaLcZ+b4M/HVx/2PA2+puO9J74s2k9/8lTWu7tnxXAn9fPNaY9a7TvwW32CWdDdwGbJE0GBH3k/qF7iz6xl+jxy8Q+hbpAOj0hUbXkg6iPhQRFwA/lfQ50kHRW4Bf9CFFxI0R8chCeWZkOwS4HBgGPkc6QPFa0l7CnmK6PyUdeH03qe/rYB4/jfIa0rn0j0TE5cCnSVeznkc6Yv1oW75Sn6ozsv0d8PYi263AEZI2FltRv690wdUVpD2Yvmcr8q0gvblOIF0/8ArScY9jSVtKdbZde7bNpLMFXku6huE36m67wnLStRfPkrSKtM4vKaZZS7vNkm9v0rr2VFIXx9GSPlJX20laRjrGdjLw8Yj43WLQMW3TrK3tZuT7WERMb3FPAkfV2XaldfAJNkzavfhb4D3FY8eSitYdRfALgU3FsHOAD7c9/3rgmOL2PqRTgI7txacSaeV9Udv900lbJm8Gvt32+FOBO4vb7yBdtvuk4vlfAp7eNu6qPmV7E+lI968Df1/8vR74X6Rz+RctW9v0Dmy7/Z9Jb6ZT6267WbL9MemUsWc2qO3eTOprfT/wVtKBtOvrbrdZ8r2PtMezqgltR9oDO3XGY6cA1zWh7ebI9/QiQ+3rXcfL0cGCTp+acwowXtxeRtqaWlncP4y0tX4saRfwYtKW1j+RPqn26Uv41Mcl+MVvt76Ax7t77iWdMzo97tcoCi1p9+mrxTh/skjZjgG+Mb3ito23jHRw+YTi/jn9zjYj5wGk4xs7gQ8X9+8FButqu1my3VXMdwVFN19dbdf2er6F1M32WuALxWP31N1uc+QbKx5rP124lvWO1EVxK7CpmP8HivpxH3BwA9a59nxXkb6Y69C617vSy1FigZ9MulDoOcX9pTOGnw+cPL0CkboezqRPRX2OlfkCHj9f/nPAR4vbv0ba43h62wtzNG1fCbxI2d7R/lhx+5Ci7Z672NnaMvwh6XzbzaR+7W8WbzjV2XYzsp1HOq3s2U1oO9JXRi8h9aFeTdoyvgl4f93r3Cz5/pl0htkLGtJ2XyHtgT2NtBV8FmnDsCnrXHu+L5Au/T+8CW3X6d900emIpE8BP4uIDcX9vUjnXb4DeA5wSpTsJ+8VSYeS+rTeFRG3Kn2h2GiRazXwnejNVWXdZHt7RNxWPHYMqVtqXZHtD+vI1k7pOoQzSW+yI0kr66HU2HZt2Y4mvdn+nHSW1UnU1HaSBkjdHPuQ2uk3SBeenE3aUj6cGtttlnyHk45P/BbpmNfLSO1Xy3qn4rTn4vbzSO/Ta0mX1Ne+zs3IdzTpSvdzSd80W9t6V0bZ89g3A+cWBxmOJK3Ex5NelLPrKuqFYyjOgZd0Bqn//2xSF9IPIp1+WXe2HxXZbietHHtIW/E31pit3b2kg4Dvi4i/k3QacHND8t1P6ie+ifS6LqO+tttDOnviUdKW+s9J6/8k8N4GtNtc+R6W9BpSwa9tvZsumoX7Sced3h8RFzag7Wbme4C08boN+K/Uu951rOwW++tIB0ofJn3j2JURcUufspUi6VrSwbXtpHNIPxQR36s1VGFGtruADQ1qt5WkLbg3kL7/fjNwXkQ8Ou8TF8Es2T4dEZvqTfXLigtPpvuy76o7z0xFvpOBz0Y666TuPPuQfnPhjaQ96r8EPhXpa7prN0u+zRHxZ/WmKqfjwi7puaTzOS8mHSzq2VfldqvYg/ggaUv485GuWmuEJmcDkLSU1P3yMClfk17XxmaDdFEX8FiU2TpaRE3OJ+lM0mm1n2va6wrNz7eQUlvsZmbWfE39aTwzM6vIhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlpn/D0QBdzhJVkBDAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -342,9 +2473,14 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, "outputs": [], "source": [] } @@ -365,7 +2501,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4,